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.     

Testing the link between species interactions and species co‐occurrence in a trophic network

Species interactions are dynamic processes that vary across environmental and ecological contexts, and operate across scale boundaries, making them difficult to quantify. Nevertheless, ecologists are increasingly interested in inferring species interactions from observational data using statistical analyses of their spatial co‐occurrence patterns. Trophic interactions present a particular challenge, as predators and prey may frequently or rarely co‐occur, depending on the spatial or temporal scale of observation. In this study, we investigate the accuracy of inferred interactions among species that both compete and trophically interact. We utilized a long‐term dataset of pond‐breeding amphibian co‐occurrences from Mt Rainier National Park (Washington, USA) and compiled a new dataset of their empirical interactions from the literature. We compared the accuracy of four statistical methods in inferring these known species interactions from spatial associations. We then used the best performing statistical method, the Markov network, to further investigate the sensitivity of interaction inference to spatial scale‐dependence and the presence of predators. We show that co‐occurrence methods are generally inaccurate when estimating trophic interactions. Further the strength and sign of inferred interactions were dependent upon the spatial scale of observation and predator presence influenced the detectability of competitive interactions among prey species. However, co‐occurrence analysis revealed new patterns of spatial association among pairs of species with known interactions. Overall, our study highlights a limiting frontier in co‐occurrence theory and the disconnect between widely implemented methodologies and their ability to accurately infer interactions in trophically‐structured communities.     

Species co-occurrence: the case of congeneric species and a causal approach to patterns of species association

Aim To test whether congeneric species are significantly associated with one another in space, either positively or negatively. Also, to provide a framework for a causal investigation of co‐occurrence patterns by a parallel comparison of interactions in geographical and ecological data matrices. Location For the analysis of congeneric species’ co‐occurrences we used 30 matrices covering a wide range of taxa and geographical areas, while for the causal investigation we used the distribution of 50 terrestrial isopod species on 20 islands and 264 sampling stations in the central Aegean archipelago, as well as a number of ecological variables for each sampling station. Methods We developed a software program (cooc ) that incorporates the species‐by‐species approach to co‐occurrence analysis using EcoSim's output of prior null model analysis of co‐occurrence. We describe this program in detail, and use it to investigate one of the most common assembly rules, namely, the decreased levels of co‐occurrence among congeneric species pairs. For the causal analysis, we proceed likewise, cross‐checking the results from the geographical and the ecological matrices. There is only one possible combination of results that can support claims for direct competition among species. Results We do not get any strong evidence for widespread competition among congeneric species, while most communities investigated do not show significant patterns of species associations. The causal analysis suggests that the principal factors behind terrestrial isopod species associations are of historical nature. Some exceptional cases are also discussed. Main conclusions Presence/absence data for a variety of taxa do not support the assembly rule that congeneric species are under more intense competition compared to less related species. Also, these same data do not suggest strong interactions among species pairs, regardless of taxonomic status. When significant species associations can be seen in such matrices, they mainly reflect the effects of history or of habitat requirements.     

Accounting for environmental variation in co‐occurrence modelling reveals the importance of positive interactions in root‐associated fungal communities

Understanding the role of interspecific interactions in shaping ecological communities is one of the central goals in community ecology. In fungal communities, measuring interspecific interactions directly is challenging because these communities are composed of large numbers of species, many of which are unculturable. An indirect way of assessing the role of interspecific interactions in determining community structure is to identify the species co‐occurrences that are not constrained by environmental conditions. In this study, we investigated co‐occurrences among root‐associated fungi, asking whether fungi co‐occur more or less strongly than expected based on the environmental conditions and the host plant species examined. We generated molecular data on root‐associated fungi of five plant species evenly sampled along an elevational gradient at a high arctic site. We analysed the data using a joint species distribution modelling approach that allowed us to identify those co‐occurrences that could be explained by the environmental conditions and the host plant species, as well as those co‐occurrences that remained unexplained and thus more probably reflect interactive associations. Our results indicate that not only negative but also positive interactions play an important role in shaping microbial communities in arctic plant roots. In particular, we found that mycorrhizal fungi are especially prone to positively co‐occur with other fungal species. Our results bring new understanding to the structure of arctic interaction networks by suggesting that interactions among root‐associated fungi are predominantly positive.     

A new integrative framework for large‐scale assessments of biodiversity and community dynamics,using littoral gastropods and crabs of British Columbia,Canada

Improving our understanding of species responses to environmental changes is an important contribution ecologists can make to facilitate effective management decisions. Novel synthetic approaches to assessing biodiversity and ecosystem integrity are needed, ideally including all species living in a community and the dynamics defining their ecological relationships. Here, we present and apply an integrative approach that links high‐throughput, multicharacter taxonomy with community ecology. The overall purpose is to enable the coupling of biodiversity assessments with investigations into the nature of ecological interactions in a community‐level data set. We collected 1195 gastropods and crabs in British Columbia. First, the General mixed Yule‐coalescent (GMYC) and the Poisson Tree Processes (PTP) methods for proposing primary species‐hypotheses based on cox1 sequences were evaluated against an integrative taxonomic framework. We then used data on the geographic distribution of delineated species to test species co‐occurrence patterns for nonrandomness using community‐wide and pairwise approaches. Results showed that PTP generally outperformed GMYC and thus constitutes a more effective option for producing species‐hypotheses in community‐level data sets. Nonrandom species co‐occurrence patterns indicative of ecological relationships or habitat preferences were observed for grazer gastropods, whereas assemblages of carnivorous gastropods and crabs appeared influenced by random processes. Species‐pair associations were consistent with current ecological knowledge, thus suggesting that applying community assembly within a large taxonomical framework constitutes a valuable tool for assessing ecological interactions. Combining phylogenetic, morphological and co‐occurrence data enabled an integrated view of communities, providing both a conceptual and pragmatic framework for biodiversity assessments and investigations into community dynamics.     

Co‐occurrence is not evidence of ecological interactions

There is a rich amount of information in co‐occurrence (presence–absence) data that could be used to understand community assembly. This proposition first envisioned by Forbes (1907) and then Diamond (1975) prompted the development of numerous modelling approaches (e.g. null model analysis, co‐occurrence networks and, more recently, joint species distribution models). Both theory and experimental evidence support the idea that ecological interactions may affect co‐occurrence, but it remains unclear to what extent the signal of interaction can be captured in observational data. It is now time to step back from the statistical developments and critically assess whether co‐occurrence data are really a proxy for ecological interactions. In this paper, we present a series of arguments based on probability, sampling, food web and coexistence theories supporting that significant spatial associations between species (or lack thereof) is a poor proxy for ecological interactions. We discuss appropriate interpretations of co‐occurrence, along with potential avenues to extract as much information as possible from such data.     

The phylogenetic structure of plant–pollinator networks increases with habitat size and isolation

Similarity among species in traits related to ecological interactions is frequently associated with common ancestry. Thus, closely related species usually interact with ecologically similar partners, which can be reinforced by diverse co‐evolutionary processes. The effect of habitat fragmentation on the phylogenetic signal in interspecific interactions and correspondence between plant and animal phylogenies is, however, unknown. Here, we address to what extent phylogenetic signal and co‐phylogenetic congruence of plant–animal interactions depend on habitat size and isolation by analysing the phylogenetic structure of 12 pollination webs from isolated Pampean hills. Phylogenetic signal in interspecific interactions differed among webs, being stronger for flower‐visiting insects than plants. Phylogenetic signal and overall co‐phylogenetic congruence increased independently with hill size and isolation. We propose that habitat fragmentation would erode the phylogenetic structure of interaction webs. A decrease in phylogenetic signal and co‐phylogenetic correspondence in plant–pollinator interactions could be associated with less reliable mutualism and erratic co‐evolutionary change.     

Invasions cause biodiversity loss and community simplification in vertebrate food webs

Global change is increasing the occurrence of perturbation events on natural communities, with biological invasions posing a major threat to ecosystem integrity and functioning worldwide. Most studies addressing biological invasions have focused on individual species or taxonomic groups to understand both, the factors determining invasion success and their effects on native species. A more holistic approach that considers multispecies communities and species’ interactions can contribute to a better understanding of invasion effects on complex communities. Here we address biological invasions on species‐rich food webs. We performed in silico experiments on empirical vertebrate food webs by introducing virtual species characterised by different ecological roles and belonging to different trophic groups. We varied a number of invasive species traits, including their diet breadth, the number of predators attacking them, and the bioenergetic thresholds below which invader and native species become extinct. We found that simpler food webs were more vulnerable to invasions, and that relatively less connected mammals were the most successful invaders. Invasions altered food web structure by decreasing species richness and the number of links per species, with most extinctions affecting poorly connected birds. Our food web approach allows identifying the combinations of trophic factors that facilitate or prevent biological invasions, and it provides testable predictions on the effects of invasions on the structure and dynamics of multitrophic communities.     

Ecological correlates of the spatial co‐occurrence of sympatric mammalian carnivores worldwide

The composition of local mammalian carnivore communities has far‐reaching effects on terrestrial ecosystems worldwide. To better understand how carnivore communities are structured, we analysed camera trap data for 108 087 trap days across 12 countries spanning five continents. We estimate local probabilities of co‐occurrence among 768 species pairs from the order Carnivora and evaluate how shared ecological traits correlate with probabilities of co‐occurrence. Within individual study areas, species pairs co‐occurred more frequently than expected at random. Co‐occurrence probabilities were greatest for species pairs that shared ecological traits including similar body size, temporal activity pattern and diet. However, co‐occurrence decreased as compared to other species pairs when the pair included a large‐bodied carnivore. Our results suggest that a combination of shared traits and top‐down regulation by large carnivores shape local carnivore communities globally.     

Rewiring and indirect effects underpin modularity reshuffling in a marine food web under environmental shifts

Species are characterized by physiological and behavioral plasticity, which is part of their response to environmental shifts. Nonetheless, the collective response of ecological communities to environmental shifts cannot be predicted from the simple sum of individual species responses, since co‐existing species are deeply entangled in interaction networks, such as food webs. For these reasons, the relation between environmental forcing and the structure of food webs is an open problem in ecology. To this respect, one of the main problems in community ecology is defining the role each species plays in shaping community structure, such as by promoting the subdivision of food webs in modules—that is, aggregates composed of species that more frequently interact—which are reported as community stabilizers. In this study, we investigated the relationship between species roles and network modularity under environmental shifts in a highly resolved food web, that is, a “weighted” ecological network reproducing carbon flows among marine planktonic species. Measuring network properties and estimating weighted modularity, we show that species have distinct roles, which differentially affect modularity and mediate structural modifications, such as modules reconfiguration, induced by environmental shifts. Specifically, short‐term environmental changes impact the abundance of planktonic primary producers; this affects their consumers’ behavior and cascades into the overall rearrangement of trophic links. Food web re‐adjustments are both direct, through the rewiring of trophic‐interaction networks, and indirect, with the reconfiguration of trophic cascades. Through such “systemic behavior,” that is, the way the food web acts as a whole, defined by the interactions among its parts, the planktonic food web undergoes a substantial rewiring while keeping almost the same global flow to upper trophic levels, and energetic hierarchy is maintained despite environmental shifts. This behavior suggests the potentially high resilience of plankton networks, such as food webs, to dramatic environmental changes, such as those provoked by global change.     

Competitive interactions change the pattern of species co‐occurrences under neutral dispersal

Non‐random patterns of species segregation and aggregation within ecological communities are often interpreted as evidence for interspecific interactions. However, it is unclear whether theoretical models can predict such patterns and how environmental factors may modify the effects of species interactions on species co‐occurrence. Here we extend a spatially explicit neutral model by including competitive effects on birth and death probabilities to assess whether competition alone is able to produce non‐random patterns of species co‐occurrence. We show that transitive and intransitive competitive hierarchies alone (in the absence of environmental heterogeneity) are indeed able to generate non‐random patterns with commonly used metrics and null models. Moreover, even weak levels of intransitive competition can increase local species richness. However, there is no simple rule or consistent directional change towards aggregation or segregation caused by competitive interactions. Instead, the spatial pattern depends on both the type of species interaction and the strength of dispersal. We conclude that co‐occurrence analysis alone may not able to identify the underlying processes that generate the patterns.     

A comprehensive framework for the study of species co‐occurrences,nestedness and turnover

Binary presence–absence matrices (rows = species, columns = sites) are often used to quantify patterns of species co‐occurrence, and to infer possible biotic interactions from these patterns. Previous classifications of co‐occurrence patterns as nested, segregated, or modular have led to contradictory results and conclusions. These analyses usually do not incorporate the functional traits of the species or the environmental characteristics of the sites, even though the outcomes of species interactions often depend on trait expression and site quality. Here we address this shortcoming by developing a method that incorporates realized functional and environmental niches, and relates them to species co‐occurrence patterns. These niches are defined from n‐dimensional ellipsoids, and calculated from the n eigenvectors and eigenvalues of the variance–covariance matrix of measured environmental or trait variables. Average niche overlap among species and the spatial distribution of niches define a triangle plot with vertices of species segregation (low niche overlap), nestedness (high niche overlap), and modular co‐occurrence (clusters of overlapping niches). Applying this framework to temperate understorey plant communities in southwest Poland, we found a consistent modular structure of species occurrences, a pattern not detected by conventional presence–absence analysis. These results suggest that, in our case study, habitat filtering is the most important process structuring understorey plant communities. Furthermore, they demonstrate how incorporating trait and environmental data into co‐occurrence analysis improves pattern detection and provides a stronger theoretical framework for understanding community structure.     

The pairwise approach to analysing species co‐occurrence

The analysis of species co‐occurrence patterns continues to be a main pursuit of ecologists, primarily because the coexistence of species is fundamentally important in evaluating various theories, principles and concepts. Examples include community assembly, equilibrium versus non‐equilibrium organization of communities, resource partitioning and ecological character displacement, the local–regional species diversity relationship, and the metacommunity concept. Traditionally, co‐occurrence has been measured and tested at the level of an entire species presence–absence matrix wherein various algorithms are used to randomize matrices and produce statistical null distributions of metrics that quantify structure in the matrix. This approach implicitly recognizes a presence–absence matrix as having some real ecological identity (e.g. a set of species exhibiting nestedness among a set of islands) in addition to being a unit of statistical analysis. An emerging alternative is to test for non‐random co‐occurrence between paired species. The pairwise approach does not analyse matrix‐level structure and thus views a species pair as the fundamental unit of co‐occurrence. Inferring process from pattern is very difficult in analyses of co‐occurrence; however, the pairwise approach may make this task easier by simplifying the analysis and resulting inferences to associations between paired species.     

Patterns of co‐occurrences in a killifish metacommunity are more related with body size than with species identity

Body size may be more important than species identity in determining species interactions and community structure. However, co‐occurrence of organisms has commonly been analysed from a taxonomic perspective and the body size is rarely taken into account. On six sampling occasions, we analysed patterns of killifish co‐occurrences in nestedness (tendency for less rich communities to be subsamples of the richest), checkerboard structure (tendency for species segregation), and modularity (tendency for groups to co‐occur more frequently than random expectation) in a pond metacommunity located in Uruguay. We contrasted co‐occurrence patterns among species and body size‐classes (individuals from different species were combined into size categories). The analysis was performed at two spatial scales: ponds (communities) and sample units within ponds. Observed nestedness was frequently smaller than the null expectation, with significantly greater deviations for body size‐classes than for species, and for sample units than for communities. At the sample unit level, individuals tended to segregate (i.e. clump into a checkerboard pattern) to a larger extent by body size rather than by taxonomy. Modularity was rarely detected, but nevertheless indicated a level of taxonomic organization not evident in nestedness or checkerboard indices. Identification of the spatial scale and organization at which ecological forces determine community structure is a basic requirement for advancement of robust theory. In our study system, these ecological forces probably structured the community by body sizes of interacting organisms rather than by species identities.     

Ecological traits influence the phylogenetic structure of bird species co‐occurrences worldwide

The extent to which species’ ecological and phylogenetic relatedness shape their co‐occurrence patterns at large spatial scales remains poorly understood. By quantifying phylogenetic assemblage structure within geographic ranges of >8000 bird species, we show that global co‐occurrence patterns are linked – after accounting for regional effects – to key ecological traits reflecting diet, mobility, body size and climatic preference. We found that co‐occurrences of carnivorous, migratory and cold‐climate species are phylogenetically clustered, whereas nectarivores, herbivores, frugivores and invertebrate eaters tend to be more phylogenetically overdispersed. Preference for open or forested habitats appeared to be independent from the level of phylogenetic clustering. Our results advocate for an extension of the tropical niche conservatism hypothesis to incorporate ecological and life‐history traits beyond the climatic niche. They further offer a novel species‐oriented perspective on how biogeographic and evolutionary legacies interact with ecological traits to shape global patterns of species coexistence in birds.     

A conceptual framework for studying the strength of plant–animal mutualistic interactions

The strength of species interactions influences strongly the structure and dynamics of ecological systems. Thus, quantifying such strength is crucial to understand how species interactions shape communities and ecosystems. Although the concepts and measurement of interaction strength in food webs have received much attention, there has been comparatively little progress in the context of mutualism. We propose a conceptual scheme for studying the strength of plant–animal mutualistic interactions. We first review the interaction strength concepts developed for food webs, and explore how these concepts have been applied to mutualistic interactions. We then outline and explain a conceptual framework for defining ecological effects in plant–animal mutualisms. We give recommendations for measuring interaction strength from data collected in field studies based on a proposed approach for the assessment of interaction strength in plant–animal mutualisms. This approach is conceptually integrative and methodologically feasible, as it focuses on two key variables usually measured in field studies: the frequency of interactions and the fitness components influenced by the interactions.     

Phylogeny and co‐occurrence of mammal species on Southeast Asian islands

Aim Islands have often been used as model systems in community ecology. The incorporation of information on phylogenetic relatedness of species in studies of island assemblage structure is still uncommon, but could provide valuable insights into the processes of island community assembly. We propose six models of island community assembly that make different predictions about the associations between co‐occurrences of species pairs on islands, phylogenetic relatedness and ecological similarity. We then test these models using data on mammals of Southeast Asian islands. Location Two hundred and forty islands of the Sundaland region of Southeast Asia. Methods We quantified the co‐occurrence of species pairs on islands, and identified pairs that co‐occur more frequently (positive co‐occurrence) or less frequently (negative co‐occurrence) than expected under null models. We then examined the distributions of these significantly deviating pairs with respect to phylogenetic relatedness and ecological differentiation, and compared these patterns with those predicted by the six community assembly models. We used permutation regression to test whether co‐occurrence patterns are predicted by relatedness, body size difference or difference in diet quality. Separate co‐occurrence matrices were analysed in this way for seven mammal families and four smaller subsets of the islands of Sundaland. Results In many matrices, average numbers of negative co‐occurrences were higher than expected under null models. This is consistent with assemblage structuring by competition, but may also result from low geographic overlap of species pairs, which contributes to negative co‐occurrences at the archipelago‐wide level. Distributions of species pairs within plots of phylogenetic distance × ecological differentiation were consistent with competition, habitat filtering or within‐island speciation models, depending on the taxon. Regressions indicated that co‐occurrence was more likely among closely related species pairs within the Viverridae and Sciuridae, but in most matrices phylogenetic distance was unrelated to co‐occurrence. Main conclusions Simple deterministic models linking co‐occurrence with phylogeny and ecology are a useful framework for interpreting distributions and assemblage structure of island species. However, island assemblages in Sundaland have probably been shaped by a complex idiosyncratic set of interacting ecological and evolutionary processes, limiting the predictive power of such models.     

Ant mosaics occur in SE Asian oil palm plantation but not rain forest and are influenced by the presence of nest‐sites and non‐native species

Interaction networks within biotic communities can be dramatically altered by anthropogenic habitat modification. Ants, an important ecological group, often interact competitively to form mosaic‐like patterns in disturbed plantation habitats, in which dominant species form mutually exclusive territories. However, the existence of these ant mosaics in pristine forests is contentious. Here we assess the relative strengths of ant competitive interactions in oil palm plantation and primary rain forest in Sabah, Malaysia, using null models of species co‐occurrence. We use two metrics: the C‐score, which measures mean degree of overall co‐occurrence, and a novel metric, the C_var‐score, which measures the variance in degree of co‐occurrence. We also investigate the role of nest sites by collecting ants from canopy and leaf litter microhabitats, and from epiphytic ferns, an important nest site for canopy ants. Furthermore, we assess whether non‐native species, which were widespread in oil palm plantation (61 occurrences vs five in rain forest) are important in driving the formation of ant mosaics. We found no evidence for ant mosaics in any primary forest microhabitat. In oil palm plantation, segregation between species was pronounced in epiphytes, weak in the rest of the canopy and absent in leaf litter communities. Intriguingly, exclusion of non‐native ant species from analyses increased the degree of negative species co‐occurrence in all three microhabitats, with species segregation in the oil palm canopy becoming statistically significant. Our results suggest that invasion of plantation habitats by non‐native species does not drive increased species segregation in ant communities. Rather, high degrees of species segregation might relate to changes in the importance of canopy nest sites, with colonies competing more strongly for these in plantations. In primary forests, weaker nest‐site limitation and the highly complex, more vertically stratified, non‐uniform canopy could lead to random co‐occurrence between ant species at the scales studied here.     

Estimating co-extinction threats in terrestrial ecosystems

The biosphere is changing rapidly due to human endeavour. Because ecological communities underlie networks of interacting species, changes that directly affect some species can have indirect effects on others. Accurate tools to predict these direct and indirect effects are therefore required to guide conservation strategies. However, most extinction-risk studies only consider the direct effects of global change—such as predicting which species will breach their thermal limits under different warming scenarios—with predictions of trophic cascades and co-extinction risks remaining mostly speculative. To predict the potential indirect effects of primary extinctions, data describing community interactions and network modelling can estimate how extinctions cascade through communities. While theoretical studies have demonstrated the usefulness of models in predicting how communities react to threats like climate change, few have applied such methods to real-world communities. This gap partly reflects challenges in constructing trophic network models of real-world food webs, highlighting the need to develop approaches for quantifying co-extinction risk more accurately. We propose a framework for constructing ecological network models representing real-world food webs in terrestrial ecosystems and subjecting these models to co-extinction scenarios triggered by probable future environmental perturbations. Adopting our framework will improve estimates of how environmental perturbations affect whole ecological communities. Identifying species at risk of co-extinction (or those that might trigger co-extinctions) will also guide conservation interventions aiming to reduce the probability of co-extinction cascades and additional species losses.     

Network structural properties mediate the stability of mutualistic communities

Key advances are being made on the structures of predator–prey food webs and competitive communities that enhance their stability, but little attention has been given to such complexity–stability relationships for mutualistic communities. We show, by way of theoretical analyses with empirically informed parameters, that structural properties can alter the stability of mutualistic communities characterized by nonlinear functional responses among the interacting species. Specifically, community resilience is enhanced by increasing community size (species diversity) and the number of species interactions (connectivity), and through strong, symmetric interaction strengths of highly nested networks. As a result, mutualistic communities show largely positive complexity–stability relationships, in opposition to the standard paradox. Thus, contrary to the commonly-held belief that mutualism's positive feedback destabilizes food webs, our results suggest that interplay between the structure and function of ecological networks in general, and consideration of mutualistic interactions in particular, may be key to understanding complexity–stability relationships of biological communities as a whole.