Research on animal colouration has grown exponentially in the last decade thanks to multidisciplinary approaches. Most studies are focused on trade-offs between communication and mimicry, which represent the two main constraints and drivers of the evolution of body colourations. Reptiles are excellent model species for investigating this field of study and lizards in particular show great variability of body colourations and their functions. We studied the lizard Podarcis siculus, analysing the variations of dorsal colour of three populations and obtained clear patterns of seasonal and ontogenetical variation of dorsal colour. According to baseline colour, males were greener and brighter than females, although no difference in saturation was recorded. According to seasonal variations, analyses showed that both sexes significantly vary in colour over the year: males reached higher peaks of hue and saturation later than females during spring, while females showed higher peaks of brightness and reached earlier similarly to hue and saturation. Ontogenetic variations were recorded only in males, which become greener, less bright and saturated with growing size. Therefore, our results suggest the occurrence of two opposing strategies in colour expression between sexes: males’ dorsal colouration plays a major role in communication, while females are more crypsis-oriented.

Research on animal colouration has been carried out intensively in the past few decades and nowadays it is a fast-growing interdisciplinary branch of zoology that involves a wide variety of biological sectors, namely genetics (Fulgione et al., 2004; Brenig et al., 2013), developmental biology (Wilson et al., 2007), reproductive biology (Wellenreuther et al., 2014) and feeding ecology (Schaefer et al., 2008), as well as all of their reciprocal interconnections. As a matter of fact, the evolution of body colouration is more often the result of the interaction between synergic and contrasting pressures rather than the outcome of one single factor (Endler, 1981). Therefore, it has been approached under various perspectives with the aim of investigating what drives the evolution and selection of specific body colourations. The task of this field of research is to understand the mechanisms involved in colour expression and perception in order to define the role of body colourations as a means of communication (Cuthill et al., 2017). Physiological studies focused largely on the mechanisms and constraints of colour expression, which is mainly determined by the innate production of pigments or, alternatively, by their acquisition from the environment via the trophic chain (Hill et al., 1994; Shawkey and D'Alba, 2017). Conspicuously coloured animals have to compromise between the communicative efficiency of a manifest colouration and the energetic cost of the production of pigments (Hill and Montgomerie, 1994). Another aspect of colour as a communication tool consists of the receiver's sensory perception that could be involved in a process of close evolutionary correlation between the colour pigments of the emitter and the visual receptors of the receiver (Lind et al., 2017; Renoult et al., 2017). Body colouration is also a highly informative trait for both inter- and intraspecific interactions as it may serve as an honest signal to warn potential threats (aposematism) (Mappes et al., 2005), attract mating partners (female choice based on the potential quality of the male) (Slagsvold and Lifjeld, 1988; LeBas and Marshall, 2000; Bajer et al., 2010) or to discriminate the sex of conspecifics (Cooper and Burns, 1987). Contrastingly, conspicuous colourations tend to increase the chance of being detected by predators causing fitness reduction. Thus, the evolution of cryptic and disruptive colour patterns, alongside deterrent or elusive behaviours, can be the adaptive response to compensate higher predation rate (Forsman, 1995; Forsman and Shine, 1995; Stuart-Fox et al., 2003).

Reptiles show noticeable variability of body colourations and patterns, and hence they are excellent model species in this research perspective. In the common adder (Vipera berus) melanism is positively correlated to higher reproductive success thanks to increased thermoregulatory efficiency, notwithstanding its negative effects on survivorship due to a higher chance to be detected by avian predators (Andrén and Nilson, 1981). Some lizards, as reported for Sceloporus grammicus, discriminate among individuals belonging to the same or different groups thanks to specific colour polymorphisms. Such variability is often characteristic of ventral parts of the body such as belly and throat, which are typically less visible than dorsal ones, thus compromising with potential predation risks (Bastiaans et al., 2014). Colour polymorphisms also influence social interactions and contests. The red males of the polymorphic Australian painted dragon Ctenophorus pictus tend to win most contests against yellow competitors; in keeping with the fact that red is commonly considered a sign of dominance, in this case, red colour has also the effect of reducing the aggressive response of the yellow counterparts (Healey et al., 2007).

The Italian wall lizard Podarcis siculus is an endemic lacertid of the Italian peninsula and the Adriatic basin and an efficient invasive species in other European and not European countries (Burke et al., 2002; Ribeiro and Sá-Sousa, 2018; Damas-Moreira et al., 2019). This lizard is dorsally green-brown for most of its activity period (March to October) and males especially show a green to azure gular region at the peak of their reproductive season (Corti et al., 2011). Its spreading success within and outside its original natural distribution range can be explained by its great adaptability in terms of trophic ecology (Zuffi and Giannelli, 2013), thermoregulatory efficiency (Kapsalas et al., 2016), and tolerance to anthropogenic environmental modifications (Mangiacotti et al., 2013). Although body colourations within the Podarcis genus have been widely investigated, most studies on Italian species have been carried out on P. muralis and its polymorphisms (Capula et al., 2009; Sacchi et al., 2013) while P. siculus has been investigated to a lesser extent under this perspective (Fulgione et al., 2008). Nonetheless, seasonal colour variation has been analysed in a recent paper, which showed that dorsal colourations match the background colour of the environment according to seasonal changes in vegetation colour, supporting the hypothesis that dorsal colouration can play a significant role in camouflage (Pellitteri-Rosa et al., 2020). However, dorsal (green) colourations in Lacertids are often involved in intersexual competition (Martín and López, 2010), where the intensity of the colouration correlates with male quality and forecasts the hierarchies of dominance in males (Olsson, 1994; Martín and López, 2009). Furthermore, darker colourations such as melanism are sometimes related to stronger immunity responses and resistance to ectoparasites (Baeckens and Van Damme, 2018). Since the breeding cycle of the Italian wall lizard matches the seasonal chromatic peak of green in local vegetation (males are most active in spring), the hypothesis that the dorsal colouration might play a role in intraspecific communication, although unlikely given that dorsal colourations are mostly camouflage-oriented, cannot be excluded. Notably, in earlier studies, neither sexual differences in dorsal colouration over the season nor the effect of body size on colour intensity have been investigated in detail (Pellitteri-Rosa et al., 2020). Therefore, this paper will focus on seasonal dorsal chromatic variations in three populations of P. siculus from northern and southern Italy, and on the possible effect of sex and body size on such variations.

The first step of the analysis consisted in calculating the mean values for hue, saturation and value and then computing their frequency distribution from the total number of pixels analysed; as shown in Fig. 1, males (m) tend to be generally greener than females (f), to show a higher amount of pigment, and to be slightly brighter.

Fig. 1.

Distribution frequencies of hue, saturation and brightness of body colouration of Italian wall lizards. Solid lines and dark grey areas are for males, dashed lines and pink areas are for females. Lines are means and area represent HDI95.

Fig. 1.

Distribution frequencies of hue, saturation and brightness of body colouration of Italian wall lizards. Solid lines and dark grey areas are for males, dashed lines and pink areas are for females. Lines are means and area represent HDI95.

The Bayesian mixed models supported the occurrence of a difference in the baseline colour (Mesor) between males and females, as the posterior distributions of the difference between sexes (m-f) deviated from zero for H and V, but not for S values (Fig. 2A). In detail, males were greener (Pm>f=0.75) and more brilliant (Pm>f=0.91) than females, but dorsal colouration was as saturated in males as in females (Pm>f=0.56).

Fig. 2.

Differential m-f posterior distribution values of the first component score of baseline colour (A), seasonal variation (B) and size (C) for each colour parameter (hue, saturation, value).

Fig. 2.

Differential m-f posterior distribution values of the first component score of baseline colour (A), seasonal variation (B) and size (C) for each colour parameter (hue, saturation, value).

The analysis of the Cosinor components on the score of the first component of hue (i.e. the seasonal effect) showed that the amplitude was not null, in both males and females (males: PA>0=0.999; females: PA>0=0.999, Table 1). This result supported the occurrence of a seasonal pattern of green expression in both sexes. However, males showed a seasonal amplitude larger than females (Pm>f=0.87, Fig. 2B), and reached the peak on 16 March (HDI95: 11 February, 12 April, Fig. 3A), whereas females did it on 21 February (HDI95: 8 December, 26 April, Fig. 3A). A further relevant effect on greenness was detected for body size, but only in males. Indeed, the hue PC1 score of males increased with body size, but it did not in females (Table 1). Consequently, the posterior probability of the differential effect of size between sexes deviated from zero (Pm>f=0.84, Fig. 2C). This result suggested a positive relation between size (age) and green expression only in males (Fig. 4A).

Fig. 3.

Effect of seasonality on hue (A), saturation (B) and value (C) differential patterns in males (solid lines) and females (dotted lines). Grey areas are 95% credibility intervals.

Fig. 3.

Effect of seasonality on hue (A), saturation (B) and value (C) differential patterns in males (solid lines) and females (dotted lines). Grey areas are 95% credibility intervals.

Fig. 4.

Effect of size (standardized SVL) on hue (A), saturation (B) and value (C) differential patterns in males (solid lines) and females (dotted lines). Grey areas are 95% credibility intervals.

Fig. 4.

Effect of size (standardized SVL) on hue (A), saturation (B) and value (C) differential patterns in males (solid lines) and females (dotted lines). Grey areas are 95% credibility intervals.

Table 1.

Fixed effect estimations (Cosinor parameters and size) on the colour spectra (scores of the first PC) of males and females of Italian wall lizard. For each parameter, the HSM (above), and HDI95 (below) are reported

Fixed effect estimations (Cosinor parameters and size) on the colour spectra (scores of the first PC) of males and females of Italian wall lizard. For each parameter, the HSM (above), and HDI95 (below) are reported
Fixed effect estimations (Cosinor parameters and size) on the colour spectra (scores of the first PC) of males and females of Italian wall lizard. For each parameter, the HSM (above), and HDI95 (below) are reported

The analysis of the Cosinor components for the PC1 of saturation confirmed that colour varied seasonally, as in this case also the amplitudes of both sexes deviated from zero (males: PA>0=0.99; females: PA>0=0.99, Table 1). Colour saturation increased during spring (Fig. 3B) and reached the maximum on 27 April in males (HDI95: 9 April, 10 May) and on 11 March in females (HDI95: 26 January, 6 June). As for hue, also the amplitude of saturation was higher in males than in females (Fig. 2B, Pm>f=0.83). Concerning the effect of size, we observed male saturation to decrease with increasing size (Table 1; Fig. 4B), while no remarkable effect was detected for females (Table 1). Consequently, the posterior probability of having a larger effect of size in males than in females was 77% (Fig. 2C).

The last analysis on the PC1 score of brightness supplied further evidence for seasonal changes in colour expression (Fig. 3C). The amplitudes of the Cosinor deviated from zero in both males and females (males: PA>0=0.99; females: PA>0=0.99, Table 1), but with an opposite pattern with respect to hue and saturation. Indeed, the amplitude was larger in females than in males (Pm>f=0.012, Fig. 2B). Males showed the maximum brightness on 28 February (HDI95: 7 February, 31 March), whereas females on 5 February (HDI95: 6 January, 12 March). As for hue and saturation, the effect of body size was observed only in males (Table 1): larger males were less bright than smaller ones (Fig. 4C). So, there was a marked difference between sexes regarding the effect of size on brightness (Pm>f=0.95, Fig. 2C).

Finally, a large portion of the variance unexplained by sex, season, and body size depended on the population, especially for saturation and brightness. These proportions were 28.7% for hue (HDI95: 9.6%-82.6%), 42.1% for saturation (HDI95: 17.2%-91.3%), and 65.4% for brightness (HDI95: 36.1%-97.6%). Notably, lizards sampled in central Italy (i.e. Tuscany) were greener, more saturated, but less brilliant than those samples in southern Italy (i.e. Calabria) (Fig. 5).

Fig. 5.

Posterior distribution values of the first component score of hue (circle), saturation (square) and brightness (triangle) for the three sampling sites (CA=Calci; BG=botanical garden; PA=Paola).

Fig. 5.

Posterior distribution values of the first component score of hue (circle), saturation (square) and brightness (triangle) for the three sampling sites (CA=Calci; BG=botanical garden; PA=Paola).

Our analysis demonstrates that the Italian wall lizard is characterised by a remarkable variability in terms of external body colouration, and that such variability is sexually and seasonally determined. Males were greener and more brilliant than females and showed a more pronounced seasonal pattern than females. In both males and females, the dorsal colouration appeared greener, more saturated and brighter during spring, whereas in autumn lizards of both sexes were browner, paler and duller. The role of ontogeny in the degree of green expression can be clearly detected only in males, which became greener, more saturated and less bright with increasing body size; while females did not show any clear pattern of colour expression concerning size.

Body colourations play a major role in the ecology of polymorphic species as they are deeply involved in a wide range of core functions and are strongly driven by selective pressures that can be summed up to the compromise between communication and crypsis (Galán, 2008; Cuthill et al., 2017). As a consequence, body colourations evolve according to the main function played and to the ecological constraints coped with. Therefore, conspicuous colourations will favour the maximization of intraspecific communication, notwithstanding the risk of being potentially detected by predators, while duller ones will be more likely to improve predator avoidance, thus relying on different ways to communicate with conspecifics (Cuthill et al., 2017). In agreement with past research on lizards that have outlined how conspicuous colourations have social significance for partner selection and male–male contests (LeBas and Marshall, 2000; Stuart-Fox et al., 2003; Husak et al., 2006; Bajer et al., 2010; Amdekar and Thaker, 2019), our data consistently indicate that males are more conspicuous than females for the whole duration of the reproductive season. This observation can be interpreted from different perspectives.

According to Pérez i de Lanuza et al. (2013) sexual size dimorphism (SSD) in lizards is largely driven by sexual selection as an outcome of male–male competition in sexually dimorphic lizards. Males of dimorphic species often show more complex, diverse and conspicuous colour patterns than females, which are instead duller and less visible; such sexual dichromatism also seems to be positively related to SSD. Hence it is reasonable to suggest that the males of P. siculus tend to be more conspicuously coloured with growing size as a potential signal of dominance towards competitors, thus favouring their potential mating success. Additionally, conspicuous colourations could also be involved in intersexual communication as male quality indicators (Cook et al., 2013; Pérez i de Lanuza et al., 2014). Coherently with these considerations, females of P. siculus are normally less colourful than males and their colour does not vary remarkably with age; this could be explained by the fact that an increase in colour expression, resulting in higher chances of being detected by predators, would reduce survivorship, and thus potential reproductive success (Bauwens and Thoen, 1981). From this perspective, it is known of other closely related European lacertids, like Psammodromus algirus and Zootoca vivipara, that male juveniles tend to show duller colourations than the adults, thus making them less detectable (Martin et al., 2013; Moreno-Rueda et al., 2021). Although it is relatively complex to determine the nature of such differences, it is reasonable to hypothesize that less conspicuous colourations are selected to favour concealment, similarly to what is displayed by females. Additionally, differential colour expression with dully coloured juveniles can also be addressed to the minimisation of predation risk, harassment by adults and to convey a signal of subordinance to adults (Lyon and Montgomerie, 1986; Hawkins et al., 2012).

When considering the effect of seasonality on colour expression patterns, it was possible to identify important patterns of periodic seasonal colour variations. Male lizards in particular showed more marked variations concerning all the three chromatic parameters of the HSV system, whereas such variations among females were mild (Table 2). The chromatic variation of males was almost in synchrony with the circannual period for H, S and V, thus indicating that colour expression (i.e. becoming greener and brighter) is tightly dependent on the phase of the activity period when specific colourations are necessary. Fig. 3A shows the effect of seasonality on colour variations for both sexes, where males tend to delay the maximum of green to the peak of the reproductive period (April–June) overlapping it with the minimum of saturation. As a consequence, males tend to be greener and brighter than females that, on the contrary, are darker throughout the whole activity period. According to these findings, it is reasonable to assume that in the Italian wall lizard two specific seasonal patterns can be identified: although both sexes show similar patterns of seasonal chromatic variations, on the one hand, males become globally more visible until the peak of reproduction, thus favouring the communicative aspect of colour expression, on the other hand females are generally more cryptic and less detectable (Marshall and Stevens, 2014; Marshall et al., 2015). This hypothesis is consistent with reported hormone seasonal variations of this species, where hormones reach their peaks simultaneously and overlap with the peak of green expression recorded in our research. Such hormones are in fact significantly involved in reproduction as both catecholamines and corticosterone (De Falco et al., 2004) and testosterone (Andò et al., 1990) are associated with aggressive behaviour that is most frequent during the reproductive period and especially in males. Therefore, it is possible that seasonal chromatic variations towards greener colourations are effectively associated with rising hormone levels in both sexes of P. siculus.

Table 2.

Number of individuals collected and mean SVL (mm) in each population

Number of individuals collected and mean SVL (mm) in each population
Number of individuals collected and mean SVL (mm) in each population

Furthermore, despite not being the focus of our research, the analyses highlighted the occurrence of variation among the studied populations. Notably, lizards from southern populations (Calabria) were brighter, but less green and saturated than those from central Italy (Tuscany). Although the cause of this discrepancy is unclear, we can hypothesize two potential reasons: it might be determined by a phylogenetic difference between two populations belonging to the separate clades P. siculus campestris (Tuscany) and P. siculus siculus (Calabria) (Senczuk et al., 2017) and resulting in morphological discrepancy in dorsal patterns (Corti et al., 2011); alternatively, darker patterns of Tuscanian lizards could be induced by environmental factors such as lower solar radiation and soil temperature (Pinna, 2017), which may affect colour expression in order to improve the efficiency of thermoregulation (Andrén and Nilson, 1981). In fact, in many colour shifting taxa (i.e., reptiles, amphibians, crustaceans and fish), temperature affects melanocyte-stimulating hormone determining melanin dispersion (Stuart-Fox and Moussalli, 2009). As reported by Norris (1967), reptiles can actively become darker or lighter and match background colouration according to environmental temperature, thus compensating lower thermal regimes with darker colourations to increase thermoregulatory efficiency. Hence, a similar adaptation can be addressed to the case of P. siculus northern populations that, under lower radiation and temperature conditions, might become darker to compensate for potentially less favourable environmental conditions.

In conclusion, our research points out that P. siculus is characterised by marked colour variations throughout the year, with distinct patterns between males and females, thus resulting in two diametrically different strategies, where the former is more communication-oriented and the latter is more crypsis-oriented. Since this is one of the first studies on the annual chromatic variations of P. siculus, further research is needed to extend our knowledge of the physiological and ecological mechanisms that could determine such variations (namely hormone secretion or pigment acquisition from the environment) and the extent to which they are perceived in the lizards’ visual spectrum. However, our experimental protocol is limited to assessing colour variations in the visible spectrum, thus indicating that this lizard matches the overall background of the environment in accordance with the findings of Pellitteri-Rosa et al. (2020). From this perspective, it is reasonable to assume that such differences in dorsal colourations are mainly camouflage-oriented to reduce detection by predators, given their ability to perceive the wavelengths within the visible spectrum. It would be important for future research to provide solid evidence supporting the assumption that lizards’ colour expression and their predators’ visual perception overlap, thus indicating that dorsal colourations could have a major role in camouflage. On the other hand, provided that lizards are known for being able to see ultra violet (UV) wavelengths as well and that ventral/lateral colourations are commonly used to communicate with conspecifics, the extant limited knowledge about UV signalling in the species we studied (Pérez i de Lanuza and Font, 2006; Stapley and Whiting, 2006) makes it still necessary to investigate in future research whether its body colourations also express UV signals and whether they play any role in intraspecific communication.

Data sampling

We collected the field data in 2019 throughout the whole activity period of the Italian wall lizard, from March to October, in three populations from central and southern Italy; the first site was in Calci, Tuscany (CA, 43°43′N, 10°31′E), the second and third were both in the proximity of Cosenza (Calabria), respectively, the botanic garden of the University of Calabria (BG, 39°21′N, 16°13′E), and the town of Paola (PA, 39°21′N, 16°2′E). Sampling surveys were performed once a month in each site by five to eight researchers for three consecutive days each time. All lizards were captured by noosing or by hand and then moved to the laboratory for measurements and photographs. Firstly, individuals were sexed and measured for snout-to-vent length (SVL) and tail length (TL) using a ruler (1 mm accuracy). Tail injury or loss was also recorded along with the breaking point (if detectable). All juveniles (SVL <50 mm; Henle, 1988) were excluded from the data collection and released after capture.

According to published data and thanks to the effectiveness and strength of this approach in studying animal colourations (Stevens et al., 2007), we used digital photography to retrieve a thorough framework of dorsal colorations of the studied populations. The method is able to provide a fine distinction, which is highly repeatable and finer than the human eye (Villafuerte and Negro, 1998). Therefore, for each individual, we took high-resolution digital images of dorsal colouration using a Nikon D50 camera at a 1.2-million-pixel resolution, equipped with a Nikkor 60 mm AF-S Micro lens, and fixed on a stand at a distance of 18 cm. Each picture was taken adjacent to a GretagMacBeth Mini ColorChecker chart (24 colour references, 5.7×8.25 cm) in a 44×44 cm lightbox illuminated with two daylight 22 W circular neon tubes (Reporter 55100 Studio-kit). After data collection, all lizards were released in the spot where they were captured. Overall, we measured and photographed 910 distinct lizards, including 419 males and 491 females (Table 2).

Colour analysis

We performed the analysis using the red-blue-green (RGB) colour system because colour vision is a common trait of diurnal vertebrates; although it has been extensively pointed out that the visual spectra of the different vertebrates can be highly variable due to the presence of an additional class of cones in birds (González-Martín-Moro et al., 2017) and reptiles (Fleishman et al., 1993), colour vision is present in lizards and in their predators, such as mammals or birds of prey (Kelber et al., 2003; Potier et al., 2018). Hence, this approach, despite not being reliable enough to assess the whole visual spectrum of both emitters (lizards) and receivers (conspecifics and predators) thoroughly, is adequate to provide partial insight about the potential constraints driving colour expression in our model species. We analysed the pictures according to the method by Bergmann and Beehner (2008) and Sacchi et al. (2013). Firstly, the Camera plug-in for Adobe Photoshop CS3 was used to generate a new colour profile that adjusted the colour in the photographs (jpeg format) to the known colour levels in each square of the ColorChecker chart. Then, for each image, we selected the region of interest (ROI) using the ‘lazoo’ tool, and pixels corresponding to point of reflected light were removed using the package ‘magick’ in R (Ooms, 2018). Subsequently, we analysed on average 59,806±44,296 pixels for each image, ranging from 2469 to 452,305 for the dorsal region. Eventually, the RGB colour values were rearranged in the hue, saturation, and value (HSV) system. It differs from the cubic geometry of the RGB colour space because it consists of a cylindrical-coordinate system in which the hue is the angle around the vertical axis and corresponds to the colour lights, the saturation is the distance from the axis, and the brightness corresponds to the distance along the axis.

Statistical analysis

Colour variation in the Italian wall lizards was assessed in a three-step analysis (Sacchi et al., 2021). First, we generated the individuals’ frequency distributions of hue, saturation, and brightness by using values computed on the whole samples of pixels selected for each individual (HSV colour spectrum). Second, three principal component analyses (PCA), one each for H, S, and V, respectively, were used on the colour spectra, and the first components, explaining 30.7%, 49.0% and 67.2% of the total variance, were used as a proxy to summarise the inter-individual variability of colouration. The hue PC score accounted for the opposite variation of the orange-yellow and green interval in the hue colour spectrum; negative scores were associated with a higher peak in the orange-yellow interval whereas positive scores were associated with a higher peak in the green interval. The saturation PC score accounted for the increase of colour saturation with an increased score, whereas the value PC score accounted for the increase of brightness with higher scores (Fig. 6). Third, the PC scores were analysed through random intercept linear mixed models (LMM) including a single-component cosinor function to model the effect of the season (Refinetti et al., 2007; Cornelissen, 2014). Cosinor models were originally developed to model circadian rhythm in physiological processes (Halberg et al., 1967), but they can be also used to model periodic variations in ecological variables (Mangiacotti et al., 2019; Sacchi et al., 2020). In Cosinor models the response variable (Y) is assumed to depend on time (t) following a regular cycle, which is incorporated in a linear model through a cosine function:
where M is the MESOR (Midline Statistic Of Rhythm, i.e. the time-corrected mean of the response), A is the amplitude (maximum absolute deviation from MESOR), τ the period of the cycle, the acrophase (i.e. the timing of highest values), and e(t) the error term (Cornelissen, 2014). The model can be linearized by rewriting the formula: Y(t)=M+βx+γz+e(t); being x=cos(2πt/τ) and z=sin (2πt/τ) the cosinor terms, and β=Acos and γ=− Asin the cosinor coefficients (Cornelissen, 2014). In our model, the two cosinor terms entered the LMM as fixed effects, with time expressed as Julian date (1= 1 January) and τ=365 to account for circannual rhythms around the time-corrected mean of H, S, V values. Additional fixed effects were sex and body size (i.e. standardized SVL). The two-way interactions sex×size, sex×x, and sex×z were also added to account for possible differential effect of sex on size and season on colour expression. Population entered the model as a random effect. The dependent variables were the PC scores of H, S and V values, and three independent models were run. LMMs were fit in a Bayesian analytical framework available in the package JAGS 4.3.0 (http://mcmc-jags.sourceforge.net/), using flat priors for coefficients and intercept (μ=0 and σ=0.001), and uninformative half-Cauchy priors (x0=0, γ=25) for both σ²error and σ²population. For all models, Markov Chain Monte Carlo parameters were set as follows: number of independent chains=three; number of iterations=34,000; burning=4000; thinning=three. Convergence was checked and results from the posterior distribution are reported as the half sample mode (HSM, Bickel and Frühwirth, 2006) plus 95% (or 50%) highest density intervals (HDI95; HDI50; Kruschke, 2010). In Bayesian statistics, the HSM is a commonly used estimator of the central tendency of posterior probability distribution robust to outliers, whereas the HDI95 defines the interval that includes the parameter with 95% probability. Parameter values in the centre of the HDI have higher credibility than parameter values at the limits of the interval. Therefore, when the HDIs of two groups do not overlap, there is a credible evidence for different group means. By contrast, to the extent of the two groups’ HDIs overlap there is no credible evidence of difference between the means. When comparing two groups (e.g. males and females), we therefore reported in Table 1 (bold) the posterior probability of their difference being different (i.e. higher or lower) from zero. Data preparation, model settings, call to JAGS, and posterior elaborations were done in R 4.0 using the package R2jags (Su and Yajima, 2015), modeest (Poncet, 2012), and HDInterval (Meredith and Kruschke, 2018).
Fig. 6.

Principal component (PC) analysis on hue, saturation and brightness of body colouration in the Italian wall lizards. Upper panels report the amount of variance associated with the first ten PCs; lower panels show the effects of the PC1 from the lowest (black) to the highest (red) scores expressed as variation with respect the mean values (i.e. the mean spectrum reported in Fig. 1).

Fig. 6.

Principal component (PC) analysis on hue, saturation and brightness of body colouration in the Italian wall lizards. Upper panels report the amount of variance associated with the first ten PCs; lower panels show the effects of the PC1 from the lowest (black) to the highest (red) scores expressed as variation with respect the mean values (i.e. the mean spectrum reported in Fig. 1).

We would like to thank the contribution of two anonymous reviewers for the remarkable and thorough suggestions that allowed a significant improvement of this research article. We thank also Cristian Marchio, Pietro Giovannelli and Armando Santoro, because they also contributed to the collection of data in Calabria. The research was performed in accordance with Italian laws on use of animals for scientific purpose and with the necessary permit from the Ministero dell'Ambiente e della Tutela del Territorio e del Mare (MATTM) and from the Istituto Superiore per la Protezione Ricerca Ambientale (ISPRA) for capture, manipulation, sample collection and release of the animals. We would also like to thank our colleague Mattia Falaschi PhD for correcting an incorrect figure label (Fig. 1, page 2).

Author contributions

Conceptualization: F.S., M.A.L.Z., S.S., R.S.; Methodology: M.A.L.Z., A.J.C., M.M., S.S., R.S.; Formal analysis: M.A.L.Z., A.J.C., M.M., F.A., R.S.; Investigation: F.S., A.J.C., L.D.V., G.G., A.E.G., F.L.L., S.L., M.M., S.S., F.A., E.S., I.T., R.S.; Resources: F.S., M.A.L.Z.; Data curation: F.S., M.A.L.Z., A.J.C., M.M., S.S., F.A., E.S., R.S.; Writing - original draft: F.S., R.S.; Writing - review & editing: M.A.L.Z., A.J.C., M.M., S.S., E.S., R.S., L.D.V., G.G., A.E.G., F.L.L., S.L., F.A., I.T.; Supervision: M.A.L.Z.; Project administration: R.S.

Funding

This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

Amdekar
,
M. S.
and
Thaker
,
M.
(
2019
).
Risk of social colours in an agamid lizard: implications for the evolution of dynamic signals
.
Biol. Lett.
15
,
20190207
.
Andò
,
S.
,
Panno
,
M. L.
,
Ciarcia
,
G.
,
Imbrogno
,
E.
,
Buffone
,
M.
,
Beraldi
,
E.
,
Sisci
,
D.
,
Angelini
,
F.
and
Botte
,
V.
(
1990
).
Plasma sex hormone concentrations during the reproductive cycle in the male lizard, Podarcis s. sicula
.
J. Reprod. Fert.
90
,
353
-
360
.
Andrén
,
C.
and
Nilson
,
G.
(
1981
).
Reproductive success and risk of predation in normal and melanistic colour morphs of the adder, Vipera berus
.
Biol. J. Linn. Soc.
15
,
235
-
246
.
Baeckens
,
S.
and
Van Damme
,
R.
(
2018
).
Immunocompetence and parasite infestation in a melanistic and normally-coloured population of the lacertid lizard, Podarcis siculus
.
Amphib-Reptilia
39
,
471
-
478
.
Bajer
,
K.
,
Molnár
,
O.
,
Török
,
J.
and
Herczeg
,
G.
(
2010
).
Female European green lizards (Lacerta viridis) prefer males with high ultraviolet throat reflectance
.
Behav. Ecol. Sociobiol.
64
,
2007
-
2014
.
Bastiaans
,
E.
,
Bastiaans
,
M. J.
,
Morinaga
,
G.
,
Gaytán
,
J. G. C.
,
Marshall
,
J. C.
,
Bane
,
B.
,
de la Cruz
,
F. M.
and
Sinervo
,
B.
(
2014
).
Female preference for sympatric vs. Allopatric male throat color morphs in the mesquite lizard (Sceloporus grammicus) species complex
.
PLoS ONE
9
,
e93197
.
Bauwens
,
D.
and
Thoen
,
C.
(
1981
).
Escape tactics and vulnerability to predation associated with reproduction in the lizard Lacerta vivipara
.
J. Anim. Ecol.
50
,
733
-
743
.
Bergman
,
T. J.
and
Beehner
,
J. C.
(
2008
).
A simple method for measuring colour in wild animals: validation and use on chest patch colour in geladas (Theropithecus gelada)
.
Biol. J. Linn. Soc.
94
,
231
-
240
.
Bickel
,
D. R.
and
Frühwirth
,
R.
(
2006
).
On a fast, robust estimator of the mode: comparisons to other robust estimators with applications
.
Comput. Stat. Data Anal.
50
,
3500
-
3530
.
Brenig
,
B.
,
Beck
,
J.
,
Floren
,
C.
,
Bornemann-Kolatzki
,
K.
,
Wiedemann
,
I.
,
Hennecke
,
S.
,
Swalve
,
H.
and
Schütz
,
E.
(
2013
).
Molecular genetics of coat colour variations in white galloway and white park cattle
.
Anim. Genet
44
,
450
-
453
.
Burke
,
R. L.
,
Hussain
,
A. A.
,
Storey
,
J. M.
and
Storey
,
K. B.
(
2002
).
Freeze tolerance and supercooling ability in the italian wall lizard, Podarcis sicula, introduced to long island, New York
.
Copeia
2002
,
836
-
842
.
Capula
,
M.
,
Rugiero
,
L.
and
Luiselli
,
L.
(
2009
).
Ecological correlates of colour pattern polymorphism along the transition zone between two morphs of the common wall lizard, Podarcis muralis
.
Ital. J. Zool.
76
,
376
-
380
.
Cook
,
E. G.
,
Murphy
,
T. G.
and
Johnson
,
M. A.
(
2013
).
Colorful displays signal male quality in a tropical anole lizard
.
Naturwissenschaften
100
,
993
-
996
.
Cooper
,
W. E.
and
Burns
,
N.
(
1987
).
Social significance of ventrolateral coloration in the fence lizard, Sceloporus undulatus
.
Anim. Behav.
35
,
526
-
532
.
Cornelissen
,
G.
(
2014
).
Cosinor-based rhythmometry
.
Theor. Biol. Med. Model.
11
,
16
.
Corti
,
C.
,
Capula
,
M.
,
Luiselli
,
L.
,
Razzetti
,
E.
and
Sindaco
,
R.
(
2011
).
Podarcis siculus (Rafinesque-Schmaltz, 1810)
. In
Fauna d'Italia, Reptilia
(ed.
C.
Corti
,
M.
Capula
,
L.
Luiselli
,
E.
Razzetti
and
R.
Sindaco
), pp.
407
-
417
.
Calderini Publisher
.
Cuthill
,
I. C.
,
Allen
,
W. L.
,
Arbuckle
,
K.
,
Caspers
,
B.
,
Chaplin
,
G.
,
Hauber
,
M. E.
,
Hill
,
G. E.
,
Jablonski
,
N. G.
,
Jiggins
,
C. D.
,
Kelber
,
A.
et al. 
(
2017
).
The biology of color
.
Science
357
,
eaan0221
.
Damas-Moreira
,
I.
,
Riley
,
J. L.
,
Harris
,
D. J.
and
Whiting
,
M. J.
(
2019
).
Can behaviour explain invasion success? A comparison between sympatric invasive and native lizards
.
Anim. Behav.
151
,
195
-
202
.
De Falco
,
M.
,
Sciarrillo
,
R.
,
Virgilio
,
F.
,
Fedele
,
V.
,
Valiante
,
S.
,
Laforgia
,
V.
and
Varano
,
L.
(
2004
).
Annual variations of adrenal gland hormones in the lizard Podarcis sicula
.
J. Comp. Physiol.
190
,
675
-
681
.
Endler
,
J. A.
(
1981
).
An overview of the relationships between mimicry and crypsis
.
Biol. J. Linn. Soc.
16
,
22
-
25
.
Fleishman
,
L. J.
,
Loew
,
E. R.
and
Leal
,
M.
(
1993
).
Ultraviolet vision in lizards
.
Nature
365
,
397
.
Forsman
,
A.
(
1995
).
Opposing fitness consequences of colour pattern in male and female snakes
.
J. Evol. Biol.
8
,
53
-
70
.
Forsman
,
A.
and
Shine
,
R.
(
1995
).
The adaptive significance of colour pattern polymorphism in the Australian scincid lizard Lampropholis delicata
.
Biol. J. Linn. Soc.
55
,
273
-
291
.
Fulgione
,
D.
,
Milone
,
M.
,
Rippa
,
D.
,
Capriglione
,
T.
,
Guarino
,
F. M.
and
Odierna
,
G.
(
2004
).
Genetic analysis in a mainland melanic lizard (Podarcis sicula) population from Cilento and Vallo di Diano national park (Campania, Italy)
.
Ital. J. Zool.
71
,
25
-
28
.
Fulgione
,
D.
,
Guglielmi
,
S.
,
Odierna
,
G.
,
Rippa
,
D.
,
Caliendo
,
M. F.
and
Rastogi
,
R. K.
(
2008
).
Morphological differentiation and genetic structure in island lizard populations
.
Zool. Sci.
25
,
465
-
474
.
Galán
,
P.
(
2008
).
Ontogenetic and sexual variation in the coloration of the lacertid lizards Iberolacerta monticola and Podarcis bocagei. Do the females prefer the greener males?
Anim. Biol.
58
,
173
-
198
.
González-Martín-Moro
,
J.
,
Hernández-Verdejo
,
J. L.
and
Clement-Corral
,
A.
(
2017
).
The visual system of diurnal raptors: updated review
.
Arch. Rev. Esp. Oftalmol.
92
,
225
-
232
.
Halberg
,
F.
,
Tong
,
Y. L.
and
Johnson
,
E. A.
(
1967
).
Circadian system phase — an aspect of temporal morphology; procedures and illustrative examples
. In
The Cellular Aspects of Biorhythms
(ed.
H.
von Meyersbach
), pp.
20
-
48
.
Berlin, Heidelberg
:
Springer
.
Hawkins
,
G. L.
,
Hill
,
G. E.
and
Mercadante
,
A.
(
2012
).
Delayed plumage maturation and delayed reproductive investment in birds
.
Biol. Rev.
87
,
257
-
274
.
Healey
,
M.
,
Uller
,
T.
and
Olsson
,
M.
(
2007
).
Seeing red: morph-specific contest success and survival rates in a colour-polymorphic agamid lizard
.
Anim. Behav.
74
,
337
-
341
.
Henle
,
K.
(
1988
).
Dynamics and ecology of three Yugoslavian populations of the Italian wall lizard (Podarcis sicula campestris De Betta) (Reptilia: Lacertidae)
.
Zool. Anz
220
,
33
-
48
.
Hill
,
G. E.
and
Montgomerie
,
R.
(
1994
).
Plumage colour signals nutritional condition in the house finch
.
Proc. R. Soc. B Biol. Sci.
258
,
47
-
52
.
Hill
,
G. E.
,
Montgomerie
,
R.
,
Inouye
,
C. Y.
and
Dale
,
J.
(
1994
).
Influence of dietary carotenoids on plasma and plumage colour in the house finch: intra- and intersexual variation
.
Funct. Ecol.
8
,
343
-
350
.
Husak
,
J. F.
,
Macedonia
,
J. M.
,
Fox
,
S. F.
and
Sauceda
,
R. C.
(
2006
).
Predation cost of conspicuous male coloration in collared lizards (Crotaphytus collaris): an experimental test using clay-covered model lizards
.
Ethology
112
,
572
-
580
.
Kapsalas
,
G.
,
Gavriilidi
,
I.
,
Adamopoulou
,
C.
,
Foufopoulos
,
J.
and
Pafilis
,
P.
(
2016
).
Effective thermoregulation in a newly established population of Podarcis siculus in Greece: a possible advantage for a successful invader
.
Acta Herpetol.
11
,
111
-
118
.
Kelber
,
A.
,
Vorobyev
,
M.
and
Osorio
,
D.
(
2003
).
Animal colour vision – behavioural tests and physiological concepts
.
Biol. Rev.
78
,
81
-
118
.
Kruschke
,
J. K.
(
2010
).
Bayesian data analysis
.
Wiley Interdiscip. Rev. Cogn. Sci.
1
,
658
-
676
.
LeBas
,
N. R.
and
Marshall
,
J. N.
(
2000
).
The role of colour in signalling and male choice in the agamid lizard Ctenophorus ornatus
.
Proc. R. Soc. B Biol. Sci.
267
,
445
-
452
.
Lind
,
O.
,
Henze
,
M. J.
,
Kelber
,
A.
and
Osorio
,
D.
(
2017
).
Coevolution of coloration and colour vision?
Philos. Trans. R. Soc. B Biol. Sci.
372
,
20160338
.
Lyon
,
B. E.
and
Montgomerie
,
R. D.
(
1986
).
Delayed plumage maturation in passerine birds: reliable signaling by subordinate males?
Evolution
40
,
605
-
615
.
Mangiacotti
,
M.
,
Scali
,
S.
,
Sacchi
,
R.
,
Bassu
,
L.
,
Nulchis
,
V.
and
Corti
,
C.
(
2013
).
Assessing the spatial scale effect of anthropogenic factors on species distribution
.
PLoS ONE
8
,
e67573
.
Mangiacotti
,
M.
,
Pezzi
,
S.
,
Fumagalli
,
M.
,
Coladonato
,
A. J.
,
d'Ettorre
,
P.
,
Leroy
,
C.
,
Bonnet
,
X.
,
Zuffi
,
M. A. L.
,
Scali
,
S.
and
Sacchi
,
R.
(
2019
).
Seasonal variations in femoral gland secretions reveals some unexpected correlations between protein and lipid components in a lacertid lizard
.
J. Chem. Ecol.
45
,
673
-
683
.
Mappes
,
J.
,
Marples
,
N.
and
Endler
,
J. A.
(
2005
).
The complex business of survival by aposematism
.
Trends Ecol. Evol.
20
,
598
-
603
.
Marshall
,
K. L. A.
and
Stevens
,
M.
(
2014
).
Wall lizards display conspicuous signals to conspecifics and reduce detection by avian predators
.
Behav. Ecol.
25
,
1325
-
1337
.
Marshall
,
K. L. A.
,
Philpot
,
K. E.
and
Stevens
,
M.
(
2015
).
Conspicuous male coloration impairs survival against avian predators in Aegean wall lizards, Podarcis erhardii
.
Ecol. Evol.
5
,
4115
-
4131
.
Martín
,
J.
and
López
,
P.
(
2009
).
Multiple color signals may reveal multiple messages in male Schreiber's green lizards, Lacerta schreiberi
.
Behav. Ecol. Sociobiol.
63
,
1743
-
1755
.
Martín
,
J.
and
López
,
P.
(
2010
).
Multimodal sexual signals in male ocellated lizards Lacerta lepida: Vitamin E in scent and green coloration may signal male quality in different sensory channels
.
Naturwissenschaften
97
,
545
-
553
.
Martin
,
M.
,
Meylan
,
S.
,
Gomez
,
D.
and
Le Gaillard
,
J.-F.
(
2013
).
Ultraviolet and carotenoid-based coloration in the viviparous lizard Zootoca vivipara (Squamata: Lacertidae) in relation to age, sex, and morphology
.
Biol. J. Linn. Soc.
110
,
128
-
141
.
Meredith
,
M.
and
Kruschke
,
J.
(
2018
).
Package ‘HDInterval’. Highest (Posterior) Density Intervals
.
Moreno-Rueda
,
G.
,
Reguera
,
S.
,
Zamora-Camacho
,
F. J.
and
Comas
,
M.
(
2021
).
Inter-individual differences in ornamental colouration in a mediterranean lizard in relation to altitude, season, sex, age, and body traits
.
Diversity
13
,
158
.
Norris
,
K. S.
(
1967
).
Colour adaptation in desert reptiles and its thermal relationships
. In
Lizard ecology – A Symposium
(ed.
W.W.
Milstead
), pp.
162
-
229
.
Columbia, MO
:
University of Missouri Press
.
Olsson
,
M.
(
1994
).
Nuptial coloration in the sand lizard, Lacerta agilis: an intra-sexually selected cue to fighting ability
.
Anim. Behav.
48
,
607
-
613
.
Ooms
,
J.
(
2020
).
magick: advanced graphics and image-processing in R. R packageversion 2.3. Available from: https://CRAN.R-project.org/package=magick (accessed 20 April 2020)
.
Pellitteri-Rosa
,
D.
,
Gazzola
,
A.
,
Todisco
,
S.
,
Mastropasqua
,
F.
and
Liuzzi
,
C.
(
2020
).
Lizard colour plasticity tracks background seasonal changes
.
Biol. Open
9
,
bio052415
.
Pérez i de Lanuza
,
G.
and
Font
,
E.
(
2006
).
Lizard blues: blue body colouration and ultraviolet polychromatism in lacertids
.
Rev. Esp. Herp.
24
,
67
-
84
.
Pérez i de Lanuza
,
G.
,
Font
,
E.
and
Monterde
,
J. L.
(
2013
).
Using visual modelling to study the evolution of lizard coloration: sexual selection drives the evolution of sexual dichromatism in lacertids
.
J. Evol. Biol.
26
,
1826
-
1835
.
Pérez i de Lanuza
,
G.
,
Carazo
,
P.
and
Font
,
E.
(
2014
).
Colours of quality: structural (but not pigment) coloration informs about male quality in a polychromatic lizard
.
Anim. Behav.
90
,
73
-
81
.
Pinna
,
S.
(
2017
).
Lineamenti di Climatologia
,
Canterano, RM
,
Italy
):
Aracne editrice
.
Poncet
,
P.
(
2012
).
Modeest: mode estimation. R package version, 2(2)
.
Potier
,
S.
,
Mitkus
,
M.
and
Kelber
,
A.
(
2018
).
High resolution of colour vision, but low contrast sensitivity in a diurnal raptor
.
Proc. R. Soc. B Biol. Sci.
285
,
20181036
.
Refinetti
,
R.
,
Cornélissen
,
G.
and
Halberg
,
F.
(
2007
).
Procedures for numerical analysis of circadian rhythms
.
Biol. Rhythm. Res.
38
,
275
-
325
.
Renoult
,
J. P.
,
Kelber
,
A.
and
Schaefer
,
H. M.
(
2017
).
Colour spaces in ecology and evolutionary biology
.
Biol. Rev.
92
,
292
-
315
.
Ribeiro
,
R.
and
Sá-Sousa
,
P.
(
2018
).
Where to live in Lisbon: urban habitat used by the introduced Italian wall lizard (Podarcis siculus)
.
Basic Appl. Herpetol.
32
,
57
-
70
.
Sacchi
,
R.
,
Pellitteri-Rosa
,
D.
,
Bellati
,
A.
,
Di Paoli
,
A.
,
Ghitti
,
M.
,
Scali
,
S.
,
Galeotti
,
P.
and
Fasola
,
M.
(
2013
).
Colour variation in the polymorphic common wall lizard (Podarcis muralis): analysis using the RGB colour system
.
Zool. Anz. A. J. Comp. Zool.
252
,
431
-
439
.
Sacchi
,
R.
,
Mangiacotti
,
M.
,
Scali
,
S.
,
Coladonato
,
A. J.
,
Pitoni
,
S.
,
Falaschi
,
M.
and
Zuffi
,
M. A. L.
(
2020
).
Statistical methodology for the evaluation of leukocyte data in wild reptile populations: a case study with the common wall lizard (Podarcis muralis)
.
PLoS ONE
15
,
e0237992
.
Sacchi
,
R.
,
Cancian
,
S.
,
Ghia
,
D.
,
Fea
,
G.
and
Coladonato
,
A.
(
2021
).
Color variation in signal crayfish Pacifastacus leniusculus
.
Curr. Zool.
67
,
35
-
43
.
Schaefer
,
H. M.
,
McGraw
,
K.
and
Catoni
,
C.
(
2008
).
Birds use fruit colour as honest signal of dietary antioxidant rewards
.
Funct. Ecol.
22
,
303
-
310
.
Senczuk
,
G.
,
Colangelo
,
P.
,
De Simone
,
E.
,
Aloise
,
G.
and
Castiglia
,
R.
(
2017
).
A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus
.
BMC Evol. Biol.
17
,
6
.
Shawkey
,
M. D.
and
D'Alba
,
L.
(
2017
).
Interactions between colour-producing mechanisms and their effects on the integumentary colour palette
.
Philos. Trans. R. Soc. B Biol. Sci.
372
,
20160536
.
Slagsvold
,
T.
and
Lifjeld
,
J. T.
(
1988
).
Plumage colour and sexual selection in the pied flycatcher Ficedula hypoleuca
.
Anim. Behav.
36
,
395
-
407
.
Stapley
,
J.
and
Whiting
,
M. J.
(
2006
).
Ultraviolet signals fighting ability in a lizard
.
Biol. Lett.
2
,
169
-
172
.
Stevens
,
M.
,
Párraga
,
C. A.
,
Cuthill
,
I. C.
,
Partridge
,
J. C.
and
Troscianko
,
T. S.
(
2007
).
Using digital photography to study animal coloration
.
Biol. J. Linn. Soc.
90
,
211
-
237
.
Stuart-Fox
,
D.
and
Moussalli
,
A.
(
2009
).
Camouflage, communication and thermoregulation: lessons from colour changing organisms
.
Phil. Trans. R. Soc. Biol. Sci.
364
,
463
-
470
.
Stuart-Fox
,
D. M.
,
Moussalli
,
A.
,
Marshall
,
N. J.
and
Owens
,
I. P. F.
(
2003
).
Conspicuous males suffer higher predation risk: visual modelling and experimental evidence from lizards
.
Anim. Behav.
66
,
541
-
550
.
Su
,
Y. S.
and
Yajima
,
S.
(
2015
)
Package ‘R2jags’. R package version 0.03-08.
http://CRAN.R-project.org/package=R2jags.
Villafuerte
,
R.
and
Negro
,
J. J.
(
1998
).
Digital imaging for colour measurement in ecological research
.
Ecol. Lett.
1
,
151
-
154
.
Wellenreuther
,
M.
,
Svensson
,
E. I.
and
Hansson
,
B.
(
2014
).
Sexual selection and genetic colour polymorphisms in animals
.
Mol. Ecol.
23
,
5398
-
5414
.
Wilson
,
D.
,
Heinsohn
,
R.
and
Endler
,
J. A.
(
2007
).
The adaptive significance of ontogenetic colour change in a tropical python
.
Biol. Lett.
3
,
40
-
43
.
Zuffi
,
M. A. L.
and
Giannelli
,
C.
(
2013
).
Trophic niche and feeding biology of the Italian wall lizard, Podarcis siculus campestris (De Betta, 1857) along western Mediterranean coast
.
Acta Herpetol.
8
,
35
-
39
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.