# Influence of intra- and interspecific variation in predator–prey body size ratios on trophic interaction strengths

Ross N Cuthbert, Ryan J Wasserman, Tatenda Dalu, Horst Kaiser, Olaf L F Weyl, Jaimie T A Dick, Arnaud Sentis, Michael W Mccoy, Mhairi E Alexander

<figure><img src="image_1.png" /></figure>

► To cite this version:

Ross N Cuthbert, Ryan J Wasserman, Tatenda Dalu, Horst Kaiser, Olaf L F Weyl, et al.. Influence of intra- and interspecific variation in predator–prey body size ratios on trophic interaction strengths. Ecology and Evolution, 2020, 10 (12), pp.5946-5962. 10.1002/ece3.6332. hal-03171607

HAL Id: hal-03171607

https://hal.inrae.fr/hal-03171607

Submitted on 17 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Influence of intra- and interspecific variation in predator-prey body size ratios on trophic interaction strengths

Ross N. Cuthbert<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID

Ryan J. Wasserman<fnref n="4" />,<fnref n="3" />

Tatenda Dalu<fnref n="5" />,<fnref n="3" /> ID

Horst Kaiser<fnref n="6" />

Olaf L. F. Weyl<fnref n="7" />,<fnref n="6" /> ID

Jaimie T. A. Dick<fnref n="2" />

Arnaud Sentis<fnref n="8" /> ID

Michael W. McCoy<fnref n="9" />

Mhairi E. Alexander<fnref n="10" />,<fnref n="11" />,<fnref n="3" />

## Correspondence

Ross N. Cuthbert, GEOMAR, Helmholtz-Zentrum für Ozeanforschung Kiel, 24105 Kiel, Germany.
Email: rossnoalcuthbert@gmail.com

## Funding information

Alexander von Humboldt Foundation; NRF-SAIAB Institutional Support system; National Research Foundation (NRF)-South

## Abstract

1. Predation is a pervasive force that structures food webs and directly influences ecosystem functioning. The relative body sizes of predators and prey may be an important determinant of interaction strengths. However, studies quantifying the combined influence of intra- and interspecific variation in predator-prey body size ratios are lacking.

2. We use a comparative functional response approach to examine interaction strengths between three size classes of invasive bluegill and largemouth bass toward three scaled size classes of their tilapia prey. We then quantify the influence of intra- and interspecific predator-prey body mass ratios on the scaling of attack rates and handling times.

3. Type II functional responses were displayed by both predators across all predator and prey size classes. Largemouth bass consumed more than bluegill at small and intermediate predator size classes, while large predators of both species were more similar. Small prey were most vulnerable overall; however, differential attack rates among prey were emergent across predator sizes. For both bluegill and largemouth bass, small predators exhibited higher attack rates toward small and intermediate prey sizes, while larger predators exhibited greater attack rates toward large prey. Conversely, handling times increased with prey size, with small bluegill exhibiting particularly low feeding rates toward medium-large prey types. Attack rates for both predators peaked unimodally at intermediate predator-prey body mass ratios, while handling times generally shortened across increasing body mass ratios.

4. We thus demonstrate effects of body size ratios on predator-prey interaction strengths between key fish species, with attack rates and handling times dependent on the relative sizes of predator-prey participants.

5. Considerations for intra- and interspecific body size ratio effects are critical for predicting the strengths of interactions within ecosystems and may drive differential ecological impacts among invasive species as size ratios shift.

<footnote n="1">¹GEOMAR, Helmholtz-Zentrum für Ozeanforschung Kiel, Kiel, Germany</footnote>

<footnote n="2">²Institute for Global Food Security, School of Biological Sciences, Queen's University Belfast, Belfast, UK</footnote>

<footnote n="3">³South African Institute for Aquatic Biodiversity, Makhanda, South Africa</footnote>

<footnote n="4">⁴Department of Zoology and Entomology, Rhodes University, Makhanda, South Africa</footnote>

<footnote n="5">⁵Department of Ecology and Resource Management, University of Venda, Thohoyandou, South Africa</footnote>

<footnote n="6">⁶Department of Ichthyology and Fisheries Science, Rhodes University, Makhanda, South Africa</footnote>

<footnote n="7">⁷DSI/NRF Research Chair in Inland Fisheries and Freshwater Ecology, South African Institute for Aquatic Biodiversity, Makhanda, South Africa</footnote>

<footnote n="8">⁸INRAE, Aix Marseille University, UMR RECOVER, Aix-en-Provence, France</footnote>

<footnote n="9">⁹Department of Biology, East Carolina University, Greenville, SC, USA</footnote>

<footnote n="10">¹⁰Institute for Biomedical and Environmental Health Research, School of Health and Life Sciences, University of the West of Scotland, Paisley, UK</footnote>

<footnote n="11">¹¹Department of Botany and Zoology, Centre for Invasion Biology, Stellenbosch University, Matieland, South Africa</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
© 2020 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.</footnote>
African Research Chairs Initiative of the
Department of Science and Innovation
(DSI), Grant/Award Number: UID 77444,
88746 and 110507; DSI-NRF Centre of
Excellence for Invasion Biology; NRF
Thuthuka, Grant/Award Number: 117700;
Natural Environment Research Council; DSI-
NRF CIB; U.S. National Science Foundation,
Grant/Award Number: 1556743

KEYWORDS

Bluegill, context-dependency, functional response, interaction strength, largemouth bass, piscivory, size-scaling, tilapia

# 1 | INTRODUCTION

Predation is an important factor in natural communities that exerts
a significant influence on both predators and prey (Dick et al., 2017;
Lima, 2002; Sih, Crowley, McPeek, Petranka, & Strohmeier, 1985).
From the perspective of the predator, the ability to capture prey
is critical in determining energy acquisition required for growth,
reproduction, and survival (Brose et al., 2019; Hempel, Neukamm,
& Thiel, 2016; Kaemingk, Graeb, & Willis, 2014; Steinhart, Stein,
& Marschall, 2004). For prey, it is an important driver in the dy-
namics that underpin many aspects of their behavior and mor-
phology allowing them to persist under predatory threats (Carlson
& Langkilde, 2014; Preisser, Bolnick, & Benard, 2005; Werner &
Peacor, 2003). As a result, a broad range of physiological, behav-
ioral, and morphological features have been shown to affect pred-
ator-prey interactions (Lagrue, Besson, & Lecerf, 2015; Portalier,
Fussmann, Loreau, & Cherif, 2019). Furthermore, predation out-
comes are often dependent on the combination of certain traits of
the species involved (Barrios-O'Neill et al., 2016; Barrios-O'Neill,
Kelly, & Emmerson, 2019; Christensen, 1996; Einfalt, Parkos, &
Wahl, 2015; Kolář, Boukal, & Sentis, 2019; Uma & Weiss, 2012).

Body size is an important trait that is strongly linked to forag-
ing outcomes. For instance, in fish, bigger predators often have
greater prey capture success owing to enhanced swimming abil-
ities, improved visual acuity, and changes to morphology such
as increased gape (Christensen, 1996; Graeb, Mangan, Jolley,
Wahl, & Dettmers, 2006; Miller, Crowder, & Rice, 1993). Similarly,
larger prey size permits greater reaction distances, faster escape
speeds, and often provides size refuges from predation (Cuthbert,
Callaghan, & Dick, 2019; Hansen & Wahl, 1981; Rodgers, Downing,
& Morrell, 2014). However, the contributions of these variables to
the interactions between predators and prey can be influenced
by the relative body size ratios of the participants, which can be
an important determinant in predator-prey dynamics (Asquith
& Vonesh, 2012; Barrios-O'Neill et al., 2016; Brose et al., 2019;
Jennings & Warr, 2003). For example, consumers may exhibit
lower attack rates toward relatively small and large resources, re-
sulting in a hump-shaped function that peaks at intermediate con-
sumer-resource ratios (Brose, 2010; McCoy, Bolker, Warkentin,
& Vonesh, 2011; Preston et al., 2018; but see Barrios-O'Neill
et al., 2016). As such, the relative differences in body size that
are observed in interacting predator and prey species can have
important implications for population structure and community
stability (Otto, Rall, & Brose, 2007; Sentis, Binzer, & Boukal, 2017).
Furthermore, the effects of body size ratio on trophic interactions

may depend on prey densities, owing to the underlying density
dependence of consumer-resource interactions (Holling, 1959).

The ability of an individual to capture resources or to avoid predation is not constant through time, but can vary as body size changes during ontogeny (e.g., McCoy et al., 2011). Many species exhibit size-specific ontogenetic shifts in resource use, which play an important role in determining growth and survival of individuals with a subsequent influence on population- and community-level processes (Reichstein, Persson, & De Roos, 2015; Rudolf & Rasmussen, 2013; Werner & Gilliam, 1984). Size-specific changes may emerge across age classes where older, larger individuals act as separate “ecological species” (Polis, 1984; Schröder, Nilsson, Persson, Kooten, & Reichstein, 2009). However, changes may also emerge in cohorts of a similar age that exhibit large variations in size (Brooks, McCoy, & Bolker, 2013; Pfister & Stevens, 2002). Likewise, this can result in discrete stages comprised of separate size classes with a range of functional roles that are based on selection of diet (De Roos, Leonardsson, Persson, & Mittelbach, 2002; De Roos, Persson, & McCauley, 2003; Huss, Persson, & Byström, 2007; Parkos & Wahl, 2010; Rudolf & Rasmussen, 2013; Urbatzka, Beeck, Van Der Velde, & Borcherding, 2008). Most species undergo such size-specific ontogenetic shifts in diet, often corresponding to external changes and internal factors, such as food supply or physiological demands (Werner & Gilliam, 1984). There is also much evidence to support the importance of relative body size to the outcome of these interactions between predatory fish and their prey (Fuiman, 1994; Lundvall, Svanbäck, Persson, & Byström, 1999; Scharf, Buckel, Juanes, & Conover, 1998). In many instances, prey size range can change during ontogeny, where maximum prey size increases while minimum prey size changes only slightly. This allows larger predators a greater prey size choice compared to smaller predators that are unable to feed on large prey and have thus a restricted prey size range (Scharf, Juanes, & Rountree, 2000). Energetic limitations, where a prey size does not provide sufficient energy compared with capture costs, may further influence prey selectivity as energetically costly prey may not be sustainable for the growth of predators (Portalier et al., 2019). However, there is a distinct lack of studies which examine effects of intraspecific size variations on the scaling of consumer-resource interactions.

Here, we assess the effect of variations in predator-prey body
size ratios on trophic interactions and quantify the influence of
both intra- and interspecific variations. We investigated whether
changes in the relative size of two fish predators (bluegill *Lepomis macrochirus* and largemouth bass *Micropterus salmoides*) and their
prey affect feeding interactions, and the relationship between prey
density and predation rate (i.e., the so-called "functional response"; Holling, 1959). Bluegill and largemouth bass are important predators in freshwaters, and both are global invaders that have been introduced extensively outside of their native range of North America. Both fishes are known to have major impacts on native communities due to predation (Almeida, Almodóvar, Nicola, Elvira, & Grossman, 2012; Ellender, Weyl, & Swartz, 2011; Godinho & Ferreira, 2000; Mittelbach, 1988). Bluegill and largemouth bass undergo ontogenic shifts in relation to resource and habitat use (Parkos & Wahl, 2010; Post, 2003; Werner & Hall, 1988), and are known to engage in piscivory (Azuma, 1992; Pelham, Pierce, & Larscheid, 2001; Taguchi, Miura, Krueger, & Sugiura, 2014). However, the window of opportunity for such changes may narrow over time because of the concurrent growth of the target fish prey (Olson, 1996).

Functional responses can be a useful means of quantifying and comparing effects of predators on prey populations under different environmental contexts (Dick et al., 2014). Both functional response form and magnitude are useful indicators of the strength of interactions between predators and their prey (Barrios-O'Neill et al., 2016; Dick et al., 2013; Dickey et al., 2020; Sheath et al., 2018). Three functional response forms are commonly characterized (Holling, 1966): Type I forms describe a linear relationship between consumption rates and resource density, and are typically associated with filter-feeding organisms that are not constrained by handling times (Jeschke, Kopp, & Tollrian, 2004); type II functional responses are hyperbolic, whereby most, if not all, resources are consumed at low densities and consumption rates are restricted at high densities by handling time; sigmoid type III functional responses allow for refuge at low densities as attack rate increases with prey density, before feeding rates again saturate at high densities. While the type II functional response is thought to be more destabilizing for prey populations and may lead to local extinctions, type III curves may stabilize resources (Juliano, 2001).

The aims of this study were to investigate intra- and interspecific resource use variations as a function of the relative predator body size of bluegill and largemouth bass toward a common fish prey species. Using functional responses, we quantified predator-prey interactions and asked whether predation rates were conserved across three size categories of both predators and prey, in a fully crossed experimental design. We predicted that larger fishes of both species would have higher maximum feeding rates of all prey sizes compared to small and intermediate predators. Similarly, we

expected predators of intermediate size to be more efficient predators on smaller and intermediate prey compared to large prey, and smaller predators to be most effective at consuming smaller prey. Further, we expected functional response attack rates to peak at intermediate predator-prey body size ratios, and for handling times to lengthen under increasing ratios. In addition, given the earlier switch to piscivory in largemouth bass and the reports of their voracity as predators (Alexander, Dick, Weyl, Robinson, & Richardson, 2014), we expected impacts on prey populations to be greater in this species.

## 2 | MATERIALS AND METHODS

### 2.1 | Animal maintenance

We examined functional responses of two fish species, the bluegill *L. macrochirus* and the largemouth bass *M. salmoides*, toward a prey fish, the Mozambique tilapia *Oreochromis mossambicus*. Young-of-year bluegill and largemouth bass were collected in March 2015 by seine netting in Mosslands (33°24′7.28″S; 26°27′6.89″E) and Yarrow (33°24′55.34″S; 26°22′34.98″E) reservoirs, near Makhanda (formerly Grahamstown), South Africa. Tilapia were supplied by AquaCulture Innovations, Makhanda. All fishes were transported to the Department of Ichthyology and Fisheries Science (DIFS), Rhodes University, Makhanda, and were housed in 600-L tanks in a closed recirculating system (water flow to the fish tanks 1 L/min; 18 ± 1°C). Fishes were allowed to acclimate to the system for at least 72 hr prior to use in predation trials.

### 2.2 | Functional response trials

Experiments were conducted on three size classes of each predator that were fully crossed with three size classes of prey. Bluegill and largemouth bass were assigned to one of three size classes based on their mass (see Table 1). Trials were completed in cages with volumes scaled to length of predatory fish (approx. 1 L 4 mm⁻¹; Table 1). Cages were scaled according to fish length in order to control for differences in predator-prey encounter rates that would naturally arise among size classes, therefore removing such differences as a confounding variable. Individual cages were constructed from 1.5 mm mesh and floated using buoyancy aids in 15 separate 300-L

<table><caption>TABLE 1 Mass (g ± SE) and length (mm ± SE) of largemouth bass, bluegill, and tilapia</caption><thead><tr><th></th><th>Size 1</th><th>Size 2</th><th>Size 3</th></tr></thead><tbody><tr><td>Largemouth bass (g, mm)</td><td>1.70 ± 0.19, 50.33 ± 0.78</td><td>4.95 ± 0.18, 73.40 ± 0.51</td><td>11.38 ± 0.47, 99.0 ± 1.62</td></tr><tr><td>Bluegill (g, mm)</td><td>1.65 ± 0.04, 45.59 ± 0.39</td><td>5.45 ± 0.11, 70.69 ± 0.61</td><td>12.29 ± 0.46, 92.33 ± 0.71</td></tr><tr><td>Tilapia (g, mm)</td><td>0.018 ± 0.001, 11.05 ± 0.05</td><td>0.035 ± 0.002, 14.65 ± 0.19</td><td>0.075 ± 0.004, 18.50 ± 0.14</td></tr><tr><td>Cage dimensions (cm, L)</td><td>24.5 × 23.0 × 23, 12.9</td><td>27.0 × 26.8 × 26.8, 19.4</td><td>29.5 × 29.5 × 29.5, 25.6</td></tr></tbody></table>

Note: Cage dimensions (height × breadth × width cm) and volume (L), scaled for length of predator, are also presented.
fiberglass tanks that were part of the same flow-through system as the holding tanks described above. This allowed five trials each of small, intermediate, and large predators to run at any one time. The order of trials was fully randomized across all treatments (predator species and size, prey size, prey density; see below). Due to the high number of trials (358 in total), fish predators were reused in the different size ratio treatments. Fish of each size were selected from a common pool of at least 60 individuals; thus, reuse occurred a maximum of two times. Fish were left for at least 2 days between trials and maintained on larval chironomids to minimize experimental prey learning. Fish predators were humanely euthanized at the end of the experiment with an overdose of clove oil.

Bluegill and largemouth bass of each size class (Table 1) were starved for 72 hr, before being randomly selected from their holding tanks and transferred to the appropriate size-scaled cage (see above) 1 hr prior to a trial. Pilot observations indicated that this was sufficient time for fish to settle in cages before consuming prey. Individual fish were then presented with one of three sizes of tilapia prey (Table 1) at five prey densities (2, 4, 8, 16, 32), with at least three replicates per density. Small tilapia prey were further presented at densities of 64 prey/cage to facilitate reaching an asymptote in consumption rates. Feeding trials were run for 1 hr, after which prey consumption was enumerated based on examination of remaining live prey. Controls were three replicates (one for each cage size) of each prey size and density in the absence of predators.

## 2.3 | Statistical analyses

Statistical analyses were conducted using R v3.4.4 (R Development Core Team, 2018). Differences in overall prey consumption among predator species, predator sizes, and prey sizes, and their two- and three-way interactions, were assessed using a generalized linear model (GLM) with a quasi-Poisson error distribution and log link, irrespective of prey density. This error family was selected owing to residual overdispersion. A backward stepwise deletion process was followed until all nonsignificant terms and interactions were removed, starting with the highest-order interaction term (Crawley, 2007). Effect sizes in the resulting model were inferred using F-tests via analysis of deviance (type III sums of squares). Prey densities of 64 were excluded as they were not present across all treatment groups (i.e., absent for medium and large prey). Significant effects in the model were analyzed with Tukey's contrast post hoc tests via least square means (Lenth, 2016). Statistical significance was inferred at the 95% confidence level in all analyses.

Binomial GLMs with logit links considering the proportion of prey consumed as a function of prey density were used to categorize functional responses. A significantly negative first-order term indicates a type II functional response, while a type III functional response is categorized by a significantly positive first-order term followed by a significantly negative second-order term (Juliano, 2001). Where evidence for a particular functional

response form was equivocal, Akaike's information criterion (AIC) was used to select models which minimized information loss, with lower values indicating a better fit (see Pritchard, Paterson, Bovy, & Barrios-O'Neill, 2017).

As prey were not replaced as they were consumed, Rogers' random predator equation was used to model functional responses (Rogers, 1972):

$$N_e = N_0 (1 - \exp(a(N_e h - T))) \quad (1)$$

where $N_e$ is the number of prey eaten, $N_0$ is the initial density of prey, $a$ is the attack rate (cage/hr), $h$ is the handling time (hr/prey), and $T$ is the experiment duration (i.e., 1 hr). The Lambert W function was used to fit the model, owing to the recursive nature of the random predator equation (Bolker, 2008; McCoy & Bolker, 2008). Multiple estimations of $a$ and $h$ parameters were generated via bootstrapping ($n = 20$ per experimental group; lower cap at 0). These parameters were next compared using Gamma GLMs (structured as before).

Bootstrapped $a$ and $h$ parameter estimates were then analyzed using polynomial regression models separately according to fish species as a function of predator-prey body mass ratios (continuous predictor; calculated from fish group mean masses). We $\log_{10}$-transformed $a$ estimates and body mass ratios prior to analyses, and $\log_{10}(x + 1)$-transformed $h$ estimates as $\log_{10}(0)$ is not defined and the bootstrap yielded a few null estimated values of handling time. Model comparisons were performed using AICc to select among linear, quadratic, and cubic polynomial models. Diagnostic plots were used to ensure residuals met the assumptions of parametric testing (Zuur, Ieno, & Elphick, 2010).

Owing to our nonreplacement experimental design, we tested whether prey depletion biased $a$ and $h$ parameter estimation from the random predator equation (Equation 1). To do so, we first simulated prey depletion during the time course of an experiment using the following population dynamic model considering a type II functional response:

$$\frac{dN}{dt} = - \frac{aNP}{1 + ahN} \quad (2)$$

where $N$ is the prey density, $P$ is the predator density, $a$ is the attack rate (cage/hr), and $h$ is the handling time (hr/prey). To generate predictions of expected prey survival in the experiment, initial values of $N_0$ and $P$ were set at the initial prey and predator densities corresponding to the experimental treatments, and the population dynamic model was integrated over the time interval of the experiment using the R package "DeSolve" (Soetaert, Petzoldt, & Setzer, 2010). To mimic our experimental data, we used the population dynamic model to simulate levels of prey depletion at four experimental durations (0.2, 0.5, 1.0, 1.5 hr) across the six experimental prey densities (2, 4, 8, 16, 32, 64). Input functional response parameters at each duration were systematically varied ($a$: 0.45, 1.00, 5.00; $h$: 0.002, 0.01, 0.30), while the opposing parameter was fixed at an intermediate value (i.e., $a = 1.00$; $h = 0.01$) ($n = 3$ simulated per density). Predator densities were maintained at one. Second, the random predator
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** Mean (±SE) prey consumed of different size classes by small, medium, and large bluegill and largemouth bass pooled across ubiquitous prey densities (up to 32 ind./cage)</figcaption></figure>

<table><caption>**TABLE 2** Linear coefficients (first-order terms) resulting from logistic regression of proportional prey consumption as a function of prey density across predator species and size classes</caption><thead><tr><th>Species</th><th>Prey size class</th><th>Predator size class</th><th>Linear coefficient, P</th></tr></thead><tbody><tr><td rowspan="9">Bluegill</td><td rowspan="3">Small</td><td>Small</td><td>-0.03, &lt;.001</td></tr><tr><td>Medium</td><td>-0.02, &lt;.001</td></tr><tr><td>Large</td><td>-0.02, &lt;.001</td></tr><tr><td rowspan="3">Medium</td><td>Small</td><td>-0.07, &lt;.001</td></tr><tr><td>Medium</td><td>-0.08, &lt;.001</td></tr><tr><td>Large</td><td>-0.11, &lt;.001</td></tr><tr><td rowspan="3">Large</td><td>Small</td><td>-0.05, 0.02</td></tr><tr><td>Medium</td><td>-0.04, 0.01</td></tr><tr><td>Large</td><td>-0.09, &lt;.001</td></tr><tr><td rowspan="6">Largemouth bass</td><td rowspan="3">Small</td><td>Small</td><td>-0.06, &lt;.001</td></tr><tr><td>Medium</td><td>-0.03, &lt;.001</td></tr><tr><td>Large</td><td>0.004, 0.35</td></tr><tr><td rowspan="3">Medium</td><td>Small</td><td>-0.15, &lt;.001</td></tr><tr><td>Medium</td><td>-0.08, &lt;.001</td></tr><tr><td>Large</td><td>-0.05, 0.003</td></tr><tr><td rowspan="3">Large</td><td>Small</td><td>-0.07, &lt;.001</td></tr><tr><td>Medium</td><td>-0.12, &lt;.001</td></tr><tr><td>Large</td><td>-0.11, &lt;.001</td></tr></tbody></table>

equation (Equation 1) was fit to these rounded consumption simulations at each duration and parameter input separately, to estimate *a* and *h*. These estimated values were then compared to the input *a* and *h* parameters based on the overlap of their standard error for each treatment, to test whether increasing experimental duration, and thus prey depletion, caused estimation bias. We calculated prey depletion for each duration and parameter scenario by dividing the total simulated number of prey eaten by the total of initial prey densities (i.e., *N*<sub>e</sub>/*N*<sub>0</sub>).

3 | RESULTS

Control prey without predators had 100% survival in all replicates, and thus, experimental deaths were attributed to consumption by predatory fish, which was also observed. Accordingly, it was not necessary to adjust functional responses for background prey mortality (see Rosenbaum & Rall, 2018). Smallest prey were most vulnerable overall ($F_{2,331} = 26.93, p < .001$; Figure 1) and experienced 158% and 73% greater mean mortality than large and medium prey, respectively (Tukey's test, both $p < .001$). Medium-sized prey were 49% more vulnerable than large prey (Tukey's test, $p = .04$). A significant interaction was found between predator species and predator size ($F_{2,331} = 7.28, p < .001$; Figure 1) due to largemouth bass consuming 235% and 122% more prey than bluegill when predators were small- and medium-sized, respectively (Tukey's test, both $p < .001$), yet only 27% more when both predators were large-sized (Tukey's test, $p = .08$). There were no other statistically clear two-way or three-way interactions, which were thus removed during the backward stepwise deletion process.

Linear coefficients considering prey consumption as a function of prey density were significantly negative, irrespective of predator species or predator and prey sizes, indicating type II functional responses (Table 2). Large largemouth bass consumption toward small prey was an exception to this owing to the lack of asymptotic prey densities (dotted curve, Figure 2b); however, the type II random predator model (Rogers, 1972) minimized information loss here as compared to the type III (Hassell 1978) and generalized models for conditions without prey replacement (both $\Delta$AIC > 2) (see Pritchard et al., 2017). Toward small prey, largemouth bass functional responses tended to be of higher magnitude than those of bluegill when predators were small- and medium-sized (Figure 2a,b). Toward medium-sized prey, largemouth bass functional response magnitudes were generally greater irrespective of their size (Figure 2c,d). Contrastingly, functional response magnitudes were more similar in height between predators toward large prey (Figure 2e,f).
FIGURE 2 Functional responses of three sizes classes of bluegill (a, c, e) and largemouth bass (b, d, f) toward small- (a, b), medium- (c, d), and large-sized (e, f) tilapia prey. Points are raw data. Note differences in axes scaling according to prey size

<figure><img src="image_4.png" /></figure>

FIGURE 3 Mean (±SE) bootstrapped attack rates toward tilapia prey of different size classes by small, medium, and large bluegill and largemouth bass

<figure><img src="image_5.png" /></figure>

Attack rates across sizes classes of bluegill and largemouth bass differed significantly according to prey size owing to a significant three-way interaction between predator species, predator size, and prey size ($F_{4,342} = 7.42, p < .001$) (Figure 3). Toward small prey, bluegill attack rates did not differ significantly among predator sizes (Tukey's test, all $p > .05$), while small and medium largemouth bass exhibited significantly higher attack rates than large largemouth bass (Tukey's test, both $p < .01$). Toward intermediate prey classes, small-sized bluegill attack rates were highest, followed by

large-sized bluegill, and both were significantly greater than medium-sized bluegill (Tukey's test, both $p < .05$). Similarly, small-sized largemouth bass attack rates were significantly greater than that observed in medium-sized largemouth bass (Tukey's test, $p < .01$). Toward large prey, large bluegill attack rates were significantly greater than the attack rates of either small- or medium-sized bluegill (Tukey's test, both $p < .001$). In turn, medium-sized bluegill attack rates were significantly higher than small-sized bluegill (Tukey's test, $p < .001$). Attack rates for largemouth bass were
<figure><img src="image_6.png" /></figure>

**FIGURE 4** Mean (±SE) bootstrapped handling times toward tilapia prey of different size classes by small, medium, and large bluegill and largemouth bass

<figure><img src="image_7.png" /></figure>

**FIGURE 5** Scaling of attack rate parameter with predator–prey body mass ratios between bluegill and largemouth bass. Polynomial regression lines are presented alongside raw data points

significantly higher by medium-sized predators compared to either small- or large-sized predators (Tukey's test, both *p* < .01). Accordingly, interactions between predator and prey sizes were displayed statistically, with small-sized predators generally exhibiting the highest attack rates toward small and intermediate prey classes, while attack rates of medium- and large-sized predators were greatest toward large prey.

A significant three-way interaction term also indicated differences in handling time among predator species and size classes toward differently sized prey (*F*_{4,342} = 3.44, *p* = .01) (Figure 4). For both predator species, handling times by small predators were always significantly longer than those of medium- and large-sized predators under all prey sizes (Tukey's test, all *p* < .001). Similarly, handling times of medium-sized bluegill were significantly longer

than those of large-sized bluegill toward small- and intermediate-sized prey (Tukey's test, both *p* < .001); this difference was not statistically significant toward large-sized prey (Tukey's test, *p* = .99). Contrasting, medium- and large-sized largemouth bass handling times were only statistically different toward small prey sizes (Tukey's test, *p* < .001). Overall, while handling times generally related negatively with predator size, we found species-specific differences, with small bluegill handling times being particularly longer than those of largemouth bass. Further, intermediate–large largemouth bass exhibited more similar prey handling capacities across prey types.

For both predator species, scaling of attack rates followed a unimodal pattern with predator–prey body mass ratios, peaking at intermediate values (Figure 5). According to their AICc values, the
best model was the cubic polynomial model for both species here (Table 3). Largemouth bass attack rates tended to reach a greater magnitude and peaked at lower predator-prey body bass ratios than bluegill.

Handling times generally decreased concurrently with increasing predator-prey body bass ratios for both fish species (Figure 6). The quadratic polynomial model was the best fit for bluegill, while the cubic model was selected in the case of largemouth bass (Table 3). Largemouth bass handling times tended to be reduced across

<table><caption>TABLE 3 Linear, quadratic, and cubic coefficients resulting from polynomial regression of bootstrapped functional response parameters as a function of predator-prey body mass ratios between predator species</caption><thead><tr><th>Species</th><th>Parameter</th><th>Term</th><th>Coefficient, P</th></tr></thead><tbody><tr><td rowspan="3">Bluegill</td><td rowspan="3">a</td><td>Linear</td><td>2.77, &lt;.001</td></tr><tr><td>Quadratic</td><td>-2.98, &lt;.001</td></tr><tr><td>Cubic</td><td>0.76, 0.11</td></tr><tr><td rowspan="3">Largemouth bass</td><td rowspan="3">a</td><td>Linear</td><td>-2.14, &lt;.001</td></tr><tr><td>Quadratic</td><td>-3.53, &lt;.001</td></tr><tr><td>Cubic</td><td>1.34, &lt;.001</td></tr><tr><td rowspan="3">Bluegill</td><td rowspan="3">h</td><td>Linear</td><td>-0.81, &lt;.001</td></tr><tr><td>Quadratic</td><td>0.15, 0.003</td></tr><tr><td>Cubic</td><td>—</td></tr><tr><td rowspan="3">Largemouth bass</td><td rowspan="3">h</td><td>Linear</td><td>-0.39, &lt;.001</td></tr><tr><td>Quadratic</td><td>0.11, &lt;.001</td></tr><tr><td>Cubic</td><td>-0.06, &lt;.001</td></tr></tbody></table>

Note: Inclusion of terms was based on ΔAkaike's Information Criterion.

predator-prey body mass ratios compared to bluegill, and particularly under low body mass ratio values.

When using the population dynamic model to simulate prey depletion across input parameter values and prey densities mimicking our experimental design, we found that estimated values of *a* and *h* were overall statistically similar to input “known” values of both *a* and *h* (duration: Figures A1 and A2; depletion: Figures A3 and A4). These included scenarios of total depletion at lowest prey densities, as per our experimental data. For instance, prey depletion was total (100%) for the scenario *a* = 5, *h* = 0.01 at prey densities 2, 4, 8, and 16 for experimental durations of 1 hr and more. Estimated attack rates and handling times always overlapped input values, with one exception at a duration shorter than our experiment (0.2 hr: *a* = 1.00, *h* = 0.30) (Figure A2e,f), and where prey depletion was low (Figure A4e,f). Accordingly, increasing experimental duration, and thus prey depletion, did not significantly bias estimates across a broad range of input *a* and *h* values.

## 4 | DISCUSSION

This study demonstrates the influence of intraspecific shifts in body size ratios during ontogeny on the strength of predator-prey interactions. Biotic processes such as predation exert profound influence upon the stability and functioning of ecosystems (Brose et al., 2019; Wasserman & Froneman, 2013), and interaction strengths between trophic groups have been shown to scale with relative body size (Barrios-O'Neill et al., 2016; Brose et al., 2019; Schröder, Kalninkat, & Arlinghaus, 2016). Given that body size ratios between predators and prey are highly variable, examinations of ecological impacts across body size variations are crucial to predict community

<figure><img src="image_8.png" /><figcaption>FIGURE 6 Scaling of handling time parameter with predator-prey body mass ratios between bluegill and largemouth bass. Polynomial regression lines are presented alongside raw data points</figcaption></figure>
outcomes (Asquith & Vonesh, 2012; Brose et al., 2019; Jennings & Warr, 2003). However, studies quantitatively examining the influence of body size variations on the strength of predator-prey interactions are lacking at the intraspecific level. As a result, there has been a relatively poor understanding of how consumer-resource participant heterogeneity can modulate predation effects (see Barrios-O'Neill et al., 2016; Gallagher, Brandl, & Stier, 2016; McCoy et al., 2011; Uiterwaal, Mares, & DeLong, 2017). This, in turn, limits insights into potential tradeoff dynamics between paired predators and prey given often asynchronous ontogenic development (e.g., Olson, 1996). The results here demonstrate that scalings of predator-prey body size can alter impacts from invasive fish toward their prey. Furthermore, they show that body size scalings can interact and manifest differently between predator species and potentially drive differential invader impacts.

We found that all size classes of both bluegill and largemouth bass exhibited a saturating type II functional response toward each size class of tilapia prey. However, there were marked differences in predation rates between the species, and predator-prey interactions differed within intraspecific body size ratios. Predation by largemouth bass was generally greater than in bluegill across all body size ratios examined here which was expected given the reported voracity of bass (e.g., Alexander et al., 2014). However, the higher predation of largemouth bass was most pronounced at small and intermediate predator sizes, and became less important where large individuals of both predator fish species were used. Indeed, there was no clear difference between feeding rates of medium- and large-sized largemouth bass, while raw consumption by bluegill always increased incrementally with body size. Although diets of many fish species change with ontogeny, for instance due to temporal variations in resource supplies or changeable physiological demands (Werner & Gilliam, 1984), the timings of shifts to piscivory can vary considerably between species. Largemouth bass can be piscivorous from as early as age-0 (Pelham et al., 2001), while shifts in bluegill piscivory may occur much later, owing to littoral zone occupancy during early ontogenic stages (Mittelbach, 1981). Thus, consumptive differences between the species presented here align with their documented feeding ecology, with interspecific differences greatest at small predator sizes. Accordingly, our results help quantitatively inform the nature of dietary shifts in empirical contexts and also help predict the ecological impacts of these invaders across size ratio differences. In particular, the high predatory impacts of largemouth bass corroborate well-established field impacts, whereby native prey populations are frequently reduced and other fish populations extirpated, while bluegill effects may be less evident (Ellender et al., 2011; Weyl, de Moor, Hill, & Weyl, 2010).

Overall, considering the effects of predator and prey body size separately, functional response magnitudes were generally highest toward smaller prey across all size classes of each predator species. Likewise, functional responses by larger predators were generally higher than in smaller body size classes indicating a higher maximum predation rate for larger predators compared to smaller

ones. While smaller prey sizes exhibit reduced reaction distances and escape speeds (Hansen & Wahl, 1981; Rodgers et al., 2014), larger predators frequently demonstrate enhanced capture rates due to improved visual acuity, increased gape, and greater swim speeds (Christensen, 1996; Graeb et al., 2006; Miller et al., 1993). Nevertheless, interactive complexity between relative predator-prey body sizes can emerge due to inefficiency in capturing prey that are relatively large or small (Brose, 2010; McCoy et al., 2011). On the other hand, sustained ecological impacts irrespective of predator size class have been shown in other studies (see Gallagher et al., 2016). Here, differences in the magnitude of functional responses of both predator species, regardless of predator size, were less marked toward large-sized prey. This may reflect general handling constraints associated with this prey size for the two predator species.

Maximum attack rates of predators are often observed at intermediate predator-prey body size ratios under certain conditions (Barrios-O'Neill et al., 2016). The present study demonstrates emergent effects of body size between predator-prey pairings on attack rates (i.e., capture/search efficiencies), with attack rates tending to scale unimodally. Small predators exhibited greater attack rates toward small- and medium-sized prey, while large bluegill and intermediate largemouth bass displayed the greatest capture efficiencies toward large prey. In turn, attack rates for both predators were lowest by large individuals toward small prey and by small predators toward large prey. This provides further evidence for the importance of predator-prey body size ratios for interaction strengths, with predators often exhibiting greater capture efficiencies toward similarly scaled prey in respect to size. Moreover, the heightened attack rate of medium-sized largemouth bass toward large prey may suggest that capture efficiencies instead peak under greater prey sizes than tested in the current study for large individuals of this species. Indeed, largemouth bass attack rates peaked unimodally under lower predator-prey body size ratios than bluegill in the present study. Given that high attack rates are conducive to greater ecological impact potential at low resource densities (Bollache, Dick, Farnsworth, & Montgomery, 2008; Cuthbert, Dickey, Coughlan, Joyce, & Dick, 2019; Dick et al., 2013), low prey density predatory impacts may be reduced toward relatively small or large prey. For large prey in particular, these patterns may elicit greater stability through refugia, given the greater range of capture efficiencies demonstrated toward this size class. Importantly, however, as with many comparative functional response studies (e.g., Sheath et al., 2018; Wasserman, Alexander, Dalu, et al., 2016), the present study employed a nonreplacement experimental design which resulted in consistently high prey depletion under certain treatments. While the statistical measures applied account for prey depletion, feeding interactions were likely constrained at low prey densities. Contrasting, replacement designs could allow for improved attack rate estimation (i.e., initial curve slope) and better discrimination between experimental treatment groups (Dick et al., 2014). Regardless, nonreplacement experimental designs such as ours provide useful comparative insights into context dependencies affecting trophic
interactions. Furthermore, our analyses through simulations indicated that increasing prey depletion through longer experimental duration did not cause bias in parameter estimation using the random predator equation, across a range of starting attack rate and handling time values. Accordingly, prey depletion likely did not impact upon the robustness of our empirical results, even where prey were totally depleted at some densities.

Patterning of handling times was more consistent in our study system than that of attack rates. We found that handling times unsurprisingly increased concurrently with prey size indicating that more time is needed to handle and digest larger prey. Handling times for both predators related negatively to increasing predator-prey body mass ratios. As prey size increased, medium- and large-sized bluegill generally exhibited greater similarity in handling times, while small bluegill displayed significant handling constraints relative to largemouth bass. Conversely, medium and large largemouth bass handling times were less consistently significantly different. These traits may again reflect a relatively early ontogenic shift to piscivory in largemouth bass (Pelham et al., 2001). As with attack rates, a greater range of handling times was displayed toward large-sized prey, likely driven by size-related handling constraints. Nonetheless, further research that presents different prey types simultaneously is required to elucidate selectivity or switching processes which may stabilize trophic interactions (e.g., Cuthbert, Dickey, McMorrow, Laverty, & Dick, 2018), given that prey types were provided singularly in the present study.

The application of functional responses offers useful insights into the density dependencies of interaction strengths comparatively (Jeschke, Kopp, & Tollrian, 2002; Wasserman, Alexander, Weyl, et al., 2016). While such comparisons can be difficult to parallel with empirical community dynamics (Vonesh, McCoy, Altwegg, Landi, & Measey, 2017a, 2017b), functional responses offer comparative insights into the influence of context dependencies on predatory impacts (Sentis & Boukal, 2018; Sheath et al., 2018; Wasserman, Alexander, Dalu, et al., 2016). Our findings demonstrate effects of body size ratios on the predatory interaction strengths between piscivorous fishes and their prey. Moreover, we show clear consumptive differences between two key invaders which may relate to differential impacts on prey across their ontogeny. While traits other than body size can influence interaction strengths (e.g., temperature: Cuthbert, Dick, Callaghan, & Dickey, 2018; Englund, Öhlund, Hein, & Diehl, 2011; Sentis, Hemptinne, & Brodeur, 2013), effects of body size have been shown to be particularly pervasive across trophic and taxonomic groups (Barrios-O'Neill et al., 2016; Brose et al., 2019). The present study furthers these findings, highlighting the importance of both inter- and intraspecific differences in predator-prey body size ratios for trophic interactions. Therefore, asynchronous ontogenic development between predators and prey likely affect community outcomes with respect to population stabilities. We thus propose that considerations for the intra- and interspecific scaling of body size in consumer-resource systems can greatly enhance our predictive capacity for community interactions, invasive species ecological impacts, and ecosystem structuring.

## ACKNOWLEDGMENTS

RNC is funded through a research fellowship from the Alexander von Humboldt Foundation. We acknowledge use of infrastructure and equipment provided by the NRF-SAIAB Research Platforms and the funding channeled through the NRF-SAIAB Institutional Support system. This study was partially funded by the National Research Foundation (NRF) (UID 77444, 88746) —South African Research Chairs Initiative of the Department of Science and Innovation (DSI) (Inland Fisheries and Freshwater Ecology, Grant No. 110507) and the DSI-NRF Centre of Excellence for Invasion Biology (CIB). TD is supported by an NRF Thuthuka Grant (Grant No. 117700). JTAD acknowledges funding from the Natural Environment Research Council. MWM acknowledges support from a fellowship from the DSI-NRF CIB, and the U.S. National Science Foundation (Grant No. 1556743). The Department of Economic Development, Environmental Affairs and Tourism (Cacadu Region) is thanked for issuing research permits (11/15CR & 12/15CR).

## CONFLICT OF INTEREST

The authors declare no conflicts of interest.

## AUTHOR CONTRIBUTIONS

**Ross N. Cuthbert:** Formal analysis (lead); validation (lead); visualization (lead); writing – original draft (lead); writing – review & editing (lead). **Ryan J. Wasserman:** Conceptualization (equal); data curation (equal); investigation (equal); methodology (equal); writing – review & editing (equal). **Tatenda Dalu:** Investigation (equal); methodology (equal); writing – review & editing (equal). **Horst Kaiser:** Methodology (equal); resources (equal); writing – review & editing (equal). **Olaf L. F. Weyl:** Conceptualization (equal); funding acquisition (equal); methodology (equal); project administration (equal); resources (equal); supervision (equal); writing – review & editing (equal). **Jaimie T. A. Dick:** Conceptualization (equal); project administration (equal); supervision (equal); writing – review & editing (equal). **Arnaud Sentis:** Formal analysis (equal); methodology (equal); software (equal); visualization (equal); writing – review & editing (equal). **Michael W. McCoy:** Conceptualization (equal); methodology (equal); writing – review & editing (equal). **Mhairi E. Alexander:** Conceptualization (equal); data curation (equal); funding acquisition (equal); investigation (equal); methodology (equal); writing – review & editing (equal).

## DATA AVAILABILITY STATEMENT

All data from this manuscript are freely available at the Dryad Digital Repository: https://doi.org/10.5061/dryad.7m0cfxppt (Cuthbert et al., 2020).

## ORCID

Ross N. Cuthbert

<figure><img src="image_9.png" /></figure>

https://orcid.org/0000-0003-2770-254X

Tatenda Dalu

<figure><img src="image_10.png" /></figure>

https://orcid.org/0000-0002-9019-7702

Olaf L. F. Weyl

<figure><img src="image_11.png" /></figure>

https://orcid.org/0000-0002-8935-3296

Arnaud Sentis

<figure><img src="image_12.png" /></figure>

https://orcid.org/0000-0003-4617-3620
REFERENCES

Alexander, M. E., Dick, J. T. A., Weyl, O. L. F., Robinson, T. B., & Richardson, D. M. (2014). Existing and emerging high impact invasive species are characterized by higher functional responses than natives. *Biology Letters, **10***, 20130946. https://doi.org/10.1098/rsbl.2013.0946

Almeida, D., Almodóvar, A., Nicola, G. G., Elvira, B., & Grossman, G. D. (2012). Trophic plasticity of invasive juvenile largemouth bass *Micropterus salmoides* in Iberian streams. *Fisheries Research, **113***, 153-158. https://doi.org/10.1016/j.fishres.2011.11.002

Asquith, C. M., & Vonesh, J. R. (2012). Effects of size and size structure on predation and inter-cohort competition in red-eyed treefrog tadpoles. *Oecologia, **170***, 629-639. https://doi.org/10.1007/s0044 2-012-2332-x

Azuma, M. (1992). Ecological release in feeding behaviour: The case of bluegills in Japan. *Hydrobiologia, **243***, 269-276. https://doi.org/10.1007/BF00007042

Barrios-O'Neill, D., Kelly, R., Dick, J. T. A., Ricciardi, A., Macisaac, H. J., & Emmerson, M. C. (2016). On the context-dependent scaling of consumer feeding rates. *Ecology Letters, **19***, 668-678. https://doi.org/10.1111/ele.12605

Barrios-O'Neill, D., Kelly, R., & Emmerson, M. E. (2019). Biomass encounter rates limit the size scaling of feeding interactions. *Ecology Letters, **22***)11, 1870-1878. https://doi.org/10.1111/ele.13380

Bolker, B. (2008). *emdbook: Ecological models and data in R*. Princeton, NJ: Princeton University Press.

Bollache, L., Dick, J. T. A., Farnsworth, K. D., & Montgomery, W. I. (2008). Comparison of the functional responses of invasive and native amphipods. *Biology Letters, **4***, 166-169. https://doi.org/10.1098/rsbl.2007.0554

Brooks, M. E., McCoy, M. W., & Bolker, B. M. (2013). A method for detecting positive growth autocorrelation without marking individuals. *PLoS ONE, **8***, e76389. https://doi.org/10.1371/journal.pone.0076389

Brose, U. (2010). Body-mass constraints on foraging behaviour determine population and food-web dynamics. *Functional Ecology, **24***, 28-34. https://doi.org/10.1111/j.1365-2435.2009.01618.x

Brose, U., Archambault, P., Barnes, A. D., Bersier, L.-F., Boy, T., Canning-Clode, J., ... Iles, A. C. (2019). Predator traits determine food-web architecture across ecosystems. *Nature Ecology and Evolution, **3***, 919-927. https://doi.org/10.1038/s41559-019-0899-x

Carlson, B. E., & Langkilde, T. (2014). Predation risk in tadpole populations shapes behavioural responses of prey but not strength of trait-mediated indirect interactions. *Oikos, **123***, 1519-1527. https://doi.org/10.1111/oik.01488

Christensen, B. (1996). Predator foraging capabilities and prey antipredator behaviours: Pre-versus postcapture constraints on size-dependent predator-prey interactions. *Oikos, **76***, 368-380. https://doi.org/10.2307/3546209

Crawley, M. J. (2007). *The R book*. Chichester, UK: Wiley.

Cuthbert, R. N., Callaghan, A., & Dick, J. T. A. (2019). Differential interaction strengths and prey preferences across larval mosquito ontogeny by a cohabiting predatory midge. *Journal of Medical Entomology, **56***, 1428-1432. https://doi.org/10.1093/jme/tjz059

Cuthbert, R. N., Dick, J. T. A., Callaghan, A., & Dickey, J. W. E. (2018). Biological control agent selection under environmental change using functional responses, abundances and fecundities; the Relative Control Potential (RCP) metric. *Biological Control, **121***, 50-57. https://doi.org/10.1016/j.biocontrol.2018.02.008

Cuthbert, R. N., Dickey, J. W. E., Coughlan, N. E., Joyce, P. W. S., & Dick, J. T. A. (2019). The Functional Response Ratio (FRR): Advancing comparative metrics for predicting the ecological impacts of invasive alien species. *Biological Invasions, **21***, 2543-2547. https://doi.org/10.1007/s10530-019-02002-z

Cuthbert, R. N., Dickey, J. W. E., McMorrow, C., Laverty, C., & Dick, J. T. A. (2018). Resistance is futile: Lack of predator switching and a preference for native prey predict the success of an invasive prey species.

Royal Society Open Science, 5, 180339. https://doi.org/10.1098/rsos.180339

Cuthbert, R. N., Wasserman, R. J., Dalu, T., Kaiser, H., Weyl, O. L. F., Dick, J. T. A., ... Alexander, M. E. (2020). Influence of intra- and interspecific variation in predator-prey body size ratios on trophic interaction strengths. *Dryad Digital Repository*. https://doi.org/10.5061/dryad.7m0cxfppt

De Roos, A. M., Leonardsson, K., Persson, L., & Mittelbach, G. G. (2002). Ontogenetic niche shifts and flexible behavior in size-structured populations. *Ecological Monographs, **72***, 271-292. https://doi.org/10.1890/0012-9615(2002)072[0271:ONSAFB]2.0.CO;2

De Roos, A. M., Persson, L., & McCauley, E. (2003). The influence of size-dependent life-history traits on the structure and dynamics of populations and communities. *Ecology Letters, **6***, 473-487. https://doi.org/10.1046/j.1461-0248.2003.00458.x

Dick, J. T. A., Alexander, M. E., Jeschke, J. M., Ricciardi, A., MacIsaac, H. J., Robinson, T. B., ... Richardson, D. M. (2014). Advancing impact prediction and hypothesis testing in invasion ecology using a comparative functional response approach. *Biological Invasions, **16***, 735-753. https://doi.org/10.1007/s10530-013-0550-8

Dick, J. T. A., Gallagher, K., Avlijas, S., Clarke, H. C., Lewis, S. E., Leung, S., ... Ricciardi, A. (2013). Ecological impacts of an invasive predator explained and predicted by comparative functional responses. *Biological Invasions, **15***), 837-846. https://doi.org/10.1007/s10530-012-0332-8

Dick, J. T. A., Laverty, C., Lennon, J. J., Barrios-O'Neill, D., Mensink, P. J., Robert Britton, J., ... Caffrey, J. M. (2017). Invader Relative Impact Potential: A new metric to understand and predict the ecological impacts of existing, emerging and future invasive alien species. *Journal of Applied Ecology, **54***), 1259-1267. https://doi.org/10.1111/1365-2664.12849

Dickey, J. W. E., Cuthbert, R. N., South, J., Britton, J. R., Caffrey, J., Chang, X., ... Dick, J. T. A. (2020). On the RIP: Using Relative Impact Potential to assess the ecological impacts of invasive alien species. *NeoBiota, **55***), 27-60. https://doi.org/10.3897/neobiota.55.49547

Einfalt, L. M., Parkos, J. J., & Wahl, D. H. (2015). Effect of predator size and prey characteristics on piscivory of juvenile largemouth bass. *Transactions of the American Fisheries Society, **144***), 682-692. https://doi.org/10.1080/00028487.2015.1027002

Ellender, B. R., Weyl, O. L. F., & Swartz, E. R. (2011). Invasion of a headwater stream by non-native fishes in the Swartkops River system, South Africa. *African Zoology, **46***), 39-46. https://doi.org/10.1080/15627020.2011.11407477

Englund, G., Öhlund, G., Hein, C. L., & Diehl, S. (2011). Temperature dependence of the functional response. *Ecology Letters, **14***), 914-921. https://doi.org/10.1111/j.1461-0248.2011.01661.x

Fuiman, L. A. (1994). The interplay of ontogeny and scaling in the interactions of fish larvae and their predators. *Journal of Fish Biology, **45***), 55-79. https://doi.org/10.1111/j.1095-8649.1994.tb01084.x

Gallagher, A. J., Brandl, S. J., & Stier, A. C. (2016). Intraspecific variation in body size does not alter the effects of mesopredators on prey. *Royal Society Open Science, **3***), 160414. https://doi.org/10.1098/rsos.160414

Godinho, F. N., & Ferreira, M. T. (2000). Composition of endemic fish assemblages in relation to exotic species and river regulation in a temperate stream. *Biological Invasions, **2***), 231-244.

Graeb, B. D. S., Mangan, M. T., Jolley, J. C., Wahl, D. H., & Dettmers, J. M. (2006). Ontogenetic changes in prey preference and foraging ability of yellow perch: Insights based on relative energetic return of prey. *Transactions of the American Fisheries Society, **135***), 1493-1498.

Hansen, M. J., & Wahl, D. H. (1981). Selection of small *Daphnia pulex* by yellow perch fry in Oneida Lake, New York. *Transactions of the American Fisheries Society, **110***), 64-71. https://doi.org/10.1577/1548-8659(1981)110<64:SO SDPB>2.0.CO;2

Hassell, M. P. (1978). *The dynamics of arthropod predator-prey systems*. Princeton, NJ: Princeton University Press.
Hempel, M., Neukamm, R., & Thiel, R. (2016). Effects of introduced round goby (*Neogobius melanostomus*) on diet composition and growth of zander (*Sander lucioperca*), a main predator in European brackish waters. *Aquatic Invasions*, 11, 167–178. https://doi.org/10.3391/ai.2016.11.2.06

Holling, C. S. (1959). Some characteristics of simple types of predation and parasitism. *Canadian Entomologist*, 91, 385-398. https://doi.org/10.4039/Ent91385-7

Holling, C. S. (1966). The functional response of invertebrate predators to prey density. *Memoirs of the Entomological Society of Canada*, 48, 1–86. https://doi.org/10.4039/entm9848fv

Huss, M., Persson, L., & Byström, P. (2007). The origin and development of individual size variation in early pelagic stages of fish. *Oecologia*, 153, 57–67. https://doi.org/10.1007/s00442-007-0719-x

Jennings, S., & Warr, K. (2003). Smaller predator-prey body size ratios in longer food chains. *Proceedings of the Royal Society of London. Series B: Biological Sciences*, 270, 1413–1417. https://doi.org/10.1098/rspb.2003.2392

Jeschke, J. M., Kopp, M., & Tollrian, R. (2002). Predator functional responses: Discriminating between handling and digesting prey. *Ecological Monographs*, 72, 95–112. https://doi.org/10.1890/0012-9615(2002)072[0095:PFRDBH]2.0.CO;2

Jeschke, J. M., Kopp, M., & Tollrian, R. (2004). Consumer-food systems: Why type I functional responses are exclusive to filter feeders. *Biological Reviews*, 79, 337–349. https://doi.org/10.1017/S1464793103006286

Juliano, S. A. (2001). Nonlinear curve fitting: Predation and functional response curves. In S. M. Scheiner, & J. Gurevitch (Eds.), *Design and analysis of ecological experiments* (pp. 159–182). New York, NY: Oxford University Press.

Kaemingk, M. A., Graeb, B. D. S., & Willis, D. W. (2014). Temperature, hatch date, and prey availability influence age-0 yellow perch growth and survival. *Transactions of the American Fisheries Society*, 143, 845–855. https://doi.org/10.1080/00028487.2014.886622

Kolář, V., Boukal, D., & Sentis, A. (2019). Predation risk and habitat complexity modify intermediate predator feeding rates and energetic efficiencies in a tri-trophic system. *Freshwater Biology*, 64, 1480–1491. https://doi.org/10.1111/fwb.13320

Lagrue, C., Besson, A. A., & Lecerf, A. (2015). Interspecific differences in antipredator strategies determine the strength of non-consumptive predator effects on stream detritivores. *Oikos*, 124, 1589–1596. https://doi.org/10.1111/oik.02272

Lenth, R. V. (2016). Least-square means: The R package lsmeans. *Journal of Statistical Software*, 69, 1–33.

Lima, S. L. (2002). Putting predators back into behavioral predator-prey interactions. *Trends in Ecology and Evolution*, 17, 70–75. https://doi.org/10.1016/S0169-5347(01)02393-X

Lundvall, D., Svanbäck, R., Persson, L., & Byström, P. (1999). Size-dependent predation in piscivores: Interactions between predator foraging and prey avoidance abilities. *Canadian Journal of Fisheries and Aquatic Sciences*, 56, 1285–1292.

McCoy, M. W., & Bolker, B. M. (2008). Trait-mediated interactions: Influence of prey size, density and experience. *Journal of Animal Ecology*, 77, 478–486. https://doi.org/10.1111/j.1365-2656.2008.01372.x

McCoy, M. W., Bolker, B. M., Warkentin, K. M., & Vonesh, J. R. (2011). Predicting predation through prey ontogeny using size-dependent functional response models. *American Naturalist*, 177, 752–766. https://doi.org/10.1086/659950

Miller, T. J., Crowder, L. B., & Rice, J. A. (1993). Ontogenetic changes in behavioural and histological measures of visual acuity in three species of fish. *Environmental Biology of Fishes*, 37, 1–8. https://doi.org/10.1007/BF00000707

Mittelbach, G. G. (1981). Foraging efficiency and body size: A study of optimal diet and habitat use. *Ecology*, 62, 1370–1386.

Mittelbach, G. G. (1988). Competition among refuging sunfishes and effects of fish density on littoral zone invertebrates. *Ecology*, 69, 614–623. https://doi.org/10.2307/1941010

Olson, M. H. (1996). Ontogenetic niche shifts in largemouth bass: Variability and consequences for first-year growth. *Ecology*, 77, 179–180. https://doi.org/10.2307/2265667

Otto, S. B., Rall, B. C., & Brose, U. (2007). Allometric degree distributions facilitate food-web stability. *Nature*, 450, 1226-1229. https://doi.org/10.1038/nature06359

Parkos, J. J., & Wahl, D. H. (2010). Intra- and intersystem variation in largemouth bass recruitment: Reproduction, prey availability, and the timing of year-class establishment. *Transactions of the American Fisheries Society*, 139, 1372–1385. https://doi.org/10.1577/T09-097.1

Pelham, M. E., Pierce, C. L., & Larscheid, J. G. (2001). Diet dynamics of the juvenile piscivorous fish community in Spirit Lake, Iowa, USA, 1997-1998. *Ecology of Freshwater Fish*, 10, 198-211. https://doi.org/10.1034/j.1600-0633.2001.100402.x

Pfister, C. A., & Stevens, F. R. (2002). The genesis of size variability in plants and animals. *Ecology*, 83, 59–72. https://doi.org/10.1890/0012-9658(2002)083[0059:TGOSVI]2.0.CO;2

Polis, G. A. (1984). Age structure component of niche width and intra-specific resource partitioning: Can age groups function as ecological species? *American Naturalist*, 123, 541–564. https://doi.org/10.1086/284221

Portalier, S. M. J., Fussmann, G. F., Loreau, M., & Cherif, M. (2019). The mechanics of predator-prey interactions: First principles of physics predict predator-prey size ratios. *Functional Ecology*, 33, 323–334. https://doi.org/10.1111/1365-2435.13254

Post, D. M. (2003). Individual variation in the timing of ontogenetic niche shifts in largemouth bass. *Ecology*, 84, 1298–1310. https://doi.org/10.1890/0012-9658(2003)084[1298:IVITTO]2.0.CO;2

Preisser, E. L., Bolnick, D. I., & Benard, M. F. (2005). Scared to death? The effects of intimidation and consumption in predator-prey interactions. *Ecology*, 86, 501–509. https://doi.org/10.1890/04-0719

Preston, D. L., Henderson, J. S., Falke, L. P., Segui, L. M., Layden, T. J., & Novak, M. (2018). What drives interaction strengths in complex food webs? A test with feeding rates of a generalist stream predator. *Ecology*, 99, 1591–1601. https://doi.org/10.1002/ecy.2387

Pritchard, D. W., Paterson, R. A., Bovy, H. C., & Barrios-O'Neill, D. (2017). Frair: An R package for fitting and comparing consumer functional responses. *Methods in Ecology and Evolution*, 8, 1528–1534.

R Development Core Team (2018). *A language and environment for statistical computing*. Vienna, Austria: R Foundation for Statistical Computing.

Reichstein, B., Persson, L., & De Roos, A. M. (2015). Ontogenetic asymmetry modulates population biomass production and response to harvest. *Nature Communications*, 6, 6441. https://doi.org/10.1038/ncomms7441

Rodgers, G. M., Downing, B., & Morrell, L. J. (2014). Prey body size mediates the predation risk associated with being "odd". *Behavioural Ecology*, 26, 242–246.

Rogers, D. (1972). Random search and insect population models. *Journal of Animal Ecology*, 41, 369–383. https://doi.org/10.2307/3474

Rosenbaum, B., & Rall, B. (2018). Fitting functional responses: Direct parameter estimation by simulating differential equations. *Methods in Ecology and Evolution*, 9, 2076–2090. https://doi.org/10.1111/2041-210X.13039

Rudolf, V. H. W., & Rasmussen, N. L. (2013). Ontogenetic functional diversity: Size structure of a keystone predator drives functioning of a complex ecosystem. *Ecology*, 94, 1046–1056. https://doi.org/10.1890/12-037.1

Scharf, F. S., Buckel, J. A., Juanes, F., & Conover, D. O. (1998). Predation by juvenile piscivorous bluefish (*Pomatomus saltatrix*): The influence of prey to predator size ratio and prey type on predator capture
success and prey profitability. *Canadian Journal of Fisheries and Aquatic Sciences*, 55, 1695-1703.

Scharf, F. S., Juanes, F., & Rountree, R. A. (2000). Predator size-prey size relationships of marine fish predators: Interspecific variation and effects of ontogeny and body size on trophic-niche breadth. *Marine Ecology Progress Series*, 208, 229-248.

Schröder, A., Kalninkat, G., & Arlinghaus, R. (2016). Individual variation in functional response parameters is explained by body size but not by behavioural types in a poeciliid fish. *Oecologia*, 182, 1129-1140.
https://doi.org/10.1007/s00442-016-3701-7

Schröder, A., Nilsson, K.A., Persson, L., Kooten, T.V., & Reichstein, V. (2009). Invasion success depends on invader body size in a size-structured mixed predation-competition community. *Journal of Animal Ecology*, 78, 1151-1162. https://doi.org/10.1111/j.1365-2656.2009.01590.x

Sentis, A., Binzer, A., & Boukal, D. S. (2017). Temperature-size responses alter food chain persistence across environmental gradients. *Ecology Letters*, 20, 852-862. https://doi.org/10.1111/ele.12779

Sentis, A., & Boukal, D. S. (2018). On the use of functional responses to quantify emergent multiple predator effects. *Scientific Reports*, 8, 11787. https://doi.org/10.1038/s41598-018-30244-9

Sentis, A., Hemptinne, J.-L., & Brodeur, J. (2013). Effects of simulated heat waves on an experimental plant-herbivore-predator food chain. *Global Change Biology*, 19, 833-842. https://doi.org/10.1111/gcb.12094

Sheath, D., Dick, J. T. A., Dickey, J. W. E., Guo, Z., Andreou, D., & Britton, J. R. (2018). Winning the arms race: Host-parasite shared evolutionary history reduces infection risks in fish final hosts. *Biology Letters*, 14, 20180363. https://doi.org/10.1098/rsbl.2018.0363

Sih, A., Crowley, P. H., McPeek, M., Petranka, J., & Strohmeier, K. (1985). Predation, competition, and prey communities - A review of field experiments. *Annual Review of Ecology and Systematics*, 16, 269-311.
https://doi.org/10.1146/annurev.es.16.110185.001413

Soetaert, K., Petzoldt, T., & Setzer, R. W. (2010). Solving differential equations in R: Package deSolve. *Journal of Statistical Software*, 33, 1-25.

Steinhart, G. B., Stein, R. A., & Marschall, E. A. (2004). High growth rate of young-of-the-year smallmouth bass in Lake Erie: A result of the round goby invasion? *Journal of Great Lakes Research*, 30, 381-389.
https://doi.org/10.1016/S0380-1330(04)70355-X

Taguchi, T., Miura, Y., Krueger, D., & Sugiura, S. (2014). Utilizing stomach content and faecal DNA analysis techniques to assess the feeding behaviour of largemouth bass *Micropterus salmoides* and bluegill *Lepomis macrochirus*. *Journal of Fish Biology*, 84, 1271-1288.

Uiterwaal, S., Mares, C., & DeLong, J. P. (2017). Body size, body size ratio, and prey type influence the functional response of damselfly nymphs. *Oecologia*, 185, 339-346. https://doi.org/10.1007/s0044 2-017-3963-8

Uma, D. B., & Weiss, M. R. (2012). Flee or fight: Ontogenetic changes in the behavior of cobweb spiders in encounters with spider-hunting wasps. *Environmental Entomology*, 41, 1474-1480. https://doi.org/10.1603/EN12126

Urbatzka, R., Beeck, P., Van Der Velde, G., & Borcherding, J. (2008). Alternative use of food resources causes intra-cohort variation in the size distribution of young-of-the-year perch (*Perca fluviatilis*). *Ecology of Freshwater Fish*, 17, 475-480.

Vonesh, J., McCoy, M. W., Altwegg, R., Landi, P., & Measey, J. (2017a). Functional responses can't unify invasion ecology. *Biological Invasions*, 19, 1673-1676. https://doi.org/10.1007/s10530-016-1356-2

Vonesh, J., McCoy, M. W., Altwegg, R., Landi, P., & Measey, J. (2017b). Rather than unifying invasion biology Dick et al's approach rests on subjective foundations. *Biological Invasions*, 19(5), 1679-1680.
https://doi.org/10.1007/s10530-016-1361-5

Wasserman, R. J., Alexander, M. E., Dalu, T., Ellender, B., Kaiser, H., & Weyl, O. L. F. (2016). Using functional responses to quantify interaction effects among predators. *Functional Ecology*, 30, 1988-1998.
https://doi.org/10.1111/1365-2435.12682

Wasserman, R. J., Alexander, M. E., Weyl, O. L. F., Barrios-O'Neill, D., Froneman, P. W., & Dalu, T. (2016). Emergent effects of structural complexity and temperature on predator-prey interactions. *Ecosphere*, 7, e01239. https://doi.org/10.1002/ecs2.1239

Wasserman, R. J., Noyon, M., Avery, T. S., & Froneman, P. W. (2013). Trophic level stability-inducing effects of predaceous early juvenile fish in an estuarine mesocosm study. *PLoS ONE*, 8, e61019. https://doi.org/10.1371/journal.pone.0061019

Werner, E. E., & Gilliam, J. F. (1984). The ontogenetic niche and species interactions in size-structured populations. *Annual Review of Ecology and Systematics*, 15, 393-425. https://doi.org/10.1146/annurev.es.15.110184.002141

Werner, E. E., & Hall, D. J. (1988). Ontogenetic habitat shifts in bluegill: The foraging rate-predation risk trade-off. *Ecology*, 69, 1352-1366.

Werner, E. E., & Peacor, S. D. (2003). A review of trait-mediated indirect interactions in ecological communities. *Ecology*, 84, 1083-1100. https://doi.org/10.1890/0012-9658(2003)084[1083:AROTI I]2.0.CO;2

Weyl, P. S. R., de Moor, F. C., Hill, M. P., & Weyl, O. L. F. (2010). The effect of largemouth bass *Micropterus salmoides* on aquatic macro-invertebrate communities in the Wit River, Eastern Cape, South Africa. *African Journal of Aquatic Science*, 35, 273-281.

Zuur, A. F., Ieno, E. N., & Elphick, C. S. (2010). A protocol for data exploration to avoid common statistical problems. *Methods in Ecology and Evolution*, 1, 3-14. https://doi.org/10.1111/j.2041-210X.2009.00001.x

**How to cite this article:** Cuthbert RN, Wasserman RJ, Dalu T, et al. Influence of intra- and interspecific variation in predator-prey body size ratios on trophic interaction strengths. *Ecol Evol.* 2020;10:5946-5962. https://doi.org/10.1002/ece3.6332
APPENDIX 1

<figure><img src="image_13.png" /><figcaption>**FIGURE A1** Functional response parameter estimates (attack rates: a, c, e; handling times, b, d, f) modelled from consumption simulations across different experimental durations and input attack rate values. Input handling times were fixed at 0.01. Estimates are ±SE and the solid lines represent the input parameter value. Note differences in axes scaling</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>FIGURE A2 Functional response parameter estimates (attack rates: a, c, e; handling times, b, d, f) modelled from consumption simulations across different experimental durations and input handling time values. Input attack rates were fixed at 1.00. Estimates are ±SE and the solid lines represent the input parameter value. Note differences in axes scaling.</figcaption></figure>
<figure><img src="image_15.png" /></figure>

**FIGURE A3** Functional response parameter estimates (attack rates: a, c, e; handling times, b, d, f) modelled from consumption simulations across different prey depletions (total prey eaten divided by total density) and input attack rate values. Input handling times were fixed at 0.01. Estimates are ±SE and the solid lines represent the input parameter value. Note differences in axes scaling.
<figure><img src="image_16.png" /><figcaption>FIGURE A4 Functional response parameter estimates (attack rates: a, c, e; handling times, b, d, f) modelled from consumption simulations across different prey depletions (total prey eaten divided by total density) and input handling time values. Input attack rates were fixed at 1.00. Estimates are ±SE and the solid lines represent the input parameter value. Note differences in axes scaling.</figcaption></figure>