Cross-Context Responses to Novelty in Rural and Urban Small Mammals

The Anthropocene is the era of urbanization. The accelerating expansion of cities occurs at the expense of natural reservoirs of biodiversity and presents animals with challenges for which their evolutionary past might not have prepared them. Cognitive and behavioral adjustments to novelty could promote animals’ persistence under these altered conditions. We investigated the structure of, and covariance between, different aspects of responses to novelty in rural and urban small mammals of two non-commensal rodent species. We ran replicated experiments testing responses to three novelty types (object, food, or space) of 47 individual common voles (Microtus arvalis) and 41 individual striped field mice (Apodemus agrarius). We found partial support for the hypothesis that responses to novelty are structured, clustering (i) speed of responses, (ii) intensity of responses, and (iii) responses to food into separate dimensions. Rural and urban small mammals did not differ in most responses to novelty, suggesting that urban habitats do not reduce neophobia in these species. Further studies investigating whether comparable response patters are found throughout different stages of colonization, and along synurbanization processes of different duration, will help illuminate the dynamics of animals’ cognitive adjustments to urban life.


INTRODUCTION
Urbanization is a defining process of the Anthropocene, stemming from the conversion of natural areas into land covered by impervious artificial surfaces (e.g., Møller, 2009;Biermann et al., 2016;Munshi-South et al., 2016). Urbanized areas are characterized by elevated habitat fragmentation, strong edge effects, high levels of disturbance and environmental heterogeneity (e.g., Van Ham et al., 2013). As the global human population rises, the expansion of urbanized areas is accelerating at the expense of natural reservoirs of biodiversity (e.g., McKinney, 2002;Grimm et al., 2008;Seto et al., 2012;Gao and O'Neill, 2020). The loss of vulnerable species intensifies biological homogenization, and results in separate city landscapes that are more similar to each other than to the natural ecosystems that they replace (e.g., McKinney, 2002;Barnosky et al., 2012;Groffman et al., 2014;Alberti, 2015;Alberti et al., 2017).
These urbanized habitats present animals with specific challenges for which their evolutionary past might not have prepared them (e.g., McDonnell and Hahs, 2015). For instance, urban dwellers have to cope with novel or altered levels of human disturbance, habitat fragmentation, resource availability, predators, and competitor and parasite communities (e.g., Sih et al., 2011). Wildlife responses to novelty and urbanization are thus a key focus of current behavioral and cognitive research (e.g., Crane et al., 2020;Mettke-Hofmann et al., 2020;Sol et al., 2020;Lee and Thornton, 2021). Since behavior and cognition largely determine how individuals interact with their surroundings, behavioral and cognitive adaptations are expected to play a major role in coping with human-induced rapid environmental change (HIREC, sensu Sih et al., 2011) (e.g., Greggor et al., 2014;Barrett et al., 2019;Goumas et al., 2020).
"Neophobia is an ecologically relevant fear behavior that arises through a cognitive assessment of novel stimuli" (Greggor et al., 2015, p. 82). Responses to novelty are classically categorized into responses to novel objects, foods, and space (e.g., Greggor et al., 2015). Responses to novelty are often inferred to be interchangeable, and responses in one context are frequently considered predictive of responses in the others (e.g., Greenberg, 2003). However, it is still not well understood whether responses to different novelty dimensions correlate across contexts, and whether they are indicative of the same underlying mechanism (e.g., Greggor et al., 2015;Crane et al., 2020). Different selective pressures might have shaped responses to novelty dimensions independently, as well as the specific type of responses (e.g., Greenberg and Mettke-Hofmann, 2001;Mettke-Hofmann, 2014;Greggor et al., 2015). High predation risk might have selected for heightened fear responses to unfamiliar objects and stimuli, while high dispersal propensity might come with elevated exploration tendency of novel spaces (e.g., Ferrari et al., 2007;Brown et al., 2014;Mettke-Hofmann, 2014). Neophobic responses toward food could have evolved in species with frequent encounters with dangerous or poisoned food (e.g., Galef, 1993;Thornton, 2008). Classifying the types of responses themselves is not always straightforward, because novelty elicits two opposite tendencies: aversion (neophobia) but also attraction (neophilia), which are often seen as a continuum (e.g., Greenberg, 2003). In many bird species neophobia and neophilia are governed by independent motivations and can vary independently (e.g., Greenberg and Mettke-Hofmann, 2001;Mettke-Hofmann et al., 2002;Mettke-Hofmann, 2014). Neophilia facilitates information gathering, whereas neophobia mainly protects the animal from unknown potential dangers by heightened caution (e.g., Greenberg, 2003;Griffin and Guez, 2014).
Responses to novelty are often found to vary consistently among individuals (e.g., Herborn et al., 2010;Yuen et al., 2016; but see e.g., Brown et al., 2013). They are commonly used as a marker of animal personality or temperament, since they are considered a stable response to challenges or risks across times or situations (e.g., Réale et al., 2007;Greggor et al., 2016a). In freeliving animals, responses to novelty are associated with reactions to other novel situations, for example innovation propensity in a novel foraging task or dietary conservatism toward new food types (e.g., Webster and Lefebvre, 2001;Thomas et al., 2004). Since positive responses to novelty facilitate interacting with aspects of the environment, they provide opportunities for learning (e.g., Herborn et al., 2010). Therefore, the propensity and intensity of novelty responses should matter in the process of adjustment to altered environmental conditions (e.g., Sol et al., 2013;McDonnell and Hahs, 2015). Neophobia has been proposed as a major predictor for the successful establishment in novel or altered habitats, and indeed several studies on urban birds report lower neophobic and/or higher neophilic responses in urban dwellers compared to conspecifics living under more natural conditions (e.g., Sol et al., 2011Sol et al., , 2013Ducatez et al., 2017; but see e.g., Audet et al., 2016;Greggor et al., 2016b). In contrast, we know very little about the role of neophobia in coping with urban environments in species more limited in dispersal. These species might be precluded from easily moving away in response to intense human disturbance. Also, the colonization of humanaltered environments could be hindered by pronounced aversion to novel stimuli in species with high predation pressure, in which individuals' avoidance of dangers, and hence survival, strongly impacts ultimate fitness (e.g., Lima and Dill, 1990;Norrdahl and Korpimäki, 1995;Ceballos, 2002;. The dangerous niche hypothesis proposes that the primary benefit of aversive responses to novelty (neophobia) is protection against danger. Therefore, neophobia should be adaptive when novel stimuli are likely dangerous (e.g., if toxic foods, traps or a high level of predation risk characterize an individual's environment; Greenberg, 2003). Novelty avoidance could then be a measure of general risk-aversion, particularly relevant in species that are heavily preyed upon, and are common (target and nontarget) victims of pest management (e.g., Brakes and Smith, 2005;Herborn et al., 2010;Crane et al., 2020).
Here, we aimed to explore the potential role of responses to novelty in coping with urban environments for small mammals. We explored the behavioral responses to novelty of wild-caught, captive-held, rural and urban small mammals, in a battery of replicated tests designed to assess responses to single novelty types. Specifically, we investigated: (i) the temporal consistency and the structure of the neophobic responses to three novelty types: object, food, space, and (ii) differences in the responses to novelty between individuals of rural and urban origin. As study species, we selected common voles (Microtus arvalis, Arvicolidae) and striped field mice (Apodemus agrarius, Muridae). These are non-commensal rodents commonly found both in rural and urban areas of northern and central Europe (e.g., Riegert et al., 2009;Łopucki et al., 2019). Common voles are fossorial, mostly herbivorous, and inhabit open grasslands (e.g., Jacob et al., 2014). Striped field mice are omnivorous, and inhabit tall grasses and bushy underwood (e.g., Babińska-Werka, 1981). Both species are known to display between-individual variation in behavioral traits under both laboratory and natural conditions (e.g., Lantová et al., 2011;Pieniążek et al., 2017;. Both species also show behavioral and cognitive adjustments to human-altered environments (e.g., Dammhahn et al., 2020;Mazza et al., 2020;Mazza and Guenther, 2021). We asked (i) whether there is one latent variable linking all the responses in the different tests; (ii) whether there is covariance between responses to the three novelty types; (iii) whether there is covariance between responses that express 'speed' (that may be read as latency to approach or neophobia) and 'intensity' (extent of the exploration or neophilia), independently of the novelty type; and (iv) whether there are common factors that underlie the specific covariance structures ('general neophobia'). We tested these alternative scenarios of phenotypic integration at the among-individual level and, thus, asked whether specific responses of individuals to novelty are correlated across contexts. Such correlations are indicative of potential genetic or permanent environmental effects on novelty responses, in contrast to correlations at the within-individual level. Specific predictions on how behavioral responses to novelty cluster or vary across novelty types are reported in Table 1. Additionally, we asked whether the behavioral responses to novelty are repeatable, and if they differ between rural and urban individuals. Based on the notions reported above we predicted that animals would show consistent between-individual differences in responses to novelty, and that urban individuals would show lower levels of neophobia and higher levels of neophilia than their rural conspecifics.

Trapping Sites
In total 48 common voles and 42 striped field mice were included in this study, and we report data for 47 common voles and 41 striped field mice. Animals originated from rural and urban wild populations of northeastern Germany. Rural common voles (N = 24; 7 females, 17 males) and striped field mice (N = 27; 13 females, 14 males) were trapped in three different sites in the region of Uckermark in 13 • 71 E, area 3,058 km 2 , Supplementary Figure 1). The area is characterized by agricultural landscapes and comparatively low human population density (for further details see Schirmer et al., 2019). Rural trapping sites were fallow lands between arable fields, composed of grassy areas streaked with nettles, bushes, and trees (ca. 0.85-1.66 ha). Urban common voles (N = 24; 11 females, 13 males) originated from four different sites, urban striped field mice (N = 15; 8 females, 7 males) from three different sites in the city of Berlin (52 • 31 N, 13 • 24 E, area 891 km 2 , Supplementary Figure 1). Urban sites were characterized by heterogeneous vegetation including grassy areas, bushes, and small trees (ca. 0.75-1.16 ha). Rural and urban trapping sites were on average 102.6 ± 3.7 km apart. Rural sites were on average 9.2 ± 6.7 km apart from each other, while urban sites were on average 13.3 ± 5.9 km apart. Common voles have an average home range size of 30-180 m 2 , an average core area of 4-25 m 2 , a daily range length of 9-49 m, and are estimated to disperse between 76 and 110 m per generation depending on sex (e.g., Jacob and Hempel, 2003;Briner et al., 2005;Gauffre et al., 2008Gauffre et al., , 2009Roos et al., 2019). Striped field mice have an average home range size of 2,737 ± 2,046 m 2 , an average core area of 600 ± 446 m 2 , and an average daily distance traveled of 697 ± 312 m (Schirmer et al., 2019. In previous capture-mark-recapture studies conducted in the same sites, it never happened that individuals trapped in one site were found in another . We therefore considered the sampled populations to be independent replicates. All sites were part of the CityScapeLabs or AgroScapeLabs experimental platforms (von der Lippe et al., 2020; https://www.bbib.org). We estimated a sealing index, i.e., the proportion of coverage of natural soil with artificial impervious surface (e.g., buildings and paved roads), as a proxy of the degree of urbanization in the animal's original environment Mazza et al., 2020). The index was calculated in a 5-km buffer around each study site within the CityScapeLabs project for the urban sites, on the basis of the biotope mapping of Berlin (Seress et al., 2014;Buchholz et al., 2018). We used the same method for the rural sites. In brief, this method scores the abundance of paved roads, buildings, and vegetation within each 100 m × 100 m cell in a defined buffer area around the focal study sites, using a GIS (v. 10/ESRI; http://www.esri.com/; Buchholz et al., 2018). Sealed surface closely corresponds to other urbanization indices, such as human population density, disturbance by humans and pets, noise, and light pollution (e.g., Corsini et al., 2017Corsini et al., , 2019Buchholz et al., 2018). Sealing indexes ranged from 2.4 to 3.7 for TABLE 1 | A priori hypotheses on the structure of behavioral responses to novelty in common voles and striped field mice.

Model Hypothesis
Model 0 Null model of response independence.
Model 1 All behavioral responses to novelty are underlined by a single common general dimension (latent variable).

Model 2
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect the speed, the frequency, and the duration of the responses.

Model 3
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect the speed, the frequency, and the duration of the responses, underlain by a common general dimension.

Model 4
Behavioral responses to novelty can be grouped into two dimensions (reflective latent variables), which reflect the speed and the intensity (frequency + duration) of the responses.
Model 5 Behavioral responses to novelty can be grouped into two dimensions (reflective latent variables), which reflect the speed and the intensity (frequency + duration) of the responses, underlain by a common general dimension.

Model 6
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect all responses to novel food and the speed and the intensity component of the responses to other novelty types (such as object and space).

Model 7
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect all responses to novel food and the speed and the intensity component of the responses to other novelty types (such as object and space), underlain by a common general dimension.

Model 8
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect the reactions to objects, foods, and spaces.

Model 9
Behavioral responses to novelty can be grouped into three dimensions (reflective latent variables), which reflect the reactions to objects, foods, and spaces, underlain by a common general dimension.
Frontiers in Ecology and Evolution | www.frontiersin.org rural and 32.0 -54.9 for urban sites, indicating a marked contrast between the level of anthropogenic disturbances the animals were exposed to in their habitats .

Animals and Housing
We captured animals during the same season, between July and September 2018, using multiple-capture live-traps (Ugglan Special Traps n. 1-2, Grahnab AB, Hillerstorp, Sweden). After a pre-baiting period of two nights, traps were activated every night between 18:00 and 19:00, and checked and de-activated every morning between 06:00 and 07:00. Trapping was performed on each site for 7-12 days and was discontinued once the capture success was <2 individuals/grid night. We are confident that prebaiting and capturing over several subsequent days minimized potential bias concerning trap avoidance. Captured animals were sexed, weighed, and checked for reproductive status. Lactating females were immediately released at the site of capture. All other animals were given additional food and transported to our animal holding facility, where they were housed individually in standard polycarbonate cages (Ehret GmbH, Germany, Type III; dimensions: 42 cm × 27 cm × 16 cm) equipped with hay, bedding, and carton shelters. Food (Ssniff NM/V1244-0) and water were provided ad libitum. The animal holding facility receives natural light from windows, and has no thermal insulation, so light, temperature and humidity follow the natural conditions occurring outside the laboratory. Testing took place between October and November 2018. All animals had been in captivity for approximately 3 months when the test battery began (common voles: 91.1 ± 16.2; striped field mice for 109.6 ± 24.1 days). The sample comprised males (28 common voles and 19 striped field mice) and females (19 common voles and 22 striped field mice), all adult and all non-reproductive. After the conclusion of the present experiment, animals were kept in captivity, involved in further behavioral tests (Mazza and Guenther, 2021), and then released at their original trapping sites with a soft-release protocol.

Experimental Procedure
At the start of each testing session, animals were taken from their cage, moved into a separate testing room, and transferred into the experimental apparatus without direct handling. Here, after a 24 h habituation period, animals were presented with a sequence of tests involving a single novelty type: novel object, novel food, and novel space, in randomized order. The habituation phase started no later than 18:00 h and lasted no less than 24 h. A test phase started between 6:00-7:00 and lasted until 18:00-19:00, when the second phase began. Each test lasted 12 h and was immediately followed by the next test of a different novelty type (either food or object or space). Light conditions in the testing room were LD 12:12 (12 h light-12 h dark), and the change was concurrent with the change in novelty type.
Voles alternate several independent short activity bouts (typically 2-3 h) and rest periods, both during daytime and at night, resulting in a polyphasic pattern (e.g., Halle and Lehmann, 1987;Halle, 2000;Wolff, 2007). Apodemus mice seem to be able to adjust their activity rhythms and phases to different factors (e.g., Wróbel and Bogdziewicz, 2015;. Our own observations in the course of a previous field study confirm that striped field mice had activity periods distributed along the 24 h (i.e., they entered the traps during both day and night) in rural as in urban habitats .
Having 12-h-long testing phases allowed us to wait for each animal to be in an active phase, and to approach the novelty in its own time. All animals had activity bouts well within this limit, and all approached the novelty within the first 4.5 h (see below for further details). State differences while testing may account for large variation in behavioral measures, which is why we opted for a protocol that allowed animals to follow their own activity rhythms and provided enough time to habituate to the new setup and housing conditions, as well as enough time to approach novelty on their own terms.
A video-recording system on the ceiling allowed simultaneous monitoring and recording of all tests. Groups of 20 conspecifics were tested simultaneously, i.e., they were in the experimental room at the same time but invisible to each other. From these videos, behavioral parameters (see below for each test) were extracted and analyzed by one observer (IC) with the program BORIS (Friard and Gamba, 2016). The overall testing session lasted 60 h (24 h of habituation and 36 h with access to the three novelty types). This allowed all animals to approach the novel stimuli in their own time, and allowed us to record their responses from the moment they made the first approach (see below). At the conclusion of the third test, animals were weighed and returned to their home cage. Tests were repeated after 15-24 days with a different object, food, space setting.
We removed all observations (N = 17) from the dataset when technical issues in any part of the test battery made the precise quantification of behaviors impossible for the observer (IC) (e.g., condensation on one of the camera lenses). The final data set comprised a total of 158 tests for 47 individual common voles and 41 individual striped field mice (mean number of test rounds per individual ±SD: 1.8 ± 0.4).

Novelty Test Battery Set-Up
The test battery was conducted in a modular set-up, consisting of two plastic boxes (IKEA "Samla" 39 cm × 28 cm × 28 cm) connected by a short tube (8 cm × 5 cm Ø) closed by a nontransparent door (Figure 1). Both boxes contained a thin layer of bedding and were closed by a lid with drilled holes to ensure ventilation. One box was used as a 'familiar space' or home box where the novel object and food were presented; the second was used for the novel space test (see below). The home box was equipped with a dark plastic tube as shelter (8 cm × 5 cm Ø). Food (Ssniff NM/V1244-0) and water were available ad libitum at all times throughout all experimental phases, including the 'novel food' test described below. At the end of each testing session, all equipment was cleaned with either ethanol or soap and water, and left to dry for 48 h.

Novel Object Test
The test started when a novel object was introduced in the home box, on the side opposite to the shelter. Different types of plastic toy animals were used as novel objects: a duck, a horse, a wolf, a shark, and a fish (each approximately 8 cm × 5 cm × 7 cm).
FIGURE 1 | Schematic representation of the experimental set-up used to test one novelty type at a time. After each test session the respective novelty type was removed from the box or, in the case of the novel space, the doors leading to it were closed again.
The type of object was randomized across species and sexes. On the second round of testing, each individual was presented with a different type of object. If the animals were not in the shelter when the novel object was introduced in the box, they were gently guided to it without direct handling. The following parameters were measured over the following 12 h: (i) latency to approach the object, defined as coming within 2 cm of it, with the head turned toward it, after leaving the shelter for the first time; (ii) number of approaches to the novel object (defined as i) during the first hour after the first approach; (iii) overall duration of approaches (defined as i). Recording latencies from the moment the novel object was placed in the cage would have confounded activity levels with responses to novelty, in case an animal had a shorter latency to approach just because it was already active when the object was placed in the cage compared to a conspecific whose activity phase started hours later. To avoid confounding effects of higher activity levels or individual differences in reaction to disturbances (e.g., experimenter interacting with the home box) or slight variations in disturbance strength, with responses to novelty, we recorded and analyzed latency to approach from the moment the animals left their shelter, i.e., when they started an activity phase. Similarly, number and duration of approaches were recorded from the moment of the first approach.
All animals participated in the test, and approached the novel object within 4 h (for further details, see the "Results" section). At the end of the test, the novel object was removed from the home box.

Novel Food Test
The test started when a tray containing an unfamiliar food was introduced in the home box, on the side opposite to the shelter. Animals were presented with two food types in randomized order in the two repeats. We used millet flakes and quinoa seeds for common voles, and standard food pellets flavored with grapefruit or "tropical" fruit syrups for striped field mice. We checked that the novel foods and flavors would be approached and eaten by animals of both species during a short pilot, conducted with animals that were not included in the present study. If the animals were not in the shelter when the novel food was introduced in the box, they were gently guided to it. The following parameters were measured over the following 12 h: (i) latency to touch the novel food, i.e., to make physical contact (touching with nose or paws, nibbling, and sniffing), after leaving the shelter for the first time; (ii) number of contacts to the novel food (defined as i) during the first hour after the first contact, (iii) overall duration of contacts with the novel food (defined as i) during the first hour after the first contact. To avoid confounding effects of higher activity levels or individual differences in reaction to disturbances (e.g., experimenter interacting with the home box) or slight variations in disturbance strength, with responses to novelty, we recorded and analyzed latency to touch the food from the moment the animals left their shelter, i.e., when they started an activity phase. Similarly, number and duration of contacts were recorded from the moment of the first contact. Given the distance between the video cameras (on the room ceiling) and the food trays, we could not distinguish between contact and actual food consumption. All animals participated in the test and approached the novel food within 4.5 h (for further details, see the "Results" section). Familiar food (standard pellet described above in section "Animals and Housing") and water were available ad libitum throughout the experiment. At the end of the test, the tray containing the novel food was removed from the home box.

Novel Space Test
The test started when the door to the second box was opened. To provide a novel space in both rounds of testing, we could not present an empty box, so we attempted to change the conformation of the available space by placing additional structures. The second box was therefore equipped with a round terracotta vase and a plastic closed pipe in the first round, and with vertical walls during the second round, creating a maze-like space (Figure 1). If the animals were not in the shelter when the door to the second box was opened, they were gently guided to it. The following parameters were measured over the following 12 h: (i) latency to visit the novel space for the first time, i.e., to enter the second box with the whole body without the tail, after leaving the shelter for the first time; (ii) number of visits to the novel space (defined as i) during the first hour after the first visit; (iii) overall duration of visits (defined as i). To avoid confounding effects of higher activity levels or individual differences in reaction to disturbances (e.g., experimenter interacting with the home box) or slight variations in disturbance strength, with responses to novelty, we recorded and analyzed latency to visit from the moment the animals left their shelter, i.e., when they started an activity phase. Similarly, number and duration of visits were recorded from the moment of the first visit.
All animals participated in the test and visited the novel space within 4 h (for further details, see the "Results" section). At the end of the test, the door leading to the novel space was closed.

Statistical Analyses
We conducted three sets of statistical analyses (described in detail below). First, we estimated adjusted repeatabilities and patterns of (co)variance among behavioral responses using a Bayesian multivariate mixed-effect model (e.g., Dingemanse et al., 2010;Moiron et al., 2019). As a second step, we fitted structural equation models (SEMs) to compare potential underlying structures of among-individual responses to novelty. We achieved this by testing the relative fit of 10 a priori hypothesized structures of phenotypic integration (Figure 2 and Table 1) (e.g., Dingemanse et al., 2010;Moiron et al., 2019). As a third step, we reduced the number of dependent variables and summarized the behaviors using a principal component analyses (PCA). We then used the summarized scores to compare rural and urban individuals. We scaled all nine variables to standard deviation units (i.e., mean centered and variance standardized) to facilitate comparison of the relative magnitudes of variance components across variables. To meet the normality requirement, latencies and durations were log-transformed prior to analyses. Data from the two species was analyzed separately, following the same protocol. Based on previous studies on small mammals' personality and cognition (e.g., Herde and Eccard, 2013;Mazza et al., 2018), we did not expect differences between the sexes, and we did not include sex in the models to avoid overfitting. The accepted significance level was α ≤ 0.05. All data analyses were conducted with R version 3.6.1 (R Core Team, 2015).

Multivariate Mixed-Effects Modeling
We estimated among-trait (co)variance at the within-and among-individual level using a Bayesian multivariate mixedeffects model [R package MCMCglmm, Markov chain Monte Carlo generalized linear mixed models (Hadfield, 2010)]. To do so, we built a random intercept model where we fitted all nine behavioral variables as response variables, and the trial number (round 1 vs. 2) and the presentation order of the tests as fixed effects. Individual identity was added as a random factor, specified as a random intercept, in each model to control for nonindependent repeated measures of the same individual. We used an uninformative prior and ran 220,000 iterations, with a burn-in of 20,000 and a thinning interval of 100. These parameter settings resulted in low temporal autocorrelation between estimates of subsequent models, which were assessed by graphical diagnostics. The resulting among-individual covariance matrices were used to test our SEM hypotheses. We also calculated adjusted individual repeatability for each behavioral variable as the proportion of the total phenotypic variance not attributable to fixed effects that was explained by individual identity and its 95% credibility intervals based on the posterior distribution (Nakagawa and Schielzeth, 2010;Dingemanse and Dochtermann, 2013).

Structural Equation Modeling
To test the relative fit of alternative hypotheses on the structure of the novelty responses, we applied SEM to 10 a priori conceived scenarios assessing how the different responses were associated among individuals (Figure 2 and Table 1). SEM provides a powerful method for syndrome structure analysis, assessing how well alternative hypotheses fit the available behavioral data. SEMs allow the detection of patterns of phenotypic trait covariance which could be missed using more classical approaches (Dochtermann and Jenkins, 2011;Dingemanse and Dochtermann, 2013). Of the 10 models, model 0 represents a statistical "null" expectation where all responses are independent from each other; model 1 represents a scenario in which all responses were underlain by a single common general dimension (N, general neophobia). Models 2-6 represent different scenarios where behavioral responses to novelty can be grouped into two or three dimensions (reflective latent variables), which reflect the speed, the frequency, and the duration of responses to novelty, and where these dimensions can be underlain by a common general dimension (N). In these models, the dimension intensity combines duration and frequency. Models 7-8 represent hypotheses of behavioral responses grouped into dimensions (reflective latent variables), which reflect the responses to novelty in the form of object, food or potential resource, and space, related to the testing set-ups and where these three dimensions can be underlain by a common general dimension (N). Models were fitted with the package lavaan (Rosseel, 2012). We then compared the models' fit and evaluated the relative support based on several indices, including the Akaike information criterion (AIC), the minimum fit function (FMIN), the normed fit index (NFI), and the parsimonious normed fit index (PNFI).

Principal Component Analysis
Since our sample sizes did not allow us to compare between rural and urban populations via SEMs, we used principal component analyses (PCA) followed by oblimin rotation to reduce the number of dependent variables (Tabachnick and Fidell, 2001). All variables were checked for suitability for a PCA by examination of the determinant of the correlation matrix, the Bartlett test, and the Kaiser-Mayer-Olkin (KMO) criterion (Field, 2012). We retained principal components with Eigenvalues > 1 (Kaiser-Guttman criterion, Kaiser, 1991). The PCA returned three main components that accounted for 62% of the variance in both species ( Table 2). We then assessed the correlation between the three components by Spearman rank correlations. We assessed repeatability of the composite scores using linear mixed effects models (LMM) with individual identity as random factor (Nakagawa and Schielzeth, 2010). We calculated 95% confidence intervals (CI) of repeatabilities for each variable by parametric bootstrapping (N = 1,000 simulation iterations) (Faraway, 2006) using the R package rptR (Stoffel et al., 2017). P-values were FIGURE 2 | Schematic representation of the a priori hypotheses on the structure of behavioral responses to novelty in common voles (Microtus arvalis) and striped field mice (Apodemus agrarius) tested with SEMs for each species separately. calculated based on 1,000 permutations. We compared each component across rural and urban populations using linear mixed effects models, fitted by restricted maximum likelihood, run with the R package lme4 (Bates et al., 2015), with individual identity as random factor and test round as a fixed factor. Time in captivity was never significant and was removed in order to select the minimum parsimonious models. More complex structures with individual identity nested in trapping site did not have a better fit (Supplementary Table 4), so we chose a simpler structure as better suited to the sample size of observations.
equipped with plenty of food and apples as water sources. We covered the traps with grass, branches, and leaves to provide shade and avoid overheating. In case of rainy nights, trapping was suspended. Lactating females were immediately released at the site of capture. All other animals were given additional oat flakes, apples, and cucumbers to provide food and water sources, and transported to our animal holding facility. If two or more animals were captured in the same trap, they were separated to avoid potentially stressful interactions during transport. Traps were firmly secured to avoid shaking and accidental opening. Temperature was kept constant at 20 • C to avoid overheating. Animals were housed individually in the animal holding facility in standard polycarbonate cages (Ehret GmbH, Germany, Type III; dimensions: 42 cm × 27 cm × 16 cm), that incorporated biologically relevant enrichment features such as natural material and refuges. Food (Ssniff NM/V1244-0) and water were provided ad libitum. Light, temperature, and humidity mirrored the natural conditions occurring outside the laboratory. We ensured that the animals learned to recognize the bottles as water sources by placing fresh apples and cucumbers directly above the bottles for the first week of captivity. Animals were never food-deprived or subjected to any form of potentially harmful manipulation. Disturbance by the experimenter was minimized by videorecording the tests, so that the experimenters' presence in the room was not required. After the conclusion of this experiment, animals were kept in captivity and involved in further behavioral tests, and later were released at the capture sites in a waterproof shelter with hay and food to allow them to slowly re-acclimatize to their former habitat.

RESULTS
All animals approached the novel object within 12,000 s, the novel food within 16,200 s and visited the novel space within 14,400 s. Descriptive statistics for measured variables are reported in Supplementary Table 5.
All behavioral variables were repeatable over time (Figure 3 and Supplementary Table 1) and showed considerable variation among individuals (Table 3). For the common voles, it was not possible to identify one model as best explaining the amongindividual covariance structure (Supplementary Table 2). While some indices to assess model fit (FMIN, NFI, and PNFI) had best values for models 8-9, support based on AIC ( AIC < 2) was similarly high also for models 1, 4, and 5. Models 8-9 postulated the integration of the behavioral variables into three dimensions (reflective latent variables), reflecting the reactions to novel objects, foods, and spaces. Model 9 posited that these three dimensions were further underlain by a common general dimension N ( Table 1). Models 1, 4-5 assumed the integration of behavioral variables into two dimensions, reflecting speed and intensity of the reactions, and the presence of a single common general dimension N.
Likewise, for the striped field mice, SEM comparisons identified models 6 and 7 as best explaining the among-individual covariance matrix based on FMIN and NFI, but model 4 was best supported based on PFNI and AIC (Supplementary Table 2). Models 6-7 posited the integration of the behavioral variables into three dimensions (reflective latent variables), reflecting responses to novel food and the speed and the intensity component of the responses to other novelty types (such as object and space). Model 7 postulated that these two dimensions were underlain by a common general dimension N. It might have been too complex, though, and warning messages indicated that "a solution could not be found." Model 4 posited the integration of the behavioral variables into two dimensions (reflective latent variables), reflecting the speed and the intensity (frequency and duration) of the responses to novelty. In both species, the null model was clearly unsupported, and so were the models grouping behavioral responses into separate clusters of speed, frequency, and duration (2-3, Supplementary Table 1).
The PCA returned three main components that accounted for 62% of the variance in both species ( Table 2). For both species, the first component best explained the variance associated with the intensity of the response, i.e., the number and duration of the interaction with the novelty type (object and space) and the second component was associated with the latencies to interact with the novel stimulus. The third component was associated with the intensity (number and duration) of interaction with the novel food in the common voles, and with the latency to approach the food and the duration of interactions with the novel object for the striped field mice. We will refer to these components as PC1(intensity), PC2(speed), and PC3(food) because their categorization, especially for the third, could be controversial. We refer to the corresponding aspect of our hypotheses in the SEM in brackets only for purposes of clarity and brevity. Repeatabilities and correlations between the three components are shown in Table 2 Table 3) and tended to inspect the food less (β = −0.44 ± 0.26, t = −1.69, P = 0.09 -Supplementary Table 3) compared to rural conspecifics. The difference in speed of approach (PC2) was only present in the second round of testing (LM second round subset: β = −0.57 ± 0.23, t = −2.49, P = 0.018). No difference was detected in any of the components for rural and urban striped field mice. PC2 (speed of approach) increased with test round in both species, while PC1 (intensity of the responses) increased for common voles and decreased for striped field mice (Supplementary Table 3).

DISCUSSION
How animals react to novelty might be an important predictor of their capability to adjust to human-altered environments. Using a battery of repeated tests involving novel object, food, and space, we found that measured behavioral responses were moderately to highly repeatable over time. Across the three novelty types and at the among-individual level, we found that responses are linked but the most likely structure of these links remains somewhat ambiguous. Although the null model was clearly unsupported, a priori hypothesized structures resembling a more general neophobic response and those of more specific response types  Correlations were estimated using a multivariate Bayesian mixed model approach, controlling for the effects of order of presentations and test round. Values in bold font indicate important correlations, owing to 95% credibility intervals not overlapping zero.
were similarly well supported. Therefore, deciphering the overall structure of responses to novelty requires further testing. Notably, urban and rural individuals did not differ in their responses concerning intensity, speed, and reaction to novel food, with the exception of common voles being faster in approaching novelty in the second round of testing compared to rural conspecifics.

Structure of the Responses to Novelty
To analyze the structure of responses to novelty, we focused on three ecologically important novelty contexts and on the among-individual level, which should reflect individual intrinsic characteristics and permanent environmental effects (e.g., Dochtermann et al., 2015). For these novelty types and at this level, we could not single out a specific a priori hypothesis as being better supported than all others, although the null model was clearly unsupported in both species, indicating the presence of a structure in the measured behavioral responses. For common voles the best supported model with most, but not all, indices indicated a structuring of the behavioral responses to novelty into three dimensions, reflecting the reactions to objects, foods, and spaces. This would be in line with the idea that the three tests illuminate different aspects of the responses to novelty in different contexts, supporting the idea that different selective pressures act on these axes separately (e.g., Greenberg and Mettke-Hofmann, 2001;Mettke-Hofmann et al., 2002;Mettke-Hofmann, 2014). Similarly, there was not one single scenario that could be selected as unequivocally better at describing the data for striped field mice. Here the a priori hypotheses that appeared better supported pointed toward the integration of the behavioral variables into three dimensions, reflecting responses to novel food, and the speed and the intensity component of the responses to objects and spaces. Specifically, the models identified a pattern where the latencies to approach the novel object and space co-varied, the frequencies and durations of approaches to these novelty types -which we summarized as describing the intensity of the interaction -co-varied on a separate axis, and responses to food grouped together in a third. Notably, the model hypothesizing a latent variable reflecting a general response, that could not be excluded for common voles, was not supported. However, in both species, some of the models did not converge, and a proper comparison among all models could not be made. Since the SEM analysis is not inferring 'causality' but rather helps to reject some a priori hypothesized structures and, thus, to provisionally accept a (set of) given models it always entails some uncertainty and we offer it as a 'work in progress result.' SEMs proved to be particularly sensitive to sample size, and although our sample sizes of observations followed the suggestion of Bentler and Chou (1987) with ca. five observations per estimated parameter, it might have limited our power to resolve structure with SEMs (e.g., Bókony et al., 2012;Wolf et al., 2013). Therefore, we supplemented our analyses with principal component analyses of behavioral variables at the phenotypic level. For both species, these analyses revealed a structure combining the most likely scenarios based on the SEMs, identifying (broadly speaking) a slow-fast continuum to approach the novelty (PC2 speed), one that described the intensity of the responses to object and space (PC1 intensity), and a third on which different aspects of the response to food were mainly represented (PC3 food), with a few differences between species.
Taken together, our results suggest that there are two common underlying factors that drive and structure different aspects of the behavioral responses to novelty. First, the speed of interaction, i.e., how long individuals will take until approaching novelty, and second its extent or intensity. This is in line with previous studies reporting an approach-avoidance conflict stemming from neophilic and neophobic reactions both being elicited by the presentation of novelty (e.g., Murphy, 1978;Mettke-Hofmann et al., 2002;Stieb et al., 2005;Stöwe et al., 2006). Here, however, we cannot tell if more intense interactions are due to higher neophilia or to the amount of information individuals need to collect on any aspect of their surroundings, therefore reflecting individual strategies of fast and shallow versus slow and thorough exploration (e.g., Carere and Locurto, 2011;Reader, 2015;Rojas-Ferrer et al., 2020). In both species the responses to food stood apart from the others. This could be interpreted in terms of behavioral reactions to novel food constituting a separate dimension and requiring a separate set of responses, that are not captured in the avoidance/attraction paradigm (e.g., Marples and Kelly, 1999).
The very weak correlations between behavioral responses at the among-individual level also seem to indicate a high degree of independence between these responses, or at least that they are reflecting separate aspects of the responses to novelty. From these results, we derive a tentative conclusion of responses to novelty being underlain by common factors expressing both the tendency to avoid novelty (neophobia) and drive to explore it (neophilia), but context-specific responses were also present. Specific responses to objects, foods, and spaces might not necessarily offer reliable predictions on the responses of small mammals in the other contexts, and while responses to objects and spaces seem to be underlain by the same factors, novel food appeared to be perceived and reacted to in a markedly different manner, perhaps indicating responses to potential unknown dangers in the first case and to potential resources in the other (e.g., Šlipogor et al., 2016).

Responses to Novelty in Rural and Urban Individuals
Urban areas differ from rural ones in spatial (e.g., habitat configurations, plant composition, and frequency of alterations), object (e.g., rubbish in the park, nesting material, or opportunities), and food (e.g., presence of human-waste, or different seed and arthropod communities) features. The types of novelty considered in this work (space, object, and food) thus encompass some of the main aspects which animals can be faced with in their environment. We did not find reduced neophobia in urban common voles and striped field mice, as most composite scores were similar across habitats. This contrasts with some previous studies reporting lower levels of neophobia in urban animals, and the general predictions made on behavioral adjustments to urban environments (e.g., Sol et al., 2011Sol et al., , 2013Lowry et al., 2013;Ducatez et al., 2017;Jarjour et al., 2019). While further work is necessary to illuminate the mechanisms maintaining comparable levels of neophobia and neophilia in rural and urban small mammal populations, these results suggest that living in the city has not altered the pressures shaping these responses in the two species. Particularly, since small mammals' fitness is strongly determined by the survival component, maintaining comparatively high life-preserving fear of novel stimuli might be the most appropriate response. Previous studies also reported comparable responses to novelty in rural and urban dwellers (e.g., Echeverría and Vassallo, 2008;Bókony et al., 2012). Although colonization and establishment in novel environments are commonly considered to promote flexible adjustment of behavioral responses to both information gathering and risk, it is still unclear whether plasticity in responses to novelty is due to a general or specific modification of fear (Greggor et al., 2016b). Neophobic responses likely evolved to protect the individual from potential harm in unknown situations (e.g., Greenberg, 2003). We have previously shown that rural and urban individuals of both species differ in consistent behavioral traits related to risk taking, with urban individuals being both bolder (i.e., more risk-prone) and more plastic in the way they express their behavioral responses Mazza et al., 2020). In this study, we report evidence of responses to novelty also being stable individual traits, as well as little to no support for differences between rural and urban populations. Taken together, these findings support the view that what is usually referred to as "boldness" indicates consistent individual differences in risk-taking behavior (Réale et al., 2007) in response to a familiar, known risk, such as being in an open unprotected empty space for rodent prey species, while neophobic responses may specifically target unknown potential dangers and be strongly maintained across populations, irrespective of habitat characteristics. The dangerous niche hypothesis posits that heightened aversion to novelty is adaptive when novel stimuli are likely dangerous, e.g., if toxic foods, traps or a high level of predation risk characterize an individual's environment (Greenberg and Mettke-Hofmann, 2001;Greenberg, 2003;Bókony et al., 2012). Differential selective pressures in different stages of colonization might explain the contrasting findings between studies. Depending on the species-specific conditions, it could prove more advantageous to express higher cautiousness during the first stages of colonization, followed by reduced neophobia and increased neophilia, allowing for higher chances of encountering new resources (e.g., Greenberg, 2003;Bókony et al., 2012). For others, especially those that are targeted by deterrents and culling efforts by humans, there might be reasons to increase aversion to novel stimuli even further (e.g., Cowan, 1977;Frid and Dill, 2002;Barrett et al., 2019;Goumas et al., 2020). This is what Greggor et al. (2016b) suggest regarding corvids, whose neophobic responses were as strong as those of their rural counterparts and much stronger than non-corvid urban birds. Similar findings were reported by Bókony et al. (2012) for house sparrows (Passer domesticus). Common voles and striped field mice are non-commensal and not the target of any management program in urban areas. Small mammals are, however, among the most common non-target victims of rodenticides intended for other species (e.g., Brakes and Smith, 2005). Urban green areas also see high densities of predators (e.g., Fischer et al., 2012).
We have no exact information on how long the two species have lived in the city of Berlin, but until ca. 250 years ago Berlin consisted mainly of villages separated by wide rural areas and therefore the process of its urbanization is relatively recent (e.g., Reulecke and Reuleke, 1977;Arandelovic and Bogunovich, 2014; http://atlasofurbanexpansion.org/cities/view/Berlin). Urban animals in our study therefore might not have been urbanized long enough to really develop major differences to their rural counterparts in this respect. Also, green spaces across the city might represent functional corridors for small mammals, eroding any local behavioral adaptations in responses to novelty. Long-term and population genetics data alongside behavioral assessments are needed to assess if and how responses to novelty vary with colonization stage, and in which direction. Since small mammals appear to avoid the road surface itself rather than traffic noise or other emissions, attention should also be given to the specific urban design, where the configuration of green areas that might be used as corridors could impact intergenerational dispersal (e.g., McGregor et al., 2008;Ascensão et al., 2016). To further illuminate the significance of responses to novelty, future work should focus on whether the novelty elicits a physiological stress response, as well as comparing responses to novelty with measures of risk responsiveness toward different potential stressors, and to the presence of humans specifically (e.g., Herborn et al., 2010;Barrett et al., 2019;Goumas et al., 2020).
Given the lack of information on individual experiences, we regard our results as rather conservative because this heterogeneity in experiences might actually erode some patterns visible only under controlled developmental conditions. Longer term studies that allow the tracking of the history of individual experiences and other developmental aspects are a necessary next step. These studies are only possible with captive-born individuals under laboratory conditions, which also have their drawbacks, for example due to reduced stimuli experiences compared to wild conditions. The present study -using wildcaught individuals of unknown developmental and experience background -is a first step toward understanding the structure of novelty responses, particularly since we combine individuals from different populations.
Since animals of both origins were kept in captivity for ca. 3 months before testing, it might be argued that their responses to novelty, as measured, are a product of habituation to the laboratory environment. Interpreting lab-obtained data to investigate evolutionary patterns can be problematic; at the same time, a shorter delay between capture and test might not have prevented a non-natural response, because the living conditions were different, and the stress of adjusting to captivity could have been a bigger confounding factor (e.g., Niemelä and Dingemanse, 2014). Several studies on various long-lived taxa report consistency in responses to novelty in captivity and in the wild (e.g., Herborn et al., 2010;Minderman et al., 2010;Pellegrini et al., 2010;van Overveld and Matthysen, 2010), throughout the lifespan (e.g., Payne et al., 2021) or across contexts (e.g., Stöwe et al., 2006;Tebbich et al., 2012;Castanheira et al., 2013). The few studies that, to our knowledge, address this issue in small mammals also report consistent responses to novelty across time and contexts in general, and across the laboratory-field context in particular, showing that neophobia measured in captivity reflected differences in neophobia in the wild (e.g., Yuen et al., 2016;Guenther and Brust, 2017;Schuster et al., 2017). Similarly, while a previous study of our group shows that urban animals were more likely to adjust their behavioral responses to the testing environment , we also showed that cognitive aspects such as innovation propensity still differed between urban and rural populations after a much longer time in captivity (approx. 1 year) (Mazza and Guenther, 2021). Further studies measuring novelty responses directly in wild populations, challenging as they are, might shed light on this aspect.
The one aspect in which we could detect differences between rural and urban common voles was that urban voles approached novelty faster compared to rural common voles. Considering the importance of survival for prey animals, we think this result could indeed be explained by the amount of different stimuli animals are exposed to in different habitats, and not to a general lower aversion to novel stimuli. Notably, these differences did not emerge until the second round of testing. We repeated the tests to assess repeatability of the response and we used variation of the stimuli to maintain the aspect of novelty. Results indicate that we could not completely avoid habituation to the test, though, as the speed of the response increased for all animals in the second round, and the intensity varied for both species: voles increased the intensity of the interaction with the novelty while striped field mice decreased it. Responses to food did not show any habituation effects, further suggesting that they form a separate axis of responses to novelty. Urban common voles habituated faster than rural conspecifics, approaching the novelty at higher speed in the second round of testing (e.g., Vincze et al., 2016). Faster habituation rates are often interpreted as a sign of higher behavioral plasticity and/or a simple form of learning (e.g., Martin and Réale, 2008;Snell-Rood, 2013). This finding is in line with a wealth of literature reporting higher levels of behavioral plasticity in urban dwellers, both at the inter-and intra-specific levels (e.g., Sol et al., 2002Sol et al., , 2013Lowry et al., 2013). In a previous study involving these individual common voles, we also found higher levels of behavioral plasticity (reversible phenotypic plasticity) in boldness and activity of urban individuals . No such difference was detected for striped field mice, although here too we could have expected higher levels of behavioral plasticity, as these individuals were also tested in a battery of problem-solving tests, often considered as an indicator of behavioral plasticity (e.g., Reader and Laland, 2003;Stamps, 2016;Audet and Lefebvre, 2017), and urban mice outperformed their rural conspecifics (Mazza and Guenther, 2021). In that case, rural and urban striped field mice alike required an extensive habituation phase (up to 5 weeks) before they could be presented with the problem-solving set-ups. In that study, our aim was to test problem solving performance unhindered by amongindividual variation in neophobia (Mazza and Guenther, 2021). The long time it took to reach the habituation criterion suggests that two repeats of the test 15-24 days apart might not be enough to detect potential differences in habituation rates between rural and urban striped field mice.

CONCLUSION
We offer this study as a first step toward understanding the structure of the responses to novelty in common voles and striped field mice. Results support the idea that responses are structured by common elements that could be broadly described as neophobia and neophilia toward objects and spaces, while reactions to novel food seem to belong to a separate dimension. Additionally, urban small mammals seem to maintain levels of neophobia and neophilia comparable to those of their rural counterparts, which could indicate that urban environments represent a dangerous niche for small mammals to live in, which rewards comparable behavioral responses to novelty. Further studies are needed to investigate the mechanisms maintaining similar levels of neophobia and neophilia in rural and urban small mammal populations, and the degree to which the constraints faced by very diverse taxa in each habitat type are comparable.

DATA AVAILABILITY STATEMENT
The data that supports the findings of this study are available as Supplementary Material of this article.

AUTHOR CONTRIBUTIONS
VM and MD conceived and designed the study. IC performed the data collection. VM analyzed the data, with help from MD and JE. VM wrote the manuscript. All authors contributed to the final draft and gave final approval for publication.

ACKNOWLEDGMENTS
This work was made possible by funding from the German Federal Ministry of Education and Research (BMBF) within the Collaborative Project 'Bridging in Biodiversity Science (BIBS)' (funding number 01LC1501A), supporting prior steps of the project. During manuscript preparation, MD was supported by a grant of the German Science Foundation (DA 1377/4-1). We acknowledge the support of the Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Potsdam. We thank Angelika Beck for preparation of the city project and Elisa Lösche and Stephanie Thiele for their help on the field during trapping. Thanks to Angela Puschmann and Elisa Lösche for helping to care for the animals in the laboratory. We thank the reviewers and handling editor for their constructive comments and suggestions.