- Research
- Open access
- Published:
Common lizard microhabitat selection varies by sex, parity mode, and colouration
BMC Ecology and Evolution volume 23, Article number: 47 (2023)
Abstract
Background
Animals select and interact with their environment in various ways, including to ensure their physiology is at its optimal capacity, access to prey is possible, and predators can be avoided. Often conflicting, the balance of choices made may vary depending on an individual’s life-history and condition. The common lizard (Zootoca vivipara) has egg-laying and live-bearing lineages and displays a variety of dorsal patterns and colouration. How colouration and reproductive mode affect habitat selection decisions on the landscape is not known. In this study, we first tested if co-occurring male and female viviparous and oviparous common lizards differ in their microhabitat selection. Second, we tested if the dorsal colouration of an individual lizard matched its basking site choice within the microhabitat where it was encountered, which could be related to camouflage and crypsis.
Results
We found that site use differed from the habitat otherwise available, suggesting lizards actively choose the composition and structure of their microhabitat. Females were found in areas with more wood and less bare ground compared to males; we speculate that this may be for better camouflage and reducing predation risk during pregnancy, when females are less mobile. Microhabitat use also differed by parity mode: viviparous lizards were found in areas with more density of flowering plants, while oviparous lizards were found in areas that were wetter and had more moss. This may relate to differing habitat preferences of viviparous vs. oviparous for clutch lay sites. We found that an individual’s dorsal colouration matched that of the substrate of its basking site. This could indicate that individuals may choose their basking site to optimise camouflage within microhabitat. Further, all individuals were found basking in areas close to cover, which we expect could be used to escape predation.
Conclusions
Our study suggests that common lizards may actively choose their microhabitat and basking site, balancing physiological requirements, escape response and camouflage as a tactic for predator avoidance. This varies for parity modes, sexes, and dorsal colourations, suggesting that individual optimisation strategies are influenced by inter-individual variation within populations as well as determined by evolutionary differences associated with life history.
Background
Within their environments, individuals of any animal species are presented with many options for where to spend their time. These microhabitat selection decisions mainly depend on the trade-offs between finding food, engaging in social interactions, and avoiding predators [1,2,3,4]. Among these, predator avoidance may play the most important role, as an incautious selection can lead to immediate death [3]. To avoid predation, an individual is usually faced with two choices: hiding or fleeing. Camouflage, the similarity between the colouration and pattern of an animal and its surroundings, is a striking example of adaptation that has become a model system within evolutionary research [5,6,7,8,9]. Studies on camouflage in wild animals suggest that prey animals inhabiting a variety of coloured microhabitats will either compromise between different microhabitats, or specialise in one microhabitat to optimise predator avoidance [7, 10]. For example, bark-resting moths such as the peppered moth (Biston betularia) show a preference for microhabitats with a similar colouration to themselves [11, 12] and will even reposition their body to enhance camouflage [13,14,15]. Camouflage may also depend on the variety of different microhabitats encountered by an individual. For example, desert lizards (Sceloporus magister) have compromised on a coloration that is intermediate between two or possibly more of the microhabitats that they use [7].
In squamates, dorsal patterning and colouration has mainly evolved to match the surrounding substrate for crypsis [16,17,18]. However, colour and patterning may also have evolved as a response to sexual signalling and environmental pressures, including temperature and protection from ultra-violet radiation [19,20,21,22,23,24]. For example, in the viper Montivipera raddei species complex, dorsal pattern colouration varies in association with thermal factors, including temperature and solar radiation [25]. However, predator avoidance tactics might change during an individual’s lifetime, depending on ontogenetic changes in colouration, sprint speed (escape performance), or changes in behaviour due to reproduction and mate acquisition [26,27,28]. Overall, studies find a complex relationship between colouration, patterning, sex, and animal personalities, suggesting that there are interdependent trade-offs in the evolution of this variation [23, 29, 30].
Providing opportunities for camouflage is an important consideration for habitat selection, but microhabitats critically also offer protection from predators [31, 32], opportunities for encountering prey [4, 33], or thermoregulatory benefits for ectotherms [2, 34]. For example, vegetation cover and short flee distance to hiding spots allow animals to avoid detection and decrease predation risk [32, 35]. On the other hand, open basking sites provide optimal opportunities to maximise exposure to sunlight and thus have thermoregulatory benefits for ectotherms [36, 37], but are more exposed to potential predators [38]. Therefore, individuals need to carefully choose the amount of time they spend in different types of microhabitat. In addition, these choices may depend on the sex of an individual and the season. For example, females are usually less mobile during pregnancy, and therefore need to adjust their behaviour as their flight speed reduces [39, 40].
While it has been well demonstrated that lizards select microhabitats with particular characteristics (e.g., for thermoregulation [41,42,43,44] or predator avoidance [1, 45]), few studies have tested whether individuals within the same species that show intraspecific variation in appearance actively choose specific backgrounds within a microhabitat. Recent studies indeed suggest that, in some species, individuals may actively choose a habitat that matches their own colouration. For example, in ground-nesting birds, individuals match their own colouration and their eggs to the habitat they nest in [46] and escape more readily when their own match to the background – but not their clutch’s – is poor [47]. In Aegean wall lizards (Podarcis erhardii) inhabiting Greek islands, individuals match their own colouration to the background, with the effect being stronger in females and in areas with higher predation risk [15]. Therefore individual decision-making on site use likely involves the local environment availability and an organism’s colour variation.
The common lizard (Zootoca vivipara) is a small, ground-dwelling lizard found across Eurasia. In addition to having the largest distributional range of any terrestrial reptile, it shows striking variation in reproductive aspects of life-history: some populations are egg-laying (oviparous), while others are live-bearing (viviparous) [48, 49]. Common lizard reproductive season starts in early spring and ends in summer, and during this time gestating females may double in weight [50]. This period is also when common lizards are most vulnerable to predators that often rely on visual cues, such as birds of prey, snakes and some mammals.
Common lizards show considerable variation in inter-individual colouration and dorsal patterning [51], with reticulated or linear patterns of dark on lighter brown background, or an intermediate pattern between reticulated and linear [23]. An individual’s dorsal patterning is established during the first year and is then stable throughout its lifetime [23]. Common lizards do not have active colour change to match or camouflage against particular substrates, as known in other lizards such as chameleons [52]. Features such as body condition and heterozygosity play important roles in common lizard mate choice [53, 54] but dorsal colour or pattern do not seem to be crucial cues. Instead the frequencies at which reticulated and linear forms occur in a population vary depending on environmental conditions and on genetic inheritance [23]. However, the association of dorsal colour and pattern on the landscape has not been investigated in detail either for males or females.
In this study we assessed microhabitat selection in syntopically occurring oviparous and viviparous common lizards during the reproductive season. We first tested if common lizard microhabitat selection was specific in the context of the environment available to them. In addition, we checked if their choice was associated with parity mode (oviparous or viviparous) and sex. Finally, we tested if common lizard dorsal colouration and patterning was associated with colouration of the basking site while accounting for size and sex differences. The hypotheses we tested were: (i) particular features of microhabitat are preferred by common lizards, (ii) oviparous and viviparous common lizards differ in microhabitat selection, (iii) pregnant females are encountered in a different microhabitat than males, (iv) dorsal colouration matches the colouration of the lizard’s basking site, and (v) dorsal patterning is associated with basking site colouration.
Results
Lizard selected microhabitat differs between sexes and parity mode
Common lizards apparently chose specific microhabitats within the sampled study area. This is shown by significant differences between the relative proportions of microhabitat categories available in the area and the microhabitat used by the lizards (Figs. 1 and 2). This was exemplified by DIC (deviance information criterion) estimates for the preferred simple lizard model (where common lizards are assumed to be a homogenous group, Table 1, rows 2 and 3) being at least 5 lower than the null model, which assumed no difference between sites with common lizards and random sites for all nine habitat types modelled (Table 1, row 1). In particular, common lizards tended to be present in areas with higher proportions of angiosperms, leaf litter, moss, water and wood, and lower proportions of cover, rocks, grass and bare ground than were randomly available in their surroundings (Table 1; Additional file 1 (Table S1); Fig. 2). Where a simple model distinguishing between lizards and randomly available sites was preferred, common lizards were present in areas with a lower proportion of grass (\({\mu }_{lizard}=0.55\) and \({\mu }_{random}=0.70\)) and higher proportions of leaf litter (\({\mu }_{lizard}=0.18\) and \({\mu }_{random}=0.10\)) than randomly available. When present, rock was found in similar proportions across groups (\({\mu }_{lizard}=0.27\) and \({\mu }_{random}=0.29\)), but was more likely to be completely absent in the vicinity of lizards (\({\alpha }_{lizard}=0.89\) and \({\alpha }_{random}=0.83\)).
For six of the nine habitat types, more complex models based on differences between lizards were preferred (Table 1). Males lizards were more likely to be observed with no bare ground (proportion = 0) in the vicinity (\({\alpha }_{random}=0.67\)) compared with both females and random availability (\({\alpha }_{random}=0.58\) and \({\alpha }_{female}=0.57\)) (Additional file 1 (Table S1). When present, the proportion of wood in the vicinity of females of both parity modes was higher (\({\mu }_{female}=0.47)\) compared with male lizards (\({\mu }_{female}=0.42)\) or random (\({\mu }_{female}=0.22\)). Viviparous common lizards were found more often in association with angiosperms (\({\alpha }_{viviparous}=0.65\)) compared with oviparous lizards (\({\alpha }_{oviparous}=0.87\)). Oviparous lizards were less likely to be found in the absence of water (\({\alpha }_{oviparous}=0.91\)) than viviparous lizards or randomly available (\({\alpha }_{viviparous}=0.99\) and \({\alpha }_{random}=0.99\)) and when present, water made up a higher proportion of an oviparous lizard’s surroundings (\({\mu }_{oviparous}=0.19\) while \({\mu }_{viviparous}\) and \({\mu }_{random}=0.14\)). Oviparous common lizards were also less likely to be found in the absence of moss (\({\alpha }_{oviparous}=0.71\)) than viviparous lizards or randomly available (\({\alpha }_{viviparous}=0.86\) and \({\alpha }_{random}=0.88\)). Finally, viviparous lizards were more likely to have a total presence or absence of cover (\({\alpha }_{viviparous}=0.12)\) with a proportion of 1 more likely than absence (\({\gamma }_{viviparous}=0.59\)), meanwhile oviparous common lizards exhibited a bimodal distribution in their preference of how much cover was available in their habitat, while viviparous had a more unimodal distribution (Fig. 2).
Dorsal colouration and patterning
Dorsal colouration was associated with parity mode, sex, body length, and body weight (Table 2). More specifically, viviparous and larger individuals tended to have a higher hue (Table 2; Additional file 2 (Table S2); Additional file 3 (Fig. S1)). Moreover, viviparous, male, larger, and lighter-weight lizards tended to have a more intense colour (model 1 (saturation); Table 2; Additional file 2 (Table S2)). Finally, oviparous had a lighter dorsal colouration compared to viviparous individuals (estimate = -3.12, t1,250 = -3.34, P < 0.001; Table 2; Additional file 2 (Table S2)).
Dorsal patterning was associated with body length and sex (model 2; Table 2; Fig. 3; Additional file 2 (Table S2)). Specifically, a larger proportion of males had a reticulate pattern compared to females (estimate = 1.068, z1,387 = 3.05, P < 0.01; Table 2; Fig. 3B) and smaller lizards were more likely to have a linear pattern (estimate = 0.071, z1,387 = 2.79, P < 0.01; Additional file 2 (Table S2)). Individuals with a reticulate pattern generally had a lighter dorsal colouration (model 3; F2,299 = 3.28, R2 = 0.015, P < 0.05; Fig. 3C; Additional file 2 (Table S2)).
Common lizard colour matches microhabitat and basking site
Common lizard dorsal hue ranged from 27 to 51 (mean: 36.8), basking sites from 20 to 69 (mean: 47.3) and habitat from 33 to 69 (mean: 54.9). Hue was positively correlated between habitat and basking site (model 4; estimate = 0.39, t1,160 = 8.75, R2 = 0.331, P < 0.001;) and more variable in basking sites than habitat (standard deviation: 10.6 vs. 7.2). Common lizard lightness ranged from 25.1 to 64.9, basking site from 12.9 to 79.4, and overall microhabitat from 17.3 to 47.6 (Additional file 4 (Fig. S2)). Basking site lightness and overall habitat lightness were significantly correlated (model 5; estimate = 0.33, t1,160 = 10.81, R2 = 0.427, P < 0.001; Table 2; Additional file 2 (Table S2)), but the distributions differ, with basking sites being lighter (mean: 36.8 vs. 30.1) and more variable in lightness than habitat (standard deviation: 12.7 vs. 6.6; Additional file 4 (Fig. S2)).
Dorsal colouration was associated with the colour of the basking microhabitat selected by the lizard (model 6). Specifically, the hue of the lizard’s dorsal area correlated with the hue of the habitat (estimate = 0.21, t1,162 = 5.29, R2 = 0.151, P < 0.001; Table 2; Fig. 4) and the basking site (estimate = 0.09, t1,162 = 3.00, R2 = 0.05, P < 0.01; Table 2; Fig. 4). Therefore, lizards generally tended to prefer habitats and basking sites that matched their own colouration. Moreover, lizards tended to choose basking sites that were slightly darker than their own colouration (estimate = -0.10, t1,162 = -2.14, R2 = 0.03, P < 0.05; Table 2; Fig. 4). In contrast to colouration, dorsal patterning, i.e. being linear or reticulated, was not associated with basking site or habitat colouration (N = 161, P > 0.1; Additional file 2 (Table S2)).
Discussion
Our findings suggest that common lizards in their natural environment actively select their microhabitat, presumably to optimise thermoregulatory and physiological requirements while minimising risk. Lizards chose all habitat components in markedly different proportions than what was available to them in the immediate area. We showed that common lizards generally avoid bare rocks and bare ground. In this way, we suggest they can avoid detection and the risk of predation in such exposed habitats, an important consideration given that they are rather slow, ground-dwelling lizards. Common lizards have relatively high humidity requirements for their environment [55, 56], and this was reflected in their preference for a wetter habitat with more moss. Moreover, leaf litter and angiosperms were preferred, presumably providing options for escaping while being less visible to predators.
As basking sites, common lizards often use wood, which was also enriched in their microhabitat relative to the immediate area. Wood generally matches well with the lizard’s background colouration and also offers crevices for quick access to escape refuges [27]. In addition, wood heats up quickly and therefore provides optimal options for thermoregulation [57]. Microhabitats with less vegetation cover available than what is randomly distributed in the environment were also preferred, presumably because common lizards require basking sites and taller vegetation reduces the amount of solar radiation that can be used for thermoregulation.
In addition to these general preferences, microhabitat selection depended on an individual’s life history: site selection differed between males and females, and oviparous and viviparous lizards. Females preferred habitats with good basking options (more wood) and less open space (‘bare ground’). During the reproductive season, females are pregnant and reduced in speed [40, 58, 59], and therefore may require relatively ‘safer’ microhabitats where they can attain their optimal body temperature quickly but also escape to a hiding spot within short distances. By staying in close proximity to a refuge and exhibiting a low flight tolerance [27], female common lizards may offset their reduced speed during pregnancy to avoid predation. By performing running performance tests in addition to microhabitat choice experiments, this potential sex-bias could be tested in the future using a larger and more balanced sample size.
Parity mode was also associated with common lizard microhabitat selection. Specifically, oviparous common lizards were found in wetter microhabitats with a higher proportion of moss and a lower proportion of flowering plants. This confirms a previous study which found that oviparous common lizards are found in wetter habitats [60]. We propose that the reason for this may be that for optimal development, egg clutches need a humid environment. Egg clutches of oviparous common lizards have been found below moss and within humid dead wood (personal observation), and in captivity egg clutches are usually deposited below moisturised moss when given the option [50]. Viviparous common lizards were found in microhabitats with a larger proportion of angiosperms, while other components did not differ between parity modes. In summary, we show that an individual’s sex and life history (male vs. female, oviparous vs. viviparous) are associated with microhabitat choice in common lizards in their natural environment, presumably to match physiological requirements and trade-offs with predation risk. This suggests long-term evolutionary differences play out in the behavioural strategies of individuals on the landscape.
Our findings further suggest that common lizards select particular microhabitats and basking sites compared to the background matrix and that this is associated with their dorsal colour. Individuals tended to choose a habitat and basking site colour closer to their own colouration (measured as hue) than other available options. This result is in agreement with other studies suggesting lizards match background to their own colouration in order to reduce detection by predators [15, 61, 62]. However, in respect to avian predators our result should be interpreted cautiously, as we did not measure the full spectrum of the lizard and its background. Both lizards and birds cover UV light in their visual spectra, which we did not measure here. In contrast to the positive correlation between the lizard’s hue and its habitat, the lightness of the lizard’s basking site tended to be darker than its own colouration. A potential explanation could be that lizards optimise their thermal environment during basking. As darker surfaces allow for more efficient thermoregulation, lighter individuals, which presumably require more time to heat up compared to darker individuals [63], may offset this difference by choosing darker surfaces while basking. In summary, our results suggest that lizard individuals are faced with a trade-off between heating efficiently and matching their background.
How common lizards choose their habitat and basking site in relation to their own colouration is an intriguing question. It is possible that lizards can visually assess their own dorsal colouration by curling, a behaviour they typically also adopt while basking. Alternatively, the locus or loci for dorsal colouration may be genetically linked to behavioural preferences for habitats that best match their colour and pattern. Another alternative may be that lizards learn to optimise their camouflaging behaviour from previous encounters with predators [64]. Basking sites differed from the otherwise immediate microhabitat in colouration and in composition of substrate and vegetation, strongly suggesting they are actively selected relative to the immediately available microhabitat.
Our findings also demonstrate the variation in common lizard colouration among individuals. In particular, colour intensity differed significantly between individuals of different parity modes, sex, body length and body weight. Oviparous individuals differed from viviparous individuals in all three dorsal colour variables: viviparous individuals had larger values in hue and saturation, and a lower value in lightness. While these differences were subtle, it confirms impressions from the field that oviparous individuals tended to be slightly lighter than viviparous. Moreover, the hue was shifted towards larger values in larger individuals, visible as a change from a darker brown towards a lighter brown. This is consistent with common lizards undergoing colour changes during ontogeny, while the correlation with body weight and overall fitness warrants further investigation. Reticulate-patterned individuals tended to be lighter in dorsal colouration. We also found that pattern is significantly associated with sex, with reticulated patterns being more common in males and linear patterns being more common in females. That dorsal patterning differs between sexes is in agreement with another study showing a higher proportion of linear patterns in females in wall lizard (Podarcis muralis) populations [24] and in a different common lizard lineage [23]. Previous research comparing dorsal patterns of common lizards and exploratory behaviour found linear individuals seem to be bolder in transitioning between habitat patches without cover, and have higher immigration rates compared to reticulate individuals [23, 29]. Therefore, there is the potential that complex behavioural traits are associated with dorsal patterning, though here we did not find this to be associated with different basking site selection.
As dorsal patterning did not have a significant effect on microhabitat selection, we suggest that the relatively darker linear morph being more abundant in females might be a thermoregulatory adaptation. A darker linear patterning may warm up the ventral column quicker during basking and improve thermoregulation [23, 24]. A higher frequency of the linear pattern in females might be explained by an increase in energy from thermoregulation aiding to escape from predators whilst gravid [24, 65, 66]. In general, the grain and geometry of a pattern, together with body colouration, may represent random samples of those visual features of the background and aid in camouflage whilst basking [61, 67]. Patterning may also cause disruptive camouflage by breaking body lines or merging body margins to match with the substrate [61, 68,69,70]. Although we did not find an association between dorsal patterning and basking site in common lizards at this relatively small spatial scale, patterning might be related to habitat on a larger scale, for example when habitats differ more substantially from each other (e.g., alpine meadows compared to peatland meadows).
Conclusions
Overall, our key findings suggest that common lizards selectively choose their microhabitat and that this is non-random relative to the immediately available sites, presumably to balance life-history dependent physiological requirements, camouflage, and predator avoidance. Common lizards are preferentially found in backgrounds that enhance matching, potentially for camouflage against predators in their natural microhabitats whilst staying close to cover to enable quick escape. Further work is needed to understand the active choice in these site selections by lizards and whether camouflaged background matching benefits survival. Other studies have suggested that genetic control alongside visual input may underlie adaptive background choices, and individuals may learn to discriminate between camouflaging and non-camouflaging backgrounds [15]. Although this study tests true individual colours (as captured by camera RGB values) and does not rely on colours being indicated by the human visual system, other studies have relied on human-oriented indices [71,72,73]. Measures of crypsis based on the human visual system may be biased or differ from that of other animals, as different groups of animals can vary greatly in their visual systems [74,75,76,77]. Finally, future studies of reptile crypsis will benefit from using spectral reflectance data covering a wider spectrum than digital photography, and linking the background and prey colouration directly to the visual system of potential predator(s) [15, 61, 78,79,80].
Materials and methods
Study area and sampling
The study was carried out in the Carinthian Alps in Austria in the Gailtal valley. The study area covers six sites that are separated by visible barriers (e.g. rivers, roads, forests) and a total area size of approximately 0.3 km2, with an altitude range from 1350 to 1550 m. The area is used for grazing cows and goats during summer months (June-September) and the habitat is dominated by alpine meadow. The study area is one of the few contact zones between egg-laying and live-bearing common lizards, two genetically differentiated and distinct phylogenetic lineages [48, 49, 60].
A total of 394 individuals were caught by hand from May to July 2016 and sexed by the presence of a penile bulge in males and its absence in females. Individuals were photographed under standardised conditions to quantify colour information accurately (see next section). Snout to vent lengths (SVL) and tail lengths (TL) were measured using digital callipers (Moore and Wright) and weights were measured using a Pesola Micro-Line Spring Scale. All female lizards were kept in individual plastic terraria (28 × 19 × 14 cm), provided with hiding places, wet moss, water and food (mealworms, crickets) ad libitum [50]. Following laying eggs or giving birth, a tail clip from each female was taken for genetic analyses and she was then released at the point of capture. Live-born young were released with their mother, while egg clutches were incubated under 24 °C and released at the point of their mother’s capture upon hatching. Parity mode was assessed by phenotype (i.e., if individuals laid a clutch of eggs or gave birth to living young) and SNP genotyping [49]. For genotyping, ddRADSeq libraries were constructed and sequenced on an Illumina HiSeq 4000 machine at 2 × 150 bp read length. Individuals were demultiplexed using STACKS [81] and mapped to the common lizard reference genome [82]. SNPs were filtered and extracted by a required presence of > 50% across individuals and a minimum allele frequency of 10%. An ADMIXTURE analysis was run on all individuals using K = 2 to account for the oviparous and viviparous phylogenetic lineages encountered in the area. If the membership value (Q) fell between 0.1 and 0.9, individuals were recorded as admixed and excluded. Higher or lower values were recorded as purely oviparous or viviparous, respectively.
Colouration and pattern data generation
Images were taken in RAW format to preserve colour information under standard settings (fixed ISO: 160, aperture: 8, exposure compensation: 0) with a Canon 70D and a 60 mm fixed lens. Each image included an ‘X-rite ColorChecker’ board to standardise colouration between images and capture accurate colour [83, 84]. First, a Digital Negative (DNG) copy of each image was made and then run through the colour checker software (X-rite ColorChecker) which establishes an accurate colour foundation in form of a profile. Next, this profile is applied to the original image before analysing to ensure standard and accurate colour results. Images were processed using Adobe Photoshop CC 2015. After applying the profile to the image, white balance was selected for the off-white square on the X-rite colour checker board and subsequently the red-green-blue (RGB) linearity for the selected square values were set to 240 nanometres (nm) [85]. Finally, RGB values were recorded for the area surrounding the dorsal stripe, where colour sampling area included the dorsum from front legs to hind legs; flanks were excluded. Hue, saturation and lightness (HSL) values were calculated to assess the colour. Hue is an approximation of the main colour (from red to blue), while saturation measures its intensity, and lightness reflects how dark or light the colour is by quantifying the amount of black or white added to the hue. Our interpretation focused mostly on hue (H) and lightness (L). We note that while our approach allows for the unbiased extraction of colours within the human visual spectrum (300–700 nm), it ignores colours within the UV spectrum. Therefore, these colour profiles are representative for mammalian predators but the spectral sensitivity range of lizards (including Zootoca vivipara [86]) and birds [75, 87] is not fully covered.
From the images, the dorsal stripe was categorised by eye as either linear, reticulated or intermediate for the anterior and posterior body part. As variation between the anterior and posterior portions of the lizard’s body has previously been reported [24], dorsal patterning was recorded separately for these two areas. The dorsal pattern of each of the 394 individuals (264 females and 130 males) was categorised and dorsal colouration quantified.
Microhabitat use
Habitat measurements were carried out in the same area where each individual was found basking and captured. The basking site was defined as the location where the lizard was first spotted and stationary. An image of the habitat was taken from above immediately after lizard capture using a Canon 70D and 60 mm lens. Each image included a 0.5 × 0.5 m quadrat square, with the basking location of approximately 0.15 × 0.15 m placed in the middle of the quadrat, and an X-rite ColorChecker board outside the quadrat. Habitat images were processed in Adobe Photoshop using the same colour-sampling methodology as for the lizard dorsal images, with profiles prepared using the X-rite ColorChecker software. The RGB values were recorded for the 0.15 × 0.15 m basking site in the centre of the quadrat. Out of all sampled individuals (N = 394), a total of 165 (120 females and 45 males) were sampled for both dorsal colouration and basking site colouration.
Habitat was categorised as made up of areas that contained grass, moss, leaf litter, angiosperm, wood, rock, bare soil, water, or ‘other’ substrate within the quadrat. Habitat category proportions were inferred using the image and a grid (25 grid squares at 3 cm x 3 cm) overlaid on the habitat image on the computer. Habitat proportions were defined as the number of grid squares of each category divided by the total number of grid squares within the quadrat. The proportion of potential cover immediately available to the lizard to escape predation was quantified as and the proportion of the quadrat which included tall grasses, wood, and other areas within the quadrat in which the individual could hide quickly. Microhabitat use could be analysed in 283 (185 females and 98 males) out of the 394 individuals.
In addition, a total of 240 random habitat images were taken throughout the study area for substrate comparison to that of basking sites. The study area can be divided into six geographical areas that are separated by dividing features, such as roads, rivers or forests [50] and random sites were determined in each of these as follows: Starting from the centre of each of the six sampling sites, a random direction was chosen by throwing a stick in the air. The stick had a designated tip, and after landing the tip of the stick determined the direction to walk. The random sampling point was reached after 20 steps, generating a nearby but separate microhabitat measure. The quadrat was then placed touching and perpendicular to the researcher’s feet and an image was taken in the same way as for the basking sites. The next random sampling point was determined by repeating the previous step, now starting from the point from which the current random sample was taken. This whole process was repeated for a total of ten random sampling points per each of the six sampling sites. Each member of the field team (N = 4) performed this step once, so that each individual sampling site contained 40 images (4 people x 10 random spots) of random habitat measures. This resulted in a total of 240 image across the whole study area (6 sampling sites x 40 images). One image was subsequently excluded due to overexposure. Habitat categories were quantified from images in the same way as for the microhabitat in which lizards were found basking.
Statistical analyses
To assess microhabitat selection in common lizards, data on the presence, absence, or proportion of various habitat types within the proximity of lizards were modelled using Bayesian zero-and-one-inflated beta (ZOIB) regression [88]. The ZOIB distribution allows modelling of proportional data that contains both zeros and ones. In such cases the response variable, \(y\), is modelled as a mixture of Bernouilli and beta distributions, from which the true zeros and ones, and the values between zero and one are generated, respectively. The probability density function is
where and \(0<\alpha ,\gamma ,\mu <1\) and \(\phi >1\). The mixture parameter, \(\alpha\), determines the extent to which the Bernoulli or beta components dominate, with higher values associated indicating greater numbers of zero or one. The parameter \(\gamma\) is the probability that \(y=1\) when \(y\) is drawn from the Bernoulli component. For a response variable where either no zeros or ones are observed \(\gamma\) is fixed as either \(\gamma =1\) or \(\gamma =0\), respectively. Habitat categories without observations \(y=1\) occurred in our data, and in these cases, the mixture parameter, \(\alpha\), determines the probability that \(y\) was zero rather than being a proportion modelled using the beta component. Continuous datapoints, \(0<y<1\), were modelled using a beta distribution parameterised with mean \(\mu\) and precision \(\phi\).
For each habitat type, nine versions of the ZOIB regression model were assessed. The first null model contained single set of parameters \(\alpha ,\gamma ,\mu\) and \(\phi\) fitted to all lizard and random sites. The next simplest model allowed \(\alpha ,\gamma ,\mu\) to vary between lizard and random sites (‘simple lizard model’). Further models allowed \(\alpha ,\gamma ,\mu\) to also vary between female and male lizards, between oviparous and viviparous, and finally by both sex and parity mode. Each model other than the null was implemented with a single \(\phi\) parameter shared by all observations and alternatively with two \(\phi\) parameters for lizard and random sites. Model comparison was done through comparison of deviance information criterion (DIC) [89], with the best model identified as that with the lowest DIC, unless a simpler model was found with a DIC score differing by less than two, in which case the simpler model was preferred.
Bayesian models were implemented with minimally informative priors for each observation class, where class may represent either random or common lizard sites in simpler models or may represent sites with common lizards of a particular sex or parity mode in more complex models. For each observation class, \(j\), prior distributions were: for the mixture parameter \(\text{l}\text{o}\text{g}\text{i}\text{t}\left({\alpha }_{j}\right)={\alpha `}_{j} \text{w}\text{h}\text{e}\text{r}\text{e}\,{\alpha `}_{j}\sim\text{N}\text{o}\text{r}\text{m}\text{a}\text{l}\left(0, 1000\right)\); for the probability in the logistic regression probability \(\text{l}\text{o}\text{g}\text{i}\text{t}\left({\gamma }_{j}\right)={\gamma `}_{j} \text{w}\text{h}\text{e}\text{r}\text{e}\,{\gamma `}_{j}\sim\text{N}\text{o}\text{r}\text{m}\text{a}\text{l}\left(0, 1\right)\); for the beta regression mean \(\text{l}\text{o}\text{g}\text{i}\text{t}\left({\mu }_{j}\right)={\mu `}_{j} \text{w}\text{h}\text{e}\text{r}\text{e}\,{\mu `}_{j}\sim\text{N}\text{o}\text{r}\text{m}\text{a}\text{l}\left(0, 1000\right)\); and beta regression precision \({\phi }_{j}={e}^{{\phi `}_{j}} \text{w}\text{h}\text{e}\text{r}\text{e}\,{\phi `}_{j}\sim\text{N}\text{o}\text{r}\text{m}\text{a}\text{l}\left(0, 100\right)\). Models were specified in the JAGS language and run using JAGS v4.3.0 [90] initiated from R (R Core Team 2013) using the package runjags, with DIC calculated using the function extract.jags [91].
Next, we assessed if common lizard dorsal colouration was dependent on any attributes of the individual (e.g., parity mode, sex, length, weight) and if it matched their background colouration. We extracted HSL (hue, saturation, lightness) values from the images of dorsal colouration, basking site, and habitat. Next, a generalised linear model (GLM) was used to identify if average dorsal colouration (measured as HSL) correlated with parity mode, sex, body length, and weight of an individual (model 1). Similarly, a GLM with a binomial distribution was used to establish if dorsal patterning (reticulated, intermediate or linear) depended on parity mode, sex, body length and weight (model 2). Using an ANOVA, we tested if dorsal colouration showed a correlation with dorsal patterning (model 3). Next, we checked if colouration of the whole microhabitat correlated with the basking site colouration (model 4). We then tested if average dorsal colouration was associated with basking site and sex (model 5), and with the colouration of the whole microhabitat and sex (model 6). Consequently, we also tested the same for dorsal patterning (model 7 and 8).
Availability of data and materials
Data are available from the Univ. of Glasgow Enlighten Repository https://doi.org/10.5525/gla.researchdata.1417. Bayesian models specified in JAGS syntax and R code required to run them on an example dataset are available at https://github.com/will-harvey/JAGSeco-evo.
References
Downes S, Shine R. Heat, safety or solitude? Using habitat selection experiments to identify a lizard’s priorities. Anim Behav. 1998;55:1387–96.
Ahnesjö J, Forsman A. Differential habitat selection by pygmy grasshopper color morphs; interactive effects of temperature and predator avoidance. Evol Ecol. 2006;20:235–57.
Gallagher AJ, Creel S, Wilson RP, Cooke SJ. Energy landscapes and the landscape of fear. Trends Ecol Evol. 2017;32:88–96.
Downes S. Trading heat and food for safety: costs of predator avoidance in a lizard. Ecol. 2001;82:2870–81.
Turner JRG. Some early works on heliconiine butterflies and their biology (Lepidoptera, Nymphalidae). J Linn Soc London Zool. 1967;46:255–66.
Linnen CR, Poh Y-P, Peterson BK, Barrett RDH, Larson JG, Jensen JD, et al. Adaptive evolution of multiple traits through multiple mutations at a single gene. Sci. 2013;339:1312–6.
Boratyński Z, Brito JC, Campos JC, Karala M, Mappes T. Large spatial scale of the phenotype-environment color matching in two cryptic species of african desert jerboas (Dipodidae: Jaculus). PLoS ONE. 2014;9(4):e94342.
Stevens M, Merilaita S. Animal camouflage: current issues and new perspectives. Philosophical Trans Royal Soc B: Biol Sci. 2009;364:423–7.
Kaufman DW. Adaptive coloration in Peromyscus polionotus: experimental selection by owls. J Mammal. 1974;55:271–83.
Merilaita S, Tuomi J, Jormalainen V. Optimization of cryptic coloration in heterogeneous habitats. Biol J Linn Soc. 1999;67:151–61.
Grant B, Howlett RJ. Background selection by the peppered moth (Biston betularia Linn.): individual differences. Biol J Linn Soc. 1988;33:217–32.
Sargent TD. Background selections of geometrid and noctuid moths. Sci. 1966;154:1674–5.
Kang CK, Moon JY, Lee SI, Jablonski PG. Camouflage through an active choice of a resting spot and body orientation in moths. J Evol Biol. 2012;25:1695–702.
Kang CK, Moon JY, Jablonski PG. Moths on tree trunks seek out more cryptic positions when their current crypticity is low. Anim Behav. 2013;86:587–94.
Marshall KLA, Philpot KE, Stevens M. Microhabitat choice in island lizards enhances camouflage against avian predators. Sci Rep. 2016;6:19815.
Rosenblum EB, Hoekstra HE, Nachman MW. Adaptive reptile color variation and the evolution of the Mc1r gene. Evol. 2004;58:1794–808.
Norris KS. Color adaptation in desert reptiles and its thermal relationships. In: Milstead WW, editor. Lizard ecology: a symposium. Columbia, MO: University of Missouri Press; 1967:162–229.
Kettlewell B. The evolution of melanism: the study of a recurring necessity with special reference to industrial melanism in the Lepidoptera. Oxford: Clarendon Press; 1973.
Porter WP, Gates DM. Thermodynamic equilibria of animals with environment. Ecol Monogr. 1969;39:227–44.
Brodie ED. Correlational selection for color pattern and antipredator behavior in the garter snake Thamnophis ordinoides. Evol. 1992;46:1284–98.
Sherbrooke WC, de Castrucci L, Hadley AM. Temperature effects on in vitro skin darkening in the mountain spiny lizard, Sceloporus jarrovi: a thermoregulatory adaptation?. Physiol Zool. 1994;67:659–72.
Kwiatkowski MA, Sullivan BK. Geographic variation in sexual selection among populations of an iguanid lizard, Sauromalus obesus (= ater). Evolution. 2002;56:2039–51.
Lepetz V, Massot M, Chaine AS, Clobert J. Climate warming and the evolution of morphotypes in a reptile. Glob Change Biol. 2009;15:454–66.
Sacchi R, Pellitteri-Rosa D, Bellati A, Scali S, Tronconi D, Galeotti P, et al. Biological correlates of two dorsal color pattern types in the common wall lizard, Podarcis muralis (Laurenti, 1768) (Squamata: Sauria: Lacertidae). Herpetozoa. 2012;25:3–11.
Rajabizadeh M, Adriaens D, Kaboli M, Sarafraz J, Ahmadi M. Dorsal colour pattern variation in eurasian mountain vipers (genus Montivipera): a trade-off between thermoregulation and crypsis. Zool Anz-A J Comp Zool. 2015;257:1–9.
Shine R. Costs of reproduction in reptiles. Oecologia. 1980;46:92–100.
Bauwens D, Thoen C. Escape tactics and vulnerability to predation associated with reproduction in the lizard Lacerta vivipara. J Anim Ecol. 1981;50:733–43.
Lima SL, Dill LM. Behavioral decisions made under the risk of predation: a review and prospectus. Can J Zool. 1990;68:619–40.
Zajitschek SR, Zajitschek F, Clobert J. The importance of habitat resistance for movement decisions in the common lizard, Lacerta vivipara. BMC Ecol. 2012;12:13.
Lorenzon P, Clobert J, Massot M. The contribution of phenotypic plasticity to adaptation in Lacerta vivipara. Evol. 2001;55:392–404.
Main KL. Predator avoidance in seagrass meadows: prey behavior, microhabitat selection, and cryptic coloration. Ecol. 1987;68:170–80.
Orrock JL. Rodent foraging is affected by indirect, but not by direct, cues of predation risk. Behav Ecol. 2004;15:433–7.
Suselbeek L, Emsens WJ, Hirsch BT, Kays R, Rowcliffe JM, Zamora-Gutierrez V, et al. Food acquisition and predator avoidance in a neotropical rodent. Anim Behav. 2014;88:41–8.
Avery RA, Bedford JD, Newcombe CP. The role of thermoregulation in lizard biology: predatory efficiency in a temperate diurnal basker. Behav Ecol Sociobiol. 1982;11:261–7.
Martín J, López P. Escape behaviour of juvenile Psammodromus algirus lizards: constraint of or compensation for limitations in body size? Behaviour. 1995;132:181–92.
Row JR, Blouin-Demers G. Thermal quality influences effectiveness of thermoregulation, habitat use, and behaviour in milk snakes. Oecologia. 2006;148:1–11.
Clifford BJ, Bauder JM, Marchand MN. Tree removal likely improves thermal quality of basking sites for an imperiled timber rattlesnake population at the northern edge of its range. Northeastern Naturalist. 2020;27:746–62.
Vanhooydonck B, Van Damme R. Relationships between locomotor performance, microhabitat use and antipredator behaviour in lacertid lizards. Funct Ecol. 2003;17:160–9.
Brodie ED. Behavioral modification as a means of reducing the cost of reproduction. Am Nat. 1989;134:225–38.
Shine R. Locomotor speeds of gravid lizards: placing “costs of reproduction” within an ecological context. Funct Ecol. 2003;7:526–33.
Quirt KC, Blouin-Demers G, Howes BJ, Lougheed SC. Microhabitat selection of five-lined skinks in northern peripheral populations. J Herpetology. 2006;40:335–42.
Asbury DA, Adolph SC. Behavioural plasticity in an ecological generalist: Microhabitat use by western fence lizards. Evol Ecol Res. 2007;9:801–15.
Gvoždík L. To heat or to save time? Thermoregulation in the lizard Zootoca vivipara (Squamata: Lacertidae) in different thermal environments along an altitudinal gradient. Can J Zool. 2002;80:479–92.
Van Damme R, Bauwens D, Verheyen RF. Thermoregulatory responses to environmental seasonality by the lizard Lacerta vivipara. Herpetologica. 1987;43:405–15.
Martin J, López P. Influence of habitat structure on the escape tactics of the lizard Psammodromus algirus. Can J Zool. 1995;73:129–32.
Stevens M, Troscianko J, Wilson-Aggarwal JK, Spottiswoode CN. Improvement of individual camouflage through background choice in ground-nesting birds. Nat Ecol Evol. 2017;1:1325–33.
Wilson-Aggarwal JK, Troscianko JT, Stevens M, Spottiswoode CN. Escape distance in ground-nesting birds differs with individual level of camouflage. Am Nat. 2016;188:231–9.
Recknagel H, Kamenos NA, Elmer KR. Common lizards break Dollo’s law of irreversibility: genome-wide phylogenomics support a single origin of viviparity and re-evolution of oviparity. Mol Phylogenet Evol. 2018;127:579–88.
Recknagel H, Carruthers M, Yurchenko AA, Nokhbatolfoghahai M, Kamenos NA, Bain MM, et al. The functional genetic architecture of egg-laying and live-bearing reproduction in common lizards. Nat Ecol Evol. 2021;5:1546–56.
Recknagel H, Elmer KR. Differential reproductive investment in co-occurring oviparous and viviparous common lizards (Zootoca vivipara) and implications for life-history trade-offs with viviparity. Oecologia. 2019;190:85–98.
Arribas OJ. Morphological variability of the Cantabro-Pyrenean populations of Zootoca vivipara with description of a new subspecies. Herpetozoa. 2009;21:123–46.
Stuart-Fox D, Whiting MJ, Moussalli A. Camouflage and colour change: antipredator responses to bird and snake predators across multiple populations in a dwarf chameleon. Biol J Linn Soc. 2006;88:437–46.
Fitze PS, Cote J, Clobert J. Mating order-dependent female mate choice in the polygynandrous common lizard Lacerta vivipara. Oecologia. 2010;162:331–41.
Laloi D, Eizaguirre C, Fédérici P, Massot M. Female choice for heterozygous mates changes along successive matings in a lizard. Behav Process. 2011;88:149–54.
Le Galliard JF, Rozen-Rechels D, Lecomte A, Demay C, Dupoue A, Meylan S. Short-term changes in air humidity and water availability weakly constrain thermoregulation in a dry-skinned ectotherm. PLoS ONE. 2021;16(2):e0247514.
Grenot CJ, Garcin L, Dao J, Hérold J-P, Fahys B, Tséré-Pagès H. How does the european common lizard, Lacerta vivipara, survive the cold of winter? Comp Biochem Physiol Part A Mol Integr Physiol. 2000;127:71–80.
Halley A. Choice of substrate and heating rate in Lacerta vivipara. Br J Herpetology. 1982;6:207–13.
Le Galliard J-F, Le Bris M, Clobert J. Timing of locomotor impairment and shift in thermal preferences during gravidity in a viviparous lizard. Funct Ecol. 2003;17:877–85.
Van Damme R, Bauwens D, Verheyen RF, Veheyen RF. Effect of relative clutch mass on sprint speed in the lizard Lacerta vivipara. J Herpetology. 1989;23:459–61.
Cornetti L, Ficetola GF, Hoban S, Vernesi C. Genetic and ecological data reveal species boundaries between viviparous and oviparous lizard lineages. Heredity. 2015;115:517–26.
Stuart-Fox DM, Moussalli A, Johnston GR, Owens IP. Evolution of color variation in dragon lizards: quantitative tests of the role of crypsis and local adaptation. Evolution. 2004;58:1549–59.
Reguera S, Zamora-Camacho FJ, Moreno-Rueda G. The lizard Psammodromus algirus (Squamata: Lacertidae) is darker at high altitudes. Biol J Linn Soc. 2014;112:132–41.
Clusella Trullas S, van Wyk JH, Spotila JR. Thermal melanism in ectotherms. J Therm Biol. 2007;32:235–45.
Stevens M, Ruxton GD. The key role of behaviour in animal camouflage. Biol Rev. 2019;94:116–34.
Andrén C, Nilson G. Reproductive success and risk of predation in normal and melanistic colour morphs of the adder, Vipera berus. Biol J Linn Soc. 1981;15:235–46.
Forsman A. Heating rates and body temperature variation in melanistic and zigzag Vipera berus: does colour make a difference? Ann Zool Fenn. 1995;32:365–74.
Endler JA. A predator’s view of animal color patterns. Evolutionary Biology. Boston, MA: Springer US; 1978. 319–64.
Lythgoe JN. The ecology of vision. Oxford: Clarendon Press; 1979.
Cott HB. Adaptive coloration in animals. London: Methuen; 1940.
Merilaita S. Crypsis through disruptive coloration in an isopod. Proc R Soc B: Biol Sci. 1998;265:1059–64.
Steward RC. Industrial and non-industrial melanism in the peppered moth, Biston betularia (L). Ecol Entomol. 1977;2:231–43.
Götmark F, Hohlfält A. Bright male plumage and predation risk in passerine birds: are males easier to detect than females? Oikos. 1995;74:475–84.
Götmark F. Conspicuous coloration in male birds is favoured by predation in some species and disfavoured in others. Proc R Soc B: Biol Sci. 1993;253:143–6.
Fleishman LJ, Loew ER, Leal M. Ultraviolet vision in lizards. Nature. 1993;365:397–7.
Cuthill IC, Partridge JC, Bennett ATD, Church SC, Hart NS, Hunt S. Ultraviolet vision in birds. Adv Study Behav. 2000;29:159–214.
Loew ER, Fleishman LJ, Foster RG, Provencio I. Visual pigments and oil droplets in diurnal lizards. J Exp Biol. 2002;205:927–38.
Cronin TW, Caldwell RL, Marshall J. Tunable colour vision in a mantis shrimp. Nature. 2001;411:547–8.
Endler JA. Evolutionary implications of the interaction between animal signals and the environment. In: Espmark Y, Amundsen T, Rosenqvist G, editors. Animal signals: signalling and signal design in animal communication. Trondheim, Norway: Tapir Academic Press; 2000:11–4.
Endler JA, Basolo AL. Sensory ecology, receiver biases and sexual selection. Trends Ecol Evol. 1998;13:415–20.
Maia R, White TE. Comparing colors using visual models. Behav Ecol. 2018;29:649–59.
Catchen J, Hohenlohe PA, Bassham S, Amores A, Cresko WA. Stacks: an analysis tool set for population genomics. Mol Ecol. 2013;22:3124–40.
Yurchenko AA, Recknagel H, Elmer KR. Chromosome-level assembly of the common lizard (Zootoca vivipara) genome. Genome Biol Evol. 2020;12:1953–60.
Stevens M, Párraga CA, Cuthill IC, Partridge JC, Troscianko TS. Using digital photography to study animal coloration. Biol J Linn Soc. 2007;90:211–37.
Bergman TJ, Beehner JC. A simple method for measuring colour in wild animals: validation and use on chest patch colour in geladas (Theropithecus gelada). Biol J Linn Soc. 2008;94:231–40.
Burgon JD, Vieites DR, Jacobs A, Weidt SK, Gunter HM, Steinfartz S, et al. Functional colour genes and signals of selection in colour-polymorphic salamanders. Mol Ecol. 2020;29:1284–99.
Martin M, Le Galliard J-F, Meylan S, Loew ER. The importance of ultraviolet and near-infrared sensitivity for visual discrimination in two species of lacertid lizards. J Exp Biol. 2015;218:468–465.
Hart NS, Hunt DM. Avian visual pigments: characteristics, spectral tuning, and evolution. Am Nat. 2007;169:7–26.
Ospina R, Ferrari SLP. Inflated beta distributions. Stat Pap. 2010;51:111–26.
Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A. Bayesian measures of model complexity and fit. J Royal Stat Society: Ser B (Statistical Methodology). 2002;64:583–639.
Plummer M. Just Another Gibbs Sampler v4.3.0 (JAGS): A program for analysis of Bayesian graphical models using Gibbs sampling. 2017. https://sourceforge.net/projects/mcmc-jags/. Accessed 8 Feb 2022.
Denwood MJ. Runjags: an R Package Providing Interface Utilities, Model Templates, parallel Computing methods and additional distributions for MCMC Models in JAGS. J Stat Softw. 2016;71:1–25.
Acknowledgements
We thank H. Leitão and M. Sutherland for assistance during fieldwork. We are grateful to Alpenferienpark Reisach for accommodation and support throughout fieldwork. We thank three anonymous referees for review comments that improved the manuscript and M. Uzzell for constructive proofreading.
Ethical guidelines (Accordance statement)
All research was done compliant with permissions for studying common lizards granted by Bezirkshauptmannschaft Hermagor (HE3-NS-959/2013).
Funding
This research was funded by NERC (NE/N003942/1) and the Univ. of Glasgow (Masters and LKAS PhD funding).
Author information
Authors and Affiliations
Contributions
KRE and HR designed the study. HR and ML collected the data. HR, WH, and ML analysed the data and generated statistical models. HR, ML and WH led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
As this study is observational of animals in their natural environment, it does not require ethical permissions for animal experimentation from the UK nor Austria. All research was done compliant with permissions for studying common lizards granted by Bezirkshauptmannschaft Hermagor (HE3-NS-959/2013).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Parameters in zero-and-one-inflated beta regression models favoured by DIC.
Additional file 2: Table S2.
Statistical models performed on lizard colouration and patterning. The effect of any significantly associated variable in a model is shown.
Additional file 3: Figure S1.
Dorsal colouration differs between parity modes and correlates with body weight. Viviparous (A) and larger (B) individuals tended to have larger hue values compared to oviparous and smaller individuals.
Additional file 4: Figure S2.
Histograms of A) the wider habitat lightness and B) basking site lightness.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Recknagel, H., Harvey, W.T., Layton, M. et al. Common lizard microhabitat selection varies by sex, parity mode, and colouration. BMC Ecol Evo 23, 47 (2023). https://doi.org/10.1186/s12862-023-02158-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12862-023-02158-2