Influence of predation risk on the sheltering behaviour of the coral-dwelling damselfish, Pomacentrus moluccensis

Predation is a key ecosystem function, especially in high diversity systems such as coral reefs. Not only is predation one of the strongest top-down controls of prey population density, but it also is a strong driver of prey behaviour and function through non-lethal effects. We ask whether predation risk influences sheltering behaviour of damselfish living in mutualism with branching corals. Host corals gain multiple advantages from the mutualistic relationship which are determined by the strength of damselfish sheltering. Distance travelled by the Lemon Damselfish Pomacentrus moluccensis away from their host colony was measured here as a proxy for sheltering strength and was expected to be shortest under highest predation risk. Predation risk, defined as a function of predator abundance and activity, turbidity and habitat complexity, was quantified at four reef slope sites in Kepulauan Seribu, Indonesia. Damselfish sheltering strength was measured using stationary unmanned video cameras. Small damselfish (< 2 cm) increased their sheltering strength under high turbidity. Predator feeding activity, but not abundance, influenced damselfish sheltering strength. Contrary to our expectations, sheltering behaviour of adult damselfish decreased under high predator activity. While these observations are in line with risk-averse behaviour by juvenile P. moluccensis, they may indicate the presence of sentinel behaviour in the adults of this species. Habitat complexity seemed to be less important as a driver of damselfish behaviour. These counterintuitive results may indicate complex social behaviour and age-specific strategies for predator avoidance.


Introduction
Predation is a key process in terrestrial and marine ecosystems (Steele 1985). Not only does predation act as a top-down regulator of prey population densities, but it is also one of the strongest drivers of prey behaviour and trophic structure of biological communities (Lima and Dill 1990;Ritchie and Johnson 2009;Estes et al. 2010;Boaden and Kingsford 2015). Predator abundance and predation intensity can vary in space and time (Bernstein et al. 1988;Bosiger et al. 2012). While the probability of being preyed upon may be a direct function of predator abundance and activity, predation risk can further be augmented by environmental factors which can influence either predators, prey, or both (De Robertis et al. 2003;Rilov et al. 2007;Wenger et al. 2013).
Predation plays a particularly important role in high-diversity ecosystems (e.g., coral reefs) by maintaining relative levels of species dominance and function (Hixon and Carr 1997). At the same time, species at higher trophic levels are disproportionally affected by human activities such as fishing, with potentially severe consequences for ecosystem functioning (Pauly et al. 1998;Estes et al. 2011).
Overfishing of predators is a common ailment of many coral reefs around the world, yet the consequences of this on ecosystem functions remain poorly understood (Mumby et al. 2012;Fenner 2014;Maire et al. 2016). The release of the corallivorous crownof-thorns starfish (Acanthaster planci, Linnaeus, 1758) from predation pressure, for example, is believed to have contributed to pervasive outbreaks and large-scale coral mortality on Indo-Pacific reefs (Jackson et al. 2001;Baird et al. 2013). Sufficient levels of predation are essential for ecosystem functioning and the maintenance of ecological processes (Lima and Dill 1990;Catano et al. 2014). Carnivorous coral reef fishes (e.g., piscivores) comprise species-rich assemblages of both long-term inhabitants of the reef habitat and transient visitors (Boaden and Kingsford 2015). Changes or spatial differences in predator biomass and predation rates are likely to have consequences for prey abundance and behaviour (Mumby et al. 2012). Quantifying the spatial and temporal variability of predation rates is therefore fundamental to understand ecosystem functioning.
Intuitively, predation risk for coral reef fishes should vary proportionally with predator abundance and activity. However, this relationship is further influenced by environmental factors such as water turbidity and topographic complexity (hereafter referred to as rugosity) (Beukers and Jones 1997;Rilov et al. 2007;Leahy et al. 2011). These factors may have complex effects on the predator's ability to detect prey and/or the level of risk perceived by the prey. High turbidity results generally from the influx of terrigenous sediment into the ocean or resuspension and implies low light availability and poor contrast between objects and the environment (De Robertis et al. 2003). Turbidity may therefore affect the visual perception of both predators and prey (Utne-Palm 2002;De Robertis et al. 2003;Wenger et al. 2013). Whether either predators or prey are at greatest disadvantage in turbid reefs remains however incompletely understood. Predator success is assumed to be higher under turbid conditions, likely due to poor risk perception by prey fishes (Wenger et al. 2013). The alteration of the predator-prey interaction may however depend on prey size, colour and sensory abilities (Utne-Palm 2002). Piscivorous fish, for instance, might be at greater advantage in turbid reefs compared to prey, given that contrast decreases more sharply for big objects than for smaller ones, and thus the detection of larger fish (predators) may be more difficult (De Robertis et al. 2003).
Coral reef rugosity can also affect piscivorous predation rates. Both, positive and negative effects of rugosity are possible for prey and predators. Predators are assumed to be at a disadvantage in highly rugose habitats, where prey detection may be difficult and physical prey refugia are abundant (Beukers and Jones 1997;Figueiredo et al. 2015). In fact, prey have higher survival rates in sites with intermediate (Crowder and Cooper 1982) or high rugosity (Connell and Jones 1991), compared to flatter areas. However, this may not necessarily be the case for territorial and site-attached prey fishes whose location is predictable and constant. For species tightly associated with benthic structures, early visual detection of predators may be obstructed by physical obstacles. Ambushing or highly mobile piscivores may therefore benefit from highly rugose reefs where they may be less likely to be perceived by prey (Rilov et al. 2007). Accordingly, higher rugosity does not always correspond to high damselfish abundance (Bergman et al. 2000).
The level of predation risk may influence predator avoidance behaviour of prey fishes (Lima and Dill 1990), as well as behaviours implicated in their ecological functions (e.g., feeding rates). High predator abundance, for instance, can lead to diminished grazing rates of herbivorous fish (Catano et al. 2014). Such behavioural patterns resulting from predator avoidance can emerge in response to either acute or chronic predation risk (Madin et al. 2010). Certain plan ktiv orous damse lfish es d isp lay stro ng mutualistic relationships with branching corals, living in groups that use the coral branches for shelter to avoid predation and currents (Liberman et al. 1995;Holbrook et al. 2000). Corals benefit from this behaviour resulting in increased growth, survival and reproduction (Liberman et al. 1995;Holbrook et al. 2008). On reefs where heavy fishing pressure has reduced predator stocks, damselfishes might display a weaker sheltering behaviour than on predator-rich reefs. As reduced sheltering effectively decreases the benefits derived by host corals from mutualistic damselfishes, an indirect link between overfishing and coral health might exist (Garcia-Herrera et al. 2017). Predator avoidance behaviours of site-attached fishes may vary in response to the environmental conditions prevailing on a reef, independently of predator abundance, and may be size-dependent (Leahy et al. 2011;Helfman and Winkelman 2015). The size-spectrum theory predicts abundant small-bodied predators and few large-bodied ones (Sogard 1997). Small-bodied predators would only have access to small prey, whereas large-bodied predators would feed on a broad range of prey sizes (Juanes 1994;Juanes and Conover 1994;Scharf et al. 2000), resulting in a higher number of potential predators for smaller prey (Sogard 1997). The combined effect of multiple environmental and intrinsic mediators of predator-prey interactions is less clear, yet crucial for predicting ecosystem dynamics.
The present study investigates whether predation risk influences the sheltering behaviour of the damselfish Pomacentrus moluccensis (Bleeker, 1853). Predation risk was defined here as the probability of damselfish being preyed upon. The latter was assumed to be a function of predator abundance, level of predatory activity, and environmental moderators such as turbidity and rugosity. The distance at which damselfish hovered away from their host colony was used as a proxy for the strength of their sheltering behaviour. We hypothesised that under high predation risk, damselfish would remain within close range of their host coral colony, and under low predation risk they would venture further away from it. We also predicted that smaller damselfish would tend to stay closer to their host colony compared to larger ones because they are susceptible to being consumed by the full predators' size spectrum. We anticipate that clarifying how predation risk affects damselfish sheltering behaviour can help predict the spatial heterogeneity of this mutualistic association, and likely infer its stability.

Study area
This study was conducted in the Kepulauan Seribu, a group of 110 islands located north of Jakarta off the coast of Java, Indonesia (Fig. 1). Most of these islands are flat, small-sized (<10 ha), and relatively close to each other (Madduppa et al. 2013;Baum et al. 2015). Four patch reef sites were selected, aiming to include a range of low-to-high predator abundance. Two of them, namely Pulau Belanda and Pulau Pelangi, were >2 km from the next densely populated island. Pulau Panjang and Pulau Kelapa were inhabited islets. All transects were deployed at the north-west side of the islets, except for Kelapa, where transects were deployed at the southwest side as conditions there were more suitable to the study design than the north side of the island. Fig. 1 Overview of the study region Kepulauan Seribu, northwest Java, Indonesia (location within Indonesia is indicated in the inset map). Study sites were located on the north-west coast of Pelangi, Belanda, Panjang and on the south-west coast for Kelapa

Study species
Our study focused on the most abundant coral-dwelling damselfish in the study area, the Lemon Damselfish, Pomacentrus moluccensis. This medium-sized planktivorous species shows low specificity for the coral host (Pratchett et al. 2012), but is preferentially associated with live corals (Wilson et al. 2008;Hall and Kingsford 2016). Both P. moluccensis and the closely related Pomacentrus amboinensis (Bleeker, 1868) have been well-studied concerning larval dispersion and settlement (Öhman et al. 1998;Holbrook et al. 2000;Munkres et al. 2007;Pratchett et al. 2012), and also used as model species to understand predator recognition strategies indicating quick adaptation to predators (Mitchell et al. , 2013Bosiger et al. 2012;Bosiger and McCormick 2014) and social learning (Bosiger et al. 2012;Manassa et al. 2013Manassa et al. , 2014Hall and Kingsford 2016). The site-attached behaviour of P. moluccensis (Hall and Kingsford 2016) and its wellstudied behavioural response to predation make this an ideal model organism to assess the influence of predation risk on shelter-seeking behaviour.

Quantifying predation risk
Two proxies for direct predation risk for P. moluccensis were quantified in this study, namely predator abundance and predator activity. Both were assessed in each site at midday to coincide with the time of highest activity by P. moluccensis. Predator abundance was quantified using underwater visual censuses along belt transects. At each of the study sites, a total of six 50 m transects were deployed at a constant depth of 4 m parallel to the coastline. A distance of 25 m between consecutive transects was used to ensure independence. All transects were inspected by the same surveyor, who counted all fishes of five free swimming piscivorous families (Labridae, Lethrinidae, Lutjanidae, Nemipteridae, Serranidae, Table S1) within an area of two metres at either side of the transect. Cave-and sanddwelling predators were not counted due to the difficulty of precisely assessing their numbers on a belt transect. Second, levels of predatory activity were quantified using feeding assays at all sites. Dried squid meat cut into circular pieces of 13 mm in diameter was offered as bait, using a modification of a protocol previously used to quantify predation rates in tropical seagrasses (Duffy et al. 2015). Along the transects used to quantify predator abundance and after conducting all fish censuses, a total of 150 baited fiberglass sticks (25 per transect) were fixed to the substratum with squid bait attached to their upper end. Sticks on each transect were separated approximately 1 m from each other (Duffy et al. 2015), which resulted in six assays stretching over a distance of about 500 m. Snorkelers counted the number of baits remaining 10, 30 and 60 min after deployment, with the number of remaining bait assumed to be inversely related to predator activity. One to two GoPro Hero 3 cameras were placed in the vicinity of each transect to record the species responsible for the removal of bait on a subset of the assay.

Quantifying the strength of sheltering behaviour in Pomacentrus moluccensis
The distance swum by damselfish away from their host coral was considered here as an adequate proxy for the strength of their sheltering behaviour, and was quantified using video surveys. Because P. moluccensis shows low host specificity (Pratchett et al. 2012), ten branching coral colonies within 25 m of the transects were selected for video surveys irrespective of their species. Corals were considered only if they were completely alive and hosted at least two damselfish, but preference was given to colonies with more than two individuals. To facilitate the visualisation of the damselfish's position with respect to their host coral, we focused on isolated colonies that provided footage with a water column background. Filmed corals were at least 15 m apart from each other to minimise potential mutual influence and exchange of individuals. The cameras were aligned with the upper edge of the corals and fixed about 1.5 m away, recording the widest side of the colony if it was not equilateral. GoPro Hero3 cameras were set to record for at least 30 min in the absence of divers or snorkelers, immediately after a 30 cm size reference was held next to the coral. Upon camera retrieval, the host coral was measured and the number of damselfishes counted, classifying them into one of three size categories: small (< 2 cm), medium (2 -4 cm), or large (> 4 cm). The host coral species as well as adjacent colonies (< 5 m away) were photographed for posterior identification.
Posterior video analysis focused on twenty minutes of diver-free video and used 'VLC media player' (VideoLAN®, Version 2.2.1, Apr 2015).
Distances swum by damselfish away from their host coral were assessed with a digital screen ruler (MB-Ruler) (Bader 2016). A single surveyor extracted the data from all videos. The first step consisted of confirming the number and sizes of damselfish present in the colony and scaling each video based on the 30 cm size reference shown initially. Data were then collected on still images captured every minute. On these images, the distances between the graphic centre of each damselfish and the nearest point of the coral, as well as the total length (TL) of each fish, were measured. Given that during video analysis, the 3D space in which damselfish move was projected onto a 2D plane (i.e. computer screen), measured distances may not have always reflected the exact distance between a fish and its host coral. This was particularly true for fish swimming towards the camera. These fish could have thus hovered at a greater distance from their coral than perceived by the observer. As a result, certain distances could have been underestimated during video analysis. Considering that this study focused on assessing relative rather than absolute distances, and that damselfish are randomly distributed around their host, the procedural bias was likely homogeneous across sites. Conclusions based on inter-site comparisons were therefore assumed to be minimally affected by this technical limitation. Furthermore, the position of the cameras filming the corals was chosen at random in relation to the spatial distribution of the fish. Even if damselfish groups did not distribute randomly in space, the relatively high sample size (20 snapshots per group and 10 groups per site) should have cancelled out this effect. Certain damselfish recorded in front of the coral were indistinguishable from others sheltering among the coral branches. For this reason, distances in both cases were measured as zeroes. Distances swum by fish that went outside the camera's field of view or behind the coral colony were not measured.

Quantifying environmental moderators of predation risk
Environmental factors likely to influence behaviour of predators and prey, and therefore mediate predation risk, were assessed at each site within a 5 m radius around the host coral colonies selected for video surveys. Rugosity was estimated visually on a scale from 0 (virtually flat) to 4 (highly rugose) based on the amount of vertical substrate, hideouts, and caves (Polunin and Roberts 1993). As a proxy of turbidity, horizontal visibility was quantified using a Secchi disk and a measuring tape deployed by two snorkelers at a depth of 1 m.

Spatial differences in predation risk
Predation risk was considered a function of predator abundance, predator activity, rugosity and turbidity. To confirm that predation risk varied spatially, we analysed each of these variables separately for a site effect. In the case of turbidity and rugosity, which were non-normally distributed, comparisons were conducted using Kruskal-Wallis tests. Post-hoc pairwise comparisons with a Bonferroni correction were conducted to identify differences between particular sites whilst avoiding Type I errors arising from multiplicity (Bretz et al. 2011). The p-value for acceptance of results as significant was thus adjusted to 0.0083. The total number of predators was calculated per transect, and differences among sites were tested using a generalised linear hypothesis test (glht) with the Tukey method for multiple comparisons (Bretz et al. 2011).
The feeding assays were assessed using survival analysis (Pyke and Thompson 1986) with the 'survival' R package (Therneau 2015). The amount of remaining baits per transect at different times was recorded to calculate survival times per transect. Baits that were recorded missing 10 and 30 min into the feeding assay were assumed to have been consumed after the snorkeler's previous observation. Kaplan-Meier curves of bait survival were then computed per site (Kalbfleisch and Prentice 2002;Rich et al. 2010). The baits remaining at the end of the assay were regarded as having an unknown survival time, i.e. the assay data contained censored data (Pyke and Thompson 1986). We therefore used a Log-Rank multiple comparison with Bonferroni correction (p-value = 0.0083) to compare the survival rates in order to test whether sites differed in predation risk (Pyke and Thompson 1986;Rich et al. 2010). As a proxy of predator feeding activity, the time point at which 50% of the bait remained unconsumed was derived for each transect from an exponential curve, which was the best fit of the Kaplan-Meier curve.

Influence of predation risk on damselfish sheltering behaviour
To determine the relative influence of predation risk and the environmental moderators of predation risk on the strength of damselfish sheltering behaviour (measured as the average distance swum by damselfish away from their host coral), a multiple regression was used. The complete model considered turbidity, rugosity, predator abundance and predator activity. We tested for collinearity among these factors by computing the variance inflation factor (VIF) of the complete model with the Bcar^package in R (Fox et al. 2017). A VIF value of 3 was chosen as a threshold for collinearity (Zuur et al. 2010). No factor exceeded the VIF threshold of 3, thus indicating there was no collinearity. Forward stepwise model selection was performed using the 'leaps' package in R (Lumley 2017), choosing parameters that decreased the AIC value of the model. A separate model was run for each damselfish size category to investigate whether prey size influences the relative roles of predation risk and the environmental factors that mediate predation risk as drivers of sheltering behaviour strength. We verified that these models conformed to normality and homoscedasticity by inspecting diagnostic residual plots (Zuur et al. 2010). The complete as well as the optimal model met the assumptions of normality and homogeneity of variances.

Intrinsic drivers of damselfish sheltering behaviour
The average strength in the sheltering behaviour of damselfish groups may be driven by the number of individuals within the group, due to damselfish performing optimal foraging, as well as by their mean total length, due to higher predation risk for smaller individuals (Townsend and Winfield 1985;Sogard 1997). To test this, we ran two separate linear models. Mean sheltering strength of all individuals per group was fitted as a function of the number of individuals within the group in a first model, and as a function of mean size category of the group in a second model. Lastly we verified that these models incurred no significant divergence from normality or heteroscedasticity (Zuur et al. 2010).

Spatial differences of predation risk
Environmental factors considered here as moderators of predation risk (i.e. turbidity and rugosity) were significantly different among sites (Fig. 2a, b). Total predator abundance was significantly lower at Pelangi than at all other sites (glht Tukey, p < 0.025). Predator feeding activity also differed among sites. Between 29.3 and 81.3% of bait was left 10 min after assay deployment, and 5.3-56% after 30 min of assay deployment. Unconsumed bait at the end of the experiment ranged from 0 to 34.7% (Fig. 2d). The highest predator feeding activity occurred at Belanda (p < 0.001), where visibility was also highest (20 m), and rugosity as well as predator abundance was moderate. Feeding activity was marginally higher at Kelapa compared to Pelangi (Log-rank, p = 0.054). The consumption of 86 baits was documented by cameras. Camera recordings show that baits were mainly consumed by Pentapodus trivittatus (Bloch, 1791), Thalassoma lunare (Linnaeus, 1758) and Cephalopholis microprion (Bleeker, 1852) (24.4%, 24.4% and 20% of baits, respectively).

Influence of predation risk on damselfish sheltering behaviour
Predator activity had a negative and significant effect on sheltering strength of entire damselfish groups (p < 0.001, Fig. 3a, Table 1), whereas turbidity had a non-significant influence (p = 0.0883). Neither rugosity nor predator abundance were included as significant predictors in any of the final models. The sheltering behaviour of large damselfish (> 4 cm) was weaker when predator activity was higher, as it was for medium-sized individuals (2-4 cm; p < 0.001 and p = 0.0275, respectively, Fig. 3b, Table 1). Sheltering behaviour was positively influenced by turbidity for smallbodied damselfishes (< 2 cm; p = 0.0227; Table 1) but unaffected by predator activity. The intercept of the model for this size category was not significantly different from zero, showing that under turbid conditions, small individuals hid among the coral branches.

Intrinsic drivers of damselfish sheltering behaviour
Sheltering behaviour of damselfish groups was strongly and negatively correlated with the mean size category of the group (p < 0.001, Fig. 4a). The largest damselfish (> 4 cm) hovered about twice as far away from the host corals compared to medium sized (2-4 cm) and small (< 2 cm) individuals. The mean number of damselfish associated with a single host coral ranged between two and 17 but had no significant influence on sheltering behaviour of the group (Fig. 4b).

Discussion
The present study investigated whether predation risk influences the strength of the sheltering behaviour of damselfish living in mutualism with branching corals. We demonstrate that sheltering behaviour is indeed influenced by predation risk, which is consistent with other studies indicating threat sensitivity and risk assessment by damselfish (Helfman 1989;McCormick and Manassa 2007;Leahy et al. 2011;Bosiger et al. 2012;Helfman and Winkelman 2015;Chivers et al. 2016). Predator activity rather than predator abundance was the main factor influencing damselfish sheltering behaviour. Similar observations were made for gobies, where a combination of visual and chemical (extract of an injured conspecific) alarm cues yielded the maximum anti-predator response (McCormick and Manassa 2007). The combination of these cues is associated to predator activity rather than mere abundance, given that an injury (or death) of a conspecific prey is required to provide the chemical cue required for the maximal response. However, the response to predation activity was the opposite of that identified in earlier studies. Instead of increasing sheltering strength, predation activity seemed to decrease it. These results are counterintuitive since we expected sheltering behaviour would be strongest where predation activity, and therefore predation risk, were highest (Liberman et al. 1995). Turbidity influenced only the smallest damselfish, whereas rugosity was less important.
Our models demonstrate that damselfish sheltering behaviour is strongly influenced by the level of predator activity. Damselfish sheltering behaviour was weakest where predator activity was highest. This was particularly true for the largest damselfish. Sheltering behaviour varied regardless of predator activity for the smallest damselfish, but was significantly and positively influenced by turbidity. This could indicate that the smallest damselfish perceive high turbidity as an enhancer of predation risk, which is consistent with studies indicating a higher mortality and decreased foraging in high turbidity conditions for other damselfish species (Leahy et al. 2011;Wenger et al. 2013).
Damselfish size had a strong effect on their sheltering behaviour. Groups of larger individuals displayed weaker sheltering behaviour than groups of smaller individuals. The size-spectrum theory predicts that smaller predators are more abundant than larger ones, which indicates that smaller prey may have a larger number of potential predators (Sogard 1997). Predator success as well as prey size selection depend on various parameters, such as mouth gape, that are largely influenced by predator size (Juanes and Conover 1994). Predators tend to target smaller prey items because this results in higher number of successful predation attempts (Juanes 1994;Juanes and Conover 1994). This further increases the potential number of predators of small-sized damselfish (Sogard 1997) and could explain the size-dependent sheltering behaviour observed here. Studies on the closely-related bicolor damselfish Stegastes partitus The direction of the correlation as well as the p-value, or else discard (Disc.) in the optimal model, are shown for predator activity, predator abundance, turbidity and rugosity, for entire damselfish groups (i.e. including all sizes) and per size category. AIC values of the model containing all parameters and the optimal model are compared Fig. 3 Influence of predator activity on damselfish sheltering behaviour. a significant relationship between the sheltering behaviour of entire damselfish groups (i.e. including all sizes) and predator activity (R 2 = 0.3172, p < 0.001); b individual sheltering behaviour of each size category of damselfish (small: triangles, medium: crosses, large: circles) in relation to predator activity (small: ns; medium: R 2 = 0.147, p = 0.0275; large: R 2 = 0.447, p < 0.001) (Poey, 1868) indicate that threat sensitivity depends on both the damselfish' and predator's body size (Helfman 1989;Helfman and Winkelman 2015). The species Stegastes planifrons (Cuvier, 1830), one of the most common damselfishes on Caribbean reefs, shows a similar, size-dependent response to predation risk. Smaller individuals display high sensitivity to threat, gradually reacting to danger, while larger individuals show a more indifferent behaviour, not always reacting to potential predators (Helfman and Winkelman 2015). According to the optimal foraging theory, larger groups of damselfish may prompt individuals to swim longer distances in the competitive pursue of plankton and optimal nutrition (Townsend and Winfield 1985). Group size seems, indeed, to act on the foraging of Dascyllus aruanus (Linnaeus, 1758) by diminishing the mean size of food items consumed and inducing an upstream swimming of highranked individuals (Forrester 1991). In the present study, the number of individuals within P. moluccensis groups did not influence shelter seeking. This might be due to the relatively small group sizes, as the mean number of individuals per group was 5.9 here compared to 7.9 in a similar study (Forrester 1991). This smaller group size might result in lower competition between individuals and thus a less pronounced optimal foraging effect. Other studies of damselfish show that no competition occurs between individuals within a group and that smaller individuals even increase energy uptake with higher density of conspecifics, which might also account for the observed results (Booth 2004). Contrary to what we expected, damselfish in general exhibited weaker sheltering behaviour under high levels of predatory activity. We expand on a number of possible explanations for this unexpected observation. Damselfish subject to high levels of predatory activity may tend to position themselves relatively high up in the water column, thus maximising their field of vision and predator detection. Coordinated vigilance has been observed for pairs of feeding rabbitfishes (Brandl and Bellwood 2015). P. moluccensis identifies usual predators rapidly and is capable of learning to recognise new predators and their activity patterns (Mitchell et al. , 2013Bosiger et al. 2012). An increasingly cautious behaviour under conditions of high predator activity would allow damselfish to quickly react to danger, and individuals exhibiting this behaviour may have a direct survival advantage. A sign of this potential learned predator response is that only large, likely more experienced individuals exhibited weak sheltering behaviour under high predation risk. Younger, thus smaller damselfish of the species Acanthochromis polyacanthus (Bleeker, 1855) have shown only a diffuse response to environmental moderators of predation risk (Leahy et al. 2011), likely because they have not yet developed predator awareness . Nonetheless, small individuals may still benefit from the behaviour of larger individuals in the group because they are capable of social learning and could react to the retreat of the adults even if they are naïve towards the predator itself, thus maximising their survival (Manassa et al. 2013). In the presence of adult conspecifics, juvenile damselfish have an increased growth rate due to higher food intake (Booth 2004). Since anti-predator behaviours (e.g. surveillance, sheltering) require a time and energy investment (Lind and Cresswell 2005), juvenile damselfish might profit from the behaviour of adult individuals by Fig. 4 Influence of mean size and number of individuals in damselfish groups on sheltering behaviour. a Significant relationship between sheltering behaviour and mean size of damselfish (R 2 = 0.3114, p < 0.001); b no significant relationship between the number of individuals in a group and sheltering behaviour allocating less time to surveillance themselves and more time to feeding. Genetic relatedness among individuals within a group has been shown for several coraldwelling damselfish (Öhman et al. 1998;Buston et al. 2009). Individuals within a group do in fact maintain a strong hierarchical and social structure (Forrester 1991). This increases even further the potential fitness gain from vigilance behaviour by adults, since they would promote their own survival as well as survival of conspecifics, contributing to a stable altruism (Foster et al. 2006). By maintaining vigilance for the group, adult individuals increase their food intake and presumably their fitness, thus resulting in carryover maternal effects during gamete production, leading to increased fitness in offspring (McCormick and Gagliano 2008;Bonin et al. 2009). Therefore, it seems reasonable to suggest that our findings support the presence of sentinel behaviour in coral reef damselfishes. Sentinel behaviour is displayed when certain individuals of cooperative animal groups concentrate on scoping the environment for predators and warn conspecifics in case of an attack. Generally, sentinel behaviour increases with the perception of predation risk, and it is displayed by elder and more experienced individuals (Ridley et al. 2010). These conditions are met in the present case.
The observed pattern may also have been caused by a confounding factor simultaneously influencing both, predator activity and sheltering behaviour. Shallow incoming plankton clouds may, for instance, prompt damselfish to decrease sheltering behaviour in favour of foraging while simultaneously increasing levels of activity of predators attracted to exposed planktivores. Data on planktonic food items was not available to discard this possibility in the present study, and should be included in future studies of damselfish sheltering behaviour.
Alternatively, larger individuals may be more prone to display courtship and mating behaviour, prompting them to swim further from their host colony (Shpigel and Fishelson 1986). However, this would not explain the observed negative relationship to predator activity.
Both hypotheses underlying the observed behavioural patterns in larger individuals (i.e. foraging versus sentinel behaviour) could be easily verified by investigating the relationship among damselfish foraging activity, sheltering behaviour and levels of predator activity. If in response to increased food availability, foraging intensity of damselfish increases along with predator activity while damselfish sheltering behaviour weakens, a confounding factor (i.e. food availability) may be responsible for the observed reduction in sheltering. Alternatively, if sheltering behaviour of damselfish decreases with predator activity without an increase in feeding when food availability is high, then damselfish likely behave cautiously rather than engaging in foraging, and the weakened sheltering behaviour is linked to sentinel behaviour. Genetic studies of the relatedness between individuals within a group could back up the results in the second case.
In our study, predator activity was mainly accounted for by mesopredator species like P. trivittatus, T. lunare and C. microprion. Although these species might appear as low risk-inducing predators, mesopredators and particularly these species do consume coral-dwelling damselfish (Bonin et al. 2009;Holmes and McCormick 2010). An absence of apex predators in the feeding assay would therefore not explain the decreased sheltering behaviour observed here for large damselfishes where predator activity was high.
Here we considered the mean distance between damselfish and their host colony as an appropriate indicator of sheltering behaviour, but there are caveats to this assumption. Similar studies used the foraging rate (Hall and Kingsford 2016), the time spent inside and outside the coral (Rilov et al. 2007), or the maximum distance covered by the damselfish (Forrester 1991). Our proxy does however provide a reasonably precise metric of sheltering behaviour since it is influenced by both the maximum distance covered and the time spent inside and outside the coral. The quality of the footage obtained here was insufficient to quantify foraging rates. Future studies would however benefit from integrating mean distance of damselfish, food availability, and foraging rates, especially if there is an interest in detecting changes in the consistency of their ecological functions. Further precision could be achieved using stereoscopic cameras. A 2D projection of a 3D image results in a certain level of underestimation of the distance between damselfish and their host coral. However, the relative distance was sufficiently informative in our study.
Whilst sheltering, damselfish move their fins constantly, thus ventilating the small spaces between the coral branches, reducing hypoxia forming as a result of coral respiration at night (Goldshmid et al. 2004). They also increase available ammonium due to excretion (Holbrook et al. 2008) and defend their host from corallivores (Chase et al. 2014). Damselfish potentially contribute to the dissipation of excess oxygen produced by the photosynthesis of the corals. Some soft corals are known to produce small-scale currents by movements of their polyps to evacuate excess oxygen (Kremien et al. 2013). Damselfish are likely to produce similar effects on hard corals through the movement of their fins, and sheltering has been shown to contribute to increased photosynthesis (Garcia-Herrera et al. 2017). Ultimately, the mutualistic relationship with damselfishes results in increased growth, survival and reproduction of the host coral (Liberman et al. 1995;Holbrook et al. 2008). The magnitude of the benefits derived by host corals from sheltering damselfish may be determined by the consistency of their mutualistic association. Benefits may be decreased when damselfish spend less time sheltering. We originally hypothesised that low predation risk would result in weakened sheltering behaviour by P. moluccensis, thus leading to an alteration of the ecological link between damselfish and their host coral. Our results, however, show that the influence of predation risk on behaviour does not always conform to this expectation, underlining the complexity of predicting how changes in predation risk will influence the relationship between damselfish and their host coral. High levels of predation are generally taken as indicators of a healthy ecosystem (Parsons 1992;Jackson et al. 2001;Shurin et al. 2002;Boaden and Kingsford 2015). Although high levels of predator activity may decrease the strength of damselfish sheltering in their host coral and potentially decrease the benefits host corals can derive from mutualistic damselfish, this unexpected link between predation risk and sheltering behaviour does not appear widespread among coral-associated damselfish. A negative effect of high levels of predation on overall ecosystem health is therefore discarded. Certainly, behavioural responses of damselfishes to predation will vary among species or even populations of the same species. Further research is thus needed to accurately predict these dynamics.