Identification of keystone taxa in rhizosphere microbial communities using different methods and their effects on compounds of the host Cinnamomum migao

Exploring keystone taxa affecting microbial community stability and host function is crucial for understanding ecosystem functions. However, identifying keystone taxa from humongous microbial communities remains challenging. We collected 344 rhizosphere and bulk soil samples from the endangered plant C. migao for 2 years consecutively. Used high-throughput sequencing 16S rDNA and ITS to obtain the composition of bacterial and fungal communities. We explored keystone taxa and the applicability and limitations of five methods (SPEC-OCCU, Zi-Pi, Subnetwork, Betweenness, and Module), as well as the impact of microbial community domain, time series, and rhizosphere boundary on the identification of keystone taxa in the communities. Our results showed that the five methods, identified abundant keystone taxa in rhizosphere and bulk soil microbial communities. However, the keystone taxa shared by the rhizosphere and bulk soil microbial communities over time decreased rapidly decrease in the five methods. Among five methods on the identification of keystone taxa in the rhizosphere community, Module identified 113 taxa, SPEC-OCCU identified 17 taxa, Betweenness identified 3 taxa, Subnetwork identified 3 taxa, and Zi-Pi identified 4 taxa. The keystone taxa are mainly conditionally rare taxa, and their ecological functions include chemoheterotrophy, aerobic chemoheterotrophy, nitrate reduction, and anaerobic photoautotrophy. The results of the random forest model and structural equation model predict that keystone taxa Mortierella and Ellin6513 may have an effects on the accumulation of 1, 4, 7, - Cycloundecatriene, 1, 5, 9, 9-tetramethyl-, Z, Z, Z-, beta-copaene, bicyclogermacrene, 1,8-Cineole in C. migao fruits, but their effects still need further evidence. Our study evidence an unstable microbial community in the bulk soil, and the definition of microbial boundary and ecologically functional affected the identification of keystone taxa in the community. Subnetwork and Module are more in line with the definition of keystone taxa in microbial ecosystems in terms of maintaining community stability and hosting function.

Improved bacterial composition and co-occurrence patterns of rhizosphere increased nutrient uptake and grain yield through cultivars mixtures in maize

Crop diversification contributes to agricultural productivity and resources efficient utilization. However, whether cultivar mixtures in maize affects soil bacterial community, nutrient uptake, plant growth and yield remains unknown. A two-year lysimetric experiment was conducted using two maize cultivars (LY16 and JS501) with different root system architectures planted in monoculture or in mixture under normal fertilization (NF), reduced fertilization (RF) or no addition of fertilizer (CK) and was assessed at the silking stages. Cultivar mixtures and monoculture of LY16 had higher shoot biomass, nutrient uptake and total root length at silking stage, and grain yield than monoculture of JS501 under NF and RF conditions. Under NF and RF conditions, cultivar mixtures and monoculture of LY16 led to an increase in bacterial diversity, significant changes in community structure, and a high abundance of Bacteroidia and biomarkers of Chitinophagaceae and Saprospiraceae (Bacteroidia). Cultivar mixtures showed specific responses from modules of the rhizosphere bacterial community co-occurrence network, and the relative abundance of keystone taxa of cultivar mixtures was higher than that of monoculture of JS501. The keystone taxa had a broad and significant positive correlation with plant nutrient accumulation and grain yield. Cultivar mixtures showed similar assembly processes of Bacteroidia with monoculture of LY16, and the increased abundance of Chitinophagaceae may lead to a healthy rhizosphere bacterial community. Overall, our findings indicate that cultivar mixtures significantly affects the assembly and composition of the rhizosphere bacterial community, and thus benefits plant nutrient acquisition and plant growth. These findings could deepen our understanding of the facilitating effect of rhizosphere functional microbial community (e.g. plant nutrition uptake or immunity)of cultivar mixtures.

NosZ I carrying microorganisms determine N2O emissions from the subtropical paddy field under elevated CO2 and strongly CO2-responsive cultivar

Elevated CO2 (eCO2) decreases N2O emissions from subtropical paddy fields, but the underlying mechanisms remain to be investigated. Herein, the response of key microbial nitrogen cycling genes to eCO2 (ambient air +200 μmol CO2 mol-1) in four rice cultivars, including two weakly CO2-responsive (W27, H5) and two strongly CO2-responsive cultivars (Y1540, L1988), was investigated. Except for nosZ I, eCO2 did not significantly alter the abundance of the other genes. NosZ I was a crucial factor governing N2O emissions, especially under eCO2 and a strongly responsive cultivar. eCO2 affected the nosZ I gene abundance (p < 0.05), for instance, the nosZ I gene abundance of cultivar W27 increased from 1.53 × 107 to 2.86 × 107 copies g-1 dw soil (p < 0.05). In the nosZ I microbial community, the known taxa were mainly Pseudomonadota (phylum) (19.74-31.72 %) and Alphaproteobacteria (class) (0.56-13.12 %). In the nosZ I community assembly process, eCO2 enhanced the role of stochasticity, increasing from 35 % to 85 % (p < 0.05), thereby inducing diffusion limitations of weakly responsive cultivars to dominate (67 %). Taken together, the increase in nosZ I gene abundance is a potential reason for the alleviation of N2O emissions from subtropical paddy fields under eCO2.

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.

Arbuscular mycorrhizal fungi and Streptomyces: brothers in arms to shape the structure and function of the hyphosphere microbiome in the early stage of interaction

BACKGROUND

Fungi and bacteria coexist in a wide variety of environments, and their interactions are now recognized as the norm in most agroecosystems. These microbial communities harbor keystone taxa, which facilitate connectivity between fungal and bacterial communities, influencing their composition and functions. The roots of most plants are associated with arbuscular mycorrhizal (AM) fungi, which develop dense networks of hyphae in the soil. The surface of these hyphae (called the hyphosphere) is the region where multiple interactions with microbial communities can occur, e.g., exchanging or responding to each other's metabolites. However, the presence and importance of keystone taxa in the AM fungal hyphosphere remain largely unknown.

RESULTS

Here, we used in vitro and pot cultivation systems of AM fungi to investigate whether certain keystone bacteria were able to shape the microbial communities growing in the hyphosphere and potentially improved the fitness of the AM fungal host. Based on various AM fungi, soil leachates, and synthetic microbial communities, we found that under organic phosphorus (P) conditions, AM fungi could selectively recruit bacteria that enhanced their P nutrition and competed with less P-mobilizing bacteria. Specifically, we observed a privileged interaction between the isolate Streptomyces sp. D1 and AM fungi of the genus Rhizophagus, where (1) the carbon compounds exuded by the fungus were acquired by the bacterium which could mineralize organic P and (2) the in vitro culturable bacterial community residing on the surface of hyphae was in part regulated by Streptomyces sp. D1, primarily by inhibiting the bacteria with weak P-mineralizing ability, thereby enhancing AM fungi to acquire P.

CONCLUSIONS

This work highlights the multi-functionality of the keystone bacteria Streptomyces sp. D1 in fungal-bacteria and bacterial-bacterial interactions at the hyphal surface of AM fungi. Video Abstract.

Fungal, but not bacterial, diversity and network complexity promote network stability during roadside slope restoration

To restore degraded roadside ecosystems, conventional methods such as revegetation and soil amendment are frequently employed. However, our understanding of the long-term effects of these restoration approaches on soil microbial diversity and network complexity across different vegetation types remains poor, which contributes to poor restoration outcomes. In this study, we explored the effects of roadside slope restoration on microbial communities across different vegetation types at varying stages of restoration. We found that restoration time had a more pronounced impact on microbial diversity than specific vegetation type. As restoration progressed, microbial network complexity and fungal diversity increased, but bacterial diversity declined, suggesting that keystone taxa may contribute to network complexity. Interestingly, bacterial network complexity increased concomitant with decreasing network modularity and robustness, which may compromise system stability. Distinct vegetation types were associated with restoration-sensitive microbial communities at different restoration stages. Leguminouse and nitrogen-fixing plants, such as Albiziak alkora, Ginkgo biloba, Rhus chinensis, Rhapis excels, and Rubia cordifolia exhibited such associations after five years of restoration. These keystone taxa included Proteobacteria, Actinobacteria, Chloroflexi, Gemmatimonadota, and Myxococcota. We also found that bacterial alpha diversity was significantly correlated with restoration time, soil pH, moisture, available phosphate, nitrate nitrogen, and plant height, while fungal diversity was primarily shaped by restoration time. Together, our findings suggest that soil properties, environmental factors, vegetation type, and dominant species can be manipulated to guide the trajectory of ecological recovery by regulating the abundance of certain microbial taxa.

Keystone taxa enhance the stability of soil bacterial communities and multifunctionality under steelworks disturbance

Continuous discharge of wastewater, emissions, and solid wastes from steelworks poses environmental risks to ecosystems. However, the role of keystone taxa in maintaining multifunctional stability during environmental disturbances remains poorly understood. To address this, we investigated the community diversity, assembly mechanisms, and soil multifunctionality of soils collected from within the steelworks (I), within 2.5 km radius from the steelworks (E), and from an undisturbed area (CK) in Jiangsu Province, China, via 16 S rRNA sequencing. Significant differences were found in the Chao1 and the richness indexes of the total taxa (p < 0.05), while the diversity of keystone taxa was not significant at each site (p > 0.05). The deterministic processes for total taxa were 42.9%, 61.9% and 47.7% in CK, E, and I, respectively. Steelworks stress increased the deterministicity of keystone taxa from 52.3% in CK to 61.9% in E and I soils. The average multifunctionality indices were 0.518, 0.506 and 0.513 for CK, E and I, respectively. Although the soil multifunctionality was positive correlated with α diversity of both the total and keystone taxa, the average degree of keystone taxa in functional network increased significantly (79.96 and 65.58, respectively), while the average degree of total taxa decreased (44.59 and 51.25, respectively) in the E and I. This suggests keystone taxa contribute to promoting the stability of ecosystems. With increasing disturbance, keystone taxa shift their function from basic metabolism (ribosome biogenesis) to detoxification (xenobiotics biodegradation, metabolism, and benzoate degradation). Here we show that keystone taxa are the most important factor in maintaining stable microbial communities and functions, providing new insights for mitigating pollution stress and soil health protection.

How redox gradient potentially influences nitrate reduction coupled with sulfur cycling: A new insight into nitrogen cycling in the hyporheic zone of effluent-dominated rivers

The coupled N and S cycling in variable redox gradients in the hyporheic zone (HZ) of the rivers receiving effluents from wastewater treatment plants is unclear. Using two representative effluent-dominated rivers as model systems, a metagenome approach was employed to explore the spatiotemporal redox zonation of the HZ and the N/S cycling processes within the system. The results manifested that nitrate reduction represented the fundamental nitrogen pathway in the HZ. Interestingly, DNRA coupled with sulfur reduction, and denitrification coupled with sulfur oxidation were respectively abundant in the oxic and anoxic zone. Lower nitrate concentration (0-2.72 mg-N/L) and more abundant genes involved in denitrification (napB, NarGHI) and sulfur oxidation (sseA, glpE) were detected in the anoxic zone. Contrarily, the nitrate concentration (0.07-4.87 mg-N/L) and the abundance of genes involved in sulfur reduction (ttrB, sudA) and DNRA (nirBD) were observed more abundant in the oxic zone. Therefore, the results verified the oxygen-limited condition did not suppress but rather facilitated the denitrification process in the presence of active S cycling. The high relative abundances of nosZ gene encoding sequence (3-5 % relative to all nitrogen-cycling processes) in both the effluent-discharging area and downstream area highly confirmed that HZ was capable of alleviating the N2O emission in the region. The functional keystone taxa were revealed through co-occurrence network analysis. The structural equation model shows that the genes of N/S cycling were positively impacted by functional keystone taxa, especially the N cycling genes. Functional keystone taxa were proven driven by the redox gradient, demonstrating their positive roles in mediating N/S cycling processes. The promoting effect on nitrate reduction coupled with sulfur cycling was clarified when redox conditions oscillated, providing a new perspective on mitigating nitrogen pollution and greenhouse gas emissions in effluent-receiving rivers.

Organic farming promotes the abundance of fungi keystone taxa in bacteria-fungi interkingdom networks

Soil bacteria-fungi interactions are essential in the biogeochemical cycles of several nutrients, making these microbes major players in agroecosystems. While the impact of the farming system on microbial community composition has been extensively reported in the literature, whether sustainable farming approaches can promote associations between bacteria and fungi is still unclear. To study this, we employed 16S, ITS, and 18S DNA sequencing to uncover how microbial interactions were affected by conventional and organic farming systems on maize crops. The Bray-Curtis index revealed that bacterial, fungal, and arbuscular mycorrhizal fungi communities were significantly different between the two farming systems. Several taxa known to thrive in healthy soils, such as Nitrosophaerales, Orbiliales, and Glomus were more abundant in the organic farming system. Constrained ordination revealed that the organic farming system microbial community was significantly correlated with the β-glucosidase activity, whereas the conventional farming system microbial community significantly correlated with soil pH. Both conventional and organic co-occurrence interkingdom networks exhibited a parallel node count, however, the former had a higher number of edges, thus being denser than the latter. Despite the similar amount of fungal nodes in the co-occurrence networks, the organic farming system co-occurrence network exhibited more than 3-fold the proportion of fungal taxa as keystone nodes than the conventional co-occurrence network. The genera Bionectria, Cercophora, Geastrum, Penicillium, Preussia, Metarhizium, Myceliophthora, and Rhizophlyctis were among the fungal keystone nodes of the organic farming system network. Altogether, our results uncover that beyond differences in microbial community composition between the two farming systems, fungal keystone nodes are far more relevant in the organic farming system, thus suggesting that bacteria-fungi interactions are more frequent in organic farming systems, promoting a more functional microbial community.

Deciphering solute transport, microbiota assembly patterns and metabolic functions in the hyporheic zone of an effluent-dominated river

We lack a clear understanding of how anthropogenic pressures, exemplified by effluent discharge from wastewater treatment plants, destabilize microbial communities in the hyporheic zone (HZ) of receiving rivers. In this study, the spatiotemporal characteristics of hydrological parameters, and the physicochemical properties of surface and subsurface water in a representative effluent-dominated river were monitored. Sequencing of 16S rRNA amplicons and metagenomes revealed the microbial community structure in the HZ of both effluent discharge area and downstream region. The keystone taxa (taxa vital in determining the composition of each microbial cluster) and the keystone functions they controlled were subsequently identified. Effluent discharge amplified the depth of the oxic/suboxic zone and the hyporheic exchange fluxes in the effluent discharge area, which was 50-120% and 40-300% higher than in the downstream region, respectively. Microbial community structure pattern analysis demonstrated an enhancement in the rate of dispersal, an increase in microbial diversity, and an improved community network complexity in the effluent discharge area. By contrast, the number of keystone taxa in the effluent discharge area was 50-70% lower than that of the downstream region, resulting in reduced community network stability and functionality. The keystone taxa controlling metabolic functions in the networks categorized to effluent discharge area were comprised of more genera related to nitrogen and sulfur cycling, e.g., Dechloromonas, Desulfobacter, Flavobacterium, Nitrosomonas, etc., highlighting a research need in monitoring species associated with nutrient element cycling in the HZ of receiving waterbodies. The results showed that the keystone taxa could contribute positively to network stability, which was negatively correlated to hyporheic exchange fluxes and redox gradients. This study provides valuable insights that will improve our understanding of how river ecosystems respond to changes in anthropogenic pressures.

Unraveling assemblage of microbial community dwelling in Dabaoshan As/Pb/Zn mine-impacted area: A typical mountain mining area of South China

Microbial community assemblage includes microorganisms from the three domains including Bacteria, Archaea, and Eukarya (Fungi), which play a crucial role in geochemical cycles of metal(loid)s in mine tailings. Mine tailings harbor vast proportions of metal(loid)s, representing a unique source of co-contamination of metal(loid)s that threaten the environment. The elucidation of the assembly patterns of microbial communities in mining-impacted ecospheres has received little attention. To decipher the microbial community assembly processes, the microbial communities from the five sites of the Dabaoshan mine-impacted area were profiled by the MiSeq sequencing of 16S rRNA (Bacteria and Archaea) genes and internal transcribed spacers (Fungi). Results indicated that the coexistence of 31 bacterial, 10 fungal, and 3 archaeal phyla, were mainly dominated by Mucilaginibacter, Cladophialophora, and Candidatus Nitrosotalea, respectively. The distribution of microorganisms was controlled by deterministic processes. The combination of Cu, Pb, and Sb was the main factor explaining the structure of microbial communities. Functional predicting analysis of bacteria and archaea based on the phylogenetic investigation of communities by reconstruction of unobserved states analyses revealed that the metabolic pathways related to arsenite transporter, arsenate reductase, and FeS cluster were important for metal detoxification. Furthermore, the ecological guilds (pathogens, symbiotrophs, and saprotrophs) of fungal communities explained 44.5 % of functional prediction. In addition, metal-induced oxidative stress may be alleviated by antioxidant enzymes of fungi communities, such as catalase. Such information provides new insights into the microbial assembly rules in co-contaminated sites.

Cattle grazing management affects soil microbial diversity and community network complexity in the Northern Great Plains

Soil microbial communities play a vital role in the biogeochemical cycling and ecological functioning of grassland, but may be affected by common land uses such as cattle grazing. Changes in microbial diversity and network complexity can affect key ecosystem functions such as nutrient cycling. However, it is not well known how microbial diversity and network complexity respond to grazing in the Northern Great Plains. Consequently, it is important to understand whether variation in grazing management alters the diversity and complexity of grassland microbial communities. We compared the effect of intensive adaptive multi-paddock (AMP) grazing and conventional grazing practices on soil microbial communities using 16S/ITS amplicon sequencing. Samples were collected from grasslands in 13 AMP ranches and 13 neighboring, conventional ranches located across the Canadian prairies. We found that AMP grazing increased fungal diversity and evenness, and led to more complex microbial associations. Acidobacteria, Actinobacteria, Gemmatimonadetes, and Bacteroidetes were keystone taxa associated with AMP grazing, while Actinobacteria, Acidobacteria, Proteobacteria, and Armatimonadetes were keystone taxa under conventional grazing. Besides overall grazing treatment effects, specific grazing metrics like cattle stocking rate and rest-to-grazing ratio affected microbial richness and diversity. Bacterial and fungal richness increased with elevated stocking rate, and fungal richness and diversity increased directly with the rest-to-grazing ratio. These results suggest that AMP grazing may improve ecosystem by enhancing fungal diversity and increasing microbial network complexity and connectivity.

Bacteroides and related species: The keystone taxa of the human gut microbiota

Microbial communities play a significant role in maintaining ecosystems in a healthy homeostasis. Presently, in the human gastrointestinal tract, there are certain taxonomic groups of importance, though there is no single species that plays a keystone role. Bacteroides spp. are known to be major players in the maintenance of eubiosis in the human gastrointestinal tract. Here we review the critical role that Bacteroides play in the human gut, their potential pathogenic role outside of the gut, and their various methods of adapting to the environment, with a focus on data for B. fragilis and B. thetaiotaomicron. Bacteroides are anaerobic non-sporing Gram negative organisms that are also resistant to bile acids, generally thriving in the gut and having a beneficial relationship with the host. While they are generally commensal organisms, some Bacteroides spp. can be opportunistic pathogens in scenarios of GI disease, trauma, cancer, or GI surgery, and cause infection, most commonly intra-abdominal infection. B. fragilis can develop antimicrobial resistance through multiple mechanisms in large part due to its plasticity and fluid genome. Bacteroidota (formerly, Bacteroidetes) have a very broad metabolic potential in the GI microbiota and can rapidly adapt their carbohydrate metabolism to the available nutrients. Gastrointestinal Bacteroidota species produce short-chain fatty acids such as succinate, acetate, butyrate, and occasionally propionate, as the major end-products, which have wide-ranging and many beneficial influences on the host. Bacteroidota, via bile acid metabolism, also play a role in in colonization-resistance of other organisms, including Clostridioides difficile, and maintenance of gut integrity.

Dominant role of rare bacterial taxa rather than abundant taxa in driving the tailing primary succession

Primary ecological succession is imperative for tailing vegetation, driven notably by microbes that enhance tailing nutrient status. Yet, the roles of abundant and rare taxa in tailing primary succession remain underexplored. This study investigates these subcommunities across three succession stages (i.e., original tailing, biological crusts, grasslands). Throughout primary succession, alpha diversity and functional gene abundances of the rare taxa (RT) group consistently rise from bare tailings to grasslands. Conversely, the abundant taxa (AT) group displays an opposing trend. Intriguingly, employing co-occurrence networks, keystone taxa, mantel tests, similarity percentage analysis, and structural equation model, the study uncovers that RT wields a more pivotal role than AT in driving tailing primary succession. Community assembly analysis reveals stochastic control of AT and deterministic control of RT. Additionally, primary succession reinforces stochastic processes in AT, while RT's deterministic process remains unaffected. By unveiling these dynamics, the research enriches our understanding of primary ecological succession in tailings. Recognition of unique diversity patterns and community assembly mechanisms for rare and abundant subcommunities advances tailing ecosystem comprehension and informs ecological restoration strategies. This study thus contributes valuable insights to the complex interplay of microbial taxa during tailing primary succession.

Understanding how various forms of phosphorus stress affect microbiome functions and boost plant disease resistance: Insights from metagenomic analysis

The plant's response to phosphorus (P) starvation suppresses its immunity and regulates rhizosphere microbial colonization. However, the impact of various P forms on plant disease resistance and microbial composition remains underreported. This paper examines the soybean rhizosphere microbiome facing co-stress from Fusarium oxysporum and diverse P forms. Macrogenomic analysis evaluates whether P addition enhances plant disease resistance and rhizosphere microbial function, and if such effects relate to P forms. Results show that different P forms mitigate F. oxysporum-induced plant inhibition by promoting P turnover. P forms predominantly affect microbial composition, followed by soil and plant properties. In soybean, the phosphate transport strategy (ugpA/Q) was selected to maintain high P to enhance immunity in the KH2PO4 treatment, while organo-P mineralization (phnH/F/W/G) was selected for superphosphate treatment. The Frankiales, a P-turnover microorganism, copiotrophic microorganisms, and indicator bacteria of plant properties, initially increase after F. oxysporum inoculation and then decrease post P addition, regardless of P forms. Additionally, the rhizosphere microbial community's metabolic activities and compounds significantly aid soybean defense against F. oxysporum, with functional types depending on P forms. Therefore, these findings establish a novel approach to enhance host defense against soil-borne diseases through P nutrition regulation to mediate host-driven metabolic activities of microbial communities.

Key taxa and mobilome-mediated responses co-reshape the soil antibiotic resistome under dazomet fumigation stress

Agrochemicals are emergingly being implicated in the widespread dissemination of antibiotic resistance genes (ARGs) in agroecosystems. However, minimal research exists on the disturbance of fumigant on soil ARGs. Focusing on a typical fumigant dazomet in a simulated soil microcosm, we characterized the dazomet-triggered timely response and longstanding dynamic of ARGs at one-fold and two-fold field recommended doses using metagenome and quantitative PCR. Dazomet treatments reduced 13.17%-69.98% of absolute abundance of 16S rRNA gene and targeted ARGs, but, awfully, boosted diversity and relative abundance of ARGs up to 1.33-1.60 and 1.62-1.90 folds, respectively. Approximately 77.28% of changes in relative abundance of ARGs could be explained by bacterial community and mobile genetic elements (MGEs). Mechanistically, primary hosts of ARGs shifted from Proteobacteria (control) to Firmicutes and Actinobacteria (treatments) accompanied with corresponding changes in their abundance by combining community analysis, host tracking analysis and antibiotic resistant bacteria assay. Meanwhile, dazomet exposure significantly increased the incidence of MGEs and stimulated the conjugation of antibiotic-resistant plasmid. In addition, absolute abundance of targeted ARGs gradually recovered in the post-fumigation stage. Collectively, our results elucidate the dazomet-triggered emergence and spread of soil ARGs and highlight the importance of navigating toward rational use of fumigant in agricultural fields.

Anthropogenic original DOM is a critical factor affecting LNA bacterial community assembly

We investigated the geographical and environmental distance-decay relationships for both of the two bacteria in the Haihe River, Tianjin, China. HNA bacteria exhibited a stronger geographical variation-dependent pattern while LNA bacteria exhibited a stronger environmental variation-dependent pattern. Variance partition analysis (VPA), Mantel test, and partial mantel test validated the discrepant impacts of geographical distance and environmental factors on their two communities. The heterogeneous selection dominated community assembly of LNA bacteria demonstrates their greater sensitivity to environmental conditions. As the deterministic environmental factor, anthropogenic original dissolved organic matter (DOM) functions exclusively on LNA bacteria, and it is the critical factor leading to the discrepant biogeographical patterns of LNA and HNA bacteria. LNA bacteria interact with HNA bacteria and mediate the DOM driving total bacteria assembly. The LNA keystone taxa, Pseudomonas, Rheinheimera, Candidatus Aquiluna, and hgcl clade are capable to compete with HNA bacteria for anthropogenic original DOM, and are potential indicators of anthropogenic pollution. Our research reveals the non-negligible effect of the LNA bacteria in regulating the ecological response of total bacteria.

Planktonic/benthic Bathyarchaeota as a "gatekeeper" enhance archaeal nonrandom co-existence and deterministic assembling in the Yangtze River

Archaea, the third proposed domain of life, mediate carbon and nutrient cycling in global natural habitats. Compared with bacteria, our knowledge about archaeal ecological modes in large freshwater environments subject to varying natural and human factors is limited. By metabarcoding analysis of 303 samples, we provided the first integrate biogeography about archaeal compositions, co-existence networks, and assembling processes within a 6000 km continuum of the Yangtze River. Our study revealed that, among the major phyla, water samples owned a higher proportion of Thaumarchaeota (62.8%), while sediments had higher proportions of Euryarchaeota (33.4%) and Bathyarchaeota (18.8%). A decline of polarization in phylum abundance profile was observed from plateau/mountain/hill to basin/plain areas, which was attributed to the increase of nutrients and metals. Planktonic and benthic Bathyarchaeota tended to co-occur with both major (e.g., methanogens or Thermoplasmata) and minor (e.g., Asgard or DPANN) taxa in the non-random networks, harboring the highest richness and abundances of keystone species and contributing the most positively to edge number, node degree, and nearest neighbor degree. Furthermore, we noted significantly positive contributions of Bathyarchaeota abundance and network complexity to the dominance of deterministic process in archaeal assembly (water: 65.3%; sediments: 92.6%), since higher carbon metabolic versatility of Bathyarchaeota would benefit archaeal symbiotic relations. Stronger deterministic assembling was identified at the lower-reach plain, and higher concentrations of ammonium and aluminum separately functioning as nutrition and agglomerator were the main environmental drivers. We lastly found that the Three Gorges Dam caused a simultaneous drop of benthic Bathyarchaeota abundance, network co-existence, and deterministic effects immediately downstream due to riverbed erosion as a local interference. These findings highlight that Bathyarchaeota are a "gatekeeper" to promote fluvial archaeal diversity, stability, and predictability under varying macroscopic and microscopic factors, expanding our knowledge about microbial ecology in freshwater biogeochemical cycling globally.

Variation of bacterial community assembly over developmental stages and midgut of Dermanyssus gallinae

Bacterial microbiota play an important role in the fitness of arthropods, but the bacterial microflora in the parasitic mite Dermanyssus gallinae is only partially explored; there are gaps in our understanding of the microbiota localization and in our knowledge of microbial community assembly. In this work, we have visualized, quantified the abundance, and determined the diversity of bacterial occupancy, not only across developmental stages of D. gallinae, but also in the midgut of micro-dissected female D. gallinae mites. We explored community assembly and the presence of keystone taxa, as well as predicted metabolic functions in the microbiome of the mite. The diversity of the microbiota and the complexity of co-occurrence networks decreased with the progression of the life cycle. However, several bacterial taxa were present in all samples examined, indicating a core symbiotic consortium of bacteria. The relatively higher bacterial abundance in adult females, specifically in their midguts, implicates a function linked to the biology of D. gallinae mites. If such an association proves to be important, the bacterial microflora qualifies itself as an acaricidal or vaccine target against this troublesome pest.

Pathogen resistance in soils associated with bacteriome network reconstruction through reductive soil disinfestation

Reductive soil disinfestation (RSD) is an effective bioremediation technique to restructure the soil microbial community and eliminate soilborne phytopathogens. Yet we still lack a comprehensive understanding of the keystone taxa involved and their roles in ecosystem functioning in degraded soils treated by RSD. In this study, the bacteriome network structure in RSD-treated soil and the subsequent cultivation process were explored. As a result, bacterial communities in RSD-treated soil developed more complex topologies and stable co-occurrence patterns. The richness and diversity of keystone taxa were higher in the RSD group (module hub: 0.57%; connector: 23.98%) than in the Control group (module hub: 0.16%; connector: 19.34%). The restoration of keystone taxa in RSD-treated soil was significantly (P < 0.01) correlated with soil pH, total organic carbon, and total nitrogen. Moreover, a strong negative correlation (r =  -0.712; P < 0.01) was found between keystone taxa richness and Fusarium abundance. Our results suggest that keystone taxa involved in the RSD network structure are capable of maintaining a flexible generalist mode of metabolism, namely with respect to nitrogen fixation, methylotrophy, and methanotrophy. Furthermore, distinct network modules composed by numerous anti-pathogen agents were formed in RSD-treated soil; i.e., the genera Hydrogenispora, Azotobacter, Sphingomonas, and Clostridium_8 under the soil treatment stage, and the genera Anaerolinea and Pseudarthrobacter under the plant cultivation stage. The study provides novel insights into the association between fungistasis and keystone or sensitive taxa in RSD-treated soil, with significant implications for comprehending the mechanisms of RSD. KEY POINTS: • RSD enhanced bacteriome network stability and restored keystone taxa. • Keystone taxa richness was negatively correlated with Fusarium abundance. • Distinct sensitive OTUs and modules were formed in RSD soil.

Effects of 2,6-di-tert-butyl-hydroxytotulene and mineral-lubricant base oils on microbial communities during lubricants biodegradation

2,6-Di-tert-butyl-hydroxytotulene (BHT) is an additive commonly used in the manufacturing of lubricants to improve their antioxidant properties. However, in this study, we found that BHT affects the biodegradation of bio-lubricants by influencing the microbial community during the degradation of bio-lubricants. Specifically, BHT was found to reduce bacterial richness in activated sludge, but it increased the relative abundance of Actinobacteria (from 21.24% to 40.89%), Rhodococcus (from 17.15% to 31.25%), Dietzia (from 0.069% to 6.49%), and Aequorivita (from 0.90% to 1.85%). LEfSe analysis and co-occurrence network analysis suggested that Actinobacteria could be potential biomarkers and keystone taxa in microbial communities. Using the MetaCyc pathway database, the study found that BHT interfered with cellular biosynthetic processes. Additionally, the study also showed that mineral-lubricant base oils, which are difficult to degrade, significantly altered the diversity and composition of the microbiome. Overall, the findings demonstrate that BHT and mineral-lubricant base oils can substantially alter bacterial richness, structure, and function, potentially contributing to the difficulty in degrading lubricants. These findings have implications for the development of more biodegradable lubricants and the management of industrial waste containing lubricants.

Organic amendments alter microbiota assembly to stimulate soil metabolism for improving soil quality in wheat-maize rotation system

Straw retention (SR) and organic fertilizer (OF) application contribute to improve soil quality, but it is unclear how the soil microbial assemblage under organic amendments mediate soil biochemical metabolism pathways to perform it. This study collected soil samples from wheat field under different application of fertilizer (chemical fertilizer, as control; SR, and OF) in North China Plain, and systematically investigated the interlinkages among microbe assemblages, metabolites, and physicochemical properties. Results showed that the soil organic carbon (SOC) and permanganate oxidizable organic carbon (LOC) in soil samples followed the trend as OF > SR > control, and the activity of C-acquiring enzymes presented significantly positive correlation with SOC and LOC. In organic amendments, bacteria and fungi community were respectively dominated by deterministic and stochastic processes, while OF exerted more selective pressure on soil microbe. Compared with SR, OF had greater potential to boost the microbial community robustness through increasing the natural connectivity and stimulating fungal taxa activities in inter-kingdom microbial networks. Altogether 67 soil metabolites were significantly affected by organic amendments, most of them belonged to benzenoids (Ben), lipids and lipid-like molecules (LL), and organic acids and derivatives (OA). These metabolites were mainly derived from lipid and amino acid metabolism pathways. A list of keystone genera such as stachybotrys and phytohabitans were identified as important to soil metabolites, SOC, and C-acquiring enzyme activity. Structural equation modeling showed that soil quality properties were closely associated with LL, OA, and PP drove by microbial community assembly and keystone genera. Overall, these findings suggested that straw and organic fertilizer might drive keystone genera dominated by determinism to mediate soil lipid and amino acid metabolism for improving soil quality, which provided new insights into understanding the microbial-mediated biological process in amending soil quality.

May microbial ecological baseline exist in continental groundwater?

BACKGROUND

Microbes constitute almost the entire biological community in subsurface groundwater and play an important role in ecological evolution and global biogeochemical cycles. Ecological baseline as a fundamental reference with less human interference has been investigated in surface ecosystems such as soils, rivers, and ocean, but the existence of groundwater microbial ecological baseline (GMEB) is still an open question so far.

RESULTS

Based on high-throughput sequencing information derived from national monitoring of 733 newly constructed wells, we find that bacterial communities in pristine groundwater exhibit a significant lateral diversity gradient and gradually approach the topsoil microbial latitudinal diversity gradient with decreasing burial depth of phreatic water. Among 74 phyla dominated by Proteobacteria in groundwater, Patescibacteria act as keystone taxa that harmonize microbes in shallower aquifers and accelerate decline in bacterial diversity with increasing well-depth. Decreasing habitat niche breadth with increasing well-depth suggests a general change in the relationship among key microbes from closer cooperation in shallow to stronger competition in deep groundwater. Unlike surface-water microbes, microbial communities in pristine groundwater are predominantly shaped by deterministic processes, potentially associated with nutrient sequestration under dark and anoxic environments in aquifers.

CONCLUSIONS

By unveiling the biogeographic patterns and mechanisms controlling the community assembly of microbes in pristine groundwater throughout China, we firstly confirm the existence of GMEB in shallower aquifers and propose Groundwater Microbial Community Index (GMCI) to evaluate anthropogenic impact, which highlights the importance of GMEB in groundwater water security and health diagnosis. Video Abstract.

Types of vegetables shape composition, diversity, and co-occurrence networks of soil bacteria and fungi in karst areas of southwest China

BACKGROUND

Microorganisms are of significant importance in soil. Yet their association with specific vegetable types remains poorly comprehended. This study investigates the composition of bacterial and fungal communities in soil by employing high-throughput sequencing of 16 S rRNA genes and ITS rRNA genes while considering the cultivation of diverse vegetable varieties.

RESULTS

The findings indicate that the presence of cultivated vegetables influenced the bacterial and fungal communities leading to discernible alterations when compared to uncultivated soil. In particular, the soil of leafy vegetables (such as cabbage and kale) exhibited higher bacterial α-diversity than melon and fruit vegetable (such as cucumber and tomato), while fungal α-diversity showed an inverse pattern. The prevailing bacterial phyla in both leafy vegetable and melon and fruit vegetable soils were Proteobacteria, Acidobacteriota, Actinobacteriota, and Chloroflexi. In leafy vegetable soil, dominant fungal phyla included Ascomycota, Olpidiomycota, Mortierellomycota, and Basidiomycota whereas in melon and fruit vegetable soil. Ascomycota, Mortierellomycota, Basidiomycota, and Rozellomycota held prominence. Notably, the relative abundance of Ascomycota was lower in leafy vegetable soil compared to melon and fruit vegetable soil. Moreover, leafy vegetable soil exhibited a more complex and stable co-occurrence network in comparison to melon and fruit vegetable soil.

CONCLUSION

The findings enhance our understanding of how cultivated soil bacteria and fungi respond to human disturbance, thereby providing a valuable theoretical basis for soil health in degraded karst areas of southwest China.

Characterization of two keystone taxa, sulfur-oxidizing, and nitrate-reducing bacteria, by tracking their role transitions in the benzo[a]pyrene degradative microbiome

BACKGROUND

Keystone taxa are drivers of microbiome structure and functioning, which may play critical roles in microbiome-level responses to recalcitrant pollution and are a key to bioremediation. However, the characterization and manipulation of such taxa is a major challenge due to the complexity of microbial communities and rapid turnover in both time and space. Here, microcosms were set up with benzo[a]-pyrene (BaP) and/or nitrate based on C-rich, S-rich, and N-limited mangrove sediments as reductive experimental models to trigger and track the turnover of keystone taxa to address this challenge.

RESULTS

Based on microbial co-occurrence network analysis, two keystone taxa, Sulfurovum and Sulfurimonas, were found to exhibit significant role transitions in different microcosms, where these two taxa played nonkeystone roles with neutral relationships in in situ mangrove sediments. However, Sulfurimonas transitioned to be keystone taxa in nitrate-replenished microcosms and formed a keystone guild with Thioalkalispira. Sulfurovum stood out in BaP-added microcosms and mutualized in a densely polycyclic aromatic hydrocarbon (PAH)-degrader-centric keystone guild with Novosphingobium and Robiginitalea, where 63.25% of added BaP was removed. Under the occurrence of nitrate and BaP, they simultaneously played roles as keystone taxa in their respective guilds but exhibited significant competition. Comparative genomics and metagenome-assembled genome (MAG) analysis was then performed to reveal the metabolic potential of those keystone taxa and to empirically deduce their functional role in keystone guilds. Sulfurimonas possesses a better sense system and motility, indicative of its aggressive role in nitrate acquisition and conversion; Sulfurovum exhibited a better ability for oxidation resistance and transporting nutrients and electrons. High-efficiency thermal asymmetric interlaced polymerase reaction (hiTAIL-PCR) and enhanced green fluorescent protein (eGFP)-labeling approaches were employed to capture and label the BaP key degrader to further experimentally verify the roles of keystone taxa Sulfurovum in the keystone guilds. Observations of the enhancement in reactive oxygen species (ROS) removal, cell growth, and degradation efficiency by co-culture of isolated keystone taxa strains experimentally demonstrated that Sulfurovum contributes to the BaP degradative microbiome against BaP toxicity.

CONCLUSIONS

Our findings suggest that the combined use of co-occurrence network analysis, comparative genomics, and co-culture of captured keystone taxa (3C-strategy) in microbial communities whose structure is strongly shaped by changing environmental factors can characterize keystone taxa roles in keystone guilds and may provide targets for manipulation to improve the function of the microbiome. Video Abstract.

Depth-dependent effects of tree species identity on soil microbial community characteristics and multifunctionality

Soil microbes play key roles that support forest ecosystem functioning, while their community characteristics are strongly determined by tree species identity. However, the majority studies primarily focus on soil microorganisms in the topsoil, resulting in limited understanding of the linkages between tree species identity and the microbial communities that inhabit deep soils. Here we investigated the diversity, structure, function, and co-occurrence networks of soil bacterial and fungal communities, as well as related soil physicochemical properties, to a depth of two meters in dryland forests dominated by either Pinus tabuliformis, a native coniferous species, Robinia pseudoacacia, an exotic broadleaf and nitrogen-fixing species, or both. Tree species identity had stronger effects on soil multifunctionality and microbial community structure in the deep layers (80-200 cm) than in the top layers (0-60 cm). In addition, fungal communities were more responsive to tree species identity, whereas bacteria were more sensitive to soil depth. Tree species identity strongly influenced microbial network stability and complexity, with higher quantities in R. pseudoacacia than the other plantations, by affecting microbial composition and their associations. The increased in microbial network complexity and the relative abundance of keystone taxa enhance the soil multifunctionality of microbial productivity, sugar and chitin degradation, and nutrient availability and cycling. Meanwhile, the relative abundance of keystone taxa was more representative of soil multifunctionality than microbial diversity. Our study highlights that tree species identity significantly influences soil microbial community characteristics and multifunctionality, especially in deep soils, which will help us understand soil nutrients processed in plantation forest ecosystem and provide a reference for tree species selection in ecological restoration.

Bacterial community dynamics and co-occurrence network patterns during different stages of biochar-driven composting

Bacterial communities were dynamically tracked at four stages of biochar-driven sheep manure pile composting, and the co-occurrence networks with keystone taxa were established. The succession of bacterial community obvious varied during the composting process, Proteobacteria predominant in initial stage (39%) then shifted into Firmicutes in thermophilic (41%) and mesophilic (27%) stages, finally the maturation stage dominant by Bacteroidota (26%). Visualizations of bacterial co-occurrence networks demonstrate more cooperative mutualism and complex interactions in the thermophilic and mesophilic phases. Noticeably, the 7.5 and 10% biochar amended composts shown highest connections (736 and 663 total links) and positive cooperation (97.37 and 97.13% positive link) as well as higher closeness centrality and betweenness centrality of keystone taxa. Overall, appropriate biochar addition alters bacterial community succession and strengthens connection between keystone taxa and other bacteria, with 7.5 and 10% biochar amended composts has intense mutualistic symbiosis among bacterial communities.

Local-Scale Damming Impact on the Planktonic Bacterial and Eukaryotic Assemblages in the upper Yangtze River

Dam construction and impoundment cause discontinuities in the natural biophysical gradients in rivers. These discontinuities may alter distinctive habitats and different microbial community assembly mechanisms upstream and downstream of dams, which reflect the potential impacts of damming on riverine aquatic ecosystems. In this study, we investigated the planktonic microbial assemblages of three large dams in the upper Yangtze River by using high-throughput sequencing. The results revealed that the alpha diversity indexes increased downstream of the dams. In addition, more eukaryotic ASVs solely occurred downstream of the dams, which indicated that a large proportion of eukaryotes appeared downstream of the dams. The nonmetric multidimensional scaling analysis indicated that there was no obvious geographic clustering of the planktonic microbial assemblages among the different locations or among the different dams. However, the dam barriers changed dam-related variables (maximum dam height and water level) and local environmental variables (water temperature, DOC, etc.) that could possibly affect the assembly of the planktonic microbial communities that are closest to the dams. A co-occurrence network analysis demonstrated that the keystone taxa of the planktonic bacteria and eukaryotes decreased downstream of the dams. In particular, the keystone taxa of the eukaryotes disappeared downstream of the dams. The robustness analysis indicated that the natural connectivity of the microbial networks decreased more rapidly upstream of the dams, and the downstream eukaryotic network was more stable. In conclusion, damming has a greater impact on planktonic eukaryotes than on bacteria in near-dam areas, and planktonic microbial assemblages were more susceptible to the environmental changes. Our study provides a better understanding of the ecological effects of river damming.

Organic Matter Decomposition in River Ecosystems: Microbial Interactions Influenced by Total Nitrogen and Temperature in River Water

Microbes contribute to the organic matter decomposition (OMD) in river ecosystems. This study considers two aspects of OMD in river ecosystems which have not been examined in scientific studies previously, and these are the microbial interactions in OMD and the influence of environmental factors on microbial interactions. Cotton strip (CS), as a substitute for organic matter, was introduced to Luanhe River Basin in China. The results of CS assay, microbial sequencing, and redundancy analysis (RDA) showed that CS selectively enriched bacterial and fungal groups related to cellulose decomposition, achieving cotton strip decomposition (CSD). Bacterial phylum Proteobacteria and fungal phyla Rozellomycota and Ascomycota were the dominant groups associated with CSD. Network analysis and Mantel test results indicated that bacteria and fungi on CS cooperatively formed an interaction network to achieve the CSD. In the network, modules 2 and 4 were significantly positively associated with CSD, which were considered as the key modules in this study. The key modules were mainly composed of phyla Proteobacteria and Ascomycota, indicating that microbes in key modules were the effective decomposers of CS. Although keystone taxa were not directly associated with CSD, they may regulate the genera in key modules to achieve the CSD, since some keystone taxa were linked with the microbial genera associated with CSD in the key modules. Total nitrogen (TN) and temperature in water were the dominant environmental factors positively influenced CSD. The key modules 2 and 4 were positively influenced by water temperature and TN in water, respectively, and two keystone taxa were positively associated with TN. This profoundly revealed that water temperature and TN influenced the OMD through acting on the keystone taxa and key modules in microbial interactions. The research findings help us to understand the microbial interactions influenced by environmental factors in OMD in river ecosystems.

Identifying the effects of cropping with different pear cultivars on microbial community composition and networks in orchard soils

The role of plant genotype in determining the assembly of soil microorganisms is widely accepted; however, the effects of cropping with different cultivars of perennial crop plants on the composition of soil microbial communities are not fully understood. In the current study, high-throughput amplicon sequencing and real-time PCR were used to investigate the major features of bacterial community composition, ecological networks, and soil physicochemical properties in three replicate pear orchards, each planted with monocultures of pear cultivars Hosui (HS) or Sucui (SC) of similar ages. A distinct difference in the composition of microbial communities was observed between soils of HS and SC orchards. A significantly greater relative abundance of Verrucomicrobia and Alphaproteobacteria whereas a significantly lower relative abundance of Betaproteobacteria were found in soils of HS cropped orchards than that in SC orchards. Sphingomonas sp., belonging to the Alphaproteobacteria, was recognized as a key species in the co-occurrence network of the microbial interactions. Moreover, redundancy analysis, the Mantel correlation test, and random forest analysis showed that soil pH was the dominant driver in determining microbial community composition in HS soils, whereas soil organic matter was the primary factor determining microbial community composition in SC soils. Altogether, we provide evidence that soils in HS orchards harbor unique microbial communities enriched with respect to microbial groups associated with nutrient cycling, whereas soils in SC orchards are dominated by a group of beneficial microbes exhibiting plant growth promotion. These findings have implications for science-based guidance for manipulation of the soil microbiome to achieve sustainable food production.

Seasonal variation characteristics of inhalable bacteria in bioaerosols and antibiotic resistance genes in Harbin

Bioaerosols have received extensive attention due to their impact on climate, ecological environment, and human health. This study aimed to reveal the driving factors that structure bacterial community composition and the transmission route of antibiotic resistance genes (ARGs) in PM_2.5. The results showed that the bacterial concentration in spring (8.76 × 10⁵ copies/m³) was significantly higher than that in summer (1.03 × 10⁵ copies/m³) and winter (4.74 × 10⁴ copies/m³). Low temperatures and air pollution in winter negatively affected bacterial concentrations. Keystone taxa were identified by network analysis. Although about 50 % of the keystone taxa had low relative abundances, the strong impact of complex interactions between keystone taxa and other taxa on bacterial community structure deserved attention. The bacterial community assembly was dominated by stochastic processes (79.3 %). Interactions between bacteria and environmental filtering together affected bacterial community composition. Vertical gene transfer played an important role in the transmission of airborne ARGs. Given the potential integration and expression of ARGs in recipients, the human exposure risk due to high concentrations of ARGs and mobile genetic elements cannot be ignored. This study highlights human exposure to inhalable bacterial pathogens and ARGs in urban areas.

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.

Soil microbial community structure dynamics shape the rhizosphere priming effect patterns in the paddy soil

Microbial community structure plays a crucial part in soil organic carbon (SOC) decomposition and variation of rhizosphere priming effects (RPEs) during plant growth. However, it is still uncertain how bacterial community structure regulates RPEs in soil and how RPE patterns respond to plant growth. Therefore, we conducted an experiment to examine the RPE response to plant growth and nitrogen (N) addition (0 (N0), 150 (N150), and 300 (N300) kg N ha^-1) using the ¹³C natural abundance method in a C₃ soil (paddy soil) - C₄ plant (maize, Zea mays L.) system; we then explored the underlying biotic mechanisms using 16S rRNA sequencing techniques. Networks were constructed to identify keystone taxa and to analyze the correlations between network functional modules of bacterial community and C decomposition. The results indicated that negative and positive RPEs occurred on Day 30 and Day 75 after maize planting, respectively. Bacterial community structure significantly changed and tended to shift from r-strategists toward K-strategists with changing labile C: N stoichiometry and soil pH during plant growth stages. The different network modules of bacterial community were aggregated in response to RPE pattern variation. Caulobacteraceae, Bacillus, and Chitinophagaceae were keystone taxa on Day 30, while Gemmatimonas, Candidatus Koribacter, and Xanthobacteraceae were keystone taxa on Day 75. Moreover, keystone taxa with different C utilization strategies were significantly different between the two growth stages and related closely to different RPE patterns. This study provides deeper insights into the network structure of bacterial communities corresponding to RPE patterns and emphasizes the significance of keystone taxa in RPE variation.

Dual regulatory effects of microplastics and heat waves on river microbial carbon metabolism

Rivers play a critical role in the global carbon cycle, but the processes can be affected by widespread microplastic (MP) pollution and the increasing frequency of heat waves (HWs) in a warming climate. However, little is known about the role of river microbes in regulating the carbon cycle under the combined action of MP pollution and HWs. Here, through seven-day MP exposure and three cycles of HW simulation experiments, we found that MPs inhibited the thermal adaptation of the microbial community, thus regulating carbon metabolism. The CO₂ release level increased, while the carbon degradation ability and the preference for stable carbon were inhibited. Metabonomic, 16 S rRNA and ITS gene analyses further revealed that the regulation of carbon metabolism was closely related to the microbial r-/K- strategy, community assembly and transformation of keystone taxa. The random forest model revealed that dissolved oxygen and ammonia-nitrogen were important variables influencing microbial carbon metabolism. The above findings regarding microbe-mediated carbon metabolism provide insights into the effect of climate-related HWs on the ecological risks of MPs.

Early warning of red tides using bacterial and eukaryotic communities in nearshore waters

Anthropogenic discharge activities have increased nutrient pollution in coastal areas, leading to algal blooms and microbial community changes. Particularly, microbial communities could easily be affected with variation in nutrient pollution, and thus offered a promising strategy to predict early red tides warning via microbial community-levels variation and their keystone taxa hysteretic responses to nutrient pollution. Herein high-throughput sequencing technology from 52 samples were used to explore the variation of microbial communities and find the significant tipping points with aggravating nutrient conditions in Xiaoping Island coastal area. Results indicated that bacterial and microeukaryote communities were generally spatial and seasonal heterogeneity and were influenced by the different nutrient conditions. Procrustes test results showed that the comprehensive index of organics polluting (OPI), total nitrogen (TN), inorganic nitrogen (DIN), and total phosphorus (TP) were significantly correlated with the composition of bacteria and microeukaryotes. A SEGMENTED analysis revealed that the threshold of TN, DIN, and NH₄-N for bacterial community were 0.23 ± 0.091 mg/L, 0.21 ± 0.084 mg/L, 0.09 ± 0.057 mg/L, respectively. Tipping points for TN, DIN, and NH₄-N agreed with the concentration during Ceratium tripos and Skeletonema costatum blooms. Co-occurrence network results found that Planktomarina, Acinetobacter, and Verrucomicrobiaceae were keystone and OPI-discriminatory taxa. The abundant changes of Planktomarina at station A1 were significantly correlated with the development of C. tripos blooms (r = 0.55, p < 0.05), and also significantly correlated with TN, DIN, and NO₃-N (r≥|0.55|, p < 0.05). The abundant changes of Acinetobacter and Verrucomicrobiaceae at station C1 were significantly correlated with the development of C. tripos blooms (r ≥ 0.77, p < 0.05), and also significantly correlated with PO₄-P (r ≥ 0.64, p < 0.05). The dynamic abundance of keystone taxa showed that the trend of rapid changes could be monitored 1.5 months before the occurrence of red tide. Therefore, this study provides an assessment method for early warning of red tide occurrence and factors that trigger red tide.

Similar assembly mechanisms but distinct co-occurrence patterns of free-living vs. particle-attached bacterial communities across different habitats and seasons in shallow, eutrophic Lake Taihu

Eutrophication due to nitrogen and phosphorus input is an increasing problem in lake ecosystems. Free-living (FL) and particle-attached (PA) bacterial communities play a primary role in mediating biogeochemical processes in these lakes and in responding to eutrophication. However, knowledge of factors governing function, assembly mechanisms, and co-occurrence patterns of these communities remain poorly understood and are key challenges in microbial ecology. To address this knowledge gap, we collected 96 samples from Lake Taihu across four seasons and investigated the bacterial community using 16S rRNA gene sequencing. Our results demonstrate that the α-diversity, β-diversity, community composition, and functional composition of FL and PA bacterial communities exhibited differing spatiotemporal dynamics. FL and PA bacterial communities displayed similar distance-decay relationships across seasons. Deterministic processes (i.e., environmental filtering and species interaction) were the primary factors shaping community assembly in both FL and PA bacteria. Similar environmental factors shaped bacterial community structure while different environmental factors drove bacterial functional composition. Habitat filtering influenced enrichment of bacteria within specific functional groups. Among them, the FL bacterial community appeared to play a critical role in methane-utilization, whereas the PA bacteria contributed more to biogeochemical cycling of carbon. FL and PA bacterial communities exhibited distinct co-occurrence pattern across different seasons. In the FL network, Methylotenera and Methylophilaceae were identified as keystone taxa, while Burkholderiaceae and the hgcI clade were keystone taxa in the PA network. The PA bacterial community appeared to possess greater stability in the face of environmental change than did FL counterparts. These results broaden our knowledge of the driving factors, co-occurrence patterns, and assembly processes in FL and PA bacterial communities in eutrophic ecosystems and provide improved insight into the underlying mechanisms responsible for these results.

Microbial co-occurrence networks driven by low-abundance microbial taxa during composting dominate lignocellulose degradation

Lignocellulose, which contains cellulose, hemicellulose and lignin, is one of the most important factors determining the rate and quality of compost decomposition, and the microbial community composition affects the rate of lignocellulose decomposition. Interactions between microbial taxa contribute significantly to ecosystem energy flow and material cycling. However, it is not clear how interactions between microbial taxa affect the degradation of lignocellulose during the composting process. For this reason we carried out aerobic co-composting experiments with maize straw and cattle manure to explore the contribution of microbial community diversity and the interaction between taxa to lignocellulosic degradation. The results showed that moisture and temperature had the greatest effect on microbial communities during composting and that lignocellulose degradation was dominated by microbial co-occurrence networks rather than microbial community diversity. Overall co-occurrence network and bacterial-fungal interactions explained 23.9-84.1 % of lignocellulosic degradation, whereas microbial diversity only accounted for 24.6-31.5 %. Interestingly, keystone taxa analysis of the microbial co-occurrence networks revealed that low-abundance taxa influenced microbial interactions driving lignocellulose degradation. Our results provide a new perspective for understanding lignocellulose degradation during composting, offering insights into important microbial interaction mechanisms for improving compost quality and efficiency.

Nanoscale zero-valent iron changes microbial co-occurrence pattern in pentachlorophenol-contaminated soil

Nanoscale zero-valent iron (nZVI) is a prominent nanomaterial for the remediation of organochlorine-contaminated soil and groundwater. However, a knowledge gap regarding the effects of the coexistence of nZVI and pollutants on soil microorganisms remains. Here, we studied the effects of nZVI on the microbial community structure, co-occurrence network, and keystone taxa in pentachlorophenol (PCP, a typical organochlorine pesticide) contaminated soils. The addition of nZVI (1000 mg/kg) had no obvious recovery effect on the microbial community structure of PCP-contaminated soil, but enhanced the connection and lowered the modularity of the microbial network. These changes were mainly present in the bacterial network rather than in the fungal or archaeal network. Moreover, the addition of nZVI increased the number of keystone taxa in the PCP-contaminated soil from 29 to 76. These keystone taxa are related to the degradation of organochlorine pollutants, carbon metabolism, and nitrogen metabolism and may thus be helpful in recovering soil ecological functions. These findings provide new insights into the interaction among nanomaterials, microorganisms, and pollutants.

Synergistic Impacts of Arsenic and Antimony Co-contamination on Diazotrophic Communities

Nitrogen (N) shortage poses a great challenge to the implementation of in situ bioremediation practices in mining-contaminated sites. Diazotrophs can fix atmospheric N₂ into a bioavailable form to plants and microorganisms inhabiting adverse habitats. Increasing numbers of studies mainly focused on the diazotrophic communities in the agroecosystems, while those communities in mining areas are still not well understood. This study compared the variations of diazotrophic communities in composition and interactions in the mining areas with different extents of arsenic (As) and antimony (Sb) contamination. As and Sb co-contamination increased alpha diversities and the abundance of nifH encoding the dinitrogenase reductase, while inhibited the diazotrophic interactions and substantially changed the composition of communities. Based on the multiple lines of evidence (e.g., the enrichment analysis of diazotrophs, microbe-microbe network, and random forest regression), six diazotrophs (e.g., Sinorhizobium, Dechloromonas, Trichormus, Herbaspirillum, Desmonostoc, and Klebsiella) were identified as keystone taxa. Environment-microbe network and random forest prediction demonstrated that these keystone taxa were highly correlated with the As and Sb contamination fractions. All these results imply that the above-mentioned diazotrophs may be resistant to metal(loid)s.

Characterizing sediment bacterial community and identifying the biological indicators in a seawater-freshwater transition zone during the wet and dry seasons

Seawater intrusion has a detrimental effect on agriculture, industry, and human health. One question of particular interest is how the microbial community responds to and reflects seawater intrusion with seasonal variation. The current study explored the seasonal changes in bacterial community composition and interaction in the vicinity of Pearl River Estuary in dry season (January) and wet season (September). Results indicated that the salinity of sediment samples obtained in dry season was higher than that in wet season. The salt stress induced a declined alpha diversity but resulted in a loosely connected and unstable biotic interaction network in the bacterial communities. Random forest prediction and redundancy analysis of bacterial community indicated that salinity substantially affected the bacterial communities. Multiple lines of evidence, including the enrichment of bacterial taxa in the high-salinity location, microbe-microbe interactions, environment-microbe interactions, and machine learning approach, demonstrated that the families Moraxellaceae and Planococcaceae were the keystone taxa and were resistant to salt stress, which suggested that both of them can be used as potential biological indicators of monitoring and controlling seawater intrusion in coastal zone areas.

Microbial diversity drives pyrene dissipation in soil

Soil microbial diversity is an essential driver of multiple ecosystem functions and services. However, the role and mechanisms of microbial diversity in the dissipation of persistent organic pollutants in soil are largely unexplored. Here, a gradient of soil microbial diversity was constructed artificially by a dilution-to-extinction approach to assess the role of soil microbial diversity in the dissipation of pyrene, a high molecular weight polycyclic aromatic hydrocarbon (PAH), in a 42-day microcosm experiment. The results showed that pyrene dissipation (98.1%) and the abundances of pyrene degradation genes (the pyrene dioxygenase gene nidA and the gram-positive PAH-ring hydroxylating dioxygenase gene PAH-RHDα GP) were highest in soils with high microbial diversity. Random-forest machine learning was combined with linear regression analysis to identify a range of keystone taxa (order level) associated with pyrene dissipation, including Sphingobacteriales, Vampirovibrionales, Blastocatellales, Myxococcales, Micrococcales and Rhodobacterales. The diversity of these keystone taxa was significantly and positively correlated with the abundance of pyrene degradation genes and the removal rate of pyrene. According to (partial) Mantel tests, keystone taxa diversity was the dominant factor determining pyrene dissipation compared with total microbial diversity. Moreover, co-occurrence network analysis revealed that diverse keystone taxa may drive pyrene dissipation via more positive interactions between keystone species and with other species in soil. Taken together, these findings provide new insights on the regulation of keystone taxa diversity to promote the dissipation of PAH in soil.

The potential role of fertilizer-derived exogenous bacteria on soil bacterial community assemblage and network formation

Manure fertilization contributes to crop production and sustainable agriculture by introducing large amounts of nutrients and exogenous microbes into soil. However, the contribution of exogenous microbes in shaping soil bacterial community and network structure after fertilization are still controversial. In this study, bacterial communities and network structure that received unsterilized (R + C) or sterilized (R + SC) manure fertilizers, as well as no fertilizer control (R), were characterized using high throughput sequencing. Results showed that the relative abundance of fertilizer-derived OTUs decreased from 10.4% to 4.6% after 90 days incubation, while the Bray-Curtis distance between the control and fertilization group (R + C and R + SC) gradually increased with the culture time. It can be supposed that manure fertilization altered soil bacterial communities by interfering the growth of indigenous bacteria rather than the colonization of fertilizer-derived bacteria. Network analysis showed that a subset of the fertilizer-derived OTUs identified as Xanthomonadales order and Promicromonospora, Constrictibacter genera acted as connectors between modules. They enhanced the interactions not only between soil-derived OTUs and fertilizer-derived OTUs, but also within indigenous bacteria, supported that the introduction of fertilizer-derived exogenous bacteria contributes large to soil bacterial network association. Moreover, fertilizer-derived OTUs presented to be positively correlated with soil pH, while majority soil-derived OTUs presented to be negatively correlated with various physicochemical variables (pH, DOC, NO₃^-, and LAP). Our study highlighted the critical role of fertilizer-derived bacteria in regulating indigenous soil microbial community and network formation after fertilization.

Variation in the diazotrophic community in a vertical soil profile contaminated with antimony and arsenic

A nitrogen (N) deficiency will usually hinder bioremediation efforts in mining-derived habitats such as occurring in mining regions. Diazotrophs can provide N to support the growth of plants and microorganisms in these environments. However, diazotrophic communities in mining areas have been not studied frequently and are more poorly understood than those in other environments, such as in agricultural soils or in the presence of legumes. The current study compares the differences in depth-resolved diazotrophic community compositions and interactions in two contrasting sites (to depths of 2 m), including a highly contaminated and a moderately contaminated site. Antimony (Sb) and arsenic (As) co-contamination induced a loosely connected biotic interaction, and a selection of deep soils by diazotrophic communities. Multiple lines of evidence, including the enrichment of diazotrophic taxa in the highly contaminated sites, microbe-microbe interactions, environment-microbe interactions, and a machine learning approach (random forests regression), demonstrated that Rhizobium was the keystone taxon within the vertical profile of contaminated soil and was resistant to the Sb and As contaminant fractions. All of these observations suggest that one diazotroph, Rhizobium, may play an important role in N fixation in the examined contaminated sites.

Elevation-related climate trends dominate fungal co-occurrence network structure and the abundance of keystone taxa on Mt. Norikura, Japan

Soil fungi play an important role in promoting nutrient cycling and maintaining ecosystem stability. Yet, there has been little understanding of how fungal co-occurrence networks differ along elevational climate gradients, a topic of interest to both macroecology and climate change studies. Based on high-throughput sequencing technology, we investigated the trend in co-occurrence network structure of soil fungal communities at 11 elevation levels along a 2300 m elevation gradient on Mt. Norikura, Japan, and identified the keystone taxa in the network, hypothesizing a progressive decline in network connectivity with elevation due to decreased plant diversity and enhanced environmental stress caused by changes in climate and soil characteristics. Our results demonstrated that network-level topological features such as network size, average degree, clustering coefficient, and modularity decreased significantly with increasing elevation, indicating that the fungal OTUs at low elevation were more closely associated and the network structure was more compact at low elevations. This conclusion was verified by the negative correlation between positive cohesion, negative cohesion and elevation. Moreover, the negative/positive cohesion ratio reached its peak value in mid-elevations with moderate environmental stress, indicating that the fungal community structure in mid-elevations was more stable than that at other elevations. We also found that the keystone taxa were more abundant at lower elevations. Furthermore, statistical analysis revealed that against a background of uniform geology, climate may play a dominant role in determining the properties and intensity of soil fungal networks, and significantly affect the abundance distribution of keystone taxa. These findings enhance understanding of the pattern and mechanism of the fungal community co-occurrence network along elevation, as well as the responses of microorganisms to climate change on a vertical scale in montane ecosystems. IMPORTANCE: Exploration of the elevational distribution of microbial networks and their driving factors and mechanisms may provide opportunities for predicting potential impacts of environmental changes, on ecosystem functions and biogeographic patterns at a broad scale. Although many studies have explored patterns of fungal community diversity and composition along various environmental gradients, it is unclear how the topological structure of co-occurrence networks shifts along elevational temperature gradients. In this study, we found that the connectivity of the fungal community decreased with increasing elevation and that climate was the dominant factor regulating co-occurrence patterns, apparently acting indirectly through soil characteristics. Our results also suggest that higher elevations on mountains have fewer keystone taxa than low elevations. These patterns may be related to the decrease of plant diversity and the increase of environmental stress along elevation gradients.

Impact of carbon/nitrogen ratio on the performance and microbial community of sequencing batch biofilm reactor treating synthetic mariculture wastewater

The differences of cultured organism species, aquaculture model and supervisor mode lead to different carbon/nitrogen ratios in mariculture wastewater. Therefore, the performance, microbial community and enzymatic activity of sequencing batch biofilm reactor were compared in treating synthetic mariculture wastewater at different chemical oxygen demand/nitrogen (COD/N) ratios. Compared with COD/N ratio of 6, the ammonia-oxidizing rate and nitrite-oxidizing rate at COD/N ratio of 5, 4 and 3 increased by 3.66 % and 3.08 %, 11.19 % and 14.95 %, and 24.50 % and 32.54 %, respectively. Similarly, the ammonia monooxygenase and nitrite oxidoreductase activities increased by 3.50 % and 6.76 %, 11.09 % and 16.22 %, and 25.43 % and 39.19 % at COD/N ratio at 5, 4 and 3, respectively. However, the denitrifying rate and denitrification enzymatic activity declined with the decrease of C/N ratio from 6 to 3. The production, protein content and polysaccharide content of loosely bound extracellular polymeric substances (LB-EPS) and tightly bound EPS (TB-EPS) reduced with the decrease of COD/N ratio from 6 to 3. The abundance of nitrifying genera increased with the decrease of COD/N ratio from 6 to 3, whereas most of denitrification genera displayed a decreasing trend. The microbial co-occurrence pattern, keystone taxa and significant difference were altered with the decrease of COD/N ratio. Among the keystone taxa, Thauera, Denitromonas, Nitrosomonas and Denitratisoma had a close link with nitrogen transformation. The present results can provide some theoretical basis for evaluating the effect of carbon/nitrogen ratio on the nitrogen removal of biological wastewater treatment systems.

Microbial community profiles in soils adjacent to mining and smelting areas: Contrasting potentially toxic metals and co-occurrence patterns

Mining and smelting activities have introduced severe potentially toxic metals (PTMs) contamination into surrounding soil settings. Influences of PTMs on microbial diversity have been widely studied. However, variations of microbial communities, network structures and community functions in different levels of PTMs contaminated soils adjacent to mining and smelting aera are still poorly investigated. In this study, microbial communities of soils around different levels of PTMs contamination were comprehensively studied by 16S rRNA gene amplicons high-throughput sequencing. Microbial interactions and module functions were also exploited to ascertain the discrepancies of PTMs concentration levels on microbial ecological functions. Results indicated that the microbial community composition was significantly distinct attributed to the phylum Protebacteria (p = 0.002) dominating in soil with high level PTMs contents but Actinobacteria (p = 0.002) in low level of PTMs-contaminated soil. Microbial α diversity was not significantly influenced by different levels of PTMs contaminations. Microorganisms proactively responded to PTMs content levels by means of strengthening network complexities and modularities among microbe-microbe interactions. The functions of main network modules were predicted associating membrane transport, amino acid metabolism, energy metabolism and carbohydrate metabolism. The PTMs detoxification and anti-oxidation were significantly strengthened at the high level of PTMs contamination. The present study demonstrated that modification of microbial community by the adaptive adjustment of microbial compositions and strengthening their network complexity and modularity.

Metagenomic insights into the effect of thermal hydrolysis pre-treatment on microbial community of an anaerobic digestion system

Thermal hydrolysis process (THP) is an effective pre-treatment method to reduce solids volume and improve biogas production during anaerobic digestion (AD) via increasing the biodegradability of waste activated sludge (WAS). However, the effects of THP pre-treated sludge on microbial diversity, interspecies interactions, and metabolism in AD systems remain largely unknown. We therefore setup and operated an anaerobic digester during a long-term period to shed light on the effect of THP pre-treatment on AD microbial ecology in comparison to conventional AD via Illumina based 16S rRNA gene amplicon sequencing and genome-centric metagenomics analysis. Results showed THP sludge significantly reduced the microbial diversity, shaped the microbial community structure, and resulted in more intense microbial interactions. Compared to WAS as the feed sludge, THP sludge shaped the core functional groups, but functional redundancy ensured the system's stability. The metabolic interactions between methanogens and syntrophic bacteria as well as the specific metabolic pathways were further elucidated. Hydrogenotrophic methanogens, Methanospirillum sp. and Methanolinea sp., were the primary contributors for methane production when treating THP and WAS, respectively, which also have potential for acetate oxidation to methane. Collectively, this study provides in-depth information on the interspecies interactions to better understand how THP pre-treatment influences AD microbial community.

Profiling of Microbial Communities in the Sediments of Jinsha River Watershed Exposed to Different Levels of Impacts by the Vanadium Industry, Panzhihua, China

The mining, smelting, manufacturing, and disposal of vanadium (V) and associated products have caused serious environmental problems. Although the microbial ecology in V-contaminated soils has been intensively studied, the impacted watershed ecosystems have not been systematically investigated. In this study, geochemistry and microbial structure were analyzed along ~30 km of the Jinsha River and its two tributaries across the industrial areas in Panzhihua, one of the primary V mining and production cities in China. Geochemical analyses showed different levels of contamination by metals and metalloids in the sediments, with high degrees of contamination observed in one of the tributaries close to the industrial park. Analyses of the V4 hypervariable region of 16S rRNA genes of the microbial communities in the sediments showed significant decrease in microbial diversity and microbial structure in response to the environmental gradient (e.g., heavy metals, total sulfur, and total nitrogen). Strong association of the taxa (e.g., Thauera, Algoriphagus, Denitromonas, and Fontibacter species) with the metals suggested selection for these potential metal-resistant and/or metabolizing populations. Further co-occurrence network analysis showed that many identified potential metal-mediating species were among the keystone taxa that were closely associated in the same module, suggesting their strong inter-species interactions but relative independence from other microorganisms in the hydrodynamic ecosystems. This study provided new insight into the microbe-environment interactions in watershed ecosystems differently impacted by the V industries. Some of the phylotypes identified in the highly contaminated samples exhibited potential for bioremediation of toxic metals (e.g., V and Cr).

Keystone taxa-mediated bacteriome response shapes the resilience of the paddy ecosystem to fungicide triadimefon contamination

The increasing input of fungicides has emerged as a global concern for agroecosystem stability and sustainability. Agroecosystem resilience has been linked to microbiome response, however, is not well understood. Focusing on a widespread triazole-class fungicide triadimefon in the paddy ecosystem, we characterized that the soils and sediments were dominant triadimefon reservoirs with the peak level at 195 μg kg^-1 and 31.3 μg kg^-1, respectively, but essential for the resilience of paddy ecosystem to triadimefon. In paddy simulation models, the half-life of triadimefon in soil-sediment was 8.4-28.9 days, while it was prolonged to 86.6-115.5 days after elimination of resident microbial community. Phospholipid fatty acid profiling and high-throughput sequencing showed that the distinctive bacterial community responses contributed to variable degradation of triadimefon in paddy soils and sediments. Sphingomonas and Xanthomonas were identified as positive responders of the keystone taxa in the responsive bacteriome, whereas Enterobacter were negative responders that declined over time. Synthetic assemblages combined with quantitative polymerase chain reaction further validated that Sphingomonas and Xanthomonas were involved in sustaining soil-sediment resilience to triadimefon contamination. Collectively, our results revealed that the shaping of soil and sediment bacteriomes was responsible for the resilience of the paddy agroecosystem to fungicide contamination.

Feedback mechanisms of periphytic biofilms to ZnO nanoparticles toxicity at different phosphorus levels

The increasing use of nanoparticles (NPs) has raised concerns about their potential environmental risks. Many researches on NPs focused on the toxicity mechanism to microorganisms, but neglect the toxicity effects in relation to nutritional conditions. Here, we evaluated the interactive effects of zinc oxide (ZnO) NPs and phosphorus (P) levels on the bacterial community and functioning of periphytic biofilms. Results showed that long-term exposure to ZnO NPs significantly reduced alkaline phosphatase activity (APA) of periphytic biofilms just in P-limited conditions. Co-occurrence network analysis indicated that ZnO NPs exposure reduced network complexity between bacterial taxa in P-limited conditions, while the opposite trend was observed in P-replete conditions. Correlation analysis and random forest modeling suggested that excessive Zn²⁺ released and high reactive oxygen species (ROS) production might be mainly responsible for the inhibition of APA induced by ZnO NPs under P-limited conditions, while adjustment of bacterial diversity and improvement of keystone taxa cooperation were the main mechanisms in maintaining APA when subjected to weak toxicity of ZnO NPs in P-replete conditions. Taken together, our results provide insights into the biological feedback mechanism involved in ZnO NPs exposure on the ecological function of periphytic biofilms in different P nutritional conditions.