Intraspeciﬁc Behavioral Variation Mediates Insect Prey Survival via Direct and Indirect E ﬀ ects

: Conspeciﬁc individuals often exhibit behavioral di ﬀ erences that inﬂuence susceptibility to predation. Yet, how such trait di ﬀ erences scale to a ﬀ ect prey population regulation and community structure remains unclear. We used an 8 day ﬁeld mesocosm experiment to explore the e ﬀ ects of intraspeciﬁc prey behavioral trait variation on survival in an herbivorous insect community. We further manipulated spider predator composition to test for top-down context-dependence of behavioral e ﬀ ects. Insect prey behavioral trait variance inﬂuenced survival through both direct (i.e., variation among conspeciﬁcs) and indirect (i.e., variation among heterospeciﬁcs) mechanisms. The behavioral variance of two prey species, Philaenus and Orchelimum , directly reduced their survival, though for Philaenus , this direct negative e ﬀ ect only occurred in the presence of a single spider predator species. In contrast, the survival of Scudderia was enhanced by the behavioral trait variance of the surrounding insect community, an indirect positive e ﬀ ect. Taken together, these results emphasize the importance of accounting for intraspeciﬁc variation in community ecology, demonstrating novel pathways by which individual-level behavioral di ﬀ erences scale to alter population and community level patterns.


Introduction
The importance of prey traits for predatory interactions spans multiple biological scales. For example, predators exert selection on prey populations by preferentially consuming individuals with certain trait values [1][2][3][4][5][6], but also modify community structure by disproportionally consuming one prey species over another [7]. This latter effect is often ascribed to mean, species-level trait differences [8][9][10]. Yet these classic empirical approaches overlook the ecological effects of trait variance within prey populations, an outcome of trait differences among individuals that can vary independently from trait means [11,12]. While it is unlikely that a single prey individual's traits will modify outcomes at higher organizational scales (but see [13]), population-level prey trait variance could be an important property linking individual traits to the community-level effects of predators. Yet, few studies to date have experimentally explored the effects of prey trait variance on food webs [14][15][16]. Here, we used a field experiment to examine the effects of insect prey behavioral trait means and behavioral trait variances on the community-structuring effects of their predators (spiders).
Prey trait variance directly impacts prey survival. For example, size-selective predators often have an optimal prey body size for which predation rates are maximized (e.g., a hump-shaped attack rate distribution [17,18]), and the same likely applies to other types of prey traits [19]. Under such a scenario, increasing variance around the optimal trait mean would increase the prevalence of defended individuals, thereby reducing predation rate [11,20]. Alternatively, if the prey population trait mean deviates from the predator's optimal prey trait value, then increasing variance could increase predation by increasing the prevalence of vulnerable individuals, though this effect would depend on the dominant direction of trait variance [11,20]. These effects of prey trait variance might further depend on predator species identity because predator species often differ in their optimal prey trait values [8,21] as well as the type of selection that they induce (directional, stabilizing, or disruptive) [22].
Trait variance within prey populations could further influence predation rates on other prey species through indirect effects, though the nature of these effects may depend upon the type of phenotypic trait under consideration. If a species' trait variance directly reduces its vulnerability to predation, then this could turn predation pressure towards other prey species within the community. Alternatively, for traits that are relatively labile, such as behavior, the trait variance of interspecific prey could lead to trait changes in a focal species that render it more vulnerable to predation [23,24]. While this latter indirect effect would not occur with prey traits that are fixed over short durations (e.g., body size), this type of indirect effect could be prevalent for traits whose expression is context-dependent, such as behavior [25]. Nevertheless, the majority of studies examine the effects of prey behavioral traits at the individual level (i.e., a natural selection context) [3], while community-level effects of behavioral trait variance remain poorly understood.
Here, we explore how the behavioral trait variances and behavioral trait means of herbivorous insect populations mediate the effects of spider predators on insect survival. We took advantage of naturally occurring boldness (behavioral trait) variation within and across prey populations and manipulated predator species composition to test for context-dependence of behavioral effects. Our study reveals both direct and indirect effects of prey boldness variance, with some effects occurring only under certain spider predator species compositions. These findings demonstrate new pathways by which behavioral trait variation within populations, often overlooked within the field of community ecology [26][27][28], could scale to regulate community assembly and food web dynamics.

Study System and Animal Collection
We investigated the effects of prey behavior on prey survival using a spider predator-insect prey community system. Our experiment was conducted in a 60 × 12 m old-field dominated by ruderal plants (e.g., Solidago spp.) at the Donald S. Wood Field Laboratory (South Shenango, PA, USA). Spiders are important mesopredators within old-field communities, feeding on a variety of herbivorous insects that in turn regulate old-field primary production [29][30][31]. Preliminary sampling revealed that the jumping spiders Phidippus clarus and Platycryptus undatus were abundant predators at our field site, while four species dominated the herbivorous insect community: meadow froghoppers Philaenus spumarius, bush katydid nymphs from the genus Scudderia, common meadow katydids Orchelimum vulgare, and silver leafhoppers Athysanus argentarius. Given their abundance and thus high likelihood of interacting in the field, we selected these two predator and four prey species for our experiment. Prior to the experiment, we verified that both jumping spider species fed readily upon these insect prey species in one-on-one feeding trials. Spiders and insects were collected by sweep-netting and maintained individually in 50 mL clear plastic vials before the start of the experiment. Spiders were fed houseflies ad libitum during this 4 day housing period, while insect vials contained a single Solidago canadensis leaf for food and moisture. Individual spiders and insects were used once in the experiment, and survivors were released.

Experimental Design
We tested the effects of insect prey behavioral trait variances and behavioral trait means on prey survival within replicate field mesocosms. Each 0.36 m 2 area mesocosm was stocked with the same insect prey community: 10 Philaenus individuals, 10 Scudderia individuals, 5 Orchelimum individuals, and 5 Athysanus individuals, reflecting the relative abundance of these species at our old-field site. Prior to the placement of insect prey in mesocosms, we assayed the boldness behavior of each individual insect using the same assay (see Prey Behavioral Assay below). Insect prey were then added to mesocosms haphazardly and blindly with respect to their individual behavior ( Figure 1) in order to generate naturally occurring behavioral variation within and across prey species' populations. Thus, mesocosm prey communities varied continuously in prey species' mean boldness and boldness variance. We manipulated spider predator species composition across mesocosms to test whether potential effects of insect prey behavior on prey survival were context dependent (i.e., only occurred under certain predator species compositions). We manipulated spider species composition using three different treatments: two Phidippus spiders per mesocosm (n = 10 mesocosm replicates), two Platycryptus spiders per mesocosm (n = 13 mesocosm replicates), or one of each spider species per mesocosm (n = 12 mesocosm replicates). Given the 0.36 m 2 area of our mesocosms, these experimental predator densities were within the range of those observed within our old-field site, as well as those reported for Phidippus within other old-field sites (4-10 spiders per m 2 [29,32]).

Prey Behavioral Assay
Prior to starting the mesocosm experiment, we measured the behavior of each insect prey individual using the same boldness assay. Boldness, defined as an individual's propensity for risk-taking, can be measured in several ways [33,34]. Following other studies [35][36][37], we measured boldness as emergence time from a refuge, where individuals that emerge (i.e., enter a potentially dangerous situation) more quickly are considered bolder than individuals that take longer to emerge.
The boldness of each prey individual was measured by placing it, using a small capillary tube, at the bottom of a 50 mL tube (11 cm in length) that opened to a larger arena that insects were unfamiliar with. We then recorded the amount of time required for each insect to climb up and out of the tube after it had begun moving within the tube. Assays were terminated if the individual did not move within 300 s. Accordingly, all boldness measurements potentially surpassing 300 s were assigned the value of 300 s. While we did not assess the repeatability of insect boldness, previous work suggests that boldness is a repeatable trait in both Philaenus and Scudderia [38], as well as for several other herbivorous insect species [39].

Mesocosm Design and Assembly
We used 60 (L) × 60 (W) × 60 cm (H) collapsible cages (1450DSV: BioQuip Products, Inc.) as field mesocosms. These cages had a lumite mesh roof and three lumite mesh walls to allow for natural rainfall and humidity variation, while the fourth wall was covered with transparent plastic. Two weeks before mesocosms were set out, we transplanted Solidago canadensis stems to lawns that bordered opposite sides of the old-field. These plants were spaced 2 m apart, forming a line around the old-field. Immediately before adding the mesocosms, we cut the plants surrounding each transplanted plant to 6 cm in height. We then placed a mesocosm over each transplanted plant and lined the bottoms of mesocosms with topsoil to prevent organism escape or incursion. The goal of this procedure was to create identical patches of old-field within mesocosms. Solidago plants further provided food and refuge habitat for insect prey during the experiment. All mesocosms were oriented with their plastic side facing north to prevent overheating during sunrise and sunset periods. Before the addition of experimental animals to mesocosms, we conducted several visual scans of all mesocosms and removed naturally occurring insects and spiders.
After this assembly process, we added insect prey to mesocosms and then spiders. Spider predator species composition treatments were assigned randomly to mesocosms, and all treatments were run concurrently. Spiders hid within refuge habitat upon being added to mesocosms, and thus all predation events occurred under natural conditions over the course of the experiment. The boldness behavior of each insect prey individual was known, allowing us to calculate both the mean boldness and boldness coefficient of variation of each prey species' population within each mesocosm (see Analysis). While individual spider behavior was assayed prior to placement in mesocosms, spider behavior did not affect prey survival [40], and spider individuals were randomly assigned to mesocosms to mimic natural predator behavioral variation in the field. After the addition of spiders, mesocosm communities were left undisturbed for eight days, allowing natural trophic dynamics to occur. At the end of this eight-day period, mesocosms were deconstructed, and surviving insects were counted within each mesocosm to measure prey survival.

Analysis
Our analysis compared the effects of insect prey behavioral trait variances and behavioral trait means on prey survival and tested whether these effects depended on spider predator species composition. We built four generalized linear models, with one model for the survival of each insect prey species (Philaenus, Scudderia, Orchelimum, Athysanus). Due to overdispersion, we modeled proportional prey survival response variables using a quasibinomial error distribution within R [41]. Boldness trait variance was calculated as the coefficient of variation (CV; the standard deviation divided by the mean), a standard measure of dispersion. To separate the direct and indirect effects of prey behavior, each model included both the focal species' (i.e., the species whose survival was being modeled) mean boldness and boldness CV (direct effects), as well as the mean boldness and boldness CV of the rest of the prey community within each mesocosm, excluding the focal species (indirect effects). We further included two-way interactions between each of these four behavioral fixed effects (focal species mean boldness, focal species boldness CV, community mean boldness, community boldness CV) and spider predator species composition (three treatments: Phidippus spiders, Platycryptus spiders, or one individual from each spider species). This approach allowed us to assess (1) direct effects of prey behavior on survival (e.g., the effect of Philaenus behavior on Philaenus survival), (2) indirect effects of prey behavior on survival (e.g., the effects of the prey community, minus Philaenus, on Philaenus survival), and (3) the dependence of both direct and indirect effects on spider predator species composition. We tested the significance of fixed effects using likelihood-ratio Chi-square tests ("Anova" function in the "Car" package [42]).

Philaenus Survival
The population-level behavioral variance of Philaenus influenced its survival (i.e., a direct effect), but this effect was dependent on spider predator species composition. Philaenus survival declined with increasing boldness variance in the presence of Phidippus jumping spiders but not in the other predator treatments, driving a significant interaction between predator treatment and Philaenus boldness variation (Table 1, Figure 2). All other factors, including measurements of mean boldness, failed to predict Philaenus survival (Table 1).  Field mesocosms were stocked with 10 Philaenus individuals, and population-level boldness variance was calculated as the coefficient of variation in boldness across these 10 individuals. Shaded areas around marginal effect regression lines represent 84% confidence intervals generated using the "ggpredict" function in the R package "ggeffects" [43]. Eighty-four percent confidence intervals approximate significant differences at a 95% confidence level [44].

Scudderia Survival
Scudderia survival increased with increasing boldness variance of the prey community (i.e., an indirect effect), and this effect was consistent across spider predator treatments (Table 1, Figure 3). We also recorded a positive effect of mean community boldness on Scudderia survival, though this effect was relatively weak (Table 1). Scudderia population-level behavioral metrics did not influence its survival (i.e., no evidence for direct effects; Table 1). Figure 3. Indirect effect of prey community boldness variance on Scudderia proportional survival across all spider predator treatments. Field mesocosms were stocked with 10 Scudderia individuals, and community boldness variance was calculated as the coefficient of variation in boldness all prey species (Philaenus, Orchelimum, Athysanus) minus the focal species Scudderia. Shaded areas around marginal effect regression line represent 84% confidence intervals generated using the "ggpredict" function in the R package "ggeffects" [43].

Orchelimum Survival
Similar to Philaenus, we recorded a negative direct effect of Orchelimum population-level behavioral variance on its survival ( Figure 4, Table 1), though this effect was consistent across spider predator treatments (Table 1). All other factors, including measurements of mean boldness, failed to predict Orchelimum survival (Table 1). Figure 4. Direct effect of Orchelimum prey boldness variance on Orchelimum proportional survival across all spider predator treatments. Field mesocosms were stocked with 5 Orchelimum individuals, and population-level boldness variance was calculated as the coefficient of variation in boldness across these 5 individuals. Shaded areas around marginal effect regression line represent 84% confidence intervals generated using the "ggpredict" function in the R package "ggeffects" [43].

Athysanus survival
None of the model factors were significant predictors of Athysanus survival (Table 1).

Discussion
While myriad studies have established the importance of prey traits for predator-prey interactions [3,17,19], these studies are typically focused at either the individual level within a natural selection context [3] or at the community level within a species selection context [45]. Such traditional approaches overlook the ecological effects of population-or group-level trait variance [11], a potential key property linking the traits of individuals to communities. Here, using a field experiment, we showed that prey behavioral trait variances have decisive effects on herbivorous insect survival as mediated by spider predators, while behavioral trait means offer relatively less explanatory power. Importantly, such effects of prey behavioral trait variance on survival arose via both direct and indirect mechanisms. Behaviorally variant populations of both Philaenus and Orchelimum experienced reduced survival compared to behaviorally homogeneous populations (i.e., direct negative effects). Yet for Philaenus, this direct negative effect was only manifested in the presence of Phidippus spiders, a context-dependent effect of behavioral trait variance. We further detected an indirect effect whereby the aggregate behavioral trait variance of the surrounding prey community, Philaenus, Orchelimum, and Athysanus, improved the survival of Scudderia. In total, these findings revealed new and intricate pathways by which prey behavioral traits scale to shape prey community structure via a shared predator.

Differences across Prey Species
The nature of direct and indirect interactions within our study system was largely prey species dependent. While Philaenus and Orchelimum boldness trait variance reduced their probability of survival, we found no evidence for direct effects in the other two prey species. Furthermore, Scudderia was the only prey species whose survival was influenced by an indirect effect. This suggests that the effects of behavioral trait variance are not universal but rather dependent on prey species' functional traits.
Direct negative effects of Philaenus and Orchelimum boldness variance on survival could be due to spiders' increased capture efficiency of behaviorally extreme prey individuals (e.g., particularly bold or shy individuals) whose availability increases in more behaviorally variant groups. Previous work by Lichtenstein et al. (2017) confirmed this direct negative effect of Philaenus behavioral variance on Philaenus survival, though this study estimated behavioral hypervolumes from three Philaenus behavioral traits: activity level, boldness, and space use [38]. If these effects of behavioral variance persist in the field, then spiders could drive stabilizing selection on Philaenus and Orchelimum boldness via the improved survival of behaviorally intermediate individuals [46]. More specifically, Philaenus experienced this direct negative effect in the presence of Phidippus spiders but not Platycryptus spiders or the treatment in which both spider species were included. Such a context-dependent direct negative effect may stem from the unique predator avoidance (i.e., escape) abilities of Philaenus froghoppers, who are better jumpers for their size than fleas [47]. We suggest it is possible that Phidippus was the only predator able to capture Philaenus individuals at either end of the boldness spectrum.
Improved Scudderia survival in the presence of behaviorally variant insect prey communities could result from Scudderia behavioral modifications over the course of the experiment. While we assayed individual insect boldness a few days before running our experiment, this does not preclude behavioral changes during the experiment. As such, it is possible that the surrounding prey community modified Scudderia behavior in such a way that enhanced their survival. Clearly, additional work is required to probe the underlying behavioral mechanisms behind these interactions and assess their generality. We believe that a better understanding of the role of functional traits in regulating behaviorally mediated interactions will increase predictability across systems.

Indirect Effect Mediated by a Shared Predator
Our study sheds new light on interactions between prey species mediated by a shared predator, a major class of indirect effects. The most well-studied of these interactions is apparent competition, which in its most basic form describes the scenario in which prey species indirectly reduce one another's densities via the positive numerical response of a shared predator [48,49]. Other work shows that agonistic interactions between prey can similarly increase predation rates, even in the absence of any predator numerical response [24]. Our study describes a related indirect effect, apparent commensalism, where the behavioral variance of the surrounding prey community increases the survival of Scudderia in the presence of a shared predator.

Conclusion and Future Directions
Individual-level trait variation is ubiquitous within natural populations, but ecologists are just beginning to appreciate the full extent of its effects on community-level processes [26,27,50]. Our study demonstrates the importance of behavioral trait variance in regulating insect survival, effectively scaling the effects of individual traits to the community level. In contrast, we found the effects of behavioral trait means to be relatively weak in magnitude (i.e., having relatively less predictive power). Our findings further suggest prey species specificity in the effects of trait variance, adding a new level of complexity to our understanding of predator-prey systems. While such complexity could be expected to diminish the predictability of predator-prey dynamics, we contend that such idiosyncratic responses might help explain the striking variation in community structure observed along natural gradients. Our study suggests that the magnitude of direct and indirect effects could change substantially across ecological gradients (e.g., trait variance within species, the presence of specific predators), thereby supporting heterogeneous landscapes of species interactions [51]. Though underexplored, we expect that future empirical and theoretical explorations of such gradient-dependence will provide new insight into the structure and functioning of natural communities.