Feeding on multiple sources: towards a universal parameterization of the functional response of a generalist predator allowing for switching

Functional Responses Shape Node and Network Level Properties of a Simplified Boreal Food Web

Ecological communities are fundamentally connected through a network of trophic interactions that are often complex and difficult to model. Substantial variation exists in the nature and magnitude of these interactions across various predators and prey and through time. However, the empirical data needed to characterize these relationships are difficult to obtain in natural systems, even for relatively simple food webs. Consequently, prey-dependent relationships and specifically the hyperbolic form (Holling’s Type II), in which prey consumption increases with prey density but ultimately becomes saturated or limited by the time spent handling prey, are most widely used albeit often without knowledge of their appropriateness. Here, we investigate the sensitivity of a simplified food web model for a natural, boreal system in the Kluane region of the Yukon, Canada to the type of functional response used. Intensive study of this community has permitted best-fit functional response relationships to be determined, which comprise linear (type I), hyperbolic (type II), sigmoidal (type III), prey- and ratio-dependent relationships, and inverse relationships where kill rates of alternate prey are driven by densities of the focal prey. We compare node- and network-level properties for a food web where interaction strengths are estimated using best-fit functional responses to one where interaction strengths are estimated exclusively using prey-dependent hyperbolic functional responses. We show that hyperbolic functional responses alone fail to capture important ecological interactions such as prey switching, surplus killing and caching, and predator interference, that in turn affect estimates of cumulative kill rates, vulnerability of prey, generality of predators, and connectance. Exclusive use of hyperbolic functional responses also affected trends observed in these metrics over time and underestimated annual variation in several metrics, which is important given that interaction strengths are typically estimated over relatively short time periods. Our findings highlight the need for more comprehensive research aimed at characterizing functional response relationships when modeling predator-prey interactions and food web structure and function, as we work toward a mechanistic understanding linking food web structure and community dynamics in natural systems.

Evolution of prudent predation in complex food webs

Prudent predators catch sufficient prey to sustain their populations but not as much as to undermine their populations' survival. The idea that predators evolve to be prudent has been dismissed in the 1970s, but the arguments invoked then are untenable in the light of modern evolution theory. The evolution of prudent predation has repeatedly been demonstrated in two-species predator-prey metacommunity models. However, the vigorous population fluctuations that these models predict are not widely observed. Here we show that in complex model food webs prudent predation evolves as a result of consumer-mediated ('apparent') competitive exclusion of resources, which disadvantages aggressive consumers and does not generate such fluctuations. We make testable predictions for empirical signatures of this mechanism and its outcomes. Then we discuss how these predictions are borne out across freshwater, marine and terrestrial ecosystems. Demonstrating explanatory power of evolved prudent predation well beyond the question of predator-prey coexistence, the predicted signatures explain unexpected declines of invasive alien species, the shape of stock-recruitment relations of fish, and the clearance rates of pelagic consumers across the latitudinal gradient and 15 orders of magnitude in body mass. Specific research to further test this theory is proposed.

Assassin snails (Anentome helena) as a biological model for exploring the effects of individual specialisation within generalist predators

Quantifying feeding behaviour of generalist predators at the population and individual levels is crucial for understanding the structure and functioning of food webs. Individual predator/consumer feeding niches can be significantly narrower than that of the population across animal taxa. In such species, the population of a generalist predator becomes essentially an ensemble of specialist individuals and this often highly affects the dynamics of the prey-predator interactions. Currently, few experimental systems exist that are both easily technically manipulated in a lab and are reliable to accurately assess effects of individual specialisation within generalist predators. Here we argue that a freshwater predaceous snail, Anentome helena (also known as an ‘assassin snail’), is a convenient and reliable experimental system to study feeding of a generalist predator on multiple food types which exhibits well-pronounced specialisation of foraging individuals. Using A. helena we experimentally test: (i) how relative prey abundances in the environment affect the feeding patterns, (ii) whether the feeding patterns are consistent over the duration of the experimental period, and (iii) compare the feeding niche breadth of individuals to that of the laboratory population. By offering four different prey snail species, at a range of relative abundances, we show that there are consistent patterns in feeding. Importantly, the consumption of each prey was independent of the relative abundance at which they were present. Individual predators showed selectivity to a particular prey, i.e. the population of assassin snails seems to be formed of individuals that specialise on different prey. Our findings would contribute to the recent revision and the ongoing debate on the classification of predator species into generalists and specialists.

Integrating disparate datasets to model the functional response of a marine predator: A case study of harbour porpoises in the southern North Sea

Quantifying consumption and prey choice for marine predator species is key to understanding their interaction with prey species, fisheries, and the ecosystem as a whole. However, parameterizing a functional response for large predators can be challenging because of the difficulty in obtaining the required data on predator diet and on the availability of multiple prey species.This study modeled a multi-species functional response (MSFR) to describe the relationship between consumption by harbour porpoises (Phocoena phocoena) and the availability of multiple prey species in the southern North Sea. Bayesian methodology was employed to estimate MSFR parameters and to incorporate uncertainties in diet and prey availability estimates. Prey consumption was estimated from stomach content data from stranded harbour porpoises. Prey availability to harbour porpoises was estimated based on the spatial overlap between prey distributions, estimated from fish survey data, and porpoise foraging range in the days prior to stranding predicted from telemetry data.Results indicated a preference for sandeels in the study area. Prey switching behavior (change in preference dependent on prey abundance) was confirmed by the favored type III functional response model. Variation in the size of the foraging range (estimated area where harbour porpoises could have foraged prior to stranding) did not alter the overall pattern of the results or conclusions.Integrating datasets on prey consumption from strandings, predator foraging distribution using telemetry, and prey availability from fish surveys into the modeling approach provides a methodological framework that may be appropriate for fitting MSFRs for other predators.

Phenotypic plasticity as a cause and consequence of population dynamics

Predicting complex species-environment interactions is crucial for guiding conservation and mitigation strategies in a dynamically changing world. Phenotypic plasticity is a mechanism of trait variation that determines how individuals and populations adapt to changing and novel environments. For individuals, the effects of phenotypic plasticity can be quantified by measuring environment-trait relationships, but it is often difficult to predict how phenotypic plasticity affects populations. The assumption that environment-trait relationships validated for individuals indicate how populations respond to environmental change is commonly made without sufficient justification. Here we derive a novel general mathematical framework linking trait variation due to phenotypic plasticity to population dynamics. Applying the framework to the classical example of Nicholson's blowflies, we show how seemingly sensible predictions made from environment-trait relationships do not generalise to population responses. As a consequence, trait-based analyses that do not incorporate population feedbacks risk mischaracterising the effect of environmental change on populations.

Building clone-consistent ecosystem models

Many ecological studies employ general models that can feature an arbitrary number of populations. A critical requirement imposed on such models is clone consistency: If the individuals from two populations are indistinguishable, joining these populations into one shall not affect the outcome of the model. Otherwise a model produces different outcomes for the same scenario. Using functional analysis, we comprehensively characterize all clone-consistent models: We prove that they are necessarily composed from basic building blocks, namely linear combinations of parameters and abundances. These strong constraints enable a straightforward validation of model consistency. Although clone consistency can always be achieved with sufficient assumptions, we argue that it is important to explicitly name and consider the assumptions made: They may not be justified or limit the applicability of models and the generality of the results obtained with them. Moreover, our insights facilitate building new clone-consistent models, which we illustrate for a data-driven model of microbial communities. Finally, our insights point to new relevant forms of general models for theoretical ecology. Our framework thus provides a systematic way of comprehending ecological models, which can guide a wide range of studies.

Hidden layers of density dependence in consumer feeding rates

Functional responses relate a consumer's feeding rates to variation in its abiotic and biotic environment, providing insight into consumer behaviour and fitness, and underpinning population and food-web dynamics. Despite their broad relevance and long-standing history, we show here that the types of density dependence found in classic resource- and consumer-dependent functional-response models equate to strong and often untenable assumptions about the independence of processes underlying feeding rates. We first demonstrate mathematically how to quantify non-independence between feeding and consumer interference and between feeding on multiple resources. We then analyse two large collections of functional-response data sets to show that non-independence is pervasive and borne out in previously hidden forms of density dependence. Our results provide a new lens through which to view variation in consumer feeding rates and disentangle the biological underpinnings of species interactions in multi-species contexts.

Systematic bias in studies of consumer functional responses

Functional responses are a cornerstone to our understanding of consumer-resource interactions, so how to best describe them using models has been actively debated. Here we focus on the consumer dependence of functional responses to evidence systematic bias in the statistical comparison of functional-response models and the estimation of their parameters. Both forms of bias are universal to nonlinear models (irrespective of consumer dependence) and are rooted in a lack of sufficient replication. Using a large compilation of published datasets, we show that - due to the prevalence of low sample size studies - neither the overall frequency by which alternative models achieve top rank nor the frequency distribution of parameter point estimates should be treated as providing insight into the general form or central tendency of consumer interference. We call for renewed clarity in the varied purposes that motivate the study of functional responses, purposes that can compete with each other in dictating the design, analysis and interpretation of functional-response experiments.

From winter to summer and back: Lessons from the parameterization of a seasonal food web model for the Białowieża forest

Dynamic food web models describe how species abundances change over time as a function of trophic and life-history parameters. They are essential to predicting the response of ecosystems to perturbations. However, they are notoriously difficult to parameterize, so that most models rely heavily either on allometric scaling of parameters or inverse estimation of biomass flows. The allometric approach makes species of comparable body mass have near-identical parameters which can generate extinctions within a trophic level. The biomass flow approach is more precise, but is restricted to steady-states, which is not appropriate for time-varying environments. Adequately parameterizing large food webs of temperate and arctic environments requires dealing both with many species of similar sizes and a strongly seasonal environment. Inspired by the rich empirical knowledge on the vertebrate food web of the Białowieża forest, we parameterize a bipartite food web model comprising 21 predators and 124 prey species. Our model is a non-autonomous coupled ordinary differential equations system that allows for seasonality in life-history and predation parameters. Birth and death rates, seasonal descriptions of diet for each species, food requirements and biomass information are combined into a seasonal parameterization of a dynamic food web model. Food web seasonality is implemented with time-varying intrinsic growth rate and interaction parameters, while predation is modelled with both type I and type II functional responses. All our model variants allow for >80% persistence in spite of massive apparent competition, and a quantitative match to observed (seasonal) biomasses. We also identify trade-offs between maximizing persistence, reproducing observed biomasses, and ensuring model robustness to sampling errors. Although multi-annual cycles are expected with type II functional responses, they are here prevented by a strong predator self-regulation. We discuss these results and possible improvements on the model. We provide a general workflow to parameterize dynamic food web models in seasonal environments, based on a real case study. This may help to better predict how biodiverse food webs respond to changing environments.

Generalist predator dynamics under kolmogorov versus non-Kolmogorov models

Ecosystems often contain multiple species across two or more trophic levels, with a variety of interactions possible. In this paper we study two classes of models for generalist predators that utilize more than one food source. These models fall into two categories: predator - two prey and predator - prey - subsidy models. For the former, we consider a generalist predator which utilizes two distinct prey species, modelled via a Kolmogorov system of equations with Type II response functions. For the latter, we consider a generalist predator which exploits both a prey population and an allochthonous resource which is provided as a subsidy to the system exogenously, again with Type II response functions. This latter class of model is no longer Kolmogorov in form, due to an exogenous forcing term modelling the input of the allochthonous resource into the system. We non-dimensionalize both models, so that their respective parameter spaces may be more easily compared, and study the dynamics possible from each type of model, which will then indicate - for specific parameter regimes - which generalist predator's preferences are more favorable to survival, including the prevalence of coexistence states. We also consider the various non-equilibrium dynamics emergent from such models, and show that the non-Kolmogorov predator - prey - subsidy model of 10 admits more regular dynamics (including steady states and one type of limit cycle), whereas the predator - two prey Kolmogorov model can feature multiple types of limit cycles, as well as multistability resulting in strong sensitivity to initial conditions (with stable limit cycles and steady states both coexisting for the same model parameters). Our results highlight several interesting differences and similarities between Kolmogorov and non-Kolmogorov models for generalist predators.

An updated perspective on spiders as generalist predators in biological control

The role of generalist predators in biological control remains controversial as they may not only reduce pest populations but also disrupt biocontrol exerted by other natural enemies. Here, we focus on spiders as a model group of generalist predators. They are among the most abundant and most diverse natural enemies in agroecosystems. We review their functional traits that influence food-web dynamics and pest suppression at organisational levels ranging from individuals to communities. At the individual and population levels, we focus on hunting strategy, body size, life stage, nutritional target, and personality (i.e., consistent inter-individual differences in behaviour). These functional traits determine the spider trophic niches. We also focus on the functional and numerical response to pest densities and on non-consumptive effects of spiders on pests. At the community level, we review multiple-predator effects and effect of alternative prey on pest suppression. Evidence for a key role of spiders in pest suppression is accumulating. Importantly, recent research has highlighted widespread non-consumptive effects and complex intraguild interactions of spiders. A better understanding of these effects is needed to optimize biocontrol services by spiders in agroecosystems.

Insecticides alter prey choice of potential biocontrol agent Philodromus cespitum (Araneae, Philodromidae)

Even though pesticides can have various sublethal effects on behaviour of biocontrol agents, no study to date has investigated the effects of pesticides on the prey choice of generalist predators. Prey choice of generalist predators is among key factors determining the predation pressure they exert on pests, because it influences fitness of predators and consequently their densities and per capita capture rate. Here, we investigated the effect of Integro (a.i. methoxyfenozide) and SpinTor 480 Sc (a.i. spinosad) on prey choice and predatory activity of the spider Philodromus cespitum, which is known significantly to reduce hemipteran and dipteran pests in fruit orchards. We compared the prey preferences of philodromids between the psyllid pest Cacopsylla pyri (Hemiptera, Psyllidae) and beneficial Theridion sp. spiders in laboratory experiments. We found that both pesticides altered the prey preferences of philodromids. While the philodromids in a control treatment preferred theridiid spiders over the psyllid pest, philodromids in the pesticide treatments exhibited no significant prey preferences. The changes in prey preferences were caused by increased predation on the psyllids, while the predation on the theridiid remained similar. We suggest that the changes in prey preferences might theoretically be due to 1) impaired sensory systems, 2) altered taste, and/or 3) altered mobility. In combination with other studies finding reduced predation on fruit flies in P. cespitum after exposure to SpinTor, our results indicate that the sublethal effects of pesticides on predatory behaviour of generalist predators can depend on prey type and/or prey community composition.

Bistability induced by generalist natural enemies can reverse pest invasions

Analytical modeling of predator-prey systems has shown that specialist natural enemies can slow, stop and even reverse pest invasions, assuming that the prey population displays a strong Allee effect in its growth. We aimed to formalize the conditions in which spatial biological control can be achieved by generalists, through an analytical approach based on reaction-diffusion equations. Using comparison principles, we obtain sufficient conditions for control and for invasion, based on scalar bistable partial differential equations. The ability of generalist predators to control prey populations with logistic growth lies in the bistable dynamics of the coupled system, rather than in the bistability of prey-only dynamics as observed for specialist predators attacking prey populations displaying Allee effects. As a consequence, prey control is predicted to be possible when space is considered in additional situations other than those identified without considering space. The reverse situation is also possible. None of these considerations apply to spatial predator-prey systems with specialist natural enemies.

Delay driven spatiotemporal chaos in single species population dynamics models

Questions surrounding the prevalence of complex population dynamics form one of the central themes in ecology. Limit cycles and spatiotemporal chaos are examples that have been widely recognised theoretically, although their importance and applicability to natural populations remains debatable. The ecological processes underlying such dynamics are thought to be numerous, though there seems to be consent as to delayed density dependence being one of the main driving forces. Indeed, time delay is a common feature of many ecological systems and can significantly influence population dynamics. In general, time delays may arise from inter- and intra-specific trophic interactions or population structure, however in the context of single species populations they are linked to more intrinsic biological phenomena such as gestation or resource regeneration. In this paper, we consider theoretically the spatiotemporal dynamics of a single species population using two different mathematical formulations. Firstly, we revisit the diffusive logistic equation in which the per capita growth is a function of some specified delayed argument. We then modify the model by incorporating a spatial convolution which results in a biologically more viable integro-differential model. Using the combination of analytical and numerical techniques, we investigate the effect of time delay on pattern formation. In particular, we show that for sufficiently large values of time delay the system's dynamics are indicative to spatiotemporal chaos. The chaotic dynamics arising in the wake of a travelling population front can be preceded by either a plateau corresponding to dynamical stabilisation of the unstable equilibrium or by periodic oscillations.

A hybrid behavioural rule of adaptation and drift explains the emergent architecture of antagonistic networks

Ecological processes that can realistically account for network architectures are central to our understanding of how species assemble and function in ecosystems. Consumer species are constantly selecting and adjusting which resource species are to be exploited in an antagonistic network. Here we incorporate a hybrid behavioural rule of adaptive interaction switching and random drift into a bipartite network model. Predictions are insensitive to the model parameters and the initial network structures, and agree extremely well with the observed levels of modularity, nestedness and node-degree distributions for 61 real networks. Evolutionary and community assemblage histories only indirectly affect network structure by defining the size and complexity of ecological networks, whereas adaptive interaction switching and random drift carve out the details of network architecture at the faster ecological time scale. The hybrid behavioural rule of both adaptation and drift could well be the key processes for structure emergence in real ecological networks.

Competence of hosts and complex foraging behavior are two cornerstones in the dynamics of trophically transmitted parasites

Multi-host trophically transmitted parasite (TTP) is a common life cycle where prey and predators are respectively intermediate and definitive hosts of the parasite. In these systems, the foraging response of the predator toward variations in prey community composition underlies the dynamic of the parasite. Therefore, modeling epidemiological dynamic of infectious diseases considering ecological predator-prey interactions is essential to understand the spreading of parasites in ecosystems. However, two important weaknesses of previous TTP models including feeding interaction can be pointed out: (i) the choice of a linear density-dependent contact rate is faintly realistic as it supposes an unlimited ingestion rate with an increase of prey density and (ii) considering only one host prey species prevents the study of host biodiversity effect due to change in the prey community composition where species have different competences to be infected and to transmit the parasite. This article attempts to address the dynamics of parasite in a context of multiple intermediate hosts differentiated by their competences and of complex foraging behavior of the predator. We present and analyze a deterministic one predator-two prey model, which is then used to explore the transmission cycle of the cestode Echinococcus multilocularis. This study examines the foraging condition for the co-existence of the prey, and then, based on the computation of the threshold measure of disease risk, R0, we show that the pattern of feeding interactions changes the relationship between disease risk and prey community composition. Finally, we disentangle the mechanism leading to the counter-intuitive observation of a decrease of disease risk while the population density of intermediate hosts increases.

Long-term transients and complex dynamics of a stage-structured population with time delay and the Allee effect

Traditionally, mathematical modeling in population ecology is mainly focused on asymptotic behavior of the model, i.e. as given by the system attractors. Recently, however, transient regimes and especially long-term transients have been recognized as playing a crucial role in the dynamics of ecosystems. In particular, long-term transients are a potential explanation of ecological regime shifts, when an apparently healthy population suddenly collapses and goes extinct. In this paper, we show that the interplay between delay in maturation and a strong Allee effect can result in long-term transients in a single species system. We first derive a simple 'conceptual' model of the population dynamics that incorporates both a strong Allee effect and maturation delay. Unlike much of the previous work, our approach is not empirical since our model is derived from basic principles. We show that the model exhibits a high complexity in its asymptotic dynamics including multi-periodic and chaotic attractors. We then show the existence of long-term transient dynamics in the system, when the population size oscillates for a long time between locally stable stationary states before it eventually settles either at the persistence equilibrium or goes extinct. The parametric space of the model is found to have a complex structure with the basins of attraction corresponding to the persistence and extinction states being of a complicated shape. This impedes the prediction of the eventual fate of the population, as a small variation in the maturation delay or the initial population size can either bring the population to extinction or ensure its persistence.

Imperfect prey selectivity of predators promotes biodiversity and irregularity in food webs

Ecological communities are often characterised by many species occupying the same trophic level and competing over a small number of vital resources. The mechanisms maintaining high biodiversity in such systems are still poorly understood. Here, we revisit the role of prey selectivity by generalist predators in promoting biodiversity. We consider a generic tri-trophic food web, consisting of a single limiting resource, a large number of primary producers and a generalist predator. We suggest a framework to describe the predator functional response, combining food selectivity for distinctly different functional prey groups with proportion-based consumption of similar prey species. Our simulations reveal that intermediate levels of prey selectivity can explain a high species richness, functional biodiversity, and variability among prey species. In contrast, perfect food selectivity or purely proportion-based food consumption leads to a collapse of prey functional biodiversity. Our results are in agreement with empirical phytoplankton rank-abundance curves in lakes.

Prey Carrying Capacity Modulates the Effect of Predation on Prey Diversity

Understanding the role of predation in regulating prey diversity is a major goal in ecology, with profound consequences for community dynamics, ecosystem structure, and conservation practice. Deterministic differential equation models predict that some predation regimes, such as prey-switching predation, should promote prey coexistence and increase prey diversity. However, such models do not capture stochastic population fluctuations that are ubiquitous in empirical study sites and nature reserves. In this article, we examine the effects of prey-switching predation on the species richness of prey communities with demographic noise. We show that in finite, discrete prey populations, the ability of prey-switching predation to promote diversity depends on the carrying capacity of the prey community and the richness of the source pool for prey. Identical predation regimes may have opposite effects on prey diversity depending on the size and productivity of the habitat or the metacommunity richness. Statistical properties of the fluctuations of prey populations determine the effect of stabilizing mechanisms on species richness. We discuss the implications of this result for empirical studies of predation in small study areas and for the management of small nature reserves.

The influence of nutritional history on the functional response of Geocoris pallidipennis to its prey, Myzus persicae

Insect artificial diets are the foundation for mass production of insect predators. Whether there is an influence of long-term rearing with artificial diet on the control ability of predators should be considered. Here, we focused on the effect of nutritional history on the functional response of Geocoris pallidipennis to Myzus persicae. The influence of nutritional history (artificial diet versus natural prey, M. persicae) on the functional response of third to fifth instar nymphs and female G. pallidipennis was examined in the laboratory. The results showed that the functional response curve of both the nymphs and the adult female of G. pallidipennis to M. persicae reflected similar trends on both nutritional histories and confirmed the type II response. Adult female G. pallidipennis reared on either M. persicae or artificial diet produced a significantly better performance than the juvenile stages tested. We estimated that adult female G. pallidipennis can consume 141.6 (artificial diet) or 131.6 (M. persicae) aphids per day, respectively. This indicated that G. pallidipennis reared on both artificial diet and M. persicae displayed high rates of predation.