The spatial scaling of species interaction networks

The propagation of disturbances in ecological networks

Despite the development of network science, we lack clear heuristics for how far different disturbance types propagate within and across species interaction networks. We discuss the mechanisms of disturbance propagation in ecological networks, and propose that disturbances can be categorized into structural, functional, and transmission types according to their spread and effect on network structure and functioning. We describe the properties of species and their interaction networks and metanetworks that determine the indirect, spatial, and temporal extent of propagation. We argue that the sampling scale of ecological studies may have impeded predictions regarding the rate and extent that a disturbance spreads, and discuss directions to help ecologists to move towards a predictive understanding of the propagation of impacts across interacting communities and ecosystems.

Dissimilarity in flea and host assemblages and their interaction networks along a spatial distance gradient: different patterns revealed by different network dissimilarity metrics

We investigated the distance-decay pattern (an increase in dissimilarity with increasing geographic distance) in regional assemblages of fleas and their small mammalian hosts, as well as their interaction networks, in four biogeographic realms. Dissimilarity of assemblages (βtotal) was partitioned into species richness differences (βrich) and species replacement (βrepl) components. Dissimilarity of networks was assessed using two metrics: (a) whole network dissimilarity (βWN) partitioned into species replacement (βST) and interaction rewiring (βOS) components and (b) D statistics, measuring dissimilarity in the pure structure of the networks, without using information on species identities and calculated for hosts-shared-by-fleas networks (Dh) and fleas-shared-by-hosts networks (Df). We asked whether the distance-decay pattern (a) occurs among interactor assemblages or their interaction networks; (b) depends on the network dissimilarity metric used; and (c) differs between realms. The βtotal and βrepl of flea and host assemblages increased with distance in all realms except for host assemblages in the Afrotropics. βrich for flea and host assemblages increased with distance in the Nearctic only. In networks, βWN and βST demonstrated a distance-decay pattern, whereas βOS was mainly spatially invariant except in the Neotropics. Correlations of Dh or Df and geographic distance were mostly non-significant. We conclude that investigations of dissimilarity in interaction networks should include both types of dissimilarity metrics (those that consider partner identities and those that consider the pure structure of networks). This will allow elucidating the predictability of some facets of network dissimilarity and the unpredictability of other facets.

Both local stability and dispersal contribute to metacommunity sensitivity to asynchronous habitat availability

The stability of isolated communities depends on the complexity of their foodwebs. However, it remains unclear how local stability interacts with dispersal in multitrophic metacommunities to shape biodiversity patterns. This lack of understanding is deeper in the more realistic frame of landscapes that exhibit non-trivial and time-varying structures. Therefore, in this study, we aim to evaluate the influence of local stabilizing factors versus dispersal in determining the sensitivity of metacommunity biodiversity to increasing asynchrony of site availability. Additionally, we assess the role of foodweb complexity and landscape structure as modulating factors. To accomplish our goals we developed a model based on random matrices for local communities, which are linked by stochastic dispersal over explicit dynamic landscapes. We ran numerical simulations and computed the effect sizes of foodweb temperature, self-limitation, dispersal ability, and all pairwise combinations, on the sensitivity of biodiversity to landscape asynchrony. In our experiments we explored gradients of species richness, foodweb connectance, number of sites, and landscape modularity. Our results showed that asynchrony among site availability periods reduced α -diversity and increased β -diversity. Asynchrony increased γ -diversity at high dispersal rates. Both local and regional stabilizing factors determined the sensitivity of metacommunities to landscape asynchrony. Local factors were more influential in landscapes with fewer sites and lower modularity, as well as in metacommunities composed of complex foodwebs. This research offers insights into the dynamics of metacommunities in dynamic landscapes, providing valuable knowledge about the interplay between local and regional factors in shaping ecological stability and species persistence.

How collectively integrated are ecological communities?

Beyond abiotic conditions, do population dynamics mostly depend on a species' direct predators, preys and conspecifics? Or can indirect feedback that ripples across the whole community be equally important? Determining where ecological communities sit on the spectrum between these two characterizations requires a metric able to capture the difference between them. Here we show that the spectral radius of a community's interaction matrix provides such a metric, thus a measure of ecological collectivity, which is accessible from imperfect knowledge of biotic interactions and related to observable signatures. This measure of collectivity integrates existing approaches to complexity, interaction structure and indirect interactions. Our work thus provides an original perspective on the question of to what degree communities are more than loose collections of species or simple interaction motifs and explains when pragmatic reductionist approaches ought to suffice or fail when applied to ecological communities.

Land use conversion increases network complexity and stability of soil microbial communities in a temperate grassland

Soils harbor highly diverse microbial communities that are critical to soil health, but agriculture has caused extensive land use conversion resulting in negative effects on critical ecosystem processes. However, the responses and adaptations of microbial communities to land use conversion have not yet been understood. Here, we examined the effects of land conversion for long-term crop use on the network complexity and stability of soil microbial communities over 19 months. Despite reduced microbial biodiversity in comparison with native tallgrass prairie, conventionally tilled (CT) cropland significantly increased network complexity such as connectivity, connectance, average clustering coefficient, relative modularity, and the number of species acting at network hubs and connectors as well as resulted in greater temporal variation of complexity indices. Molecular ecological networks under CT cropland became significantly more robust and less vulnerable, overall increasing network stability. The relationship between network complexity and stability was also substantially strengthened due to land use conversion. Lastly, CT cropland decreased the number of relationships between network structure and environmental properties instead being strongly correlated to management disturbances. These results indicate that agricultural disturbance generally increases the complexity and stability of species "interactions", possibly as a trade-off for biodiversity loss to support ecosystem function when faced with frequent agricultural disturbance.

Geographic range size and species morphology determines the organization of sponge host-guest interaction networks across tropical coral reefs

Sponges are widely spread organisms in the tropical reefs of the American Northwest-Atlantic Ocean, they structure ecosystems and provide services such as shelter, protection from predators, and food sources to a wide diversity of both vertebrates and invertebrates species. The high diversity of sponge-associated fauna can generate complex networks of species interactions over small and large spatial-temporal gradients. One way to start uncovering the organization of the sponge host-guest complex networks is to understand how the accumulated geographic area, the sponge morphology and, sponge taxonomy contributes to the connectivity of sponge species within such networks. This study is a meta-analysis based on previous sponge host-guest literature obtained in 65 scientific publications, yielding a total of 745 host-guest interactions between sponges and their associated fauna across the Caribbean Sea and the Gulf of Mexico. We analyzed the sponge species contribution to network organization in the Northwest Tropical Atlantic coral reefs by using the combination of seven complementary species-level descriptors and related this importance with three main traits, sponge-accumulated geographic area, functional sponge morphology, and sponges' taxonomy bias. In general, we observed that sponges with a widespread distribution and a higher accumulated geographic area had a greater network structural contribution. Similarly, we also found that Cup-like and Massive functional morphologies trend to be shapes with a greater contribution to the interaction network organization compared to the Erect and Crust-like morphos. Lastly, we did not detect a taxonomy bias between interaction network organization and sponges' orders. These results highlight the importance of a specific combination of sponge traits to promote the diversity of association between reef sponges and their guest species.

Habitat area and edges affect the length of trophic chains in a fragmented forest

The food chain length represents how much energy reaches different trophic levels in food webs. Environmental changes derived from human activities have the potential to affect chain length. We explore how habitat area and edges affect chain length through: (1) a bottom-up effect of abundance ('pyramid hypothesis'); (2) the truncation of the highest trophic level ('trophic-rank hypothesis'); and (3) changes in species connectivity patterns ('connectivity hypothesis'). We built plant-leaf miner-parasitoid food webs in 19 remnants of a fragmented Chaco forest from central Argentina. On each remnant, we constructed food webs from different locations at the forest interior and edges. For each food web, we registered the abundance of species, the species richness of each trophic level, estimated the connectivity of their networks, and the average food chain length. We used structural equation models to evaluate the direct and indirect effects of habitat area and edge/interior location on food chain length mediated by species richness, abundance and connectivity. We found no direct effects of habitat area on chain length but chains were longer at forest edges than at their interior. The three mechanisms were supported by our results, although they showed different strengths. First, we found that the interior favours a bottom-up abundance effect from herbivores to parasitoids that positively affected chain length; second, we found that the forest area positively affects plant richness, which has a strong effect on the number of resources used by consumers, with a positive effect on chain length. Third, the remnant area and interior position favoured plant richness with a negative effect on the abundance of parasitoids, which had a positive effect on chain length. In general, the strongest effects on chain length were detected through changes in abundance rather than species richness although abundance was less affected by habitat fragmentation. We evaluated for the first time the effects of human-driven habitat fragmentation on the length of trophic chains in highly diverse plant-herbivore-parasitoid networks. Despite the loss of species, small habitat fragments and edges embedded in the agricultural matrix can support interaction networks, making them conservation targets in managed landscapes.

Strong succession in prokaryotic association networks and community assembly mechanisms in an acid mine drainage-impacted riverine ecosystem

Acid mine drainage (AMD) serves as an ideal model system for investigating microbial ecology, interaction, and assembly mechanism in natural environments. While previous studies have explored the structure and function of microbial communities in AMD, the succession patterns of microbial association networks and underlying assembly mechanisms during natural attenuation processes remain elusive. Here, we investigated prokaryotic microbial diversity and community assembly along an AMD-impacted river, from the extremely acidic, heavily polluted headwaters to the nearly neutral downstream sites. Microbial diversity was increased along the river, and microbial community composition shifted from acidophile-dominated to freshwater taxa-dominated communities. The complexity and relative modularity of the microbial networks were also increased, indicating greater network stability during succession. Deterministic processes, including abiotic selection of pH and high contents of sulfur and iron, governed community assembly in the headwaters. Although the stochasticity ratio was increased downstream, manganese content, microbial negative cohesion, and relative modularity played important roles in shaping microbial community structure. Overall, this study provides valuable insights into the ecological processes that govern microbial community succession in AMD-impacted riverine ecosystems. These findings have important implications for in-situ remediation of AMD contamination.

Spatial scaling of soil microbial co-occurrence networks in a fragmented landscape

Habitat loss has been a primary threat to biodiversity. However, species do not function in isolation but often associate with each other and form complex networks. Thus, revealing how the network complexity and stability scale with habitat area will give us more insights into the effects of habitat loss on ecosystems. In this study, we explored the relationships between the island area and the network complexity and stability of soil microbes. We found that the complexity and stability of soil microbial co-occurrence networks scale positively with island area, indicating that habitat loss will potentially simplify and destabilize soil microbial networks.

Colonization-competition dynamics of basal species shape food web complexity in island metacommunities

Exploring how food web complexity emerges and evolves in island ecosystems remains a major challenge in ecology. Food webs assembled from multiple islands are commonly recognized as highly complex trophic networks that are dynamic in both space and time. In the context of global climate change, it remains unclear whether food web complexity will decrease in a monotonic fashion when undergoing habitat destruction (e.g., the inundation of islands due to sea-level rise). Here, we develop a simple yet comprehensive patch-dynamic framework for complex food web metacommunities subject to the competition-colonization tradeoff between basal species. We found that oscillations in food web topological complexity (characterized by species diversity, mean food chain length and the degree of omnivory) emerge along the habitat destruction gradient. This outcome is robust to changing parameters or relaxing the assumption of a strict competitive hierarchy. Having oscillations in food web complexity indicates that small habitat changes could have disproportionate negative effects on species diversity, thus the success of conservation actions should be evaluated not only on changes in biodiversity, but also on system robustness to habitat alteration. Overall, this study provides a parsimonious mechanistic explanation for the emergence of food web complexity in island ecosystems, further enriching our understanding of metacommunity assembly.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00167-0.

Shortcomings of reusing species interaction networks created by different sets of researchers

Given the requisite cost associated with observing species interactions, ecologists often reuse species interaction networks created by different sets of researchers to test their hypotheses regarding how ecological processes drive network topology. Yet, topological properties identified across these networks may not be sufficiently attributable to ecological processes alone as often assumed. Instead, much of the totality of topological differences between networks-topological heterogeneity-could be due to variations in research designs and approaches that different researchers use to create each species interaction network. To evaluate the degree to which this topological heterogeneity is present in available ecological networks, we first compared the amount of topological heterogeneity across 723 species interaction networks created by different sets of researchers with the amount quantified from non-ecological networks known to be constructed following more consistent approaches. Then, to further test whether the topological heterogeneity was due to differences in study designs, and not only to inherent variation within ecological networks, we compared the amount of topological heterogeneity between species interaction networks created by the same sets of researchers (i.e., networks from the same publication) with the amount quantified between networks that were each from a unique publication source. We found that species interaction networks are highly topologically heterogeneous: while species interaction networks from the same publication are much more topologically similar to each other than interaction networks that are from a unique publication, they still show at least twice as much heterogeneity as any category of non-ecological networks that we tested. Altogether, our findings suggest that extra care is necessary to effectively analyze species interaction networks created by different researchers, perhaps by controlling for the publication source of each network.

Direct and indirect viral associations predict coexistence in wild plant virus communities

Viruses are a vastly underestimated component of biodiversity that occur as diverse communities across hierarchical scales from the landscape level to individual hosts. The integration of community ecology with disease biology is a powerful, novel approach that can yield unprecedented insights into the abiotic and biotic drivers of pathogen community assembly. Here, we sampled wild plant populations to characterize and analyze the diversity and co-occurrence structure of within-host virus communities and their predictors. Our results show that these virus communities are characterized by diverse, non-random coinfections. Using a novel graphical network modeling framework, we demonstrate how environmental heterogeneity influences the network of virus taxa and how the virus co-occurrence patterns can be attributed to non-random, direct statistical virus-virus associations. Moreover, we show that environmental heterogeneity changed virus association networks, especially through their indirect effects. Our results highlight a previously underestimated mechanism of how environmental variability can influence disease risks by changing associations between viruses that are conditional on their environment.

Changes in the bacterial rare biosphere after permanent application of composted tannery sludge in a tropical soil

Composted tannery sludge (CTS) promotes shifts in soil chemical properties, affecting microbial communities. Although the effect of CTS application on the bacterial community has been studied, it is unclear whether this impact discriminates between the dominant and rare species. This present study investigated how the dominant and rare bacterial communities respond over time to different concentrations of CTS application (0, 2.5, 5, 10, and 20 tons/ha) for 180 days. The richness of operational taxonomic units (OTU) was 30-fold higher in the rare than in the dominant biosphere. While some phyla shifted their relative abundance differently in the dominant and rare biosphere, some genera increased their relative abundance under higher CTS concentrations, such as Nocardioides (∼100%), Rubrobacter (∼300%), and Nordella (∼400%). Undominated processes largely governed the dominant biosphere (76.97%), followed by homogeneous (12.51%) and variable (8.03%) selection, and to a lesser extent, the dispersal limitation (2.48%). The rare biosphere was driven by the CTS application as evidenced by the exclusively homogeneous selection (100%). This study showed that the rare biosphere was more sensitive to changes in soil chemical parameters due to CTS application, which evidences the importance explore this portion of the bacterial community for its biotechnological use in contaminated soils.

Fungal diversities and community assembly processes show different biogeographical patterns in forest and grassland soil ecosystems

Soil fungal community has been largely explored by comparing their natural diversity. However, there is a relatively small body of literature concerned with fungal community assembly processes and their co-occurrence network correlations carried out across large spatial-temporal scales with complex environmental gradients in natural ecosystems and different habitats in China. Thus, soil fungal community assembly processes were assessed to predict changes in soil function in 98 different forest and grassland sites from the Sichuan, Hubei, and Hebei Provinces of China using high-throughput sequencing of nuclear ribosomal internal transcribed spacer 2 (ITS-2). The 10 most abundant fungal phyla results showed that Ascomycota was the most abundant phylum in forests from Sichuan province (64.42%) and grassland habitats from Hebei province (53.46%). Moreover, core fungal taxa (487 OTUs) represented 0.35% of total fungal OTUs. We observed higher fungal Shannon diversity and richness (the Chao1 index) from diverse mixed forests of the Sichuan and Hubei Provinces than the mono-cultured forest and grassland habitats in Hebei Province. Although fungal alpha and beta diversities exhibited different biogeographical patterns, the fungal assembly pattern was mostly driven by dispersal limitation than selection in different habitats. Fungal co-occurrence analyses showed that the network was more intense at Saihanba National Forest Park (SNFP, Hebei). In contrast, the co-occurrence network was more complex at boundaries between forests and grasslands at SNFP. Additionally, the highest number of positive (co-presence or co-operative) correlations of fungal genera were inferred from grassland habitat, which led fungal communities to form commensalism relationships compared to forest areas with having higher negative correlations (mutual exclusion or competitive). The generalized additive model (GAM) analysis showed that the association of fungal Shannon diversity and richness indices with geographical coordinates did not follow a general pattern; instead, the fluctuation of these indices was restricted to local geographical coordinates at each sampling location. These results indicated the existence of a site effect on the diversity of fungal communities across our sampling sites. Our observation suggested that higher fungal diversity and richness of fungal taxa in a particular habitat are not necessarily associated with more complex networks.

Plant-frugivore network simplification under habitat fragmentation leaves a small core of interacting generalists

Habitat fragmentation impacts seed dispersal processes that are important in maintaining biodiversity and ecosystem functioning. However, it is still unclear how habitat fragmentation affects frugivorous interactions due to the lack of high-quality data on plant-frugivore networks. Here we recorded 10,117 plant-frugivore interactions from 22 reservoir islands and six nearby mainland sites using the technology of arboreal camera trapping to assess the effects of island area and isolation on the diversity, structure, and stability of plant-frugivore networks. We found that network simplification under habitat fragmentation reduces the number of interactions involving specialized species and large-bodied frugivores. Small islands had more connected, less modular, and more nested networks that consisted mainly of small-bodied birds and abundant plants, as well as showed evidence of interaction release (i.e., dietary expansion of frugivores). Our results reveal the importance of preserving large forest remnants to support plant-frugivore interaction diversity and forest functionality.

Smaller communities, such as those on islands, under ecological network simplification reduce interactions between specialist organisms.

Island and Mountain Ecosystems as Testbeds for Biological Control in the Anthropocene

For centuries, islands and mountains have incited the interest of naturalists, evolutionary biologists and ecologists. Islands have been the cradle for biogeography and speciation theories, while mountain ranges have informed how population adaptation to thermal floors shapes the distribution of species globally. Islands of varying size and mountains’ altitudinal ranges constitute unique “natural laboratories” where one can investigate the effects of species loss or global warming on ecosystem service delivery. Although invertebrate pollination or seed dispersal processes are steadily being examined, biological control research is lagging. While observations of a wider niche breadth among insect pollinators in small (i.e., species-poor) islands or at high (i.e., colder) altitudes likely also hold for biological control agents, such remains to be examined. In this Perspective piece, we draw on published datasets to show that island size alone does not explain biological control outcomes. Instead, one needs to account for species’ functional traits, habitat heterogeneity, host community make-up, phenology, site history or even anthropogenic forces. Meanwhile, data from mountain ranges show how parasitism rates of Noctuid moths and Tephritid fruit flies exhibit species- and context-dependent shifts with altitude. Nevertheless, future empirical work in mountain settings could clarify the thermal niche space of individual natural enemy taxa and overall thermal resilience of biological control. We further discuss how global databases can be screened, while ecological theories can be tested, and simulation models defined based upon observational or manipulative assays in either system. Doing so can yield unprecedented insights into the fate of biological control in the Anthropocene and inform ways to reinforce this vital ecosystem service under global environmental change scenarios.

How network size strongly determines trophic specialisation: A technical comment on Luna et al. (2022)

Luna et al. (2022) concluded that the environment contributes to explaining specialisation in open plant-pollinator networks. When reproducing their study, we instead found that network size alone largely explained the variation in their specialisation metrics. Thus, we question whether empirical network specialisation is driven by the environment.

A patch-dynamic metacommunity perspective on the persistence of mutualistic and antagonistic bipartite networks

The structure of interactions between species within a community plays a key role in maintaining biodiversity. Previous studies have found that the effects of these structures might substantially differ depending on interaction type, for example, a highly connected and nested architecture stabilizes mutualistic communities, while the stability of antagonistic communities is enhanced in modular and weakly connected structures. Here we show that, when network dynamics are modelled using a patch-dynamic metacommunity framework, the qualitative differences between antagonistic and mutualistic systems disappear, with nestedness and modularity interacting to promote metacommunity persistence. However, the interactive effects are significantly weaker in antagonistic metacommunities. Our model also predicts an increase in connectance, nestedness and modularity over time in both types of interaction, except in antagonistic networks where nestedness declines. At steady state, we find a strong negative correlation between nestedness and modularity in both mutualistic and antagonistic metacommunities. These predictions are consistent with the structural trends found in a large dataset of real-world antagonistic and mutualistic communities.

Biotic filtering by species' interactions constrains food-web variability across spatial and abiotic gradients

Despite intensive research on species dissimilarity patterns across communities (i.e. β-diversity), we still know little about their implications for variation in food-web structures. Our analyses of 50 lake and 48 forest soil communities show that, while species dissimilarity depends on environmental and spatial gradients, these effects are only weakly propagated to the networks. Moreover, our results show that species and food-web dissimilarities are consistently correlated, but that much of the variation in food-web structure across spatial, environmental, and species gradients remains unexplained. Novel food-web assembly models demonstrate the importance of biotic filtering during community assembly by (1) the availability of resources and (2) limiting similarity in species' interactions to avoid strong niche overlap and thus competitive exclusion. This reveals a strong signature of biotic filtering processes during local community assembly, which constrains the variability in structural food-web patterns across local communities despite substantial turnover in species composition.

Ecological network complexity scales with area

Larger geographical areas contain more species-an observation raised to a law in ecology. Less explored is whether biodiversity changes are accompanied by a modification of interaction networks. We use data from 32 spatial interaction networks from different ecosystems to analyse how network structure changes with area. We find that basic community structure descriptors (number of species, links and links per species) increase with area following a power law. Yet, the distribution of links per species varies little with area, indicating that the fundamental organization of interactions within networks is conserved. Our null model analyses suggest that the spatial scaling of network structure is determined by factors beyond species richness and the number of links. We demonstrate that biodiversity-area relationships can be extended from species counts to higher levels of network complexity. Therefore, the consequences of anthropogenic habitat destruction may extend from species loss to wider simplification of natural communities.

Dispersal Limitation Plays Stronger Role in the Community Assembly of Fungi Relative to Bacteria in Rhizosphere Across the Arable Area of Medicinal Plant

Understanding the ecological patterns of rhizosphere microbial communities is critical for propelling sustainable agriculture and managing ecosystem functions by exploiting microorganisms. However, this knowledge is still unclear, especially under host-associated large-scale and regarding the comparison between bacteria and fungi. We examined community assembly processes and community characters including environmental thresholds and co-occurrence patterns across the cultivatable area of Panax notoginseng for bacteria and fungi. Both are vital members of the rhizosphere but differ considerably in their life history and dispersal potentiality. Edaphic factors drove the parallel variations of bacterial and fungal communities. Although bacterial and fungal communities exhibited similar biogeographic patterns, the assembly of fungi was more driven by dispersal limitation than selection compared with bacteria. This finding supported the 'size-dispersal' hypothesis. pH and total nitrogen respectively mediated the relative importance of deterministic and stochastic processes in shaping bacterial and fungal communities. In addition, fungal communities exhibited potentially broader environmental thresholds and more modular co-occurrence patterns than bacteria (bacteria: 0.67; fungi: 0.78). These results emphasized the importance of dispersal limitation in structuring rhizosphere microbiota and shaping community features of ecologically distinct microorganisms. This study provides insights into the improved prediction and management of the key functions of rhizosphere microbiota.

Effect of aerobic/anoxic duration on the performance, microbial activity and microbial community of sequencing batch biofilm reactor treating synthetic mariculture wastewater

The effect of aerobic/anoxic duration on the performance, microbial community and enzymatic activity of sequencing batch biofilm reactor (SBBR) were investigated in treating mariculture wastewater. The microbial oxygen uptake rate and nitrifying rate gradually decreased with the aerobic/anoxic duration from 120/210 to 30/300 min, whereas the nitrite reducing rate and nitrate reducing rate had the opposite results. The activities of dehydrogenase, ammonia monooxygenase and nitrite oxidoreductase gradually decreased with the aerobic/anoxic duration from 120/210 to 30/300 min, but the activities of nitrate reductase and nitrite reductase had a gradual increment. The microbial nitrogen removal rates had similar changing trends to their corresponding enzymatic activities at different aerobic/anoxic duration. The variation of aerobic/anoxic duration obviously affected the microbial richness and diversity of SBBR. The co-occurrence, keystone taxa and significant difference of microbial community had some changes with the aerobic/anoxic duration from 120/210 to 30/300 min.

The spatial scaling of food web structure across European biogeographical regions

The species-area relationship (SAR) is one of the most well-established scaling patterns in ecology. Its implications for understanding how communities change across spatial gradients are numerous, including the effects of habitat loss on biodiversity. However, ecological communities are not mere collections of species. They are the result of interactions between these species forming complex networks that tie them together. Should we aim to grasp the spatial scaling of biodiversity as a whole, it is fundamental to understand the changes in the structure of interaction networks with area. In spite of a few empirical and theoretical studies that address this challenge, we still do not know much about how network structure changes with area, or what are the main environmental drivers of these changes. Here, using the meta-network of potential interactions between all terrestrial vertebrates in Europe (1140 species and 67 201 feeding interactions), we analysed network-area relationships (NARs) that summarize how network properties scale with area. We do this across ten biogeographical regions, which differ in environmental characteristics. We found that the spatial scaling of network complexity strongly varied across biogeographical regions. However, once the variation in SARs was accounted for, differences in the shape of NARs vanished. On the other hand, the proportion of species across trophic levels remained remarkably constant across biogeographical regions and spatial scales, despite the great variation in species richness. Spatial variation in mean annual temperature and habitat clustering were the main environmental determinants of the shape of both SARs and NARs across Europe. Our results suggest new avenues in the exploration of the effects of environmental factors on the spatial scaling of biodiversity. We argue that NARs can provide new insights to analyse and understand ecological communities.

Network ecology in dynamic landscapes

Network ecology is an emerging field that allows researchers to conceptualize and analyse ecological networks and their dynamics. Here, we focus on the dynamics of ecological networks in response to environmental changes. Specifically, we formalize how network topologies constrain the dynamics of ecological systems into a unifying framework in network ecology that we refer to as the 'ecological network dynamics framework'. This framework stresses that the interplay between species interaction networks and the spatial layout of habitat patches is key to identifying which network properties (number and weights of nodes and links) and trade-offs among them are needed to maintain species interactions in dynamic landscapes. We conclude that to be functional, ecological networks should be scaled according to species dispersal abilities in response to landscape heterogeneity. Determining how such effective ecological networks change through space and time can help reveal their complex dynamics in a changing world.

Species-area and network-area relationships in host-helminth interactions

The scaling relationship observed between species richness and the geographical area sampled (i.e. the species-area relationship (SAR)) is a widely recognized macroecological relationship. Recently, this theory has been extended to trophic interactions, suggesting that geographical area may influence the structure of species interaction networks (i.e. network-area relationships (NARs)). Here, we use a global dataset of host-helminth parasite interactions to test existing predictions from macroecological theory. Scaling between single locations to the global host-helminth network by sequentially adding networks together, we find support that geographical area influences species richness and the number of species interactions in host-helminth networks. However, species-area slopes were larger for host species relative to their helminth parasites, counter to theoretical predictions. Lastly, host-helminth network modularity-capturing the tendency of the network to form into separate subcommunities-decreased with increasing area, also counter to theoretical predictions. Reconciling this disconnect between existing theory and observed SAR and NAR will provide insight into the spatial structuring of ecological networks, and help to refine theory to highlight the effects of network type, species distributional overlap, and the specificity of trophic interactions on NARs.

Island Colonization and Environmental Sustainability in the Postglacial Mediterranean

Island environments present challenges to human colonization, but we have a poor understanding of how environmental difference drives heterogeneous patterns of insular settlement. In this paper, we assess which environmental and geographic variables positively or negatively affect the long-term sustainability of human settlement on islands. Using the postglacial Mediterranean basin as a case study, we assess the impact of area, isolation index, species richness, and net primary productivity (NPP) on patterns of island occupation for both hunter-gatherer and agropastoral populations. We find that models involving area most effectively accounts for sustainability in hunter-gatherer island settlement. The agropastoral data are noisier, perhaps due to culturally specific factors responsible for the distribution of the data; nonetheless, we show that area and NPP exert profound influence over sustainability of agropastoral island settlement. We conclude by suggesting that this relates to the capacity of these variables to impact demographic robusticity directly.

Relative effects of anthropogenic pressures, climate, and sampling design on the structure of pollination networks at the global scale

Pollinators provide crucial ecosystem services that underpin to wild plant reproduction and yields of insect-pollinated crops. Understanding the relative impacts of anthropogenic pressures and climate on the structure of plant-pollinator interaction networks is vital considering ongoing global change and pollinator decline. Our ability to predict the consequences of global change for pollinator assemblages worldwide requires global syntheses, but these analytical approaches may be hindered by variable methods among studies that either invalidate comparisons or mask biological phenomena. Here we conducted a synthetic analysis that assesses the relative impact of anthropogenic pressures and climatic variability, and accounts for heterogeneity in sampling methodology to reveal network responses at the global scale. We analyzed an extensive dataset, comprising 295 networks over 123 locations all over the world, and reporting over 50,000 interactions between flowering plant species and their insect visitors. Our study revealed that anthropogenic pressures correlate with an increase in generalism in pollination networks while pollinator richness and taxonomic composition are more related to climatic variables with an increase in dipteran pollinator richness associated with cooler temperatures. The contrasting response of species richness and generalism of the plant-pollinator networks stresses the importance of considering interaction network structure alongside diversity in ecological monitoring. In addition, differences in sampling design explained more variation than anthropogenic pressures or climate on both pollination networks richness and generalism, highlighting the crucial need to report and incorporate sampling design in macroecological comparative studies of pollination networks. As a whole, our study reveals a potential human impact on pollination networks at a global scale. However, further research is needed to evaluate potential consequences of loss of specialist species and their unique ecological interactions and evolutionary pathways on the ecosystem pollination function at a global scale.

Trait positions for elevated invasiveness in adaptive ecological networks

Our ability to predict the outcome of invasion declines rapidly as non-native species progress through intertwined ecological barriers to establish and spread in recipient ecosystems. This is largely due to the lack of systemic knowledge on key processes at play as species establish self-sustaining populations within the invaded range. To address this knowledge gap, we present a mathematical model that captures the eco-evolutionary dynamics of native and non-native species interacting within an ecological network. The model is derived from continuous-trait evolutionary game theory (i.e., Adaptive Dynamics) and its associated concept of invasion fitness which depicts dynamic demographic performance that is both trait mediated and density dependent. Our approach allows us to explore how multiple resident and non-native species coevolve to reshape invasion performance, or more precisely invasiveness, over trait space. The model clarifies the role of specific traits in enabling non-native species to occupy realised opportunistic niches. It also elucidates the direction and speed of both ecological and evolutionary dynamics of residing species (natives or non-natives) in the recipient network under different levels of propagule pressure. The versatility of the model is demonstrated using four examples that correspond to the invasion of (i) a horizontal competitive community; (ii) a bipartite mutualistic network; (iii) a bipartite antagonistic network; and (iv) a multi-trophic food web. We identified a cohesive trait strategy that enables the success and establishment of non-native species to possess high invasiveness. Specifically, we find that a non-native species can achieve high levels of invasiveness by possessing traits that overlap with those of its facilitators (and mutualists), which enhances the benefits accrued from positive interactions, and by possessing traits outside the range of those of antagonists, which mitigates the costs accrued from negative interactions. This ‘central-to-reap, edge-to-elude’ trait strategy therefore describes the strategic trait positions of non-native species to invade an ecological network. This model provides a theoretical platform for exploring invasion strategies in complex adaptive ecological networks.

Landscape modification and nutrient-driven instability at a distance

Almost 50 years ago, Michael Rosenzweig pointed out that nutrient addition can destabilise food webs, leading to loss of species and reduced ecosystem function through the paradox of enrichment. Around the same time, David Tilman demonstrated that increased nutrient loading would also be expected to cause competitive exclusion leading to deleterious changes in food web diversity. While both concepts have greatly illuminated general diversity-stability theory, we currently lack a coherent framework to predict how nutrients influence food web stability across a landscape. This is a vitally important gap in our understanding, given mounting evidence of serious ecological disruption arising from anthropogenic displacement of resources and organisms. Here, we combine contemporary theory on food webs and meta-ecosystems to show that nutrient additions are indeed expected to drive loss in stability and function in human-impacted regions. Our models suggest that destabilisation is more likely to be caused by the complete loss of an equilibrium due to edible plant species being competitively excluded. In highly modified landscapes, spatial nutrient transport theory suggests that such instabilities can be amplified over vast distances from the sites of nutrient addition. Consistent with this theoretical synthesis, the empirical frequency of these distant propagating ecosystem imbalances appears to be growing. This synthesis of theory and empirical data suggests that human modification of the Earth is strongly connecting distantly separated ecosystems, causing rapid, expansive and costly nutrient-driven instabilities over vast areas of the planet. Similar to existing food web theory, the corollary to this spatial nutrient theory is that slowing down spatial nutrient pathways can be a potent means of stabilising degraded ecosystems.

The Structure of Ecological Networks Across Levels of Organization

Interactions connect the units of ecological systems, forming networks. Individual-based networks characterize variation in niches among individuals within populations. These individual-based networks merge with each other, forming species-based networks and food webs that describe the architecture of ecological communities. Networks at broader spatiotemporal scales portray the structure of ecological interactions across landscapes and over macroevolutionary time. Here, I review the patterns observed in ecological networks across multiple levels of biological organization. A fundamental challenge is to understand the amount of interdependence as we move from individual-based networks to species-based networks and beyond. Despite the uneven distribution of studies, regularities in network structure emerge across scales due to the fundamental architectural patterns shared by complex networks and the interplay between traits and numerical effects. I illustrate the integration of these organizational scales by exploring the consequences of the emergence of highly connected species for network structures across scales.

Revisiting the Links-Species Scaling Relationship in Food Webs

Predicting the number of interactions among species in a food web is an important task. These trophic interactions underlie many ecological and evolutionary processes, ranging from biomass fluxes, ecosystem stability, resilience to extinction, and resistance against novel species. We investigate and compare several ways to predict the number of interactions in food webs. We conclude that a simple beta-binomial model outperforms other models, with the added desirable property of respecting biological constraints. We show how this simple relationship gives rise to a predicted distribution of several quantities related to link number in food webs, including the scaling of network structure with space and the probability that a network will be stable.

Towards an integrative understanding of soil biodiversity

Soil is one of the most biodiverse terrestrial habitats. Yet, we lack an integrative conceptual framework for understanding the patterns and mechanisms driving soil biodiversity. One of the underlying reasons for our poor understanding of soil biodiversity patterns relates to whether key biodiversity theories (historically developed for aboveground and aquatic organisms) are applicable to patterns of soil biodiversity. Here, we present a systematic literature review to investigate whether and how key biodiversity theories (species-energy relationship, theory of island biogeography, metacommunity theory, niche theory and neutral theory) can explain observed patterns of soil biodiversity. We then discuss two spatial compartments nested within soil at which biodiversity theories can be applied to acknowledge the scale-dependent nature of soil biodiversity.

Scaling-up biodiversity-ecosystem functioning research

A rich body of knowledge links biodiversity to ecosystem functioning (BEF), but it is primarily focused on small scales. We review the current theory and identify six expectations for scale dependence in the BEF relationship: (1) a nonlinear change in the slope of the BEF relationship with spatial scale; (2) a scale‐dependent relationship between ecosystem stability and spatial extent; (3) coexistence within and among sites will result in a positive BEF relationship at larger scales; (4) temporal autocorrelation in environmental variability affects species turnover and thus the change in BEF slope with scale; (5) connectivity in metacommunities generates nonlinear BEF and stability relationships by affecting population  synchrony at local and regional scales; (6) spatial scaling in food web structure and diversity will generate scale dependence in ecosystem functioning. We suggest directions for synthesis that combine approaches in metaecosystem and metacommunity ecology and integrate cross‐scale feedbacks. Tests of this theory may combine remote sensing with a generation of networked experiments that assess effects at multiple scales. We also show how anthropogenic land cover change may alter the scaling of the BEF relationship. New research on the role of scale in BEF will guide policy linking the goals of managing biodiversity and ecosystems.

A multitrophic perspective on biodiversity-ecosystem functioning research

Concern about the functional consequences of unprecedented loss in biodiversity has prompted biodiversity-ecosystem functioning (BEF) research to become one of the most active fields of ecological research in the past 25 years. Hundreds of experiments have manipulated biodiversity as an independent variable and found compelling support that the functioning of ecosystems increases with the diversity of their ecological communities. This research has also identified some of the mechanisms underlying BEF relationships, some context-dependencies of the strength of relationships, as well as implications for various ecosystem services that mankind depends upon. In this paper, we argue that a multitrophic perspective of biotic interactions in random and non-random biodiversity change scenarios is key to advance future BEF research and to address some of its most important remaining challenges. We discuss that the study and the quantification of multitrophic interactions in space and time facilitates scaling up from small-scale biodiversity manipulations and ecosystem function assessments to management-relevant spatial scales across ecosystem boundaries. We specifically consider multitrophic conceptual frameworks to understand and predict the context-dependency of BEF relationships. Moreover, we highlight the importance of the eco-evolutionary underpinnings of multitrophic BEF relationships. We outline that FAIR data (meeting the standards of findability, accessibility, interoperability, and reusability) and reproducible processing will be key to advance this field of research by making it more integrative. Finally, we show how these BEF insights may be implemented for ecosystem management, society, and policy. Given that human well-being critically depends on the multiple services provided by diverse, multitrophic communities, integrating the approaches of evolutionary ecology, community ecology, and ecosystem ecology in future BEF research will be key to refine conservation targets and develop sustainable management strategies.

Reshaping our understanding of species' roles in landscape-scale networks

In network ecology, landscape-scale processes are often overlooked, yet there is increasing evidence that species and interactions spill over between habitats, calling for further study of interhabitat dependencies. Here, we investigate how species connect a mosaic of habitats based on the spatial variation of their mutualistic and antagonistic interactions using two multilayer networks, combining pollination, herbivory and parasitism in the UK and New Zealand. Developing novel methods of network analysis for landscape-scale ecological networks, we discovered that few plant and pollinator species acted as connectors or hubs, both within and among habitats, whereas herbivores and parasitoids typically have more peripheral network roles. Insect species' roles depend on factors other than just the abundance of taxa in the lower trophic level, exemplified by larger Hymenoptera connecting networks of different habitats and insects relying on different resources across different habitats. Our findings provide a broader perspective for landscape-scale management and ecological community conservation.

The geographical variation of network structure is scale dependent: understanding the biotic specialization of host-parasitoid networks

Research on the structure of ecological networks suggests that a number of universal patterns exist. Historically, biotic specialization has been thought to increase towards the Equator. Yet, recent studies have challenged this view showing non-conclusive results. Most studies analysing the geographical variation in biotic specialization focus, however, only on the local scale. Little is known about how the geographical variation of network structure depends on the spatial scale of observation (i.e., from local to regional spatial scales). This should be remedied, as network structure changes as the spatial scale of observation changes, and the magnitude and shape of these changes can elucidate the mechanisms behind the geographical variation in biotic specialization. Here we analyse four facets of biotic specialization in host-parasitoid networks along gradients of climatic constancy, classifying the networks according to their spatial extension (local or regional). Namely, we analyse network connectance, consumer diet overlap, consumer diet breadth, and resource vulnerability at both local and regional scales along the gradients of both current climatic constancy and historical climatic change. While at the regional scale none of the climatic variables are associated to biotic specialization, at the local scale, network connectance, consumer diet overlap, and resource vulnerability decrease with current climatic constancy, whereas consumer generalism increases (i.e., broader diet breadths in tropical areas). Similar patterns are observed along the gradient of historical climatic change. We provide an explanation based on different beta-diversity for consumers and resources across the geographical gradients. Our results show that the geographical gradient of biotic specialization is not universal. It depends on both the facet of biotic specialization and the spatial scale of observation.

The three regimes of spatial recovery

An enduring challenge for ecology is identifying the drivers of ecosystem and population stability. In a spatially explicit context, key features to consider are landscape spatial structure, local interactions, and dispersal. Substantial work has been done on each of these features as a driver of stability, but little is known on the interplay between them. Missing has been a more integrative approach, able to map and identify different dynamical regimes, predicting a system's response to perturbations. Here we first consider a simple scenario, i.e., the recovery of a homogeneous metapopulation from a single localized pulse disturbance. The analysis of this scenario reveals three fundamental recovery regimes: Isolated Regime when dispersal is not significant, Rescue Regime when dispersal mediates recovery, and Mixing Regime when perturbations spread throughout the system. Despite its simplicity, our approach leads to remarkably general predictions. These include the qualitatively different outcomes of various scenarios of habitat fragmentation, the surprising benefits of local extinctions on population persistence at the transition between regimes, and the productivity shifts of metacommunities in a changing environment. This study thus provides context to known results and insight into future directions of research.

Modular structure in fish co-occurrence networks: A comparison across spatial scales and grouping methodologies

Network modules are used for diverse purposes, ranging from delineation of biogeographical provinces to the study of biotic interactions. We assess spatial scaling effects on modular structure, using a multi-step process to compare fish co-occurrence networks at three nested scales. We first detect modules with simulated annealing and use spatial clustering tests (interspecific distances among species' range centroids) to determine if modules consist of species with broadly overlapping ranges; strong spatial clustering may reflect environmental filtering, while absence of spatial clustering may reflect positive interspecific relationships (commensalism or mutualism). We then use non-hierarchical, multivariate cluster analysis as an alternative method to identify fish subgroups, we repeat spatial clustering tests for the multivariate clusters, then compare spatial clustering results among modules and clusters. Next, we compare species lists within modules and clusters, and estimate congruence as the proportion of species assigned to the same groups by the two methods. Finally, we use a well-documented nest associate complex (fishes that deposit eggs in the gravel nests of a common host) to assess whether strong within-group associations may, in fact, reflect positive interspecific relationships. At each scale, 2-4 network modules were detected but a consistent relationship between scale and the number of modules was not observed. Significant spatial clustering was detected at all scales for network modules and multivariate clusters but was less prevalent at smaller scales. Congruence between modules and clusters was always < 90% and generally decreased as the number of groups increased. At all scales, the complete nest associate complex was completely preserved within a single network module, but not within a single multivariate cluster. Collectively, our results suggest that network modules are promising tools for studying positive interactions and that smaller scales may be preferable in this research.

Towards a multi-trophic extension of metacommunity ecology

Metacommunity theory provides an understanding of how spatial processes determine the structure and function of communities at local and regional scales. Although metacommunity theory has considered trophic dynamics in the past, it has been performed idiosyncratically with a wide selection of possible dynamics. Trophic metacommunity theory needs a synthesis of a few influential axis to simplify future predictions and tests. We propose an extension of metacommunity ecology that addresses these shortcomings by incorporating variability among trophic levels in 'spatial use properties'. We define 'spatial use properties' as a set of traits (dispersal, migration, foraging and spatial information processing) that set the spatial and temporal scales of organismal movement, and thus scales of interspecific interactions. Progress towards a synthetic predictive framework can be made by (1) documenting patterns of spatial use properties in natural food webs and (2) using theory and experiments to test how trophic structure in spatial use properties affects metacommunity dynamics.

Network Analyses Can Advance Above-Belowground Ecology

An understanding of above-belowground (AG-BG) ecology is important for evaluating how plant interactions with enemies, symbionts, and decomposers affect species diversity and will respond to global changes. However, research questions and experiments often focus on only a limited number of interactions, creating an incomplete picture of how entire communities may be involved in AG-BG community ecology. Therefore, a pressing challenge is to formulate hypotheses of AG-BG interactions when considering communities in their full complexity. Here we discuss how network analyses can be a powerful tool to progress AG-BG research, link across scales from individual to community and ecosystem, visualize community interactions between the two (AG and BG) subsystems, and develop testable hypotheses.