Introduction
Murid rodents are among the most abundant granivores in tropical forests and can influence seed fate by consuming or destroying seeds and by dispersing them to different microhabitats (Wells & Bagchi Reference Wells and Bagchi2005). Murids can influence species-specific seed survival (Hautier et al. Reference Hautier, Saner, Philipson, Bagchi, Ong and Hector2010) if foraging decisions are associated with dietary preferences based on seed morphology (Dirzo & Mendoza Reference Dirzo and Mendoza2007, Cramer Reference Cramer2014). Predator-induced changes in granivory could therefore affect seed survival (Laundré et al. Reference Laundré, Hernández, Medina, Campanella, López-Portillo, González-Romero, Grajales-Tam, Burke, Gronemeyer and Browning2014, Bestion et al. Reference Bestion, Cucherousset, Teyissier and Cote2015) with disproportionate effects on seeds with particular morphological traits (Xiao et al. Reference Xiao, Wang, Harris and Zhang2006, Forget & Jansen Reference Forget and Jansen2007, Bricker et al. Reference Bricker, Pearson and Maron2010). Most of what we know about these processes, however, comes from temperate systems; the indirect effects of predators on seed survival via impact on granivorous rodents in tropical forests are much less well known.
Predator scent can elicit physiological and behavioural changes in prey (Apfelbach et al. Reference Apfelbach, Blanchard, Blanchard, Hayes and McGregor2005, Takahashi et al. Reference Takahashi, Nakashima, Hong and Watanabe2005, Sündermann et al. Reference Sündermann, Scheumann and Zimmermann2008). Rodents often use scents (e.g., from urine) to gauge predation risk (Rosell Reference Rosell2001, Apfelbach et al. Reference Apfelbach, Blanchard, Blanchard, Hayes and McGregor2005, Ferrero et al. Reference Ferrero, Lemon, Fluegge, Pashkovski, Korzan, Datta, Spehr, Fendt and Liberles2011, Bytheway et al. Reference Bytheway, Carthey and Banks2013) and may spend less time eating seeds or exhibit increased vigilance in riskier microhabitats (Lima & Dill Reference Lima and Dill1990). In some cases, the indirect effects of predation risk may be as strong as direct reductions in prey abundance for plant survival (Hernández & Laundré Reference Hernández and Laundré2005, Preisser et al. Reference Preisser, Bolnick and Benard2005). Though the effects of predator scent on prey foraging behaviour have been documented, whether predator-induced changes in foraging translate into altered seed removal, and survival has not been well explored in tropical systems.
Methods
We used motion-triggered, infrared cameras (Reconyx HC500) and field experiments to test whether simulated predator presence affected visitation and seed removal rates of murid rodents in Danum Valley Conservation Area (DVCA; N5.10189°, E117.688°), Sabah, Malaysian Borneo. DVCA (438 km2) is the largest primary dipterocarp forested area in Sabah and is a Class 1 (‘Totally protected’) Forest Reserve with prohibited resource extraction (Hazebroek et al. Reference Hazebroek, Adlin and Sinun2012). Murid rodents (> 27 species) are the main mammal granivores, mostly from the genera Rattus and Maxomys (Phillips & Phillips Reference Phillips and Phillips2016). Bearded pigs (Sus barbatus), lesser (Tragulus napu) and greater (T. javanicus) chevrotains, and red (Muntiacus muntjac) and Bornean yellow (M. atherodes) muntjac also occur in the system and eat seeds. All vertebrate species detected in our study area are listed in Supplementary Material Figure S1.
We simulated predator presence using bobcat (Lynx rufus) urine gel (Bobcat P-Gel, PMart, Sandy Point, ME) at 10 paired experimental stations ca. 500 m apart (> the average daily movement of rodents in our site; Wells et al. Reference Wells, Pfeiffer, Lakim and Kalko2006). Each station contained two 1 × 1 m quadrats (treatment and control) spaced 50 m apart (i.e., within broadly similar microhabitat conditions), with one camera trap overlooking each quadrat. We placed 5 g of urine gel in the treatment quadrats at several locations (Figure 1) and left the control quadrats unmodified. Cameras were active 24 hours per day and set to take 10 photos in rapid-fire succession upon triggering, with high trigger sensitivity and no quiet period between triggers.
At each quadrat, we placed four seeds from each of four plant species, for a total of 16 seeds, ∼ 20 cm from plot edges and from adjacent seeds. Seed order was random, but conspecifics were not placed next to each other. Seeds were set out for four days, after which removal was determined by our field team counting the number of seeds missing from each plant taxon, with removal attributed to all granivores. We used camera trap photograph to determine seed removal by rodents in particular, with removal events considered independent when a photograph showed a rodent handling or removing seed and subsequent photographs showed a rodent with no seeds in its feet or mouth.
We used seeds from four plant genera varying in seed size: pumpkin (Cucurbita pepo; mean seed mass = 0.25 g), longan (Dimocarpus longan; 0.78 g), rambutan (Nephelium lappaceum; 1.00 g), and oil palm (Elaeis guineensis; 6.61 g). E. guineensis is native to Africa, but rodents are known to prey on seeds of this species in Southeast Asia (Buckle et al. Reference Buckle, Chia, Fenn and Visvalingam1997). C. pepo is native to North America, but Cucurbita spp. are known to be eaten by murids including Rattus spp. and Mus spp., both of which occur in our system (de Guia & Quibod Reference de Guia and Quibod2014). Dimocarpus and Nephelium are in the Sapindaceae family, whose seeds are eaten by murid rodents in Borneo (Blackham & Corlett Reference Blackham and Corlett2015).
Data analysis
We used a generalised linear mixed-effects models (GLMM) with Poisson error distributions to compare the number of independent murid detections between treatment and control sites. We used camera trap station ID as a random intercept (a factor with 10 levels representing each treatment and control pair) to account for microhabitat differences among camera trap stations that may have affected small-scale habitat selection (Bowers & Dooley Reference Bowers and Dooley1993, Mayor et al. Reference Mayor, Schneider, Schaefer and Mahoney2009). The number of independent murid detections was estimated from camera trap photos; detections were considered independent if they were separated by any non-murid species or were >5 minutes apart.
Next, we used a GLMM with binomial error distributions to model overall seed removal (by all granivorous taxa) as a function of two explanatory variables: the number of murid rodent detections and treatment. Here, the fate of each seed at the end of our study period was categorised as present or absent (i.e. seed removal, determined by our field team), and paired camera station ID was modelled as the random intercept.
To determine if seed removal by rodents in particular was influenced by predator scent, we used GLMMs with a binomial error distribution to model the overall proportion of seeds removed by rodents as a function of experimental treatment. We assessed overall removal as we could not distinguish between seed species being removed or handled by rodents in camera trap photographs. As before, camera station ID was used as the random effect. The presence of large-bodied granivores (i.e., bearded pigs) could have affected murid visitation to seed quadrats and seed removal (Keesing Reference Keesing1998), so we also ran a separate GLMM with the same model structure but excluded data from all quadrats visited by bearded pigs.
Finally, we examined if seed removal varied among seed species. We tested whether removal (by all granivores) differed among seed taxa (Cucurbita, Dimocarpus, Nephelium, and Elaeis) and whether predator scent had any effect on removal. First, we assessed seed removal (a binary response: seed present or missing from the quadrat at the end of the study) as a function of seed species and treatment. We used a post hoc Tukey test to compare seed removal between seed species. Finally, we ran separate GLMMs for each plant taxon, modeling individual seed removal caused by any animal (missing versus present) for each species, as a function of treatment and of murid rodent detections. For all models using species-specific removal as the response, we modelled overall removal (i.e., due to murids or non-murids), rather than seed removal by rodents per se because we were unable to distinguish seed species in the camera trap photos. For both types of analysis of species-specific seed removal, paired camera station ID was included in models as the random intercept.
All regression models were run using the lme4 package (Bates & Maechler Reference Bates and Maechler2011) in R 3.3.2 (R Core Team 2018). The R package Performance was used to check for multicollinearity (Lüdecke et al. Reference Lüdecke, Ben-Shacar, Patil, Waggoner and Makowski2021). Only covariates with low or moderate correlation were retained in analyses. Comparisons of seed removal between species were performed using the emmeans package in R (Lenth Reference Lenth2020).
Though our model, predator species is not native to the study region, bobcats are potent predators of rodents elsewhere (Neale & Sacks Reference Neale and Sacks2001, Thornton et al. Reference Thornton, Sunquist and Main2004), and their urine is a valid proxy for rodent predation risk for several reasons. Rather than responses conditioned to specific predators, prey species typically have innate, generalized responses towards predators in general (Hayes et al. Reference Hayes, Nahrung and Wilson2006). Predator scent recognition is innate in murid rodents, and bobcat urine has been used in experiments measuring changes in rodent foraging (Orrock et al. Reference Orrock, Danielson and Brinkerhoff2004) and behaviour (Fendt Reference Fendt2006). Completely predator-naive laboratory rats (which, like the species we studied, are from the Muridae family but not sympatric with bobcats anywhere) have been shown to respond defensively to bobcat urine despite having spent their whole lives in captivity (Fendt Reference Fendt2006). Lab rat responses may differ somewhat from those of murids in nature, but studies with captive animals suggest they exhibit broadly similar responses to those under natural conditions (Apfelbach et al. Reference Apfelbach, Blanchard, Blanchard, Hayes and McGregor2005). Non-sympatric predator scents elicit behavioural changes in other murid species as well (Bramley & Waas Reference Bramley and Waas2001, Dielenberg & McGregor Reference Dielenberg and McGregor2001, Carthey & Banks Reference Carthey and Banks2016). Therefore, we were satisfied that bobcat urine could serve as a valid ‘model predator’ scent cue, given that urine from predators native to the study area (e.g., leopard cats, Prionailurus bengalensis) was not available.
Results
Murid visits to seed plots were similar between control and treatment sites (βtreatment = 0.270, SE = 0.247, df = 17, P = 0.261). Overall, seed removal by all granivores (i.e., murids and non-murids) was higher in predator urine quadrats (βtreatment = 0.662, SE = 0.263, df = 317, P = 0.036
On average, murids removed 1.7 seeds (SD = 2.70) across all quadrats, but the proportion of seeds removed was not significantly influenced by predator scent (βtreatment = 0.652, SE = 0.411, df = 17, P = 0.113; Figure S2). In the absence of bearded pigs, average seed removal (by all other taxa, not just rodents) was 2.54 seeds (SD = 3.33) per site and removal by murids did not differ between treatment and control sites (β = 0.690, SE = 0.475, df = 8, P = 0.146).
Seed removal by all granivores (including camera sites visited by bearded pigs) varied across plant taxa. Removal of seeds of all other plant taxa was more frequent than that of Elaeis (β Cucurbita = 5.298, SE = 0.752, df = 313, P < 0.001; β Dimocarpus = 3.177, SE = 0.517, df = 313, P < 0.001; β Nephelium = 1.973, SE = 0.467, df = 313, P < 0.001). Results for comparisons of removal between species are shown in Table S1. Removal of Cucurbita seeds was higher than that of all other species (Figure 2). Rodent visits were not significantly related to removal rates in any of the genera (P > 0.06; Table S1). Seed removal in Elaeis was marginally significantly associated with rodent visits (βMurid = 0.512 SE = 0.264, df = 76, P = 0.052, Table S2) but, after excluding an outlier, there was no apparent effect (βMurid = 0.269, SE = 0.375, df = 72, P = 0.474). Dimocarpus was the only taxon for which seed removal was affected by predator scent; removal was higher at sites with predator urine (βtreatment = 2.002, SE = 0.878, df = 76, P = 0.023, Table S2).
Discussion
Predator scent did not influence murid rodent visits to camera stations. The lack of difference in rodent visitation between treatments suggests that the perceived level of predation risk associated with our deployment of bobcat urine gel was insufficient to cause granivores to avoid the sites. Other studies have demonstrated murid foraging responses to predator scent, including from non-native predators (Wolff Reference Wolff2004, Apfelbach et al. Reference Apfelbach, Blanchard, Blanchard, Hayes and McGregor2005, Ramp et al. Reference Ramp, Russell and Croft2005, Carthey & Banks Reference Carthey and Banks2016). The lack of treatment effects that we observed could simply demonstrate that the perceived level of predation risk was not sufficient to alter foraging behaviour, such that food rewards outweighed potential risks of foraging in ‘predator’ plots. Trait-mediated predation effects are known in some systems (Schmitz et al. Reference Schmitz, Krivan and Ovadia2004, Preisser et al. Reference Preisser, Bolnick and Benard2005) but may be far from ubiquitous. Indeed, predators often have no detectable influence on prey distributions (Brodie & Giordano Reference Brodie and Giordano2013). Even when prey do avoid predators spatially, avoidance movements are often very temporary, rendering changes in overall foraging patterns negligible (Kauffman et al. Reference Kauffman, Brodie, Jules and Url2013, Brodie et al. Reference Brodie, Aslan, Rogers, Redford, Maron, Bronstein and Groves2014). It is possible that rodents only respond to predator scent immediately after application and that the perceived risk decreases as the urine is masked by other scents in the environment (McFrederick et al. Reference McFrederick, Fuentes, Roulston, Kathilankal and Lerdau2009, Bytheway et al. Reference Bytheway, Carthey and Banks2013). However, rodents visited our stations throughout the study period (including immediately after gel application), suggesting that any temporal changes in urine gel potency did not influence rodent behaviour. Changes in granivore responses to predator odors could be context-dependent, varying with extrinsic factors such as temperature and weather (Herman & Valone Reference Herman and Valone2000, Orrock & Danielson Reference Orrock and Danielson2009). Finally, food availability (i.e., seed density) could influence animal willingness to forage as well as how much food they would consume in habitats perceived as risky (Brown et al. Reference Brown, Morgan and Dow1992). Overall, rodent food availability in our system may have been high as our study took place during a dipterocarp tree masting event providing known rodent food (Phillips & Phillips Reference Phillips and Phillips2016).
Seed removal by rodents was unaffected by predator urine, but removal by all granivores combined varied with seed species and may have been influenced by seed traits. Granivores may prefer certain seed types depending on energy content (Xiao et al. Reference Xiao, Wang, Harris and Zhang2006), handling time, body-to-seed-size ratio (Muñoz & Bonal Reference Muñoz and Bonal2008), toxicity, or morphology (Myster & Pickett Reference Myster and Pickett1993, Hulme & Benkman Reference Hulme, Benkman, Herrera and Pellmyr2002). We found that removal rates varied among plant genera, suggesting that seed traits might influence granivore dietary preference. Although larger seeds provide more energy (Charnov Reference Charnov1976, Mack Reference Mack1998, Brewer Reference Brewer2001), the largest seeds in our experiment (Elaeis) were removed least often, while the smallest (Cucurbita) were removed most often. Therefore, granivores in DVCA might prefer smaller seeds because of reduced handling times (Dirzo & Mendoza Reference Dirzo and Mendoza2007, Muñoz & Bonal Reference Muñoz and Bonal2008, Wang et al. Reference Wang, Ye and Chen2013). If removal was driven by seed size, low Elaeis removal could point to a dietary size threshold for seed predators, whereby granivory has a greater negative effect on seeds below a certain weight (Dirzo & Mendoza Reference Dirzo and Mendoza2007, Perez-Ramos et al. Reference Perez-Ramos, Garcia-De La Cruz and Gomez-Aparicio2017). Also, granivores may have perceived the added time spent handling large Elaeis seeds in high-risk treatment quadrats as too high relative to energetic gains (Lima & Bednekoff Reference Lima and Bednekoff1999, Dirzo & Mendoza Reference Dirzo and Mendoza2007). Foraging for small seeds may be perceived as less risky, even in the presence of predators.
Conclusion
The use of predator scent to simulate predation risk can reveal how prey species perceive risk in the environment and how any subsequent changes in foraging behaviour might affect plants. Large-scale experiments could provide more insight into how widespread loss of top predators might indirectly affect plant communities. Future studies should use experimental exclosures to identify prey-specific changes in foraging behaviour. We also recommend placing seeds at a range of densities to investigate the influence of food availability on foraging decisions in high risk areas. Finally, evaluating prey responses towards a range of predator species could help determine whether rodents show similar responses towards other non-native predators.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S0266467422000050
Acknowledgements
This study was conducted with permission from Yayasan Sabah, the Sabah Biodiversity Council, and the Danum Valley Management Committee. This work was part of the Southeast Asian Rainforest Research Partnership (SEARRP) and was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) post-graduate doctoral award (PGS-D) to AG and NSERC Discovery and Canadian Foundation for Innovation grants to JB. We are grateful for logistical and field assistance from G. Reynolds, A. Karolus, S. Espinola, S. Mang, T. Rachinski, S. Havalic, and T. Ferrel. We thank two anonymous reviewers for their feedback.
Data Availability Statement
Data are available on Figshare.com (https://figshare.com/authors/Alys_Granados/4108249).
Competing Interests
The authors have none to declare.