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.     

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 meta‐analysis of nestedness and turnover components of beta diversity across organisms and ecosystems

Aim

The number of studies investigating the nestedness and turnover components of beta diversity has increased substantially, but our general understanding of the drivers of turnover and nestedness remains elusive. Here, we examined the effects of species traits, spatial extent, latitude and ecosystem type on the nestedness and turnover components of beta diversity.

Location

Global.

Time period

1968–2017.

Major taxa studied

From bacteria to mammals.

Methods

From the 99 studies that partition total beta diversity into its turnover and nestedness components, we assembled 269 and 259 data points for the pairwise and multiple site beta‐diversity metrics, respectively. Our data covered a broad variation in species dispersal type, body size and trophic position. The data were from freshwater, marine and terrestrial realms, and encompassed geographical areas from the tropics to near polar regions. We used linear modelling as a meta‐regression tool to analyse the data.

Results

Pairwise turnover, multiple site turnover and total beta diversity all decreased significantly with latitude. In contrast, multiple site nestedness showed a positive relationship with latitude. Beta‐diversity components did not generally differ among the realms. The turnover component and total beta diversity increased with spatial extent, whereas nestedness was scale invariant for pairwise metrics. Multiple site beta‐diversity components did not vary with spatial extent. Surprisingly, passively dispersed organisms had lower turnover and total beta diversity than flying organisms. Body size showed a relatively weak relationship with beta diversity but had important interactions with trophic position, thus also affecting beta diversity via interactive effects. Producers had significantly higher average pairwise turnover and total beta diversity than carnivores.

Main conclusions

The present results provide evidence that species turnover, being consistently the larger component of total beta diversity, and nestedness are related to the latitude of the study area and intrinsic organismal features. We showed that two beta‐diversity components had generally opposing patterns with regard to latitude. We highlight that beta‐diversity partition may give additional insights into the underlying causes of spatial variability in biotic communities compared with total beta diversity alone.     

A multiple-site dissimilarity measure for species presence/absence data and its relationship with nestedness and turnover

Multiple-site dissimilarity may be caused by two opposite processes of meta-community organization, such as species nestedness and turnover. Therefore, discriminating among these contributions is necessary for linking multiple-site dissimilarity to ecosystem functioning. This paper introduces a measure of multiple-site dissimilarity or beta diversity for presence/absence data that is based on information on species absences from the species × sites matrix. It is also shown that the newly proposed dissimilarity index can be additively partitioned into species nestedness and turnover.     

Limited niche differentiation within remarkable co‐occurrences of congeneric species: Monomorium ants in the Australian seasonal tropics

Niche theory predicts that few closely related species can co‐occur because such species tend to be ecologically similar and niche differentiation is required to avoid competitive exclusion. We analyse the co‐occurrence of a remarkable 10–15 species of the ant genus Monomorium occurring within single 10 × 10 m plots in a tropical savanna of northern Australia. Most of the species are undescribed, so we use genetic analysis to validate our species demarcations. We document nest dispersion patterns, and investigate differentiation in the three primary niche dimensions: space, time and food. We also examine species differences in competitive abilities, by describing rates of foraging activity, foraging ranges, worker aggression, and levels of behavioural dominance. Analyses of nest and forager distributions showed very limited evidence of spatial segregation within plots. The great majority of species foraged either exclusively or primarily during daylight hours. Body size and isotopic analyses indicated very limited dietary differentiation. Such limited niche partitioning occurred despite the species differing markedly in their competitive abilities as measured by rates of resource discovery, recruitment and monopolization. Our findings defy the traditional assumption that multiple closely related and ecologically similar species of highly interactive taxa cannot co‐occur. It seems very likely that species coexistence in our study system is determined to a very large degree by stochastic processes relating to dispersal and establishment, as predicted by neutral theory. However, neutral theory assumes competitive equivalence, whereas we found very marked differences in the competitive abilities of our co‐occurring species. We suggest that competitive exclusion is prevented by the modular nature of ant colonies, with competition limiting colony performance but not preventing occurrence. We conclude that other factors that allow species persistence, and not just competitive equivalence, can allow dispersal and establishment processes to drive species coexistence.     

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.     

The relationship between species replacement,dissimilarity derived from nestedness,and nestedness

Aim Beta diversity can be partitioned into two components: dissimilarity due to species replacement and dissimilarity due to nestedness ( Baselga, 2010 , Global Ecology and Biogeography, 19 , 134–143). Several contributions have challenged this approach or proposed alternative frameworks. Here, I review the concepts and methods used in these recent contributions, with the aim of clarifying: (1) the rationale behind the partitioning of beta diversity into species replacement and nestedness‐resultant dissimilarity, (2) how, based on this rationale, numerators and denominators of indices have to match, and (3) how nestedness and nestedness‐resultant dissimilarity are related but different concepts. Innovation The rationale behind measures of species replacement (turnover) dictates that the number of species that are replaced between sites (numerator of the index) has to be relativized with respect to the total number of species that could potentially be replaced (denominator). However, a recently proposed partition of Jaccard dissimilarity fails to do this. In consequence, this partition underestimates the contribution of species replacement and overestimates the contribution of richness differences to total dissimilarity. I show how Jaccard dissimilarity can be partitioned into meaningful turnover and nestedness components, and extend these new indices to multiple‐site situations. Finally the concepts of nestedness and nestedness‐resultant dissimilarity are discussed. Main conclusions Nestedness should be assessed using consistent measures that depend both on paired overlap and matrix filling, e.g. NODF, whereas beta‐diversity patterns should be examined using measures that allow the total dissimilarity to be separated into the components of dissimilarity due to species replacement and dissimilarity due to nestedness. In the case of multiple‐site dissimilarity patterns, averaged pairwise indices should never be used because the mean of the pairwise values is unable to accurately reflect the multiple‐site attributes of dissimilarity.     

Partitioning the turnover and nestedness components of beta diversity

Aim  Beta diversity (variation of the species composition of assemblages) may reflect two different phenomena, spatial species turnover and nestedness of assemblages, which result from two antithetic processes, namely species replacement and species loss, respectively. The aim of this paper is to provide a unified framework for the assessment of beta diversity, disentangling the contribution of spatial turnover and nestedness to beta-diversity patterns.
Innovation  I derive an additive partitioning of beta diversity that provides the two separate components of spatial turnover and nestedness underlying the total amount of beta diversity. I propose two families of measures of beta diversity for pairwise and multiple-site situations. Each family comprises one measure accounting for all aspects of beta diversity, which is additively decomposed into two measures accounting for the pure spatial turnover and nestedness components, respectively. Finally, I provide a case study using European longhorn beetles to exemplify the relevance of disentangling spatial turnover and nestedness patterns.
Main conclusion  Assigning the different beta-diversity patterns to their respective biological phenomena is essential for analysing the causality of the processes underlying biodiversity. Thus, the differentiation of the spatial turnover and nestedness components of beta diversity is crucial for our understanding of central biogeographic, ecological and conservation issues.     

Linking microbial co‐occurrences to soil ecological processes across a woodland‐grassland ecotone

Ecotones between distinct ecosystems have been the focus of many studies as they offer valuable insights into key drivers of community structure and ecological processes that underpin function. While previous studies have examined a wide range of above‐ground parameters in ecotones, soil microbial communities have received little attention. Here we investigated spatial patterns, composition, and co‐occurrences of archaea, bacteria, and fungi, and their relationships with soil ecological processes across a woodland‐grassland ecotone. Geostatistical kriging and network analysis revealed that the community structure and spatial patterns of soil microbiota varied considerably between three habitat components across the ecotone. Woodland samples had significantly higher diversity of archaea while the grassland samples had significantly higher diversity of bacteria. Microbial co‐occurrences reflected differences in soil properties and ecological processes. While microbial networks were dominated by bacterial nodes, different ecological processes were linked to specific microbial guilds. For example, soil phosphorus and phosphatase activity formed the largest clusters in their respective networks, and two lignolytic enzymes formed joined clusters. Bacterial ammonia oxidizers were dominant over archaeal oxidizers and showed a significant association (p < 0.001) with potential nitrification (PNR), with the PNR subnetwork being dominated by Betaproteobacteria. The top ten keystone taxa comprised six bacterial and four fungal OTUs, with Random Forest Analysis revealing soil carbon and nitrogen as the determinants of the abundance of keystone taxa. Our results highlight the importance of assessing interkingdom associations in soil microbial networks. Overall, this study shows how ecotones can be used as a model to delineate microbial structural patterns and ecological processes across adjoining land‐uses within a landscape.     

A network approach for inferring species associations from co‐occurrence data

Positive and negative associations between species are a key outcome of community assembly from regional species pools. These associations are difficult to detect and can be caused by a range of processes such as species interactions, local environmental constraints and dispersal. We integrate new ideas around species distribution modeling, covariance matrix estimation, and network analysis to provide an approach to inferring non‐random species associations from local‐ and regional‐scale occurrence data. Specifically, we provide a novel framework for identifying species associations that overcomes three challenges: 1) correcting for indirect effects from other species, 2) avoiding spurious associations driven by regional‐scale distributions, and 3) describing these associations in a multi‐species context. We highlight a range of research questions and analyses that this framework is able to address. We show that the approach is statistically robust using simulated data. In addition, we present an empirical analysis of > 1000 North American tree communities that gives evidence for weak positive associations among small groups of species. Finally, we discuss several possible extensions for identifying drivers of associations, predicting community assembly, and better linking biogeography and community ecology.     

Fifty years of co‐evolution and beyond: integrating co‐evolution from molecules to species

Fifty years after Ehrlich and Raven's seminal paper, the idea of co‐evolution continues to grow as a key concept in our understanding of organic evolution. This concept has not only provided a compelling synthesis between evolutionary biology and community ecology, but has also inspired research that extends beyond its original scope. In this article, we identify unresolved questions about the co‐evolutionary process and advocate for the integration of co‐evolutionary research from molecular to interspecific interactions. We address two basic questions: (i) What is co‐evolution and how common is it? (ii) What is the unit of co‐evolution? Both questions aim to explore the heart of the co‐evolutionary process. Despite the claim that co‐evolution is ubiquitous, we argue that there is in fact little evidence to support the view that reciprocal natural selection and coadaptation are common in nature. We also challenge the traditional view that co‐evolution only occurs between traits of interacting species. Co‐evolution has the potential to explain evolutionary processes and patterns that result from intra‐ and intermolecular biochemical interactions within cells, intergenomic interactions (e.g. nuclear‐cytoplasmic) within species, as well as intergenomic interactions mediated by phenotypic traits between species. Research that bridges across these levels of organization will help to advance our understanding of the importance of the co‐evolutionary processes in shaping the diversity of life on Earth.     

Multisite and multispecies measures of overlap,co‐occurrence,and co‐diversity

Several indices measure the association or segregation between two species and the similarity or differentiation between two sets of species. These indices are based on the overlap in the distribution of species (measured with the number of co‐occurrences) or on the overlap in species composition of sites (measured with the number of species that are shared between two sites). This paper shows that when evaluating more than two species the number of overlaps and the number of pairwise co‐occurrences are not equal, as it is the case for two species. Equivalently, when comparing more than two species assemblages, the number of overlaps differ from the number of instances of species sharing by pairs of sites (the ‘co‐diversities’). Considering this distinction, two different types of multispecies and multisite indices can be derived: indices of general overlap and indices of co‐occurrence or co‐diversity. Here I present a complete series of the two types of indices that correspond to the popular Jaccard, Sørensen, and Simpson two‐species or two‐site indices. Indices of general overlap are defined by three parameters (the total number of species, the total number of sites, and the total number of occurrences), whereas indices of co‐occurrence or co‐diversity depend on those parameters plus an additional one that is defined by the values of species richness or range size. Consequently, the two types of indices respond differently to null models, depending on the parameters that are fixed or randomized. All indices correlate well with the mean of the traditional indices calculated pair by pair, and the correspondence is extremely close for the new indices of co‐occurrence and co‐diversity. These properties should be useful in clarifying some of the confusion that exists in the current discussion about the advantages and disadvantages of pairwise vs community‐wide approaches in the analysis of diversity.     

Effects of dispersal mode on the environmental and spatial correlates of nestedness and species turnover in pond communities

Advances in metacommunity theory have made a significant contribution to understanding the drivers of variation in biological communities. However, there has been limited empirical research exploring the expression of metacommunity theory for two fundamental components of beta diversity: nestedness and species turnover. In this paper, we examine the influence of local environmental and a range of spatial variables (hydrological connectivity, proximity and overall spatial structure) on total beta diversity and the nestedness and turnover components of beta diversity for the entire macroinvertebrate community and active and passively dispersing taxa within pond habitats. High beta diversity almost entirely reflects patterns of species turnover (replacement) rather than nestedness (differences in species richness) in our dataset. Local environmental variables were the main drivers of total beta diversity, nestedness and turnover when the entire community was considered and for both active and passively dispersing taxa. The influence of spatial processes on passively dispersing taxa, total beta diversity and nestedness was significantly greater than for actively dispersing taxa. Our results suggest that species sorting (local environmental variables) operating through niche processes was the primary mechanism driving total beta diversity, nestedness and turnover for the entire community and active and passively dispersing taxa. In contrast, spatial factors (hydrological connectivity, proximity and spatial eigenvectors) only exerted a secondary influence on the nestedness and turnover components of beta diversity.     

An outcome model for human bladder cancer: A comprehensive study based on weighted gene co‐expression network analysis

The precision evaluation of prognosis is crucial for clinical treatment decision of bladder cancer (BCa). Therefore, establishing an effective prognostic model for BCa has significant clinical implications. We performed WGCNA and DEG screening to initially identify the candidate genes. The candidate genes were applied to construct a LASSO Cox regression analysis model. The effectiveness and accuracy of the prognostic model were tested by internal/external validation and pan‐cancer validation and time‐dependent ROC. Additionally, a nomogram based on the parameter selected from univariate and multivariate cox regression analysis was constructed. Eight genes were eventually screened out as progression‐related differentially expressed candidates in BCa. LASSO Cox regression analysis identified 3 genes to build up the outcome model in E‐MTAB‐4321 and the outcome model had good performance in predicting patient progress free survival of BCa patients in discovery and test set. Subsequently, another three datasets also have a good predictive value for BCa patients' OS and DFS. Time‐dependent ROC indicated an ideal predictive accuracy of the outcome model. Meanwhile, the nomogram showed a good performance and clinical utility. In addition, the prognostic model also exhibits good performance in pan‐cancer patients. Our outcome model was the first prognosis model for human bladder cancer progression prediction via integrative bioinformatics analysis, which may aid in clinical decision‐making.     

Distinguishing between turnover and nestedness in the quantification of biotic homogenization

Compositional changes through local extinction and colonization are inherent to natural communities, but human activities are increasingly influencing the rate and nature of the species being lost and gained. Biotic homogenization refers to the process by which the compositional similarity of communities increases over time through a non-random reshuffling of species. Despite the extensive conceptual development of the homogenization framework, approaches to quantify patterns of homogenization are scarcely developed. Most studies have used classical dissimilarity indices that actually quantify two components of compositional variation: turnover and nestedness. Here we demonstrate that a method that partitions those two components reveals patterns of homogenization that are otherwise obscured using traditional techniques. The forest understorey vegetation of an unmanaged reserve was recorded in permanent plots in 1979 and 2009. In only thirty years, the local species richness significantly decreased and the variation in the species composition from site to site shifted towards a structure with reduced true species turnover and increased dissimilarity due to nestedness. A classic analysis masked those patterns. In summary, we illustrated the need to move beyond the simple quantification of homogenization using classical indices and advocate integration of the multitude of ways to quantify community similarity into the homogenization framework.     

Boreal predator co‐occurrences reveal shared use of seismic lines in a working landscape

Interspecific interactions are an integral aspect of ecosystem functioning that may be disrupted in an increasingly anthropocentric world. Industrial landscape change creates a novel playing field on which these interactions take place, and a key question for wildlife managers is whether and how species are able to coexist in such working landscapes. Using camera traps deployed in northern Alberta, we surveyed boreal predators to determine whether interspecific interactions affected occurrences of black bears (Ursus americanus), coyotes (Canis latrans), and lynx (Lynx canadensis) within a landscape disturbed by networks of seismic lines (corridors cut for seismic exploration of oil and gas reserves). We tested hypotheses of species interactions across one spatial‐only and two spatiotemporal (daily and weekly) scales. Specifically, we hypothesized that (1) predators avoid competition with the apex predator, gray wolf (Canis lupus), (2) they avoid competition with each other as intraguild competitors, and (3) they overlap with their prey. All three predators overlapped with wolves on at least one scale, although models at the daily and weekly scale had substantial unexplained variance. None of the predators showed avoidance of intraguild competitors or overlap with prey. These results show patterns in predator space use that are consistent with both facilitative interactions or shared responses to unmeasured ecological cues. Our study provides insight into how predator species use the working boreal landscape in relation to each other, and highlights that predator management may indirectly influence multiple species through their interactions.     

Genetic diversity,population structure and sex‐biased dispersal in three co‐evolving species

Genetic diversity and spatial structure of populations are important for antagonistic coevolution. We investigated genetic variation and population structure of three closely related European ant species: the social parasite Harpagoxenus sublaevis and its two host species Leptothorax acervorum and Leptothorax muscorum. We sampled populations in 12 countries and analysed eight microsatellite loci and an mtDNA sequence. We found high levels of genetic variation in all three species, only slightly less variation in the host L. muscorum. Using a newly introduced measure of differentiation (Jost’s D_est ), we detected strong population structuring in all species and less male‐biased dispersal than previously thought. We found no phylogeographic patterns that could give information on post‐glacial colonization routes – northern populations are as variable as more southern populations. We conclude that conditions for Thompson’s geographic mosaic of coevolution are ideal in this system: all three species show ample genetic variation and strong population structure.