Plant-mycorrhizal fungus co-occurrence network lacks substantial structure

The root microlandscape of arbuscular mycorrhizal fungi

Understanding the drivers of assemblages of arbuscular mycorrhizal fungi (AMF) is essential to leverage the benefits of AMF for plant growth and health. Arbuscular mycorrhizal fungi are heterogeneously distributed in space even at small scale. We review the role of plant distribution in driving AMF assemblages (the passenger hypothesis), using a transposition of the conceptual framework of landscape ecology. Because rooting systems correspond to habitat patches with limited carrying capacity that differ in quality due to host-preference effects, we suggest considering plant communities as mosaics of AMF microhabitats. We review how predictions from landscape ecology apply to plant community effects on AMF, and the existing evidence that tests these predictions. Although many studies have been conducted on the effect of plant compositional heterogeneity on AMF assemblages, they mostly focused on the effect of plant richness, while only a few investigated the effect of configurational heterogeneity, plant connectivity or plant community temporal dynamics. We propose key predictions and future prospects to fill these gaps. Considering plant communities as landscapes extends the passenger hypothesis by including a spatially explicit dimension and its associated ecological processes and may help understand and manipulate AMF assemblages at small spatial scales.

Organization of the macroinvertebrate community in a tropical annual agroecosystem into modules

The structure of macroinvertebrate communities in agroecosystems has been assumed to be modular and organized around key herbivore pests. We characterized the macroinvertebrate community in the annual organic brassica agroecosystem in tropical central Brazil to determine if the community was a random assemblage of independent populations or was organized into repeatable multi-species components. We sampled 36 macroinvertebrate taxa associated with six organic brassica farms at biweekly intervals during the dry season during two years in the Distrito Federal, Brazil. We used an unconstrained ordination based on latent variable modeling (boral) with negative binomial population counts to analyze community composition independent of variation in sample abundance. We evaluated observed community structure by comparing it with randomized alternatives. We found that the community was not a random assemblage and consistently organized itself into two modules based around the major herbivores; one with lepidoptera and whiteflies and their associated natural enemies which was gradually replaced during the season by one with brassica aphids, aphid parasitoids and coccinellids. This analysis suggests that the historical and present-day focus on pest herbivores and their associated species in agroecosystems may be justified based on community structure.

Stochastic processes dominate soil arbuscular mycorrhizal fungal community assembly along an elevation gradient in central Japan

Arbuscular mycorrhizal (AM) fungi play an important role in facilitating ecosystem function and stability. Yet, their community response patterns and ecological assembly processes along elevational gradients which cross a range of climates and soil conditions remain elusive. We used Illumina MiSeq sequencing to examine trends in soil AM fungal community along an elevational gradient from 100 m to 2300 m in central Japan. A total of 750 operational taxonomic units (OTUs) affiliated to 12 AM fungal genera were identified from soil samples, and the AM fungal community composition differed strongly with elevation, with variance explained more by climate, followed by soil and plant factors. The AM fungal α-diversity, network connectivity and complexity between AM fungal taxa and also with plant communities all exhibited a maximum at the mid-elevation of 800 m and then declined, principally influenced by soil pH and precipitation. Stochastic processes dominated AM fungal community assembly across the whole elevation gradient, with homogenizing dispersal being the main process. Only when AM fungal communities were contrasted across a relatively broad range of elevations, did variable selection (deterministic process) became significant, and even then in a mixed role with stochasticity. While OTUs of AM fungi are clearly adapted to particular environmental ranges, stochasticity due to rapid dispersal has a major role in determining their occurrence, suggesting that AM fungi may possess generalized and interchangeable niches, and can adjust their distribution rapidly - at least on the scale of a single mountain. This finding emphasizes that the roles of AM fungi in plant ecology may be non-specific and easily substituted, and furthermore that there is rapid local scale dispersal, which may allow plants to maintain effective AM associations under environmental change.

The Effects of Species Abundance, Spatial Distribution, and Phylogeny on a Plant-Ectomycorrhizal Fungal Network

Plant and root fungal interactions are among the most important belowground ecological interactions, however, the mechanisms underlying pairwise interactions and network patterns of rhizosphere fungi and host plants remain unknown. We tested whether neutral process or spatial constraints individually or jointly best explained quantitative plant-ectomycorrhizal fungal network assembly in a subtropical forest in southern China. Results showed that the observed plant-ectomycorrhizal fungal network had low connectivity, high interaction evenness, and an intermediate level of specialization, with nestedness and modularity both greater than random expectation. Incorporating information on the relative abundance and spatial overlap of plants and fungi well predicted network nestedness and connectance, but not necessarily explained other network metrics such as specificity. Spatial overlap better predicted pairwise species interactions of plants and ectomycorrhizal fungi than species abundance or a combination of species abundance and spatial overlap. There was a significant phylogenetic signal on species degree and interaction strength for ectomycorrhizal fungal but not for plant species. Our study suggests that neutral processes (species abundance matching) and niche/dispersal-related processes (implied by spatial overlap and phylogeny) jointly drive the shaping of a plant-ectomycorrhizal fungal network.

Impacts of Urbanization Undermine Nestedness of the Plant-Arbuscular Mycorrhizal Fungal Network

Cities are prone to ecological problems, yet the impacts of rapid global urbanization on the feedback between above- and belowground subsystems remain largely unknown. We sampled the roots of 8 common herbaceous plants within the Fifth Ring (urban areas) and in Jiufeng National Forest Park (rural areas) in Beijing (China) to assess the impacts of urbanization on the network of plant-arbuscular mycorrhizal (AM) fungal associations. Using Illumina MiSeq sequencing, 81 AM fungal OTUs were identified in 78 herb root samples. The Shannon, Simpson, and Pielou indices of root AM fungi in urban areas were significantly higher than those in rural areas. In this study, a significantly nested mycorrhizal association network was observed in rural areas (NODF = 64.68), whereas a non-nested pattern was observed in urban areas (NODF = 55.50). The competition index C-score (0.0769) of AM fungi in urban areas was slightly lower than that in rural areas (0.1431), and the species specialization (d’) of 8 host plants and fungal dissimilarity among 8 host plants in urban areas were significantly lower than those in rural areas. Convergent associations among hosts may be an important factor influencing this non-nested pattern of the plant-AM fungi network in urban areas. Generalists, rather than specialists, were enhanced during the establishment of mycorrhizal associations in urban areas. Our results suggest that reduced selectivity of host plants, and generalist promotion and specialist reduction of AM fungi during urbanization may contribute to the non-nested network of plant-AM fungal associations.

Defining non-indigenous fish assemblage types in Mediterranean rivers: Network analysis and management implications

The design and implementation of appropriate management actions to tackle the spread and negative impacts of non-indigenous fish species (NIFS) in freshwater ecosystems still remains a complex task. In an effort to address the limitations of current managerial approaches, our study develops and applies a classification framework to define non-indigenous fish assemblage types (FATs) in Mediterranean riverine ecosystems and identifies the linkage with various regional, local, biotic and abiotic environmental factors. This framework contributes to a pre-invasion stage screening and to the design of effective type-specific post-invasion management actions dealing with specific NIFS assemblages. A network analysis was applied in order to extract NIFS community level information from 393 samples obtained from 51 river basins of Greece, while a multivariate ordination analysis was conducted to detect the factors best explaining the structure and distribution of FATs. The association patterns of the study's inputs and outcomes are illustrated through an alluvial diagram, providing insights across different spatial scales. In total, five FATs were generated where the major NIFS contributors of average similarity within each modularity class revealed the key indicator species (Gambusia-FAT; Carassius/Lepomis-FAT; Pseudorasbora-FAT; Salmonids-FAT and Carp-FAT). Overall, the identified FATs varied spatially, indicating different community structure, mainly based on the diverse habitat preferences and life-history traits of indicator species. Alien FATs were mainly distributed within large and transboundary rivers, while Translocated and Salmonids FATs mostly occupied ecoregions with relatively depauperate faunas and often in biodiversity hotspots. The results of this study can identify conservation priorities within FATs, inform specific-type post-invasion management actions tackling NIFS, while in addition may provide valuable information for protecting high-priority water bodies before invasion.

A methodological framework to analyse determinants of host-microbiota networks, with an application to the relationships between Daphnia magna's gut microbiota and bacterioplankton

The past 30 years have seen both a surge of interest in assessing ecological interactions using tools borrowed from network theory and an explosion of data on the occurrence of microbial symbionts thanks to next-generation sequencing. Given that classic network methods cannot currently measure the respective effects of different environmental and biological drivers on network structure, we here present two methods to elucidate the determinants of bipartite interaction networks. The first method is based on classifications and compares communities within networks to the grouping of nodes by treatment or similar controlling groups. The second method assesses the link between multivariate explanatory variables and network structure using redundancy analyses after singular value decomposition. In both methods, the significance of effects can be gauged through two randomizations. Our methods were applied to experimental data on Daphnia magna and its interactions with gut microbiota and bacterioplankton. The whole network was affected by Daphnia's diet (algae and/or cyanobacteria) and sample type, but not by Daphnia genotype. At coarse grains, bacterioplankton and gut microbiota communities were different. At this scale, the structure of the gut microbiota-based network was not linked to any explanatory factors, while the bacterioplankton-based network was related to both Daphnia's diet and genotype. At finer grains, Daphnia's diet and genotype affected both microbial networks, but the effect of diet on gut microbiota network structure was mediated solely by differences in microbial richness. While no reciprocal effect between the microbial communities could be found, fine-grained analyses presented a more nuanced picture, with bacterioplankton likely affecting the composition of the gut microbiota. Our methods are widely applicable to bipartite networks, can elucidate both controlled and environmental effects in experimental setting using a large amount of sequencing data and can tease apart reciprocal effects of networks on one another. The twofold approach we propose has the advantage of being able to tease apart effects at different scales of network structure, thus allowing for detailed assessment of reciprocal effects of linked networks on one another. As such, our network methods can help ecologists understand huge datasets reporting microbial co-occurrences within different hosts.

Initial soil community drives heathland fungal community trajectory over multiple years through altered plant-soil interactions

Dispersal limitation, biotic interactions, and environmental filters interact to drive plant and fungal community assembly, but their combined effects are rarely investigated. This study examines how different heathland plant and fungal colonization scenarios realized via three biotic treatments - addition of mature heathland-derived sod, addition of hay, and no additions - affect soil fungal community development over 6 yr along a manipulated pH gradient in a large-scale experiment starting from an agricultural, topsoil removed state. Our results show that both biotic and abiotic (pH) treatments had a persistent influence on the development of fungal communities, but that sod additions diminished the effect of abiotic treatments through time. Analysis of correlation networks between soil fungi and plants suggests that the reduced effect of pH in the sod treatment, where both soil and plant propagules were added, might be due to plant-fungal interactions since the sod additions caused stronger, more specific, and more consistent connections compared with the no addition treatment. Based on these results, we suggest that the initial availability of heathland fungal and plant taxa, which reinforce each other, can significantly steer further fungal community development to an alternative configuration, overriding the otherwise prominent effect of abiotic (pH) conditions.

Species composition of arbuscular mycorrhizal communities changes with elevation in the Andes of South Ecuador

Arbuscular mycorrhizal fungi (AMF) are the most prominent mycobionts of plants in the tropics, yet little is known about their diversity, species compositions and factors driving AMF distribution patterns. To investigate whether elevation and associated vegetation type affect species composition, we sampled 646 mycorrhizal samples in locations between 1000 and 4000 m above sea level (masl) in the South of Ecuador. We estimated diversity, distribution and species compositions of AMF by cloning and Sanger sequencing the 18S rDNA (the section between AML1 and AML2) and subsequent derivation of fungal OTUs based on 99% sequence similarity. In addition, we analyzed the phylogenetic structure of the sites by computing the mean pairwise distance (MPD) and the mean nearest taxon difference (MNTD) for each elevation level. It revealed that AMF species compositions at 1000 and 2000 masl differ from 3000 and 4000 masl. Lower elevations (1000 and 2000 masl) were dominated by members of Glomeraceae, whereas Acaulosporaceae were more abundant in higher elevations (3000 and 4000 masl). Ordination of OTUs with respect to study sites revealed a correlation to elevation with a continuous turnover of species from lower to higher elevations. Most of the abundant OTUs are not endemic to South Ecuador. We also found a high proportion of rare OTUs at all elevations: 79-85% of OTUs occurred in less than 5% of the samples. Phylogenetic community analysis indicated clustering and evenness for most elevation levels indicating that both, stochastic processes and habitat filtering are driving factors of AMF community compositions.

Integrating plant species contribution to mycorrhizal and seed dispersal mutualistic networks

Mutualistic interactions like those established between plants and mycorrhizal fungi or seed dispersers are key drivers of plant population dynamics and ecosystem functioning; however, these interactions have rarely been explored together. We assembled a tripartite fungi-plant-disperser network in the Gorongosa National Park-Mozambique, to test (1) if diversity and importance of plant mutualists above- and belowground are correlated, and (2) whether biotically and abiotically dispersed plants are associated with distinct arbuscular mycorrhizal fungi (AMF). We quantified seed dispersal by animals for 1 year and characterized the AMF of 26 common plant species. Sixteen plant species were dispersed by 15 animals and colonized by 48 AMF virtual taxa (VT), while the remaining 10 plant species were not dispersed by animals and associated with 34 AMF VT. We found no evidence for a correlation between the number of plant partners above- and belowground or on plant specialization on both types of partners. We also found no evidence for differentiation of AMF communities between biotically and abiotically dispersed plants. Our results suggest that the establishment of plant interactions with seed dispersers and mycorrhizal fungi is largely independent and that both biotically and abiotically dispersed plants seem to associate with similar communities of AMF.

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.

Symmetric assembly and disassembly processes in an ecological network

The processes whereby ecological networks emerge, persist and decay throughout ecosystem development are largely unknown. Here we study networks of plant and arbuscular mycorrhizal fungal (AMF) communities along a 120 000 year soil chronosequence, as they undergo assembly (progression) and then disassembly (retrogression). We found that network assembly and disassembly were symmetrical, self-reinforcing processes that together were capable of generating key attributes of network architecture. Plant and AMF species that had short indirect paths to others in the community (i.e. high centrality), rather than many direct interaction partners (i.e. high degree), were best able to attract new interaction partners and, in the case of AMF species, also to retain existing interactions with plants during retrogression. We then show using simulations that these non-random patterns of attachment and detachment promote nestedness of the network. These results have implications for predicting extinction sequences, identifying focal points for invasions and suggesting trajectories for restoration.

Host preference and network properties in biotrophic plant-fungal associations

Analytical methods can offer insights into the structure of biological networks, but mechanisms that determine the structure of these networks remain unclear. We conducted a synthesis based on 111 previously published datasets to assess a range of ecological and evolutionary mechanisms that may influence the plant-associated fungal interaction networks. We calculated the relative host effect on fungal community composition and compared nestedness and modularity among different mycorrhizal types and endophytic fungal guilds. We also assessed how plant-fungal network structure was related to host phylogeny, environmental and sampling properties. Orchid mycorrhizal fungal communities responded most strongly to host identity, but the effect of host was similar among all other fungal guilds. Community nestedness, which did not differ among fungal guilds, declined significantly with increasing mean annual precipitation on a global scale. Orchid and ericoid mycorrhizal fungal communities were more modular than ectomycorrhizal and root endophytic communities, with arbuscular mycorrhizal fungi in an intermediate position. Network properties among a broad suite of plant-associated fungi were largely comparable and generally unrelated to phylogenetic distance among hosts. Instead, network metrics were predominantly affected by sampling and matrix properties, indicating the importance of study design in properly inferring ecological patterns.

Temperature-mediated local adaptation alters the symbiotic function in arbuscular mycorrhiza

Variation in the symbiotic function of arbuscular mycorrhizal fungi (AM fungi) has been demonstrated among distinct biotic and abiotic interactions. However, there is little knowledge on how local temperature conditions influence the functional divergence of AM symbionts in alpine ecosystems. Here, we conduct a reciprocal inoculation experiment to explore the three-way interactions among plants, AM fungal inoculum and temperature at sites of contrasting elevation. Evidence of local adaptation of plant growth was found only under low temperature conditions, with no consistent local versus foreign effect found in AM fungal performance. The origin of either the plant or the inoculum relative to the temperature was important in explaining symbiotic function. Specifically, when inoculum and temperature were sympatric but allopatric to the plant, poor adaptation by the plant to the novel environment was clearly found under both temperature conditions. Further analysis found that the symbiotic function was inversely related to fungal diversity under high temperature conditions. These results suggest that local adaptation represents a powerful factor in the establishment of novel combinations of plant, inoculum and temperature, and confirms the importance of taking into account both biotic and abiotic interactions in the prediction of the response of symbionts to global environmental change.

Multiple Assembly Rules Drive the Co-occurrence of Orthopteran and Plant Species in Grasslands: Combining Network, Functional and Phylogenetic Approaches

Understanding the factors underlying the co-occurrence of multiple species remains a challenge in ecology. Biotic interactions, environmental filtering and neutral processes are among the main mechanisms evoked to explain species co-occurrence. However, they are most often studied separately or even considered as mutually exclusive. This likely hampers a more global understanding of species assembly. Here, we investigate the general hypothesis that the structure of co-occurrence networks results from multiple assembly rules and its potential implications for grassland ecosystems. We surveyed orthopteran and plant communities in 48 permanent grasslands of the French Jura Mountains and gathered functional and phylogenetic data for all species. We constructed a network of plant and orthopteran species co-occurrences and verified whether its structure was modular or nested. We investigated the role of all species in the structure of the network (modularity and nestedness). We also investigated the assembly rules driving the structure of the plant-orthopteran co-occurrence network by using null models on species functional traits, phylogenetic relatedness and environmental conditions. We finally compared our results to abundance-based approaches. We found that the plant-orthopteran co-occurrence network had a modular organization. Community assembly rules differed among modules for plants while interactions with plants best explained the distribution of orthopterans into modules. Few species had a disproportionately high positive contribution to this modular organization and are likely to have a key importance to modulate future changes. The impact of agricultural practices was restricted to some modules (3 out of 5) suggesting that shifts in agricultural practices might not impact the entire plant-orthopteran co-occurrence network. These findings support our hypothesis that multiple assembly rules drive the modular structure of the plant-orthopteran network. This modular structure is likely to play a key role in the response of grassland ecosystems to future changes by limiting the impact of changes in agricultural practices such as intensification to some modules leaving species from other modules poorly impacted. The next step is to understand the importance of this modular structure for the long-term maintenance of grassland ecosystem structure and functions as well as to develop tools to integrate network structure into models to improve their capacity to predict future changes.