Fluvial network organization imprints on microbial co-occurrence networks

Impact of extreme rainfall and flood events on harmful cyanobacterial communities and ecological safety in the Baiyangdian Lake Basin, China

Globally, climate change has intensified extreme rainfall events, leading to substantial hydrological changes in aquatic ecosystems. These changes, in turn, have increased the frequency of harmful algal blooms, particularly those of cyanobacteria. This study examines cyanobacterial community dynamics in the Baiyangdian Lake Basin, China, after heavy rainfall and flooding events. The aim was to clarify how such extreme hydrological events affect cyanobacterial populations in floodplain ecosystems and assess related ecological risks. The results demonstrated a significant increase in cyanobacterial diversity, exemplified by an increase of the Shannon diversity index from an average of 1.72 to 2.1 (p < 0.05). Following heavy rainfall and subsequent flooding, the average relative abundance of cyanobacteria in the microbial community increased from 7.59 % to 9.62 %, along a notable rise in the abundance of harmful cyanobacteria. The community structure exhibited notable differences after flooding, showing an increase in species richness, but a decrease in community tightness and clustering, as well as a reduction in niche overlap among harmful cyanobacteria. Environmental factors such as dissolved oxygen, water temperature, and pH were identified as crucial predictors of harmful cyanobacterial community differences and abundance variations resulting from flooding. These findings provide a critical framework for predicting ecological risks associated with the expansion of bloom-forming cyanobacteria in large shallow lake basins, particularly under intensified rainfall and flooding events. This insight is essential to anticipate potential ecological disruptions in sensitive aquatic ecosystems.

Stronger biogeographical pattern of bacterioplankton communities than biofilm communities along a riverine ecosystem: a local scale study of the Kaidu river in the arid and semi-arid northwest of China

Although the biogeographical pattern and mechanisms underlying microbial assembly have been well-explored in lentic ecosystems, the relevant scenarios in lotic ecosystems remain poorly understood. By sequencing the bacterial communities in bacterioplankton and biofilm, our study detected their distance-decay relationship (DDR), and the balance between deterministic and stochastic processes, along the Kaidu river in an arid and semi-arid region of northwest China. Our results revealed that bacterioplankton and biofilm had significantly contrasting community structures. The bacterioplankton communities showed a gradually decreasing trend in alpha-diversity from the headwater to the river mouth, contrasting with the alpha-diversity of biofilm communities which was constant along the river length. Both bacterioplankton and biofilm showed significant DDRs along the 500-km river corridor with the slope of the bacterioplankton DDR being steeper than that of the biofilm DDR, which implies a stronger biogeography of bacterioplankton than biofilm. Relative to biofilm communities, the species interactions formed a denser and more complex network in the bacterioplankton communities than in the biofilm communities. Our results also revealed that there was a transition of community assembly from deterministic to stochastic processes upstream to downstream, although both the bacterioplankton and biofilm communities were mainly regulated by deterministic processes within the entire river. All these empirical results expand our knowledge of microbial ecology in an arid and semi-arid lotic ecosystem.

The shifts in microbial interactions and gene expression caused by temperature and nutrient loading influence Raphidiopsis raciborskii blooms

Climate change and the trophic status of water bodies are important factors in global occurrence of cyanobacterial blooms. The aim of this study was to explore the cyanobacteria‒bacterial interactions that occur during Raphidiopsis raciborskii (R. raciborskii) blooms by conducting microcosm simulation experiments at different temperatures (20 °C and 30 °C) and with different phosphorus concentrations (0.01 mg/L and 1 mg/L) using an ecological model of microbial behavior and by analyzing microbial self-regulatory strategies using weighted gene coexpression network analysis (WGCNA). Three-way ANOVA revealed significant effects of temperature and phosphorus on the growth of R. raciborskii (P < 0.001). The results of a metagenomics-based analysis of bacterioplankton revealed that the synergistic effects of both climate and trophic changes increased the ability of R. raciborskii to compete with other cyanobacteria for dominance in the cyanobacterial community. The antagonistic effects of climate and nutrient changes favored the occurrence of R. raciborskii blooms, especially in eutrophic waters at approximately 20 °C. The species diversity and richness indices differed between the eutrophication treatment group at 20 °C and the other treatment groups. The symbiotic bacterioplankton network revealed the complexity and stability of the symbiotic bacterioplankton network during blooms and identified the roles of key species in the network. The study also revealed a complex pattern of interactions between cyanobacteria and non-cyanobacteria dominated by altruism, as well as the effects of different behavioral patterns on R. raciborskii bloom occurrence. Furthermore, this study revealed self-regulatory strategies that are used by microbes in response to the dual pressures of temperature and nutrient loading. These results provide important insights into the adaptation of microbial communities in freshwater ecosystems to environmental change and provide useful theoretical support for aquatic environmental management and ecological restoration efforts.

Vertical heterogeneity enhances network complexity and stability of co-occurrence microbes in the eastern Indian Ocean

Microbes are core to driving biogeochemical cycles and differ between sun-drenched surface and relatively dark deep oceans. However, their distinct contributions to the organization and association of communities are still remaining elusive. Here, their assembly and co-occurrence stability are systematically researched along the surface and vertical gradients in the eastern Indian Ocean. The distribution of surface microbes was grouped tightly with closer phylogenetic distance and broader niche breadth, and separately from those vertical samples. Clear distance-decay of community similarity was observed in surface microbes with lower richness, while more diverse microeukaryotes and prokaryotes were observed in surface and vertical environments, respectively. Co-occurrence microbes along vertical gradients had a more complex network that was dominated by prokaryotes, while exhibited a lower modularity compared to the surface network. Microbial associations along vertical gradients were more stable and resilient, with lower robustness, higher vulnerability, and a relatively consistent fragmentation. Moreover, prokaryotes contribute greatly to the network topology and stability compared to microeukaryotes in surface environments, emphasizing their distinct functions and survival strategies in maintaining community stability across spatial variations. Environmental selection and community differentiation led to the divergence in organization and potential function of microbes. This study shed light on new perspectives on how marine microbes were associated with and influenced by spatial heterogeneity and their distinct roles in community organization in the face of environmental fluctuations.

Anthropogenic activity shapes the assemble and co-occurrence pattern of microbial communities in fishing harbors around the Bohai economic circle

This study aimed to elucidate the effects of coastal environmental stress on the composition of sediment bacterial communities and their cooccurrence patterns in fishing harbors around the Bohai Economic Circle, China. Compared with the natural sea area, fishing harbors contained higher levels of organic pollution (organic pollution index = 0.12 ± 0.026) and considerably reduced bacterial richness and evenness. The distributions of sediment microbial communities clustered along the pollutant concentration gradients across fishing harbors. Betaproteobacteria dominated (76%) organically polluted fishing harbors, which were mostly disturbed by anthropogenic activities. However, the harbors also revealed the absence of numerous pathogenic (Coxiella and Legionella) and photosynthetic (Synechococcus and Leptolyngbya) bacteria. Abundant genera, including Thiobacillus and Arenimonas, exhibited a positive correlation with total phosphorus and a negative correlation with total nitrogen in sediments. Meanwhile, Sulfurovum, Psychrobacter, and Woeseia showed the opposite trend. Pollutant accumulation and anthropogenic activities caused the decrease in the sediment microbial diversity and dispersal ability and promoted convergent evolution. Severely polluted harbors with simplified cooccurrence networks revealed the presence of destabilized microbial communities. In addition, the modularity of bacterial networks decreased with organic pollution. Our results provide important insights into the adjustment mechanism of microbial communities to community organization and functions under environmental pollution stress. Overall, this study enhanced our understanding of how microbial communities in coastal sediments adapted and survived amidst anthropogenic activities like oily effluent discharges from large ships, wash water, domestic sewage, garbage, and fisheries wastes. It also examined their resilience to future contamination.

Networks as tools for defining emergent properties of microbiomes and their stability

The potential promise of the microbiome to ameliorate a wide range of societal and ecological challenges, from disease prevention and treatment to the restoration of entire ecosystems, hinges not only on microbiome engineering but also on the stability of beneficial microbiomes. Yet the properties of microbiome stability remain elusive and challenging to discern due to the complexity of interactions and often intractable diversity within these communities of bacteria, archaea, fungi, and other microeukaryotes. Networks are powerful tools for the study of complex microbiomes, with the potential to elucidate structural patterns of stable communities and generate testable hypotheses for experimental validation. However, the implementation of these analyses introduces a cascade of dichotomies and decision trees due to the lack of consensus on best practices. Here, we provide a road map for network-based microbiome studies with an emphasis on discerning properties of stability. We identify important considerations for data preparation, network construction, and interpretation of network properties. We also highlight remaining limitations and outstanding needs for this field. This review also serves to clarify the varying schools of thought on the application of network theory for microbiome studies and to identify practices that enhance the reproducibility and validity of future work. Video Abstract.

Community organization and network complexity and stability: contrasting strategies of prokaryotic versus eukaryotic microbiomes in the Bohai Sea and Yellow Sea

Unraveling the effects of spatial gradients on microbiome assembly and association is a challenging topic that remains understudied in the coastal ecosystem. Here, we aimed to investigate the effects of spatial variation on the network complexity and stability of plankton microbiomes in the Bohai Sea and Yellow Sea. These seas serve as spawning and nursery grounds for economically important fisheries valued at billions of dollars annually. Environmental heterogeneity structures microbial communities into distinct spatial patterns, leading to complex direct/indirect relationships and broader ecological niches of bacterioplankton compared to microeukaryotic communities. Interestingly, salinity gradients positively influenced the richness of rare subgroups of bacterioplankton, while the rare microeukaryotic subgroups showed an opposite trend. Abundant subgroups of prokaryotic/eukaryotic microbiomes exhibited greater environmental niche breadth and lower phylogenetic distance compared to the rare subgroups. Stochastic processes contributed greatly to microbiome dynamics, and deterministic processes governed the bacterioplankton organization with a lower phylogenetic turnover rate. Compared to microeukaryotes, bacterioplankton exhibit higher network modularity, complexity, and robustness and lower fragmentation, and vulnerability. These observations offer vital insights into the anti-interference ability and resistance of plankton microbiomes in response to environmental gradients in terms of organization and survival strategy as well as their adaptability to environmental disturbances.IMPORTANCEAn in-depth understanding of community organization and stability of coastal microbiomes is crucial to determining the sustainability of marine ecosystems, such as the Bohai Sea and Yellow Sea. Distinct responses between prokaryotic and eukaryotic microbiomes to spatial heterogeneity were observed in terms of geographical distribution, phylogenetic distance, niche breadth, and community assembly process. Environmental variations are significantly correlated with the dynamics of rare eukaryotic plankton subcommunities compared to prokaryotic plankton subcommunities. Deterministic processes shaped prokaryotic plankton community organization with a lower phylogenic turnover rate. Rare subgroups had noticeably higher phylogenetic distance and lower niche breadth than the corresponding abundant subgroups. Prokaryotic microbiomes had higher molecular network complexity and stability compared to microeukaryotes. Results presented here show how environmental gradients alter both the geographical characteristics of the microbial organization in coastal seas and also their co-occurrence network complexity and stability and thus have critical implications for nutrient and energy cycling.

Dispersal limitation determines the ecological processes that regulate the seasonal assembly of bacterial communities in a subtropical river

Bacteria play a crucial role in pollutant degradation, biogeochemical cycling, and energy flow within river ecosystems. However, the underlying mechanisms governing bacterial community assembly and their response to environmental factors at seasonal scales in subtropical rivers remain poorly understood. In this study, we conducted 16S rRNA gene amplicon sequencing on water samples from the Liuxi River to investigate the composition, assembly processes, and co-occurrence relationships of bacterial communities during the wet season and dry season. The results demonstrated that seasonal differences in hydrochemistry significantly influenced the composition of bacterial communities. A more heterogeneous community structure and increased alpha diversity were observed during the dry season. Water temperature emerged as the primary driver for seasonal changes in bacterial communities. Dispersal limitation predominantly governed community assembly, however, during the dry season, its contribution increased due to decreased immigration rates. Co-occurrence network analysis reveals that mutualism played a prevailing role in shaping bacterial community structure. Compared to the wet season, the network of bacterial communities exhibited higher modularity, competition, and keystone species during the dry season, resulting in a more stable community structure. Although keystone species displayed distinct seasonal variations, Proteobacteria and Actinobacteria were consistently abundant keystone species maintaining network structure in both seasons. Our findings provide insights into how bacterial communities respond to seasonal environmental changes, uncovering underlying mechanisms governing community assembly in subtropical rivers, which are crucial for the effective management and conservation of riverine ecosystems.

Comparing Sediment Bacterial Communities of Volcanic Lakes and Surrounding Rivers in Inner Mongolia Autonomous Region, Northeastern China

Volcanic lakes originate from a volcanic crater or caldera, and were a crucial component of aquatic ecosystems. Sediment bacteria play an important role in the nutrient cycling of aquatic ecosystems; however, their patterns distribution in volcanic lakes and the surrounding river habitats are unknown. In this study, we compare the sediment bacterial communities and their co-occurrence networks between these two habitats in the Inner Mongolia Autonomous Region, Northeastern China (the Arxan UNESCO Global Geopark), using 16S rRNA gene amplicon sequencing. The results revealed that there were significant variations in the physicochemical parameters of the sediment between these two habitats. The bacterial α-diversity, β-diversity, and community composition of the sediment also significantly differed between these two habitats. Network analysis showed that the co-occurrence patterns and keystone taxa in the sediment differed between these two habitats. The sediment bacterial communities in the river habitats were more stable than those in the lake habitats in the face of environmental change. Canonical correspondence analysis demonstrated that both physical (pH and MC) and nutrition-related factors (TN, TP, LOI, and TOC) were the most important environmental factors shaping the variations of bacterial community composition (BCC) in the sediment between these two habitats. This work could greatly improve our understanding of the sediment BCC of the sediment from aquatic ecosystems in the UNESCO Global Geopark.

Forest management positively reshapes the phyllosphere bacterial community and improves community stability

Research has shown that forest management can improve the post-drought growth and resilience of Qinghai spruce in the eastern Qilian Mountains, located on the northeastern Tibetan Plateau. However, the impact of such management on the tree-associated phyllosphere microbiome is not yet fully understood. This study provides new evidence of positive forest management effects on the phyllosphere microbiome after extreme drought, from the perspectives of community diversity, structure, network inference, keystone species, and assembly processes. In managed Qinghai spruce forest, the α-diversity of the phyllosphere bacterial communities increased, whereas the β-diversity decreased. In addition, the phyllosphere bacterial community became more stable and resistant, yet less complex, following forest management. Keystone species inferred from a bacterial network also changed under forest management. Furthermore, forest management mediated changes in community assembly processes, intensifying the influence of determinacy, while diminishing that of stochasticity. These findings support the hypothesis that management can re-assemble the phyllosphere bacterial community, enhance community stability, and ultimately improve tree growth. Overall, the study highlights the importance of forest management on the phyllosphere microbiome and furnishes new insights into forest conservation from the perspective of managing microbial processes and effects.

Strategies for regulating the intensity of different cyanobacterial blooms: Insights from the dynamics and stability of bacterioplankton communities

The occurrence of cyanobacterial blooms is increasing in frequency and magnitude due to climate change and human activities, which poses a direct threat to drinking water security. The impacts of abiotic and biotic factors on the development of blooms have been well studied; however, control strategies for different bloom intensities have rarely been explored from the perspective of the dynamics and stability of bacterioplankton communities. Here, a network analysis was used to investigate the interactions and stability of microbial communities during different periods of R. raciborskii bloom in an inland freshwater lake. The abundance and diversity of rare taxa were significantly higher than that of abundant taxa throughout the bloom cycle. At the pre-bloom (PB) stage, microbial interactions among the different bacterial groups were weak but strongly negatively correlated, indicating low robustness and weak disturbance resistance within the community. However, community stability was better, and microbial interactions became more complicated at the high-bloom (HB) and low-bloom (LB) stages. Interestingly, rare taxa were significantly responsible for community stability and connectivity despite their low relative abundance. The Mantel test revealed that Secchi depth (SD), orthophosphate (PO43--P), and dissolved oxygen (DO) were significantly positively correlated with abundant taxa, rare taxa and PB. DO was significantly positively correlated with HB, intermediate taxa, and rare taxa, while water temperature (WT), N/P and total nitrogen (TN) were significantly positively correlated with LB, abundant taxa, intermediate taxa, and rare taxa. These findings suggest that reducing the PO43--P concentration at the PB stage may be an effective approach to preventing the development of R. raciborskii blooms, while regulating rare taxa at the HB and LB stages may be a key factor in controlling R. raciborskii blooms.

Towards (better) fluvial meta-ecosystem ecology: a research perspective

Rivers are an important component of the global carbon cycle and contribute to atmospheric carbon exchange disproportionately to their total surface area. Largely, this is because rivers efficiently mobilize, transport and metabolize terrigenous organic matter (OM). Notably, our knowledge about the magnitude of globally relevant carbon fluxes strongly contrasts with our lack of understanding of the underlying processes that transform OM. Ultimately, OM processing en route to the oceans results from a diverse assemblage of consumers interacting with an equally diverse pool of resources in a spatially complex network of heterogeneous riverine habitats. To understand this interaction between consumers and OM, we must therefore account for spatial configuration, connectivity, and landscape context at scales ranging from local ecosystems to entire networks. Building such a spatially explicit framework of fluvial OM processing across scales may also help us to better predict poorly understood anthropogenic impacts on fluvial carbon cycling, for instance human-induced fragmentation and changes to flow regimes, including intermittence. Moreover, this framework must also account for the current unprecedented human-driven loss of biodiversity. This loss is at least partly due to mechanisms operating across spatial scales, such as interference with migration and habitat homogenization, and comes with largely unknown functional consequences. We advocate here for a comprehensive framework for fluvial networks connecting two spatially aware but disparate lines of research on (i) riverine metacommunities and biodiversity, and (ii) the biogeochemistry of rivers and their contribution to the global carbon cycle. We argue for a research agenda focusing on the regional scale-that is, of the entire river network-to enable a deeper mechanistic understanding of naturally arising biodiversity-ecosystem functioning coupling as a major driver of biogeochemically relevant riverine carbon fluxes.

Microbial community structure and co-occurrence network stability in seawater and microplastic biofilms under prometryn pollution in marine ecosystems

Prometryn has been extensively detected in marine environment because of its widespread usage in agriculture and aquaculture and has been concerns since its serious effects on aquatic organisms. However, its impact on the microbial community in the marine ecosystem including seawater and biofilm is still unclear. Therefore, a short-term indoor microcosm experiment of prometryn exposure was conducted. This study found that prometryn had a more significant impact on the structure and stability of the microbial community in seawater compared to microplastic biofilms. Additionally, we observed that the assembly of the microbial community in biofilms was more affected by stochastic processes than in seawater under the exposure of prometryn. Our study provided evidence for the increasing impact of the microbial communities under the stress of prometryn and microplastics.

The role of agricultural drainage, storm-events, and natural filtration on the biogeochemical cycling capacity of aquatic and sediment environments in Lake Erie's drainage basin

Lake Erie is the most at risk of the Great Lakes for degraded water quality due to non-point source pollution caused by agricultural activities in the lake's watershed. The extent and temporal patterns of nutrient loading from these agricultural activities is influenced by the timing of agronomic events, precipitation events, and water flow through areas of natural filtration within the watershed. Downstream impacts of these nutrient loading events may be moderated by the co-loading of functionally relevant biogeochemical cycling microbial communities from agricultural soils. This study quantified loading patterns of these communities from tile drain sources, assessed whether functional communities from agricultural sources influenced downstream microbial functionality, and investigated how distance from agricultural sources, storm events, and areas of natural filtration altered nutrient cycling and nutrient fluxes in aquatic and sediment environments. Water and sediment samples were collected in the Wigle Creek watershed in Ontario, from tile drains through to Lake Erie, from May to November 2021, and microbial nitrogen (N) and phosphorous (P) cycling capacity (quantitative PCR), and nutrient levels were evaluated. Results showed that N and P functional groups were co-loaded with nutrients, with increased loading occurring during storm events and during agricultural activities including fertilization and harvest. Overall functional capacity in the aquatic environment decreased with distance from the agricultural sources and as water transited through natural filtration areas. In contrast, the sediment environment was more resilient to both agricultural disturbances and abiotic factors. This study expands our understanding of when and where different stages of N and P cycling occurs in agriculturally impacted watersheds, and identifies both seasons and regions to target with nutrient mitigation strategies.

River connectivity determines microbial assembly processes and leads to alternative stable states in river networks

River network is a common form of lotic ecosystems. Variances in river connection modes would form networks with significantly different structures, and further affect aquatic organisms. Microbial communities are vital organisms of river networks, they participate in numerous biogeochemical processes. Identifying associations between microbial community and structural features of river networks are essential for maintaining environmental quality. Thus, dendritic (DRN) and trellised river networks (TRN) were studied by combining molecular biological tools, ecological theory and hydrodynamic calculation. Results illustrated that river connectivity, a vital structural feature exhibiting mass transport ability of river network, increased relative importance of homogeneous selection processes in microbial assembly, which would further shape community with alternative stable states. Between the two researched river networks, DRN possessed higher connectivity, which made homogeneous selection as the driving force in community assembly. The microbial communities in DRN were consisted of species occupying similar ecological niche, and exhibited two alternative stable states, which can decrease influences of environmental disturbance on community composition. On the contrary, lower connectivity of TRN decreased proportions of homogeneous selection in community assembly, which further led to species occupying varied ecological niche. The microbial community exhibited only one stable state, and environmental disturbance would cause loss of ecological niche and significantly alter community composition. This study could provide useful information for the optimization of river connection engineering.

Hydrological disturbances enhance stochastic assembly processes and decrease network stability of algae communities in a highland floodplain system

Floodplains are hotspots for biodiversity research and conservation worldwide. Hydrological disturbances can profoundly influence the ecological processes and functions of floodplain systems by altering key biological groups such as algae communities. However, the impacts of flood disturbance on the assembly processes and co-occurrence patterns of algae communities in floodplain ecosystems are still unclear. To ascertain the response patterns of algae communities to flood disturbance, we characterized planktonic and benthic algae communities in 144 water and sediment samples collected from the Tibetan floodplain during non-flood and flood periods based on 23S ribosomal RNA gene sequencing. Results showed that planktonic algae exhibited higher diversity and greater compositional variations compared with benthic communities after flood disturbance. Flooding promoted algae community homogenization at horizontal (rivers vs. oxbow lakes) and vertical levels (water vs. sediment). Stochastic processes governed the assembly of distinct algae communities, and their ecological impacts were enhanced in response to flooding. In the non-flood period, dispersal limitation (81.78 %) was the primary ecological process driving algae community assembly. In the flood period, the relative contribution of ecological drift (72.91 %) to algae community assembly markedly increased, with dispersal limitation (22.61 %) being less important. Flooding reduced the interactions among algae taxa, resulting in lower network complexity and stability. Compared with the planktonic algae subnetworks, the benthic subnetworks showed greater stability in the face of flooding. Findings of this study broaden our understanding of how algae communities respond to hydrological disturbances from an ecological perspective and could be useful for the management of highland floodplain ecosystems.

Environmental gradients shape microbiome assembly and stability in the East China sea

Microbiomes play a key role in marine ecosystem functioning and sustainability. Their organization and stability in coastal areas, particularly in anthropogenic-influenced regions, however, remains unclear compared with an understanding of how microbial community shifts respond to marine environmental gradients. Here, the assembly and community associations across vertical and horizontal gradients in the East China Sea are systematically researched. The seawater microbial communities possessed higher robustness and lower fragmentation and vulnerability compared to the sediment microbiomes. Spatial gradients act as a deterministic filtering factor for microbiome organization. Microbial communities had lower phylogenetic distance and higher niche breadth in the nearshore and offshore areas compared to intermediate areas. The phylogenetic distance of microbiomes decreased from the surface to the bottom but the niche breadth was enhanced in surface and bottom environments. Vertical gradients destabilized microbial associations, while the community diversity was enhanced. Multivariate regression tree analysis and canonical correspondence analysis indicated that depth, distance from shore, nutrient availability, temperature, salinity, and chlorophyll a, affected the distribution and co-occurrence of microbial groups. Our results highlight the crucial roles of environmental gradients in determining microbiome association and stability. These results improve our understanding of the survival strategies/adaptive mechanisms of microbial communities in response to environmental variation and provide new insights for protecting the ecosystems and maintaining the sustainability of ecological functions.

Distribution heterogeneity of sediment bacterial community in the river-lake system impacted by nonferrous metal mines: Diversity, composition and co-occurrence patterns

Metal(loid) pollution caused by mining activities can affect microbial communities. However, knowledge of the diversity, composition, and co-occurrence patterns of bacterial communities in aquatic systems impacted by nonferrous metal mines. Here, the metal(loid) contents and bacterial communities in sediments from the Zijiang River (tributary to mainstream) to Dongting Lake were investigated by geochemical and molecular biology methods. The results indicated that the river sediments had lower pH and higher ecological risk of metal(loid)s than the lake sediment. The diversity and composition of bacterial communities in river sediments significantly (p < 0.05) differed from those in lake sediments, showing distributional heterogeneity. The biomarkers of tributary, mainstream, and lake sediments were mainly members of Deltaproteobacteria, Firmicutes, and Nitrospirae, respectively, reflecting species sorting in different habitats. Multivariate statistical analysis demonstrated that total and bioavailable Sb, As, and Zn were positively correlated with bacterial community richness. pH, TOC, TN, and Zn were crucial factors in shaping the distribution difference of bacterial communities. Environment-bacteria network analysis indicated that pH, SO42-, and total and bioavailable As and Sb greatly influenced the bacterial composition at the genus level. Bacteria-bacteria network analysis manifested that the co-occurrence network in mainstream sediments with a higher risk of metal(loid) pollution exhibited higher modularity and connectivity, which might be the survival mechanism for bacterial communities adapted to metal(loid) pollution. This study can provide a theoretical basis for understanding the ecological status of aquatic systems.

Distribution and co-occurrence networks of the bacterial community in sediment cores from the subtropical Daya Bay, China

The bacterial community plays an important role in biogeochemical cycles in marine sediment. However, little is known about the vertical profiles and co-occurrence patterns of bacterial community in sediment cores from the marine environment. In this study, five sediment cores were taken from a subtropical bay in China, heavily impacted by anthropogenic activities. The bacterial composition in sediment cores was investigated by using high-throughput sequencing of the 16S rRNA gene. A principal coordinates analysis and an adonis analysis of the operational taxonomic unit (OTU) compositions showed that spatial variation, rather than vertical variation, determined the bacterial structure in sediment cores. The bacterial complexity varied greatly across the five sediment cores, and the rare taxa played an important role in supporting the stability of the bacterial network. This study revealed that sediment properties and anthropogenic activities may induce a shift in the bacterial composition in sediment cores of a subtropical bay.

Global patterns and drivers of lead concentration in inland waters

Water bodies are important carriers for lead (Pb) biogeochemical cycling, which is a key pathway of Pb transport. Although existing studies on Pb loading in inland waters have developed rapidly, a quantitative assessment of the distribution patterns and drivers of Pb concentration in inland waters at the global scale remains unclear. Here, by analyzing 1790 observations collected from 386 independent publications, we assessed the spatial distribution and drivers of Pb concentration in inland waters worldwide. We found that (1) globally, the median of Pb concentration in inland waters was 5.81 μg L-1; (2) among different inland water types, Pb concentration was higher in rivers, and the highest Pb concentration was in industrial land in terms of land use type; (3) Pb concentration in inland waters were positively driven by potential evapotranspiration, elevation and road density; and (4) Pb concentration showed a negative relationship with absolute latitude, decreasing from tropic to boreal regions. Overall, our global assessment of the patterns and drivers of Pb concentration in inland waters contributed to a better understanding of the natural and anthropogenic attributions of Pb in the inland hydrological cycling.

Dynamics of dominant rhizospheric microbial communities responsible for trichlorfon absorption and translocation in maize seedlings

In this study, the available phosphorus (AP) and TCF concentrations in soils and maize (Zea mays) seedling tissues were measured in response to escalating TCF concentrations during 216 hr of culture. Maize seedlings growth considerably enhanced soil TCF degradation, reaching the highest of 73.2% and 87.4% at 216 hr in 50 and 200 mg/kg TCF treatments, respectively, and increased AP contents in all the seedling tissues. Soil TCF was majorly accumulated in seedling roots, reaching maximum concentration of 0.017 and 0.076 mg/kg in TCF-50 and TCF-200, respectively. The hydrophilicity of TCF might hinder its translocation to the aboveground shoot and leaf. Using bacterial 16 S rRNA gene sequencing, we found that TCF addition drastically lessened bacterial community interactions and hindered the complexity of their biotic networks in rhizosphere than in bulk soils, leading to the homogeneity of bacterial communities that were resistant or prone to TCF biodegradation. Mantel test and redundancy analysis suggested a significant enrichment of dominant species Massilia belonging to Proteobacteria phyla, which in turn affecting TCF translocation and accumulation in maize seedling tissues. This study provided new insight into the biogeochemical fate of TCF in maize seedling and the responsible rhizobacterial community in soil TCF absorption and translocation.

Copper and temperature shaped abundant and rare community assembly respectively in the Yellow River

Untangling assembly and microbial interaction of abundant and rare microbiota in aquatic ecosystem is pivotal for understanding how community assembly respond to environmental variables and co-occurrence patterns. Here, we explored the assembly mechanisms, their drivers, and species co-occurrence of abundant and rare microbiomes in the Yellow River using 16S rRNA gene sequencing in Lanzhou, China. Here, abundant community was ubiquitous across all sites, whereas rare community was uneven distributed. The richness and community dissimilarity of rare taxa were significantly greater than those of abundant ones. Stochastic processes structured the rare community assembly in spring and winter, while deterministic processes shaped the abundant and rare community assembly in other seasons and all sites. Copper and water temperature mediated the balance between deterministic and stochastic processes of abundant and rare community, respectively. A few abundant taxa with closer relationships frequently occupied central positions and had a great effect on other co-occurrences in the network, while the majority of keystone microbiota were rare microbiome and played a considerable part in maintaining the network structure. Our study provides some ecological proposals for water quality management and ecological stability of the Yellow River. KEY POINTS: • Deterministic process dominated abundant and rare community assembly. • Cu and TW mediated the balance of abundant and rare community assembly respectively. • Abundant taxa had a greater effect on other co-occurrences in the network.

Inter-basin water diversion homogenizes microbial communities mainly through stochastic assembly processes

Inter-basin water transfer is an effective manner to achieve the optimal allocation of water resources, while accompanied by some ecological effects. The responses of microorganisms to water diversion and the ecological processes in regulating the community assembly are still unclear. Taking the eastern route of South-to-North Water Diversion Project as the study area, we investigated the microbial community patterns and the underlying assemblage processes in habitats with different hydrological connectivity, including isolated lakes, connected lakes and man-made canal. The results showed that microbial communities in the canal had higher diversity, lower dissimilarity, weaker compositional variation, and stronger co-occurrence patterns compared with that in the connected and isolated lakes. These findings suggested that the increase of connectivity among natural aquatic habitats due to water diversion can homogenize microbial communities and reduce microbial heterogeneity. The neutral and null models demonstrated the importance of stochastic processes in shaping microbial community assembly. Dispersal limitation and variable selection were the predominant mechanisms structuring microbial communities in the isolated lakes. Due to the homogenized environmental condition and the enhanced hydrologic connectivity in the canal and the connected lakes, microbial communities had higher dispersal capability and ecological drift occurred more frequently in these lotic habitats. The variations in microbial community structure were mainly driven by biotic ecological succession than abiotic factors, with positive and negative cohesion explained 63% and 25% of variability, respectively. Six taxa were considered as the potential introduced microorganisms, which may favor the nutrient biogeochemical cycling and the organic matter degradation, but may also bring ecological risks. Overall, this study provides a deeper understanding of the ecological consequences of inter-basin water diversion, and helps the regulation and management of these projects.

Spatial and temporal change determined co-occurrence networks stability and community assembly processes of epipelagic seawater microbial community in the Nordic Sea

The Nordic Sea has a vital impact on the global climate change, occupies a significant status in the physical oceanography research, on account of its intersection of complex ocean currents. To explore the influence of seasonal and spatial heterogeneity in its epipelagic seawater on the microbial community structure, a total of 54 seawater samples from 18 stations and 3 water layers (0 m, 50 m, 100 m) were collected in the summer of 2017 and the autumn of 2018 from the Norwegian Sea, the Greenland Sea and the vicinity of Jan Mayen Island in the Nordic Sea. Alpha- and Beta- diversity analysis showed that significant differences were found between characteristic bacterial groups in detached or mixed currents of corresponding seasons, as endemic OTUs with seasonal and ocean current characteristics which revealed the existence of spatiotemporal patterns of microbial communities in the Nordic Sea. Moreover, co-occurrence networks were conducted to show different degree of complexity and stability of microbial community response to spatiotemporal dynamic changes. Furthermore, the flow and collision between ocean currents do have an impact on the community assembly processes by affecting the migration and dispersal of microbial communities. This study reflects the response of microbial communities to the spatiotemporal dynamics and reveals the microbial community assembly mechanisms under complex hydrological condition represented in the Nordic Sea.

Assembly and Network Stability of Planktonic Microorganisms under the Influence of Salinity Gradient: an Arctic Case Study from the Lena River Estuary to the Laptev Sea

The diversity and primary productivity in the Arctic ecosystem are rapidly changing due to global warming. Microorganisms play a vital role in biogeochemical cycling. However, the diversity of planktonic microorganism communities in the Laptev Sea, one of the most important marginal seas of the Western Arctic Ocean, have not been studied sufficiently in depth. The diversity and community structure of the planktonic microorganisms in the surface water were investigated at 20 stations on the Lena River flowing into the Laptev Sea. Multivariate statistical analyses demonstrated clear spatial patterns in the α diversity and community structure for microorganisms under different salinity levels. Co-occurrence networks of microbial communities revealed that spatial variation promoted differentiation of the characteristics and stability of microbial networks in the Laptev Sea. Contrary to expectations, abundant taxa were found to not have a large influence on the stability and resilience of microbial interactions in the region. On the contrary, less-abundant taxa were found to have far greater influence. The stability and resilience of the prokaryotic and microeukaryotic networks in the Lena River estuary and the continental shelf provided valuable insights into the impact of freshwater and land inflow disturbances on microbial assemblage. Overall, these results enhance our understanding of the composition of microbial communities and provide insights into how spatial changes of abundant versus rare species alter the nature and stability of microbial networks from the Lena River estuary to the Laptev Sea. In addition, this study explored microbial interactions and their ability to resist future disturbances. IMPORTANCE The regime of the Laptev Sea depends closely on the runoff of the Lena River. Microorganisms are essential components of aquatic food webs and play a significant role in polar ecosystems. In this study, we provided a basic microbial data set as well as new insights into the microbial networks from the Lena River estuary to the Laptev Sea, while exploring their potential to resist future disturbances. A comprehensive and systematic study of the community structure and function of the planktonic microorganisms in the Laptev Sea would greatly enhance our understanding of how polar microbial communities respond to the salinity gradient under climate warming.

Seasonal succession of microbial community co-occurrence patterns and community assembly mechanism in coal mining subsidence lakes

Coal mining subsidence lakes are classic hydrologic characteristics created by underground coal mining and represent severe anthropogenic disturbances and environmental challenges. However, the assembly mechanisms and diversity of microbial communities shaped by such environments are poorly understood yet. In this study, we explored aquatic bacterial community diversity and ecological assembly processes in subsidence lakes during winter and summer using 16S rRNA gene sequencing. We observed that clear bacterial community structure was driven by seasonality more than by habitat, and the α-diversity and functional diversity of the bacterial community in summer were significantly higher than in winter (p < 0.001). Canonical correspondence analysis indicated that temperature and chlorophyll-a were the most crucial contributing factors influencing the community season variations in subsidence lakes. Specifically, temperature and chlorophyll-a explained 18.26 and 14.69% of the community season variation, respectively. The bacterial community variation was driven by deterministic processes in winter but dominated by stochastic processes in summer. Compared to winter, the network of bacterial communities in summer exhibited a higher average degree, modularity, and keystone taxa (hubs and connectors in a network), thereby forming a highly complex and stable community structure. These results illustrate the clear season heterogeneity of bacterial communities in subsidence lakes and provide new insights into revealing the effects of seasonal succession on microbial assembly processes in coal mining subsidence lake ecosystems.

Crop diversity promotes the recovery of fungal communities in saline-alkali areas of the Western Songnen Plain

Introduction

Phytoremediation is an effective strategy for saline land restoration. In the Western Songnen Plain, northeast China, soil fungal community recovery for saline phytoremediation has not been well documented among different cropping patterns. In this study, we tested how rotation, mixture, and monoculture cropping patterns impact fungal communities in saline-alkali soils to assess the variability between cropping patterns.

Methods

The fungal communities of the soils of the different cropping types were determined using Illumina Miseq sequencing.

Results

Mixture and rotation promoted an increase in operational taxonomic unit (OTU) richness, and OTU richness in the mixture system decreased with increasing soil depth. A principal coordinate analysis (PCoA) showed that cropping patterns and soil depths influenced the structure of fungal communities, which may be due to the impact of soil chemistry. This was reflected by soil total nitrogen (TN) and electrical conductivity (EC) being the key factors driving OTU richness, while soil available potassium (AK) and total phosphorus (TP) were significantly correlated with the relative abundance of fungal dominant genus. The relative abundance of Leptosphaerulina, Alternaria, Myrothecium, Gibberella, and Tetracladium varied significantly between cropping patterns, and Leptosphaerulina was significantly associated with soil chemistry. Soil depth caused significant differences in the relative abundance of Fusarium in rotation and mixture soils, with Fusarium more commonly active at 0-15 cm deep soil. Null-model analysis revealed that the fungal community assembly of the mixture soils in 0-15 cm deep soil was dominated by deterministic processes, unlike the other two cropping patterns. Furthermore, fungal symbiotic networks were more complex in rotation and mixture than in monoculture soils, reflected in more nodes, more module hubs, and connectors. The fungal networks in rotation and mixture soils were more stable than in monoculture soils, and mixture networks were obviously more connected than rotations. FUNGuild showed that the relative proportion of saprotroph in rotation and mixture was significantly higher than that in monocultures. The highest proportion of pathotroph and symbiotroph was exhibited in rotation and mixture soils, respectively.

Discussion

Overall, mixture is superior to crop rotation and monocultures in restoring fungal communities of the saline-alkali soils of the Western Songnen Plain, northeast China.

Temporal Succession of Bacterial Community Structure, Co-occurrence Patterns, and Community Assembly Process in Epiphytic Biofilms of Submerged Plants in a Plateau Lake

In shallow macrophytic lakes, epiphytic biofilms are formed on the surface of submerged plant stems and leaves because of algae and bacterial accumulation. Epiphytic biofilms significantly impact the health of the host vegetation and the biogeochemical cycling of lake elements. However, community diversity, species interactions, and community assembly mechanisms in epiphytic bacterial communities (EBCs) of plants during different growth periods are not well understood. We investigated the successional dynamics, co-occurrence patterns, and community assembly processes of epiphytic biofilm bacterial communities of submerged plants, Najas marina and Potamogeton lucens, from July to November 2020. The results showed a significant seasonal variation in EBC diversity and richness. Community diversity and richness increased from July to November, and the temperature was the most important driving factor for predicting seasonal changes in EBC community structure. Co-occurrence network analysis revealed that the average degree and graph density of the network increased from July to November, indicating that the complexity of the EBC network increased. The bacterial community co-occurrence network was limited by temperature, pH, and transparency. The phylogeny-based null model analysis showed that deterministic processes dominated the microbial community assembly in different periods, increasing their contribution. In addition, we found that as the dominance of deterministic processes increased, the microbial co-occurrence links increased, and the potential interrelationships between species became stronger. Thus, the findings provide insights into the seasonal variability of EBC assemblage and co-occurrence patterns in lacustrine ecosystems.

Soil bacterial communities of paddy are dependent on root compartment niches but independent of growth stages from Mollisols of Northeast China

Introduction

Deep insights into adhering soil of root zones (rhizosphere and rhizoplane) microbial community could provide a better understanding of the plant-microbe relationship. To better understand the dynamics of these microbial assemblies over the plant life cycle in rhizodeposition along rice roots.

Methods

Here, we investigated bacterial distribution in bulk, rhizosphere, and rhizoplane soils at tillering, heading, and mature stage, from rice (Oryza sativa) fields of the Northeast China.

Results and Discussion

Our results revealed that soil bacterial α-diversity and community composition were significantly affected by root compartment niches but not by temporal change. Compared to rhizoplane soils in the same period, bulk in the heading and rhizosphere in the mature had the largest increase in Shannon's index, with 11.02 and 14.49% increases, respectively. Proteobacteria, Chloroflexi, Bacteroidetes, and Acidobacteria are predominant across all soil samples, bulk soil had more phyla increased across the growing season than that of root related-compartments. Deterministic mechanisms had a stronger impact on the bacterial community in the compartments connected to the roots, with the relative importance of the bulk soil, rhizoplane and rhizosphere at 83, 100, and 56%, respectively. Because of ecological niche drivers, the bacterial networks in bulk soils exhibit more complex networks than rhizosphere and rhizoplane soils, reflected by more nodes, edges, and connections. More module hub and connector were observed in bulk (6) and rhizoplane (5) networks than in rhizosphere (2). We also detected shifts from bulk to rhizoplane soils in some functional guilds of bacteria, which changed from sulfur and nitrogen utilization to more carbon and iron cycling processes. Taken together, our results suggest distinct bacterial network structure and distribution patterns among rhizosphere, rhizoplane, and bulk soils, which could possibly result in potential functional differentiation. And the potential functional differentiation may be influenced by plant root secretions, which still needs to be further explored.

SCNIC: Sparse correlation network investigation for compositional data

Microbiome studies are often limited by a lack of statistical power due to small sample sizes and a large number of features. This problem is exacerbated in correlative studies of multi-omic datasets. Statistical power can be increased by finding and summarizing modules of correlated observations, which is one dimensionality reduction method. Additionally, modules provide biological insight as correlated groups of microbes can have relationships among themselves. To address these challenges, we developed SCNIC: Sparse Cooccurrence Network Investigation for compositional data. SCNIC is open-source software that can generate correlation networks and detect and summarize modules of highly correlated features. Modules can be formed using either the Louvain Modularity Maximization (LMM) algorithm or a Shared Minimum Distance algorithm (SMD) that we newly describe here and relate to LMM using simulated data. We applied SCNIC to two published datasets and we achieved increased statistical power and identified microbes that not only differed across groups, but also correlated strongly with each other, suggesting shared environmental drivers or cooperative relationships among them. SCNIC provides an easy way to generate correlation networks, identify modules of correlated features and summarize them for downstream statistical analysis. Although SCNIC was designed considering properties of microbiome data, such as compositionality and sparsity, it can be applied to a variety of data types including metabolomics data and used to integrate multiple data types. SCNIC allows for the identification of functional microbial relationships at scale while increasing statistical power through feature reduction.

Bacterial communities in a subarctic stream network: Spatial and seasonal patterns of benthic biofilm and bacterioplankton

Water-column bacterial communities are assembled by different mechanisms at different stream network positions, with headwater communities being controlled by mass effects (advection of bacteria from terrestrial soils) while downstream communities are mainly driven by environmental sorting. Conversely, benthic biofilms are colonized largely by the same set of taxa across the entire network. However, direct comparisons of biofilm and bacterioplankton communities along whole stream networks are rare. We used 16S rRNA gene amplicon sequencing to explore the spatiotemporal variability of benthic biofilm (2 weeks old vs. mature biofilm) and water-column communities at different network positions of a subarctic stream from early summer to late autumn. Amplicon sequence variant (ASV) richness of mature biofilm was about 2.5 times higher than that of early biofilm, yet the pattern of seasonality was the same, with the highest richness in midsummer. Biofilm bacterial richness was unrelated to network position whereas bacterioplankton diversity was negatively related to water residence time and distance from the source. This pattern of decreasing diversity along the network was strongest around midsummer and diminished greatly as water level increased towards autumn. Biofilm communities were phylogenetically clustered at all network positions while bacterioplankton assemblages were phylogenetically clustered only at the most downstream site. Both early and mature biofilm communities already differed significantly between upstream (1st order) and midstream (2nd order) sections. Network position was also related to variation in bacterioplankton communities, with upstream sites harbouring substantially more unique taxa (44% of all upstream taxa) than midstream (20%) or downstream (8%) sites. Some of the taxa that were dominant in downstream sections were already present in the upmost headwaters, and even in riparian soils, where they were very rare (relative abundance <0.01%). These patterns in species diversity and taxonomic and phylogenetic community composition of the riverine bacterial metacommunity were particularly strong for water-column communities, whereas both early and mature biofilm exhibited weaker spatial patterns. Our study demonstrated the benefits of studying bacterioplankton and biofilm communities simultaneously to allow testing of ecological hypotheses about biodiversity patterns in freshwater bacteria.

Bacterioplankton dispersal and biogeochemical function across Alaskan Arctic catchments

In Arctic catchments, bacterioplankton are dispersed through soils and streams, both of which freeze and thaw/flow in phase, seasonally. To characterize this dispersal and its potential impact on biogeochemistry, we collected bacterioplankton and measured stream physicochemistry during snowmelt and after vegetation senescence across multiple stream orders in alpine, tundra, and tundra-dominated-by-lakes catchments. In all catchments, differences in community composition were associated with seasonal thaw, then attachment status (i.e. free floating or sediment associated), and then stream order. Bacterioplankton taxonomic diversity and richness were elevated in sediment-associated fractions and in higher-order reaches during snowmelt. Families Chthonomonadaceae, Pyrinomonadaceae, and Xiphinematobacteraceae were abundantly different across seasons, while Flavobacteriaceae and Microscillaceae were abundantly different between free-floating and sediment-associated fractions. Physicochemical data suggested there was high iron (Fe⁺ ) production (alpine catchment); Fe⁺ production and chloride (Cl^- ) removal (tundra catchment); and phosphorus (SRP) removal and ammonium (NH₄ ⁺ ) production (lake catchment). In tundra landscapes, these 'hot spots' of Fe⁺ production and Cl^- removal accompanied shifts in species richness, while SRP promoted the antecedent community. Our findings suggest that freshet increases bacterial dispersal from headwater catchments to receiving catchments, where bacterioplankton-mineral relations stabilized communities in free-flowing reaches, but bacterioplankton-nutrient relations stabilized those punctuated by lakes.

Legacy Effects of Biochar and Compost Addition on Arbuscular Mycorrhizal Fungal Community and Co-Occurrence Network in Black Soil

Compost and biochar are beneficial soil amendments which derived from agricultural waste, and their application was proven to be effective practices for promoting soil fertility. Arbuscular mycorrhizal (AM) fungi form symbiotic associations with most crop plant species, and are recognized as one group of the most important soil microorganisms to increase food security in sustainable agriculture. To understand the legacy effects of compost and biochar addition on AM fungal communities, a field study was conducted on the Songnen Plain, Northeast China. Two years after application, compost addition improved soil aggregate stability, but we did not detect a legacy effect of compost addition on AM fungal community. Our results indicated that AM fungal Shannon diversity and Pielou evenness indices were significantly increased by one-time biochar addition, but unaffected by compost addition after two year’s application. PERMANOVA analysis also revealed a legacy effect of biochar addition on AM fungal community. Network analysis revealed a dramatically simplified AM fungal co-occurrence network and small network size in biochar added soils, demonstrated by their topological properties (e.g., low connectedness and betweenness). However, AM fungal community did not differ among aggregate fractions, as confirmed by the PERMANOVA analysis as well as the fact that only a small number of AM fungal OTUs were shared among aggregate fractions. Consequently, the current study highlights a stronger legacy effect of biochar than compost addition on AM fungi, and have implications for agricultural practices.

HiLi-chip: A high-throughput library construction chip for comprehensive profiling of environmental microbial communities

Investigating the contribution and associations of environmental microbes to ecological health and human well-being is in great demand with the goal of One Health proposed. To achieve the goal, there is an urgent need for accurate approaches to obtaining a large amount of high-resolution molecular information from various microbes. In this study, we developed a high-throughput library construction chip (HiLi-Chip) for profiling environmental microbial communities and evaluated its performance. The HiLi-Chip showed high conformity with the conventional Pacbio method in terms of α-diversity, community composition of abundant bacteria (>83%), as well as rare taxa (>84%) and human pathogens detection (>67%), indicating its advantages of accuracy, high-throughput, cost-efficiency, and broad practicability. It is suggested that the optimal strategy of the HiLi-Chip was a 2.4 μL PCR mixture per sample (∼2.4 ng DNA) with a 216-sample × 24-replicate format. We have successfully applied the HiLi-Chip to the Jiulongjiang River and identified 51 potential human bacterial pathogens with a total relative abundance of 0.22%. Additionally, under limited nutrients and similar upstream environments, bacteria tended to impose competitive pressures, resulting in a more connected network at the downstream river confluence (RC). Whereas narrow niche breadth of bacteria and upstream environmental heterogeneity probably promoted niche complementary and environment selection leading to fewer links at RC in the midsection of the river. Core bacteria might represent the entire bacterial community and enhance network stability through synergistic interactions with other core bacteria. Collectively, our results demonstrate that the HiLi-Chip is a robust tool for rapid comprehensive profiling of microbial communities in environmental samples and has significant implications for a profound understanding of environmental microbial interactions.

Tillage Practices and Residue Management Manipulate Soil Bacterial and Fungal Communities and Networks in Maize Agroecosystems

Tillage practices and residue management are highly important agricultural practices. However, very few studies have examined the influence of tillage practices and residue management on both bacterial and fungal communities and network patterns in consecutive years. We examined the effects of different tillage practices, including no tillage, rotary tillage, and deep tillage, on soil bacterial and fungal communities and co-occurrence networks following residue removal and residue retention in 2017 and 2018. This study showed that both bacterial and fungal communities were unaffected by tillage practices in 2017, but they were significantly impacted in 2018. Soil fungal operational taxonomic unit (OTU) richness was significantly enhanced by deep tillage compared with no tillage in 2018, while bacterial OTU richness was unaffected in either year. Tillage practices had differing effects on soil microbial co-occurrence networks, with rotary and deep tillage increasing the complexity of bacterial networks but simplifying fungal networks. However, residue retention only induced a shift in the fungal community and simplified soil bacterial and fungal networks in 2018. This study highlights the dissimilar responses of bacterial and fungal networks to tillage practices and emphasizes that tillage practice is more important than residue management in shaping soil microbial communities.

Beyond the snapshot: identification of the timeless, enduring indicator microbiome informing soil fertility and crop production in alkaline soils

BACKGROUND

Microorganisms are known to be important drivers of biogeochemical cycling in soil and hence could act as a proxy informing on soil conditions in ecosystems. Identifying microbiomes indicative for soil fertility and crop production is important for the development of the next generation of sustainable agriculture. Earlier researches based on one-time sampling have revealed various indicator microbiomes for distinct agroecosystems and agricultural practices as well as their importance in supporting sustainable productivity. However, these microbiomes were based on a mere snapshot of a dynamic microbial community which is subject to significant changes over time. Currently true indicator microbiomes based on long-term, multi-annual monitoring are not available.

RESULTS

Here, using samples from a continuous 20-year field study encompassing seven fertilization strategies, we identified the indicator microbiomes ecophysiologically informing on soil fertility and crop production in the main agricultural production base in China. Among a total of 29,184 phylotypes in 588 samples, we retrieved a streamlined consortium including 2% of phylotypes that were ubiquitously present in alkaline soils while contributing up to half of the whole community; many of them were associated with carbon and nutrient cycling. Furthermore, these phylotypes formed two opposite microbiomes. One indicator microbiome dominated by Bacillus asahii, characterized by specific functional traits related to organic matter decomposition, was mainly observed in organic farming and closely associated with higher soil fertility and crop production. The counter microbiome, characterized by known nitrifiers (e.g., Nitrosospira multiformis) as well as plant pathogens (e.g., Bacillus anthracis) was observed in nutrient-deficit chemical fertilizations. Both microbiomes are expected to be valuable indictors in informing crop yield and soil fertility, regulated by agricultural management.

CONCLUSIONS

Our findings based on this more than 2-decade long field study demonstrate the exciting potential of employing microorganisms and maximizing their functions in future agroecosystems. Our results report a "most-wanted" or "most-unwanted" list of microbial phylotypes that are ready candidates to guide the development of sustainable agriculture in alkaline soils.

Microeukaryotic Community Shifting Along a Lentic-Lotic Continuum

As an important regulator of ecosystem functions in river systems, microeukaryotes play an important role in energy and material conversion, yet little is known about the shift along a lentic-lotic continuum. In this study, the 18S rRNA genes sequencing was used to identify the microeukaryotic communities at 82 sites along a lentic-lotic continuum with the aim of understanding the impact of upstream inlet river on microeukaryotic communities in Baiyang Lake (BYD) and its downstream. Our results showed that the upstream inlet river affected the diversity and community composition of microeukaryotes in BYD and downstream rivers, and environmental variables greatly affected the composition of microeukaryotic community. The community composition in BYD had lower variabilities. Co-occurrence network analysis revealed that the network was non-random and clearly parsed into three modules, and different modules were relatively more abundant to a particular area. As keystone taxa, some nodes of the upstream microeukaryotic network played an important role in structuring network and maintaining the stability of the ecosystem. In BYD and downstream, the microeukaryotic network was highly fragmented, and the loss of keystone taxa would have an adverse impact on the integrity and function of the microeukaryotic community. Microeukaryotes had strong tendencies to co-occur, which may contribute to the stability and resilience of microeukaryotic communities. Overall, these findings extend the current understanding of the diversity and community composition of microeukaryotic along a lentic-lotic continuum.

Compositional Changes and Co-Occurrence Patterns of Planktonic Bacteria and Microeukaryotes in a Subtropical Estuarine Ecosystem, the Pearl River Delta

Planktonic microorganisms in aquatic ecosystems form complex assemblages of highly interactive taxa and play key roles in biogeochemical cycles. However, the microbial interactions within bacterial and microeukaryotic communities, and the mechanisms underpinning the responses of abundant and rare microbial taxa to environmental disturbances in the river estuary remain unknown. Here, 16S and 18S rRNA gene sequencing were used to investigate the compositional changes and the co-occurrence patterns of bacterial and microeukaryotic communities. The results showed that the rare taxa in the bacterial communities were more prevalent than those in the microeukaryotic communities and may influence the resilience and resistance of microorganisms to environmental variations in estuarine ecosystems. The environmental variations had strong effects on the microeukaryotic communities and their assembly mechanisms but not on the bacterial communities in our studied area. However, based on co-occurrence network analyses, the bacterial communities had stronger links and more complex interactions than microeukaryotic communities, suggesting that bacterial networks may help improve the buffering capacities of the estuarine ecosystem against environmental change. The keystone taxa of bacteria mainly belonged to rare subcommunities, which further illustrates that rare taxa may play fundamental roles in network persistence. Overall, these results provide insights into the microbial responses of aquatic ecosystems to environmental heterogeneity.

Rare-earth metal oxide nanoparticles decouple the linkage between soil bacterial community structure and function by selectively influencing potential keystone taxa

Excessive production and application of rare-earth metal oxide nanoparticles warrants assessment of their environmental risks. Little is known about the impact of these nanoparticles on soil bacterial communities. We quantified the effects of nano-Gd₂O₃ and nano-La₂O₃, at the different concentrations and exposure regimes, on soil bacterial community structure and function as well as the structure-function relationship. Further, we constructed and analyzed a co-occurrence network to identify and characterize potential keystone taxa that were related to the enzyme activities and responded to the increasing concentrations of nanoparticles. Both nano-Gd₂O₃ and nano-La₂O₃ significantly altered the bacterial community structure and function in a concentration-dependent manner; however, these negative effects were observed on day 1 or day 7 but not on day 60, indicating that these effects were transient and the bacterial communities can mitigate the effect of these nanoparticles over time. Interestingly, the nanoparticle exposure decoupled the relationship between the structure and function of the soil bacterial communities. The decoupling was due to changes in the composition and relative abundances of potential keystone taxa related to bacterial community functions. Altogether, we provide insights into the interactions between the rare-earth metal oxide nanoparticles and soil bacterial communities. Our results facilitate the environmental risk assessment and safe usage of rare-earth metal oxide nanoparticles.

Spatio-temporal patterns of multi-trophic biodiversity and food-web characteristics uncovered across a river catchment using environmental DNA

Accurate characterisation of ecological communities with respect to their biodiversity and food-web structure is essential for conservation. However, combined empirical study of biodiversity and multi-trophic food webs at a large spatial and temporal resolution has been prohibited by the lack of appropriate access to such data from natural systems. Here, we assessed biodiversity and food-web characteristics across a 700 km² riverine network over seasons using environmental DNA. We found contrasting biodiversity patterns between major taxonomic groups. Local richness showed statistically significant, season-dependent increases and decreases towards downstream location within the catchment for fish and bacteria, respectively. Meanwhile, invertebrate richness remained spatially unchanged but varied across seasons. The structure of local food webs, such as link density and nestedness, also varied across space and time. However, these patterns did not necessarily mirror those observed for biodiversity and functional feeding characteristics. Our results suggest that biodiversity patterns and food-web dynamics are not directly scalable to each other even at the same spatial and temporal scales. In order to conserve species diversity as well as the functional trophic integrity of communities, patterns of biodiversity and food-web characteristics must thus be jointly studied.

By sampling environmental DNA across a large riverine network over multiple seasons, the varied dynamics between biodiversity and food-web dynamics are revealed.

Vertical and temporal variations of soil bacterial and archaeal communities in wheat-soybean rotation agroecosystem

Soil microbes are an essential component of terrestrial ecosystems and drive many biogeochemical processes throughout the soil profile. Prior field studies mainly focused on the vertical patterns of soil microbial communities, meaning their temporal dynamics have been largely neglected. In the present study, we investigated the vertical and temporal patterns of soil bacterial and archaeal communities in a wheat-soybean rotation agroecosystem at a depth of millions of sequences per sample. Our results revealed different vertical bacterial and archaeal richness patterns: bacterial richness was lowest in the deep soil layer and peaked in the surface or middle soil layer. In contrast, archaeal richness did not differ among soil layers. PERMANOVA analysis indicated that both bacterial and archaeal community compositions were significantly impacted by soil depth but unaffected by sampling time. Notably, the proportion of rare bacteria gradually decreased along with the soil profile. The rare bacterial community composition was the most important indicator for soil nutrient fertility index, as determined by random forest analysis. The soil prokaryotic co-occurrence networks of the surface and middle soil layers are more connected and harbored fewer negative links than that of the deep soil layer. Overall, our results highlighted soil depth as a more important determinant than temporal variation in shaping the soil prokaryotic community and interspecific interactions and revealed a potential role of rare taxa in soil biogeochemical function.

Environmental effects of nanoparticles on the ecological succession of gut microbiota across zebrafish development

The environmental stresses could significantly affect the structure and functions of microbial communities colonized in the gut ecosystem. However, little is known about how engineered nanoparticles (ENPs), which have recently become a common pollutant in the environment, affect the gut microbiota across fish development. Based on the high-throughput sequencing of the 16S rRNA gene amplicon, we explored the ecological succession of gut microbiota in zebrafish exposed to nanoparticles for three months. The nanoparticles used herein including titanium dioxide nanoparticles (nTiO₂, 100 μg/L), zinc oxide nanoparticles (nZnO, 100 μg/L), and selenium nanoparticles (nSe, 100 μg/L). Our results showed that nanoparticles exposure reduced the alpha diversity of gut microbiota at 73-90 days post-hatching (dph), but showed no significant effects at 14-36 dph. Moreover, nTiO₂ significantly (p < 0.05) altered the composition of the gut microbial communities at 73-90 dph (e.g., decreasing abundance of Cetobacterium and Vibrio). Moreover, we found that homogeneous selection was the major process (16.6-57.8%) governing the community succession of gut microbiota. Also, nanoparticles exposure caused topological alterations to microbial networks and led to increased positive interactions to destabilize the gut microbial community. This study reveals the environmental effects of nanoparticles on the ecological succession of gut microbiota across zebrafish development, which provides novel insights to understand the gut microbial responses to ENPs over the development of aquatic animals.

Bacterioplankton Diversity and Distribution in Relation to Phytoplankton Community Structure in the Ross Sea Surface Waters

Primary productivity in the Ross Sea region is characterized by intense phytoplankton blooms whose temporal and spatial distribution are driven by changes in environmental conditions as well as interactions with the bacterioplankton community. However, the number of studies reporting the simultaneous diversity of the phytoplankton and bacterioplankton in Antarctic waters are limited. Here, we report data on the bacterial diversity in relation to phytoplankton community structure in the surface waters of the Ross Sea during the Austral summer 2017. Our results show partially overlapping bacterioplankton communities between the stations located in the Terra Nova Bay (TNB) coastal waters and the Ross Sea Open Waters (RSOWs), with a dominance of members belonging to the bacterial phyla Bacteroidetes and Proteobacteria. In the TNB coastal area, microbial communities were characterized by a higher abundance of sequences related to heterotrophic bacterial genera such as Polaribacter spp., together with higher phytoplankton biomass and higher relative abundance of diatoms. On the contrary, the phytoplankton biomass in the RSOW were lower, with relatively higher contribution of haptophytes and a higher abundance of sequences related to oligotrophic and mixothrophic bacterial groups like the Oligotrophic Marine Gammaproteobacteria (OMG) group and SAR11. We show that the rate of diversity change between the two locations is influenced by both abiotic (salinity and the nitrogen to phosphorus ratio) and biotic (phytoplankton community structure) factors. Our data provide new insight into the coexistence of the bacterioplankton and phytoplankton in Antarctic waters, suggesting that specific rather than random interaction contribute to the organic matter cycling in the Southern Ocean.

Hydrological properties predict the composition of microbial communities cycling methane and nitrogen in rivers

Sediment microbial communities drive the biogeochemical cycles that make rivers globally important sources and sinks of carbon (C) and nitrogen (N). The structure of these communities is strongly determined by the local physico-chemical environment. However, we currently lack an understanding of the factors that determine microbial community structures at the catchment scale. Here, we show that the contribution of groundwater to total river flow (quantified as base flow index; BFI) predicts the structure and diversity of the different microbial functional groups that cycle N and C across nine UK rivers, spanning a geological BFI gradient from 0.23 (clay sediment) to 0.95 (chalk gravel sediment). Furthermore, the GC-content (percentage of guanine-cytosine bases in a DNA sequence) and codon-usage bias of ammonia monooxygenase DNA sequences, and the hydrophobicity and net-charge of the corresponding amino acid sequences, were all strongly correlated with BFI, likely reflecting physiological adaptations to different riverbed sediment structure along the BFI gradient. Our results offer an opportunity to overcome the "paradox of scales" that has seen microbial ecologists focus on small- rather than large-scale environmental variables, enabling us to scale-up our understanding of microbial biogeochemistry to the catchment and beyond.

Microbial co-occurrence network topological properties link with reactor parameters and reveal importance of low-abundance genera

Operational factors and microbial interactions affect the ecology in anaerobic digestion systems. From 12 lab-scale reactors operated under distinct engineering conditions, bacterial communities were found driven by temperature, while archaeal communities by both temperature and substrate properties. Combining the bacterial and archaeal community clustering patterns led to five sample groups (ambient, mesophilic low-solid-substrate, mesophilic, mesophilic co-digestion and thermophilic) for co-occurrence network analysis. Network topological properties were associated with substrate characteristics and hydrolysis-methanogenesis balance. The hydrolysis efficiency correlated (p < 0.05) with clustering coefficient positively and with normalized betweenness negatively. The influent particulate COD ratio and the relative differential hydrolysis-methanogenesis efficiency (D_efficiency) correlated negatively with the average path length (p < 0.05). Individual genera’s topological properties showed more connector genera in thermophilic network, representing stronger inter-module communication. Individual genera’s normalized degree and betweenness revealed that lower-abundance genera (as low as 0.1%) could perform central hub roles and communication roles, maintaining the stability and functionality of the microbial community.

Spatial Variability in Streambed Microbial Community Structure across Two Watersheds

Both spatial and temporal variability are key attributes of sedimentary microbial communities, and while spatial effects on beta-diversity appear to dominate at larger distances, the character of spatial variability at finer scales remains poorly understood, especially for headwater stream communities. We investigated patterns of microbial community structure (MCS) in biofilms attached to streambed sediments from two watersheds across spatial scales spanning <1 m within a single stream to several hundred kilometers between watersheds. Analyses of phospholipid fatty acid (PLFA) profiles indicated that the variations in MCS were driven by increases in the relative abundance of microeukaryotic photoautotrophs and their contribution to total microbial biomass. Furthermore, streams within watersheds had similar MCS, underscoring watershed-level controls of microbial communities. Moreover, bacterial community structure assayed as either PCR-denaturing gradient gel electrophoresis (PCR-DGGE) fingerprints or PLFA profiles edited to remove microeukaryotes indicated a distinct watershed-level biogeography. No distinct stream order-level distributions were identified, although DGGE analyses clearly indicated that there was greater variability in community structure among 1st-order streams than among 2nd- and 3rd-order streams. Longitudinal gradients in microbial biomass and structure showed that the greatest variations were associated with 1st-order streams within a watershed, and 68% of the variation in total microbial biomass was explained by sediment atomic carbon-to-nitrogen ratio (C:N ratio), percent carbon, sediment surface area, and percent water content. This study confirms a distinct microbial biogeography for headwater stream communities driven by environmental heterogeneity across distant watersheds and suggests that eukaryotic photoautotrophs play a key role in structuring bacterial communities on streambed sediments. IMPORTANCE Microorganisms in streams drive many biogeochemical reactions of global significance, including nutrient cycling and energy flow; yet, the mechanisms responsible for the distribution and composition of streambed microbial communities are not well known. We sampled sediments from multiple streams in two watersheds (Neversink River [New York] and White Clay Creek [WCC; Pennsylvania] watersheds) and measured microbial biomass and total microbial and bacterial community structures using phospholipid and molecular methods. Microbial and bacterial community structures displayed a distinct watershed-level biogeography. The smallest headwater streams within a watershed showed the greatest variation in microbial biomass, and 68% of that variation was explained by the atomic carbon-to-nitrogen ratio (C:N ratio), percent carbon, sediment surface area, and percent water content. Our study revealed a nonrandom distribution of microbial communities in streambeds, and showed that microeukaryotic photoautotrophs, environmental heterogeneity, and geographical distance influence microbial composition and spatial distribution.

Estuarine gradients dictate spatiotemporal variations of microbiome networks in the Chesapeake Bay

BACKGROUND

Annually reoccurring microbial populations with strong spatial and temporal variations have been identified in estuarine environments, especially in those with long residence time such as the Chesapeake Bay (CB). However, it is unclear how microbial taxa cooccurr and how the inter-taxa networks respond to the strong environmental gradients in the estuaries.

RESULTS

Here, we constructed co-occurrence networks on prokaryotic microbial communities in the CB, which included seasonal samples from seven spatial stations along the salinity gradients for three consecutive years. Our results showed that spatiotemporal variations of planktonic microbiomes promoted differentiations of the characteristics and stability of prokaryotic microbial networks in the CB estuary. Prokaryotic microbial networks exhibited a clear seasonal pattern where microbes were more closely connected during warm season compared to the associations during cold season. In addition, microbial networks were more stable in the lower Bay (ocean side) than those in the upper Bay (freshwater side). Multivariate regression tree (MRT) analysis and piecewise structural equation modeling (SEM) indicated that temperature, salinity and total suspended substances along with nutrient availability, particulate carbon and Chl a, affected the distribution and co-occurrence of microbial groups, such as Actinobacteria, Bacteroidetes, Cyanobacteria, Planctomycetes, Proteobacteria, and Verrucomicrobia. Interestingly, compared to the abundant groups (such as SAR11, Saprospiraceae and Actinomarinaceae), the rare taxa including OM60 (NOR5) clade (Gammaproteobacteria), Micrococcales (Actinobacteria), and NS11-12 marine group (Bacteroidetes) contributed greatly to the stability of microbial co-occurrence in the Bay. Modularity and cluster structures of microbial networks varied spatiotemporally, which provided valuable insights into the 'small world' (a group of more interconnected species), network stability, and habitat partitioning/preferences.

CONCLUSION

Our results shed light on how estuarine gradients alter the spatiotemporal variations of prokaryotic microbial networks in the estuarine ecosystem, as well as their adaptability to environmental disturbances and co-occurrence network complexity and stability.

The Metabolic Regimes at the Scale of an Entire Stream Network Unveiled Through Sensor Data and Machine Learning

Streams and rivers form dense networks that drain the terrestrial landscape and are relevant for biodiversity dynamics, ecosystem functioning, and transport and transformation of carbon. Yet, resolving in both space and time gross primary production (GPP), ecosystem respiration (ER) and net ecosystem production (NEP) at the scale of entire stream networks has been elusive so far. Here, combining Random Forest (RF) with time series of sensor data in 12 reach sites, we predicted annual regimes of GPP, ER, and NEP in 292 individual stream reaches and disclosed properties emerging from the network they form. We further predicted available light and thermal regimes for the entire network and expanded the library of stream metabolism predictors. We found that the annual network-scale metabolism was heterotrophic yet with a clear peak of autotrophy in spring. In agreement with the River Continuum Concept, small headwaters and larger downstream reaches contributed 16% and 60%, respectively, to the annual network-scale GPP. Our results suggest that ER rather than GPP drives the metabolic stability at the network scale, which is likely attributable to the buffering function of the streambed for ER, while GPP is more susceptible to flow-induced disturbance and fluctuations in light availability. Furthermore, we found large terrestrial subsidies fueling ER, pointing to an unexpectedly high network-scale level of heterotrophy, otherwise masked by simply considering reach-scale NEP estimations. Our machine learning approach sheds new light on the spatiotemporal dynamics of ecosystem metabolism at the network scale, which is a prerequisite to integrate aquatic and terrestrial carbon cycling at relevant scales.

Terrestrial connectivity, upstream aquatic history and seasonality shape bacterial community assembly within a large boreal aquatic network

During transit from soils to the ocean, microbial communities are modified and re-assembled, generating complex patterns of ecological succession. The potential effect of upstream assembly on downstream microbial community composition is seldom considered within aquatic networks. Here, we reconstructed the microbial succession along a land-freshwater-estuary continuum within La Romaine river watershed in Northeastern Canada. We captured hydrological seasonality and differentiated the total and reactive community by sequencing both 16 S rRNA genes and transcripts. By examining how DNA- and RNA-based assemblages diverge and converge along the continuum, we inferred temporal shifts in the relative importance of assembly processes, with mass effects dominant in spring, and species selection becoming stronger in summer. The location of strongest selection within the network differed between seasons, suggesting that selection hotspots shift depending on hydrological conditions. The unreactive fraction (no/minor RNA contribution) was composed of taxa with diverse potential origins along the whole aquatic network, while the majority of the reactive pool (major RNA contribution) could be traced to soil/soilwater-derived taxa, which were distributed along the entire rank-abundance curve. Overall, our findings highlight the importance of considering upstream history, hydrological seasonality and the reactive microbial fraction to fully understand microbial community assembly on a network scale.

Bacterial Communities Present Distinct Co-occurrence Networks in Sediment and Water of the Thermokarst Lakes in the Yellow River Source Area

Thermokarst lakes are a ubiquitous and important landscape feature in cold regions and are changing tremendously due to the accelerated climate change. In thermokarst lakes, sediment and water are two distinct but highly interconnected habitats, harboring different bacterial communities in terms of taxonomic composition. However, the co-occurrence networks of these bacterial communities remain unclear. Here, we investigate the co-occurrence ecological networks of sediment and water bacterial communities for thermokarst lakes in the Yellow River Source Area on the Qinghai-Tibet Plateau. The results show that the bacterial communities construct distinct co-occurrence networks in sediment and water. The metacommunity network was parsed into four major modules formed by the operational taxonomic units (OTUs) enriched in sediment or water independently, and water-enriched OTUs exhibited much closer interconnections than sediment-enriched OTUs. When considering the sediment and water bacterial networks separately, different topological properties and modular patterns present: the sediment bacterial network was more clustered while the modules less responded to the environmental variables. On the contrary, the water bacterial network was more complex with the OTUs more interconnected and its modules more responded to the environmental variables. Moreover, the results of the structural equation model suggest that, by the influence of environmental variations on individual modules, the water bacterial communities would be more vulnerable under the fact of accelerating climate change. This study provides insights beyond a conventional taxonomic perspective, adding our knowledge of the potential mechanisms structuring bacterial community assembly and improving our prediction of the responses of this fast-changing ecosystem to future climate change.