Selection against instability: stable subgraphs are most frequent in empirical food webs

Food web structure can be characterized by the particular frequencies of subgraphs found within them. Although there are thirteen possible configurations of three species subgraphs, some are consistently over‐represented in empirical food webs. This is a robust pattern that is found across marine, freshwater or terrestrial environments. The preferential elimination of unstable subgraphs during the assembly of the food web can explain the observed pattern. It follows from this hypothesis that there should be differences in the stability of different subgraphs, and that stability should be positively correlated to their frequency in food webs. Using 50 food webs collected from a variety of databases I determined the frequency of each of the thirteen possible subgraphs with respect to randomized webs. Then by numerical simulation I determined the quasi sign stability (QSS) of each subgraph. My results clearly show a positive correlation between QSS and over‐representation of the different subgraphs in empirical food webs.     

Emergence of non-random structure in local food webs generated from randomly structured regional webs

Previous studies have shown that high-resolution, empirical food webs possess a non-random network structure, typically characterized by uniform or exponential degree distributions. However, the empirical food webs that have been investigated for their structural properties represent local communities that are only a subset of a larger pool of regionally coexisting species. Here, we use a simple model to investigate the effects of regional food web structure on local food webs that are assembled by two simple processes: random immigration of species from a source web (regional food web), and random extinction of species within the local web. The model shows that local webs with non-random degree distributions can arise from randomly structured source webs. A comparison of local webs assembled from randomly structured source webs with local webs assembled from source webs generated by the niche model shows that the former have higher species richness at equilibrium, but have a nonlinear response to changing extinction rates. These results imply that the network structure of regional food webs can play a significant role in the assembly and dynamics of local webs in natural ecosystems. With natural landscapes becoming increasingly fragmented, understanding such structure may be a necessary key to understanding the maintenance and stability of local species diversity.     

Persistence increases with diversity and connectance in trophic metacommunities

Background

We are interested in understanding if metacommunity dynamics contribute to the persistence of complex spatial food webs subject to colonization-extinction dynamics. We study persistence as a measure of stability of communities within discrete patches, and ask how do species diversity, connectance, and topology influence it in spatially structured food webs.

Methodology/Principal Findings

We answer this question first by identifying two general mechanisms linking topology of simple food web modules and persistence at the regional scale. We then assess the robustness of these mechanisms to more complex food webs with simulations based on randomly created and empirical webs found in the literature. We find that linkage proximity to primary producers and food web diversity generate a positive relationship between complexity and persistence in spatial food webs. The comparison between empirical and randomly created food webs reveal that the most important element for food web persistence under spatial colonization-extinction dynamics is the degree distribution: the number of prey species per consumer is more important than their identity.

Conclusions/Significance

With a simple set of rules governing patch colonization and extinction, we have predicted that diversity and connectance promote persistence at the regional scale. The strength of our approach is that it reconciles the effect of complexity on stability at the local and the regional scale. Even if complex food webs are locally prone to extinction, we have shown their complexity could also promote their persistence through regional dynamics. The framework we presented here offers a novel and simple approach to understand the complexity of spatial food webs.     

Connectance determines invasion success via trophic interactions in model food webs

Human induced global change has greatly altered the structure and composition of food webs through the invasion of non‐native species and the extinction of native species. Much attention has been paid to the effects of species deletions on food web structure and stability. However, recent empirical evidence suggests that for most taxa local species richness has increased as successful invasions outpace extinctions at this scale. This pattern suggests that food webs, which represent feeding interactions at the local scale, may be increasing in species richness. Knowledge of how food web structure relates to invasive species establishment and the effect of successful invaders on subsequent food web structure remains an unknown but potentially important aspect of global change. Here we explore the effect of food web topology on invasion success in model food webs to develop hypotheses about how the distribution of biodiversity across trophic levels affects the success of invasion at each trophic level. Our results suggest a connectance (C) based framework for predicting invasion success in food webs due to the way that C constrains the number of species at each trophic level and thus the number of potential predators and prey for an invader at a given trophic level. We use the relationship between C and the proportion of species at each trophic level in 14 well studied food webs to make the following predictions; 1) the success of basal invaders will increase as C increases due to the decrease in herbivores in high C webs, 2) herbivore invasion success will decrease as C increases due to the decrease in the proportion of basal species and increase in intermediate species and omnivores in high C webs. 3) Top predator invasion success will increase as C increases due to the increase in intermediate prey species. However, it is not clear how the relative influence of trophic structure compares to empirically known predictors of invasion success such as invader traits, propagule pressure, and resource availability.     

Food-web complexity emerging from ecological dynamics on adaptive networks

Food webs are complex networks describing trophic interactions in ecological communities. Since Robert May's seminal work on random structured food webs, the complexity-stability debate is a central issue in ecology: does network complexity increase or decrease food-web persistence? A multi-species predator-prey model incorporating adaptive predation shows that the action of ecological dynamics on the topology of a food web (whose initial configuration is generated either by the cascade model or by the niche model) render, when a significant fraction of adaptive predators is present, similar hyperbolic complexity-persistence relationships as those observed in empirical food webs. It is also shown that the apparent positive relation between complexity and persistence in food webs generated under the cascade model, which has been pointed out in previous papers, disappears when the final connection is used instead of the initial one to explain species persistence.     

Evidence for the existence of a robust pattern of prey selection in food webs

Food webs aim to provide a thorough representation of the trophic interactions found in an ecosystem. The complexity of empirical food webs, however, is leading many ecologists to focus dynamic ecosystem studies on smaller microcosm or mesocosm studies based upon community modules, which comprise three to five species and the interactions likely to have ecological relevance. We provide here a structural counterpart to community modules. We investigate food-web 'motifs' which are n-species connected subgraphs found within the food web. Remarkably, we find that the over- and under-representation of three-species motifs in empirical food webs can be understood through comparison to a static food-web model, the niche model. Our result conclusively demonstrates that predation upon species with some 'characteristic' niche value is the prey selection mechanism consistent with the structural properties of empirical food webs.     

Current food web models cannot explain the overall topological structure of observed food webs

Topological food webs illustrating “who eats whom” in different systems exhibit similar, non‐random, structures suggesting that general rules govern food web structure. Current food web models correctly predict many measures of food web topology from knowledge of species richness and connectance (fraction of possible predator–prey links that actually occur), together with assumptions about the ecological rules governing “who eats whom”. However, current measures are relatively insensitive to small changes in topology. Improvement of, and discrimination among, current models requires development of new measures of food web structure. Here I examine whether current food web models (cascade, niche, and nested hierarchy models, plus a random null model) can predict a new measure of food web structure, structural stability. Structural stability complements other measures of food web topology because it is sensitive to changes in topology that other measures often miss. The cascade and null models respectively over‐ and underpredict structural stability for a set of 17 high‐quality food webs. While the niche and nested hierarchy models provide unbiased predictions on average, their 95% confidence intervals frequently fail to include the observed data. Observed structural stabilities for all models are overdispersed compared to model predictions, and predicted and observed structural stabilities are uncorrelated, indicating that important sources of variation in structural stability are not captured by the models. Crucially, poor model performance arises because observed variation in structural stability is unrelated to variation in species richness and connectance. In contrast, almost all other measures of food web topology vary with species richness and connectance in natural webs. No model that takes species richness and connectance as the only input parameters can reproduce observed variation in structural stability. Further progress in predicting and explaining food web topology will require fundamentally new models based on different input parameters.     

Predicting food‐web structure with metacommunity models

Synthesis Metacommunity theory aims to elucidate the relative influence of local and regional‐scale processes in generating diversity patterns across the landscape. Metacommunity research has focused largely on assemblages of competing organisms within a single trophic level. Here, we test the ability of metacommunity models to predict the network structure of the aquatic food web found in the leaves of the northern pitcher plant Sarracenia purpurea. The species‐sorting and patch‐dynamics models most accurately reproduced nine food web properties, suggesting that local‐scale interactions play an important role in structuring Sarracenia food webs. Our approach can be applied to any well‐resolved food web for which data are available from multiple locations. The metacommunity framework explores the relative influence of local and regional‐scale processes in generating diversity patterns across the landscape. Metacommunity models and empirical studies have focused mostly on assemblages of competing organisms within a single trophic level. Studies of multi‐trophic metacommunities are predominantly restricted to simplified trophic motifs and rarely consider entire food webs. We tested the ability of the patch‐dynamics, species‐sorting, mass‐effects, and neutral metacommunity models, as well as three hybrid models, to reproduce empirical patterns of food web structure and composition in the complex aquatic food web found in the northern pitcher plant Sarracenia purpurea. We used empirical data to determine regional species pools and estimate dispersal probabilities, simulated local food‐web dynamics, dispersed species from regional pools into local food webs at rates based on the assumptions of each metacommunity model, and tested their relative fits to empirical data on food‐web structure. The species‐sorting and patch‐dynamics models most accurately reproduced nine food web properties, suggesting that local‐scale interactions were important in structuring Sarracenia food webs. However, differences in dispersal abilities were also important in models that accurately reproduced empirical food web properties. Although the models were tested using pitcher‐plant food webs, the approach we have developed can be applied to any well‐resolved food web for which data are available from multiple locations.     

The dynamics of spatially coupled food webs

The dynamics of ecological systems include a bewildering number of biotic interactions that unfold over a vast range of spatial scales. Here, employing simple and general empirical arguments concerning the nature of movement, trophic position and behaviour we outline a general theory concerning the role of space and food web structure on food web stability. We argue that consumers link food webs in space and that this spatial structure combined with relatively rapid behavioural responses by consumers can strongly influence the dynamics of food webs. Employing simple spatially implicit food web models, we show that large mobile consumers are inordinately important in determining the stability, or lack of it, in ecosystems. More specifically, this theory suggests that mobile higher order organisms are potent stabilizers when embedded in a variable, and expansive spatial structure. However, when space is compressed and higher order consumers strongly couple local habitats then mobile consumers can have an inordinate destabilizing effect. Preliminary empirical arguments show consistency with this general theory.     

Dynamics of compartmented and reticulate food webs in relation to energetic flows

Using simple food webs, we address how the interactions of food web structure and energetic flows influence dynamics. We examine the effect of food web topologies with equivalent energetics (i.e., trophic interactions are equivalent at each trophic level), following which we vary energetic flows to include weak and strong interactions or nonequivalent energetics. In contrast to some work (Pimm 1979), we find that compartmented webs are more stable than reticulate webs. However, we find that nonequivalent energetics can stabilize previously unstable reticulate structures. It is not only weak flows that can be stabilizing but also the arrangement of the flows that emphasizes stabilizing mechanisms. We find that the main stabilizing mechanism is asynchrony, where structures and energetic arrangements that decrease synchrony such as internal segregation or competition will stabilize dynamics. Since compartments allow prey dynamics to behave somewhat independently, compartmentation readily promotes stability. In addition, these results can be scaled from simple food webs to more complex webs with many interacting subsystems so that linking weak subsystems to strong ones can stabilize dynamics. We show that food web dynamics are determined not only by topology but also the arrangement of weak and strong energetic flows.     

Complex networks: two ways to be robust?

Recent studies of biological networks have focused on the distribution of the number of links per node. However, the connectivity distribution does not uncover all the complexity of their topology. Here, we analyse the relation between the connectivity of a species and the average connectivity of its nearest neighbours in three of the most resolved community food webs. We compare the pattern arising with the one recently reported for protein networks and for a simple null model of a random network. Whereas two highly connected nodes are unlikely to be connected between each other in protein networks, the reverse happens in food webs. We discuss this difference in organization in relation to the robustness of biological networks to different types of perturbation.     

Network topology: patterns and mechanisms in plant-herbivore and host-parasitoid food webs

1.?Biological communities are organized in complex interaction networks such as food webs, which topology appears to be non-random. Gradients, compartments, nested subsets and even combinations of these structures have been shown in bipartite networks. However, in most studies only one pattern is tested against randomness and mechanistic hypotheses are generally lacking. 2.?Here we examined the topology of regional, coexisting plant-herbivore and host-parasitoid food webs to discriminate between the mentioned network patterns. We also evaluated the role of species body size, local abundance, regional frequency and phylogeny as determinants of network topology. 3.?We found both food webs to be compartmented, with interaction range boundaries imposed by host phylogeny. Species degree within compartments was mostly related to their regional frequency and local abundance. Only one compartment showed an internal nested structure in the distribution of interactions between species, but species position within this compartment was unrelated to species size or abundance. 4.?These results suggest that compartmentalization may be more common than previously considered, and that network structure is a result of multiple, hierarchical, non-exclusive processes.     

Biodiversity lessens the risk of cascading extinction in model food webs

Due to the complex interactions between species in food webs, the extinction of one species could lead to a cascade of further extinctions and hence cause dramatic changes in species composition and ecosystem processes. We found that the risk of additional species extinction, following the loss of one species in model food webs, decreases with the number of species per functional group. For a given number of species per functional group, the risk of further extinctions is highest when an autotroph is removed and lowest when a top predator is removed. In addition, stability decreases when the distribution of interaction strengths in the webs is changed from equal to skew (few strong and many weak links). We also found that omnivory appears to stabilize model food webs. Our results indicate that high biodiversity may serve as an insurance against radical ecosystem changes.     

Motif-Role-Fingerprints: The Building-Blocks of Motifs,Clustering-Coefficients and Transitivities in Directed Networks

Complex networks are frequently characterized by metrics for which particular subgraphs are counted. One statistic from this category, which we refer to as motif-role fingerprints, differs from global subgraph counts in that the number of subgraphs in which each node participates is counted. As with global subgraph counts, it can be important to distinguish between motif-role fingerprints that are ‘structural’ (induced subgraphs) and ‘functional’ (partial subgraphs). Here we show mathematically that a vector of all functional motif-role fingerprints can readily be obtained from an arbitrary directed adjacency matrix, and then converted to structural motif-role fingerprints by multiplying that vector by a specific invertible conversion matrix. This result demonstrates that a unique structural motif-role fingerprint exists for any given functional motif-role fingerprint. We demonstrate a similar result for the cases of functional and structural motif-fingerprints without node roles, and global subgraph counts that form the basis of standard motif analysis. We also explicitly highlight that motif-role fingerprints are elemental to several popular metrics for quantifying the subgraph structure of directed complex networks, including motif distributions, directed clustering coefficient, and transitivity. The relationships between each of these metrics and motif-role fingerprints also suggest new subtypes of directed clustering coefficients and transitivities. Our results have potential utility in analyzing directed synaptic networks constructed from neuronal connectome data, such as in terms of centrality. Other potential applications include anomaly detection in networks, identification of similar networks and identification of similar nodes within networks. Matlab code for calculating all stated metrics following calculation of functional motif-role fingerprints is provided as S1 Matlab File.     

Shrinking by numbers: landscape context affects the species composition but not the quantitative structure of local food webs

1. With habitat fragmentation spreading around the world, there is a pressing need to understand its impacts on local food webs. To date, few studies have examined the effects of landscape context on multiple local communities in a quantitative, spatially realistic setting. 2. To examine how the isolation of a food web affects its structure, we construct local food webs of specialist herbivores and their natural enemies on 82 individual oaks (Quercus robur) growing in different landscape contexts. 3. Across this set of webs, we find that communities in isolated habitat patches not only contained fewer species than did well-connected ones, but also differed in species composition. 4. Surprisingly, the effects observed in terms of species composition were not reflected in the quantitative interaction structure of local food webs: landscape context had no detectable effect on either the interaction evenness, linkage density, connectance, generality or vulnerability of local webs. 5. We conclude that the quantitative structure of food webs may be stable in the face of habitat fragmentation, despite clear-cut impacts on individual species. This finding offers hope-inspiring news for conservation, but should clearly be verified by empirical studies across both naturally and more recently fragmented systems.     

Food webs: Ordering species according to body size yields high degree of intervality

Food webs, the networks describing “who eats whom” in an ecosystem, are nearly interval, i.e. there is a way to order the species so that almost all the resources of each consumer are adjacent in the ordering. This feature has important consequences, as it means that the structure of food webs can be described using a single (or few) species' traits. Moreover, exploiting the quasi-intervality found in empirical webs can help build better models for food web structure. Here we investigate which species trait is a good proxy for ordering the species to produce quasi-interval orderings. We find that body size produces a significant degree of intervality in almost all food webs analyzed, although it does not match the maximum intervality for the networks. There is also a great variability between webs. Other orderings based on trophic levels produce a lower level of intervality. Finally, we extend the concept of intervality from predator-centered (in which resources are in intervals) to prey-centered (in which consumers are in intervals). In this case as well we find that body size yields a significant, but not maximal, level of intervality. These results show that body size is an important, although not perfect, trait that shapes species interactions in food webs. This has important implications for the formulation of simple models used to construct realistic representations of food webs.     

Using trophic hierarchy to understand food web structure

Link arrangement in food webs is determined by the species' feeding habits. This work investigates whether food web topology is organized in a gradient of trophic positions from producers to consumers. To this end, we analyzed 26 food webs for which the consumption rate of each species was specified. We computed the trophic positions and the link densities of all species in the food webs. Link density measures how much each species contributes to the distribution of energy in the system. It is expressed as the number of links species establish with other nodes, weighted by their magnitude. We computed these two metrics using various formulations developed in the ecological network analysis framework. Results show a positive correlation between trophic position and link density across all the systems, regardless the specific formulas used to measure the two quantities. We performed the same analysis on the corresponding binary matrices (i.e. removing information about rates). In addition, we investigated the relation between trophic position and link density in: a) simulated binary webs with same connectance as the original ones; b) weighted webs with constant topology but randomized link strengths and c) weighted webs with constant connectance where both topology and link strengths are randomized. The correlation between the two indices attenuates, vanishes or becomes negative in the case of binary food webs and simulated data (weighted and unweighted).
According to our analysis, link density in food webs decreases with trophic position so that it is greatly reduced toward the top of the trophic hierarchy. This outcome, that seems to challenge previous conclusions based on null models, strongly depends on link quantification. Including interaction strengths may improve substantially our understanding of food web organization, and possibly contradict results based on the analysis of binary webs.     

Compilation and network analyses of cambrian food webs

A rich body of empirically grounded theory has developed about food webs—the networks of feeding relationships among species within habitats. However, detailed food-web data and analyses are lacking for ancient ecosystems, largely because of the low resolution of taxa coupled with uncertain and incomplete information about feeding interactions. These impediments appear insurmountable for most fossil assemblages; however, a few assemblages with excellent soft-body preservation across trophic levels are candidates for food-web data compilation and topological analysis. Here we present plausible, detailed food webs for the Chengjiang and Burgess Shale assemblages from the Cambrian Period. Analyses of degree distributions and other structural network properties, including sensitivity analyses of the effects of uncertainty associated with Cambrian diet designations, suggest that these early Paleozoic communities share remarkably similar topology with modern food webs. Observed regularities reflect a systematic dependence of structure on the numbers of taxa and links in a web. Most aspects of Cambrian food-web structure are well-characterized by a simple “niche model,” which was developed for modern food webs and takes into account this scale dependence. However, a few aspects of topology differ between the ancient and recent webs: longer path lengths between species and more species in feeding loops in the earlier Chengjiang web, and higher variability in the number of links per species for both Cambrian webs. Our results are relatively insensitive to the exclusion of low-certainty or random links. The many similarities between Cambrian and recent food webs point toward surprisingly strong and enduring constraints on the organization of complex feeding interactions among metazoan species. The few differences could reflect a transition to more strongly integrated and constrained trophic organization within ecosystems following the rapid diversification of species, body plans, and trophic roles during the Cambrian radiation. More research is needed to explore the generality of food-web structure through deep time and across habitats, especially to investigate potential mechanisms that could give rise to similar structure, as well as any differences.