Roll with the fear: environment and state dependence of pill bug (Armadillidium vulgare) personalities

Most studies on animal personality evaluate individual mean behaviour to describe individual behavioural strategy, while often neglecting behavioural variability on the within-individual level. However, within-individual behavioural plasticity (variation induced by environment) and within-individual residual variation (regulatory behavioural precision) are recognized as biologically valid components of individual behaviour, but the evolutionary ecology of these components is still less understood. Here, we tested whether behaviour of common pill bugs (Armadillidium vulgare) differs on the among- and within-individual level and whether it is affected by various individual specific state-related traits (sex, size and Wolbachia infection). To this aim, we assayed risk-taking in familiar vs. unfamiliar environments 30 times along 38 days and applied double modelling statistical technique to handle the complex hierarchical structure for both individual-specific trait means and variances. We found that there are significant among-individual differences not only in mean risk-taking behaviour but also in environment- and time-induced behavioural plasticity and residual variation. Wolbachia-infected individuals took less risk than healthy conspecifics; in addition, individuals became more risk-averse with time. Residual variation decreased with time, and individuals expressed higher residual variation in the unfamiliar environment. Further, sensitization was stronger in females and in larger individuals in general. Our results suggest that among-individual variation, behavioural plasticity and residual variation are all (i) biologically relevant components of an individual’s behavioural strategy and (ii) responsive to changes in environment or labile state variables. We propose pill bugs as promising models for personality research due to the relative ease of getting repeated behavioural measurements.


Introduction
Behaviour is one of the most flexible traits of animals (West-Eberhard 2003), yet some level of repeatability in behaviour across time and ecological situations (i.e. animal personality) exists (Bell et al. 2009;. Intuitively, the presence of non-random among-individual behavioural variation should constrain behavioural plasticity (Niemelä et al. 2013). This is true to a certain extent, but individuals still preserve the ability to alter their behaviour in response to changing environment, while their behaviour relative to each other remains different (Biro et al. 2010;Dingemanse et al. 2010;Briffa et al. 2013;Mathot and Dingemanse 2014). Further, it seems that in addition to non-random variation in mean behaviour, individuals can also show variation in their reaction to environmental change (within-individual behavioural plasticity) Westneat et al. 2011;Dingemanse and Wolf 2013;Mitchell and Biro 2017). Finally, biological validity and importance of withinindividual behavioural variation not induced by environmental change, or in other words, the 'rigidity' of an individual's behaviour type in a certain environment (within-individual residual variation), were recognized recently (Stamps et al. 2012;Biro and Adriaenssens 2013;Briffa 2013;Briffa et al. 2013). Hence, within-individual behavioural plasticity (hereafter: behavioural plasticity) and within-individual residual variation (hereafter: residual variation) should be considered as potentially independent components of individual behavioural strategy Kralj-Fišer and Schneider 2012;Briffa 2013;Dingemanse and Wolf 2013;Westneat et al. 2013Westneat et al. , 2015Mitchell et al. 2016). However, background mechanisms affecting emergence of individual variation in behavioural plasticity and residual variation are less understood. In addition, it is still not entirely clear whether these components can evolve independently or individual differences in within-individual behavioural variation are related to personality (Niemelä et al. 2013;Mathot and Dingemanse 2014;Stamps 2016).
Recently, a growing body of studies suggests that even short-term variation in ecological conditions and inherently labile state-linked traits could create stable differences in behavioural strategies Lichtenstein et al. 2016;Horváth et al. 2017). Parasites are among the most important environmental factors known to create stable behavioural differences (Barber and Dingemanse 2010;Kortet et al. 2010;Poulin 2013). Direct or indirect negative effects of infections (e.g. low body condition) may result in differences in individual state, which may eventuate the emergence of individual behavioural strategies in order to cope with these disadvantages (DiRienzo et al. 2015Horváth et al. 2016). In arthropods, Wolbachia are important parasitic bacteria (Hilgenboecker et al. 2008;Werren et al. 2008) with remarkable effects on hosts' physiology, including parthenogenesis, reproductive incompatibility, feminization and male killing (see Werren and Windsor 2000;Werren et al. 2008;Le Clec'h et al. 2012, 2013. Although behavioural impact of Wolbachia is less documented, it seems that infection generally affects mating behaviour of males (see Ming et al. 2015;Moreau et al. 2001;Vala et al. 2004;Zhao et al. 2013). Also, in a parasitic wasp (Leptipolinia heterotoma), it was found that Wolbachia reduces locomotor activity of both sexes (Fleury et al. 2000). Generally speaking, Wolbachia is expected either to decrease behavioural activity by impairing physiological performance of the hosts or to increase it by host manipulation or inducing some sort of terminal investment (Sicard et al. 2010;Chevalier et al. 2011).
Conglobation is a special type of tonic immobility and is a common defensive behaviour for various arthropod taxa (Tuf et al. 2015). Conglobation can be used as a proxy of risktaking behaviour (Carter et al. 2012;Beckmann and Biro 2013). Species of pill bugs (family: Armadillidiidae, order: Isopoda, subphylum: Crustacea) are capable to roll their body into an uninterrupted sphere, hiding their vulnerable posterior appendages (uropods), legs and antennae. Conglobation protects common pill bug (Armadillidium vulgare) from most invertebrate and small-sized vertebrate predators, while larger vertebrate predators easily overlook small, immobile prey (Matsuno and Moriyama 2012;Tuf et al. 2015). Considering how easy it is to measure conglobation time, pill bugs might be excellent models for animal personality research, where the current statistical approaches are data hungry and the necessary number of within-individual repeated measures is challenging to reach with most species .
In the present paper, we studied risk-taking of adult A. vulgare by performing 30 repeated behavioural assays in two different environments (familiar vs. unfamiliar). We studied the effects of environment and various state variables (Wolbachia infection, body size, sex) on individual behaviour on different levels. First, we were interested whether individuals differed in mean risk-taking (i.e. among-individual variation), in their reaction norms (i.e. behavioural plasticity) and residual variation. Second, we tested whether these components are affected by environment and individual state. We expected lower risk-taking in the unfamiliar environment. We had no prediction regarding the effects of sex or size. Regarding residual variation, theory predicts increased within-individual variability under predation risk (Hugie 2003) as prey animals may reduce the probability of capture by predators by displaying unpredictable behaviour (Humphries and Driver 1970;Jones et al. 2011). However, empirical data are somewhat contradictory (see Briffa 2013;Velasque and Briffa 2016;Urszán et al. 2018). Displaying the highest residual variation may not be the best antipredator strategy; also, level of behavioural rigidity may depend on prevailing environmental conditions (Richardson et al. 2018), but more importantly, on development (see Bierbach et al. 2017). Thus, although we did not form a directional hypotheses regarding how different environments would affect residual variation of risk-taking in pill bugs, we expect that the level of residual variation indeed differs between familiar vs. unfamiliar environments. Further, we expected both environmentally induced plasticity in the form of lower risk-taking in the unfamiliar, potentially dangerous environment and plasticity along time in the form of habituation to the laboratory conditions. The latter was expected to be stronger in the unfamiliar environment.

Study animals
relatively dense vegetation and high cover of leaf litter. We searched for animals under leaves and decaying wood at different sites. Only one individual was collected at a certain spot, and at least 50 m was left between spots in order to reduce the chance of sampling individuals from the same family. Individuals were transported to the facilities of Eötvös Loránd University, where they were housed individually in white opaque plastic boxes (15.5 cm × 11 cm × 12 cm, length, width, height, respectively) with a 1-cm-deep substrate consisting of mixture of coconut fibre and soil. Humidity was maintained by spraying the substrate with distilled water twice a day, and fresh carrot was provided as food weekly. The experimental animals have always eaten from the food, but never fully consumed it; hence, the experimental animals were fed ad libitum. We provided 12-h light period per day during the experiment. Dim light was provided by Repti Glo 2.0 Full Spectrum Terrarium Lamps (Exo Terra, Rolf V. Hagen Inc., Holm, Germany), which do not emit considerable heat but mimics the full spectrum natural light. At the end of the experiment, we sexed the individuals and measured the body size (diameter of conglobated individuals; to the nearest 0.01 mm) and specimens were conserved in 96% ethanol for Wolbachia screening (see below). All individuals were screened for Wolbachia (see Supplementary Material for details).

Behavioural assays
Individual behaviour was evaluated 30 times during 38 days between 12 May and 19 June. Days without measurements were distributed randomly. Each individual's behaviour was evaluated in two different environments (familiar vs. unfamiliar) daily, resulting in 60 repeats per individual. Behavioural assays were carried out between 9.00 and 12.00 a.m. (UTC + 02.00), and 1-h break was provided between the two assays. The order of individuals within and between environments was randomized daily.
Risk-taking was estimated by latency to restart activity (time spent immobile in conglobation) after a simulated predator attack. Animals were removed from their home boxes, after which the experimenter gently squeezed the animal to trigger conglobation, then dropped it to the surface depending on treatment from a standard (10 cm) height. This treatment is similar to manipulation by larger vertebrate predators (e.g. birds and lizards) (see Tuf et al. 2015). In the familiar environment treatment, individuals were elevated from and dropped back to their home boxes, while in the unfamiliar environment treatment, individuals were dropped to a white plastic sheet illuminated directly (with 40 × G4 Halogen Light Bulb, 10 W, 12 V). Animals were considered to restart activity when they fully stretched their body and started to move their legs in an attempt to escape. If an individual did not restart activity in 15 min, the assay was stopped and the individual was assigned maximum score (900 s). This happened only in seven cases, including four individuals; hence, we used the maximum score in the subsequent analyses. For unknown reasons, more than half of the collected individuals died soon after being transported to the laboratory and one individual had to be removed from the analyses because of its extremely outlying size. Thus, we used data (60 risk-taking measurements) from 25 individuals, 11 males and 14 females.

Statistical analysis
Latency data were log-transformed to achieve normal distribution of model residuals. To be able to fit reaction norms, environment was treated as a continuous measure by assigning − 1 to the 'familiar' situation and 1 to the 'unfamiliar.' Continuous variables were centred for the analyses by bringing them to scale with a zero mean and unit variance. Dummy variables were created for the categorical variables to use them in the Bayesian modelling (see below). To describe the hierarchically structured behavioural data, we relied on a framework based on linear mixed modelling (LMMs) according to the following equation: where Y ijk is the latency for individual i measured in day j and at the kth observation within a day, β 0 is the population mean latency, β 1 -β 6 are the mean level parameters that describe the effect of covariates x 1 -x 6 (environment, day, sex, size, infection status and order within a day, respectively), ind 0i is the random effect term capturing the deviation of individual i from the population mean, while day 0j is the random effect term depicting the deviation caused by day-specific effects. The model also considers random slopes ind 1i and ind 2i to deal with among-individual differences in plasticity with respect to novel environment and with respect to day (habituation), respectively. The model requires the following assumptions: Accordingly, the random terms are assumed to be normally distributed with a mean of zero and a respective variance, i.e. σ 2 ind 0 (among-individual variance in mean latency), σ 2 day 0 (among-day variance in mean latency), σ 2 ind 1 (among-individual variance in plasticity) and σ 2 ind 2 (among-individual variance in habituation). The error term is assumed to rely on a common residual variance σ 2 res (within-individual variance). We defined the above starting model based on a list of considerations. First, all of the fixed predictors are biologically relevant and can be hypothesized to affect the focal behavioural trait either at the among-individual or at the within-individual level. Second, the random effects were chosen to describe the hierarchical structure of the data arising from the design of the study, as well as to accommodate our main predictions concerning individual differences in plasticity and habituation. Accordingly, with this model, we could test whether (i) risk-taking behaviour was linked to individual characteristics like sex, size or health status and whether (ii) environment-or timeinduced (i.e. habituation) behavioural plasticity was present. For simplicity and to avoid too many parameters to be estimated, we did not define interaction terms among the fixed predictors. Similarly, we have not considered covariance between random slopes and intercepts. To verify that the random part of the above model is appropriate for further evaluation, we defined alternative models in the lme4 R package (Bates et al. 2015) and examined their goodness of fit relative to the model described in Eq. (1) by using likelihood ratio test. These investigations revealed that both random intercept terms are significant (P < 0.001 for both ind 0i and day 0j ), that random intercept and slope models offer better fit to the data than the random intercept only models (P = 0.018 for ind 1i and P < 0.001 for ind 2i ) and that allowing correlation between random intercepts and slopes does not imply further improvement statistically (P = 0.052 for the correlation between ind 0i and ind 1i , P = 0.373 for the correlation between ind 0i and ind 2i ).
Given Eq. (6), the above model assumes that residuals have a homogeneous variance, which is a strict assumption and does not accommodate a possibility for testing amongindividual differences in predictability, which corresponds to one of the main hypotheses of this study. To account for heterogeneous residual structure, we adopted an approach based on double hierarchical general linear models, which allows fitting the main and the dispersion parts of an LMM within the same statistical framework Nelder 1996, 2006), with the latter capturing the essence of residual variation (Westneat et al. 2013;Cleasby et al. 2015;Mitchell and Biro 2017). Accordingly, we kept the above model and added another model for the standard deviation (SD) of risk taking, as follows. First, we assumed heterogeneous residual variance by replacing Eq. (6) with which permits distinct residuals for each observation that can be further described as In this equation, γ 0 is the mean log residual SD, γ 1 -γ 6 are parameters describing the effect of covariates (environment, day, sex, size, infection status and order within a day, respectively) on predictability. The random terms ind σ0i and day σ0j reflect individual-and day-specific deviations, respectively, from the population-specific mean log SD. For these random terms, we assumed that yielding that different individuals and days can be characterized by different residual variation, and the among-individual and among-day variance of residual variation can be estimated as σ 2 σind and σ 2 σday . We were also interested in whether individual-specific behavioural plasticity and habituation can be linked to individual characteristics. To this end, we further extended the modelling framework by adding linear regressions that describe individual-specific plasticity and habituation. Therefore, we replaced Eqs. (4) and (5) with in which each individual can be described by a specific mean plasticity and habituation value depending on their individual characteristics as specified by Here, δ 1 -δ 3 and φ 1 -φ 3 stand for parameters that link the main individual-specific attributes (x 3 -x 5 , sex, size and infection status) to plasticity and habituation, respectively. Note that the regressions are forced through the origin (intercept is zero); thus, the individual-specific mean is fixed to be 0.
The parameters of the above models are estimated iteratively and depend on one another allowing the test of our main predictions in a single statistical framework (see graphical representation of the whole model structure in electronic supplementary material). Such inferences typically require Bayesian approaches based on Markov chain Monte Carlo (MCMC) processes (Gelman et al. 2004). For this purpose, we applied procedures available in program JAGS (Plummer 2003) which we controlled from within the R statistical environment (R Developmental Core Team 2018) using the package rjags (Plummer 2014). For each model, we defined three MCMC chains from different starting values and with 10,000 iterations of burn in period before sampling the posterior distribution. The posterior sample relied on the subsequent 100,000 iterations that were combined across chains. Before interpreting the results, we applied the conventional model diagnostic procedures for each chain to verify convergence and the lack of autocorrelation (Gelman et al. 2004). Bayesian methods require priors to be defined for each parameter estimated. We set these values to have a minimal influence on the posterior distribution, as we had no previous knowledge about them (uninformative, flat priors). The means of the posterior distributions were used for further interpretation by also considering the associated 95% credible intervals (95% CrI). If the generated 95% CrI of the posterior distribution for a parameter did not include zero, the parameter was considered to have a significant effect. Model codes are available in the electronic supplementary material.

Discussion
Here, we demonstrated substantial differences in risk-taking of A. vulgare at several hierarchical level of behavioural variation. Pill bugs showed significant between-individual differences not just in mean risk-taking (i.e. among-individual variation), but in the degree to which they adjust their behaviour to previous environmental conditions (i.e. behavioural plasticity), and how consistently express their behaviour in any given environment (i.e. residual variation). These patterns add to prior studies indicating that both behavioural plasticity and residual variation can be seen as potentially independent components of animal behavioural variation Stamps et al. 2012;Biro and Adriaenssens 2013;Briffa 2013;Briffa et al. 2013;Dingemanse and Wolf 2013;Stamps 2016;Chang et al. 2017;Guayasamin et al. 2017;Lichtenstein et al. 2017). Further, in line with several prior studies, we found that individual variation in residual variation was influenced by individual state or environmental differences (Briffa 2013;Bridger et al. 2015;Westneat et al. 2015). Our results provide empirical support to the notion that inter-individual differences in within-individual behavioural variation may be the outcome of adaptive processes rather than reflecting non-functional variation (Biro and Adriaenssens 2013). Here, we discuss how these effects on Fig. 1 Differences in risk-taking over time (habituation) in common pill bug (Armadillidium vulgare). Note that risk-taking is a latency variable, i.e. lower values represent higher risk-taking among-individual variation, residual variation and behavioural plasticity fit to the existing theories.

Among-individual variation
A. vulgare individuals infected with Wolbachia took less risk than their uninfected conspecifics. It has been suggested that infection with mild effects on the hosts' fitness is coupled with higher rates of behavioural activity, as infection often increases energetic needs (Lafferty and Morris 1996;García-Longoria et al. 2014;Gyuris et al. 2016). On the other hand, parasites with severe negative effects on their hosts' fitness are more likely to reduce behavioural activity, due to heavily reduced state (e.g. low body condition; Ferguson et al. 2011;Hammond-Tooke et al. 2012;Poulin 2013), which may eventuate the emergence of individual behavioural strategies in order to cope with these disadvantages (Barber and Dingemanse 2010;Kortet et al. 2010;Horváth et al. 2016).
Pathological effects of Wolbachia infection are rarely studied; however, it is known that the parasite is able to avoid or even manipulate the host's immune system and affect senescence processes directly (Braquart-Varnier et al. 2008;Le Clec'h et al. 2012). In A. vulgare, different strains have differently severe effects on the host. For instance, wVulC is a widespread and invasive feminizing strain, inducing low haemocyte level and intense septicemia, reducing the host lifespan considerably (Braquart-Varnier et al. 2008;Sicard et al. 2010;Chevalier et al. 2011;Le Clec'h et al. 2012, 2013. On the other hand, wDil has no proven direct effects on A. vulgare fitness, but does induce cytoplasmic incompatibility which may generate indirect costs (Sicard et al. 2010;Valette et al. 2013). Here, lowered risk-taking in infected pill bugs could be a result of compensation for lowered body condition and physiological performance. However, considering the low prevalence of Wolbachia infection in or sample, this finding should be interpreted with caution (Bell et al. 2010). Thus, our conclusions regarding strength and true mechanisms behind this pattern here are rather tentative. Nevertheless, we believe that this finding at least warrants more targeted research on the potential effect of Wolbachia infection on hosts' behaviour.

Residual variation
Residual variation in behaviour was found to decrease across days, in correspondence with human psychology literature as well as observations on various vertebrate and invertebrate taxa showing residual variation decreasing with increasing experience (Stamps et al. 2012;Stamps and Krishnan 2014). On the contrary, in a recent study performed on guppies (Poecilia reticulata), Mitchell et al. (2016) report no change in residual variation across a timespan similar to ours. It is known that adult individuals may acclimate by their novel environment quicker and express lower residual variability within shorter periods (Biro 2012); nevertheless, ontogenetic effects are rather implausible in our case. A more likely possibility is that residual variation in risk-taking decreased due to continued acclimation to our experimental procedure (Biro and Adriaenssens 2013;Mitchell et al. 2016).
Pill bugs express significantly higher residual variation in the unfamiliar than in the familiar environment. This finding is consistent with recent reports from both vertebrate and invertebrate taxa Stamps et al. 2012;Briffa 2013;Nakayama et al. 2016), suggesting that potentially risky environments decrease predictability of behaviour (but see Urszán et al. 2018). Thus, high residual variation Fig. 2 Individual behavioural reaction norms across a time (sensitization) and b environments in common pill bug (Armadillidium vulgare). Note that risk-taking is a latency variable, i.e. lower values represent higher risk-taking could be seen as an antipredator response (Biro and Adriaenssens 2013;Briffa et al. 2013). Predation pressure substantially affects behavioural actions of individuals in order to avoid potentially risky encounters (Bell and Sih 2007;Kortet et al. 2010;Luttbeg and Sih 2010;Engqvist et al. 2015;Sih et al. 2015). However, taking that substantial differences in residual variation were present irrespective of environmental factors, other most likely internal factors should also affect residual variation (Sih et al. 2004(Sih et al. , 2015Briffa 2013;Bierbach et al. 2017). Again, we cannot be sure regarding exact background mechanisms of this pattern. It is known that in A. vulgare, mating is linked to moulting cycle, during which individuals are more vulnerable to predators (Beauché and Richard 2013). It is likely that high residual variation in emergence from conglobation during the reproductive season helps secure survival and thus future reproductive success of the individuals.  4 Association between behavioural plasticity of risk-taking over time (sensitization) and size in common pill bug (Armadillidium vulgare). Individual sensitization is represented by the slope of the individual behavioural reaction norm in response to time Fig. 3 Differences in risk-taking residual variation over time in common pill bug (Armadillidium vulgare). Estimates were obtained from the statistical model (see Table 1)

Behavioural plasticity
We found significant decrease of risk-taking with time. As habituation is assumed to reduce unnecessary antipredator responses (i.e. length of conglobation) (Rodríguez-Prieto et al. 2010Vincze et al. 2016), our pattern rather suggests a reverse response that is known as sensitization; an internal mechanism intensifies behavioural response to constant stimulation (Bee 2001;Martin and Réale 2008;Stamps et al. 2012;Osborn and Briffa 2017). Theory predicts that sensitization eventually will fade and habituation becomes the main pattern of behavioural change over time, but empirical studies provide limited support for this (Bee 2001;Osborn and Briffa 2017).
Our data indicate no sign of habituation (i.e. lowered risktaking) during 30 days of experiment in either the familiar or in the unfamiliar environment. High level of sensitization is assumed to be linked to high stimulus rate and intensity (Bee 2001), and since our assays were conducted mostly on a daily basis, the result is somewhat in line with this prediction. Individual-specific behavioural plasticity was not affected by individual state. On the other hand, among-individual variance in slope was influenced by sex and size, i.e. females and larger specimens (irrespective of sex) became shier with time. Empirical data suggest that individual-specific differences in behavioural plasticity are the outcome of variability of inherently stable or labile state variables (Luttbeg and Sih 2010;Wolf and Weissing 2010;Mathot and Dingemanse 2014;Araya-Ajoy and Dingemanse 2017;Mitchell and Biro 2017). However, based on current correlative data, we cannot reconstruct the exact biological mechanism in the background of this pattern, especially if we take into account that differences in state are affected by both genetic and environmental variation (Mathot and Dingemanse 2014). The most plausible explanation is that high behavioural plasticity likely secures future reproductive success of females and large pill bugs. In general, the patterns reported in this subsection have also added to the growing number of data indicating that despite behaving in a consistent way, individuals still maintain the capability to adjust their behaviour to changing environmental conditions (i.e. being behaviourally plastic) (DeWitt et al. 1998;Dingemanse and Wolf 2013;Mathot and Dingemanse 2014).

Conclusions
Taken together, we found components of behavioural variation (intercept, slope and residuals) to exhibit amongindividual variation and to be sensitive to different variables related to individual state. These results suggest that all three components are integral parts of an individual's behavioural strategy and that individuals are indeed plastic upon environmental challenge and within-population behavioural variation can be at least partly explained by variation in fixed and labile state variables. We recommend studying behavioural variation in an integrative approach and along longer observational periods, as animal personality sensu lato, or in other words, individual behavioural strategy seems to be indeed more than just variation in mean behaviour.