Phage strategies facilitate bacterial coexistence under environmental variability.

Bacterial communities are often exposed to temporal variations in resource availability, which exceed bacterial generation times and thereby affect bacterial coexistence. Bacterial population dynamics are also shaped by bacteriophages, which are a main cause of bacterial mortality. Several strategies are proposed in the literature to describe infections by phages, such as "Killing the Winner", "Piggyback the loser" (PtL) or "Piggyback the Winner" (PtW). The two temperate phage strategies PtL and PtW are defined by a change from lytic to lysogenic infection when the host density changes, from high to low or from low to high, respectively. To date, the occurrence of different phage strategies and their response to environmental variability is poorly understood. In our study, we developed a microbial trophic network model using ordinary differential equations (ODEs) and performed 'in silico' experiments. To model the switch from the lysogenic to the lytic cycle, we modified the lysis rate of infected bacteria and their growth was turned on or off using a density-dependent switching point. We addressed whether and how the different phage strategies facilitate bacteria coexistence competing for limiting resources. We also studied the impact of a fluctuating resource inflow to evaluate the response of the different phage strategies to environmental variability. Our results show that the viral shunt (i.e. nutrient release after bacterial lysis) leads to an enrichment of the system. This enrichment enables bacterial coexistence at lower resource concentrations. We were able to show that an established, purely lytic model leads to stable bacterial coexistence despite fluctuating resources. Both temperate phage models differ in their coexistence patterns. The model of PtW yields stable bacterial coexistence at a limited range of resource supply and is most sensitive to resource fluctuations. Interestingly, the purely lytic phage strategy and PtW both result in stable bacteria coexistence at oligotrophic conditions. The PtL model facilitates stable bacterial coexistence over a large range of stable and fluctuating resource inflow. An increase in bacterial growth rate results in a higher resilience to resource variability for the PtL and the lytic infection model. We propose that both temperate phage strategies represent different mechanisms of phages coping with environmental variability. Our study demonstrates how phage strategies can maintain bacterial coexistence in constant and fluctuating environments.

Thermal acclimation increases the stability of a predator-prey interaction in warmer environments.

Global warming over the next century is likely to alter the energy demands of consumers and thus the strengths of their interactions with their resources. The subsequent cascading effects on population biomasses could have profound effects on food web stability. One key mechanism by which organisms can cope with a changing environment is phenotypic plasticity, such as acclimation to warmer conditions through reversible changes in their physiology. Here, we measured metabolic rates and functional responses in laboratory experiments for a widespread predator-prey pair of freshwater invertebrates, sampled from across a natural stream temperature gradient in Iceland (4-18℃). This enabled us to parameterize a Rosenzweig-MacArthur population dynamical model to study the effect of thermal acclimation on the persistence of the predator-prey pairs in response to warming. Acclimation to higher temperatures either had neutral effects or reduced the thermal sensitivity of both metabolic and feeding rates for the predator, increasing its energetic efficiency. This resulted in greater stability of population dynamics, as acclimation to higher temperatures increased the biomass of both predator and prey populations with warming. These findings indicate that phenotypic plasticity can act as a buffer against the impacts of environmental warming. As a consequence, predator-prey interactions between ectotherms may be less sensitive to future warming than previously expected, but this requires further investigation across a broader range of interacting species.

The biggest losers: habitat isolation deconstructs complex food webs from top to bottom.

Habitat fragmentation threatens global biodiversity. To date, there is only limited understanding of how the different aspects of habitat fragmentation (habitat loss, number of fragments and isolation) affect species diversity within complex ecological networks such as food webs. Here, we present a dynamic and spatially explicit food web model which integrates complex food web dynamics at the local scale and species-specific dispersal dynamics at the landscape scale, allowing us to study the interplay of local and spatial processes in metacommunities. We here explore how the number of habitat patches, i.e. the number of fragments, and an increase of habitat isolation affect the species diversity patterns of complex food webs (α-, ß-, γ-diversities). We specifically test whether there is a trophic dependency in the effect of these two factors on species diversity. In our model, habitat isolation is the main driver causing species loss and diversity decline. Our results emphasize that large-bodied consumer species at high trophic positions go extinct faster than smaller species at lower trophic levels, despite being superior dispersers that connect fragmented landscapes better. We attribute the loss of top species to a combined effect of higher biomass loss during dispersal with increasing habitat isolation in general, and the associated energy limitation in highly fragmented landscapes, preventing higher trophic levels to persist. To maintain trophic-complex and species-rich communities calls for effective conservation planning which considers the interdependence of trophic and spatial dynamics as well as the spatial context of a landscape and its energy availability.

Consistent temperature dependence of functional response parameters and their use in predicting population abundance.

Global warming is one of the greatest threats to the persistence of populations: increased metabolic demands should strengthen pairwise species interactions, which could destabilize food webs at the higher organizational levels. Quantifying the temperature dependence of consumer-resource interactions is thus essential for predicting ecological responses to warming. We explored feeding interactions between different predator-prey pairs in controlled-temperature chambers and in a system of naturally heated streams. We found consistent temperature dependence of attack rates across experimental settings, though the magnitude and activation energy of attack rate were specific to each predator, which varied in mobility and foraging mode. We used these parameters along with metabolic rate measurements to estimate energetic efficiency and population abundance with warming. Energetic efficiency accurately estimated field abundance of a mobile predator that struggled to meet its metabolic demands, but was a poor predictor for a sedentary predator that operated well below its energetic limits. Temperature effects on population abundance may thus be strongly dependent on whether organisms are regulated by their own energy intake or interspecific interactions. Given the widespread use of functional response parameters in ecological modelling, reconciling outcomes from laboratory and field studies increases the confidence and precision with which we can predict warming impacts on natural systems.

Predator traits determine food-web architecture across ecosystems.

Predator-prey interactions in natural ecosystems generate complex food webs that have a simple universal body-size architecture where predators are systematically larger than their prey. Food-web theory shows that the highest predator-prey body-mass ratios found in natural food webs may be especially important because they create weak interactions with slow dynamics that stabilize communities against perturbations and maintain ecosystem functioning. Identifying these vital interactions in real communities typically requires arduous identification of interactions in complex food webs. Here, we overcome this obstacle by developing predator-trait models to predict average body-mass ratios based on a database comprising 290 food webs from freshwater, marine and terrestrial ecosystems across all continents. We analysed how species traits constrain body-size architecture by changing the slope of the predator-prey body-mass scaling. Across ecosystems, we found high body-mass ratios for predator groups with specific trait combinations including (1) small vertebrates and (2) large swimming or flying predators. Including the metabolic and movement types of predators increased the accuracy of predicting which species are engaged in high body-mass ratio interactions. We demonstrate that species traits explain striking patterns in the body-size architecture of natural food webs that underpin the stability and functioning of ecosystems, paving the way for community-level management of the most complex natural ecosystems.

Bridging Scales: Allometric Random Walks Link Movement and Biodiversity Research.

Integrating mechanistic models of movement and behavior into large-scale movement ecology and biodiversity research is one of the major challenges in current ecological science. This is mainly due to a large gap between the spatial scales at which these research lines act. Here, we propose to apply trait-based movement models to bridge this gap and generalize movement trajectories across species and ecosystems. We show how to use species traits (e.g., body mass) to generate allometric random walks and illustrate in two worked examples how this facilitates general predictions of species-interaction traits, meta-community structures, and biodiversity patterns. Thereby, allometric random walks foster a closer integration of movement ecology and biodiversity research by scaling up from small-scale mechanistic measurements to a predictive understanding of movement and biodiversity patterns in different landscapes.

Applying generalized allometric regressions to predict live body mass of tropical and temperate arthropods.

The ecological implications of body size extend from the biology of individual organisms to ecosystem-level processes. Measuring body mass for high numbers of invertebrates can be logistically challenging, making length-mass regressions useful for predicting body mass with minimal effort. However, standardized sets of scaling relationships covering a large range in body length, taxonomic groups, and multiple geographical regions are scarce. We collected 6,212 arthropods from 19 higher-level taxa in both temperate and tropical locations to compile a comprehensive set of linear models relating live body mass to a range of predictor variables. We measured live weight (hereafter, body mass), body length and width of each individual and conducted linear regressions to predict body mass using body length, body width, taxonomic group, and geographic region. Additionally, we quantified prediction discrepancy when using parameters from arthropods of a different geographic region. Incorporating body width into taxon- and region-specific length-mass regressions yielded the highest prediction accuracy for body mass. Using regression parameters from a different geographic region increased prediction discrepancy, causing over- or underestimation of body mass depending on geographical origin and whether body width was included. We present a comprehensive range of parameters for predicting arthropod body mass and provide guidance for selecting optimal scaling relationships. Given the importance of body mass for functional invertebrate ecology and the paucity of adequate regressions to predict arthropod body mass from different geographical regions, our study provides a long-needed resource for quantifying live body mass in invertebrate ecology research.

A general scaling law reveals why the largest animals are not the fastest.

Speed is the fundamental constraint on animal movement, yet there is no general consensus on the determinants of maximum speed itself. Here, we provide a general scaling model of maximum speed with body mass, which holds across locomotion modes, ecosystem types and taxonomic groups. In contrast to traditional power-law scaling, we predict a hump-shaped relationship resulting from a finite acceleration time for animals, which explains why the largest animals are not the fastest. This model is strongly supported by extensive empirical data (474 species, with body masses ranging from 30 µg to 100 tonnes) from terrestrial as well as aquatic ecosystems. Our approach unravels a fundamental constraint on the upper limit of animal movement, thus enabling a better understanding of realized movement patterns in nature and their multifold ecological consequences.

How patch size and refuge availability change interaction strength and population dynamics: a combined individual- and population-based modeling experiment.

Knowledge on how functional responses (a measurement of feeding interaction strength) are affected by patch size and habitat complexity (represented by refuge availability) is crucial for understanding food-web stability and subsequently biodiversity. Due to their laborious character, it is almost impossible to carry out systematic empirical experiments on functional responses across wide gradients of patch sizes and refuge availabilities. Here we overcame this issue by using an individual-based model (IBM) to simulate feeding experiments. The model is based on empirically measured traits such as body-mass dependent speed and capture success. We simulated these experiments in patches ranging from sizes of petri dishes to natural patches in the field. Moreover, we varied the refuge availability within the patch independently of patch size, allowing for independent analyses of both variables. The maximum feeding rate (the maximum number of prey a predator can consume in a given time frame) is independent of patch size and refuge availability, as it is the physiological upper limit of feeding rates. Moreover, the results of these simulations revealed that a type III functional response, which is known to have a stabilizing effect on population dynamics, fitted the data best. The half saturation density (the prey density where a predator consumes half of its maximum feeding rate) increased with refuge availability but was only marginally influenced by patch size. Subsequently, we investigated how patch size and refuge availability influenced stability and coexistence of predator-prey systems. Following common practice, we used an allometric scaled Rosenzweig-MacArthur predator-prey model based on results from our in silico IBM experiments. The results suggested that densities of both populations are nearly constant across the range of patch sizes simulated, resulting from the constant interaction strength across the patch sizes. However, constant densities with decreasing patch sizes mean a decrease of absolute number of individuals, consequently leading to extinction of predators in the smallest patches. Moreover, increasing refuge availabilities also allowed predator and prey to coexist by decreased interaction strengths. Our results underline the need for protecting large patches with high habitat complexity to sustain biodiversity.

Predicting the consequences of species loss using size-structured biodiversity approaches.

Understanding the consequences of species loss in complex ecological communities is one of the great challenges in current biodiversity research. For a long time, this topic has been addressed by traditional biodiversity experiments. Most of these approaches treat species as trait-free, taxonomic units characterizing communities only by species number without accounting for species traits. However, extinctions do not occur at random as there is a clear correlation between extinction risk and species traits. In this review, we assume that large species will be most threatened by extinction and use novel allometric and size-spectrum concepts that include body mass as a primary species trait at the levels of populations and individuals, respectively, to re-assess three classic debates on the relationships between biodiversity and (i) food-web structural complexity, (ii) community dynamic stability, and (iii) ecosystem functioning. Contrasting current expectations, size-structured approaches suggest that the loss of large species, that typically exploit most resource species, may lead to future food webs that are less interwoven and more structured by chains of interactions and compartments. The disruption of natural body-mass distributions maintaining food-web stability may trigger avalanches of secondary extinctions and strong trophic cascades with expected knock-on effects on the functionality of the ecosystems. Therefore, we argue that it is crucial to take into account body size as a species trait when analysing the consequences of biodiversity loss for natural ecosystems. Applying size-structured approaches provides an integrative ecological concept that enables a better understanding of each species' unique role across communities and the causes and consequences of biodiversity loss.

Analyzing pathogen suppressiveness in bioassays with natural soils using integrative maximum likelihood methods in R.

The potential of soils to naturally suppress inherent plant pathogens is an important ecosystem function. Usually, pathogen infection assays are used for estimating the suppressive potential of soils. In natural soils, however, co-occurring pathogens might simultaneously infect plants complicating the estimation of a focal pathogen's infection rate (initial slope of the infection-curve) as a measure of soil suppressiveness. Here, we present a method in R correcting for these unwanted effects by developing a two pathogen mono-molecular infection model. We fit the two pathogen mono-molecular infection model to data by using an integrative approach combining a numerical simulation of the model with an iterative maximum likelihood fit. We show that in presence of co-occurring pathogens using uncorrected data leads to a critical under- or overestimation of soil suppressiveness measures. In contrast, our new approach enables to precisely estimate soil suppressiveness measures such as plant infection rate and plant resistance time. Our method allows a correction of measured infection parameters that is necessary in case different pathogens are present. Moreover, our model can be (1) adapted to use other models such as the logistic or the Gompertz model; and (2) it could be extended by a facilitation parameter if infections in plants increase the susceptibility to new infections. We propose our method to be particularly useful for exploring soil suppressiveness of natural soils from different sites (e.g., in biodiversity experiments).

Animal diversity and ecosystem functioning in dynamic food webs.

Species diversity is changing globally and locally, but the complexity of ecological communities hampers a general understanding of the consequences of animal species loss on ecosystem functioning. High animal diversity increases complementarity of herbivores but also increases feeding rates within the consumer guild. Depending on the balance of these counteracting mechanisms, species-rich animal communities may put plants under top-down control or may release them from grazing pressure. Using a dynamic food-web model with body-mass constraints, we simulate ecosystem functions of 20,000 communities of varying animal diversity. We show that diverse animal communities accumulate more biomass and are more exploitative on plants, despite their higher rates of intra-guild predation. However, they do not reduce plant biomass because the communities are composed of larger, and thus energetically more efficient, plant and animal species. This plasticity of community body-size structure reconciles the debate on the consequences of animal species loss for primary productivity.

Interactive effects of warming, eutrophication and size structure: impacts on biodiversity and food-web structure.

Warming and eutrophication are two of the most important global change stressors for natural ecosystems, but their interaction is poorly understood. We used a dynamic model of complex, size-structured food webs to assess interactive effects on diversity and network structure. We found antagonistic impacts: Warming increases diversity in eutrophic systems and decreases it in oligotrophic systems. These effects interact with the community size structure: Communities of similarly sized species such as parasitoid-host systems are stabilized by warming and destabilized by eutrophication, whereas the diversity of size-structured predator-prey networks decreases strongly with warming, but decreases only weakly with eutrophication. Nonrandom extinction risks for generalists and specialists lead to higher connectance in networks without size structure and lower connectance in size-structured communities. Overall, our results unravel interactive impacts of warming and eutrophication and suggest that size structure may serve as an important proxy for predicting the community sensitivity to these global change stressors.

Unifying elemental stoichiometry and metabolic theory in predicting species abundances.

While metabolic theory predicts variance in population density within communities depending on population average body masses, the ecological stoichiometry concept relates density variation across communities to varying resource stoichiometry. Using a data set including biomass densities of 4959 populations of soil invertebrates across 48 forest sites we combined these two frameworks. We analyzed how the scaling of biomass densities with population-averaged body masses systematically interacts with stoichiometric variables. Simplified analyses employing either only body masses or only resource stoichiometry are highly context sensitive and yield variable and often misleading results. Our findings provide strong evidence that analyses of ecological state variables should integrate allometric and stoichiometric variables to explain deviations from predicted allometric scaling and avoid erroneous conclusions. In consequence, our study provides an important step towards unifying two prominent ecological theories, metabolic theory and ecological stoichiometry.

Body masses, functional responses and predator-prey stability.

The stability of ecological communities depends strongly on quantitative characteristics of population interactions (type-II vs. type-III functional responses) and the distribution of body masses across species. Until now, these two aspects have almost exclusively been treated separately leaving a substantial gap in our general understanding of food webs. We analysed a large data set of arthropod feeding rates and found that all functional-response parameters depend on the body masses of predator and prey. Thus, we propose generalised functional responses which predict gradual shifts from type-II predation of small predators on equally sized prey to type-III functional-responses of large predators on small prey. Models including these generalised functional responses predict population dynamics and persistence only depending on predator and prey body masses, and we show that these predictions are strongly supported by empirical data on forest soil food webs. These results help unravelling systematic relationships between quantitative population interactions and large-scale community patterns.

Universal temperature and body-mass scaling of feeding rates.

Knowledge of feeding rates is the basis to understand interaction strength and subsequently the stability of ecosystems and biodiversity. Feeding rates, as all biological rates, depend on consumer and resource body masses and environmental temperature. Despite five decades of research on functional responses as quantitative models of feeding rates, a unifying framework of how they scale with body masses and temperature is still lacking. This is perplexing, considering that the strength of functional responses (i.e. interaction strengths) is crucially important for the stability of simple consumer-resource systems and the persistence, sustainability and biodiversity of complex communities. Here, we present the largest currently available database on functional response parameters and their scaling with body mass and temperature. Moreover, these data are integrated across ecosystems and metabolic types of species. Surprisingly, we found general temperature dependencies that differed from the Arrhenius terms predicted by metabolic models. Additionally, the body-mass-scaling relationships were more complex than expected and differed across ecosystems and metabolic types. At local scales (taxonomically narrow groups of consumer-resource pairs), we found hump-shaped deviations from the temperature and body-mass-scaling relationships. Despite the complexity of our results, these body-mass- and temperature-scaling models remain useful as a mechanistic basis for predicting the consequences of warming for interaction strengths, population dynamics and network stability across communities differing in their size structure.

The dynamics of food chains under climate change and nutrient enrichment.

Warming has profound effects on biological rates such as metabolism, growth, feeding and death of organisms, eventually affecting their ability to survive. Using a nonlinear bioenergetic population-dynamic model that accounts for temperature and body-mass dependencies of biological rates, we analysed the individual and interactive effects of increasing temperature and nutrient enrichment on the dynamics of a three-species food chain. At low temperatures, warming counteracts the destabilizing effects of enrichment by both bottom-up (via the carrying capacity) and top-down (via biological rates) mechanisms. Together with increasing consumer body masses, warming increases the system tolerance to fertilization. Simultaneously, warming increases the risk of starvation for large species in low-fertility systems. This effect can be counteracted by increased fertilization. In combination, therefore, two main drivers of global change and biodiversity loss can have positive and negative effects on food chain stability. Our model incorporates the most recent empirical data and may thus be used as the basis for more complex forecasting models incorporating food-web structure.

Climate change effects on macrofaunal litter decomposition: the interplay of temperature, body masses and stoichiometry.

Macrofauna invertebrates of forest floors provide important functions in the decomposition process of soil organic matter, which is affected by the nutrient stoichiometry of the leaf litter. Climate change effects on forest ecosystems include warming and decreasing litter quality (e.g. higher C : nutrient ratios) induced by higher atmospheric CO(2) concentrations. While litter-bag experiments unravelled separate effects, a mechanistic understanding of how interactions between temperature and litter stoichiometry are driving decomposition rates is lacking. In a laboratory experiment, we filled this void by quantifying decomposer consumption rates analogous to predator-prey functional responses that include the mechanistic parameters handling time and attack rate. Systematically, we varied the body masses of isopods, the environmental temperature and the resource between poor (hornbeam) and good quality (ash). We found that attack rates increased and handling times decreased (i) with body masses and (ii) temperature. Interestingly, these relationships interacted with litter quality: small isopods possibly avoided the poorer resource, whereas large isopods exhibited increased, compensatory feeding of the poorer resource, which may be explained by their higher metabolic demands. The combination of metabolic theory and ecological stoichiometry provided critically important mechanistic insights into how warming and varying litter quality may modify macrofaunal decomposition rates.

Warming effects on consumption and intraspecific interference competition depend on predator metabolism.

1. Model analyses show that the stability of population dynamics and food web persistence increase with the strength of interference competition. Despite this critical importance for community stability, little is known about how external factors such as the environmental temperature affect intraspecific interference competition. 2. We aimed to fill this void by studying the functional responses of two ground beetle species of different body size, Pterostichus melanarius and Poecilus versicolor. These functional response experiments were replicated across four predator densities and two temperatures to address the impact of temperature on intraspecific interference competition. 3. We generally expected that warming should increase the speed of movement, encounter rates and in consequence interference among predator individuals. In our experiment, this expectation was supported by the results obtained for the larger predator, P. melanarius, whereas the opposite pattern characterized the interference behaviour of the smaller predator P. versicolor. 4. These results suggest potentially nontrivial implications for the effects of environmental temperature on intraspecific interference competition, for which we propose an explanation based on the different sensitivity to warming of metabolic rates of both species. As expected, increasing temperature led to stronger interference competition of the larger species, P. melanarius, which exhibited a weaker increase in metabolic rate with increasing temperature. The stronger increase in the metabolic rate of the smaller predator, P. versicolor, had to be compensated by increasing searching activity for prey, which did not leave time for increasing interference. 5. Together, these results suggest that any generalization how interference competition responds to warming should also take the species' metabolic response to temperature increases into account.

Phylogenetic grouping, curvature and metabolic scaling in terrestrial invertebrates.

For more than a century, the scaling of animal metabolic rates with individual body masses and environmental temperature has predominantly been described by power-law and exponential relationships respectively. Many theories have been proposed to explain these scaling relationships, but were challenged by empirically documented curvatures on double-logarithmic scales. In the present study, we present a novel data set comprising 3661 terrestrial (mainly soil) invertebrate respiration rates from 192 independent sources across a wide range in body masses, environmental temperatures and phylogenetic groups. Although our analyses documented power-law and exponential scaling with body masses and temperature, respectively, polynomial models identified curved deviations. Interestingly, complex scaling models accounting for phylogenetic groups were able to remove curvatures except for a negative curvature at the highest temperatures (>30 °C) indicating metabolic down regulation. This might indicate that the tremendous differences in invertebrate body architectures, ecology and physiology may cause severely different metabolic scaling processes.