Interdomain ecological networks between plants and microbes

Environmental DNA Metabarcoding Reveals the Influence of Environmental Heterogeneity on Microeukaryotic Plankton in the Offshore Waters of East China Sea

Microeukaryotic plankton are essential to marine food webs and biogeochemical cycles, with coastal seas playing a critical role in aquatic ecosystems. Understanding the diversity of microeukaryotic plankton, deciphering their community structure and succession patterns, and identifying the key factors influencing these dynamics remain central challenges in coastal ecology. In this study, we examine patterns of biodiversity, community structure, and co-occurrence using environmental DNA (eDNA)-based methods. Our results show a linear correlation between α-diversity and distance from the shore, with nutrient-related factors, especially inorganic nitrogen, being the primary determinants of the spatial distribution of plankton communities. Alternation of coastal habitat have shifted the succession patterns of coastal eukaryotic plankton communities from stochastic to deterministic processes. Additionally, our observations indicate that the topology and structure of eukaryotic plankton symbiotic patterns and networks are significantly influenced by environmental heterogeneity such as nutrients, which increase the vulnerability and decrease the stability of offshore ecological networks. Overall, our study demonstrates that the distribution of microeukaryotic plankton communities is influenced by factors related to environmental heterogeneity.

Forest-to-Cropland Conversion Reshapes Microbial Hierarchical Interactions and Degrades Ecosystem Multifunctionality at a National Scale

Conversion from natural lands to cropland, primarily driven by agricultural expansion, could significantly alter soil microbiome worldwide; however, influences of forest-to-cropland conversion on microbial hierarchical interactions and ecosystem multifunctionality have not been fully understood. Here, we examined the effects of forest-to-cropland conversion on intratrophic and cross-trophic microbial interactions and soil ecosystem multifunctionality and further disclosed their underlying drivers at a national scale, using Illumina sequencing combined with high-throughput quantitative PCR techniques. The forest-to-cropland conversion significantly changed the structure of soil microbiome (including prokaryotic, fungal, and protistan communities) while it did not affect its alpha diversity. Both intrakingdom and interkingdom microbial networks revealed that the intratrophic and cross-trophic microbial interaction patterns generally tended to be more modular to resist environmental disturbance introduced from forest-to-cropland conversion, but this was insufficient for the cross-trophic interactions to maintain stability; hence, the protistan predation behaviors were still disturbed under such conversion. Moreover, key soil microbial clusters were declined during the forest-to-cropland conversion mainly because of the increased soil total phosphorus level, and this drove a great degradation of the ecosystem multifunctionality (by 207%) in cropland soils. Overall, these findings comprehensively implied the negative effects of forest-to-cropland conversion on the agroecosystem, from microbial hierarchical interactions to ecosystem multifunctionality.

The development of the biological soil crust regulates the fungal distribution and the stability of fungal networks

The heterogeneous composition of fungi plays an indispensable role in the foundation of the multifunctionalities of ecosystems within drylands. The precise mechanisms that govern fluctuations in soil fungal assemblages in dryland ecosystems remain incompletely elucidated. In this study, biological soil crusts (biocrusts) at different successional stages in the Gurbantunggut Desert were used as substrates to examine the characteristics and driving factors that influence fungal abundance and community dynamics during biocrust development using qPCR and high-throughput sequencing of the ITS2 region. The findings showed that the physicochemical properties changed significantly with the development of biocrusts. In particular, total nitrogen increased 4.8 times, along with notable increases in ammonium, total phosphorus (2.1 times) and soil organic carbon (6.5 times). Initially, there was a rise in fungal abundance, which was subsequently followed by a decline as the biocrust developed, with the highest abundance detected in lichen crust (2.66 × 107 copies/g soil) and the lowest in bare sand (7.98 × 106 copies/g soil). Ascomycetes and Basidiomycetes emerged as dominant phyla, collectively forming 85% of the fungal community. As the biocrust developed, noticeable alterations occurred in fungal community compositions, resulting from changes in the relative proportions of Dothideomycetes, Lecanoromycetes and unclassified ascomycetes. Nitrogen, phosphorus, organic carbon content, and pH of biocrusts were identified as direct or indirect regulators of fungal abundance and community structure. The complexity of fungal networks increased as biocrusts developed as revealed by network analysis, but reduced in the stability of fungal communities within algal and lichen crusts. Keystone species within the fungal community also underwent changes as biocrust developed. These results suggested that shifts in interspecies relationships among fungi could further contribute to the variation in fungal communities during the development of biocrusts.

Interaction patterns and keystone taxa of bacterial and eukaryotic communities during sulfamethoxazole mineralization in lake sediment

Sulfamethoxazole (SMX) is a common antibiotic pollutant in aquatic environments, which is highly persistent under various conditions and significantly contributes to the spread of antibiotic resistance. Biodegradation is the major pathway to eliminate antibiotics in the natural environment. The roles of bacteria and eukaryotes in the biodegradation of antibiotics have received considerable attention; however, their successions and co-occurrence patterns during the biodegradation of antibiotics remain unexplored. In this study, 13C-labled SMX was amended to sediment samples from Zhushan Bay (ZS), West Shore (WS), and Gonghu Bay (GH) in Taihu Lake to explore the interplay of bacterial and eukaryotic communities during a 30-day incubation period. The cumulative SMX mineralization on day 30 ranged from 5.2 % to 19.3 %, which was the highest in WS and the lowest in GH. The bacterial community showed larger within-group interactions than between-group interactions, and the positive interactions decreased during incubation. However, the eukaryotic community displayed larger between-group interactions than within-group interactions, and the positive interactions increased during incubation. The proportion of negative interactions between bacteria and eukaryotes increased during incubation. Fifty genera (including 46 bacterial and 4 eukaryotic genera) were identified as the keystone taxa due to their dominance in the co-occurrence network and tolerance to SMX. The cumulative relative abundance of these keystone taxa significantly increased during incubation and was consistent with the SMX mineralization rate. These taxa closely cooperated and played vital roles in co-occurrence networks and microbial community interactions, signifying their crucial role in SMX mineralization. These findings broadened our understanding of the complex interactions of microorganisms under SMX exposure and their potential functions during SMX mineralization, providing valuable insights for in situ bioremediation.

High levels of cadmium altered soil archaeal activity, assembly, and co-occurrence network in volcanic areas

Soil microbial communities are essential to biogeochemical cycles. However, the responses of microorganisms in volcanic soil with high heavy metal levels remain poorly understood. Here, two areas with high levels of cadmium (Cd) from the same volcano were investigated to determine their archaeal composition and assembly. In this study, the Cd concentrations (0.32-0.38 mg/ kg) in the volcanic soils exceeded the standard risk screening values (GB15618-2018) and correlated with archaeal communities strongly (P < 0.05). Moreover, the area with elevated levels of Cd (periphery) exhibited a greater diversity of archaeal species, albeit with reduced archaeal activity, compared to the area with lower levels of Cd (center). Besides, stochastic processes mainly governed the archaeal communities. Furthermore, the co-occurrence network was simplest in the periphery. The proportion of positive links between taxa increased positively with Cd concentration. Moreover, four keystone taxa (all from the family Nitrososphaeraceae) were identified from the archaeal networks. In its entirety, this study has expanded our comprehension of the variations of soil archaeal communities in volcanic areas with elevated cadmium levels and serves as a point of reference for the agricultural development of volcanic soils in China.

Relationships among protozoa, bacteria and fungi in polycyclic aromatic hydrocarbon-contaminated soils

Soil bacterial and fungal communities play key roles in the degradation of organic contaminants, and their structure and function are regulated by bottom-up and top-down factors. Microbial ecological effects of polycyclic aromatic hydrocarbons (PAHs) and trophic interactions among protozoa and bacteria/fungi in PAH-polluted soils have yet to be determined. We investigated the trophic interactions and structure of the microbiome in PAH-contaminated wasteland and farmland soils. The results indicated that the total concentration of the 16 PAHs (∑PAHs) was significantly correlated with the Shannon index, NMDS1 and the relative abundances of bacteria, fungi and protozoa (e.g., Pseudofungi) in the microbiome. Structural equation modelling and linear fitting demonstrated cascading relationships among PAHs, protozoan and bacterial/fungal communities in terms of abundance and diversity. Notably, individual PAHs were significantly correlated with microbe-grazing protozoa at the genus level, and the abundances of these organisms were significantly correlated with those of PAH-degrading bacteria and fungi. Bipartite networks and linear fitting indicated that protozoa indirectly modulate PAH degradation by regulating PAH-degrading bacterial and fungal communities. Therefore, protozoa might be involved in regulating the microbial degradation of PAHs by predation in contaminated soil.

Protists regulate microbially mediated organic carbon turnover in soil aggregates

Soil protists, the major predator of bacteria and fungi, shape the taxonomic and functional structure of soil microbiome via trophic regulation. However, how trophic interactions between protists and their prey influence microbially mediated soil organic carbon turnover remains largely unknown. Here, we investigated the protistan communities and microbial trophic interactions across different aggregates-size fractions in agricultural soil with long-term fertilization regimes. Our results showed that aggregate sizes significantly influenced the protistan community and microbial hierarchical interactions. Bacterivores were the predominant protistan functional group and were more abundant in macroaggregates and silt + clay than in microaggregates, while omnivores showed an opposite distribution pattern. Furthermore, partial least square path modeling revealed positive impacts of omnivores on the C-decomposition genes and soil organic matter (SOM) contents, while bacterivores displayed negative impacts. Microbial trophic interactions were intensive in macroaggregates and silt + clay but were restricted in microaggregates, as indicated by the intensity of protistan-bacterial associations and network complexity and connectivity. Cercozoan taxa were consistently identified as the keystone species in SOM degradation-related ecological clusters in macroaggregates and silt + clay, indicating the critical roles of protists in SOM degradation by regulating bacterial and fungal taxa. Chemical fertilization had a positive effect on soil C sequestration through suppressing SOM degradation-related ecological clusters in macroaggregate and silt + clay. Conversely, the associations between the trophic interactions and SOM contents were decoupled in microaggregates, suggesting limited microbial contributions to SOM turnovers. Our study demonstrates the importance of protists-driven trophic interactions on soil C cycling in agricultural ecosystems.

Changes in bacterial communities during rice cultivation remove phenolic constraints on peatland carbon preservation

Northern peatlands contain ~30% of terrestrial carbon (C) stores, but in recent decades, 14% to 20% of the stored C has been lost because of conversion of the peatland to cropland. Microorganisms are widely acknowledged as primary decomposers, but the keystone taxa within the bacterial community regulating C loss from cultivated peatlands remain largely unknown. In this study, we investigated the bacterial taxa driving peat C mineralization during rice cultivation. Cultivation significantly decreased concentrations of soil organic C, dissolved organic C (DOC), carbohydrates, and phenolics but increased C mineralization rate (CMR). Consistent with the classic theory that phenolic inhibition creates a "latch" that reduces peat C decomposition, phenolics were highly negatively correlated with CMR in cultivated peatlands, indicating that elimination of inhibitory phenolics can accelerate soil C mineralization. Bacterial communities were significantly different following peatland cultivation, and co-occurrence diagnosis analysis revealed substantial changes in network clusters of closely connected nodes (modules) and bacterial keystone taxa. Specifically, in cultivated peatlands, bacterial modules were significantly negatively correlated with phenolics, carbohydrates, and DOC. While keystone taxa Xanthomonadales, Arthrobacter, and Bacteroidetes_vadinHA17 can regulate bacterial modules and promote carbon mineralization. Those observations indicated that changes in bacterial modules can promote phenolic decomposition and eliminate phenolic inhibition of labile C decomposition, thus accelerating soil organic C loss during rice cultivation. Overall, the study provides deeper insights into microbe-driven peat C loss during rice cultivation and highlights the crucial role of keystone bacterial taxa in the removal of phenolic constraints on peat C preservation.

Geographical, climatic, and soil factors control the altitudinal pattern of rhizosphere microbial diversity and its driving effect on root zone soil multifunctionality in mountain ecosystems

Shifts in rhizosphere soil microorganisms of dominant plants' response to climate change profoundly impact mountain soil ecosystem multifunctionality; relatively little is known about the relationship between them and how they depend on long-term environmental drivers. Here, we conducted analyses of rhizosphere microbial altitudinal pattern, community assembly, and co-occurrence network of 6 dominant plants in six typical vegetation zones ranging from 1350 to 2900 m (a.s.l.) in Helan Mountains by absolute quantitative sequencing technology, and finally related the microbiomes to root zone soil multifunctionality ('soil multifunctionality' hereafter), the environmental dependence of the relationship was explored. It was found that the altitudinal pattern of rhizosphere soil bacterial and fungal diversities differed significantly. Higher co-occurrence and more potential interactions of Stipa breviflora and Carex coninux were found at the lowest and highest altitudes. Bacterial α diversity, the identity of some dominant bacterial and fungal taxa, had significant positive or negative effects on soil multifunctionality. The effect sizes of positive effects of microbial diversity on soil multifunctionality were greater than those of negative effects. These results indicated that the balance of positive and negative effects of microbes determines the impact of microbial diversity on soil multifunctionality. As the number of microbes at the phylum level increases, there will be a net gain in soil multifunctionality. Our study reveals that geographical and climatic factors can directly or modulate the effects of soil properties on rhizosphere microbial diversity, thereby affecting the driving effect of microbial diversity on soil multifunctionality, and points to the rhizosphere bacterial diversity rather than the fungi being strongly associated with soil multifunctionality. This work has important ecological implications for predicting how multiple environment-plant-soil-microorganisms interactions in mountain ecosystems will respond to future climate change.

Affinity difference determines the assembly and interaction mode of anammox community reconstructed by siderophores

Anammox community is the core of anammox process. The constancy of the anammox community determines the stability of the anammox process and the ability of withstand environmental impact. Community stability is influenced by the assembly and interaction mode of the community. This study aimed to explore the assembly, interaction mode, and stability of anammox community influenced by two siderophores (enterobactin and putrebactin) specific for Ca. Brocadia and Ca. Kuenenia as produced in our previous research. Siderophores improved the stability of the anammox community, among which vulnerability dropped by 30.02 % and 72.53 % respectively. Enterobactin and putrebactin altered the succession speed and assembly pattern of communities, with a respective increase of 9.77 % and 80.87 % in the deterministic process of anammox community assembly, respectively. Enterobactin and putrebactin reduced the dependence of Ca. Brocadia and Ca. Kuenenia on companion bacteria by 60 items and 27 items respectively. The affinity of different siderophore-Fe with bacterial membrane receptors caused variations in community reconstruction, with Ca. Brocadia and Ca. Kuenenia exhibiting the highest affinity with enterobactin-Fe (-11.4 kcal/mol) and putrebactin-Fe (-9.0 kcal/mol), respectively. This study demonstrated how siderophores can enhance the stability of anammox process by regulating assembly and interaction mode of anammox community, while also revealing the underlying molecular mechanisms.

Interkingdom network analyses reveal microalgae and protostomes as keystone taxa involved in nutrient cycling in large freshwater lake sediment

Few studies have explored the role of interkingdom interactions between bacteria and microeukaryotes in nutrient cycling in lake ecosystems. We conducted sediment sampling from 40 locations covering Hongze Lake and analyzed their chemical properties. Intra- and interkingdom networks were constructed using 16S and 18S rRNA gene amplicon sequencing. Microeukaryotic intranetworks were more complex in spring than in autumn, while no clear variation in the complexity of bacterial intranetworks was found between autumn and spring. Larger and more complex bacterial-microeukaryotic bipartite networks emerged in spring than in autumn, correlated with lower carbon, nitrogen, and phosphorus levels in spring, likely resulting in intense microbial competition. Bacteria and microeukaryotes played different topological roles in interkingdom networks, with microeukaryotes contributing to the networks' greater complexity. Seven keystone modules were identified in spring and autumn nutrient cycling. Importantly, keystone taxa in these modules belonged to photoautotrophic microalgae or predatory protostomes, indicating that these organisms are key drivers in lake sediment nutrient cycling. Our results suggested that nutrient content variation in autumn and spring changes interkingdom networks' topological structure between bacteria and microeukaryotes. Microalgae and protostomes are essential in freshwater lake nutrient cycling and may be targeted to modulate nutrient cycling in large freshwater ecosystems.

Pedogenesis of typical zonal soil drives belowground bacterial communities of arable land in the Northeast China Plain

Belowground bacterial communities play essential roles in maintaining ecosystem multifunction, while our understanding of how and why their distribution patterns and community compositions may change with the distinct pedogenetic conditions of different soil types is still limited. Here, we evaluated the roles of soil physiochemical properties and biotic interactions in driving belowground bacterial community composition across three typical zonal soil types, including black calcium soil (QS), typical black soil (HL) and dark brown soil (BQL), with distinct pedogenesis on the Northeast China Plain. Changes in soil bacterial diversity and community composition in these three zonal soil types were strongly correlated with soil pedogenetic features. SOC concentrations in HL were higher than in QS and BQL, but bacterial diversity was low, and the network structure revealed greater stability and connectivity. The composition of the bacterial community correlated significantly with soil pH in QS but with soil texture in BQL. The bacterial co-occurrence network of HL had higher density and clustering coefficients but lower edges, and different keystone species of networks were also detected. This work provides a basic understanding of the driving mechanisms responsible for belowground bacterial biodiversity and distribution patterns over different pedogenetic conditions in agroecosystems.

Inoculation with thermophiles enhanced the food waste bio-drying and complicated interdomain ecological networks between bacterial and fungal communities

Bio-drying is a practical approach for treating food waste (FW). However, microbial ecological processes during treatment are essential for improving the dry efficiency, and have not been stressed enough. This study analyzed the microbial community succession and two critical periods of interdomain ecological networks (IDENs) during FW bio-drying inoculated with thermophiles (TB), to determine how TB affects FW bio-drying efficiency. The results showed that TB could rapidly colonize in the FW bio-drying, with the highest relative abundance of 5.13%. Inoculating TB increased the maximum temperature, temperature integrated index and moisture removal rate of FW bio-drying (55.7 °C, 219.5 °C, and 86.11% vs. 52.1 °C, 159.1 °C, and 56.02%), thereby accelerating the FW bio-drying efficiency by altering the succession of microbial communities. The structural equation model and IDEN analysis demonstrated that TB inoculation complicated the IDENs between bacterial and fungal communities by significantly and positively affecting bacterial communities (b = 0.39, p < 0.001) and fungal communities (b = 0.32, p < 0.01), thereby enhancing interdomain interactions between bacteria and fungi. Additionally, inoculation TB significantly increased the relative abundance of keystone taxa, including Clostridium sensu stricto, Ochrobactrum, Phenylobacterium, Microvirga and Candida. In conclusion, the inoculation of TB could effectively improve FW bio-drying, which is a promising technology for rapidly reducing FW with high moisture content and recovering resources from it.

Synergistic effect of biotic and abiotic factors drives microbiota succession and assembly in medium-temperature Daqu

BACKGROUND

The feasibility of fortification techniques to improve the quality of medium-temperature Daqu (MTD) by inoculation functional isolates has been demonstrated. However, it is unclear what is the effect of inoculation on the controllability during the MTD fermentation process. Here, inoculated a single strain of Bacillus licheniformis, and the microbiota composed of Bacillus velezensis and Bacillus subtilis, were used to investigate the synergistic effect of biotic and abiotic factors on the succession and assembly of the MTD microbiota during the process.

RESULTS

The biotic factors promoted the proliferation of microorganisms that arrived early at the MTD. Subsequently, this alteration might inhibit microorganisms that colonized later in the MTD microecosystem, thereby assembling a different but more stable microbial community. Moreover, the biotic factors making bacterial community assembly were dominated by variable selection earlier, whereas the fungal community assembly was dominated mainly by extreme abiotic factors rather than biotic factors. Interestingly, fermentation temperature and moisture were significantly associated with the succession and assembly of the fortified MTD community. Meanwhile, the effect of the environmental variables on endogenous variables was also significant. Thus, changes in endogenous variables could be mitigated by adjusting environmental variables to regulate the process of MTD fermentation.

CONCLUSION

Biotic factors cause rapid changes of the microbiota during the MTD fermentation process, which could be controlled indirectly by regulating environmental variables. Meanwhile, a more stable MTD ecological network might be beneficial for enhancing the stability of MTD quality. © 2023 Society of Chemical Industry.

Resource-dependent biodiversity and potential multi-trophic interactions determine belowground functional trait stability

BACKGROUND

For achieving long-term sustainability of intensive agricultural practices, it is pivotal to understand belowground functional stability as belowground organisms play essential roles in soil biogeochemical cycling. It is commonly believed that resource availability is critical for controlling the soil biodiversity and belowground organism interactions that ultimately lead to the stabilization or collapse of terrestrial ecosystem functions, but evidence to support this belief is still limited. Here, we leveraged field experiments from the Chinese National Ecosystem Research Network (CERN) and two microcosm experiments mimicking high and low resource conditions to explore how resource availability mediates soil biodiversity and potential multi-trophic interactions to control functional trait stability.

RESULTS

We found that agricultural practice-induced higher resource availability increased potential cross-trophic interactions over 316% in fields, which in turn had a greater effect on functional trait stability, while low resource availability made the stability more dependent on the potential within trophic interactions and soil biodiversity. This large-scale pattern was confirmed by fine-scale microcosm systems, showing that microcosms with sufficient nutrient supply increase the proportion of potential cross-trophic interactions, which were positively associated with functional stability. Resource-driven belowground biodiversity and multi-trophic interactions ultimately feedback to the stability of plant biomass.

CONCLUSIONS

Our results indicated the importance of potential multi-trophic interactions in supporting belowground functional trait stability, especially when nutrients are sufficient, and also suggested the ecological benefits of fertilization programs in modern agricultural intensification. Video Abstract.

Electrostimulation triggers an increase in cross-niche microbial associations toward enhancing organic nitrogen wastewater treatment

As an efficient wastewater pretreatment biotechnology, electrostimulated hydrolysis acidification (eHA) has been used to accelerate the removal of refractory pollutants, which is closely related to the effects of electrostimulation on microbial interspecies associations. However, the ecological processes underpinning such linkages remain unresolved, especially for the microbial communities derived from different niches, such as the electrode surface and plankton. Herein, the principles of cross-niche microbial associations and community assembly were investigated using molecular ecological network and phylogenetic bin-based null model analysis (iCAMP) based on 16S rRNA gene sequences. The electrostimulated planktonic sludge and electrode biofilm displayed significantly (P < 0.05) 1.67 and 1.53 times higher organic nitrogen pollutant (azo dye Alizarin Yellow R) degradation efficiency than non-electrostimulation group, and the corresponding microbial community composition and structure were significantly (P < 0.05) changed. Electroactive bacteria and functional degraders were enriched in the electrode biofilm and planktonic sludge, respectively. Notably, electrostimulation strengthened the synergistic microbial associations (1.8 times more links) between sludge and biofilm members. Additionally, both electrostimulation and cross-niche microbial associations induced greater importance of deterministic assembly. Overall, this study highlights the specificity of cross-electrode surface microbial associations and ecological processes with electrostimulation and advances our understanding of the manipulation of sludge microbiomes in engineered wastewater treatment systems.

Deciphering microeukaryotic-bacterial co-occurrence networks in coastal aquaculture ponds

Microeukaryotes and bacteria are key drivers of primary productivity and nutrient cycling in aquaculture ecosystems. Although their diversity and composition have been widely investigated in aquaculture systems, the co-occurrence bipartite network between microeukaryotes and bacteria remains poorly understood. This study used the bipartite network analysis of high-throughput sequencing datasets to detect the co-occurrence relationships between microeukaryotes and bacteria in water and sediment from coastal aquaculture ponds. Chlorophyta and fungi were dominant phyla in the microeukaryotic-bacterial bipartite networks in water and sediment, respectively. Chlorophyta also had overrepresented links with bacteria in water. Most microeukaryotes and bacteria were classified as generalists, and tended to have symmetric positive and negative links with bacteria in both water and sediment. However, some microeukaryotes with high density of links showed asymmetric links with bacteria in water. Modularity detection in the bipartite network indicated that four microeukaryotes and twelve uncultured bacteria might be potential keystone taxa among the module connections. Moreover, the microeukaryotic-bacterial bipartite network in sediment harbored significantly more nestedness than that in water. The loss of microeukaryotes and generalists will more likely lead to the collapse of positive co-occurrence relationships between microeukaryotes and bacteria in both water and sediment. This study unveils the topology, dominant taxa, keystone species, and robustness in the microeukaryotic-bacterial bipartite networks in coastal aquaculture ecosystems. These species herein can be applied for further management of ecological services, and such knowledge may also be very useful for the regulation of other eutrophic ecosystems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00159-6.

Grass-microbial inter-domain ecological networks associated with alpine grassland productivity

Associations between grasses and soil microorganisms can strongly influence plant community structures. However, the associations between grass productivity and diversity and soil microbes, as well as the patterns of co-occurrence between grass and microbes remain unclear. Here, we surveyed grass productivity and diversity, determined soil physicochemical, and sequenced soil archaea, bacteria and fungi by metabarcoding technology at 16 alpine grasslands. Using the Distance-decay relationship, Inter-Domain Ecological Network (IDEN), and Mantel tests, we investigated the relationship between grass productivity, diversity and microbial diversity, and the patterns of co-occurrence between grass and microbial inter-domain network in alpine grassland. We found the archaea richness, bacteria richness and Shannon, and fungi α-diversity were significantly negatively correlation with grass diversity, but archaea and bacteria diversity were positively correlation with grass productivity. Moreover, an increase in microbial β-diversity was observed along with increased discrepancy in grass diversity and productivity and soil variables. Variance partitioning analysis suggested that the contribution of grass productivity on microbial community was higher than that of soil variables and grass diversity, which implies that microbial community was more related to grass productivity. Inter-Domain Ecological Network showed that the grass species formed complex and stable ecological networks with some bacterial, archaeal, and fungal species, and the grass-fungal ecological networks showed the highest robustness, which indicated that soil fungi could better co-coexist with aboveground grass in alpine grasslands. Besides, the connectivity degrees of the grass-microbial network were significantly positively correlated with grass productivity, suggesting that the coexistence pattern of grasses and microbes had a positive feedback effect on the grass productivity. The results are important for establishing the regulatory mechanisms between plants and microorganisms in alpine grassland ecosystems.

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.

Fungi and cercozoa regulate methane-associated prokaryotes in wetland methane emissions

Wetlands are natural sources of methane (CH₄) emissions, providing the largest contribution to the atmospheric CH₄ pool. Changes in the ecohydrological environment of coastal salt marshes, especially the surface inundation level, cause instability in the CH₄ emission levels of coastal ecosystems. Although soil methane-associated microorganisms play key roles in both CH₄ generation and metabolism, how other microorganisms regulate methane emission and their responses to inundation has not been investigated. Here, we studied the responses of prokaryotic, fungal and cercozoan communities following 5 years of inundation treatments in a wetland experimental site, and molecular ecological networks analysis (MENs) was constructed to characterize the interdomain relationship. The result showed that the degree of inundation significantly altered the CH₄ emissions, and the abundance of the pmoA gene for methanotrophs shifted more significantly than the mcrA gene for methanogens, and they both showed significant positive correlations to methane flux. Additionally, we found inundation significantly altered the diversity of the prokaryotic and fungal communities, as well as the composition of key species in interactions within prokaryotic, fungal, and cercozoan communities. Mantel tests indicated that the structure of the three communities showed significant correlations to methane emissions (p < 0.05), suggesting that all three microbial communities directly or indirectly contributed to the methane emissions of this ecosystem. Correspondingly, the interdomain networks among microbial communities revealed that methane-associated prokaryotic and cercozoan OTUs were all keystone taxa. Methane-associated OTUs were more likely to interact in pairs and correlated negatively with the fungal and cercozoan communities. In addition, the modules significantly positively correlated with methane flux were affected by environmental stress (i.e., pH) and soil nutrients (i.e., total nitrogen, total phosphorus and organic matter), suggesting that these factors tend to positively regulate methane flux by regulating microbial relationships under inundation. Our findings demonstrated that the inundation altered microbial communities in coastal wetlands, and the fungal and cercozoan communities played vital roles in regulating methane emission through microbial interactions with the methane-associated community.

Freshwater trophic status mediates microbial community assembly and interdomain network complexity

Freshwater microorganisms and their interactions are important drivers of nutrient cycling that are in turn affected by nutrient status, causing shifts in microbial community diversity, composition, and interactions. However, the impact of water trophic status on bacterial-archaeal interdomain interactions remains poorly understood. This study focused on the impact of trophic status, as characterized by trophic state index (TSI), on the interdomain interactions of freshwater microbial communities from 45 ponds in Hangzhou. Our results showed that the mesotrophic wetland bordering on lightly eutrophic (Hemu: TSI of 49; lightly eutrophic is defined as 50 ≤ TSI <60) harbored a much more complex bacterial-archaeal interdomain network, which showed significantly (P < 0.05) higher connectivity than the wetlands with lower (TSI of 38) or higher (TSI of 57) trophic levels. Notably, light eutrophication strengthened the network modules' negative associations with organic carbon through some network hubs, which could trigger carbon loss in wetlands. We also detected a non-linear response of interdomain network complexity to the increasing of nutrients with a turning point of approximately TSI 50. Quantitative estimates of community assembly processes and structural equation modelling analysis indicated that chlorophyll-a, total nitrogen, and total phosphorus could regulate interdomain network complexity (50% of the variation explanation rate) by driving microbial community assembly. This study demonstrates that microbial interdomain network complexity could be used as a bioindicator for ecological changes, which would helpful for improving ecological assessment of the freshwater eutrophication.

Community response of soil microorganisms to combined contamination of polycyclic aromatic hydrocarbons and potentially toxic elements in a typical coking plant

Both polycyclic aromatic hydrocarbons (PAHs) and potentially toxic elements (PTEs) of coking industries impose negative effects on the stability of soil ecosystem. Soil microbes are regarded as an essential moderator of biochemical processes and soil remediation, while their responses to PAHs-PTEs combined contamination are largely unknown. In the present study, soil microbial diversity and community composition in the typical coking plant under the chronic co-exposure of PAHs and PTEs were investigated and microbial interaction networks were built to reveal microbial co-occurrence patterns. The results indicated that the concentrations of PAHs in the soil inside the coking plant were significantly higher than those outside the plant. The mean concentration of ∑16PAHs was 2894.4 ng·g^-1, which is 5.58 times higher than that outside the plant. The average Hg concentration inside the coking plant was 22 times higher than the background value of Hebei province. The soil fungal community inside the coking plant showed lower richness compared with that of outside community, and there are significant difference in the bacterial and fungal community composition between inside and outside of coking plant (p < 0.01). Predicted contribution of different environmental factors to each dominant species based on random forest identified 20 and 25 biomarkers in bacteria and fungi, respectively, that were highly sensitive to coking plant soil in operation, such as Betaproteobacteria，Sordariomycetes and Dothideomycetes. Bacterial and fungal communities were shaped by the soil chemical properties (pH), PTEs (Hg), and PAHs together in the coking plant soils. Furthermore, the bacterial and fungal interaction patterns were investigated separately or jointly by intradomain and interdomain networks. Competition is the main strategy based on the co-exclusion pattern in fungal community, and the competitive relationship inside the coking plant is more complex than that outside the plant. In contrast, cooperation is the dominant strategy in bacterial networks based on the co-occurrence pattern. The present study provided insights into microbial response strategies and the interactions between bacteria and fungi under long-term combined contamination.

Microbial network complexity drives non-linear shift in biodiversity-nutrient cycling in a saline urban reservoir

Aquatic biodiversity is important in mediating ecosystem functioning, contributing to ecosystem sustainability and human wellbeing. However, how microbial network complexity affects the biodiversity-nutrient cycling relationship in saline freshwater ecosystems remains underexplored. Using high-resolution time-series data, we examined the relationships between microeukaryotic-bacterial community network complexity, biodiversity and multi-nutrient cycling in an urban reservoir undergoing a freshwater salinization-desalinization cycle. We found that low microbial diversity enhanced ecosystem multi-nutrient cycling under high salinity stress. In addition, multi-nutrient cycling declined with increased network complexity. Further, we found a non-linear relationship between salinity-induced shifts in the complexity of the microbial network and biodiversity-nutrient cycling (BNC) relationship of keystone taxa, i.e. the strength of the BNC relationship first became weak and then strong with increased network complexity. Together, these results highlighted the significant insight that there is not always positive relationship between biodiversity/network complexity and multi-nutrient cycling, even between network complexity and BNC relationship in real-world ecosystems, suggesting that preserving microbial association is important in aquatic health managing and evaluating the freshwater salinization problem.

The Stability of Phyto-Zooplanktonic Networks Varied with Zooplanktonic Sizes in Chinese Coastal Ecosystem

The linkages between phytoplankton and zooplankton are crucial for the stability of complex food webs and the flow of energy within the marine ecosystem. Despite body size exhibiting multiple effects on the planktonic community assembly and the dispersal scale, its role in determining the stability of phyto-zooplanktonic co-occurrence patterns remains unclear. Here, we focused on more than 13,000 kilometers of the Chinese coast to study the diatom-dominated plankton ecosystem and to report the significant negative effects of zooplanktonic body sizes on the topological properties of phyto-zooplanktonic networks (PZNs) by using more than 500 species from 251 samples taken along the coastline. PZNs tended to be more complex and stable when phytoplankton associated with smaller zooplankton. Particularly, the subnetworks of dominant phytoplankton displayed differences with different zooplanktonic body sizes. The zooplankton with larger and smaller body sizes tended to interact with dinoflagellates and diatoms, respectively. Additionally, abiotic factors (i.e., water temperature, pH, salinity, and metal concentrations) displayed significant effects on PZNs via the shifting of zooplanktonic composition, and the zooplanktonic body sizes altered the network modules' associations with different environmental factors. Our study elucidated the general relationship between zooplanktonic body sizes and the stability of PZNs, which provides new insights into marine food webs. IMPORTANCE Body size is a key life trait of aquatic plankton that affects organisms' metabolic rates and ecological functions; however, its specific effects on interactions between phytoplankton and zooplankton are poorly understood. We collected planktonic species and their body size data along more than 13,000 kilometers of coastline to explore the role of zooplanktonic body size in maintaining the stability of phyto-zooplanktonic networks (PZNs). We found that zooplankton play a more important role in maintaining PZN stability than do phytoplankton as well as that the PZN would be more complex and stable with smaller zooplankton. Furthermore, this work revealed that body size significantly determined the relationships between environmental factors and network structure. Overall, these findings lay a general relationship between zooplanktonic body sizes and the stability of PZNs, which helps us further explore the micro food web of coastal ecosystems.

Soil microbial community assembly and stability are associated with potato (Solanum tuberosum L.) fitness under continuous cropping regime

Continuous cropping obstacles caused by the over-cultivation of a single crop trigger soil degradation, yield reduction and the occurrence of plant disease. However, the relationships among stability, complexity and the assembly process of soil microbial community with continuous cropping obstacles remains unclear. In this study, molecular ecological networks analysis (MENs) and inter-domain ecological networks analysis (IDENs), and a new index named cohesion tools were used to calculate the stability and complexity of soil microbial communities from eight potato cultivars grown under a continuous cropping regime by using the high-throughput sequencing data. The results showed that the stability (i.e., robustness index) of the bacterial and fungal communities for cultivar ZS5 was significantly higher, and that the complexity (i.e., cohesion values) was also significantly higher in the bacterial, fungal and inter-domain communities (i.e., bacterial-fungal community) of cultivar ZS5 than other cultivars. Network analysis also revealed that Actinobacteria and Ascomycota were the dominant phyla within intra-domain networks of continuous cropping potato soil communities, while the phyla Proteobacteria and Ascomycota dominated the correlation of the bacterial-fungal network. Infer community assembly mechanism by phylogenetic-bin-based null model analysis (iCAMP) tools were used to calculate the soil bacterial and fungal communities' assembly processes of the eight potato cultivars under continuous cropping regime, and the results showed that the bacterial community was mainly dominated by deterministic processes (64.19% - 81.31%) while the fungal community was mainly dominated by stochastic processes (78.28% - 98.99%), indicating that the continuous-cropping regime mainly influenced the potato soil bacterial community assembly process. Moreover, cultivar ZS5 possessed a relatively lower homogeneous selection, and a higher TP, TN, AP and yield than other cultivars. Our results indicated that the soil microbial network stability and complexity, and community assemble might be associated with yield and soil properties, which would be helpful in the study for resistance to potato continuous cropping obstacles.

Microbial cooperation promotes humification to reduce antibiotic resistance genes abundance in food waste composting

Antibiotic resistance genes (ARGs) fate in a full-scale Food waste (FW) facility was investigated. Results showed that with the changes in ARGs, microbial networks could be naturally divided into two clusters, named as the ARGs increasing group (AI group) and the ARGs decreasing group (AD group). The significant difference between two groups (i.e. stronger microbial competition in the AI group and stronger microbial cooperation in the AD group) implied that the variation in ARGs over time were caused by a switch between competition and cooperation. These results indicated that microbial competition might increase ARGs abundance, while cooperation might reduce it. Meanwhile, structural-equation-model (SEM model) showed that humification indexes (e.g. GI value) was an indicator for characterizing microbial interactions and ARGs. The results of the linear model further confirmed that mature compost (GI values > 92.6 %) could reduce the risk of ARGs.

Core fungal species strengthen microbial cooperation in a food-waste composting process

In ecosystem engineering research, the contribution of microbial cooperation to ecosystem function has been emphasized. Fungi are one of the predominant decomposers in composting, but thus far, less attention has been given to fungal than to bacterial cooperation. Therefore, network and cohesion analyses were combined to reveal the correlation between fungal cooperation and organic matter (OM) degradation in ten composting piles. Positive cohesion, reflecting the cooperation degree, was positively linked to the degradation rate of OM. From the community perspective, core species (i.e., Candida tropicalis, Issatchenkia orientails, Kazachstania exigua, and Dipodascus australiensis) with high occurrence frequency and abundance were the key in regulating positive cohesion. These species were highly relevant to functional genera associated with OM degradation in both fungal and bacterial domains. Therefore, focusing on these core fungal species might be an appropriate strategy for targeted regulation of functional microbes and promotion of degradation rates.

The wilt pathogen induces different variations of root-associated microbiomes of plant

Root-associated compartments, including the rhizosphere, rhizoplane, and endosphere, live with diverse microbial communities which profoundly affect plant growth and health. However, a systematic understanding of the microbiome assembly across the rhizosphere, rhizoplane, and endosphere under pathogen invasion remains elusive. Using 16S high-throughput sequencing, we studied how bacterial wilt disease affected the variation and assembly of the three continuous root-associated microbiomes of tobacco. The results indicated that microorganisms were gradually filtered from the rhizosphere to the endosphere. With the pathogen invasion, the rhizosphere, rhizoplane and endosphere microbiomes selected and recruited different beneficial bacterial taxa. Some recruited bacteria were also identified as keystone members in networks (i.e., Bosea in the endosphere). The microbiomes of endosphere and rhizoplane were more sensitive to plant disease than the rhizosphere microbiome. Still, response strategies of the rhizoplane and endosphere to disease were obviously different. Microbial networks of the rhizoplane became complex in diseased samples and genes involved in sporulation formation and cell cycle were enriched. However, microbial networks of the diseased endosphere were disrupted, and functional genes related to nitrogen utilization and chemotaxis were significantly increased, indicating the importance of nitrogen resources supply of plants for the endosphere microbiome under pathogen invasion. Our results provide novel insights for understanding the different responses of the root-associated microbiomes to plant disease.

Succession of function, assembly, and interaction of microbial community in sequencing biofilm batch reactors under selenite stress

The mechanism of interaction between selenite, a toxic substance, and the microbial community in wastewater is still not well understood. Herein, nine sequencing biofilm batch reactors were used to systematically investigate the response of the microbial community to the continuous selenite stress. The results showed that selenite affected the reactor performance and reduced the biofilm mass. Also, it increased the proportion of the living cells, and changed the protein and polysaccharide composition of the biofilm as well as cellular secretions. Selenite facilitated the removal of NO₃-N, according to water-quality and bioinformatics analyses. As such, the selenite was converted into selenium nanoparticles. α-diversity analysis further revealed that 20 μM selenite enhanced the microbial community resilience, while 200 μM selenite had the reverse effect. Community composition analysis showed that Variovorax, Rhizobium, and Simkania had positive correlations with selenite (P < 0.05). Functional prediction suggested that selenite changed the C, N, and S cycle functions. Furthermore, determinism dominated the community assembly process, and the deterministic proportion increased with the increase of selenite concentration. Network analysis showed that selenite improved the stability and positive correlation ratio of the overall microbial network, and accelerated the communication between microorganisms. However, when compared with the 20 μM selenite, the 200 μM selenite boosted the competition and parasitism/predation among microorganisms. Low-abundance genera played a key role in the network of selenite-reducing microbial community. In addition, under selenite stress, biofilm network exhibited better stability and faster information exchange than suspended network, and the positive association between biofilm and suspended microorganisms increased. All in all, this research sheds light on the interaction between selenite and microbial community, as well as provides crucial information on selenium-containing wastewater.

Succession of soil microbial community in a developing mid-channel bar: The role of environmental disturbance and plant community

Succession of microbial and plant communities is crucial for the development and the stability of soil ecological functions. The relative role of plant communities and environmental disturbance in shaping the microbial community in a newly established habitat remains unclear. In this study, a mid-channel bar (MCB) exposed to an environmental disturbance gradient in the Yangtze River was studied to explore the effects of such disturbance and plant community traits on the succession of the soil microbial community. Bulk and rhizospheric soils were collected from the MCB and classified according to their level of exposure to environmental disturbance: head, central and tail. These subsequently underwent high-throughput sequencing and interdomain ecological network (IDEN) analysis to identify and characterize the predominant microbial groups present in the soils at each disturbance level. Furthermore, at each site, the presence and distribution of the plant community was also noted. The present study demonstrated that both bulk soil nutrients and plant community exhibited significant spatial distribution dependent on the level of disturbance and this influenced the composition of the microbial community. In less eroded parts of the MCB, i.e., the central, nutrients accumulated, promoting growths of plants. This in turn encouraged a more diverse microbial community, dominated by the bacterial genus Pseudarthrobacter. Plant showed a stronger association with bulk soil microbial communities compared to rhizosphere soil microbial communities. Particularly, Triarrhena sacchariflora and Hemarthria altissima, present in sites of low disturbance, exhibiting a more extensive plant-microbe association. They thus played a key role in shaping the soil microbial community. In general, however, plant species did not directly determine the composition of the bacterial community, but instead altered the nutritive state of the soil to promote microbial growth. Such findings are of significant value for conservation practices of newly formed ecosystems, which requires an integrated understanding of the role of environmental disturbance and plants on soil microbial community assemblage.

Terrigenous dissolved organic matter persists in the energy-limited deep groundwaters of the Fennoscandian Shield

The deep terrestrial biosphere encompasses the life below the photosynthesis-fueled surface that perseveres in typically nutrient and energy depleted anoxic groundwaters. The composition and cycling of this vast dissolved organic matter (DOM) reservoir relevant to the global carbon cycle remains to be deciphered. Here we show that recent Baltic Sea-influenced to ancient pre-Holocene saline Fennoscandian Shield deep bedrock fracture waters carried DOM with a strong terrigenous signature and varying contributions from abiotic and biotic processes. Removal of easily degraded carbon at the surface-to-groundwater transition and corresponding microbial community assembly processes likely resulted in the highly similar DOM signatures across the notably different water types that selected for a core microbiome. In combination with the aliphatic character, depleted δ¹³C signatures in DOM indicated recent microbial production in the oldest, saline groundwater. Our study revealed the persistence of terrestrially-sourced carbon in severely energy limited deep continental groundwaters supporting deep microbial life.

Tracking the diversity and interaction of methanogens in the energy recovery process of a full-scale wastewater treatment plant

Methanogens have been significant for the achievement of carbon neutrality in wastewater treatment plants due to their crucial roles in the anaerobic digestion of sludge. Nevertheless, the phylogenetic diversity of methanogens and their versatile metabolism have been continuously investigated, the current scientific knowledge regarding these microbes appears inadequate and requires more evaluations. This study is considered an endeavor in which functional genes sequencing was used to reveal the diversity of methanogens in the sludge process of the wastewater treatment plant. The information obtained was substantially more than that employing 16s sequencing. The methanogenic microbial resources were appropriate to sustain a self-inoculated energy recovery with a potential ability to boost methane production. A constancy was observed in 16 S rRNA gene and mcrA gene sequencing results, where the bacterial or Methanosaeta concilii dominated community of DS (digest sludge) was distinct from the inoculum sources TS (total sludge), CTS (concentrated total sludge), and HTS (hydrolysis total sludge), indicating the independent development of DS. A quantitative cross-network was constructed by coupling the absolute quantify of 16 S rRNA and mcrA sequences. The Methanobacterium petrolearium actively interacted with bacteria in the DS community rather than the dominant species (Methanosaeta concilii). Moreover, the unclassified methanogens were identified to be significantly prevalent in all communities, suggesting that unknown methanogenic taxa might be imperative in accomplishing community functions. Collectively, the findings of this research study will shed light on the comprehensive knowledge of microbial communities, especially the methanogenic microbiota. This will further enhance the exploration of the phylogenetic diversity of methanogens and their corresponding impacts in energy recovery from wastewater treatment plants.

Interkingdom Plant–Soil Microbial Ecological Network Analysis under Different Anthropogenic Impacts in a Tropical Rainforest

Plants and their associated soil microorganisms interact with each other and form complex relationships. The effects of slash-and-burn agriculture and logging on aboveground plants and belowground microorganisms have been extensively studied, but research on plant–microbial interkingdom ecological networks is lacking. In this study, using old growth forest as a control, we used metagenomic data (ITS and 16S rRNA gene amplified sequences) and plant data to obtain interdomain species association patterns for three different soil disturbance types (slash-and-burn, clear cutting and selective cutting) in a tropical rainforest based on interdomain ecological network (IDEN) analysis. Results showed that the soil bacterial–fungal and plant–microbe ecological networks had different topological properties among the three forest disturbance types compared to old growth forest. More nodes, links, higher modularity and negative proportion were found in the selective cutting stand, indicating higher stability with increasing antagonistic relationships and niche differentiation. However, the area of slash-and-burn forest yield opposite results. Network module analysis indicated that different keystone species were found in the four forest types, suggesting alternative stable states among them. Different plant species had more preferential associations with specific fungal taxa than bacterial taxa at the genus level and plant–microbe associations lagged behind bacterial–fungal associations. Overall, compared with old growth forests, the bacterial–fungal and plant–microbe ecological networks in the slash-and-burn and clear cutting stands were simpler, while the network in the selective cutting stand was more complex. Understanding the relationships between aboveground plants and belowground microorganisms under differing disturbance patterns in natural ecosystems will help in better understanding the surrounding ecosystem functions of ecological networks.

The large-scale spatial patterns of ecological networks between phytoplankton and zooplankton in coastal marine ecosystems

Although autotrophic phytoplankton and heterotrophic zooplankton both play important roles in the food web of marine ecosystem, their comprehensive interactions and spatial patterns at continental scale remain poorly studied. Here, we collected 251 seawater samples along 13,000 km of Chinese coastline, and microscopically investigated the latitudinal gradients of planktonic diversities. In total, 307 phytoplanktonic and 311 zooplanktonic species were visually identified. Using the newly developed Inter-Domain Ecological Networks (IDENs) approach, the phytoplankton-zooplankton interaction networks were constructed. We found that the phyto-zooplankton network structure was varied across three regions, more complex and numerous connections along the southern coast than in the north. In addition, some particular associations between zooplanktonic and phytoplanktonic groups were found to be localized in specific regions. Furthermore, the seawater temperature and salinity were the major driving force for shaping planktonic interaction networks. These results provide a deeper understanding of planktonic biogeography and phytoplankton-zooplankton interaction patterns.

iNAP: An integrated network analysis pipeline for microbiome studies

Integrated network analysis pipeline (iNAP) is an online analysis pipeline for generating and analyzing comprehensive ecological networks in microbiome studies. It is implemented in two sections, that is, network construction and network analysis, and integrates many open-access tools. Network construction contains multiple feasible alternatives, including correlation-based approaches (Pearson's correlation and Spearman's rank correlation along with random matrix theory, and sparse correlations for compositional data) and conditional dependence-based methods (extended local similarity analysis and sparse inverse covariance estimation for ecological association inference), while network analysis provides topological structures at different levels and the potential effects of environmental factors on network structures. Considering the full workflow, from microbiome data set to network result, iNAP contains the molecular ecological network analysis pipeline and interdomain ecological network analysis pipeline (IDENAP), which correspond to the intradomain and interdomain associations of microbial species at multiple taxonomic levels. Here, we describe the detailed workflow by taking IDENAP as an example and show the comprehensive steps to assist researchers to conduct the relevant analyses using their own data sets. Afterwards, some auxiliary tools facilitating the pipeline are introduced to effectively aid in the switch from local analysis to online operations. Therefore, iNAP, as an easy-to-use platform that provides multiple network-associated tools and approaches, can enable researchers to better understand the organization of microbial communities. iNAP is available at http://mem.rcees.ac.cn:8081 with free registration.

The Coexistence Relationship Between Plants and Soil Bacteria Based on Interdomain Ecological Network Analysis

The relationship between plants and their associated soil microbial communities plays a crucial role in maintaining ecosystem processes and function. However, identifying these complex relationships is challenging. In this study, we constructed an interdomain ecology network (IDEN) of plant–bacteria based on SparCC pairwise associations using synchronous aboveground plant surveys and belowground microbial 16S rRNA sequencing among four different natural forest types along the climate zones in China. The results found that a total of 48 plants were associated with soil bacteria among these four sites, and soil microbial group associations with specific plant species existed within the observed plant–bacteria coexistence network. Only 0.54% of operational taxonomy units (OTUs) was shared by the four sites, and the proportion of unique OTUs for each site ranged from 43.08 to 76.28%, which occupied a large proportion of soil bacterial community composition. The plant–bacteria network had a distinct modular structure (p < 0.001). The tree Acer tetramerum was identified as the network hubs in the warm temperate coniferous and broad-leaved mixed forests coexistence network and indicates that it may play a key role in stabilizing of the community structure of these forest ecosystems. Therefore, IDEN of plant–bacteria provides a novel perspective for exploring the relationships of interdomain species, and this study provides valuable insights into understanding coexistence between above-ground plants and below-ground microorganisms.

Pollution alters methanogenic and methanotrophic communities and increases dissolved methane in small ponds

Small ponds have become a hotspot of greenhouse gas emissions, but our understanding of methane (CH₄) cycling and its biological regulation in small polluted ponds remains limited. To assess how pollution affects CH₄ content, we investigated dissolved CH₄ concentrations, water and sediments properties, methanogenic and methanotrophic communities in two types of small polluted ponds. Compared with low pollution (LP) ponds, high pollution (HP) ponds showed significantly (P < 0.05) higher dissolved CH₄ in water. Sequencing of methyl coenzyme M reductase (mcrA) and particulate methane monooxygenase (pmoA) genes showed that HP led to significant (P < 0.05) shifts of CH₄-cycling microbial communities, with increased Shannon index of sediment methanogenic communities and water methanotrophic communities. There were also strong negative associations (P < 0.05) between dissolved CH₄ concentrations and interdomain methanogen-methanotroph network connectivity in water and sediments, respectively. The partial least squares path modeling indicated that dissolved oxygen, total organic carbon, ammonium nitrogen and nitrate nitrogen of water, and total nitrogen and total carbon of sediment, and CH₄-cycling microbes could regulate the CH₄ content. This study clarified the effects of environmental deterioration on CH₄ cycling in small ponds, highlighting the use of methanogen-methanotroph network connectivity to assess the CH₄ production.

Warming reshaped the microbial hierarchical interactions

Global warming may alter microbially mediated ecosystem functions through reshaping of microbial diversity and modified microbial interactions. Here, we examined the effects of 5-year experimental warming on different microbial hierarchical groups in a coastal nontidal soil ecosystem, including prokaryotes (i.e., bacteria and archaea), fungi, and Cercozoa, which is a widespread phylum of protists. Warming significantly altered the diversity and structure of prokaryotic and fungal communities in soil and additionally decreased the complexity of the prokaryotic network and fragmented the cercozoan network. By using the Inter-Domain Ecological Network approach, the cross-trophic interactions among prokaryotes, fungi, and Cercozoa were further investigated. Under warming, cercozoan-prokaryotic and fungal-prokaryotic bipartite networks were simplified, whereas the cercozoan-fungal network became slightly more complex. Despite simplification of the fungal-prokaryotic network, the strengthened synergistic interactions between saprotrophic fungi and certain prokaryotic groups, such as the Bacteroidetes, retained these phyla within the network under warming. In addition, the interactions within the fungal community were quite stable under warming conditions, which stabilized the interactions between fungi and prokaryotes or protists. Additionally, we found the microbial hierarchical interactions were affected by environmental stress (i.e., salinity and pH) and soil nutrients. Interestingly, the relevant microbial groups could respond to different soil properties under ambient conditions, whereas under warming these two groups tended to respond to similar soil properties, suggesting network hub species responded to certain environmental changes related to warming, and then transferred this response to their partners through trophic interactions. Finally, warming strengthened the network modules' negative association with soil organic matters through some fungal hub species, which might trigger soil carbon loss in this ecosystem. Our study provides new insights into the response and feedback of microbial hierarchical interactions under warming scenario.

A review of root exudates and rhizosphere microbiome for crop production

Increasing crop yields and ensuring food security is a major global challenge. In order to increase crop production, chemical fertilizers and pesticides are excessively used. However, the significance of root exudates is understudied. Beneficial interactions between plant and rhizosphere microbiome are critical for plant fitness and health. In this review, we discuss the application and progress of current research methods and technologies in terms of root exudates and rhizosphere microbiome. We summarize how root exudates promote plant access to nitrogen, phosphorus, and iron, and how root exudates strengthen plant immunity to cope with biotic stress by regulating the rhizosphere microbiome, and thereby reducing dependence on fertilizers and pesticides. Optimizing these interactions to increase plant nutrient uptake and resistance to biotic stresses offers one of the few untapped opportunities to confront sustainability issues in food security. To overcome the limitations of current research, combination of multi-omics, imaging technology together with synthetic communities has the potential to uncover the interaction mechanisms and to fill the knowledge gap for their applications in agriculture to achieve sustainable development.

Fungi-Bacteria Associations in Wilt Diseased Rhizosphere and Endosphere by Interdomain Ecological Network Analysis

In the plant rhizosphere and endosphere, some fungal and bacterial species regularly co-exist, however, our knowledge about their co-existence patterns is quite limited, especially during invasion by bacterial wilt pathogens. In this study, the fungal communities from soil to endophytic compartments were surveyed during an outbreak of tobacco wilt disease caused by Ralstonia solanacearum. It was found that the stem endophytic fungal community was significantly altered by pathogen invasion in terms of community diversity, structure, and composition. The associations among fungal species in the rhizosphere and endosphere infected by R. solanacearum showed more complex network structures than those of healthy plants. By integrating the bacterial dataset, associations between fungi and bacteria were inferred by Inter-Domain Ecological Network (IDEN) approach. It also revealed that infected samples, including both the rhizosphere and endosphere, had more complex interdomain networks than the corresponding healthy samples. Additionally, the bacterial wilt pathogenic Ralstonia members were identified as the keystone genus within the IDENs of both root and stem endophytic compartments. Ralstonia members was negatively correlated with the fungal genera Phoma, Gibberella, and Alternaria in infected roots, as well as Phoma, Gibberella, and Diaporthe in infected stems. This suggested that those endophytic fungi may play an important role in resisting the invasion of R. solanacearum.

High-Throughput Single-Cell Technology Reveals the Contribution of Horizontal Gene Transfer to Typical Antibiotic Resistance Gene Dissemination in Wastewater Treatment Plants

The spread of antibiotic resistance genes (ARGs) has gained much attention worldwide, while the contribution of vertical gene transfer (VGT) and horizontal gene transfer (HGT) is still elusive. Here, we improved an emerging high-throughput single-cell-based technology, emulsion, paired isolation, and concatenation polymerase chain reaction (epicPCR), by lengthening the sequence of ARG in the fused ARG-16S rRNA fragments to cover the variance of both ARG and its hosts. The improved epicPCR was applied to track the hosts of a widely detected ARG, sul1 gene, in five urban wastewater treatment plants (UWTPs) during two seasons. The sul1 host bacteria were highly diverse and mostly classified as Proteobacteria and Bacteroidetes. Clear seasonal divergence of α-diversity and interaction networks were present in the host community. The consensus phylogenetic trees of the sul1 gene and their host demonstrated incorrespondence on the whole and regularity on abundant groups, suggesting the important role of both HGT and VGT, respectively. The relative importance of these two ways was further measured; HGT (54%) generally played an equal or even more important role as VGT (46%) in UWTPs. The application of the improved epicPCR technology provides a feasible approach to quantify the relative contributions of VGT and HGT in environmental dissemination of ARGs.

How Mycorrhizal Associations Influence Orchid Distribution and Population Dynamics

Orchid distribution and population dynamics are influenced by a variety of ecological factors and the formation of holobionts, which play key roles in colonization and ecological community construction. Seed germination, seedling establishment, reproduction, and survival of orchid species are strongly dependent on orchid mycorrhizal fungi (OMF), with mycorrhizal cheating increasingly observed in photosynthetic orchids. Therefore, changes in the composition and abundance of OMF can have profound effects on orchid distribution and fitness. Network analysis is an important tool for the study of interactions between plants, microbes, and the environment, because of the insights that it can provide into the interactions and coexistence patterns among species. Here, we provide a comprehensive overview, systematically describing the current research status of the effects of OMF on orchid distribution and dynamics, phylogenetic signals in orchid-OMF interactions, and OMF networks. We argue that orchid-OMF associations exhibit complementary and specific effects that are highly adapted to their environment. Such specificity of associations may affect the niche breadth of orchid species and act as a stabilizing force in plant-microbe coevolution. We postulate that network analysis is required to elucidate the functions of fungal partners beyond their effects on germination and growth. Such studies may lend insight into the microbial ecology of orchids and provide a scientific basis for the protection of orchids under natural conditions in an efficient and cost-effective manner.

Soil microbial communities influencing organic phosphorus mineralization in a coastal dune chronosequence in New Zealand

The Haast chronosequence in New Zealand is an ∼6500-year dune formation series, characterized by rapid podzol development, phosphorus (P) depletion and a decline in aboveground biomass. We examined bacterial and fungal community composition within mineral soil fractions using amplicon-based high-throughput sequencing (Illumina MiSeq). We targeted bacterial non-specific acid (class A, phoN/phoC) and alkaline (phoD) phosphomonoesterase genes and quantified specific genes and transcripts using real-time PCR. Soil bacterial diversity was greatest after 4000 years of ecosystem development and associated with an increased richness of phylotypes and a significant decline in previously dominant taxa (Firmicutes and Proteobacteria). Soil fungal communities transitioned from predominantly Basidiomycota to Ascomycota along the chronosequence and were most diverse in 290- to 392-year-old soils, coinciding with maximum tree basal area and organic P accumulation. The Bacteria:Fungi ratio decreased amid a competitive and interconnected soil community as determined by network analysis. Overall, soil microbial communities were associated with soil changes and declining P throughout pedogenesis and ecosystem succession. We identified an increased dependence on organic P mineralization, as found by the profiled acid phosphatase genes, soil acid phosphatase activity and function inference from predicted metagenomes (PICRUSt2).

Fungal-Bacterial Cooccurrence Patterns Differ between Arbuscular Mycorrhizal Fungi and Nonmycorrhizal Fungi across Soil Niches

Soil bacteria and fungi are known to form niche-specific communities that differ between actively growing and decaying roots. Yet almost nothing is known about the cross-kingdom interactions that frame these communities and the environmental filtering that defines these potentially friendly or competing neighbors. We explored the temporal and spatial patterns of soil fungal (mycorrhizal and nonmycorrhizal) and bacterial cooccurrence near roots of wild oat grass, Avena fatua, growing in its naturalized soil in a greenhouse experiment. Amplicon sequences of the fungal internal transcribed spacer (ITS) and bacterial 16S rRNA genes from rhizosphere and bulk soils collected at multiple plant growth stages were used to construct covariation-based networks as a step toward identifying fungal-bacterial associations. Corresponding stable-isotope-enabled metagenome-assembled genomes (MAGs) of bacteria identified in cooccurrence networks were used to inform potential mechanisms underlying the observed links. Bacterial-fungal networks were significantly different in rhizosphere versus bulk soils and between arbuscular mycorrhizal fungi (AMF) and nonmycorrhizal fungi. Over 12 weeks of plant growth, nonmycorrhizal fungi formed increasingly complex networks with bacteria in rhizosphere soils, while AMF more frequently formed networks with bacteria in bulk soils. Analysis of network-associated bacterial MAGs suggests that some of the fungal-bacterial links that we identified are potential indicators of bacterial breakdown and consumption of fungal biomass, while others intimate shared ecological niches.IMPORTANCE Soils near living and decomposing roots form distinct niches that promote microorganisms with distinctive environmental preferences and interactions. Yet few studies have assessed the community-level cooccurrence of bacteria and fungi in these soil niches as plant roots grow and senesce. With plant growth, we observed increasingly complex cooccurrence networks between nonmycorrhizal fungi and bacteria in the rhizosphere, while mycorrhizal fungal (AMF) and bacterial cooccurrence was more pronounced in soil further from roots, in the presence of decaying root litter. This rarely documented phenomenon suggests niche sharing of nonmycorrhizal fungi and bacteria, versus niche partitioning between AMF and bacteria; both patterns are likely driven by C substrate availability and quality. Although the implications of species cooccurrence are fiercely debated, MAGs matching the bacterial nodes in our networks possess the functional potential to interact with the fungi that they are linked to, suggesting an ecological significance of fungal-bacterial cooccurrence patterns.

Successional Change of the Fungal Microbiome Pine Seedling Roots Inoculated With Tricholoma matsutake

The pine mushroom (Tricholoma matsutake; Agaricales, Tricholomataceae) is an ectomycorrhizal fungus that produces a commercially valuable, edible mushrooms. Attempts to artificially cultivate T. matsutake has so far been unsuccessful. One method used to induce T. matsutake to produce fruiting bodies of in the wild is shiro (mycelial aggregations of T. matsutake) transplantation. In vitro ectomycorrhization of T. matsutake with seedlings of Pinus densiflora has been successful, but field trials showed limited production of fruiting bodies. Few studies have been done to test what happens after transplantation in the wild, whether T. matsutake persists on the pine seedling roots or gets replaced by other fungi. Here, we investigated the composition and the interaction of the root fungal microbiome of P. densiflora seedlings inoculated with T. matsutake over a 3 year period after field transplantation, using high-throughput sequencing. We found a decline of T. matsutake colonization on pine roots and succession of mycorrhizal fungi as P. densiflora seedlings grew. Early on, roots were colonized by fast-growing, saprotrophic Ascomycota, then later replaced by early stage ectomycorrhiza such as Wilcoxina. At the end, more competitive Suillus species dominated the host roots. Most of the major OTUs had negative or neutral correlation with T. matsutake, but several saprotrophic/plant pathogenic/mycoparasitic species in genera Fusarium, Oidiodendron, and Trichoderma had positive correlation with T. matsutake. Four keystone species were identified during succession; two species (Fusarium oxysporum, and F. trincintum) had a positive correlation with T. matsutake, while the other two had a negative correlation (Suillus granulatus, Cylindrocarpon pauciseptatum). These findings have important implications for further studies on the artificial cultivation of T. matsutake.

Host phylogeny, diet, and habitat differentiate the gut microbiomes of Darwin's finches on Santa Cruz Island

Darwin's finches are an iconic example of an adaptive radiation with well-characterized evolutionary history, dietary preferences, and biogeography, offering an unparalleled opportunity to disentangle effects of evolutionary history on host microbiome from other factors like diet and habitat. Here, we characterize the gut microbiome in Darwin's finches, comparing nine species that occupy diverse ecological niches on Santa Cruz island. The finch phylogeny showed moderate congruence with the microbiome, which was comprised mostly of the bacterial phyla Firmicutes, Actinobacteria, and Proteobacteria. Diet, as measured with stable isotope values and foraging observations, also correlated with microbiome differentiation. Additionally, each gut microbial community could easily be classified by the habitat of origin independent of host species. Altogether, these findings are consistent with a model of microbiome assembly in which environmental filtering via diet and habitat are primary determinants of the bacterial taxa present with lesser influence from the evolutionary history between finch species.