Exploratory behaviour divergence between surface populations, cave colonists and a cave population in the water louse, Asellus aquaticus

Behaviour is considered among the most important factors in colonising new habitats. While population divergence in behaviour is well-documented, intraspecific variation in exploratory behaviour in species with populations successfully colonising and adapting to extreme (compared to the ‘typical’) habitats is less understood. Here, by studying surface- vs. cave-adapted populations of water louse (Asellus aquaticus), we tested whether (i) adaptation to the special, ecologically isolated cave habitat includes a decrease in explorativeness and (ii) recent, surface-type cave colonists are more explorative than their surface conspecifics from the source population. We repeatedly tested dispersal related novel area exploration and dispersal speed in both the presence and absence of light. We found that surface populations showed higher behavioural activity in dark than in light, and they were more explorative and dispersed faster than their cave conspecifics. Recent colonists showed a trend of higher dispersal speed compared to their source surface population. We suggest that extreme and isolated habitats like caves might work as ‘dispersal traps’ following successful colonisation, because adaptation to these habitats includes the reduction of explorativeness. Furthermore, we suggest that individuals with higher explorativeness are likely to colonise markedly new environments. Finally, we provide experimental evidence about surface A. aquaticus moving more in dark than in light. Environmental conditions in caves are differing drastically from those of the surface. Consequently, animals colonising subterranean habitats are subject to different selective forces than those experienced by the ancestral surface-living population. Behaviour is believed to be a key factor in successful colonisation to novel habitats; however, intraspecific behavioural variation in species with both surface- and cave-adapted populations is less known. Here, we compared dispersal related novel area exploration and dispersal speed across surface and cave-adapted populations of the freshwater crustacean Asellus aquaticus. Our results show that cave-adapted A. aquaticus are significantly less explorative and disperse slower than surface-type populations, indicating that caves may act as ‘dispersal traps’, where adaptation includes the loss of explorativeness. Also, recent cave colonists show a trend to be faster dispersers than peers from the surface source population, suggesting that individuals with higher explorativeness are likely to colonise markedly different new environments.


Introduction
Heritable variation in behavioural traits within species and populations is an important factor for successful colonisation of novel habitats or withstanding environmental variation (Morse 1980). Consistent between-individual variation in behaviour over time and across ecological situations within population (i.e. animal personality) is a common and widespread phenomenon across the animal kingdom (for reviews and meta-analyses see e.g. Sih et al. 2004;Réale et al. 2007;Bell et al. 2009;Dingemanse and Wolf 2010; Communicated by N. A Dochtermann Garamszegi et al. 2013). Furthermore, a growing number of theoretical and empirical works suggest that individuals settle in habitats, which suits their personality (here: individual behavioural configuration) the best, implying the existence of personality matching habitat choice (Edelaar et al. 2008;Jacob et al. 2015;Saltz et al. 2018;Saltz 2019). Between-individual variability in dispersal is associated with differences in various behavioural traits: dispersers (or colonisers) seem to be more risk-taking, more explorative, asocial or aggressive than non-dispersers (Blumstein et al. 2009;Cote et al. 2010a, b;Spiegel et al. 2017). Betweenindividual behavioural variation may be especially important in the colonisation of markedly different new habitats. By 'markedly different habitat', we mean a new habitat that has fundamentally distinct environmental properties compared to the habitats occupied by the given species in general. In other words, such habitat differs in one or more major environmental parameters (e.g. a marked shift in mean temperature or composition of food sources, presence of significant new predator/competitor species) from that of the colonising population's habitat. Such marked habitat shifts are not frequent for the majority of species, but a good example can be the colonisation of caves by surface species.
Caves and related subterranean habitats are characterised by the absence of light, food scarcity and simplified communities, and are strongly buffered against daily, seasonal and yearly environmental variations (Romero 2009;Borowsky and Cohen 2013;Culver and Pipan 2019). Adaptation to the cave environment, often summarised as troglomorphism (e.g. eye degeneration, depigmentation, lowered metabolism, extended life cycles, absence of circadian rhythm; see Porter and Crandall 2003;Mejía-Ortíz et al. 2006), makes it almost impossible for a troglomorphic population to leave the cave and successfully (re)colonise surface habitats. In accordance with this notion, recent empirical studies indeed suggest that dispersal ability of troglobionts (obligate cave dwellers) is, if not zero, extremely limited (see Lefébure et al. 2007;Trontelj et al. 2009;Stern et al. 2017;Balogh et al. 2020). Hence, markedly different and isolated habitats like caves might work as a 'dispersal trap' after colonisation by surface populations. Two main hypotheses are explaining the transition from surface to subterranean habitats. First, the 'adaptive shift' hypothesis suggests that populations of a surface (epigean) species invade subterranean (hypogean) habitats to exploit novel resources (Howarth 1980). Second, the 'climatic relict' hypothesis suggests that a species may be forced to colonise underground habitats to avoid uninhabitable environmental conditions on the surface, for instance, cave colonisation can be triggered by glaciation events (Peck and Finston 1993;Rivera et al. 2002;Danielopol and Rouch 2005;Juan and Emerson 2010). However, not only climate change can induce cave colonisations, but also other environmental challenges, like predation pressure, can cause the same effect (Romero 1985;Tobler 2009).
At any rate, the evolution of dispersal-related behaviours following the successful colonisation of markedly different and isolated habitats acting as dispersal traps is rarely studied. Generally, cave-adapted organisms are expected to have lowered metabolic rates and decreased movement activity (see Hüppop 2000;Hervant et al. 2001), an extreme example being the olm (Proteus anguinus), which shows a surprisingly low movement activity, revealed by a recent capture-mark-recapture study . Theory also suggest that environments with high temporal stability and predictability should favour reduced exploration (Sih et al. 2004;Careau et al. 2009). In line with this, Mettke-Hofmann et al. (2002) found that parrot species that live in complex habitats, such as forest edges, are more explorative than species found in simpler habitats, implying that habitat complexity favours high explorativeness. Hence, dispersalrelated exploratory behaviour is expected to show a clear decrease in cave-adapted organisms; however, empirical data are scarce at best, and results seem to be inconclusive in this regard. For example, the only laboratory experiment to our knowledge that directly studied dispersal-related novel area exploration in surface and cave-adapted populations of the freshwater crustacean Asellus aquaticus found no difference between the morphotypes (Brengdahl 2016). Furthermore, we do not know whether personality variation within the source population plays a role in the colonisation, i.e. whether the boldest, most explorative individuals enter the markedly different new habitat, or whether colonisation happens randomly due to environmental pressures explained by either or both of the adaptive shift and climatic relict hypotheses.
Here, we aimed to study the link between exploratory behaviour and cave colonisation on different levels. First, we were interested in how exploratory behaviour varies between surface and cave populations. Since we assumed that caves act as dispersal traps and the environment within caves has low variability, we predicted that explorativeness shows a significant decrease in cave populations compared to surface ones. Second, we tested whether individuals actually dispersing into caves are different from the source population's mean in their explorativeness. We predicted that the most explorative individuals from a population are the ones entering the markedly different new habitat, i.e. potential colonists express higher explorativeness than the average in their source population. Our model was the water louse, Asellus aquaticus (Linnaeus 1758). The species is a widespread freshwater crustacean that colonised some caves in Europe, where its populations exhibit 'troglomorphic' phenotypes (lack of pigmentation, loss of vision, enlarged sensory, ambulatory appendages, etc.) (Prevorčnik et al. 2004;Konec et al. 2015).

Study system and sampling
We studied three surface populations and one cave population. Molnár János Cave is a hypogene (water in the cave is not coming from the surface) cave filled with thermal water of constant temperature of 23-24 °C (Erőss et al. 2006;Bodor et al. 2015). There is no exogenous food in the cave, only endogenous bacteria-forming mats (see Herczeg et al. 2020 for details). The cave community is extremely simple, and there are no predators of A. aquaticus present. The cave is connected to a surface pond formed by the water outflow right at the cave entrance (47.518277° N, 19.035999° E), called Malom Lake, and there is no physical barrier between the habitats. Despite this fact, the troglomorphic population in Molnár János Cave is genetically isolated from the surrounding surface populations (including Malom Lake) for at least 60 000 years (Pérez-Moreno et al. 2017). The water temperature in Malom Lake is identical to that of the cave and constant all year round. However, Malom Lake is subjected to the natural surface light regime and can be seen otherwise as a typical surface habitat, including the presence of fish predators, i.e. guppies (Poecilia reticulata) that were introduced here during the twentieth century (Berczik 1956) and occur in extremely high density as native fish are absent.
We sampled Molnár János Cave, Malom Lake, recent cave colonists (surface type individuals coming from Malom Lake, collected in Molnár János Cave ca. 250 m from the cave entrance), and two surface populations from the vicinity, Gőtés Lake (47.354357° N, 19.22980° E) and the Dunakeszi Peat-moor (47.615613° N, 19.126392° E) between 16 and 18 May 2019. Samples were collected by hand sorting with a mesh net, except for the animals from the Molnár János Cave, where a modified Sket bottle was used (Chevaldonné et al. 2008) and cave diving was necessary. We considered collected individuals larger than 4 mm as adults (Hasu et al. 2007;Bloor 2011). As gravid females of A. aquaticus form a brood pouch (marsupium) to carry offspring (see Lafuente et al. 2021) and likely display different behavioural activity, we only used non-gravid females to avoid any bias in this respect. Gőtés Lake and Dunakeszi Peat-moor populations experience natural surface light regime and temperature fluctuations of normal, non-thermal freshwaters typical to the region. These populations represent typical surface habitats widespread in the region, and we chose them randomly. Surface type individuals recently entering the cave (hereafter 'potential colonists') from Malom Lake most likely do not form a population, but for simplicity's sake, we will refer to all compared groups (four populations and the colonists) as 'populations'.
All animals were transported immediately after capture to the aquacultural facilities of the Eötvös Loránd University (Budapest, Hungary) and housed individually in 90 × 25 mm (diameter and height, respectively) plastic Petri dishes, with sandpaper-coarsened bottoms aiding the animals' normal movement (Fišer et al. 2019). Cave and colonist individuals were kept in constant darkness, while surface populations in a 16-h light:8-h dark daily light cycle (for the methodology, see the next section). The temperature in the lab was set to 23-24 °C, which is the constant water temperature in the Molnár János Cave and Malom Lake, and within the natural range of the other surface populations at this time of sampling. Water collected from the natural habitats was used for keeping and testing throughout the laboratory period. We had the following number of adult individuals in the tests reported here: Molnár János Cave: 10 males (M)/8 females (F), colonists: 8 M/5 F, Malom Lake: 17 M/13 F, Gőtés Lake: 14 M/13 F; Dunakeszi Peat-moor: 14 M/15 F.

Behavioural assays
We used custom-built chambers for both keeping the study animals under the above-detailed light settings and videorecording their behaviour. The chambers had the following dimensions: 100 cm length × 55 cm width × 105 cm height. The chambers' sides and top were covered by black plastic sheets to block any incoming light. The chambers were equipped with two light sources. On the top, we installed LEDs imitating daylight (colour temperature = 4500 K, colour rendering index > 90), while on the bottom, we installed infrared LEDs (wavelength = 920 nm). This infrared wavelength is outside of the visible range of A. aquaticus (Dember and Richman 2012). The lights could be switched on/ off from the outside. Daylights were used to produce the planned daily light regime, while infrared LEDs were used to make the video-records in either light or dark. On the bottom of the chambers, we placed an opaque plexiglass sheet to diffuse the infrared light evenly and to provide a platform for the Petri dishes. In all chambers, we mounted four webcams (Logitech C920 Full HD; Logitech, Lausanne, Switzerland) to the top that were modified for infrared videorecording. We used the OBS Studio software (OBS Studio Contributors) to capture videos (5 frames per second, HD [1280 × 720 pixels] resolution).
The assays took place between 21 May and 1 June on 12 consecutive days. Animals were allowed to acclimate under their natural light regime without any disturbance. Potential colonists were acclimated in dark, since they were collected from the cave. During acclimation and testing, we provided no food for the test individuals. The tests started approximately at 11.00 am (UTC + 2.00) on each day. Individual behaviour was tested in chambers described above. To test novel area exploration and dispersal speed, we used a cluster of 18 open-top mazes made of Plexiglas (see Fig. 1). The lateral walls were black, while the bottom was transparent and sandpaper-coarsened to enable animals' normal movement (Fišer et al. 2019). A long division wall was installed in the middle of each maze (5.4 cm × 44.25 cm × 3 cm, width × length × height, respectively) that resulted in a U-shaped compartment. In each branch of the maze, we put seven removable obstacles (2 cm wide). This, considering the tip of the division wall too, resulted in 15 obstacles an animal had to pass around to get from one to the other end of the maze. Each population was randomly split into two groups: One of them was tested in the presence of light, while the other group in darkness. After 48 h, treatments were changed; thus, each subgroup was measured in light and dark alternately. Altogether, every individual was tested six times (three times in light + three times in dark), resulting in a total of 702 observations. We have to note here that the light treatment is completely unnatural for the cave population and can only be used to test for light avoidance (based on genetic-studies, cave-dwelling A. aquaticus probably retained the ability to detect light; Pérez-Moreno et al. 2018). However, both treatments are biologically relevant for the surface populations. It is noteworthy that even though the species is intuitively treated as diurnal, there is a field study indicating that A. aquaticus might move more during the night than during the day (Andrikovics 1981). We could test a maximum of three groups daily and testing order was randomised each day. Arenas were filled with water from the habitat of the tested population. We used a pipette to transport individuals from their Petri dish to a transparent cylinder placed at the starting point of each maze. The maze used for each individual was randomly selected beforehand. Once all individuals were placed in the maze, we waited 5 min and started recording (lasted for 60 min) and removed the cylinders to allow free movement.
We extracted two dispersal-related variables from the videos. Total number of obstacles crossed in any directions (an animal was considered to cross an obstacle when its whole body crossed) correlated highly with the farthest obstacle crossed (Spearman's rho = 0.87; 95% confidence interval (CI) = 0.84-0.89), and we chose the former to describe 'novel area exploration' because it involves an element of thoroughness. The second variable was 'dispersal speed' (for individuals that reached the farthest possible obstacle [15], the sequential number of this was divided by the time [s] required to reach it, while in cases when individuals did not reach the farthest possible obstacle, we divided the sequential number of the farthest obstacle crossed by 3600 s, i.e. the entire duration of the experiment). We emphasise that both variables are describing exploration sensu Réale et al. (2007), i.e. both variables are describing movement behaviour in a new (potentially risky) situation; hence, they are clearly different from activity sensu Réale et al. (2007), which must be measured in a non-novel and non-risky situation.  (bottom); note that the animal on the picture was magnified for better visibility. Asellus aquaticus were individually tested for crossing the obstacles and exploring the maze; total number of obstacles crossed in any direction represented 'novel area exploration' (an animal was considered to cross an obstacle when its whole body crossed) and 'dispersal speed' (for individuals that reached the farthest possible obstacle [15], the sequential number of this was divided by the actual time [s] required to reach it, while in cases when individuals did not reach the farthest possible obstacle, we divided the sequential number of the farthest obstacle crossed by 3600 s, i.e. the entire duration of the experiment)

Statistical analyses
To analyse population, treatment and sex effects on novel area exploration, a generalised linear mixed model (GLMM) with negative binomial distribution and logit link function was used. We chose this method, as transformation (log or square root) of count data is generally contraindicated (O'Hara and Kotze 2010). We added 'population' (the four sampled populations and the potential colonists), 'treatment' (light vs. dark), sex and their interactions as fixed effects. Given that the interpretation of three-way factorial interactions can be problematic, and to avoid overparameterisation, we included only the two-way interactions. To control for habituation to the test setup, we added the standardised (mean = 0, sd = 1) order of trials (hereafter, 'time') as a single fixed effect. We added individual identity as a random intercept. We also added random slopes (individual × time) as a random term, but left it in the final model only if it improved model fit. Error distribution and link function applied in the GLMMs were chosen after inspection of Q-Q plots of the model residuals. Fixed effects were tested by Wald's chi-squared tests and random effects by likelihood ratio tests. P values for the likelihood ratio tests were calculated following Zuur et al. (2009). We extracted the model's estimated marginal means using the emmeans package (Lenth 2019). To compare groups, we looked for the presence/absence of overlaps between 85% CIs following Payton et al. (2003), who demonstrated that the lack of overlap in 83-84% CIs is analogous to a P-value < 0.05. We report the proportion of explained variance by the fixed factors (marginal R 2 ) and by both fixed and random factors (conditional R 2 ) available in the MuMIn package (Barton 2009) based on the method of Nakagawa and Schielzeth (2013). Dispersal speed data were square root transformed to achieve normal distribution of the model residuals. We analysed this behavioural variable using a linear mixed model (LMM) built in the same way as described for the GLMM above. We built both the GLMM and the LMM with the R packages lme4 and lmerTest (Bates et al. 2015;Kuznetsova et al. 2016, respectively) in R 4.1.0 (R Developmental Core Team 2021).
The rptR add-on package (Stoffel et al. 2017) was used to calculate repeatability (a statistical test for the presence and 'strength' of animal personality, i.e. consistent among-individual differences over time or across ecological contexts; see Réale et al. 2007;Sih et al. 2012;Niemelä and Dingemanse 2018). Enhanced agreement repeatability (hereafter: eaR) for novel area exploration and dispersal speed in the pooled sample (i.e. all populations combined) was calculated separately for light regimes. This method allows us to fit improved models, in which the variance explained by fixed effects is calculated by the variance in the linear predictor, including the fixed effects' variance in the denominator (see Stoffel et al. 2017). Models were parameterised as described above; nevertheless, as we fitted them separately for different light regimes, treatment effect and its interactions were omitted. We built GLMMs for the negative binomially distributed novel area exploration data, following the methods of Nakagawa and Schielzeth (2010), which utilise a multiplicative overdispersion GLMM with a logit-link and using penalised quasilikelihood (PQL) estimation for repeatability on the original scale. Significance of eaR estimates (i.e. for random terms) was provided by randomisation tests, giving robust measures of statistical significance in the case of non-Gaussian data (Nakagawa and Schielzeth 2010). 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). Quantification of uncertainty for the variance explained by fixed effects (as for other variance components) was provided by parametric bootstrapping. LMMs were run to estimate repeatability for dispersal speed. CIs were calculated by nonparametric bootstrapping, while significance for eaR estimates is provided by likelihood ratio test, both sampled at each 1000th iteration.

Results
Repeatability estimates for novel area exploration and dispersal speed are given in Table S1. For novel area exploration, estimates in the pooled sample indicated significant, moderate (compared to Bell et al. 2009) repeatability in both treatments (light: eaR = 0.29 (95% CI = 0.13-0.35); dark: 0.22 (0.09-0.28)). Repeatability for dispersal speed in the pooled sample was present only in the presence of light, showing moderate strength (eaR = 0.21 (0.11-0.33)).
Population, treatment and their interaction all affected novel area exploration (Table 1). All other populations in both treatments showed higher novel area exploration than cave individuals in dark (Fig. 2). The cave population showed increased novel area exploration comparable to those of the surface populations in the -for them unnatural -light treatment. Surface populations explored more in dark than in light, the trend being weak in Dunakeszi Peat-moor. Treatment had no effect on the potential colonists' exploration. Colonists had higher exploration in light, than surface populations, including their source population (Malom Lake), but their exploration in dark was similar to the surface populations. Furthermore, we found a significant sex × treatment interaction effect (Fig. 3): Males were more active than females in dark, but not in light (see Table S3). Habituation was also significant: individuals became less active by time (Fig S1a). Finally, individual differences in novel area exploration and habituation were both significant. The fixed effects explained 36.9% of the total variance, and the full models explained 74.2%.
We found significant population and population × environment interaction effects on dispersal speed (Table 1). All populations in dark and Dunakeszi Peat-moor and potential colonists in light dispersed faster than the cave population in dark (Fig. 4). Cave A. aquaticus showed comparable dispersal speed to the surface populations in the -for them unnatural -light treatment. Potential colonists from the Malom Lake dispersed faster than their source population (Malom Lake) in light and tended to do so in in dark (see Fig. 4). We also found significant habituation in dispersal speed, with a decrease by time (Fig S1b). Individual differences in dispersal speed were significant. The fixed effects explained 12.2% of the total variance, and the full models explained 27.2%.

Population divergence in behavioural types
Our understanding of the evolution of dispersal-related behaviours in populations adapted to isolated and specialised habitats, i.e. dispersal traps, is rather incomplete. We expected generalist, surface populations of A. aquaticus to be more explorative (i.e. showing higher novel area exploration and being faster dispersers) than conspecifics adapted to the Molnár János Cave. This expectation is supported by the results. All surface populations showed higher exploration, irrespective of light treatment, than cave-adapted individuals in the dark (the natural condition for cave individuals). The dispersal speed patterns were similar, when surface populations were tested in the dark; however, the divergence was weak in Gőtés Lake and Malom Lake individuals tested in light. A. aquaticus in the Molnár János cave are adapted to an isolated, stable, predictable and (compared to the surface habitats) homogenous environment with no predators and high quantity of probably low-quality endogenous food (see Herczeg et al. 2020;. Thus, this special environment most probably acts as a dispersal trap, i.e. the presumably photophobic cave-adapted individuals are unlikely to leave the cave. Increased behavioural activity, comparable to surface populations, could be observed in the Molnár János Cave population in the presence of light. In environments with adverse stimuli, it is expected that novel area exploration is increased as it enhances the chance of finding a more suitable environmental patch, while slowing down can indicate the onset of favourable conditions (Fraenkel and Gunn 1961;Breed and Moore 2021). In our case, increased activity is likely caused by photophobia in the cave-adapted individuals (Janzer and Ludwig 1952;Fišer et al. 2016). Note that despite the marked eye reduction of the cave-adapted individuals, they still preserved the ability to detect light (Pérez-Moreno et al. 2018).
Because of the lack of large predators in most subterranean habitats, behaviour of cave-adapted individuals is virtually not affected by trade-offs between foraging profitability and predation risk. Furthermore, food availability in the majority of caves is limited compared to surface habitats (Culver and Pipan 2019). Therefore, increased activity in caves is expected to find food fulfilling energetic demands (Culver and Poulson 1971;Hüppop 2000). On the other hand, as increased movement activity itself leads to higher energy expenditure, prolonged food shortage was shown to reduce locomotor activity (together with reduced metabolic and respiratory rates) in several hypogean taxa (e.g. Niphargus sp. and Stenasellus virei Hervant et al. 1997;Hervant and Renault 2002;Proteus anguinus Hervant et al. 2001). Although endogenous bacterial mats in the Molnár János Cave can be found in large amounts, our recent findings indicate that cave-adapted A. aquaticus in the Molnár János Cave not just maintained the ability to feed on decaying leaf litter, but actually, this food is preferred over bacterial mats despite the latter are the only available source of food in the cave . Regarding the nutritional content of Molnár János Cave bacterial mats, we possess no exact information, but it was shown recently that potentially toxic metals (e.g. As, Hg, Pb, Sn, Sr, Zn) may accumulate in the biofilm (see Dobosy et al. 2016;Enyedi et al. 2019). This might result in poor quality diet and, along with our previous behavioural observations , indirectly indicates that bacterial mats in this cave might be seen as an obstacle for colonisation (surface populations unconnected to the cave avoided bacteria almost entirely) and that the Molnár János Cave population might be somewhat food deprived after all. In a previous study (Berisha et al. 2022), we found elevated movement activity in the Molnár János cave population compared to surface populations in a familiar environment, while in the present study, we found decreased explorativeness in the cave population in a novel environment. Based on these results, we suggest that cave-adapted individuals moves more in a familiar, perceived risk-free situation, probably as an adaptation to the lack of predators to maximise foraging success, but they are less active in a novel situation. Fig. 4 Dispersal speed in the absence vs. presence of light in the four tested populations of Asellus aquaticus (significant population × environment interaction). White background indicates surface type individuals, while light grey indicates surface type individuals from the Malom Lake found in the cave, and dark grey indicates cave-adapted individuals. DM, Dunakeszi Peat-moor; GL, Gőtés Lake; ML, Malom Lake; C(ML), Colonists from the Malom Lake; MJC, Molnár János Cave. Estimated marginal means and 85% confidence interval (CI) are shown. Note that lack of overlap in 83-84% CIs is analogous to a P value < 0.05 (Payton et al. 2003)

Light-induced behavioural plasticity and sexual dimorphism
In addition to the detected photophobia in cave-adapted A. aquaticus (see above), we were curious about effect of light regime on the behaviour of surface populations too. While these populations are intuitively expected to be diurnal, and in fact, they are active during the day, Andrikovics (1981) showed that surface-dwelling A. aquaticus are three times more likely to be found in funnel traps during the night than during the day. This suggests an increased night-time activity of A. aquaticus, plausibly as an outcome of low predation rate in darkness. Our results are supporting Andrikovics's (1981) observation: surface populations showed increased novel area exploration in dark compared to light, showing similar trends in dispersal speed. Based on these results, we have to reject the notion of the primarily diurnal activity of the species. Increased behavioural activity in darkness could also explain why A. aquaticus is successful in colonising subterranean habitats.
We also found male A. aquaticus to be more explorative than females in darkness. Male A. aquaticus were previously shown to seek actively for mates during the mating season, unlike females (Bertin et al. 2002), and males are performing precopula or mate guarding (Thompson and Manning 1981), which in the significant majority of cases followed by fertilisation (Eroukhmanoff et al. 2009). Therefore, the higher behavioural activity of males observed in our study can be explained by reproductive behaviour. Male-biased sexual dimorphism in exploration, expressed only in dark, further emphasises the importance of night activity for the species. However, this question definitely needs further studies.

Behaviour of potential cave colonists
Personality related dispersal is a well-documented phenomenon (e.g. Blumstein et al. 2009;Cote et al. 2010a;Spiegel et al. 2017). However, less is known about personality related colonisation of markedly different habitats (but see studies about colonising urban environments, e.g. Schuett et al. 2018;Baxter-Gilbert et al. 2019). Cave colonisation by surface populations is explained by two hypotheses: the first assumes that the driving forces are the new, yet unexploited resources provided by the caves, while the second states that caves serve as refuges from the harsh surface environments (Howarth 1980;Peck and Finston 1993;Rivera et al 2002;Danielopol and Rouch 2005;Juan and Emerson 2010). Irrespective of the environmental driver, a subset of surface individuals must enter the new cave environment for the start of the colonisation process. Even if we consider the high behavioural activity of A. aquaticus in total darkness (see above) as a sign of exaptation to the cave life (see Pérez-Moreno et al. 2017), aquatic cave environments are still different (simple and stable biotic and abiotic environment) from surface water bodies. This is especially true for the Molnár János Cave, where the only food sources are endogenous bacterial mats absent from the surface habitats. These bacterial mats are unpreferred compared to surface food or unnatural, but nutritious food even by the cave population (see Herczeg et al. 2020;. Theory predicts that individual dispersal within the same habitat type is linked to personality; as intuitively expected, bolder (i.e. more active, explorative, risk-taking) individuals disperse further than their shier conspecifics (Dingemanse et al. 2003;Cote et al. 2010aCote et al. , b, 2017Myles-Gonzalez et al. 2015). Colonisation of new habitats (in an urbanisation context) was also linked to personality in a similar way (Atwell et al. 2012;Bókony et al. 2012;Schuett et al. 2018;Baxter-Gilbert et al. 2019). Based on this, we predicted that surface A. aquaticus individuals from the Malom Lake that were found in the cave ('potential colonists') are more explorative than the average Malom Lake individual. Note that 'accidental' entrance to the cave is highly unlikely, as the Boltív Spring is flowing from the Molnár János Cave to the Malom Lake, which means that potential colonists have to move actively against the current to enter the cave. We have to also note that we have no information on how much time potential colonists have spent in the cave; it is equally plausible that they entered recently (i.e. within days) or that their parents entered the cave and they were born there.
Regarding novel area exploration, potential colonists showed higher explorativeness under both light treatments, than the other surface populations (including Malom Lake) in light, while they did not differ in dark. The difference was similar for dispersal speed in light, but colonists also tended to disperse faster in dark. Hence, it seems plausible to suggest that more explorative surface individuals are more likely to colonise a markedly different habitat. However, we have to consider alternative explanations too. We found recently that potential colonists from the Malom Lake highly prefer surface food (decaying leaves) over cave food (bacterial mats); in fact, they spent the most time with feeding on surface food compared to other populations , implying that they were food-deprived ; hence, it is likely that they are starving in the cave. Therefore, it is possible that the detected trend for high behavioural activity is simply a result of starving individuals searching for the preferred surface food.
Unlike in other surface-type populations, behavioural activity of potential colonists from the Malom Lake did not decrease in the presence of light. As mentioned above, we possess no information regarding the time potential colonists spent in the cave; further, our experimental design did not allow us to test for genetic vs. environmental effects. Yet, it is very unlikely that this behavioural difference between colonists and other surface populations is the result of evolutionary (genetic) adaptation. Phenotypic plasticity (i.e. a genotype's ability to develop or express alternative phenotypes in different environments; West-Eberhard, 2003) is a more plausible explanation. This behavioural pattern is likely not reflecting genetic adaptation to the cave environment, which should trigger a behavioural response similar to that found in the cave-adapted population.

Conclusions
Taken together, we found support for various links between exploratory behaviour and the colonisation of markedly new habitats. First, we found that following successful colonisation of a highly specialised and isolated cave habitat by A. aquaticus, mean exploratory behaviour of the cave population decreased. The result suggests that highly specialised and island-like habitats, like caves, can act as dispersal traps, where the benefits of explorativeness are negligible. Second, our results suggest that more explorative individuals are likely to enter the cave from the source surface population. This result implies that exploratory personality not only is linked to dispersal, but also determines the probability of colonising new habitats. Finally, we found an additional pattern that provides important details to the biology of A. aquaticus: surface A. aquaticus are not clearly diurnal. They are not only active in light, but also active in total darkness, and actually more active in dark than at light. This finding strengthens the field observation about high night activity of A. aquaticus (Andrikovics 1981) and also suggests that the species is exapted to cave life, explaining the species' repeated successful colonisations of various caves in Europe (Prevorčnik et al. 2004;Konec et al. 2015). We have to note here that our results originate from a single, unique cave system; furthermore, we practically compared four populations (plus potential cave colonists from a surface population) and not two habitats. Thus, we cannot make generalisations about A. aquaticus behaviour in caves and it is reasonable to expect that different selection regimes in other cave systems might cause behaviour to evolve in different directions; further, we cannot exclude the possibility of stochastic events in trait evolution via founder effects or genetic drift (Wessel et al. 2013;Miller et al. 2020). However, as the Molnár János Cave population is isolated for at least 60,000 years from the closest surface populations (Pérez-Moreno et al. 2017), including the one in the directly connected Malom Lake, we still expect that any patterns where the cave population was clearly divergent from the three surface populations would result from adaptation to the cave environment.
Funding Open access funding was provided by Eötvös Loránd University. The study was funded by the National Research, Development and Innovation Office for international cooperation (SNN 125627). GeH also gained support from the Postdoctoral research grant of the National Research, Development and Innovation Office (PD 132041). GB was supported by the ÚNKP-20-4 New National Excellence Program of the Ministry of Innovation and Technology from the source of the National Research, Development and Innovation Office.
Data availability Data are available in the Electronic Supplementary Material.

Declarations
Ethical approval Experiments were performed according to the guidelines of the Hungarian Act of Animal Care and Experimentation (1998, XXVIII, Sect. 243/1998, which conforms to the regulation of animal experiments by the European Union. At the end of the experiment, animals were euthanised and preserved in RNAlater for further research.

Conflict of interest The authors declare no competing interests.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.