From diversity to complexity: Microbial networks in soils.

Network analysis has been used for many years in ecological research to analyze organismal associations, for example in food webs, plant-plant or plant-animal interactions. Although network analysis is widely applied in microbial ecology, only recently has it entered the realms of soil microbial ecology, shown by a rapid rise in studies applying co-occurrence analysis to soil microbial communities. While this application offers great potential for deeper insights into the ecological structure of soil microbial ecosystems, it also brings new challenges related to the specific characteristics of soil datasets and the type of ecological questions that can be addressed. In this Perspectives Paper we assess the challenges of applying network analysis to soil microbial ecology due to the small-scale heterogeneity of the soil environment and the nature of soil microbial datasets. We review the different approaches of network construction that are commonly applied to soil microbial datasets and discuss their features and limitations. Using a test dataset of microbial communities from two depths of a forest soil, we demonstrate how different experimental designs and network constructing algorithms affect the structure of the resulting networks, and how this in turn may influence ecological conclusions. We will also reveal how assumptions of the construction method, methods of preparing the dataset, and definitions of thresholds affect the network structure. Finally, we discuss the particular questions in soil microbial ecology that can be approached by analyzing and interpreting specific network properties. Targeting these network properties in a meaningful way will allow applying this technique not in merely descriptive, but in hypothesis-driven research. Analysing microbial networks in soils opens a window to a better understanding of the complexity of microbial communities. However, this approach is unfortunately often used to draw conclusions which are far beyond the scientific evidence it can provide, which has damaged its reputation for soil microbial analysis. In this Perspectives Paper, we would like to sharpen the view for the real potential of microbial co-occurrence analysis in soils, and at the same time raise awareness regarding its limitations and the many ways how it can be misused or misinterpreted.

Plant facilitation through mycorrhizal symbiosis is stronger between distantly related plant species.

The tendency of closely related plant species to share natural enemies has been suggested to limit their co-occurrence and performance, but we lack a deep understanding on how mutualistic interactions such as the mycorrhizal symbiosis affect plant-plant interactions depending on the phylogenetic relatedness of the interacting plants. We hypothesise that the effect of the mycorrhizal symbiosis on plant-plant facilitative interactions depends on the phylogenetic distance between the nurse and facilitated plants. A recently published meta-analysis compiled the strength of plant facilitative interactions in the presence or absence (or reduced abundance) of mycorrhizal fungi. We use phylogenetically informed Bayesian linear models to test whether the effect size is influenced by the phylogenetic distance between the plant species involved in each plant facilitative interaction. Conspecific facilitative interactions are more strongly enhanced by mycorrhizal fungi than interactions between closely related species. In heterospecific interactions, the effect of the mycorrhizal symbiosis on plant facilitation increases with the phylogenetic distance between the nurse and facilitated plant species. Our result showing that the effect of mycorrhizal symbiosis on the facilitation interactions between plants depends on their phylogenetic relatedness provides new mechanisms to understand how facilitation is assembling ecological communities.

Unexpected bird-feather mite associations revealed by DNA metabarcoding uncovers a dynamic ecoevolutionary scenario.

The high relevance of host-switching for the diversification of highly host-specific symbionts (i.e., those commonly inhabiting a single host species) demands a better understanding of host-switching dynamics at an ecological scale. Here, we used DNA metabarcoding to study feather mites on passerine birds in Spain, sequencing mtDNA (COI) for 25,540 individual mites (representing 64 species) from 1,130 birds (representing 71 species). Surprisingly, 1,228 (4.8%) mites from 84 (7.4%) birds were found on host species that were not the expected to be a host according to a recent bird-feather mite associations catalog. Unexpected associations were widespread across studied mite (40.6%) and bird (43.7%) species and showed smaller average infrapopulation sizes than typical associations. Unexpected mite species colonized hosts being distantly related to the set of their usual hosts, but with similar body size. The network of bird-mite associations was modular (i.e., some groups of bird and mite species tended to be more associated with each other than with the others), with 75.9% of the unexpected associations appearing within the module of the typical hosts of the mite species. Lastly, 68.4% of mite species found on unexpected hosts showed signatures of genetic differentiation, and we found evidence for reproduction or the potential for it in many of the unexpected associations. Results show host colonization as a common phenomenon even for these putatively highly host-specific symbionts. Thus, host-switching by feather mites, rather than a rare phenomenon, appears as a relatively frequent phenomenon shaped by ecological filters such as host morphology and is revealed as a fundamental component for a dynamic coevolutionary and codiversification scenario.

Nurse plants transfer more nitrogen to distantly related species.

Plant facilitative interactions enhance co-occurrence between distant relatives, partly due to limited overlap in resource requirements. We propose a different mechanism for the coexistence of distant relatives based on positive interactions of nutrient sharing. Nutrients move between plants following source-sink gradients driven by plant traits that allow these gradients to establish. Specifically, nitrogen (N) concentration gradients can arise from variation in leaf N content across plants species. As many ecologically relevant traits, we hypothesize that leaf N content is phylogenetically conserved and can result in N gradients promoting N transfer among distant relatives. In a Mexican desert community governed by facilitation, we labelled nurse plants (Mimosa luisana) with 15 N and measured its transfer to 14 other species in the community, spanning the range of phylogenetic distances to the nurse plant. Nurses established steeper N source-sink gradients with distant relatives, increasing 15 N transfer toward these species. Nutrient sharing may provide long-term benefits to facilitated plants and may be an overlooked mechanism maintaining coexistence and increasing the phylogenetic diversity of plant communities.

Soil Characteristics Driving Arbuscular Mycorrhizal Fungal Communities in Semiarid Mediterranean Soils.

UNLABELLED: We investigated communities of arbuscular mycorrhizal fungi (AMF) in the roots and the rhizosphere soil of Brachypodium retusum in six different natural soils under field conditions. We explored phylogenetic patterns of AMF composition using indicator species analyses to find AMF associated with a given habitat (root versus rhizosphere) or soil type. We tested whether the AMF characteristics of different habitats or contrasting soils were more closely related than expected by chance. Then we used principal-component analysis and multivariate analysis of variance to test for the relative contribution of each factor in explaining the variation in fungal community composition. Finally, we used redundancy analysis to identify the soil properties that significantly explained the differences in AMF communities across soil types. The results pointed out a tendency of AMF communities in roots to be closely related and different from those in the rhizosphere soil. The indicator species analyses revealed AMF associated with rhizosphere soil and the root habitat. Soil type also determined the distribution of AMF communities in soils, and this effect could not be attributed to a single soil characteristic, as at least three soil properties related to microbial activity, i.e., pH and levels of two micronutrients (Mn and Zn), played significant roles in triggering AMF populations. IMPORTANCE: Communities of arbuscular mycorrhizal fungi (AMF) are main components of soil biota that can determine the productivity of ecosystems. These fungal assemblages vary across host plants and ecosystems, but the main ecological processes that shape the structures of these communities are still largely unknown. A field study in six different soil types from semiarid areas revealed that AMF communities are significantly influenced by habitat (soil versus roots) and soil type. In addition, three soil properties related to microbiological activity (i.e., pH and manganese and zinc levels) were the main factors triggering the distribution of AMF. These results contribute to a better understanding of the ecological factors that can shape AMF communities, an important soil microbial group that affects multiple ecosystem functions.

Fungal phylogenetic diversity drives plant facilitation.

Plant-plant facilitation is a crucial ecological process, as many plant species (facilitated) require the presence of an established individual (nurse) to recruit. Some plant facilitative interactions disappear during the ontogenetic development of the facilitated plant but others persist, even when the two plants are adults. We test whether the persistence of plant facilitative interactions is explained by the phylogenetic diversity of mutualistic and non-mutualistic fungi that the nurse and the facilitated species add to the shared rhizosphere. We classify plant facilitative interactions as persistent and non-persistent interactions and quantify the phylogenetic diversity of mutualistic and non-mutualistic fungi added by the plant species to the shared rhizosphere. Our results show that the facilitated species add less phylogenetic diversity of non-mutualistic fungi when plant facilitative interactions persist than when they do not persist. However, persistent and non-persistent facilitative interactions did not differ in the phylogenetic diversity of mutualistic fungi added by the facilitated species to the shared rhizosphere. Finally, the fungal phylogenetic diversity added by the nurse to the shared rhizosphere did not differ between persistent and non-persistent interactions. This study suggests that considering the fungal associates of the plant species involved in facilitative interactions can shed light on the mechanisms of persistence for plant-plant interactions.

Modularity reveals the tendency of arbuscular mycorrhizal fungi to interact differently with generalist and specialist plant species in gypsum soils.

Patterns in plant-soil biota interactions could be influenced by the spatial distribution of species due to soil conditions or by the functional traits of species. Gypsum environments usually constitute a mosaic of heterogeneous soils where gypsum and nongypsum soils are imbricated at a local scale. A case study of the interactions of plants with arbuscular mycorrhizal fungi (AMF) in gypsum environments can be illustrative of patterns in biotic interactions. We hypothesized that (i) soil characteristics might affect the AMF community and (ii) there are differences between the AMF communities (modules) associated with plants exclusive to gypsum soils (gypsophytes) and those associated with plants that show facultative behavior on gypsum and/or marly-limestone soils (gypsovags). We used indicator species and network analyses to test for differences between the AMF communities harbored in gypsophyte and gypsovag plants. We recorded 46 operational taxonomic units (OTUs) belonging to nine genera of Glomeromycota. The indicator species analysis showed two OTUs preferentially associating with gypsum soils and three OTUs preferentially associating with marly-limestone soils. Modularity analysis revealed that soil type can be a major factor shaping AMF communities, and some AMF groups showed a tendency to interact differently with plants that had distinct ecological strategies (gypsophytes and gypsovags). Characterization of ecological networks can be a valuable tool for ascertaining the potential influence of above- and below-ground biotic interactions (plant-AMF) on plant community composition.

Plant phylodiversity enhances soil microbial productivity in facilitation-driven communities.

The classical relationship between biodiversity and ecosystem functioning can be better understood when the phylogenetic component of biodiversity is considered. We linked plant phylodiversity and ecosystem functioning in a water-limited gypsum ecosystem driven by plant facilitation. We tested whether (1) plant facilitation relaxes the abiotic filter imposed by gypsum, allowing the establishment of non-gypsophyte plant species, and consequently increasing plant phylodiversity, and (2) plant phylodiversity influences soil microbial productivity. Our data revealed that the gypsophyte Ononis tridentata spatially determines a macrophytic mosaic, ameliorates the microenvironment, and maximizes plant richness and phylodiversity through facilitating non-gypsophyte species. Beyond the direct effect of the nurse plant on soil microbial biomass, activity, and respiration, the analyses suggest a direct effect of plant phylodiversity (MPD) on these general indicators of soil microbial productivity. Plant diversity (Shannon index) neither correlated with the mentioned parameters nor with specific indicators of C, N and P cycling. This is the first report of a relationship between producer phylodiversity and decomposer productivity, which supports phylogenetic diversity as a relevant player of the ecosystem functioning.

Plant metabolic response to stress in an arid ecosystem is mediated by the presence of neighbors.

Plant neighbors in arid environments can ameliorate abiotic stress by reducing insolation, but they also attract herbivores and pathogens, especially when neighbors are close relatives that share similar antagonists. Plants' metabolic profiles provide a chemical fingerprint of the physiological processes behind plant responses to different environmental stresses. For example, abscisic acid and proline, mainly involved in stomatal closure and osmotic adjustment, can induce plant responses to abiotic stress, while jasmonic acid and salicylic acid primarily regulate plant defense to herbivory or pathogens. Neighbor plants can generate contrasting ecological contexts, modulating plant responses to abiotic and biotic stresses. We hypothesize that plant metabolic profile is modulated by its neighbors in a vegetation patch, expecting a higher investment in metabolites related to biotic-stress tolerance (i.e., herbivory or pathogens) when growing associated with other plants, especially to phylogenetically close relatives, compared to plants growing alone. We show that plants from five species growing with neighbors invest more in biotic-stress tolerance while their conspecifics, growing alone, invest more in abiotic-stress tolerance. This tendency in plants' metabolic profiles was not affected by the phylogenetic diversity of their neighborhood. Linking physiological snapshots with community processes can contribute to elucidating metabolic profiles derived from plant-plant interactions.

Phylogenetic and functional constraints of plant facilitation rewiring.

Facilitative interactions bind community species in intricate ecological networks, preserving species that would otherwise be lost. The traditional understanding of ecological networks as static components of biological communities overlooks the fact that species interactions in a network can fluctuate. Analyzing the patterns that cause those shifts can reveal the principles that govern the identity of pairwise interactions and whether they are predictable based on the traits of the interacting species and the local environmental contexts in which they occur. Here we explore how abiotic stress and phylogenetic and functional affinities constrain those shifts. Specifically, we hypothesize that rewiring the facilitative interactions is more limited in stressful than in mild environments. We present evidence of a distinct pattern in the rewiring of facilitation-driven communities at different stress levels. In highly stressful environments with a firm reliance on facilitation, rewiring is limited to growing beneath nurse species with traits to overcome harsh stressful conditions. However, when environments are milder, rewiring is more flexible, although it is still constrained to nurses that are close relatives. Understanding the ability of species to rewire their interactions is crucial for predicting how communities may respond to the unprecedented rate of perturbations on Earth.

RecruitNet: A global database of plant recruitment networks.

Plant recruitment interactions (i.e., what recruits under what) shape the composition, diversity, and structure of plant communities. Despite the huge body of knowledge on the mechanisms underlying recruitment interactions among species, we still know little about the structure of the recruitment networks emerging in ecological communities. Modeling and analyzing the community-level structure of plant recruitment interactions as a complex network can provide relevant information on ecological and evolutionary processes acting both at the species and ecosystem levels. We report a data set containing 143 plant recruitment networks in 23 countries across five continents, including temperate and tropical ecosystems. Each network identifies the species under which another species recruits. All networks report the number of recruits (i.e., individuals) per species. The data set includes >850,000 recruiting individuals involved in 118,411 paired interactions among 3318 vascular plant species across the globe. The cover of canopy species and open ground is also provided. Three sampling protocols were used: (1) The Recruitment Network (RN) protocol (106 networks) focuses on interactions among established plants ("canopy species") and plants in their early stages of recruitment ("recruit species"). A series of plots was delimited within a locality, and all the individuals recruiting and their canopy species were identified; (2) The paired Canopy-Open (pCO) protocol (26 networks) consists in locating a potential canopy plant and identifying recruiting individuals under the canopy and in a nearby open space of the same area; (3) The Georeferenced plot (GP) protocol (11 networks) consists in using information from georeferenced individual plants in large plots to infer canopy-recruit interactions. Some networks incorporate data for both herbs and woody species, whereas others focus exclusively on woody species. The location of each study site, geographical coordinates, country, locality, responsible author, sampling dates, sampling method, and life habits of both canopy and recruit species are provided. This database will allow researchers to test ecological, biogeographical, and evolutionary hypotheses related to plant recruitment interactions. There are no copyright restrictions on the data set; please cite this data paper when using these data in publications.

Both abundant and rare fungi colonizing Fagus sylvatica ectomycorrhizal root-tips shape associated bacterial communities.

Ectomycorrhizal fungi live in close association with their host plants and form complex interactions with bacterial/archaeal communities in soil. We investigated whether abundant or rare ectomycorrhizal fungi on root-tips of young beech trees (Fagus sylvatica) shape bacterial/archaeal communities. We sequenced 16S rRNA genes and fungal internal transcribed spacer regions of individual root-tips and used ecological networks to detect the tendency of certain assemblies of fungal and bacterial/archaeal taxa to inhabit the same root-tip (i.e. modularity). Individual ectomycorrhizal root-tips hosted distinct fungal communities associated with unique bacterial/archaeal communities. The structure of the fungal-bacterial/archaeal association was determined by both, dominant and rare fungi. Integrating our data in a conceptual framework suggests that the effect of rare fungi on the bacterial/archaeal communities of ectomycorrhizal root-tips contributes to assemblages of bacteria/archaea on root-tips. This highlights the potential impact of complex fine-scale interactions between root-tip associated fungi and other soil microorganisms for the ectomycorrhizal symbiosis.

Arabidopsis thaliana populations show clinal variation in a climatic gradient associated with altitude.

• Understanding the adaptive basis of life history variation is a central goal in evolutionary ecology. The use of model species enables the combination of molecular mechanistic knowledge with ecological and evolutionary questions, but the study of life history variation in natural environments is required to merge these disciplines. • Here, we tested for clinal variation in life history and associated traits along an environmental and altitudinal gradient in the model species Arabidopsis thaliana. Seventeen natural populations of A. thaliana were geo-referenced in north-eastern Spain on a gradient in which precipitation increases but maximum spring temperature and minimum winter temperature decrease with altitude. • One hundred and eighty-nine genotypes from the 17 populations were grown under uniform controlled conditions. Variations in traits related to biomass allocation, fecundity, phenology and vegetative growth were tested for relationships with the altitude and climatic variables associated with the home sites. Above-ground mass, number of rosette leaves at bolting, developmental time and seed weight increased with the home site's altitude. Root allocation, vegetative growth during winter and number of seeds decreased with altitude. • We suggest that the differences among home sites provide clues to the variation in adaptive strategies associated with the climatic gradient. We compared these results with adaptations and clinal relationships reported for other species and with molecular mechanisms described in Arabidopsis.

Temporal variation in genetic diversity and effective population size of Mediterranean and subalpine Arabidopsis thaliana populations.

Currently, there exists a limited knowledge on the extent of temporal variation in population genetic parameters of natural populations. Here, we study the extent of temporal variation in population genetics by genotyping 151 genome-wide SNP markers polymorphic in 466 individuals collected from nine populations of the annual plant Arabidopsis thaliana during 4 years. Populations are located along an altitudinal climatic gradient from Mediterranean to subalpine environments in NE Spain, which has been shown to influence key demographic attributes and life cycle adaptations. Genetically, A. thaliana populations were more variable across space than over time. Common multilocus genotypes were detected several years in the same population, whereas low-frequency multilocus genotypes appeared only 1 year. High-elevation populations were genetically poorer and more variable over time than low-elevation populations, which might be caused by a higher overall demographic instability at higher altitudes. Estimated effective population sizes were low but also showed a significant decreasing trend with increasing altitude, suggesting a deeper impact of genetic drift at high-elevation populations. In comparison with single-year samplings, repeated genotyping over time captured substantially higher amount of genetic variation contained in A. thaliana populations. Furthermore, repeated genotyping of populations provided novel information on the genetic properties of A. thaliana populations and allowed hypothesizing on their underlying mechanisms. Therefore, including temporal genotyping programmes into traditional population genetic studies can significantly increase our understanding of the dynamics of natural populations.

Functional trait effects on ecosystem stability: assembling the jigsaw puzzle.

Under global change, how biological diversity and ecosystem services are maintained in time is a fundamental question. Ecologists have long argued about multiple mechanisms by which local biodiversity might control the temporal stability of ecosystem properties. Accumulating theories and empirical evidence suggest that, together with different population and community parameters, these mechanisms largely operate through differences in functional traits among organisms. We review potential trait-stability mechanisms together with underlying tests and associated metrics. We identify various trait-based components, each accounting for different stability mechanisms, that contribute to buffering, or propagating, the effect of environmental fluctuations on ecosystem functioning. This comprehensive picture, obtained by combining different puzzle pieces of trait-stability effects, will guide future empirical and modeling investigations.

Phylogenetic rewiring in mycorrhizal-plant interaction networks increases community stability in naturally fragmented landscapes.

Although ecological networks are usually considered a static representation of species' interactions, the interactions can change when the preferred partners are absent (rewiring). In mutualistic networks, rewiring with non-preferred partners can palliate extinction cascades, contributing to communities' stability. In spite of its significance, whether general patterns can shape the rewiring of ecological interactions remains poorly understood. Here, we show a phylogenetic constraint in the rewiring of mycorrhizal networks, so that rewired interactions (i.e., with non-preferred hosts) tend to involve close relatives of preferred hosts. Despite this constraint, rewiring increases the robustness of the fungal community to the simulated loss of their host species. We identify preferred and non-preferred hosts based on the probability that, when the two partners co-occur, they actually interact. Understanding general patterns in the rewiring of interactions can improve our predictions of community responses to interactions' loss, which influences how global changes will affect ecosystem stability.

Incorporating phylogenetic metrics to microbial co-occurrence networks based on amplicon sequences to discern community assembly processes.

Co-occurrence network analysis based on amplicon sequences is increasingly used to study microbial communities. Patterns of co-existence or mutual exclusion between pairs of taxa are often interpreted as reflecting positive or negative biological interactions. However, other assembly processes can underlie these patterns, including species failure to reach distant areas (dispersal limitation) and tolerate local environmental conditions (habitat filtering). We provide a tool to quantify the relative contribution of community assembly processes to microbial co-occurrence patterns, which we applied to explore soil bacterial communities in two dry ecosystems. First, we sequenced a bacterial phylogenetic marker in soils collected across multiple plots. Second, we inferred co-occurrence networks to identify pairs of significantly associated taxa, either co-existing more (aggregated) or less often (segregated) than expected at random. Third, we assigned assembly processes to each pair: patterns explained based on spatial or environmental distance were ascribed to dispersal limitation (2%-4%) or habitat filtering (55%-77%), and the remaining to biological interactions. Finally, we calculated the phylogenetic distance between taxon pairs to test theoretical expectations on the linkages between phylogenetic patterns and assembly processes. Aggregated pairs were more closely related than segregated pairs. Furthermore, habitat-filtered aggregated pairs were closer relatives than those assigned to positive interactions, consistent with phylogenetic niche conservatism and cooperativism among distantly related taxa. Negative interactions resulted in equivocal phylogenetic signatures, probably because different competitive processes leave opposing signals. We show that microbial co-occurrence networks mainly reflect environmental tolerances and propose that incorporating measures of phylogenetic relatedness to networks might help elucidate ecologically meaningful patterns.

Correction: Community structure informs species geographic distributions.

[This corrects the article DOI: 10.1371/journal.pone.0197877.].

Community structure informs species geographic distributions.

Understanding what determines species' geographic distributions is crucial for assessing global change threats to biodiversity. Measuring limits on distributions is usually, and necessarily, done with data at large geographic extents and coarse spatial resolution. However, survival of individuals is determined by processes that happen at small spatial scales. The relative abundance of coexisting species (i.e. 'community structure') reflects assembly processes occurring at small scales, and are often available for relatively extensive areas, so could be useful for explaining species distributions. We demonstrate that Bayesian Network Inference (BNI) can overcome several challenges to including community structure into studies of species distributions, despite having been little used to date. We hypothesized that the relative abundance of coexisting species can improve predictions of species distributions. In 1570 assemblages of 68 Mediterranean woody plant species we used BNI to incorporate community structure into Species Distribution Models (SDMs), alongside environmental information. Information on species associations improved SDM predictions of community structure and species distributions moderately, though for some habitat specialists the deviance explained increased by up to 15%. We demonstrate that most species associations (95%) were positive and occurred between species with ecologically similar traits. This suggests that SDM improvement could be because species co-occurrences are a proxy for local ecological processes. Our study shows that Bayesian Networks, when interpreted carefully, can be used to include local conditions into measurements of species' large-scale distributions, and this information can improve the predictions of species distributions.

Network structure embracing mutualism-antagonism continuums increases community robustness.

Theory predicts that contrasting properties of mutualistic and antagonistic networks differentially promote community resilience to species loss. However, the outcome of most ecological interactions falls within a continuum between mutualism and antagonism, and we ignore the extent to which this interactions' continuum might influence community stability. Using a large data set of interactions, we compared co-extinction cascades that either consider or ignore the mix of beneficial and detrimental actions that parrots exert on plants. When the antagonism-mutualism continuum was considered, a combination of the properties that separately enhance community stability in ecological networks emerged. This combination of properties led to an overall increase of the parrot community robustness to face plant species loss. Our results highlight that the conditional outcomes of interactions can influence the structure of ecological networks, thus affecting our predictions of community stability against eventual changes.