Energetics, Patterns of Interaction Strengths, and Stability in Real Ecosystems

The long-term changes in food web structure and ecosystem functioning of a shallow lake: Implications for the lake management

The food web structure (FWS) and ecosystem functioning (EF) of lakes worldwide are impacted by multiple disturbances. The historical evolution of the FWS and EF are not well understood due to the lack of sufficient long-term records of biotic variations. This study reconstructed the food web models in the 1950s, 1980s, 1990s, 2000s, and 2010s for Baiyangdian Lake (BYDL), the largest shallow lake in northern China, using the Linear Inverse Modeling (LIM) and investigated EF in different periods. Our results confirmed that the FWS and EF of BYDL have undergone significant changes. The biomass of phytoplankton continuously increased, and the primary productivity of phytoplankton began to replace the primary productivity of submerged macrophytes in the 2000s and became the largest energy flow in the food web. Changes in the energy flow of primary producers are transmitted to high-trophic functional groups, which affects the diet composition of fish. Based on the ecological network analysis indices and food web stability indicators, it was concluded that after a turning point in the 1990s, the ecosystem showed initial stability and then gradually became unstable. Water level fluctuations and nutrient enrichment may be the key driving factors for changes in ecosystem state. Therefore, to maintain a good state of the ecosystem, we recommend implementing comprehensive management measures of hydrological management, nutrient-loading reduction, and biomanipulation. Furthermore, this study extended LIM to lake ecosystems, which may provide a new method for lake ecological environment management.

Ecological network assembly: how the regional metaweb influences local food webs

1. Local food webs result from a sequence of colonisations and extinctions by species from the regional pool or metaweb, i.e., the assembly process. Assembly is theorised to be a selective process: whether or not certain species or network structures can persist is partly determined by local processes including habitat filtering and dynamical constraints. Consequently, local food web structure should reflect these processes. 2. The goal of this study was to test evidence for these selective processes by comparing the structural properties of real food webs to the expected distribution given the metaweb. We were particularly interested in ecological dynamics; if the network properties commonly associated with dynamical stability are indeed the result of stability constraints, then they should deviate from expectation in the direction predicted by theory. 3. To create a null expectation, we used the novel approach of randomly assembling model webs by drawing species and interactions from the empirical metaweb. The assembly model permitted colonisation and extinction, and required a consumer species to have at least one prey, but had no habitat type nor population dynamical constraints. Three data sets were used: (1) the marine Antarctic metaweb, with 2 local food-webs; (2) the 50 lakes of the Adirondacks; and (3) the arthropod community from Florida Keys' classic defaunation experiment. 4. Contrary to our expectations, we found that there were almost no differences between empirical webs and those resulting from the null assembly model. Few empirical food webs showed significant differences with network properties, motif representations and topological roles. Network properties associated with stability did not deviate from expectation in the direction predicted by theory. 5. Our results suggest that - for the commonly used metrics we considered - local food web structure is not strongly influenced by dynamical nor habitat restrictions. Instead, the structure is inherited from the metaweb. This suggests that the network properties typically attributed as causes or consequences of ecological stability are instead a by-product of the assembly process (i.e., spandrels), and may potentially be too coarse to detect the true signal of dynamical constraint.

Microbial Diversity Characteristics of Areca Palm Rhizosphere Soil at Different Growth Stages

The rhizosphere microflora are key determinants that contribute to plant health and productivity, which can support plant nutrition and resistance to biotic and abiotic stressors. However, limited research is conducted on the areca palm rhizosphere microbiota. To further study the effect of the areca palm’s developmental stages on the rhizosphere microbiota, the rhizosphere microbiota of areca palm (Areca catechu) grown in its main producing area were examined in Wanning, Hainan province, at different vegetation stages by an Illumina Miseq sequence analysis of the 16S ribosomal ribonucleic acid and internal transcribed spacer genes. Significant shifts of the taxonomic composition of the bacteria and fungi were observed in the four stages. Burkholderia-Caballeronia-Paraburkholderia were the most dominant group in stage T1 and T2; the genera Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium were decreased significantly from T1 to T2; and the genera Acidothermus and Bacillus were the most dominant in stage T3 and T4, respectively. Meanwhile, Neocosmospora, Saitozyma, Penicillium, and Trichoderma were the most dominant genera in the stage T1, T2, T3, and T4, respectively. Among the core microbiota, the dominant bacterial genera were Burkholderia-Caballeronia-Paraburkholderia and Bacillus, and the dominant fungal genera were Saitozyma and Trichoderma. In addition, we identified five bacterial genera and five fungal genera that reached significant levels during development. Finally, we constructed the OTU (top 30) interaction network of bacteria and fungi, revealed its interaction characteristics, and found that the bacterial OTUs exhibited more extensive interactions than the fungal OTUs. Understanding the rhizosphere soil microbial diversity characteristics of the areca palm could provide the basis for exploring microbial association and maintaining the areca palm’s health.

Community RNA-Seq: multi-kingdom responses to living versus decaying roots in soil

Roots are a primary source of organic carbon input in most soils. The consumption of living and detrital root inputs involves multi-trophic processes and multiple kingdoms of microbial life, but typical microbial ecology studies focus on only one or two major lineages. We used Illumina shotgun RNA sequencing to conduct PCR-independent SSU rRNA community analysis ("community RNA-Seq") and simultaneously assess the bacteria, archaea, fungi, and microfauna surrounding both living and decomposing roots of the annual grass, Avena fatua. Plants were grown in 13CO2-labeled microcosms amended with 15N-root litter to identify the preferences of rhizosphere organisms for root exudates (13C) versus decaying root biomass (15N) using NanoSIMS microarray imaging (Chip-SIP). When litter was available, rhizosphere and bulk soil had significantly more Amoebozoa, which are potentially important yet often overlooked top-down drivers of detritusphere community dynamics and nutrient cycling. Bulk soil containing litter was depleted in Actinobacteria but had significantly more Bacteroidetes and Proteobacteria. While Actinobacteria were abundant in the rhizosphere, Chip-SIP showed Actinobacteria preferentially incorporated litter relative to root exudates, indicating this group's more prominent role in detritus elemental cycling in the rhizosphere. Our results emphasize that decomposition is a multi-trophic process involving complex interactions, and our methodology can be used to track the trajectory of carbon through multi-kingdom soil food webs.

A roadmap towards predicting species interaction networks (across space and time)

Networks of species interactions underpin numerous ecosystem processes, but comprehensively sampling these interactions is difficult. Interactions intrinsically vary across space and time, and given the number of species that compose ecological communities, it can be tough to distinguish between a true negative (where two species never interact) from a false negative (where two species have not been observed interacting even though they actually do). Assessing the likelihood of interactions between species is an imperative for several fields of ecology. This means that to predict interactions between species-and to describe the structure, variation, and change of the ecological networks they form-we need to rely on modelling tools. Here, we provide a proof-of-concept, where we show how a simple neural network model makes accurate predictions about species interactions given limited data. We then assess the challenges and opportunities associated with improving interaction predictions, and provide a conceptual roadmap forward towards predictive models of ecological networks that is explicitly spatial and temporal. We conclude with a brief primer on the relevant methods and tools needed to start building these models, which we hope will guide this research programme forward. This article is part of the theme issue 'Infectious disease macroecology: parasite diversity and dynamics across the globe'.

Local stability properties of complex, species-rich soil food webs with functional block structure

Ecologists have long debated the properties that confer stability to complex, species-rich ecological networks. Species-level soil food webs are large and structured networks of central importance to ecosystem functioning. Here, we conducted an analysis of the stability properties of an up-to-date set of theoretical soil food web models that account both for realistic levels of species richness and the most recent views on the topological structure (who is connected to whom) of these food webs. The stability of the network was best explained by two factors: strong correlations between interaction strengths and the blocked, nonrandom trophic structure of the web. These two factors could stabilize our model food webs even at the high levels of species richness that are typically found in soil, and that would make random systems very unstable. Also, the stability of our soil food webs is well-approximated by the cascade model. This result suggests that stability could emerge from the hierarchical structure of the functional organization of the web. Our study shows that under the assumption of equilibrium and small perturbations, theoretical soil food webs possess a topological structure that allows them to be complex yet more locally stable than their random counterpart. In particular, results strongly support the general hypothesis that the stability of rich and complex soil food webs is mostly driven by correlations in interaction strength and the organization of the soil food web into functional groups. The implication is that in real-world food web, any force disrupting the functional structure and distribution pattern of interaction strengths (i.e., energy fluxes) of the soil food webs will destabilize the dynamics of the system, leading to species extinction and major changes in the relative abundances of species.

Multiple phenotypic traits as triggers of host attacks towards ant symbionts: body size, morphological gestalt, and chemical mimicry accuracy

Background

Ant colonies are plagued by a diversity of arthropod guests, which adopt various strategies to avoid or to withstand host attacks. Chemical mimicry of host recognition cues is, for example, a common integration strategy of ant guests. The morphological gestalt and body size of ant guests have long been argued to also affect host hostility, but quantitative studies testing these predictions are largely missing. We here evaluated three guest traits as triggers of host aggression—body size, morphological gestalt, and accuracy in chemical mimicry—in a community of six Eciton army ant species and 29 guest species. We quantified ant aggression towards 314 guests in behavioral assays and, for the same individuals, determined their body size and their accuracy in mimicking ant cuticular hydrocarbon (CHC) profiles. We classified guests into the following gestalts: protective, myrmecoid, staphylinid-like, phorid-like, and larval-shaped. We expected that (1) guests with lower CHC mimicry accuracy are more frequently attacked; (2) larger guests are more frequently attacked; (3) guests of different morphological gestalt receive differing host aggression levels.

Results

Army ant species had distinct CHC profiles and accuracy of mimicking these profiles was variable among guests, with many species showing high mimicry accuracy. Unexpectedly, we did not find a clear relationship between chemical host similarity and host aggression, suggesting that other symbiont traits need to be considered. We detected a relationship between the guests’ body size and the received host aggression, in that diminutive forms were rarely attacked. Our data also indicated that morphological gestalt might be a valuable predictor of host aggression. While most ant-guest encounters remained peaceful, host behavior still differed towards guests in that ant aggression was primarily directed towards those guests possessing a protective or a staphylinid-like gestalt.

Conclusion

We demonstrate that CHC mimicry accuracy does not necessarily predict host aggression towards ant symbionts. Exploitation mechanisms are diverse, and we conclude that, besides chemical mimicry, other factors such as the guests’ morphological gestalt and especially their body size might be important, yet underrated traits shaping the level of host hostility against social insect symbionts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00427-8.

Cascading effects of moth outbreaks on subarctic soil food webs

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

Sulfur deposition changed the community structure of soil nematodes by affecting omnivores-predators

Nematodes generally occupy multiple trophic levels in detrital food webs, which play a vital role in energy flow, material conversion and nematodes community structure stability in the underground ecosystem. Sulfur (S) is one of the important soil nutrients, and it plays an important role in the nutrient cycle of grassland ecosystem. However, the impacts of S on soil fauna and subsurface detrital food webs in grassland ecosystems were rarely studied. Accordingly, to investigate the effects of sulfur deposition on soil nematodes and detrital food webs, we conducted a S addition experiment with distinct intensities from 0 to 50 g S m^-2 yr^-1 (S 0, S 1, S 2, S 5, S 10, S 15, S 20, and S 50) to simulated sulfur deposition in a meadow steppe of northern China. We documented a significant effect of S addition on the diversity and richness of nematodes, and the species richness of soil nematodes was high in the study site. But S addition had no significant effect on the total abundance and dominant species of nematodes (Cervidellus and Aphelenchus). Results of correlation analysis and structural equation modeling consistently indicated that omnivores-predators were significantly affected by sulfur addition. A significant increase in the Structural Index (which indicates food web structure) suggested increased top-down forces and changed community structure, although bacterivores, fungivores, plant parasites did not significantly. The present results suggest that sulfur deposition would change the composition of nematode community, affect the stability of nematode community structure, and increase the disturbance to the underground ecosystem. The study provides that the detailed information of the response of nematode to S deposition can be used to analyze the process of global change affecting the underground ecosystem.

Energetic constraints imposed on trophic interaction strengths enhance resilience in empirical and model food webs

Food web stability and resilience are at the heart of understanding the structure and functioning of ecosystems. Previous studies show that models of empirical food webs are substantially more stable than random ones, due to a few strong interactions embedded in a majority of weak interactions. Analyses of trophic interaction loops show that in empirical food webs the patterns of the interaction strengths prevent the occurrence of destabilizing heavy loops and thereby enhances resilience. Yet, it is still unexplored which biological mechanisms cause these patterns that enhance food web resilience. We quantified food web resilience using the real part of the maximum eigenvalue of the Jacobian matrix of the food web from a seagrass bed in the Yellow River Delta (YRD) wetland, that could be parametrized by the empirical data of the food web. We found that the empirically based Jacobian matrix of the YRD food web indicated a much higher resilience than random matrices with the same element values but arranged in random ways. Investigating the trophic interaction loops revealed that the high resilience was due to a negative correlation between the negative and positive interaction strengths (per capita top-down and bottom-up effects, respectively) within positive feedback loops with three species. The negative correlation showed that when the negative interaction strengths were strong the positive was weak, and vice versa. Our invented reformulation of loop weight in terms of biomasses and specific production rates showed that energetic properties of the trophic groups in the loop and mass-balance constraints, for example, the food uptake has to balance all losses, created the negative correlation between the interaction strengths. This result could be generalized using a dynamic intraguild predation model, which delivered the same pattern for a wide range of model parameters. Our results shed light on how energetic constraints at the trophic group and food web level create a pattern of interaction strengths within trophic interaction loops that enhances food web resilience.

Nematodes as Drivers of Plant Performance in Natural Systems

Nematodes form an important part of soil biodiversity as the most abundant and functionally diverse animals affecting plant performance. Most studies on plant-nematode interactions are focused on agriculture, while plant-nematode interactions in nature are less known. Here we highlight that nematodes can contribute to vegetation dynamics through direct negative effects on plants, and indirect positive effects through top-down predation on plant-associated organisms. Global change alters these interactions, of which better understanding is rapidly needed to better predict functional consequences. By expanding the knowledge of plant-nematode interactions in natural systems, an increase in basic understanding of key ecological topics such as plant-soil interactions and plant invasion dynamics will be obtained, while also increasing the insights and potential biotic repertoire to be applicable in sustainable plant management.

Seasonal Stabilities of Soil Nematode Communities and Their Relationships with Environmental Factors in Different Temperate Forest Types on the Chinese Loess Plateau

The bottom-up effects of vegetation have been documented to be strong drivers of the soil food web structure and functioning in temperate forests. However, how the forest type affects the stability of the soil food web is not well known. In the Ziwuling forest region of the Loess Plateau, we selected three typical forests, Pinus tabuliformis Carrière (PT), Betula platyphylla Sukaczev (BP), and Quercus liaotungensis Koidz. (QL), to investigate the soil nematode community characteristics in the dry (April) and rainy (August) season, and analyzed their relationships with the soil properties. The results showed that the characteristics of the soil nematode communities and their seasonal variations differed markedly among the forest types. Compared to P. tabuliformis (PT), the B. platyphylla (BP) and Q. liaotungensis (QL) forests had higher plant diversity and more easily decomposed litters, which were more effective for improving the soil resource availability, thus, leading to more beneficial effects on the soil nematode community. In both the dry and rainy season, the soil nematode abundance was the highest in the BP forest. The Shannon–Wiener diversity index (H’), Pielou’s evenness index (J’), and nematode channel ratio index (NCR) were higher, while the Simpson dominance index (λ) and plant parasite index (PPI) were lower, in the BP and QL forests compared with in the PT forest. From the dry to rainy season, the total nematode abundance and the abundance of fungivores, bacterivores, and omnivore-predators, significantly increased in the QL and PT forests, and the values of the Wasilewska index (WI), maturity index (MI), H’, J’, λ, and NCR showed the most significant seasonal variability in the PT forest, which were mainly driven by changes in the soil labile C and N and the moisture content between the two seasons. Generally, the seasonal stability of the soil nematode communities was the highest in the BP forest and the poorest in the PT forest, probably due to variations in the plant diversity. Our results suggest the importance of tree species and diversity as bottom-up regulating factors of the soil food web structure, function, and seasonal stability, which has important implications for sustainable forest management in the Loess Plateau and other temperate regions.

Influence of sea ice dynamics on population energetics of Western Hudson Bay polar bears

The Arctic marine ecosystem has experienced extensive changes in sea ice dynamics, with significant effects on ice-dependent species such as polar bears (Ursus maritimus). We used annual estimates of the numbers of bears onshore in the core summering area, age/sex structure and body condition data to estimate population energy density and storage energy in Western Hudson Bay polar bears from 1985 to 2018. We examined intra-population variation in energetic patterns, temporal energetic trends and the relationship between population energetics and sea ice conditions. Energy metrics for most demographic classes declined over time in relation to earlier sea ice breakup, most significantly for solitary adult females and subadult males, suggesting their greater vulnerability to nutritional stress than other age/sex classes. Temporal declines in population energy metrics were related to earlier breakup and longer lagged open-water periods, suggesting multi-year effects of sea ice decline. The length of the open-water period ranged from 102 to 166 days and increased significantly by 9.9 days/decade over the study period. Total population energy density and storage energy were significantly lower when sea ice breakup occurred earlier and the lagged open-water period was longer. At the earliest breakup and a lagged open-water period of 180 days, population energy density was predicted to be 33% lower than our minimum estimated energy density and population storage energy was predicted to be 40% lower than the minimum estimated storage energy. Consequently, over the study, the total population energy density declined by 53% (mean: 3668 ± 386 MJ kg-1/decade) and total population storage energy declined by 56% (mean: 435900 ± 46770 MJ/decade). This study provides insights into ecological mechanisms linking population responses to sea ice decline and highlights the significance of maintaining long-term research programs.

Stability of ecosystems enhanced by species-interaction constraints

Ecosystem stability is a central question both in theoretical and applied biology. Dynamical systems theory can be used to analyze how growth rates, carrying capacities, and patterns of species interactions affect the stability of an ecosystem. The response to increasing complexity has been extensively studied and the general conclusion is that there is a limit. While there is a complexity limit to stability at which global destabilisation occurs, the collapse rarely happens suddenly if a system is fully viable (no species is extinct). In fact, when complexity is successively increased, we find that the generic response is to go through multiple single-species extinctions before a global collapse. In this paper we demonstrate this finding via both numerical simulations and elaborations of theoretical predictions. We explore more biological interaction patterns, and, perhaps most importantly, we show that constrained interaction structures-a constant row sum in the interaction matrix-prevent extinctions from occurring. This makes an ecosystem more robust in terms of allowed complexity, but it also means singles-species extinctions do not precede or signal collapse-a drastically different behavior compared to the generic and commonly assumed case. We further argue that this constrained interaction structure-limiting the total interactions for each species-is biologically plausible.

Seasonal food webs with migrations: multi-season models reveal indirect species interactions in the Canadian Arctic tundra

Models incorporating seasonality are necessary to fully assess the impact of global warming on Arctic communities. Seasonal migrations are a key component of Arctic food webs that still elude current theories predicting a single community equilibrium. We develop a multi-season model of predator-prey dynamics using a hybrid dynamical systems framework applied to a simplified tundra food web (lemming-fox-goose-owl). Hybrid systems models can accommodate multiple equilibria, which is a basic requirement for modelling food webs whose topology changes with season. We demonstrate that our model can generate multi-annual cycling in lemming dynamics, solely from a combined effect of seasonality and state-dependent behaviour. We compare our multi-season model to a static model of the predator-prey community dynamics and study the interactions between species. Interestingly, including seasonality reveals indirect interactions between migrants and residents not captured by the static model. Further, we find that the direction and magnitude of interactions between two species are not necessarily accurate using only summer time-series. Our study demonstrates the need for the development of multi-season models and provides the tools to analyse them. Integrating seasonality in food web modelling is a vital step to improve predictions about the impacts of climate change on ecosystem functioning. This article is part of the theme issue 'The changing Arctic Ocean: consequences for biological communities, biogeochemical processes and ecosystem functioning'.

Trophic interactions as determinants of the arbuscular mycorrhizal fungal community with cascading plant-promoting consequences

BACKGROUND

The soil mycobiome is composed of a complex and diverse fungal community, which includes functionally diverse species ranging from plant pathogens to mutualists. Among the latter are arbuscular mycorrhizal fungi (AMF) that provide phosphorous (P) to plants. While plant hosts and abiotic parameters are known to structure AMF communities, it remains largely unknown how higher trophic level organisms, including protists and nematodes, affect AMF abundance and community composition.

RESULTS

Here, we explored the connections between AMF, fungivorous protists and nematodes that could partly reflect trophic interactions, and linked those to rhizosphere P dynamics and plant performance in a long-term manure application setting. Our results revealed that manure addition increased AMF biomass and the density of fungivorous nematodes, and tailored the community structures of AMF, fungivorous protists, and nematodes. We detected a higher abundance of AMF digested by the dominant fungivorous nematodes Aphelenchoides and Aphelenchus in high manure treatments compared to no manure and low manure treatments. Structural equation modeling combined with network analysis suggested that predation by fungivorous protists and nematodes stimulated AMF biomass and modified the AMF community composition. The mycorrhizal-fungivore interactions catalyzed AMF colonization and expression levels of the P transporter gene ZMPht1;6 in maize roots, which resulted in enhanced plant productivity.

CONCLUSIONS

Our study highlights the importance of predation as a key element in shaping the composition and enhancing the biomass of AMF, leading to increased plant performance. As such, we clarify novel biological mechanism of the complex interactions between AMF, fungivorous protists, and nematodes in driving P absorption and plant performance. Video Abstract.

Comparison of bacterial communities in soil samples with and without tomato bacterial wilt caused by Ralstonia solanacearum species complex

BACKGROUND

Ralstonia solanacearum is one of the most notorious soil-borne phytopathogens. It causes a severe wilt disease with deadly effects on many economically important crops. The microbita of disease-suppressive soils are thought that they can contribute to the disease resistance of crop plants, thus, evaluation of the microbial community and their interaction characteristics between suppressive soil (SS) and conducive soil (CS) will help to understand resistance mechanism. To do this, the bacterial community structure, correlation analysis with soil chemical properties, interaction network of SS (nearly no disease in three years), and CS (suffered heavy bacterial wilt disease) were analyzed.

RESULTS

A higher bacterial community diversity index was found in SS, the relative abundance of Nocardioides, Gaiella and norank_f_Anaerolineaceae were significantly more than that of the CS. Moreover, the relative abundance of main genera Bacillus, norank_o_Gaiellales, Roseiflexus, and norank_o_Gemmatimonadaceae were significantly more than that of the CS. Redundancy analysis at the genus level indicated that the available phosphate played a key role in the bacterial community distribution, and its role was negatively correlated with soil pH, organic matter content, alkali-hydrolyzable nitrogen, and available potassium contents. Interaction network analysis further demonstrated that greater diversity at the genus level existed in the SS network and formed a stable network. Additionally, the species of Mycobacterium, Cyanobacteria, and Rhodobiaceae are the key components that sustain the network stability. Seven clusters of orthologous groups exhibited significant differences between SS and CS. Moreover, 55 bacterial strains with distinct antagonistic activities to R. solancearum were isolated and identified from the healthy tomato plant rhizosphere soil of the CS.

CONCLUSIONS

Our findings indicate that the bacterial diversity and interaction network differed between the CS and SS samples, providing a good foundation in the study of bacterial wilt.

From winter to summer and back: Lessons from the parameterization of a seasonal food web model for the Białowieża forest

Dynamic food web models describe how species abundances change over time as a function of trophic and life-history parameters. They are essential to predicting the response of ecosystems to perturbations. However, they are notoriously difficult to parameterize, so that most models rely heavily either on allometric scaling of parameters or inverse estimation of biomass flows. The allometric approach makes species of comparable body mass have near-identical parameters which can generate extinctions within a trophic level. The biomass flow approach is more precise, but is restricted to steady-states, which is not appropriate for time-varying environments. Adequately parameterizing large food webs of temperate and arctic environments requires dealing both with many species of similar sizes and a strongly seasonal environment. Inspired by the rich empirical knowledge on the vertebrate food web of the Białowieża forest, we parameterize a bipartite food web model comprising 21 predators and 124 prey species. Our model is a non-autonomous coupled ordinary differential equations system that allows for seasonality in life-history and predation parameters. Birth and death rates, seasonal descriptions of diet for each species, food requirements and biomass information are combined into a seasonal parameterization of a dynamic food web model. Food web seasonality is implemented with time-varying intrinsic growth rate and interaction parameters, while predation is modelled with both type I and type II functional responses. All our model variants allow for >80% persistence in spite of massive apparent competition, and a quantitative match to observed (seasonal) biomasses. We also identify trade-offs between maximizing persistence, reproducing observed biomasses, and ensuring model robustness to sampling errors. Although multi-annual cycles are expected with type II functional responses, they are here prevented by a strong predator self-regulation. We discuss these results and possible improvements on the model. We provide a general workflow to parameterize dynamic food web models in seasonal environments, based on a real case study. This may help to better predict how biodiverse food webs respond to changing environments.

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.

Patterns of local, intercontinental and interseasonal variation of soil bacterial and eukaryotic microbial communities

Although ongoing research has revealed some of the main drivers behind global spatial patterns of microbial communities, spatio-temporal dynamics of these communities still remain largely unexplored. Here, we investigate spatio-temporal variability of both bacterial and eukaryotic soil microbial communities at local and intercontinental scales. We compare how temporal variation in community composition scales with spatial variation in community composition, and explore the extent to which bacteria, protists, fungi and metazoa have similar patterns of temporal community dynamics. All soil microbial groups displayed a strong correlation between spatial distance and community dissimilarity, which was related to the ratio of organism to sample size. Temporal changes were variable, ranging from equal to local between-sample variation, to as large as that between communities several thousand kilometers apart. Moreover, significant correlations were found between bacterial and protist communities, as well as between protist and fungal communities, indicating that these microbial groups change in tandem, potentially driven by interactions between them. We conclude that temporal variation can be considerable in soil microbial communities, and that future studies need to consider temporal variation in order to reliably capture all drivers of soil microbiome changes.

Soil microbial community responses to climate extremes: resistance, resilience and transitions to alternative states

A major challenge for advancing our understanding of the functional role of soil microbial communities is to link changes in their structure and function under climate change. To address this challenge requires new understanding of the mechanisms that underlie the capacity of soil microbial communities to resist and recover from climate extremes. Here, we synthesize emerging understanding of the intrinsic and extrinsic factors that influence the resistance and resilience of soil microbial communities to climate extremes, with a focus on drought, and identify drivers that might trigger abrupt changes to alternative states. We highlight research challenges and propose a path for advancing our understanding of the resistance and resilience of soil microbial communities to climate extremes, and of their vulnerability to transitions to alternative states, including the use of trait-based approaches. We identify a need for new approaches to quantify resistance and resilience of soil microbial communities, and to identify thresholds for transitions to alternative states. We show how high-resolution time series coupled with gradient designs will enable detecting response patterns to interacting drivers. Finally, to account for extrinsic factors, we suggest that future studies should use environmental gradients to track soil microbial community responses to climate extremes in space and time. This article is part of the theme issue 'Climate change and ecosystems: threats, opportunities and solutions'.

Foreseeing the future of mutualistic communities beyond collapse

Changing conditions may lead to sudden shifts in the state of ecosystems when critical thresholds are passed. Some well-studied drivers of such transitions lead to predictable outcomes such as a turbid lake or a degraded landscape. Many ecosystems are, however, complex systems of many interacting species. While detecting upcoming transitions in such systems is challenging, predicting what comes after a critical transition is terra incognita altogether. The problem is that complex ecosystems may shift to many different, alternative states. Whether an impending transition has minor, positive or catastrophic effects is thus unclear. Some systems may, however, behave more predictably than others. The dynamics of mutualistic communities can be expected to be relatively simple, because delayed negative feedbacks leading to oscillatory or other complex dynamics are weak. Here, we address the question of whether this relative simplicity allows us to foresee a community's future state. As a case study, we use a model of a bipartite mutualistic network and show that a network's post-transition state is indicated by the way in which a system recovers from minor disturbances. Similar results obtained with a unipartite model of facilitation suggest that our results are of relevance to a wide range of mutualistic systems.

First study of food webs in a large glacial river: the trophic role of invasive trout

Abstract The aim of this study was to determine the food webs structure of a large Patagonian river in two river sections (Upstream and Midstream) and to evaluate isotopic overlap between native and introduced species. We used stable isotope analyses of δ15N and δ13C and stomach content. The Upstream section had a more complex food webs structure with a greater richness of macroinvertebrates and fish species than Midstream. Upstream basal resources were dominated by filamentous algae. Lake Trout were found to have a higher trophic position than all other fish species in that area although, the most abundant fish species, were Rainbow Trout. Depending on the life stage, Rainbow Trout shifted from prey to competitor/predator. In the Midstream section, the base of the food webs was dominated by coarse particulate organic matter, and adult Rainbow Trout had the highest trophic level. Isotopic values changed among macroinvertebrates and fish for both areas. The two most abundant native and invasive species — Puyen and Rainbow Trout — showed an isotopic separation in Midstream but did not in Upstream areas. The presence of invasive fish that occupy top trophic levels can have a significant impact on native fish populations that have great ecological importance in the region.

Metabarcoding data allow for reliable biomass estimates in the most abundant animals on earth

Microscopic organisms are the dominant and most diverse organisms on Earth. Nematodes, as part of this microscopic diversity, are by far the most abundant animals and their diversity is equally high. Molecular metabarcoding is often applied to study the diversity of microorganisms, but has yet to become the standard to determine nematode communities. As such, the information metabarcoding provides, such as in terms of species coverage, taxonomic resolution and especially if sequence reads can be linked to the abundance or biomass of nematodes in a sample, has yet to be determined. Here, we applied metabarcoding using three primer sets located within ribosomal rRNA gene regions to target assembled mock-communities consisting of 18 different nematode species that we established in 9 different compositions. We determined abundances and biomass of all species added to examine if relative sequence abundance or biomass can be linked to relative sequence reads. We found that nematode communities are not equally represented by the three different primer sets and we found that relative read abundances almost perfectly correlated positively with relative species biomass for two of the primer sets. This strong biomass-read number correlation suggests that metabarcoding reads can reveal biomass information even amongst more complex nematode communities as present in the environment and possibly can be transferred to better study other groups of organisms. This biomass-read link is of particular importance for more reliably assessing nutrient flow through food-webs, as well as adjusting biogeochemical models through user-friendly and easily obtainable metabarcoding data.

Nonlinear population dynamics are ubiquitous in animals

Nonlinear dynamics, where a change in the input is not proportional to a change in the output, are often found throughout nature, for example in biochemical kinetics. Because of the complex suite of interacting abiotic and biotic variables present in ecosystems, animal population dynamics are often thought to be driven in a nonlinear, state-dependent fashion. However, so far these have only been identified in model organisms and some natural systems. Here we show that nonlinear population dynamics are ubiquitous in nature. We use nonlinear forecasting to analyse 747 datasets of 228 species to find that insect population trends were highly nonlinear (74%), followed by mammals (58%), bony fish (49%) and birds (35%). This indicates that linear, equilibrium-based model assumptions may fail at predicting population dynamics across a wide range of animal taxa. We show that faster-reproducing animals are more likely to have nonlinear and high-dimensional dynamics, supporting past ecological theory. Lastly, only a third of time series were predictable beyond two years; therefore, the ability to predict animal population trends using these methods may be limited. Our results suggest that the complex dynamics necessary to cause regime shifts and other transitions may be inherent in a wide variety of animals.

Climate change and land use induce functional shifts in soil nematode communities

Land-use intensification represents one major threat to the diversity and functioning of terrestrial ecosystems. In the face of concurrent climate change, concerns are growing about the ability of intensively managed agroecosystems to ensure stable food provisioning, as they may be particularly vulnerable to climate extreme-induced harvest losses and pest outbreaks. Extensively managed systems, in contrast, were shown to mitigate climate change based on plant diversity-mediated effects, such as higher functional redundancy or asynchrony of species. In this context, the maintenance of soils is essential to sustain key ecosystem functions such as nutrient cycling, pest control, and crop yield. Within the highly diverse soil fauna, nematodes represent an important group as their trophic spectrum ranges from detritivores to predators and they allow inferences to the overall state of the ecosystem (bioindicators). Here, we investigated the effects of simulated climate change and land-use intensity on the diversity and abundance of soil nematode functional groups and functional indices in two consecutive years. We revealed that especially land use induced complex shifts in the nematode community with strong seasonal dynamics, while future climate led to weaker effects. Strikingly, the high nematode densities associated with altered climatic conditions and intensive land use were a consequence of increased densities of opportunists and potential pest species (i.e., plant feeders). This coincided with a less diverse and less structured community with presumably reduced capabilities to withstand environmental stress. These degraded soil food web conditions represent a potential threat to ecosystem functioning and underline the importance of management practices that preserve belowground organisms.

Non-trophic interactions strengthen the diversity-functioning relationship in an ecological bioenergetic network model

Ecological communities are undeniably diverse, both in terms of the species that compose them as well as the type of interactions that link species to each other. Despite this long recognition of the coexistence of multiple interaction types in nature, little is known about the consequences of this diversity for community functioning. In the ongoing context of global change and increasing species extinction rates, it seems crucial to improve our understanding of the drivers of the relationship between species diversity and ecosystem functioning. Here, using a multispecies dynamical model of ecological communities including various interaction types (e.g. competition for space, predator interference, recruitment facilitation in addition to feeding), we studied the role of the presence and the intensity of these interactions for species diversity, community functioning (biomass and production) and the relationship between diversity and functioning.Taken jointly, the diverse interactions have significant effects on species diversity, whose amplitude and sign depend on the type of interactions involved and their relative abundance. They however consistently increase the slope of the relationship between diversity and functioning, suggesting that species losses might have stronger effects on community functioning than expected when ignoring the diversity of interaction types and focusing on feeding interactions only.

The question of how species diversity contributes to the functioning of ecological communities has intrigued ecologists for decades, and is especially relevant in the current context of species extinctions. Ecological communities are not only diverse in terms of the species that compose them but also in terms of the way they interact with each other: for example, species compete for space and for food, eat and facilitate each other. The diversity of ways species interact has rarely been taken into account in the study of ecological communities, although widely acknowledged. Here we show that the diversity of interaction types matters: it affects species diversity, community functioning and the relationship between them by strengthening this relationship. This means that when the diversity of interaction types is taken into account, species losses have stronger impacts on the functioning of ecological communities. Our results therefore suggest that species loss may have more important consequences than expected based on classical models that do not take the diversity of interaction types into account.

Conversion of rainforest to oil palm and rubber plantations alters energy channels in soil food webs

In the last decades, lowland tropical rainforest has been converted in large into plantation systems. Despite the evident changes above ground, the effect of rainforest conversion on the channeling of energy in soil food webs was not studied. Here, we investigated community-level neutral lipid fatty acid profiles in dominant soil fauna to track energy channels in rainforest, rubber, and oil palm plantations in Sumatra, Indonesia. Abundant macrofauna including Araneae, Chilopoda, and Diplopoda contained high amounts of plant and fungal biomarker fatty acids (FAs). Lumbricina had the lowest amount of plant, but the highest amount of animal-synthesized C20 polyunsaturated FAs as compared to other soil taxa. Mesofauna detritivores (Collembola and Oribatida) contained high amounts of algal biomarker FAs. The differences in FA profiles between taxa were evident if data were analyzed across land-use systems, suggesting that soil fauna of different size (macro- and mesofauna) are associated with different energy channels. Despite that, rainforest conversion changed the biomarker FA composition of soil fauna at the community level. Conversion of rainforest into oil palm plantations enhanced the plant energy channel in soil food webs and reduced the bacterial energy channel; conversion into rubber plantations reduced the AMF-based energy channel. The changes in energy distribution within soil food webs may have significant implications for the functioning of tropical ecosystems and their response to environmental changes. At present, these responses are hard to predict considering the poor knowledge on structure and functioning of tropical soil food webs.

Latitudinal variation in soil nematode communities under climate warming-related range-expanding and native plants

Current climate change has led to latitudinal and altitudinal range expansions of numerous species. During such range expansions, plant species are expected to experience changes in interactions with other organisms, especially with belowground biota that have a limited dispersal capacity. Nematodes form a key component of the belowground food web as they include bacterivores, fungivores, omnivores and root herbivores. However, their community composition under climate change-driven intracontinental range-expanding plants has been studied almost exclusively under controlled conditions, whereas little is known about actual patterns in the field. Here, we use novel molecular sequencing techniques combined with morphological quantification in order to examine nematode communities in the rhizospheres of four range-expanding and four congeneric native species along a 2,000 km latitudinal transect from South-Eastern to North-Western Europe. We tested the hypotheses that latitudinal shifts in nematode community composition are stronger in range-expanding plant species than in congeneric natives and that in their new range, range-expanding plant species accumulate fewest root-feeding nematodes. Our results show latitudinal variation in nematode community composition of both range expanders and native plant species, while operational taxonomic unit richness remained the same across ranges. Therefore, range-expanding plant species face different nematode communities at higher latitudes, but this is also the case for widespread native plant species. Only one of the four range-expanding plant species showed a stronger shift in nematode community composition than its congeneric native and accumulated fewer root-feeding nematodes in its new range. We conclude that variation in nematode community composition with increasing latitude occurs for both range-expanding and native plant species and that some range-expanding plant species may become released from root-feeding nematodes in the new range.

Food-web interaction strength distributions are conserved by greater variation between than within predator-prey pairs

Species interactions in food webs are usually recognized as dynamic, varying across species, space, and time because of biotic and abiotic drivers. Yet food webs also show emergent properties that appear consistent, such as a skewed frequency distribution of interaction strengths (many weak, few strong). Reconciling these two properties requires an understanding of the variation in pairwise interaction strengths and its underlying mechanisms. We estimated stream sculpin feeding rates in three seasons at nine sites in Oregon to examine variation in trophic interaction strengths both across and within predator-prey pairs. Predator and prey densities, prey body mass, and abiotic factors were considered as putative drivers of within-pair variation over space and time. We hypothesized that consistently skewed interaction strength distributions could result if individual interaction strengths show relatively little variation, or alternatively, if interaction strengths vary but shift in ways that conserve their overall frequency distribution. Feeding rate distributions remained consistently and positively skewed across all sites and seasons. The mean coefficient of variation in feeding rates within each of 25 focal species pairs across surveys was less than half the mean coefficient of variation seen across species pairs within a survey. The rank order of feeding rates also remained conserved across streams, seasons and individual surveys. On average, feeding rates on each prey taxon nonetheless varied by a hundredfold, with some feeding rates showing more variation in space and others in time. In general, feeding rates increased with prey density and decreased with high stream flows and low water temperatures, although for nearly half of all species pairs, factors other than prey density explained the most variation. Our findings show that although individual interaction strengths exhibit considerable variation in space and time, they can nonetheless remain relatively consistent, and thus predictable, compared to the even larger variation that occurs across species pairs. These results highlight how the ecological scale of inference can strongly shape conclusions about interaction strength consistency and help reconcile how the skewed nature of interaction strength distributions can persist in highly dynamic food webs.

Microbial Network and Soil Properties Are Changed in Bacterial Wilt-Susceptible Soil

Bacterial wilt disease is a devastating disease of crops, which leads to huge economic loss worldwide. It is hypothesized that the occurrence of bacterial wilt may be related to changes in soil chemical properties and microbial interactions. In this study, we compared the soil chemical properties and microbial network structures of a healthy soil (HS) and a bacterial wilt-susceptible soil (BWS). The contents of available nitrogen, potassium, and phosphorus and the soil pH in the BWS were significantly lower than those in the HS. BWS showed nutrient deficiency and acidification in comparison with the HS. The structure and composition of the BWS network were quite different from those of the HS network. The BWS network had fewer modules and edges and lower connectivity than the HS network. The HS network contained more interacting species, more key microorganisms, and better high-order organization and thus was more complex and stable than the BWS network. Most nodes and module memberships were unshared by the two networks, while the ones that were shared showed different topological roles. Some generalists in the HS network became specialists in the BWS network, indicating that the topological roles of microbes were changed and key microorganisms were shifted in the BWS. In summary, the composition and structure of the microbial network of the BWS were different from that of the HS. Many microbial network connections were missing in the BWS, which most likely provided conditions leading to higher rates of bacterial wilt disease.IMPORTANCE Bacterial wilt disease is caused by the pathogen Ralstonia solanacearum and is a widespread devastating soilborne disease leading to huge economic losses worldwide. The soil microbial community is crucial to the capacity of soils to suppress soilborne diseases through complex interactions. Network analysis can effectively explore these complex interactions. In this study, we used a random matrix theory (RMT)-based network approach to investigate the changes in microbial network and associated microbial interactions in a bacterial wilt-susceptible soil (BWS) in comparison to a healthy soil (HS). We found that the structure and composition of the microbial network in BWSs were quite different from those of the HS. The BWS network had fewer modules, edges, and key microorganisms and lower connectivity than the HS network. In the BWSs, apparently the topological role of microbes was changed and key microorganisms were shifted to specialists.

Competition and Predation in Soil Fungivorous Microarthropods Using Stable Isotope Ratio Mass Spectrometry

The soil food web is often described as having three main energy channels: root, bacterial and fungal. Here we provide quantitative data using a sensitive stable isotope ratio mass spectrometry procedure with microcosms on species interactions in the fungal pathway. We measured 15N and 13C enrichment in microarthropods through grazing rare isotope enriched fungal mycelia. Experimental treatments were various combinations of 1, 2, 3, 4 microarthropods species. We used three fungivores (the collembolan Lepidocyrtus curvicollis, the Astigmata Tyrophagus putrescentiae, the Oribatida Oribatula tibialis), and the Mesostigmata predator Hypoaspis acquilifer. We collected individuals of each species separately, as well as their feces, and molt where available. All three fungivorous microarthropods consumed significantly more than their own body weight per day. The three fungivores differed in their consumption of the mycelium as it was not equally palatable to each. The Mesostigmata predator Hypoaspis also differed in its microarthropod prey preference. In multiple species combinations microarthropod behavioral interactions modified consumption and predation rates. Our selection of mites of different sizes, with varied preference for the mycelium, combined with differing predation rates on each mite, demonstrate that even three trophic level interactions with only five interacting species are not predictably simple. The interpretation of the stable isotope results and consumed-excreted weights indicate that: (a) behavior and microscopic observations should not be ignored in competition-predation interactions, and (b) functional guilds can take advantage of more diverse food opportunities. The reality of mixed diets complicates functional guild assignments that are reflected in 15N and 13C isotope levels at natural abundances in the environment. Microcosm experiments with this sensitive technique can help decipher the interpretation of rare isotope natural abundance values, as well as providing measured consumption, growth, and excretion rate values for modeling soil food web interactions.

Using Network Theory to Understand and Predict Biological Invasions

Understanding and predicting biological invasions is challenging because of the complexity of many interacting players. A holistic approach is needed with the potential to simultaneously consider all relevant effects and effectors. Using networks to describe the relevant anthropogenic and ecological factors, from community-level to global scales, promises advances in understanding aspects of invasion from propagule pressure, through establishment, spread, and ecological impact of invaders. These insights could lead to development of new tools for prevention and management of invasions that are based on species' network characteristics and use of networks to predict the ecological effects of invaders. Here, we review the findings from network ecology that show the most promise for invasion biology and identify pressing needs for future research.

Horizontal and vertical diversity jointly shape food web stability against small and large perturbations

The biodiversity of food webs is composed of horizontal (i.e. within trophic levels) and vertical diversity (i.e. the number of trophic levels). Understanding their joint effect on stability is a key challenge. Theory mostly considers their individual effects and focuses on small perturbations near equilibrium in hypothetical food webs. Here, we study the joint effects of horizontal and vertical diversity on the stability of hypothetical (modelled) and empirical food webs. In modelled food webs, horizontal and vertical diversity increased and decreased stability, respectively, with a stronger positive effect of producer diversity on stability at higher consumer diversity. Experiments with an empirical plankton food web, where we manipulated horizontal and vertical diversity and measured stability from species interactions and from resilience against large perturbations, confirmed these predictions. Taken together, our findings highlight the need to conserve horizontal biodiversity at different trophic levels to ensure stability.

Species-level predation network uncovers high prey specificity in a Neotropical army ant community

Army ants are among the top arthropod predators and considered keystone species in tropical ecosystems. During daily mass raids with many thousand workers, army ants hunt live prey, likely exerting strong top-down control on prey species. Many tropical sites exhibit a high army ant species diversity (>20 species), suggesting that sympatric species partition the available prey niches. However, whether and to what extent this is achieved has not been intensively studied yet. We therefore conducted a large-scale diet survey of a community of surface-raiding army ants at La Selva Biological Station in Costa Rica. We systematically collected 3,262 prey items from eleven army ant species (genera Eciton, Nomamyrmex and Neivamyrmex). Prey items were classified as ant prey or non-ant prey. The prey nearly exclusively consisted of other ants (98%), and most booty was ant brood (87%). Using morphological characters and DNA barcoding, we identified a total of 1,103 ant prey specimens to the species level. One hundred twenty-nine ant species were detected among the army ant prey, representing about 30% of the known local ant diversity. Using weighted bipartite network analyses, we show that prey specialization in army ants is unexpectedly high and prey niche overlap very small. Besides food niche differentiation, we uncovered a spatiotemporal niche differentiation in army ant raid activity. We discuss competition-driven multidimensional niche differentiation and predator-prey arms races as possible mechanisms underlying prey specialization in army ants. By combining systematic prey sampling with species-level prey identification and network analyses, our integrative approach can guide future research by portraying how predator-prey interactions in complex communities can be reliably studied, even in cases where morphological prey identification is infeasible.

Effects of plant functional group removal on structure and function of soil communities across contrasting ecosystems

Loss of plant diversity has an impact on ecosystems worldwide, but we lack a mechanistic understanding of how this loss may influence below-ground biota and ecosystem functions across contrasting ecosystems in the long term. We used the longest running biodiversity manipulation experiment across contrasting ecosystems in existence to explore the below-ground consequences of 19 years of plant functional group removals for each of 30 contrasting forested lake islands in northern Sweden. We found that, against expectations, the effects of plant removals on the communities of key groups of soil organisms (bacteria, fungi and nematodes), and organic matter quality and soil ecosystem functioning (decomposition and microbial activity) were relatively similar among islands that varied greatly in productivity and soil fertility. This highlights that, in contrast to what has been shown for plant productivity, plant biodiversity loss effects on below-ground functions can be relatively insensitive to environmental context or variation among widely contrasting ecosystems.

Exotic earthworms maintain soil biodiversity by altering bottom-up effects of plants on the composition of soil microbial groups and nematode communities

Bottom-up effects of plants on soil communities can be modified by the activity of exotic earthworms, by altering resource availability for soil food webs through feeding, burrowing, and casting activities. The present study explored effects of plants (planting of shrubs) on soil micro-food webs (composition of soil microbial and nematode communities), and whether these effects were altered by the activity of exotic earthworms (exotic earthworms addition). Planted shrubs resulted in a non-significant increase of bacterial biomass and significantly increased the abundance of different nematode trophic groups and total nematode biomass, indicating that planted shrubs had significant bottom-up effects on soil bacteria and nematodes. Planted shrubs decreased nematode diversity, evenness, and richness, but increased nematode dominance in the plots where the abundance of exotic earthworms was not amended. By contrast, these effects of shrub presence on soil biodiversity were not found in the plots that received exotic earthworms. In addition, planted shrubs increased the total energy flux to the nematode community. By contrast, the elevated activity of exotic earthworms mitigated the increase in total energy flux to nematodes in the presence of shrubs, and increased the ratio of fungal to bacterial PLFAs. Both of these changes indicate reduced energy flux in the plots with added exotic earthworms. Nematode diversity decreased, while nematode dominance increased with increasing total energy flux to nematodes, probably because few species benefited from high energy flux. Our study indicates that exotic earthworms can maintain soil biodiversity by reducing the energy flux through soil food webs.

Protist communities are more sensitive to nitrogen fertilization than other microorganisms in diverse agricultural soils

BACKGROUND

Agricultural food production is at the base of food and fodder, with fertilization having fundamentally and continuously increased crop yield over the last decades. The performance of crops is intimately tied to their microbiome as they together form holobionts. The importance of the microbiome for plant performance is, however, notoriously ignored in agricultural systems as fertilization disconnects the dependency of plants for often plant-beneficial microbial processes. Moreover, we lack a holistic understanding of how fertilization regimes affect the soil microbiome. Here, we examined the effect of a 2-year fertilization regime (no nitrogen fertilization control, nitrogen fertilization, and nitrogen fertilization plus straw amendment) on entire soil microbiomes (bacteria, fungi, and protist) in three common agricultural soil types cropped with maize in two seasons.

RESULTS

We found that the application of nitrogen fertilizers more strongly affected protist than bacterial and fungal communities. Nitrogen fertilization indirectly reduced protist diversity through changing abiotic properties and bacterial and fungal communities which differed between soil types and sampling seasons. Nitrogen fertilizer plus straw amendment had greater effects on soil physicochemical properties and microbiome diversity than nitrogen addition alone. Moreover, nitrogen fertilization, even more together with straw, increased soil microbiome network complexity, suggesting that the application of nitrogen fertilizers tightened soil microbiomes interactions.

CONCLUSIONS

Together, our results suggest that protists are the most susceptible microbiome component to the application of nitrogen fertilizers. As protist communities also exhibit the strongest seasonal dynamics, they serve as the most sensitive bioindicators of soil changes. Changes in protist communities might have long-term effects if some of the key protist hubs that govern microbiome complexities as top microbiome predators are altered. This study serves as the stepping stone to promote protists as promising agents in targeted microbiome engineering to help in reducing the dependency on exogenous unsustainably high fertilization and pesticide applications.

Plant-soil feedbacks promote coexistence and resilience in multi-species communities

Both ecological theory and empirical evidence suggest that negative frequency dependent feedbacks structure plant communities, but integration of these findings has been limited. Here we develop a generic model of frequency dependent feedback to analyze coexistence and invasibility in random theoretical and real communities for which frequency dependence through plant-soil feedbacks (PSFs) was determined empirically. We investigated community stability and invasibility by means of mechanistic analysis of invasion conditions and numerical simulations. We found that communities fall along a spectrum of coexistence types ranging from strict pair-wise negative feedback to strict intransitive networks. Intermediate community structures characterized by partial intransitivity may feature "keystone competitors" which disproportionately influence community stability. Real communities were characterized by stronger negative feedback and higher robustness to species loss than randomly assembled communities. Partial intransitivity became increasingly likely in more diverse communities. The results presented here theoretically explain why more diverse communities are characterized by stronger negative frequency dependent feedbacks, a pattern previously encountered in observational studies. Natural communities are more likely to be maintained by strict negative plant-soil feedback than expected by chance, but our results also show that community stability often depends on partial intransitivity. These results suggest that plant-soil feedbacks can facilitate coexistence in multi-species communities, but that these feedbacks may also initiate cascading effects on community diversity following from single-species loss.

Protists: Puppet Masters of the Rhizosphere Microbiome

The rhizosphere microbiome is a central determinant of plant performance. Microbiome assembly has traditionally been investigated from a bottom-up perspective, assessing how resources such as root exudates drive microbiome assembly. However, the importance of predation as a driver of microbiome structure has to date largely remained overlooked. Here we review the importance of protists, a paraphyletic group of unicellular eukaryotes, as a key regulator of microbiome assembly. Protists can promote plant-beneficial functions within the microbiome, accelerate nutrient cycling, and remove pathogens. We conclude that protists form an essential component of the rhizosphere microbiome and that accounting for predator-prey interactions would greatly improve our ability to predict and manage microbiome function at the service of plant growth and health.