*Article* **Inter-Individual Differences in Ornamental Colouration in a Mediterranean Lizard in Relation to Altitude, Season, Sex, Age, and Body Traits**

**Gregorio Moreno-Rueda 1, Senda Reguera 1, Francisco J. Zamora-Camacho <sup>2</sup> and Mar Comas 1,\***


**\*** Correspondence: marcomas@ugr.es

**Abstract:** Animals frequently show complex colour patterns involved in social communication, which attracts great interest in evolutionary and behavioural ecology. Most researchers interpret that each colour in animals with multiple patches may either signal a different bearer's trait or redundantly convey the same information. Colour signals, moreover, may vary geographically and according to bearer qualities. In this study, we analyse different sources of colour variation in the eastern clade of the lizard *Psammodromus algirus*. Sexual dichromatism markedly differs between clades; both possess lateral blue eyespots, but whereas males in the western populations display strikingly colourful orange-red throats during the breeding season, eastern lizards only show some commissure pigmentation and light yellow throats. We analyse how different colour traits (commissure and throat colouration, and the number of blue eyespots) vary according to body size, head size (an indicator of fighting ability), and sex along an elevational gradient. Our findings show that blue eyespots function independently from colour patches in the commissure and throat, which were interrelated. Males had more eyespots and orange commissures (which were yellow or colourless in females). Throat colour saturation and the presence of coloured commissures increased in older lizards. The number of eyespots, presence of a coloured commissure, and throat colour saturation positively related to head size. However, while the number of eyespots was maximal at lowlands, throat colour saturation increased with altitude. Overall, our results suggest that this lizard harbours several colour signals, which altitudinally differ in their importance, but generally provide redundant information. The relevance of each signal may depend on the context. For example, all signals indicate head size, but commissure colouration may work well at a short distance and when the lizard opens the mouth, while both throat and eyespots might work better at long distance. Meanwhile, throat colouration and eyespots probably work better in different light conditions, which might explain the altitudinal variation in the relative importance of each colour component.

**Keywords:** colouration; social signals; *Psammodromus algirus*; lizards; altitudinal gradient

#### **1. Introduction**

Colour ornaments are frequent in the animal kingdom, typically involved in social communication [1]. However, understanding the evolution of colour patterns in animals is challenging because it often arises from the interaction of concomitant selective pressures. Sexual selection favours sexual dichromatism and colourful patterns [2], whereas natural selection selects for dull and cryptic colourations [3]. Besides, animals frequently show complex and contrasting colour patterns consisting of multiple colour patches [4]. Multiple ornamental colour patches within individuals may be the result of different selective pressures on each patch [5,6], and so each component of colouration may be related to different individual traits (multiple message hypothesis; e.g., [7–9]). Alternatively, different

**Citation:** Moreno-Rueda, G.; Reguera, S.; Zamora-Camacho, F.J.; Comas, M. Inter-Individual Differences in Ornamental Colouration in a Mediterranean Lizard in Relation to Altitude, Season, Sex, Age, and Body Traits. *Diversity* **2021**, *13*, 158. https://doi.org/ 10.3390/d13040158

Academic Editor: Luc Legal

Received: 8 March 2021 Accepted: 5 April 2021 Published: 6 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

colour patches may act as redundant signals providing similar information, in this way increasing the reliability of the signal (backup hypothesis, [10,11]).

In reptiles, colouration is the result of a mix of pigments (carotenoids, pterins, and melanin) and structural layers including crystalline platelets [12]. In lizards, the size and spectral characteristics of colour patches have been related to different traits [13]. For example, different characteristics of colour patches can reflect fighting ability [14–16], as well as reproductive status in both females [17–19] and males [20–22], among several other traits.

Understanding the processes that generate intraspecific phenotypic variation is essential in evolutionary ecology. Colouration may vary geographically if selective pressures also vary [23–25]. In this sense, elevational gradients offer a valuable study framework encompassing a considerable environmental variation in a relatively short spatial range [26]. For example, Badyaev [27] described an elevational pattern in birds, in which the strength of sexual selection on colouration decreased with altitude. Badyaev [27] proposed several explanations for this pattern, some of which apply only to species with prolonged parental care. However, it is unknown whether such a pattern is also applicable to reptiles, studies on lizards providing mixed results [28–32].

In the present study, we investigated several sources of inter-individual variation in social colouration of the Mediterranean lizard *Psammodromus algirus* in Sierra Nevada Mountain (Spain). *Psammodromus algirus* is a medium-large lacertid (53–80 mm snoutvent length, SVL, in our study area) that inhabits shrubby habitats in the Mediterranean region of western North Africa, the Iberian Peninsula, and southern France [33]. In the Iberian Peninsula, this species is split into two phylogeographic clades [34], which differ in male colouration during the breeding season (see Figures 1 and 2 in [33]). In the western clade, adult males typically show orange-red colour on the head and throat during the breeding season [19,30,35–37], while young males only exhibit an orange spot in their mouth commissures [38]. These colour patches are absent in females. Meanwhile, social colouration in the eastern clade has been much less studied. In the eastern clade, adult males present an orange spot in their mouth commissures, but not orange colouration in the head and throat [39], just like young males in the western clade. Both adult males and females may exhibit a yellow patch on their throats during the breeding season [39]. Besides, *P. algirus* from both clades display a variable number of blue-ultraviolet eyespots in their flanks, which are more numerous in males than in females, and show little seasonal variation [39,40].

The aim of this study is to examine the sources of inter-individual variation in the colour variables involved in social communication of a southern population within the eastern clade of *P. algirus*. We were interested in the diversity of colour patches in this species and its marked geographical variation. We related several colour traits (throat lightness, chroma, and hue, commissure colour and patch size, and the number of flank eyespots) with:


Our final goal is to add to the knowledge of the evolution of lizard colouration. Concretely, our main purposes are to understand whether such a variety of signals provide different information on lizard quality or well provide redundant information, as well as to give insights on the sources of geographic variation in social colouration.

**Figure 1.** The three-dimensional map (bottom panel) is a representation of the Sierra Nevada Mountain and displays the location of the six sites sampled during this study along the elevational gradient: 300 (1), 700 (2), 1200 (3), 1700 (4), 2200 (5), and 2500 (6) m asl. The location of Sierra Nevada in the Iberian Peninsula (**top**, **left**) and an image of the lizard (**top**, **right**) are also shown.

**Figure 2.** Representation of procedures for quantifying the coloured commissure area, using the software Image J (version 1.60). (**A**) First, we performed lateral photographs of the lizard with the mouth opens; (**B**) Second, we selected the coloured area with the "Colour threshold" tool; (**C**) Third, we measured the area with the "Analyse particles" tool.

#### **2. Material and Methods**

#### *2.1. Sampling and General Procedures*

Fieldwork was performed in the Sierra Nevada mountain system (SE Spain), where *P. algirus* inhabits from 200 to 2800 m asl, with permission of the Andalusian government and National Park of Sierra Nevada (references GMN/GyB/JMIF and ENSN/JSG/JEGT/MCF). We sampled six localities sited at 300, 700, 1200, 1700, 2200, and 2500 m asl (Figure 1) with a similar structure of vegetation (more details in [44]). Lizards were captured by hand and transported to the lab in cotton bags. Sampling occurred during their activity season in Sierra Nevada, spanning from March to September [44], during the years 2010–2013. We captured 482 adult lizards (males/females per year, 2010: 43|58, 2011: 49|50, 2012: 65|60, 2013: 82|75). We measured their SVL with a metal ruler (accuracy 1 mm), body mass with a digital scale (accuracy 0.01 g) and head width with a digital calliper (accuracy 0.01 mm). Sex was determined according to femoral-pore development (more developed in males [29]). Only adults were considered because we were interested in inter- and intra-sexual communication. Given that SVL did not differ between sexes in our study area [45], we considered adults lizards with the minimum body size for which we found gravid females. Notice that body size varies with altitude [45], so this minimum SVL was estimated for each altitude (at 300 m: 53 mm; 700–1700 m: 55 mm; 2200 m: 62 mm; 2500 m: 63 mm). Because lizards were part of a long-term study, they were marked by toe clipping and resampled lizards (~5%) were excluded to avoid pseudoreplication. Toe clipping is frequently used to mark lizards with little impact on welfare and survival [46]. Toes of a subsample (*n* = 118) were conserved in ethanol and used for age determination using phalanx skeletochronology (detailed methods in [47]). Skeletochronology was a technique widely used to estimate age in reptiles (at the accuracy of a year), which has proven to be very accurate [48].

#### *2.2. Quantification of Colour Patches*

Along the four years of the study, we measured throat (gular region) colour with a colourimeter (Minolta CM-2600d). The colourimeter, placed on the lizard's skin, projected three beams of light through a 3-mm-diameter hole. As a result, it took three measures of reflectance and automatically gave the average of each colour component [49]. These components correspond to the L\*a\*b\* colour-space of the Commission Internationale d'Eclairage (CIE 1976 [50]). The device measured the reflectance of the surface on the spectrum of light of the visible range, from 400 to 700 nm wavelengths. We did not study skin reflectance in the ultra-violet (UV) range of the spectrum because the throat has no reflectance peak in the UV range. L\*a\*b\* colour space is a 3-dimensional rectangular colour space. L\* axis represents lightness (0 is black, 100 is white); a\* axis represents red-green gradient (positive values are red, negative values are green); b\* axis represents blue-yellow gradient (positive values are yellow, negative values are blue). From L\*a\*b\* values we determined *chroma* (saturation or purity) as C\* = [(a\*)<sup>2</sup> + (b\*)2] 1/2 (measured as the percentage distance from the centre [0] of the colour space to its circumference [100] where pure spectral colours are represented); and hue angle (the "colour" in common parlance) as H\* = tan<sup>−</sup>1(b\*/a\*) [51].

In 2011–2013, we additionally recorded the presence or absence, and colour (orange or yellow) if present, of a patch in the mouth commissures, and counted the total number of blue eyespots in the lizard flanks. For those lizards presenting a well-differentiated colour patch in commissures, we took photographs from the right side of the head (consistently with the mouth open; Figure 2). We used a Canon Power Shot SX200 IS digital camera and a graph paper background for size reference. Then, we measured the area of this pigmented patch using the software Image J (version 1.60 [52]). First, we scaled the photos using the graph paper (in mm) and the "Set scale" tool. Then, we adjusted the area of the patch using the "Colour threshold" tool. We afterwards measured the area of each coloured patch with the "Analyse particles" tool (Figure 2). The software measured all the patches delineated and gave the total selected area in mm2.

#### *2.3. Statistical Analyses*

Sample sizes differed among the variables measured (Appendix A). Eyespot number and commissure colour and size were not recorded in 2010. The size of the commissure patch was measured only in lizards in which it was present (50% of lizards). Age was estimated in a random subsample of 118 lizards. For diverse reasons, some data for different variables were lost in several individuals. The fact that not all variables were available in every individual conditioned the statistical analyses we could perform minimising the loss of sample size.

Firstly, we checked for possible outliers in every variable by using Cleveland plots [53]. Secondly, we graphically checked the normality and homoscedasticity of the variables [53]. Throat hue and area of the commissure were log-transformed to match these assumptions. Throat lightness was arcsin-transformed. SVL, body mass and head width were logtransformed to meet homoscedasticity and linearity [54]. Throat chroma and the number of eyespots were not transformed. Although all analyses were carried out with the variables transformed, raw data are shown in graphics and when providing mean values (with SE).

The colour variables measured (number of eyespots, throat lightness, throat chroma, throat hue, presence and colour of commissure patch, and area of the commissure patch) may not be independent among themselves. For this reason, in a first analysis, we examined the bivariate Pearson product-moment correlation between continuous variables, and we tested with Anova whether the value of the colour variables differed among categories of commissure colour (colourless, orange, yellow).

In a second analysis, we carried out a linear model for every colour variable with a continuous distribution, and a multinomial model linked to a logit function for commissure colour. In these models, the predictor variables were SVL, body mass, head width (all log-transformed, continuous), altitude (six levels: 300, 700, 1200, 1700, 2200, 2500 m), season (date when lizards were captured, arranged into three categories: March–May, June–July, August–September), sex (males, females), and year (2010, 2011, 2012, 2013; data for 2010 only available for throat lightness, chroma, and hue). All analyses were carried out with

R 3.6.1 [55]. Multinomial models were performed using the function "multinom" of the package "nnet" v. 7.3–14 [56]. In addition, we applied a model selection approach based on the Akaike Information Criterion (AIC) with the "MuMIn" package [57]. In this analysis, models with a value of ΔAIC less than 2 were chosen [58]. Moreover, given that some colour variables showed some correlation (see below), we repeated the best models including the other variables of colouration as covariates to check the soundness of the findings. Normality and homoscedasticity of every model's residuals were checked [53]. Only significant effects supported by the three statistical approaches (full model, model selection, and best models controlling for other colour variables) were considered as sound. We also checked for interactions between the predictor variables. However, most interactions were non-significant or lacked soundness, so we decided not to show them for the sake of clarity.

Lastly, lacertid lizards typically have indeterminate growth. Consequently, older lizards are also larger. In this way, age and body size may be confounded. To test the effect of age, we repeated the analyses with a subsample of known-age individuals, by including the age estimated with skeletochronology in the models. Given that these models included SVL and age, the effect of the two variables can be disentangled.

#### **3. Results**

#### *3.1. Relationship among Colour Parameters*

The number of blue eyespots was unrelated to the other colouration variables (throat chroma, lightness and hue, and area of the commissure), but lizards with orange commissures had significantly more eyespots (*F*2, 271 = 31.54, *p* < 0.001; Table 1). This finding seems to be linked to the fact that males have more eyespots and a higher probability to present orange commissures (see below). Indeed, when sex was introduced as a predictor, the relationship between number of eyespots and colour of commissures disappeared (*F*2, 270 = 1.49, *p* = 0.23; effect of sex: *F*1, 270 = 39.15, *p* < 0.001). Meanwhile, throat colouration parameters were interrelated; chroma and lightness were negatively correlated, while hue was positively correlated with lightness and so negatively related to chroma (Table 1). Nonetheless, correlation coefficients were relatively low (|*r*| ≤ 0.51). The size of commissure patch was positively correlated with throat chroma, and negatively with hue and lightness (Table 1). In addition, individuals with a commissure patch (orange or yellow) had more saturated throats (*F*2, 274 = 23.31, *p* < 0.001; Table 1). The reverse occurred for throat lightness (*F*2, 274 = 17.96, *p* < 0.001; Table 1). Throat hue was unrelated to commissure colour (*F*2, 274 = 0.24, *p* = 0.79). Lizards with orange commissures tended to have larger commissure patches than lizards with yellow commissures (0.74 ± 0.07 vs. 0.48 ± 0.16 mm2; *<sup>F</sup>*1, 111 = 3.14, *<sup>p</sup>* = 0.079). Therefore, mouth commissure and throat colourations were seemingly interrelated.

#### *3.2. Correlates of Colouration*

The full model (Table 2) showed that the number of eyespots was maximal at a low elevation and minimal at mid-elevation (Figure 3a). Males had more eyespots than females (Figure 3b). Moreover, lizards with larger heads showed more eyespots (Figure 3c). The number of eyespots did not significantly vary with season or body size (SVL or mass). According to a model selection approach, the best model was that including altitude, sex, and head size (Table 3). The second-best model included those variables plus SVL (ΔAIC = 1.79), which, however, failed to significantly explain variation in the number of eyespots.

Throat lightness significantly varied with altitude (Table 2). Concretely, lizards from localities at 2200 and 2500 m had darker throat than those from lower altitudes (Figure 4a). Lizard throat was darker at the end of the breeding season (Figure 4b). Lizards with larger SVL and heads showed darker throats than smaller lizards (Table 2; Figure 4c; data not shown for SVL). Besides, throat lightness varied among years (Table 2). The four best models (ΔAIC < 2) included altitude, season, year, and head size as predictors (Table 3). Given that throat lightness was correlated with throat chroma and hue, we repeated the

best model including these two variables as predictors to control for them. The resulting model was qualitatively very similar to the best model (data not shown), with the exception that SVL was no longer significantly related to throat lightness. SVL was not included in the fourth best model (Table 3), so there was little support for this variable being related to throat lightness.

Throat chroma increased with altitude (Figure 5a) and decreased with the advance of the season (Figure 5b; Table 2). Throat chroma increased with SVL and head size (Figure 5c, data not shown for head size). Besides, throat chroma showed interannual variation (Table 2). With the model selection, three models had ΔAIC < 2, all including altitude, season, year, SVL, and head size as significant predictors of throat chroma (Table 3). When controlling for throat lightness and hue, results were similar, but head size was no longer significant, and a marginally significant effect of body mass emerged (data not shown).

**Table 1.** Correlations among the continuous variables (provided sample size -in the subscript-, correlation coefficient and *p*-value are provided), and average values (±SE) for each category of commissure colour (no colour, yellow, or orange). Sample sizes for each category of commissure colour between brackets (for eyespots of lizards with orange commissures, the sample size was 103). Raw data are shown, but statistical tests were carried out with transformed variables when necessary. In bold significant relationships or differences. Different superscripts indicate significant differences according to an unequal N HSD post hoc test.


**Table 2.** Full models for each colour variable. In bold, predictors with a significant effect. Degree of freedom (df) as well as *F*-value for lineal models and *χ*<sup>2</sup> for multinomial model are shown. \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.


**Figure 3.** The average number of eyespots (with SE) in the flanks of the lizard *P. algirus*, according to altitude (**A**) and sex (**B**), and the relationship between the number of eyespots and head width (**C**). Notice that raw data are shown, but statistical analyses were performed with transformed data when necessary.



**Figure 4.** Average values of throat lightness (with SE) in the lizard *P. algirus*, according to altitude (**A**) and season (**B**), and the relationship between throat lightness and head width (**C**). Notice that raw data are shown, but statistical analyses were performed with transformed data when necessary.

**Figure 5.** Average values of throat chroma (with SE) in the lizard *P. algirus*, according to altitude (**A**) and season (**B**), and the relationship between throat chroma and SVL (**C**). Notice that raw data are shown, but statistical analyses were performed with transformed data when necessary.

According to the full model, throat hue showed covariation with altitude, season, and year, but not with sex or lizard morphology (Table 2). Throat hue tended to decrease with altitude and increased with the advance of the season (data not shown for simplicity). The model selection provided similar results, but, the best models also included sex, which significantly explained part of the variation in throat hue (Table 3). The inclusion of throat lightness and chroma in the model did not alter significantly the results (not shown).

The colour of the commissure differed between sexes, males usually having orange commissures, while females had yellow or no coloured commissure (Figure 6a; Table 2). The frequency of lizards with orange commissures decreased with the advance of the breeding season (Figure 6b; Table 2). Lastly, individuals with coloured commissure (either orange or yellow) had larger heads than individuals without commissure (Figure 6c, Table 2). No other variable was significantly related to commissure colour (Table 2). Model selection approach selected for two models including sex, season, and head size as significant predictors (Table 3).

**Figure 6.** Frequency of *P. algirus* lizards with colourless, yellow or orange commissure according to sex (**A**) and season (**B**), and head width (with SE) of lizards according to commissure colour (**C**).

The area of the coloured commissure varied significantly with season and year (Table 2), tending to decrease with the advance of the breeding season (Figure 7). Eight models were selected by model selection. All models included season and year as significant predictors. No model included altitude. The remaining variables were included in five models, being significant in some, but not in others (Table 3). A model including throat lightness, chroma, and hue did not alter significantly these results (not shown).

#### *3.3. The Effect of Age*

As expected, SVL increased with age (*F*4, 113 = 14.06, *p* < 0.001; Appendix B). We found that older individuals were more likely to have coloured commissures (*χ*<sup>2</sup> <sup>6</sup> = 12.83, *p* = 0.046; Figure 8a). Moreover, throat chroma increased with age (Figure 8b). This last result was confirmed by model selection. Two models were selected including the variables age + body mass (AICc = 678.33) and age + body mass + head size (AICc = 679.14). Notice

that the two models included age, but not SVL, as a predictor of throat chroma. Previous analyses showed that throat chroma was the only colour variable that increased with SVL (Table 4). However, the present findings suggest that throat chroma is influenced by age rather than by body size.

**Figure 7.** Average values of commissure area (with SE) in the lizard *P. algirus*, according to season.

**Figure 8.** Frequency of lizards with colourless, yellow or orange commissure according to age (**A**) and average (with SE) throat chroma according to age (**B**).

**Table 4.** Summary of the main results. In bold are the clear results, whereas somewhat dubious results are not in bold.


#### **4. Discussion**

Table 4 summarises the conclusions based on the main results. The findings in this study suggest that blue eyespots function independently from colour patches in the commissure and throat, while throat and commissure colouration were interrelated. Moreover, eyespots and both commissure and throat varied in different ways with altitude

and season. Eyespots varied with altitude following a U-shaped pattern, while throat colour increased in saturation (chroma) with altitude. Also, throat colouration was darker in the highest elevations, where lizards are typically darker [59]. While eyespots were permanent colour patches (as in the western clade [40]), commissure and throat colouration varied with the season. At the beginning of the breeding season, we found more lizards with coloured commissures and throat colour was darker and more saturated. Similar findings were reported for a northern population of the same phylogenetic clade [39], which suggests that throat and commissure patches are involved in communication during breeding. Similarly, in the western clade, orange head colouration in adult males and coloured commissures in young males are present during the breeding season [19,35,36,38].

We found some degree of sexual dichromatism; males typically had more eyespots and orange commissures (which were yellow or colourless in females). Hence, although sexual dichromatism was not as pronounced as in the western clade, a slight dichromatism is still present in the eastern clade. Throat saturation and the presence of coloured commissures was indicative of age. Meanwhile, the number of eyespots, presence of a coloured commissure, and throat saturation, all were indicators of head size, and thus of fighting ability. Therefore, although colour patches changed with altitude and season in different ways (probably reflecting different subjacent colour production costs and mechanisms and/or selective pressures), they seem, in general, to be redundant indicating the same traits of lizards: sex, fighting ability, and age. One possibility is that these signals are used in different contexts, indicating the same traits at different distances or in different light conditions.

#### *4.1. Colouration as Indicators of Fighting Ability, Sex, and Age*

All colour patches measured (eyespots, commissure and throat) seem indicators of head size (Table 4), which is known as a good indicator of fighting ability [60]. Lizards with greater fighting ability are more successful in defending their territories and hence obtain higher mating success [37]. Several studies suggest that blue-ultraviolet colouration, as that present in eyespots, is related to fighting ability [16,61,62], but pigmentary colourations may also serve as indicators of fighting ability [63]. During a contest, signals correlated with fighting ability may help individuals to assess the relative competitive ability of rivals and so to avoid being involved in a costly physical combat [2]. To be informative, such signals should be honest indicators of fighting ability. Colouration may be a reliable indicator of fighting ability when maintained by social interactions [64]. Indeed, in the western clade of *P. algirus*, orange head colouration is a good indicator of social dominance, but orange-headed lizards are also more often involved in fighting, so subordinate lizards would pay a cost if vividly coloured [41]. Moreover, throat and commissure colouration might act as amplifiers of head size. This may be especially important in commissure colouration, which is displayed only when the mouth is open, hence allegedly showing the intention of biting and amplifying the perceived mouth size [65]. Furthermore, in the western clade, several costs have been associated with male nuptial colouration, such as reduced immune capacity, increased risk of ectoparasitism and ultimately reduced survival, which could serve to maintain honesty [30,35,36].

While lizards from western populations are strongly dichromatic, sexual dichromatism in our study population was reduced to commissure colour (orange in males, yellow or colourless in females) and the number of eyespots (more numerous in males). Therefore, commissure colour seems to intervene in sex recognition. The fact that females display yellow throat and blue eyespots similarly to males suggests that females also use colour patches in social communication. Female lizards frequently display colour patches, usually related to receptiveness to breed [17–19]. However, alternative explanations, typically poorly explored, are possible. Females could use throat colouration (and yellow commissures) to indicate fighting ability to rivals. Females could also indicate some type of individual quality to potential mates [66]. However, colouration in females might simply be the result of correlated selection in males [67]. Anyway, the reduced sexual dichromatism

in the eastern clade suggests that sexual selection on male traits is weaker in the eastern than in the western clade.

In Sierra Nevada, older individuals expressed both the yellow patch on the throat and the orange commissures more frequently than young lizards. Given that this is a trans-sectional study, we cannot disentangle whether colourful lizards lived longer (duller lizards selectively disappearing from the population as they aged), or lizards invested more in colouration as they aged. In either case, these colour patterns would indicate higher survival capability of the bearer. Younger lizards may be selected for concealment, mimicking females, hence not developing the yellow patch or the orange commissure until they grow large enough to compete with older males [38,41]. However, dominant males may detect female-mimicking males using chemosensorial cues [68]. Alternatively, colour cues indicating lizard age could be used to evaluate the survival prospects of potential mates [43]. In order for signals of survival ability to be honest, they should be costly to produce or maintain for their bearers [69]. For instance, in the western clade, conspicuous orange head colourations make individuals more visible to predators [70], so only high-quality lizards may survive older.

#### *4.2. Altitudinal Variation in Colouration*

We report that eyespots and throat colouration followed different altitudinal patterns; the number of blue eyespots was highest at low elevation and then followed a U-shaped trend with altitude, while saturation of throat colouration linearly increased with elevation. These discrepant patterns are hard to explain on the only basis of altitudinal variation in sexual selection pressure. Moreover, altitudinal colour variation was similar in males and females, supporting the idea that it is provoked by natural selection, not by sexual selection [71]. Our findings contrast with those reported for a population of the western clade in central Spain, where *P. algirus* lizards have less saturated throats at higher altitude [29,30]. In our study population, lizard throat was darker at high elevations, which has also been reported in other lizard species [28,32]. Darker colouration with ascending elevation is a likely consequence of high dermal melanin [59].

This pattern of altitudinal variation in lizard colour signals and, especially, why eyespots and throat colouration covaried differently with elevation, requires an explanation in which both colour signals are differentially affected by selective pressures. Blue eyespots are structural colourations, mainly produced by the combination of a thick and well-arranged layer of iridophores and basal eumelanin. Meanwhile, yellow and orange colourations are produced by pigments such as carotenoids and pterins. We discuss several hypotheses that could explain the altitudinal patterns reported.


could partially explain the highest number of eyespots at low elevations, hardly could it explain why investment in throat colouration increases with ascending altitude.


#### *4.3. Comparison with the Western Clade*

In the western clade, testosterone produces orange heads in adult males, but only an orange commissure in young males [38]. To produce or have an orange head is costly, as coloured individuals are more implied in fights, and colouration increases parasite susceptibility and reduces survival [30,35,36,41]. Why do males in the eastern clade not have orange heads during the breeding season? One possibility is paedomorphism, if adult lizards are retaining their youthful characteristics. This may occur because environmental conditions in the zones inhabited by the eastern clade favour an augment in the costs and/or decrease in the benefits associated with orange colouration. In the western clade, highly coloured males have larger home ranges that overlap with more females and so are more successful in mate acquisition [35,37,41]. The eastern clade inhabits more arid zones, with less plant cover and probably less food availability [42,95]. This might influence the social behaviour of this lizard (which is unstudied in the eastern clade), needing larger territories in which the control of several female territories is more difficult, and the contact with other males is rarer. Therefore, both intra- and intersexual selection could be lower in the eastern clade and so the benefits of an orange head could decrease. It is also possible that lizards in the eastern clade are more exposed to predators, given that they live in an opener environment [95]. Higher risk of predation would select against colourful males [96]. These questions remain open and show that *P. algirus* lizards offer a valuable opportunity to study and understand the evolution of nuptial colouration in lizards.

**Author Contributions:** Conceptualization, G.M.-R.; methodology, G.M.-R., S.R. and M.C.; resources: G.M.-R.; investigation: S.R., F.J.Z.-C. and M.C.; data curation: S.R., M.C. and G.M.-R.; formal analysis, G.M.-R. and S.R.; writing, original draft, G.M.-R. and S.R.; writing, additional inputs, F.J.Z.-C. and M.C.; supervision, G.M.-R.; project administration, G.M.-R.; funding acquisition, G.M.-R. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was funded by the Spanish government (Ministerio de Ciencia e Innovación) with funds from the European Union (project CGL2009-13185). FJZC (AP2009-3505) and SR (AP2009- 1325) were supported by two pre-doctoral grants from the Ministerio de Ciencia e Innovación (FPU programme).

**Institutional Review Board Statement:** The research was conducted under both Junta de Andalucía and National Park of Sierra Nevada research permits (references GMN/GyB/JMIF, ENSN/JSG/JEGT/MCF, ENSN/JSG/BRL/MCF, SGMN/GyB/JMIF, and SSA/SI/MD/ps) issued to the authors.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are openly available in FigShare at https://doi.org/10.6084/m9.figshare.14376698.v1, accessed on 6 April 2021.

**Acknowledgments:** We thank the personnel from the Espacio Natural de Sierra Nevada for their constant support.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**


**Table A1.** Sample sizes per sex (males|females) and per altitude for each variable.

#### **Appendix B**

#### **References**

