Soil prokaryotic and fungal biome structures associated with crop disease status across the Japan Archipelago.

Archaea, bacteria, and fungi in the soil are increasingly recognized as determinants of agricultural productivity and sustainability. A crucial step for exploring soil microbiomes with important ecosystem functions is to perform statistical analyses on the potential relationship between microbiome structure and functions based on comparisons of hundreds or thousands of environmental samples collected across broad geographic ranges. In this study, we integrated agricultural field metadata with microbial community analyses by targeting 2,903 bulk soil samples collected along a latitudinal gradient from cool-temperate to subtropical regions in Japan (26.1-42.8 °N). The data involving 632 archaeal, 26,868 bacterial, and 4,889 fungal operational taxonomic units detected across the fields of 19 crop plant species allowed us to conduct statistical analyses (permutational analyses of variance, generalized linear mixed models, randomization analyses, and network analyses) on the relationship among edaphic factors, microbiome compositions, and crop disease prevalence. We then examined whether the diverse microbes form species sets varying in potential ecological impacts on crop plants. A network analysis suggested that the observed prokaryotes and fungi were classified into several species sets (network modules), which differed substantially in association with crop disease prevalence. Within the network of microbe-to-microbe coexistence, ecologically diverse microbes, such as an ammonium-oxidizing archaeon, an antibiotics-producing bacterium, and a potentially mycoparasitic fungus, were inferred to play key roles in shifts between crop-disease-promotive and crop-disease-suppressive states of soil microbiomes. The bird's-eye view of soil microbiome structure will provide a basis for designing and managing agroecosystems with high disease-suppressive functions.IMPORTANCEUnderstanding how microbiome structure and functions are organized in soil ecosystems is one of the major challenges in both basic ecology and applied microbiology. Given the ongoing worldwide degradation of agroecosystems, building frameworks for exploring structural diversity and functional profiles of soil microbiomes is an essential task. Our study provides an overview of cropland microbiome states in light of potential crop-disease-suppressive functions. The large data set allowed us to explore highly functional species sets that may be stably managed in agroecosystems. Furthermore, an analysis of network architecture highlighted species that are potentially used to cause shifts from disease-prevalent states of agroecosystems to disease-suppressive states. By extending the approach of comparative analyses toward broader geographic ranges and diverse agricultural practices, agroecosystem with maximized biological functions will be further explored.

Interaction network rewiring and species' contributions to community-scale flexibility.

The architecture of species interaction networks is a key factor determining the stability of ecological communities. However, the fact that ecological network architecture can change through time is often overlooked in discussions on community-level processes, despite its theoretical importance. By compiling a time-series community dataset involving 50 spider species and 974 Hexapoda prey species/strains, we quantified the extent to which the architecture of predator-prey interaction networks could shift across time points. We then developed a framework for finding species that could increase the flexibility of the interaction network architecture. Those "network coordinator" species are expected to promote the persistence of species-rich ecological communities by buffering perturbations in communities. Although spiders are often considered as generalist predators, their contributions to network flexibility vary greatly among species. We also found that detritivorous prey species can be cores of interaction rewiring, dynamically interlinking below-ground and above-ground community dynamics. We further found that the predator-prey interactions between those network coordinators differed from those highlighted in the standard network-analytical framework assuming static topology. Analyses of network coordinators will add a new dimension to our understanding of species coexistence mechanisms and provide platforms for systematically prioritizing species in terms of their potential contributions in ecosystem conservation and restoration.

Deterministic and stochastic processes generating alternative states of microbiomes.

The structure of microbiomes is often classified into discrete or semi-discrete types potentially differing in community-scale functional profiles. Elucidating the mechanisms that generate such "alternative states" of microbiome compositions has been one of the major challenges in ecology and microbiology. In a time-series analysis of experimental microbiomes, we here show that both deterministic and stochastic ecological processes drive divergence of alternative microbiome states. We introduced species-rich soil-derived microbiomes into eight types of culture media with 48 replicates, monitoring shifts in community compositions at six time points (8 media × 48 replicates × 6 time points = 2304 community samples). We then confirmed that microbial community structure diverged into a few state types in each of the eight medium conditions as predicted in the presence of both deterministic and stochastic community processes. In other words, microbiome structure was differentiated into a small number of reproducible compositions under the same environment. This fact indicates not only the presence of selective forces leading to specific equilibria of community-scale resource use but also the influence of demographic drift (fluctuations) on the microbiome assembly. A reference-genome-based analysis further suggested that the observed alternative states differed in ecosystem-level functions. These findings will help us examine how microbiome structure and functions can be controlled by changing the "stability landscapes" of ecological community compositions.

Leaf, root, and soil microbiomes of an invasive plant, Ardisia crenata, differ between its native and exotic ranges.

Introduction: Ecological underpinnings of the invasion success of exotic plants may be found in their interactions with microbes, either through the enemy release hypothesis and the enhanced mutualism hypothesis. Whereas recent high-throughput sequencing techniques have significantly expanded our understanding of plant-associated microbiomes and their functional guilds, few studies to date have used these techniques to compare the microbiome associated with invasive plants between their native and exotic ranges. Methods: We extracted fungal and bacterial DNA within leaf endosphere, root endosphere and soil of an invasive plant, Ardisia crenata, sampled from their native range Japan and exotic range Florida, USA. Using Illumina sequencing data, we compared microbial community compositions and diversity between the native and exotic ranges, and tested whether abundance of pathogenic or mutualistic microbes differ between the native or exotic ranges in accordance to the enemy release hypothesis or the enhanced mutualism hypothesis. Results: Fungal and bacterial community compositions differed among leaves, roots and soil, and between the native and exotic ranges. Despite a higher microbial diversity in the soil in the exotic range than in the native range, the microbial diversity within leaf and root was lower in the exotic range compared to the native range. In addition, leaves in the native range harbored a greater number of plant pathogenic fungi compared to those in the exotic range. Discussion: These patterns suggest plant controls over what microbes become associated with leaves and roots. The higher abundance of leaf pathogenic fungi, including the pathogen which is known to cause specific disease in A. crenata in the exotic range than in the native range, support the enemy release hypothesis and highlighted potential importance of examining microbial communities both above- and below-ground.

Metagenomic analysis of ecological niche overlap and community collapse in microbiome dynamics.

Species utilizing the same resources often fail to coexist for extended periods of time. Such competitive exclusion mechanisms potentially underly microbiome dynamics, causing breakdowns of communities composed of species with similar genetic backgrounds of resource utilization. Although genes responsible for competitive exclusion among a small number of species have been investigated in pioneering studies, it remains a major challenge to integrate genomics and ecology for understanding stable coexistence in species-rich communities. Here, we examine whether community-scale analyses of functional gene redundancy can provide a useful platform for interpreting and predicting collapse of bacterial communities. Through 110-day time-series of experimental microbiome dynamics, we analyzed the metagenome-assembled genomes of co-occurring bacterial species. We then inferred ecological niche space based on the multivariate analysis of the genome compositions. The analysis allowed us to evaluate potential shifts in the level of niche overlap between species through time. We hypothesized that community-scale pressure of competitive exclusion could be evaluated by quantifying overlap of genetically determined resource-use profiles (metabolic pathway profiles) among coexisting species. We found that the degree of community compositional changes observed in the experimental microbiome was correlated with the magnitude of gene-repertoire overlaps among bacterial species, although the causation between the two variables deserves future extensive research. The metagenome-based analysis of genetic potential for competitive exclusion will help us forecast major events in microbiome dynamics such as sudden community collapse (i.e., dysbiosis).

Dynamics of species-rich predator-prey networks and seasonal alternations of core species.

In nature, entangled webs of predator-prey interactions constitute the backbones of ecosystems. Uncovering the network architecture of such trophic interactions has been recognized as the essential step for exploring species with great impacts on ecosystem-level phenomena and functions. However, it has remained a major challenge to reveal how species-rich networks of predator-prey interactions are continually reshaped through time in the wild. Here, we show that dynamics of species-rich predator-prey interactions can be characterized by remarkable network structural changes and alternations of species with greatest impacts on community processes. On the basis of high-throughput detection of prey DNA from 1,556 spider individuals collected in a grassland ecosystem, we reconstructed dynamics of interaction networks involving, in total, 50 spider species and 974 prey species and strains through 8 months. The networks were compartmentalized into modules (groups) of closely interacting predators and prey in each month. Those modules differed in detritus/grazing food chain properties, forming complex fission-fusion dynamics of belowground and aboveground energy channels across the seasons. The substantial shifts of network structure entailed alternations of spider species located at the core positions within the entangled webs of interactions. These results indicate that knowledge of dynamically shifting food webs is crucial for understanding temporally varying roles of 'core species' in ecosystem processes.

Facilitative interaction networks in experimental microbial community dynamics.

Facilitative interactions between microbial species are ubiquitous in various types of ecosystems on the Earth. Therefore, inferring how entangled webs of interspecific interactions shift through time in microbial ecosystems is an essential step for understanding ecological processes driving microbiome dynamics. By compiling shotgun metagenomic sequencing data of an experimental microbial community, we examined how the architectural features of facilitative interaction networks could change through time. A metabolic modeling approach for estimating dependence between microbial genomes (species) allowed us to infer the network structure of potential facilitative interactions at 13 time points through the 110-day monitoring of experimental microbiomes. We then found that positive feedback loops, which were theoretically predicted to promote cascade breakdown of ecological communities, existed within the inferred networks of metabolic interactions prior to the drastic community-compositional shift observed in the microbiome time-series. We further applied "directed-graph" analyses to pinpoint potential keystone species located at the "upper stream" positions of such feedback loops. These analyses on facilitative interactions will help us understand key mechanisms causing catastrophic shifts in microbial community structure.

Alternative stable states, nonlinear behavior, and predictability of microbiome dynamics.

BACKGROUND: Microbiome dynamics are both crucial indicators and potential drivers of human health, agricultural output, and industrial bio-applications. However, predicting microbiome dynamics is notoriously difficult because communities often show abrupt structural changes, such as "dysbiosis" in human microbiomes. METHODS: We integrated theoretical frameworks and empirical analyses with the aim of anticipating drastic shifts of microbial communities. We monitored 48 experimental microbiomes for 110 days and observed that various community-level events, including collapse and gradual compositional changes, occurred according to a defined set of environmental conditions. We analyzed the time-series data based on statistical physics and non-linear mechanics to describe the characteristics of the microbiome dynamics and to examine the predictability of major shifts in microbial community structure. RESULTS: We confirmed that the abrupt community changes observed through the time-series could be described as shifts between "alternative stable states" or dynamics around complex attractors. Furthermore, collapses of microbiome structure were successfully anticipated by means of the diagnostic threshold defined with the "energy landscape" analysis of statistical physics or that of a stability index of nonlinear mechanics. CONCLUSIONS: The results indicate that abrupt microbiome events in complex microbial communities can be forecasted by extending classic ecological concepts to the scale of species-rich microbial systems. Video Abstract.

Core species and interactions prominent in fish-associated microbiome dynamics.

BACKGROUND: In aquatic ecosystems, the health and performance of fish depend greatly on the dynamics of microbial community structure in the background environment. Nonetheless, finding microbes with profound impacts on fish's performance out of thousands of candidate species remains a major challenge. METHODS: We examined whether time-series analyses of microbial population dynamics could illuminate core components and structure of fish-associated microbiomes in the background (environmental) water. By targeting eel-aquaculture-tank microbiomes as model systems, we reconstructed the population dynamics of the 9605 bacterial and 303 archaeal species/strains across 128 days. RESULTS: Due to the remarkable increase/decrease of constituent microbial population densities, the taxonomic compositions of the microbiome changed drastically through time. We then found that some specific microbial taxa showed a positive relationship with eels' activity levels even after excluding confounding effects of environmental parameters (pH and dissolved oxygen level) on population dynamics. In particular, a vitamin-B12-producing bacteria, Cetobacterium somerae, consistently showed strong positive associations with eels' activity levels across the replicate time series of the five aquaculture tanks analyzed. Network theoretical and metabolic modeling analyses further suggested that the highlighted bacterium and some other closely-associated bacteria formed "core microbiomes" with potentially positive impacts on eels. CONCLUSIONS: Overall, these results suggest that the integration of microbiology, ecological theory, and network science allows us to explore core species and interactions embedded within complex dynamics of fish-associated microbiomes.  Video Abstract.

Synergistic and Offset Effects of Fungal Species Combinations on Plant Performance.

In natural and agricultural ecosystems, survival and growth of plants depend substantially on residing microbes in the endosphere and rhizosphere. Although numerous studies have reported the presence of plant-growth promoting bacteria and fungi in below-ground biomes, it remains a major challenge to understand how sets of microbial species positively or negatively affect plants' performance. By conducting a series of single- and dual-inoculation experiments of 13 plant-associated fungi targeting a Brassicaceae plant species (Brassica rapa var. perviridis), we here systematically evaluated how microbial effects on plants depend on presence/absence of co-occurring microbes. The comparison of single- and dual-inoculation experiments showed that combinations of the fungal isolates with the highest plant-growth promoting effects in single inoculations did not have highly positive impacts on plant performance traits (e.g., shoot dry weight). In contrast, pairs of fungi with small/moderate contributions to plant growth in single-inoculation contexts showed the greatest effects on plants among the 78 fungal pairs examined. These results on the offset and synergistic effects of pairs of microbes suggest that inoculation experiments of single microbial species/isolates can result in the overestimation or underestimation of microbial functions in multi-species contexts. Because keeping single-microbe systems under outdoor conditions is impractical, designing sets of microbes that can maximize performance of crop plants is an important step for the use of microbial functions in sustainable agriculture.

Aboveground herbivores drive stronger plant species-specific feedback than belowground fungi to regulate tree community assembly.

Ectomycorrhizal (EcM) tree species often become more dominant than arbuscular mycorrhizal (AM) tree species in temperate forests, but they generally coexist. Theory predicts that ecological feedback mediated by aboveground herbivory and/or belowground microbes could explain these dominance/coexistence patterns. An experimental test of how aboveground/belowground organisms associated with AM/EcM trees mediate ecological feedbacks has been lacking at the community-level. By establishing AM and EcM tree sapling assemblages in mesocosms and then introducing seedlings of each type in a reciprocal planting experiment, we compared seedling performance under varying sapling species (conspecifics, heterospecifics within the same and different mycorrhizal types), using traits that reflect either aboveground herbivory-mediated feedback or belowground fungal-mediated feedback or both. When examining seedling traits that reflect aboveground herbivory-mediated feedbacks (i.e., foliar damage), AM plants tended to experience less foliar damage and EcM plants more damage under conspecific versus heterospecific saplings within the same mycorrhizal types, and aboveground herbivory-mediated feedback was species-specific rather than mycorrhizal type-specific. Conversely, when examining traits that reflect belowground fungal-mediated feedbacks, both AM and EcM plant species often exhibited mycorrhizal type-specific feedbacks (e.g., greater aboveground biomass under the same versus different mycorrhizal-type saplings) rather than species-specific feedbacks. Furthermore, tree species affected by herbivory-mediated feedback were less affected by belowground feedback, indicating that the relative importance of the feedbacks varied among plant species. Analysis of plant-associated organisms verified that the feedback outcomes corresponded with species accumulation of belowground fungi (but not of aboveground herbivores). Thus, aboveground herbivores drive stronger plant species-specific feedback than belowground fungi to regulate temperate tree diversity.

Scoring Species for Synthetic Community Design: Network Analyses of Functional Core Microbiomes.

Constructing biological communities is a major challenge in both basic and applied sciences. Although model synthetic communities with a few species have been constructed, designing systems consisting of tens or hundreds of species remains one of the most difficult goals in ecology and microbiology. By utilizing high-throughput sequencing data of interspecific association networks, we here propose a framework for exploring "functional core" species that have great impacts on whole community processes and functions. The framework allows us to score each species within a large community based on three criteria: namely, topological positions, functional portfolios, and functional balance within a target network. The criteria are measures of each species' roles in maximizing functional benefits at the community or ecosystem level. When species with potentially large contributions to ecosystem-level functions are screened, the framework also helps us design "functional core microbiomes" by focusing on properties of species groups (modules) within a network. When embedded into agroecosystems or human gut, such functional core microbiomes are expected to organize whole microbiome processes and functions. An application to a plant-associated microbiome dataset actually highlighted potential functional core microbes that were known to control rhizosphere microbiomes by suppressing pathogens. Meanwhile, an example of application in mouse gut microbiomes called attention to poorly investigated bacterial species, whose potential roles within gut microbiomes deserve future experimental studies. The framework for gaining "bird's-eye" views of functional cores within networks is applicable not only to agricultural and medical data but also to datasets produced in food processing, brewing, waste water purification, and biofuel production.

Consortia of anti-nematode fungi and bacteria in the rhizosphere of soybean plants attacked by root-knot nematodes.

Cyst and root-knot nematodes are major risk factors of agroecosystem management, often causing devastating impacts on crop production. The use of microbes that parasitize or prey on nematodes has been considered as a promising approach for suppressing phytopathogenic nematode populations. However, effects and persistence of those biological control agents often vary substantially depending on regions, soil characteristics and agricultural practices: more insights into microbial community processes are required to develop reproducible control of nematode populations. By performing high-throughput sequencing profiling of bacteria and fungi, we examined how root and soil microbiomes differ between benign and nematode-infected plant individuals in a soybean field in Japan. Results indicated that various taxonomic groups of bacteria and fungi occurred preferentially on the soybean individuals infected by root-knot nematodes or those uninfected by nematodes. Based on a network analysis of potential microbe-microbe associations, we further found that several fungal taxa potentially preying on nematodes (Dactylellina (Orbiliales), Rhizophydium (Rhizophydiales), Clonostachys (Hypocreales), Pochonia (Hypocreales) and Purpureocillium (Hypocreales)) co-occurred in the soybean rhizosphere at a small spatial scale. This study suggests how 'consortia' of anti-nematode microbes can derive from indigenous (resident) microbiomes, providing basic information for managing anti-nematode microbial communities in agroecosystems.

Factors Influencing Leaf- and Root-Associated Communities of Bacteria and Fungi Across 33 Plant Orders in a Grassland.

In terrestrial ecosystems, plants interact with diverse taxonomic groups of bacteria and fungi in the phyllosphere and rhizosphere. Although recent studies based on high-throughput DNA sequencing have drastically increased our understanding of plant-associated microbiomes, we still have limited knowledge of how plant species in a species-rich community differ in their leaf and root microbiome compositions. In a cool-temperate semi-natural grassland in Japan, we compared leaf- and root-associated microbiomes across 137 plant species belonging to 33 plant orders. Based on the whole-microbiome inventory data, we analyzed how sampling season as well as the taxonomy, nativeness (native or alien), lifeform (herbaceous or woody), and mycorrhizal type of host plants could contribute to variation in microbiome compositions among co-occurring plant species. The data also allowed us to explore prokaryote and fungal lineages showing preferences for specific host characteristics. The list of microbial taxa showing significant host preferences involved those potentially having some impacts on survival, growth, or environmental resistance of host plants. Overall, this study provides a platform for understanding how plant and microbial communities are linked with each other at the ecosystem level.

Leaf-associated microbiomes of grafted tomato plants.

Bacteria and fungi form complex communities (microbiomes) in above- and below-ground organs of plants, contributing to hosts' growth and survival in various ways. Recent studies have suggested that host plant genotypes control, at least partly, plant-associated microbiome compositions. However, we still have limited knowledge of how microbiome structures are determined in/on grafted crop plants, whose above-ground (scion) and below-ground (rootstock) genotypes are different with each other. By using eight varieties of grafted tomato plants, we examined how rootstock genotypes could determine the assembly of leaf endophytic microbes in field conditions. An Illumina sequencing analysis showed that both bacterial and fungal community structures did not significantly differ among tomato plants with different rootstock genotypes: rather, sampling positions in the farmland contributed to microbiome variation in a major way. Nonetheless, a further analysis targeting respective microbial taxa suggested that some bacteria and fungi could be preferentially associated with particular rootstock treatments. Specifically, a bacterium in the genus Deinococcus was found disproportionately from ungrafted tomato individuals. In addition, yeasts in the genus Hannaella occurred frequently on the tomato individuals whose rootstock genotype was "Ganbarune". Overall, this study suggests to what extent leaf microbiome structures can be affected/unaffected by rootstock genotypes in grafted crop plants.

Timing of evolutionary innovation: scenarios of evolutionary diversification in a species-rich fungal clade, Boletales.

Acquisition of mutualistic symbiosis could provide hosts and/or symbionts with novel ecological opportunities for evolutionary diversification. Such a mechanism is one of the major components of coevolutionary diversification. However, whether the origin of mycorrhizal symbiosis promotes diversification in fungi still requires clarification. Here, we aimed to reveal evolutionary diversification in a clade comprising ectomycorrhizal (ECM) fungi. Based on a phylogenic tree inferred from the sequences of 87 single-copy genes, we reconstructed the origins of ECM symbiosis in a species-rich basidiomycetous order, Boletales. High-resolution phylogeny of Boletales revealed that ECM symbiosis independently evolved from non-ECM states at least four times in the group. Among them, only the second most recent event, occurring in the clade of Boletaceae, was inferred to involve an almost synchronous rapid diversification and rapid transition from non-ECM to ECM symbiosis. Our results contradict the hypothesis of evolutionary priority effect, which postulates the greatest ecological opportunities in the oldest lineages. Therefore, the novel resources that had not been pre-empted by the old ECM fungal lineages - supposedly the coevolving angiosperm hosts - could be available for the young ECM fungal lineages, which resulted in evolutionary diversification occurring only in the young ECM fungal lineages.

Mycorrhizal fungi mediate the direction and strength of plant-soil feedbacks differently between arbuscular mycorrhizal and ectomycorrhizal communities.

Plants influence their soil environment, which affects the next generation of seedlings that can be established. While research has shown that such plant-soil feedbacks occur in the presence of mycorrhizal fungi, it remains unclear when and how mycorrhizal fungi mediate the direction and strength of feedbacks in tree communities. Here we show that arbuscular mycorrhizal and ectomycorrhizal fungal guilds mediate plant-soil feedbacks differently to influence large-scale patterns such as tree species coexistence and succession. When seedlings are grown under the same mycorrhizal type forest, arbuscular mycorrhizal plant species exhibit negative or neutral feedbacks and ectomycorrhizal plant species do neutral or positive feedbacks. In contrast, positive and neutral feedbacks dominate when seedlings are grown in associations within the same versus different mycorrhizal types. Thus, ectomycorrhizal communities show more positive feedbacks than arbuscular mycorrhizal communities, potentially explaining why most temperate forests are ectomycorrhizal.

Structural diversity across arbuscular mycorrhizal, ectomycorrhizal, and endophytic plant-fungus networks.

BACKGROUND: Below-ground linkage between plant and fungal communities is one of the major drivers of terrestrial ecosystem dynamics. However, we still have limited knowledge of how such plant-fungus associations vary in their community-scale properties depending on fungal functional groups and geographic locations. METHODS: By compiling a high-throughput sequencing dataset of root-associated fungi in eight forests along the Japanese Archipelago, we performed a comparative analysis of arbuscular mycorrhizal, ectomycorrhizal, and saprotrophic/endophytic associations across a latitudinal gradient from cool-temperate to subtropical regions. RESULTS: In most of the plant-fungus networks analyzed, host-symbiont associations were significantly specialized but lacked "nested" architecture, which has been commonly reported in plant-pollinator and plant-seed disperser networks. In particular, the entire networks involving all functional groups of plants and fungi and partial networks consisting of ectomycorrhizal plant and fungal species/taxa displayed "anti-nested" architecture (i.e., negative nestedness scores) in many of the forests examined. Our data also suggested that geographic factors affected the organization of plant-fungus network structure. For example, the southernmost subtropical site analyzed in this study displayed lower network-level specificity of host-symbiont associations and higher (but still low) nestedness than northern localities. CONCLUSIONS: Our comparative analyses suggest that arbuscular mycorrhizal, ectomycorrhizal, and saprotrophic/endophytic plant-fungus associations often lack nested network architecture, while those associations can vary, to some extent, in their community-scale properties along a latitudinal gradient. Overall, this study provides a basis for future studies that will examine how different types of plant-fungus associations collectively structure terrestrial ecosystems.

Publisher Correction: Core microbiomes for sustainable agroecosystems.

Owing to a technical error, this Perspective was originally published without its received and accepted dates; the dates "Received: 31 December 2017; Accepted: 23 March 2018" have now been included in all versions.