Competition and evolution along environmental gradients: patterns, boundaries and sympatric divergence

The present study examined how competitive interactions and environmental conditions generate species boundaries and determine species distributions. A spatially explicit, quantitative genetic, two-species competition model was used to manipulate the strengths of competition, gene flow and local adaptation along environmental gradients. This allowed us to assess the long-term persistence of each species and whether the ranges they inhabited had boundaries in space or were unlimited. We found that a species boundary arises along less steep environmental gradients when the strength of stabilizing selection and diversifying selection are similar. We also found that a species boundary may arise along shallow environmental gradients if interspecific competition is more intense than intraspecific, which relaxes previous requirements for steep gradients for generating range limits. We determined an analytical form for the critical environmental gradient as a function of ecological and genetic parameters at which a species boundary is expected to arise by competition. Results suggest an alternative to resource competition as an explanation for phenotypic divergence between sympatric competitors. Competitors sharing a trait that is under stabilizing selection along an environmental gradient may segregate spatially and evolve in different regions, with phenotypic sympatric divergence reflecting the resulting clines. Along various types of environmental gradients, variation in stabilizing selection intensities could lead to contrasting patterns in the distribution of species. For stabilizing selection strengths in accord with field data estimates, this study predicts that the level of sympatric character divergence would be limited along environmental gradients.     

Shared environmental responses drive co‐occurrence patterns in river bird communities

Positive or negative patterns of co‐occurrence might imply an influence of biotic interactions on community structure. However, species may co‐occur simply because of shared environmental responses. Here, we apply two complementary modelling methodologies – a probabilistic model of significant pairwise associations and a hierarchical multivariate probit regression model – to 1) attribute co‐occurrence patterns in 100 river bird communities to either shared environmental responses or to other ecological mechanisms such as interaction with heterospecifics, and 2) examine the strength of evidence for four alternative models of community structure. Species co‐occurred more often than would be expected by random community assembly and the species composition of bird communities was highly structured. Co‐occurrence patterns were primarily explained by shared environmental responses; species’ responses to the environmental variables were highly divergent, with both strong positive and negative environmental correlations occurring. We found limited evidence for behaviour‐driven assemblage patterns in bird communities at a large spatial scale, although statistically significant positive associations amongst some species suggested the operation of facilitative mechanisms such as heterospecific attraction. This lends support to an environmental filtering model of community assembly as being the principle mechanism shaping river bird community structure. Consequently, species interactions may be reduced to an ancillary role in some avifaunal communities, meaning if shared environmental responses are not quantified studies of co‐occurrence may overestimate the role of species interactions in shaping community structure.     

Biodiversity and species interactions: extending Lotka–Volterra community theory

A new analysis of the nearly century‐old Lotka–Volterra theory allows us to link species interactions to biodiversity patterns, including: species abundance distributions, estimates of total community size, patterns of community invasibility, and predicted responses to disturbance. Based on a few restrictive assumptions about species interactions, our calculations require only that the community is sufficiently large to allow a mean‐field approximation. We develop this analysis to show how an initial assemblage of species with varying interaction strengths is predicted to sort out into the final community based on the species’ predicted target densities. The sorting process yields predictions of covarying patterns of species abundance, community size, and species interaction strengths. These predictions can be tested using enrichment experiments, examination of latitudinal and productivity gradients, and features of community assembly.     

Predator diversity and environmental change modify the strengths of trophic and nontrophic interactions

Understanding the dependence of species interaction strengths on environmental factors and species diversity is crucial to predict community dynamics and persistence in a rapidly changing world. Nontrophic (e.g. predator interference) and trophic components together determine species interaction strengths, but the effects of environmental factors on these two components remain largely unknown. This impedes our ability to fully understand the links between environmental drivers and species interactions. Here, we used a dynamical modelling framework based on measured predator functional responses to investigate the effects of predator diversity, prey density, and temperature on trophic and nontrophic interaction strengths within a freshwater food web. We found that (i) species interaction strengths cannot be predicted from trophic interactions alone, (ii) nontrophic interaction strengths vary strongly among predator assemblages, (iii) temperature has opposite effects on trophic and nontrophic interaction strengths, and (iv) trophic interaction strengths decrease with prey density, whereas the dependence of nontrophic interaction strengths on prey density is concave up. Interestingly, the qualitative impacts of temperature and prey density on the strengths of trophic and nontrophic interactions were independent of predator identity, suggesting a general pattern. Our results indicate that taking multiple environmental factors and the nonlinearity of density‐dependent species interactions into account is an important step towards a better understanding of the effects of environmental variations on complex ecological communities. The functional response approach used in this study opens new avenues for (i) the quantification of the relative importance of the trophic and nontrophic components in species interactions and (ii) a better understanding how environmental factors affect these interactions and the dynamics of ecological communities.     

Individual phenotypic variation reduces interaction strengths in a consumer–resource system

Natural populations often show variation in traits that can affect the strength of interspecific interactions. Interaction strengths in turn influence the fate of pairwise interacting populations and the stability of food webs. Understanding the mechanisms relating individual phenotypic variation to interaction strengths is thus central to assess how trait variation affects population and community dynamics. We incorporated nonheritable variation in attack rates and handling times into a classical consumer–resource model to investigate how variation may alter interaction strengths, population dynamics, species persistence, and invasiveness. We found that individual variation influences species persistence through its effect on interaction strengths. In many scenarios, interaction strengths decrease with variation, which in turn affects species coexistence and stability. Because environmental change alters the direction and strength of selection acting upon phenotypic traits, our results have implications for species coexistence in a context of habitat fragmentation, climate change, and the arrival of exotic species to native ecosystems.     

Dynamics of species interaction strength in space, time and with developmental stage

Quantifying species interaction strengths enhances prediction of community dynamics, but variability in the strength of species interactions in space and time complicates accurate prediction. Interaction strengths can vary in response to density, indirect effects, priority effects or a changing environment, but the mechanism(s) causing direction and magnitudes of change are often unclear. We designed an experiment to characterize how environmental factors influence the direction and the strength of priority effects between sessile species. We estimated per capita non-trophic effects of barnacles (Semibalanus balanoides) on newly settled germlings of the fucoid, Ascophyllum nodosum, in the presence and absence of consumers in experiments on rocky shores throughout the Gulf of Maine, USA. Per capita effects on germlings varied among environments and barnacle life stages, and these interaction strengths were largely unaltered by changing consumer abundance. Whereas previous evidence shows adult barnacles facilitate fucoids, here, we show that recent settlers and established juveniles initially compete with germlings. As barnacles mature, they switch to become facilitators of fucoids. Consumers caused variable mortality of germlings through time comparable to that from competition. Temporally variable effects of interactors (e.g. S. balanoides), or spatial variation in their population structure, in different regions differentially affect target populations (e.g. A. nodosum). This may affect abundance of critical stages and the resilience of target species to environmental change in different geographical regions.     

Local and collective transitions in sparsely-interacting ecological communities

Interactions in natural communities can be highly heterogeneous, with any given species interacting appreciably with only some of the others, a situation commonly represented by sparse interaction networks. We study the consequences of sparse competitive interactions, in a theoretical model of a community assembled from a species pool. We find that communities can be in a number of different regimes, depending on the interaction strength. When interactions are strong, the network of coexisting species breaks up into small subgraphs, while for weaker interactions these graphs are larger and more complex, eventually encompassing all species. This process is driven by the emergence of new allowed subgraphs as interaction strength decreases, leading to sharp changes in diversity and other community properties, and at weaker interactions to two distinct collective transitions: a percolation transition, and a transition between having a unique equilibrium and having multiple alternative equilibria. Understanding community structure is thus made up of two parts: first, finding which subgraphs are allowed at a given interaction strength, and secondly, a discrete problem of matching these structures over the entire community. In a shift from the focus of many previous theories, these different regimes can be traversed by modifying the interaction strength alone, without need for heterogeneity in either interaction strengths or the number of competitors per species.     

Patterns of species interaction strength in assembled theoretical competition communities

Strength of interactions between species may be an important tool in our effort to understand community structure. Recent theoretical and empirical findings suggest that despite the presence of some strong interactions, weak interactions prevail in communities. Here, we examine how mean interaction strengths change as theoretical competition communities assemble and what the distribution of interaction coefficients is in the communities that are formed during the assembly process. Our results show that the mean competition strengths fall as assembly progresses and that most interactions in the communities formed are weak. Communities that are invulnerable to further invasions are those where interspecific interactions are weaker than the average interaction strength between species in the pool. If these results can be generalized to more than one trophic level, implications for management and conservation of natural communities are substantial.     

Global change alters the stability of food webs

Recent research has generally shown that a small change in the number of species in a food web can have consequences both for community structure and ecosystem processes. However ‘change’ is not limited to just the number of species in a community, but might include an alteration to such properties as precipitation, nutrient cycling and temperature. How such changes might affect species interactions is important, not just through the presence or absence of interactions, but also because the patterning of interaction strengths among species is intimately associated with community stability. Interaction strengths encompass such properties as feeding rates and assimilation efficiencies, and encapsulate functionally important information with regard to ecosystem processes. Interaction strengths represent the pathways and transfer of energy through an ecosystem. We review the best empirical data available detailing the frequency distribution of interaction strengths in communities. We present the underlying (but consistent) pattern of species interactions and discuss the implications of this patterning. We then examine how such a basic pattern might be affected given various scenarios of ‘change’ and discuss the consequences for community stability and ecosystem functioning.     

Canopy Openness Enhances Diversity of Ant–Plant Interactions in the Brazilian Amazon Rain Forest

In closed‐canopy tropical forest understory, light availability is a significant determinant of habitat diversity because canopy structure is highly variable in most tropical forests. Consequently, variation in canopy cover affects the composition and distribution of plant species via creating variable light environments. Nevertheless, little is known about how variation in canopy openness structures patterns of plant–animal interactions. Because of the great diversity and dominance of ants in tropical environments, we used ant–plant interactions as a focal network to evaluate how variation in canopy cover influences patterns of plant–insect interactions in the Brazilian Amazon rain forest. We observed that small increases in canopy openness are associated with increased diversity of ant–plant interactions in our study area, and this change is independent of plant or ant species richness. Additionally, we found smaller niche overlap for both ants and plants associated with greater canopy openness. We hypothesize that enhanced light availability increases the breadth of ant foraging sources because variation in light availability gives rise to plant resources of different quality and amounts. Moreover, greater light availability promotes vegetative growth in plants, creating ant foraging ‘bridges’ between plants. In sum, our results highlight the importance of environmental heterogeneity as a determinant of ant–plant interaction diversity in tropical environments.     

The impact of environmental variability and species composition on the stability of experimental microbial populations and communities

Understanding how environmental fluctuations affect the stability of populations and communities is complex, for example, because direct effects of environmental variability on populations may be modified and propagated across communities by species interactions. One way to explore and further understand these complexities is via a factorial manipulation of community composition and environmental conditions. Using laboratory based aquatic microcosms we manipulated environmental fluctuation by creating two environments; one with variable light and one with constant light. Within these environments, community composition was manipulated by constructing communities from all possible combinations of three species that vary in their reliance on light for growth (an autotroph: a diatom completely reliant on light, a heterotroph: a Paramecium species not reliant on light, and a mixotroph: a Paramecium species somewhat reliant on light). Community composition was predicted to affect populations and communities by introducing and altering competitive interactions between species and affecting the degree of niche differentiation between species. We found that population stability was predominantly influenced by an interaction between community composition and environmental variability, whereby the effect of environmental variability synergistically combined with effects of community composition to reduce population stability. Covariance of populations was determined by an interaction between community composition and environmental variability, though this did not result from the effect of niche differentiation between species. Species interactions drove correlations between population biomass and the environment which otherwise did not exist. Our results demonstrate the complex and interrelated effects of abiotic and biotic factors on population and community stability, and suggest the need to consider aspects of community composition when predicting the impact of environmental fluctuations.     

Redefining the community: a species-based approach

We propose that effective community size can be defined on the basis of the web of indirect interactions experienced on average by each individual species. Indirect interaction chains are composed of links provided by direct interactions. We analyzed previously published data on 20 assemblages of species. Chain strengths were estimated by the weakest link and by the product of link strength. The average strength of the interaction chain decreased with increasing numbers of links with both models. Positive indirect interactions in chains with an even number of links offset negative direct interactions. We set the community size by the chain length where 95% of the indirect interactions are weaker than 10% of the mean of the absolute value of direct interaction strength. Using the multiplicative model, seven assemblages had a community size (web of interaction length) of three links, one of four links, and the remainder of communities were too small to set community size. The analysis suggests that effective communities of size are rarely investigated in ecological experiments.     

Decay of interspecific avian flock networks along a disturbance gradient in Amazonia

Our understanding of how anthropogenic habitat change shapes species interactions is in its infancy. This is in large part because analytical approaches such as network theory have only recently been applied to characterize complex community dynamics. Network models are a powerful tool for quantifying how ecological interactions are affected by habitat modification because they provide metrics that quantify community structure and function. Here, we examine how large-scale habitat alteration has affected ecological interactions among mixed-species flocking birds in Amazonian rainforest. These flocks provide a model system for investigating how habitat heterogeneity influences non-trophic interactions and the subsequent social structure of forest-dependent mixed-species bird flocks. We analyse 21 flock interaction networks throughout a mosaic of primary forest, fragments of varying sizes and secondary forest (SF) at the Biological Dynamics of Forest Fragments Project in central Amazonian Brazil. Habitat type had a strong effect on network structure at the levels of both species and flock. Frequency of associations among species, as summarized by weighted degree, declined with increasing levels of forest fragmentation and SF. At the flock level, clustering coefficients and overall attendance positively correlated with mean vegetation height, indicating a strong effect of habitat structure on flock cohesion and stability. Prior research has shown that trophic interactions are often resilient to large-scale changes in habitat structure because species are ecologically redundant. By contrast, our results suggest that behavioural interactions and the structure of non-trophic networks are highly sensitive to environmental change. Thus, a more nuanced, system-by-system approach may be needed when thinking about the resiliency of ecological networks.     

Genes under weaker stabilizing selection increase network evolvability and rapid regulatory adaptation to an environmental shift

Regulatory networks play a central role in the modulation of gene expression, the control of cellular differentiation, and the emergence of complex phenotypes. Regulatory networks could constrain or facilitate evolutionary adaptation in gene expression levels. Here, we model the adaptation of regulatory networks and gene expression levels to a shift in the environment that alters the optimal expression level of a single gene. Our analyses show signatures of natural selection on regulatory networks that both constrain and facilitate rapid evolution of gene expression level towards new optima. The analyses are interpreted from the standpoint of neutral expectations and illustrate the challenge to making inferences about network adaptation. Furthermore, we examine the consequence of variable stabilizing selection across genes on the strength and direction of interactions in regulatory networks and in their subsequent adaptation. We observe that directional selection on a highly constrained gene previously under strong stabilizing selection was more efficient when the gene was embedded within a network of partners under relaxed stabilizing selection pressure. The observation leads to the expectation that evolutionarily resilient regulatory networks will contain optimal ratios of genes whose expression is under weak and strong stabilizing selection. Altogether, our results suggest that the variable strengths of stabilizing selection across genes within regulatory networks might itself contribute to the long‐term adaptation of complex phenotypes.     

Non-neutral community dynamics: empirical predictions for ecosystem function and diversity from linearized consumer–resource interactions

A general model of linearized species interactions, essentially Lotka–Volterra theory, applied to questions of biodiversity has previously been shown to be a powerful tool for understanding local species–abundance patterns and community responses to environmental change for a single trophic level. Here this approach is extended to predict community composition and responses to environmental changes in trophically structured systems. We show how resource and consumer species richness and their relative abundances vary with the means and variances in enrichment level and strengths of intra- and interspecific interactions. Also demonstrated are the responses of local resource and consumer species richness to the global species pools at both trophic levels, as well as the covariation with net resource productivity. These predictions for resource and consumer specific responses to changes in environmental enrichment and global biodiversity are directly testable.     

Using experimental indices to quantify the strength of species interactions

Few methods for demonstrating the effects of species interactions rival that of the manipulative experiment ( Kareiva and Levin 2002 ). The now commonly performed removal or addition of predators, competitors, or mutualists to experimentally replicated populations of recipient species has irrefutably shown that species can have important effects on each other's populations, and that the strengths of these effects can vary considerably across environmental contexts and species identities. A large body of research is directed towards understanding this variation to forecast community dynamics and dissect how species interactions regulate community structure ( Chesson 2000 , McCann 2000 , Stachowicz 2001 , Duffy 2002 ). It has frequently been pointed out that progress in this field will require more explicit connections between ecological theory and the realities of nature ( Laska and Wootton 1998 , Parker et al. 1999 , Berlow et al. 2004 , Agrawal et al. 2007 ). This requires that interaction strengths measured experimentally be appropriate to the biology of our study systems and the mathematical abstractions we ascribe to them. Here we clarify some of the assumptions made in the application of the two most commonly used indices for measuring the strengths of species interactions with manipulative removal/addition experiments: Paine's index and the dynamic index. We explain how the values these indices are intended to measure – the per capita interaction strength between two species – are typically not estimated in common currencies. We then introduce extensions to these indices that alleviate a subset of previously made assumptions and limitations. These include a reformulation of Paine's index appropriate for an open‐recruitment system, and an extension of the dynamic index applicable to interactions of a particular nonlinear form. While we focus our language on predator–prey interactions, our discussions are pertinent to the measurement of species effects in other types of interactions as well ( Mitchell and Wass 1996 , Freckleton et al. 2009 ).     

Landscape diversity promotes stable food-web architectures in large rivers

Uncovering relationships between landscape diversity and species interactions is crucial for predicting how ongoing land-use change and homogenization will impact the stability and persistence of communities. However, such connections have rarely been quantified in nature. We coupled high-resolution river sonar imaging with annualized energetic food webs to quantify relationships among habitat diversity, energy flux, and trophic interaction strengths in large-river food-web modules that support the endangered Pallid Sturgeon. Our results demonstrate a clear relationship between habitat diversity and species interaction strengths, with more diverse foraging landscapes containing higher production of prey and a greater proportion of weak and potentially stabilizing interactions. Additionally, rare patches of large and relatively stable river sediments intensified these effects and further reduced interaction strengths by increasing prey diversity. Our findings highlight the importance of landscape characteristics in promoting stabilizing food-web architectures and provide direct relevance for future management of imperilled species in a simplified and rapidly changing world.     

Beyond connectedness: why pairwise metrics cannot capture community stability

The connectedness of species in a trophic web has long been a key structural characteristic for both theoreticians and empiricists in their understanding of community stability. In the past decades, there has been a shift from focussing on determining the number of interactions to taking into account their relative strengths. The question is: How do the strengths of the interactions determine the stability of a community? Recently, a metric has been proposed which compares the stability of observed communities in terms of the strength of three‐ and two‐link feedback loops (cycles of interaction strengths). However, it has also been suggested that we do not need to go beyond the pairwise structure of interactions to capture stability. Here, we directly compare the performance of the feedback and pairwise metrics. Using observed food‐web structures, we show that the pairwise metric does not work as a comparator of stability and is many orders of magnitude away from the actual stability values. We argue that metrics based on pairwise‐strength information cannot capture the complex organization of strong and weak links in a community, which is essential for system stability.     

Stabilizing mechanisms in a food web with an introduced omnivore

Intraguild predation (IGP) is an omnivorous food web configuration in which the top predator consumes both a competitor (consumer) and a second prey that it shares with the competitor. This omnivorous configuration occurs frequently in food webs, but theory suggests that it is unstable unless stabilizing mechanisms exist that can decrease the strength of the omnivore and consumer interaction. Although these mechanisms have been documented in native food webs, little is known about whether they operate in the context of an introduced species. Here, we study a marine mussel aquaculture system where the introduction of omnivorous mussels should generate an unstable food web that favors the extinction of the consumer, yet it persists. Using field and laboratory approaches, we searched for stabilizing mechanisms that could reduce interaction strengths in the food web. While field zooplankton counts suggested that mussels influence the composition and abundance of copepods, stable isotope results indicated that life‐history omnivory and cannibalism facilitated the availability of prey refugia, and reduced competition and the interaction strength between the mussel omnivore and zooplankton consumers. In laboratory experiments, however, we found no evidence of adaptive feeding which could weaken predator–consumer interactions. Our food web study suggests that the impact of an introduced omnivore may not only depend on its interaction with native species but also on the availability of stabilizing mechanisms that alter the strength of those interactions.     

Patterns in intraspecific interaction strengths and the stability of food webs

A common approach to analyse stability of biological communities is to calculate the interaction strength matrix. Problematic in this approach is defining intraspecific interaction strengths, represented by diagonal elements in the matrix, due to a lack of empirical data for these strengths. Theoretical studies have shown that an overall increase in these strengths enhances stability. However, the way in which the pattern in intraspecific interaction strengths, i.e. the variation in these strengths between species, influences stability has received little attention. We constructed interaction strength matrices for 11 real soil food webs in which four patterns for intraspecific interaction strengths were chosen, based on the ecological literature. These patterns included strengths that were (1) similar for all species, (2) trophic level dependent, (3) biomass dependent, or (4) death rate dependent. These four patterns were analysed for their influence on (1) ranking food webs by their stability and (2) the response in stability to variation of single interspecific interaction strengths. The first analysis showed that ranking the 11 food webs by their stability was not strongly influenced by the choice of diagonal pattern. In contrast, the second analysis showed that the response of food web stability to variation in single interspecific interaction strengths was sensitive to the choice of diagonal pattern. Notably, stability could increase using one pattern and decrease using another. This result asks for deliberate approaches to choose diagonal element values in order to make predictions on how particular species, interactions, or other food web parameters affect food web stability.