Ain’t going down without a fight: state-and environment-dependence of antipredator defensive aggressive personalities in Carpetan rock lizard

Aggression is one of the most frequently studied behavioural traits across a wide range of taxa; however, most studies evaluate aggressive behaviour in a social context, in which aggressive interactions between conspecifics are motivated by resource control (offensive or social aggression). However, in an antipredator context, the primary role of aggression is defence (defensive or antipredator aggression). Although the neuroendocrinology of antipredator aggression is often studied in domesticated and laboratory animals, how environment and individual state affect this behavioural trait in the wild is largely unknown. Here, by conducting a manipulative experiment, we tested whether (i) consistent between-individual differences (i.e. animal personality) are present in antipredator aggression in adult male Carpetan rock lizards (Iberolacerta cyreni) and (ii) short-term environmental changes (presence vs. absence of predator cues) and differences in individual state (body length, head size, hind limb length) affect individual mean behaviour (i.e. behavioural type). We found moderate-high repeatability in antipredator aggression (willingness to bite a human), indicating the presence of animal personality in this behavioural trait. Lizards were on average more defensive in the presence of predator cues; furthermore, short-legged males showed higher antipredator aggression than long-legged males in the presence of predator cues, probably as an attempt to balance their decreased escape speed. Larger (~ older) males were more defensive than smaller ones, probably due to their increased fighting ability. We conclude that antipredator aggression is an important part of an individual’s behavioural repertoire and its expression is driven by both environmental situation and individual state. Antipredator/defensive aggression is not the primary antipredatory response; however, when other ways of escape are not possible, actually hurting the predator could be the only way of survival. While this behaviour obviously has substantial effects on fitness, it is severely understudied compared to social/offensive aggression. In a manipulative experiment, we found that there are consistent between-individual differences in antipredator aggression (i.e. willingness to bite during handling) of adult male Carpetan rock lizards (Iberolacerta cyreni), supporting the presence of animal personality and suggesting that this behavioural trait might respond to natural selection. Furthermore, short-term environmental variation (i.e. presence vs. absence of predator cues) in interaction with individual state affected antipredator aggression of individuals, emphasising the ecological and evolutionary relevance of this behaviour.


Introduction
Aggression is broadly defined as any overt fighting behaviour or signal of imminent behaviour with the capacity of harm (Moyer 1968;Huntingford 1976;Nelson 2006;Cain et al. 2011). Evolutionary behavioural ecology of consistent individual behavioural differences across time and ecological situations in a population (i.e. animal personality) became intensively studied during the last decades (for reviews and meta-Communicated by S. Joy Downes analyses see Réale et al. 2007;Sih and Bell 2008;Bell et al. 2009;Sih et al. 2012;Garamszegi et al. 2013;Niemelä and Dingemanse 2018), and aggression is among the most frequently studied behavioural traits in this field across several taxa (Mafli et al. 2011;Sih et al. 2015;Santostefano et al. 2016;Michelangeli et al. 2017;Szász et al. 2019). However, aggression is usually evaluated in a social context, where it is motivated by resource control, including territorial and dominance disputes (often termed 'offensive aggression'; see Nelson 2006). As members of the same species tend to utilise the same habitats and resources, the most common form of offensive aggression is the agonistic interaction between conspecifics (e.g. male-male combat; see McEvoy et al. 2013;Lee et al. 2014;Briffa et al. 2015); however, in a wider sense, the challenger does not necessarily has to be a conspecific (e.g. agonistic encounters between lions and hyenas over each other's prey; Watts and Holekamp 2008). Importantly, offensive aggressive encounters rarely result in death and the outcome of the interactions highly depends on the fighting capabilities of the opponents (Nelson 2006). In contrast, aggression in an antipredator context is primarily and directly motivated by imminent danger of death (often termed as 'defensive aggression'; see Nelson 2006;Cain et al. 2011). However, using the terms offensive vs. defensive aggression can be sometimes misleading, because (i) the weaker party in an offensive aggressive fight might actually defend itself via aggressive actions, while (ii) defensive aggression might include potentially dangerous prey behaving offensively towards a (yet) passive predator as a preventive measure. Therefore, we prefer and use the terms social (instead of offensive) and antipredator (instead of defensive) aggression in our paper. Antipredator aggression is typically triggered by nearcontact proximity of a potential predator. Antipredator attack (e.g. bite) is not a primary antipredatory response, as it occurs only if fleeing is not possible, or defensive threats, such as vocalisation, defensive posture or displaying teeth and/or claws are ineffective and the only feasible way to escape is trying to actually hurt the predator (Cowlishaw 1994;Nelson 2006;Blanchard et al. 2008). Hence, antipredator aggression could be seen as a 'last resort' against predation, similarly to another, yet completely opposite antipredator behavioural trait, tonic immobility (i.e. the lack of movement; see Gallup and Rager 1996;Edelaar et al. 2012;Horváth et al. 2019).
Neuroendocrine control of antipredator aggression and its relationship with stress is well studied, although mainly in domesticated and laboratory animals (e.g. rodents; see Rodgers and Depaulis 1982;Popova 2005;Pesce et al. 2011), which are often selected to show relatively little defence to handling. In contrast, how environment and individual state affect this behavioural trait in the wild is less studied (Cury de Barros et al. 2010;Baxter-Gilbert and Riley 2018). In behavioural ecology, an individual's state reflects any feature that affects the cost and benefits of its behavioural actions (Houston and McNamara 1999). These features may involve both inherently stable (e.g. sex differences, morphology) and labile features (e.g. energy reserves, health state, reproductive value). In addition, labile differences may include not only internal characteristics but the social environment (e.g. local density, behaviour of social partners) and ecological environment (e.g. predator, competitors, parasites) as well (see Sih et al. 2015). As mentioned above, predatory imminence is the main stimulus activating antipredator aggression; however, some stimuli may represent different range of imminence of threat. For instance, presence of active predators usually triggers high level of defensive behaviour, while passive predators do not (Blanchard et al. 2008). Similarly, one would expect that stimuli representing high level of risk (e.g. visual) should elicit stronger antipredator responses and defensiveness than cues representing the potential presence of predators without immediate risk (e.g. chemical cues) (Amo et al. 2006). Nonetheless, predator chemical cues (kairomones) were repeatedly shown to be enough to trigger glucocorticoid stress hormone release (e.g. corticosterone; Cockrem and Silverin 2002;Thomas et al. 2006;Trompeter and Langkilde 2011), and even short-term exposure to predator odour affects antipredator behaviour substantially (de Paula et al. 2005;Teyssier et al. 2014;Lloren et al. 2019) and may trigger defensive behaviour (Kalynchuk et al. 2004). Additionally, individual differences in antipredator aggression depending on labile intrinsic state variables are expected to change fast, as the underlying mechanisms are highly variable in time. For instance, in aggressive encounters between conspecifics, the outcome of aggressive interactions is often affected by body size and weight (see Rowland 1989;Nowbahari et al. 1999), and thus larger prey individuals might rely more on antipredator aggression than smaller ones. Furthermore, weaponry can be a predictor of fight success as well. In lizards, head size is an important proxy of fighting ability (Gvozdik and Van Damme 2003), and thus large-headed individuals likely defend themselves more effectively from predators than smallheaded conspecifics. On the other hand, individuals better at escaping might employ an antipredator strategy where defensiveness is less prominent. In lizards, relative limb length was shown to positively predict locomotor performance (Bauwens et al. 1995); hence, individuals with relatively shorter limbs might rely more on antipredator aggression than their conspecifics with longer limbs.
Here, we tested (i) whether antipredator aggression personality is present in adult male Carpetan rock lizards (Iberolacerta cyreni), and (ii) whether short-term environmental changes (presence vs. absence of predator cues) and differences in individual state (body size, head size, and hind limb length) affect antipredator aggression behavioural type (individual mean behaviour). To this end, we performed a manipulative experiment in which we repeatedly tested each individual's defensiveness both in the presence and absence of predator chemical cues (i.e. scent of smooth snake, Coronella austriaca). As detecting predator cues connotes a stressful situation, we predicted lizards to show stronger antipredator aggression after being exposed to smooth snake scent. According to the state-dependent safety principle, individuals with higher state are expected to show higher behavioural activity for its benefits, because their high state allows them to deal with the increased costs (Luttbeg and Sih 2010). Hence, we expected males with higher state regarding antipredator aggression (larger individuals with larger heads; Gvozdik and Van Damme 2003) to be more defensive than lower-state conspecifics. We also expected males with higher state regarding escape performance (individuals with longer limbs; Bauwens et al. 1995;Vanhooydonck et al. 2001) to be less defensive than males with lower expected escape abilities.

Study animals
We noosed 25 adult male I. cyreni between 15 and 17 May of 2016 (coinciding with the peak of the mating period of this species) at the 'Alto del Telégrafo' peak (Sierra de Guadarrama, Madrid prov., Spain, 1900 m asl). The animals were transported to the 'El Ventorrillo' field station (MNCN-CSIC), approximately 5 km from the capture site, where they were housed individually outdoor in grey opaque boxes (57 cm × 37 cm × 30 cm; length, width, height, respectively). In the boxes, we used a layer of coconut fibres as substrate (2-3 cm thick, approximately), and we provided shelters (one per box) made of fibreboards (20 cm × 15 cm × 1 cm; length, width, height, respectively). The shelters were made without bottom and were open from the front, providing a suitable hiding place for the lizards, but the shelters could be also quickly removed with minimal disturbance to the animals. Between the experimental trials, we covered the boxes with fine metal mesh to protect animals from bird predators. Individuals spent 10-12 days in the boxes used for the experiment to acclimate before the treatments started. Water and food (house crickets, Acheta domestica) were provided ad libitum. At the end of the experiment, all lizards were released at their original capture sites, without any sign of injury.

Individual traits
We measured snout-vent length (SVL; 67.39 ± 3.66; mean ± standard deviation (SD)), head size, and the length of the limbs of lizards using a digital calliper to the nearest 0.01 mm. To characterise head size, we ran a principal component analysis on four head measures (head length, head width, jaw width and head height) that produced a single principal component with strong positive loadings (61% variation explained; all factor loadings > 0.67), which we used in our analyses as a head size variable. To characterise limb length, we measured the left and right femurs and tibias, and then summing the mean femur (19.97 ± 2.06; mean ± SD) and tibia lengths (18.19 ± 1.19; mean ± SD) for every individual. Because during the experiment all animals received food ad libitum, we did not analyse body weight. Note that adding body weight measured at various stages during the experiment to our models never changed the results qualitatively (data not shown).

Behavioural assays
Behavioural assays took place on five consecutive days (between May 28th and June 2nd; we skipped an observation on 29 May because of thick cloud cover and rain). Individuals were tested in their home cages, as this approach more likely resembles how animals react in their natural home ranges (Beckmann and Biro 2013). We aimed to test individual behaviour in two, ecologically relevant situations. Thus, we placed a filter paper (4 cm × 22 cm) impregnated with the scent of a smooth snake (Coronella austriaca) inside the terraria of half of the lizards (Van Damme and Quick 2001; Amo et al. 2006). This snake is a saurophagous specialist, the main predator of I. cyreni (Martín 2015). The donor snake came from the same area where we had captured the lizards. Impregnated filter papers were obtained by placing them in the bottom of the snake's terrarium for 24 h before the experiments. Conversely, we placed a scentless clean filter paper of the same size into the terraria of the remaining half of lizards. Individuals were repeatedly tested in both situations (on different days). The order of exposure to the treatments was semi-randomised such that each individual received a sequence of treatments (Briffa 2013) (see Supplementary Text for more information). Filter papers were placed at 9.00 am (UTC + 2.00) and at the same time, shelters were removed. Tests for antipredator aggression took place the same day at 12.00 pm. We note that defensive responses of lizards are strongly dependent on temperature (Hertz et al. 1982;Crowley and Pietruszka 1983). As all behavioural data reported here were obtained on low wind, sunny days, and environmental variation between cages should have been minimal and random, we are certain that all experimental lizards could reach their preferred body temperature. Therefore, observed differences in the experiment should reflect individual differences in antipredator aggression rather than temperature effects. Following a similar standardised procedure in all trials, the same experimenter (GeH) caught the lizard by hand and recorded individual attempts of antipredator defensive bites while being held in the hand for 1 min. We treated antipredator aggression as a binomial variable, giving score 1 to individuals that attempt to bite the experimenter and score 0 to those lizards that did not. We note here that besides antipredator aggression, we also tested for social aggression of each individual towards an 'intruder' male placed into the focal (resident) male's home cage between 10.30 and 11.00 am. Results from these tests are not going to be discussed in this paper. However, we controlled for the potential effects of these tests on the antipredator aggression tests discussed here (see 'Statistical analyses'). Behavioural assessment was not blind regarding the test animals' identity; however, it does not pose a problem since the subjectivity in our methods was minimal.

Statistical analyses
To estimate repeatability for the binomial antipredator aggression data, we applied a generalised linear mixed model (GLMM) with a binomial distribution and logit link function on the pooled sample. In this model, antipredator aggression was the response variable. We also built separate GLMMs to examine the effect of treatments on repeatability. We used the rptR add-on package following the methods of Nakagawa and Schielzeth (2010). This method uses a multiplicative overdispersion GLMM with a logit link and using penalised quasi likelihood (PQL) estimation for repeatability on the original scale. Significance was estimated by randomisation tests. However, we report repeatabilities estimated on the underlying latent (link) scale as most original-scale repeatabilities are conditional for non-Gaussian data (Nakagawa and Schielzeth 2010). To control for variance explained by size differences between focal males, we added SVL as fixed effect to GLMMs otherwise similar to the previous ones and calculated enhanced agreement repetabilities (see Stoffel et al. 2017). Variance explained by the fixed effects are calculated by the variance on the link scale.
To test whether the treatments affected the mean antipredator aggression, we ran a binomial GLMM with logit link. In our model, antipredator aggression was the response variable, while predation treatment, SVL, relative head length, relative hind limb length and the two-way interactions of treatment and all morphology traits as fixed effects. 'Relative' length variables were residuals from trait-SVL linear regressions. We used them instead of the original raw variables to avoid multicollinearity. In a pilot GLMM, we tested the potential effects of preceding social aggression tests on antipredator aggression. Due to logistic constraints, we had to use intruder males (used only as stimuli) both smaller and larger than the focal individuals in the social aggression tests. Therefore, to test for the effects of social aggression tests, we added the SVL difference between resident and intruder males (integer variable, negative numbers indicating large intruder while positive numbers indicating small intruder males) as a correction variable; however, as this variable did not affect antipredator aggression (see Supplementary Table 1), we excluded this variable from our final model. Fixed effects were tested by Wald's chi-square tests, while random effects were estimated using likelihood ratio tests (LRT). The models were run with the lme4 package (Bates et al. 2015). We tested potential habituation effects by including the z-transformed (standardised to mean = 0, standard deviation = 1) order of trials (hereafter: time) both as a single-fixed effect and random slopes (i.e. the interaction with individual) in our mixed-effect models. Random-intercept (random effect: individual) and random slope (random effects: individual, individual × time) models were compared; we decided to leave the random slope term in the final model only if it improved the model fit. We report marginal and conditional R 2 estimations for our models based on the method of Nakagawa and Schielzeth (2013). All analyses were performed using R 4.0.2 (R Developmental Core Team 2020).

Behavioural repeatability
Randomisation tests indicated significant repeatability of antipredator aggression across all assays (Table 1). As randomisation test gives robust measures of statistical significance in the case of non-Gaussian data (Nakagawa and Schielzeth 2010), we consider antipredator aggression personality being present in all treatments and also in the pooled sample. Size differences explained a substantial amount of phenotypic variation (Table 1), reducing the relative contribution of individual differences. Nevertheless, there was still significant moderate-high repeatability present at the individual level regarding both treatments (Table 1). Since confidence intervals highly overlapped between the treatments, the difference between repeatability estimates cannot be considered as significant.

Behavioural type
The GLMM indicated significant effect of perceived predation risk on antipredator aggression (Table 2, Fig. 1): lizards were on average more defensive in the presence of predator cues. Furthermore, we found significant predator treatment effect and predator treatment × relative hind limb length interaction (see Table 2). To interpret the interaction, we ran two separate GLMMs for the presence and absence of predator treatments. We found a negative relationship between relative hind limb length and antipredator aggression in the presence (χ 2 = 3.92, df = 1, P = 0.047; Fig. 2a), but not in the absence of predator (χ 2 = 0.72, df = 1, P = 0.39; Fig. 2b). Moreover, SVL had a significant effect on antipredator aggression (Table 2, Fig. 3): larger males were more defensive than smaller conspecifics. Time effect was also significant (Table 2), individuals became more defensive by time; however, individual trends did not differ ( Table 2). The fixed effects explained 54%, while the whole model 84% of the total variance, which can be seen as good explanatory power for behavioural variables.

Discussion
Antipredator aggression can be seen as prey's last resort when escape from the predators failed and prey is being caught. This behaviour is obviously affecting prey's fitness; however, it is severely understudied in natural populations compared to social aggression (but see Cury de Barros et al. 2010;Markó et al. 2013; Baxter-Gilbert and Riley 2018). Here, by conducting a manipulative experiment with adult male I. cyreni, we asked whether animal personality (consistent between-individual difference) is present in antipredator aggression and also tested whether ecological conditions (presence vs. absence of perceived predation risk before the trials) and individual state (body size, head size, limb length) affect the expression of this behaviour. We found strong indication for antipredator aggression personality, both in the presence and absence of perceived predation risk before the mimicked attack, emphasising that this trait is an important aspect of the behavioural repertoire of I. cyreni and might be under natural selection. We found lizards exposed to chemical stimuli from their snake predator, and larger males in general to show increased antipredator aggression. We also found an environment-dependent individual state effect: individuals with relatively short legs were more defensive compared to their long-legged conspecifics, but only after exposure to snake predator cues. Individuals became more defensive by  time; however, habituation is assumed to reduce unnecessary antipredator responses (i.e. probability of bites) (Rodríguez-Prieto et al. 2010Vincze et al. 2016); thus, our pattern suggests a reverse response. Such response is known as sensitization: an internal mechanism that intensifies behavioural response to constant stimulation (Bee 2001;Martin and Réale 2008;Stamps et al. 2012;Osborn and Briffa 2017). Below, we discuss first the relevance and consequences of animal personality being present in antipredator aggression, and second, the effects of external environment and internal state on this behaviour.

Behavioural repeatability
Antipredator aggression was significantly repeatable within individuals. We are aware that the wide confidence intervals of our repeatability estimates indicate that our sample size is somewhat low, especially in terms of repeated measures. However, we note that all of the enhanced agreement repeatability estimate presented here can be seen as moderate-high (see Bell et al. 2009); hence, we think that the pattern we present regarding the presence of behavioural repeatability (a statistical test for animal personality) is robust. While body size variation explained a substantial proportion of behavioural variation (larger lizards showing higher antipredator aggression), there was still significant between-individual behavioural difference present in the examined population after controlling for body size. Thus, the most important question our result raises is what are the costs of antipredator aggression, i.e. why do some individuals not defend themselves actively in a supposedly life-threatening situation? Unarguably, being captured by a predator rarely holds any positive outcomes for the prey regardless its behaviour; yet, if survival is on the stake, a defensive attack may initially provide some advantage (Blanchard et al. 2008). However, there are other passive alternative antipredator strategies in many animals. In particular, tonic immobility, when the prey remains motionless when captured by a predator, or in other situations of extreme fear  is also widespread (Gallup and Rager 1996;Miyatake et al. 2008). It has been argued that tonic immobility can increase prey survival because the lack of movement makes the predator believe it already killed the prey and stop the attack, releasing the prey to start ingestion, allowing the prey to escape (Thompson et al. 1981). Even though we did not quantify 'struggling', we observed that lizards either vigorously struggled and bit or they were passive. Therefore, the absence of defensive behaviour in some of the lizards in our experiment might be considered as a sign of tonic immobility. Tonic immobility strategy has been previously considered in personality studies (Edelaar et al. 2012;Santos et al. 2015;Horváth et al. 2019). Our present findings hint that between-individual differences in antipredatory last resort behavioural traits might be interpreted along a tonic immobility-antipredator aggression continuum. This hypothesis warrants further studies.

Behavioural type
In line with our expectations, lizards exposed to chemical cues from their snake predator were more defensive in general than males from the control treatment. Previous experimental results show that olfactory predator cues are able to induce strong antipredator behavioural responses in various amphibian and reptile taxa (e.g. Rana dalmatina tadpoles, see Hettyey et al. 2015;Urszán et al. 2015aUrszán et al. , b, 2018; R. latastei tadpoles, see Scribano et al. 2020; Anguis fragilis, see Cabido et al. 2004; Podarcis muralis, see Amo et al. 2006; P. tliguerta, P. sicula and Lacerta bedriagae, see Van Damme and Quick 2001); nevertheless, studies rarely test the effect of these treatments on antipredator attacks. Similar to our present findings, López and Martín (2001) reported the fossorial amphisbaenian Blanus cinereus to show increased willingness to bite towards predator-scent (Coronella girondica and Scolopendra sp.) impregnated swabs than control stimuli. These results suggest that presence of chemical cues from predators connotes a stressful situation which leads to increased antipredator aggression. However, as predator presence vs. absence was changed daily in our experiment, the detected pattern here could be seen as short-term responses, representing a form of activational behavioural plasticity (see Snell-Rood 2013). Furthermore, partly in line with our expectations, males with shorter hind legs showed increased antipredator aggression, but only after exposure to snake predator cues. Limb length was showed to be strongly positively correlated with sprint speed (i.e. locomotor performance; Bauwens et al. 1995;Vanhooydonck et al. 2001), which is an important component of life history trade-offs and suitable proxy of individual quality (Garland 1984;Le Galliard et al. 2004;Irschick et al. 2008;Husak et al. 2016;Winchell et al. 2018). Hence, males with longer hind legs could be seen as higher-, while short-legged males as lower-state individuals. A comparable pattern was reported by Hertz et al. (1982), who found a negative relationship between body temperature (high body temperature increases locomotor performance) and defensiveness in the lizards Agama savignyi and A. pallida. In line with these previous results, it is plausible that longer-legged males have higher chance to outrun and out-manoeuver potential predators, thus rely less on defensive attacks to escape. Furthermore, as predator cues in our experiment indicate the presence of an ambush snake predator, an active escape might trigger a quick predator attack. As the ability to escape highly depends on sprint speed, short-legged males are at a disadvantage compared to long-legged conspecifics, thus, holding their ground and confronting the predator actively might be the best strategy for them. Increased antipredator aggression under constrained locomotor performance in tegu lizards (Tupinambis merianae) was shown previously by Cury de Barros et al. (2010). Also in line with the predictions, larger (õ lder) males were more defensive in general than smaller (ỹ ounger) conspecifics. In a previous study, large I. cyreni individuals were found to be more prone to flight (i.e. more riskaverse) than smaller ones (Martín and López 2003). This pattern corresponds to state-dependent safety principle (Luttbeg and Sih 2010): although larger lizards are more conspicuous to predators and thus they flee earlier than smaller ones (Martín and López 2003;Baxter-Gilbert and Riley 2018), when cornered, they can confront the predator with higher chance than smaller conspecifics due to their more powerful bites (Gvozdik and Van Damme 2003;Cury de Barros et al. 2010).

Conclusion
Taken together, we found antipredator aggression showing consistent between-individual differences in the studied population of I. cyreni males, suggesting that this behavioural trait (i) should be considered in animal personality studies and (ii) might respond to natural selection assuming that the betweenindividual variation has a genetic component. We are aware that significant repeatability during a short period cannot be seen as a sign of animal personality in the classical sense, as stable individual differences throughout the life or several years. However, we think that such differences are still informative in the given ecological context, i.e. the short and highly synchronised mating season (e.g. Horváth et al. 2016Horváth et al. , 2017. Antipredator aggression was positively associated with the presence of predator cues and body size (regardless of the presence of predator cues), while relatively short-legged males showed increased defensiveness only when predator cues were present, probably as an antipredatory strategy compensating limited escape speed. Results of our experiment indicate that differences in inherently stable state variables (size, hind limb length) affect individual antipredator aggression, further emphasising the ecological and evolutionary relevance of this behaviour. Nevertheless, as other inherently labile traits, such as levels of testosterone are proved to substantially affect intraspecific social aggression in a various of taxa (Veiga et al. 1998;Weiss and Moore 2004;O'Connor et al. 2014;Szász et al. 2019), further studies are needed to reveal potential background mechanisms behind variation in antipredator aggression. Linking antipredator-and social aggression in a behavioural syndrome framework (Sih et al. 2004;Dingemanse et al. 2012;Dochtermann and Dingemanse 2013), as well as understanding the inheritance of antipredator aggression also warrants future studies.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Ethical approval All applicable international, national, and institutional guidelines for the care and use of animals were followed. No lizard was injured during the handling related to the study and after completion of the experiments, animals were returned in good health at the exact site of capture. The study was performed under license (permit number: 10/ 056780.9/16) from the Environmental Agency of Madrid Government ("Consejería de Medio Ambiente de la Comunidad de Madrid", Spain), and in accordance with the national animal welfare standards and protocols supervised by the Bioethical Committee of the Spanish Research Council (CSIC).
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/.