Coherent dynamics and association networks among lake bacterioplankton taxa

Bacterial composition in Swedish raw drinking water reveals three major interacting ubiquitous metacommunities

BACKGROUND

Surface raw water used as a source for drinking water production is a critical resource, sensitive to contamination. We conducted a study on Swedish raw water sources, aiming to identify mutually co-occurring metacommunities of bacteria, and environmental factors driving such patterns.

METHODS

The water sources were different regarding nutrient composition, water quality, and climate characteristics, and displayed various degrees of anthropogenic impact. Water inlet samples were collected at six drinking water treatment plants over 3 years, totaling 230 samples. The bacterial communities of DNA sequenced samples (n = 175), obtained by 16S metabarcoding, were analyzed using a joint model for taxa abundance.

RESULTS

Two major groups of well-defined metacommunities of microorganisms were identified, in addition to a third, less distinct, and taxonomically more diverse group. These three metacommunities showed various associations to the measured environmental data. Predictions for the well-defined metacommunities revealed differing sets of favored metabolic pathways and life strategies. In one community, taxa with methanogenic metabolism were common, while a second community was dominated by taxa with carbohydrate and lipid-focused metabolism.

CONCLUSION

The identification of ubiquitous persistent co-occurring bacterial metacommunities in freshwater habitats could potentially facilitate microbial source tracking analysis of contamination issues in freshwater sources.

A Preliminary Study of Bacterioplankton Community Structure in the Taiyangshan Wetland in Ningxia and Its Driving Factors

The Taiyangshan Wetland, a valuable wetland resource in the arid zone of central Ningxia, is critical for flood storage and drought resistance, climate regulation, and biodiversity protection. Nevertheless, the community structure and diversity of bacterioplankton in the Taiyangshan Wetland remains unclear. High-throughput sequencing was used to analyze the differences in bacterioplankton structure and major determinants in the Taiyangshan Wetland from April to October 2020. The composition and diversity of the bacterioplankton community varied significantly in different sampling periods but showed negligible differences across lake regions. Meanwhile, the relative abundances of bacterioplankton Bacteroidetes, Actinobacteria, Firmicutes, Chloroflexi, Tenericutes, Epsilonbacteraeota, and Patescibacteria were significantly different in different sampling periods, while the relative abundances of Cyanobacteria in different lake regions were quite different. Network analysis revealed that the topological attributes of co-occurrence pattern networks of bacterioplankton were high, and bacterioplankton community compositions were complicated in the month of July. A mantel test revealed that the bacterioplankton community in the entire wetland was affected by water temperature, electrical conductivity, dissolved oxygen, salinity, total nitrogen, ammonia nitrogen, chemical oxygen demand, fluoride, and sulfate. The bacterioplankton community structure was affected by ten environmental parameters (e.g., water temperature, dissolved oxygen, salinity, and permanganate index) in April, while the bacterioplankton community was only related to 1~2 environmental parameters in July and October. The bacterioplankton community structure in Lake Region IV was related to seven environmental parameters, including dissolved oxygen, pH, total nitrogen, and chemical oxygen demand, whereas the bacterioplankton community structures in the other three lake regions were related to two environmental parameters. This study facilitates the understanding of the bacterioplankton community in wetlands in arid areas and provides references to the evaluation of aquatic ecological management of the Taiyangshan Wetland.

Bacterioplankton seasonality in deep high-mountain lakes

Due to global warming, shorter ice cover duration might drastically affect the ecology of lakes currently undergoing seasonal surface freezing. High-mountain lakes show snow-rich ice covers that determine contrasting conditions between ice-off and ice-on periods. We characterized the bacterioplankton seasonality in a deep high-mountain lake ice-covered for half a year. The lake shows a rich core bacterioplankton community consisting of three components: (i) an assemblage stable throughout the year, dominated by Actinobacteria, resistant to all environmental conditions; (ii) an ice-on-resilient assemblage dominating during the ice-covered period, which is more diverse than the other components and includes a high abundance of Verrucomicrobia; the deep hypolimnion constitutes a refuge for many of the typical under-ice taxa, many of which recover quickly during autumn mixing; and (iii) an ice-off-resilient assemblage, which members peak in summer in epilimnetic waters when the rest decline, characterized by a dominance of Flavobacterium, and Limnohabitans. The rich core community and low random elements compared to other relatively small cold lakes can be attributed to its simple hydrological network in a poorly-vegetated catchment, the long water-residence time (ca. 4 years), and the long ice-cover duration; features common to many headwater deep high-mountain lakes.

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

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

Data-driven causal analysis of observational biological time series

Complex systems are challenging to understand, especially when they defy manipulative experiments for practical or ethical reasons. Several fields have developed parallel approaches to infer causal relations from observational time series. Yet, these methods are easy to misunderstand and often controversial. Here, we provide an accessible and critical review of three statistical causal discovery approaches (pairwise correlation, Granger causality, and state space reconstruction), using examples inspired by ecological processes. For each approach, we ask what it tests for, what causal statement it might imply, and when it could lead us astray. We devise new ways of visualizing key concepts, describe some novel pathologies of existing methods, and point out how so-called 'model-free' causality tests are not assumption-free. We hope that our synthesis will facilitate thoughtful application of methods, promote communication across different fields, and encourage explicit statements of assumptions. A video walkthrough is available (Video 1 or https://youtu.be/AIV0ttQrjK8).

Identification of key water parameters and microbiological compositions triggering intensive N2O emissions during landfill leachate treatment process

Landfill leachate treatment processes tend to emit more N₂O compared to domestic wastewater treatment. This discrepancy may be ascribed to leachate water characteristics such as high refractory COD, ammonium (NH₄⁺) content, and salinity. In this work, the leachate influent was varied to examine the N₂O emission scenarios. NH₄⁺-N, COD, and Cl^- concentrations ranged between 1000-2500, 1000-10,000, and 500-3000 mg L^-1, respectively. Simultaneously, we attempted to combine statistical analysis with high-throughput sequencing to understand the microbial mechanism with regards to N₂O emission. Results show that the strong N₂O emissions occur in the nitrifying tank due to the intensive aeration. The system receiving the lowest COD shows the maximum N₂O emission factor of 42.7% of the removed nitrogen. Both redundancy analysis and a structural equation model verify that insufficient degradable organics are the key water parameter triggering intensive N₂O emission within the designed influent limits. Furthermore, two essential but non-abundant functional bacteria, Flavobacterium (acting as a denitrifier) and Nitrosomonas (acting as a nitrifier), are identified as the core functional species that dramatically influence N₂O emissions. An increase in influent COD promotes the proliferation of Flavobacterium and inhibits Nitrosomonas, which in turn reduce N₂O release. Meanwhile, two keystone species of Castellaniella and Saprospiraceae unclassified are recognized. They may supply a suitable niche and integrity of the microbial community for N-cycle functional bacteria. These findings reveal the essential role of non-abundant species in microbial community, and expand the current understanding of microbial interactions underlying N₂O dynamics in leachate treatment systems.

Seasonal Succession and Coherence Among Bacteria and Microeukaryotes in Lake Baikal

Microorganisms exhibit seasonal succession governed by physicochemical factors and interspecies interactions, yet drivers of this process in different environments remain to be determined. We used high-throughput sequencing of 16S rRNA and 18S rRNA genes to study seasonal dynamics of bacterial and microeukaryotic communities at pelagic site of Lake Baikal from spring (under-ice, mixing) to autumn (direct stratification). The microbial community was subdivided into distinctive coherent clusters of operational taxonomic units (OTUs). Individual OTUs were consistently replaced during different seasonal events. The coherent clusters change their contribution to the microbial community depending on season. Changes of temperature, concentrations of silicon, and nitrates are the key factors affected the structure of microbial communities. Functional prediction revealed that some bacterial or eukaryotic taxa that switched with seasons had similar functional properties, which demonstrate their functional redundancy. We have also detected specific functional properties in different coherent clusters of bacteria or microeukaryotes, which can indicate their ability to adapt to seasonal changes of environment. Our results revealed a relationship between seasonal succession, coherency, and functional features of freshwater bacteria and microeukaryotes.

Microbial interactions enhanced environmental fitness and expanded ecological niches under dibutyl phthalate and cadmium co-contamination

Co-contamination of organic pollutants and heavy metals is universal in the natural environment. Dibutyl phthalate (DBP), a typical plasticizer, frequently coexists with cadmium (Cd) in nature. However, little attention has been given to the impacts of co-contamination by DBP and Cd on microbial communities or the responses of microbes. To address this, a microcosm experiment was conducted by supplying the exogenous DBP-degrading bacterium Glutamicibacter nicotianae ZM05 to investigate the interplay among DBP-Cd co-contamination, the exogenous DBP-degrading bacterium G. nicotianae ZM05, and indigenous microorganisms. To adapt to co-contamination stress, microbial communities adjust their diversity, interactions, and functions. The stability of the microbial community decreased under co-contamination, as evidenced by lower diversity, simpler network, and fewer ecological niches. Microbial interactions were strengthened, as evidenced by enriched pathways related to microbial communications. Meanwhile, interactions between microorganisms enhanced the environmental fitness of the exogenous DBP-degrading bacterium ZM05. Based on co-occurrence network prediction and coculture experiments, metabolic interactions between the non-DBP-degrading bacterium Cupriavidus metallidurans ZM16 and ZM05 were proven. Strain ZM16 utilized protocatechuic acid, a DBP downstream metabolite, to relieve acid inhibition and adsorbed Cd to relieve toxic stress. These findings help to explain the responses of bacterial and fungal communities to DBP-Cd co-contamination and provide new insights for the construction of degrading consortia for bioremediation.

Bio-fertilizer Amendment Alleviates the Replanting Disease under Consecutive Monoculture Regimes by Reshaping Leaf and Root Microbiome

Replanting disease is a growing problem in intensive agricultural systems. Application of bio-fertilizer containing beneficial microbes contributes to disease suppression and is a promising strategy to control replanting disease. However, the effect of both replanting disease and bio-fertilizer amendment on the assembly of crop microbiota in leaves and roots and their relationships to crop yield and quality remains elusive. In these experiments, roots and leaves of Radix pseudostellariae were collected from different consecutive monoculture and bio-fertilizer amended fields, and the associated microbiota were characterized by bacterial 16S rRNA gene sequencing and quantitative PCR. Consecutive monoculture altered the bacterial community structure and composition and significantly increased the abundance of potential pathogenic Ralstonia and Fusarium oxysporum in leaves and roots. Furthermore, bio-fertilizer application alleviated replanting disease by decreasing the pathogen load, increasing the potential beneficial genera Pseudomonas, Streptomyces, Paenibacillus, and Bradyrhizobium. The proportion of positive correlations in the co-occurrence network of bio-fertilizer application was the highest, implying that bio-fertilizer potentially enhanced ecological commensalism or mutualism of the bacterial community across the two compartments. Structural equation models indicated that bio-fertilizer had a positive and indirect effect on both yield and quality by shaping the leaf microbiota and the root microbiota. Our findings highlight the role of leaf and root microbiota on replanting disease, showing that bio-fertilizer contributes to alleviating replanting disease by improving microbe-microbe interactions.

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

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

Host Species Influence the Gut Microbiota of Endemic Cold-Water Fish in Upper Yangtze River

The fish gut microbiome plays an important role in nutrition absorption and energy metabolism. Studying the gut microbes of cold-water fish is important to understand the dietary adaptation strategies in extreme environments. In this study, the gut samples of Schizothorax wangchiachii (SW, herbivorous), Schizothorax kozlovi (SK, omnivorous), and Percocypris pingi (PP, carnivorous) in the upper Yangtze River were collected, and we sequenced 16S rRNA amplicon to study the potential relationship between gut microbes and host species. The results showed that gut microbial composition and diversity were significantly different between the three cold-water fishes. These fishes had different key taxa in their gut microbes, including bacteria involved in the breakdown of food (e.g., Cetobacterium, Aeromonas, and Clostridium sensu stricto 10). The highest alpha diversity indices (e.g., Chao 1 index) were identified in the herbivore (SW), followed by the carnivore (PP), and the lowest in the omnivore (SK). Non-metric multidimensional scaling (NMDS) results revealed that the gut microbial community of these species was different between host species. The neutral community model (NCM) showed that the microbial community structure of SW was shaped by stochastic processes, and the highest species dispersal was found in SW, followed by PP, and the lowest in SK. The results of niche breadth agreed with these findings. Our results demonstrated that host species influenced the gut microbiome composition, diversity, and microbial community assembly processes of the three cold-water fishes. These findings implied that the variation of gut microbiome composition and function plays a key role in digesting and absorbing nutrients from different foods in cold-water fish.

Time-lagged interspecies interactions prevail during biofilm development in moving bed biofilm reactor

Clarifying the essential succession dynamics of interspecies interactions during biofilm development is crucial for the regulation and application of biofilm-based processes. In this study, regular and time-series phylogenetic molecular ecological networks (pMENs) were constructed to investigate ordinary and time-lagged interspecies interactions during biofilm development in a moving bed biofilm reactor (MBBR). Positive interactions dominated both regular (89.78%) and time-series (77.04%) ecological networks, suggesting that extensive cooperative behaviors facilitated biofilm development. The pronounced directional interactions (72.52%) in the time-series network further indicated that time-lagged interspecies interactions prevailed in the biofilm development process. Specifically, the proportion of directional negative interactions was higher than that of positive interactions, implying that interspecific competition preferred to be time-lagged. The time-series network revealed that module hubs exhibited extensive time-lagged positive interactions with their neighbors, and most of them exhibited altruistic behaviors. Keystone species possessing more positive interactions were positively correlated with biofilm biomass, NO₃ ^- -N concentrations, and the removal efficiencies of NH₄ ⁺ -N and COD. However, keystone species and peripherals that were negatively targeted by their neighbors showed positive correlations with the concentrations of NO₂ ^- -N, polysaccharides, and proteins in the soluble microbial products. The data highlight that the time-series network can provide directional microbial interactions along with the biofilm development process, which would help to predict the tendency of community shifts and propose efficient strategies for the regulation of biofilm-based processes. This article is protected by copyright. All rights reserved.

High stability of autochthonous dissolved organic matter in karst aquatic ecosystems: Evidence from fluorescence

Biological carbon pump (BCP) in karst areas has received intensive attention for years due to their significant contribution to the global missing carbon sink. The stability of autochthonous dissolved organic matter (Auto-DOM) produced by BCP in karst aquatic ecosystems may play a critical role in the missing carbon sink. However, the source of dissolved organic matter (DOM) in inland waters and its consumption by planktonic bacteria have not been thoroughly examined. Recalcitrant dissolved organic matter (RDOM) may exist in karst aquatic ecosystem as in the ocean. Through the study of the chromophoric dissolved organic matter (CDOM) and the interaction between CDOM and the planktonic bacterial community under different land uses at the Shawan Karst Water-carbon Cycle Test Site, SW China, we found that C2, as the fluorescence component of Auto-DOM mineralised by planktonic bacteria, may have some of the characteristics of RDOM and is an important DOM source in karst aquatic ecosystems. The stability ratio (Fmax_(C2/(C1+C2))) of Auto-DOM reached 89.6 ± 6.71% in winter and 64.1 ± 7.19% in spring. Moreover, correlation-based network analysis determined that the planktonic bacterial communities were controlled by different fluorescence types of CDOM, of which C1 (fresh Auto-DOM), C3 (conventional allochthonous DOM (Allo-DOM)) and C4 (the Allo-DOM mineralised by bacteria) were clustered in one module together with prevalent organic-degrading planktonic bacteria; C2 was clustered in another tightly combined module, suggesting specific microbial utilization strategies for the C2 component. In addition, some important planktonic bacterium and functional genes (including chemotrophic heterotrophs and photosynthetic bacteria) were found to be affected by high Ca²⁺ and dissolved inorganic carbon (DIC) concentrations in karst aquatic ecosystems. Our research showed that Auto-DOM may be as an important carbon sink as the Allo-DOM in karst ecosystems, the former generally being neglected based on a posit that it is easily and first mineralized by planktonic bacteria.

Seagrasses and local environment control the bacterial community structure and carbon substrate utilization in brackish sediments

Seagrasses are complex benthic coastal ecosystems that play a crucial role in organic matter cycling and carbon sequestration. However, little is known about how seagrasses influence the structure and carbon utilization potential of benthic bacterial communities. This study examined the bacterial communities in monospecific and mixed meadows of seagrasses and compared with bulk (unvegetated) sediments from Chilika, a brackish water coastal lagoon of India. High-throughput sequencing of 16S rRNA genes revealed a vegetation effect in terms of differences in benthic bacterial community diversity, composition, and abundances in comparison with bulk sediments. Desulfobacterales, Chromatiales, Enterobacteriales, Clostridiales, Vibrionales, and Acidimicrobiales were major taxa that contributed to differences between seagrass and bulk sediments. Seagrasses supported ∼5.94 fold higher bacterial abundances than the bulk due to rich organic carbon stock in their sediments. Co-occurrence network demonstrated much stronger potential interactions and connectedness in seagrass bacterial communities compared to bulk. Chromatiales and Acidimicrobiales were identified as the top two keystone taxa in seagrass bacterial communities, whereas, Dehalococcoidales and Rhizobiales were in bulk communities. Seagrasses and local environmental factors, namely, water depth, water pH, sediment salinity, redox potential, total organic carbon, available nitrogen, sediment texture, sediment pH, and sediment core depth were the major drivers of benthic bacterial community composition. Carbon metabolic profiling revealed that heterotrophic bacteria in seagrass sediments were much more metabolically diverse and active than bulk. The utilization of carbon substrate guilds, namely, amino acids, amines, carboxylic acids, carbohydrates, polymers, and phenolic compounds was enhanced in seagrass sediments. Metabolic mapping predicted higher prevalence of sulfate-reducer and N₂ fixation metabolic functions in seagrass sediments. Overall, this study showed that seagrasses control benthic bacterial community composition and diversity, enhance heterotrophic carbon substrate utilization, and play crucial roles in organic matter cycling including degradation of hydrocarbon and xenobiotics in coastal sediments.

Animal Activities of the Key Herbivore Plateau Pika (Ochotona curzoniae) on the Qinghai-Tibetan Plateau Affect Grassland Microbial Networks and Ecosystem Functions

Plateau pikas (Ochotona curzoniae) are high-altitude model animals and famous "ecosystem engineers" on the Qinghai-Tibet Plateau. Pika activities may accelerate the degradation of alpine meadows. Nevertheless, little is known about the responses of bacterial, fungal, and archaeal communities, and ecosystem multifunctionality to pika perturbations. To address this question, we studied the impacts of only pika disturbance and combined disturbance (pika disturbance and grazing) on ecological networks of soil microbial communities and ecosystem multifunctionality. Our results demonstrated that Proteobacteria, Ascomycota, and Crenarchaeota were dominant in bacteria, fungi, and archaea, respectively. Bacteria, fungi, and archaea were all influenced by the combined disturbance of grazing and pika. Most fungal communities became convergent, while bacterial and archaeal communities became differentiated during the succession of surface types. In particular, the bacterial and fungal networks were less stable than archaeal networks. In response to the interference, cross-domain cooperation between bacterial and fungal communities increased, while competitive interactions between bacterial and archaeal communities increased. Pika disturbance at high intensity significantly reduced the ecosystem multifunctionality. However, the mixed effects of grazing and pika weakened such influences. This study revealed how pika activities affected microbial networks and ecosystem multifunctionality. These results provide insights to designing reasonable ecological management strategies for alpine grassland ecosystems.

Biological Microbial Interactions from Cooccurrence Networks in a High Mountain Lacustrine District

A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed coexisting. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, where biotic interactions might be required to make the most of an extreme environment. We studied a high-throughput gene data set of alpine lakes (>220 Pyrenean lakes) with cooccurrence network analysis to infer potential biotic interactions, using the combination of a probabilistic method for determining significant cooccurrences and coexclusions between pairs of species and a conceptual framework for classifying the nature of the observed cooccurrences and coexclusions. This computational approach (i) determined and quantified the importance of environmental variables and spatial distribution and (ii) defined potential interacting microbial assemblages. We determined the properties and relationships between these assemblages by examining node properties at the taxonomic level, indicating associations with their potential habitat sources (i.e., aquatic versus terrestrial) and their functional strategies (i.e., parasitic versus mixotrophic). Environmental variables explained fewer pairs in bacteria than in microbial eukaryotes for the alpine data set, with pH alone explaining the highest proportion of bacterial pairs. Nutrient composition was also relevant for explaining association pairs, particularly in microeukaryotes. We identified a reduced subset of pairs with the highest probability of species interactions (“interacting guilds”) that significantly reached higher occupancies and lower mean relative abundances in agreement with the carrying capacity hypothesis. The interacting bacterial guilds could be more related to habitat and microdispersal processes (i.e., aquatic versus soil microbes), whereas for microeukaryotes trophic roles (osmotrophs, mixotrophs, and parasitics) could potentially play a major role. Overall, our approach may add helpful information to guide further efforts for a mechanistic understanding of microbial interactions in situ.

IMPORTANCE A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed to coexist. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, in which biotic interactions might be required to make the most of an extreme environment. Microbial metacommunities are too often only studied in terms of their environmental niches and geographic barriers since they show inherent difficulties to quantify biological interactions and their role as drivers of ecosystem functioning. Our study highlights that telling apart potential interactions from both environmental and geographic niches may help for the initial characterization of organisms with similar ecologies in a large scope of ecosystems, even when information about actual interactions is partial and limited. The multilayered statistical approach carried out here offers the possibility of going beyond taxonomy to understand microbiological behavior in situ.

Responses of mesophilic anaerobic sludge microbiota to thermophilic conditions: Implications for start-up and operation of thermophilic THP-AD systems

Anaerobic digestion (AD) has been widely employed for wastewater and organic waste treatment, in which methanogenesis is highly driven by close microbial interactions among intricate microbial communities. However, the ecological processes underpinning the community assembly that support methanogenesis in such engineered ecosystems remain largely unknown, especially when exposed to challenging circumstances (e.g., high temperature, ammonium content). Here, eight AD bioreactors were seeded with four different inocula (two from full-scale mesophilic AD systems and the other two from lab-scale mesophilic AD systems), and were operated under thermophilic conditions (55 °C) for treating thermal hydrolysis process (THP) pre-treated waste activated sludge to investigate how mesophilic community responds to thermophilic conditions during the long-term cultivation. Results showed that the inocula collected from the full-scale systems were more resilient than that from the lab-scale systems, which may be primarily attributed to indigenous robust methanogens. As a result, the former efficiently generated methane which was predominantly contributed by Methanothermobacter and Methanosarcina (healthy AD ecosystem), while methanogenic activity was remarkably prohibited in the latter (dysfunctional AD ecosystem). Thermophilic environment was a strong selection force, resulting in the convergence of microbial communities in both the healthy and dysfunctional AD ecosystems. Deterministic processes predominated the community assembly regardless of AD ecosystem function, but stronger influences of stochastic processes were observed in dysfunctional AD ecosystems, which was likely attributable from the stronger effect of immigrants from the feedstock. As indicated by molecular ecological network analysis, the microbial network structures in the healthy AD ecosystems were more stable than those in the dysfunctional AD ecosystems. Although keystone taxa were different among the bioreactors, most of which played vital roles in organic hydrolysis/fermentation. To sum up, this study greatly improved our understanding of the relationships between microbiological traits and AD ecosystem function under thermophilic conditions, which could provide useful information to guide thermophilic AD (e.g., THP-AD) start-up and health diagnosis during operation.

Per-, poly-fluoroalkyl substances (PFASs) and planktonic microbiomes: Identification of biotic and abiotic regulations in community coalescence and food webs

The importance of per-, poly-fluoroalkyl substances (PFASs) effects on riverine microbiomes is receiving increased recognition in the environmental sciences. However, few studies have explored how PFASs affect microbiomes across trophic levels, specifically through predator-prey interactions. This study examined the community profiles of planktonic archaea, bacteria, fungi, algae, protozoa, and metazoa in a semi-industrial and agricultural river alongside their interactions with 15 detected PFASs. As abiotic factors, PFASs affected community coalescence more than biogenic substances (p < 0.05). For biotic regulations, sub-communities in rare biospheres (including always rare taxa-ART and critically rare taxa-CRT) contributed to spatial community coalescence more than sub-communities in abundant biospheres (always abundant taxa-AAT and critically abundant taxa-CAT) (p < 0.05). Metazoa-bacteria (Modularity = 1.971) and protozoa-fungi (1.723) were determined to be the most stable predator-prey networks. Based on pathway models, short-chain PFBA (C4) was shown to weaken the trophic transfer efficiencies from heterotrophic bacteria (HB) to heterotrophic flagellates (HF) (p < 0.05). Long-chain PFTeDA (C14) promoted HB to amoeba (p < 0.05), which we postulate is the pathway for PFTeDA to enter the microbial food chain. Our preliminary results elucidated the influence of PFASs on planktonic microbial food webs and highlighted the need to consider protecting and remediating riverine ecosystems containing PFASs.

The counteraction of anammox community to long-term nitrite stress: Crucial roles of rare subcommunity

Understanding the temporal dynamics and recovery of anammox community under nitrite stress is critical for successful application of anammox-related processes. Here, the response behaviors of anammox community were investigated to characterize the reactor performance and ecological function under varied levels of nitrite stress (changing from 0, 50, 100, 200 to 0 mg-N/L) across a large temporal scale (588 days). The nitrogen removal rates decreased from 0.51 ± 0.02 to 0.16 ± 0.04 kg-N/(m³·d) under nitrite stress from 0 to 200 mg-N/L, while it was recovered to 0.29 ± 0.06 kg-N/(m³·d) as nitrite stress terminated. A strong community succession was driven by the initial nitrite stress of 50 mg-N/L, while the community dissimilarity mainly resulted from the increased beta diversity of rare subcommunity. Meanwhile, the rare subcommunity with high functional redundancy likely warranted the functional resilience of anammox community across the nitrite stress gradients. Moreover, the increased positive interactions between anammox bacteria and side populations supported the resilience of anammox after discontinuing nitrite stress, which facilitated the recovery of nitrogen removal efficiency. This study deciphers the interspecies interactions and functional redundancy of rare subcommunity in shaping the robustness and resilience of anammox-related processes when treating nitrite fluctuated wastewater.

Properties of sediment dissolved organic matter respond to eutrophication and interact with bacterial communities in a plateau lake

Sediment dissolved organic matter (DOM) in inland waters is commonly affected by environmental changes. However, knowledge about how sediment DOM responds to eutrophication and the associations between sediment DOM and bacterial communities requires further investigation. We selected a sediment core from Dianchi Lake (China) that was dated from 1864 to 2019 by the activity of radionuclides (²¹⁰Pb and ¹³⁷Cs). δ¹³C_DOC changes fit well with the historical record that heavy eutrophic status in Dianchi Lake were observed since 1980s. Large amounts of dissolved organic carbon (DOC), chromophoric (CDOM) and fluorescent (FDOM) DOM accumulated at the top of the sediments during the eutrophication period (1982-present). The additional algae sources with a higher degradation rate altered the composition, aromaticity and humification of DOM. After long-term mineralization, the remaining DOM became more and more recalcitrant and kept a relatively stable level at older sediments. A co-occurrence network analysis revealed that Proteobacteria, Chloroflexi, Acidobacteriota, Bacteroidota and Desulfobacterota were the most abundant species at the phylum level and clustered into three primary modules. Different microbes shared unique preferences for niches, causing a heterogeneous bacterial distribution at different depths. We conducted Spearman's correlation and redundancy analysis (RDA) to explore potential interactions between bacterial community and sediment DOM. The richness and diversity of bacterial communities were positively related to DOM content, suggesting abundant DOM can produce more available resources for bacteria. RDA results showed some specific species might modify DOM composition and structure. This study suggests that sediment DOM properties were regulated by source transformation during eutrophication, and emphasizes the importance of microbial role on sediment biogeochemical process.

Bacterioplankton community assembly in migratory fish habitat: a case study of the southern East China Sea

The health conditions of fish habitats, which affect fish health, can be reflected by the structure and ecological function of the bacterioplankton community to some extent. However, studies on bacterioplankton in the whole habitat of migratory fish, which can be divided into different functional types, are still limited. To fill this gap, we investigated the characteristics of bacterioplankton communities in three habitat types in a typical migratory fish habitat, the southern East China Sea, using 16S rRNA gene amplicon sequencing. Our study showed that the structure of the bacterioplankton community was significantly divided according to habitat type. Dispersal limitation and heterogeneous selection both contributed to the bacterioplankton community assembly through estimation of β nearest taxon index (βNTI), and redundancy analysis (RDA) further explained that the water temperature, salinity, and nutrients were deterministic factors responsible for differences in the bacterioplankton community. Additionally, different ecological functional modules dominated by functional bacterioplankton in different habitat types were identified by co-occurrence network analysis, including a hydrocarbon-degrading module dominated by Psychrobacter and health-related modules containing Ascidiaceihabitans and Pseudoalteromonas. Based on the composition of environmental bacterioplankton, our findings provide a theoretical basis for understanding the distribution of different habitat types in the southern East China Sea during the breeding period of migratory fish.

Rare soil bacteria are more responsive in desertification restoration than abundant bacteria

Soil microbes play key roles in ecosystem functions, especially in the recovery of ecosystems from disturbance, and exploring community assembly under changing environments has long been a central theme in microbial ecology. The response of abundant and rare bacteria in desertified land to restoration is still unclear. Here, we investigated the effects of vegetation restoration on the assemblage patterns of abundant and rare bacteria in soil across the four sandy lands (Hulunbeir, Horqin, Otindag, and Mu Us) in northern China. Our results revealed that abundant bacteria maintained a relatively stable state under restoration, whereas rare taxa were more responsive, indicating the higher resilience of the rare community to change. Our network analysis also showed that restoration promoted destabilizing properties in rare, but not in abundant, bacterial co-occurrence networks in soil. Environmental selection played a key role in abundant and rare community assembly under restoration. Of the two, the rare subcommunity was mainly affected by environmental filtering. The variations in the abundant and rare communities at the sampling sites under restoration were controlled mainly by plant species richness, and stronger effects were observed in the rare taxa. Overall, these results provide new insight into the mechanisms controlling bacterial community assembly in response to vegetation restoration.

Time after time: Detecting annual patterns in stream bacterial biofilm communities

To quantify the major environmental drivers of stream bacterial population dynamics, we modelled temporal differences in stream bacterial communities to quantify community shifts, including those relating to cyclical seasonal variation and more sporadic bloom events. We applied Illumina MiSeq 16S rRNA bacterial gene sequencing of 892 stream biofilm samples, collected monthly for 36‐months from six streams. The streams were located a maximum of 118 km apart and drained three different catchment types (forest, urban and rural land uses). We identified repeatable seasonal patterns among bacterial taxa, allowing their separation into three ecological groupings, those following linear, bloom/trough and repeated, seasonal trends. Various physicochemical parameters (light, water and air temperature, pH, dissolved oxygen, nutrients) were linked to temporal community changes. Our models indicate that bloom events and seasonal episodes modify biofilm bacterial populations, suggesting that distinct microbial taxa thrive during these events including non‐cyanobacterial community members. These models could aid in determining how temporal environmental changes affect community assembly and guide the selection of appropriate statistical models to capture future community responses to environmental change.

Can we shape microbial communities to enhance biological activated carbon filter performance?

A new focus on biofiltration has emerged that aims to shape microbial communities to improve treatment efficacy. It is therefore necessary to understand the linkages between microbial community structure and biofilter function. However, the assembly and interaction of microbial communities in biological activated carbon (BAC) filters are unknown. In this study, we selected one coal-based granular activated carbon (GAC), GAC-13, with simultaneously developed micropore and micro-level macropore volume used for a bench-scale BAC column experiment, and compared it with other coal-based GACs and wood-based GAC in terms of the dissolved organic carbon (DOC) removal and microbial community characteristics. The results showed that there was no difference between the DOC removal efficiency of BAC-13 and the other two coal-based BAC filters with high iodine value in the period dominated by adsorption, while the DOC removal efficiency of BAC-13 (64.7±0.6%) was significantly higher than that of other BAC filters (36.3±0.8-54.1±0.4%) with a difference of 0.3-0.7 mg/L in DOC during the steady state. The bacterial communities were strongly assembled by deterministic rather than stochastic factors, where the surface polarity of GAC had a greater effect on the microbial communities than its physical properties. The corresponding co-occurrence network revealed that microbes in the BAC filter may be more cooperative than competitive. The keystone bacterium Hyphomicrobium, which had a relatively low abundance, contributed 0.3-1% more to the most abundant functions and produced 5-21 proteins/(g·GAC) more than the dominant bacterium Sphingobium. The metaproteomic-based approach could provide more accurate information regarding the contributions of different species to metabolic functions. The pore size distribution of GAC was found to be an important factor in determining BAC filter performance; the most important pore sizes were micropores and micro-level macropores (0.2-10 μm and >100 μm in diameter), and the latter impacted the abundance of keystone species. Overall, our findings provide new insights into shaping microbial communities by optimizing pore size structure to improve BAC performance, especially the abundance of keystone species.

Network-based analysis reveals differences in plant assembly between the native and the invaded ranges

Associated with the introduction of alien species in a new area, interactions with other native species within the recipient community occur, reshaping the original community and resulting in a unique assemblage. Yet, the differences in community assemblage between native and invaded ranges remain unclear. Mediterranean grasslands provide an excellent scenario to study community assembly following transcontinental naturalisation of plant species. Here, we compared the community resemblance of plant communities in Mediterranean grasslands from both the native (Spain) and invaded (Chile) ranges. We used a novel approach, based on network analysis applied to co-occurrence analysis in plant communities, allowing us to study the co-existence of native and alien species in central Chile. This useful methodology is presented as a step forward in invasion ecology studies and conservation strategies. We found that community structure differed between the native and the invaded range, with alien species displaying a higher number of connections and, therefore, acting as keystones to sustain the structure within the invaded community. Alien species acting like keystones within the Chilean grassland communities might exacerbate the threat posed by biological invasions for the native biodiversity assets. Controlling the spread of the alien species identified here as keystones should help managing potential invasion in surrounding areas. Network analyses is a free, easy-to-implement and straightforward visual tool that can be widely used to reveal shifts in native communities and elucidate the role of multiple invaders into communities.

Varied interspecies interactions between anammox and denitrifying bacteria enhanced nitrogen removal in a single-stage simultaneous anammox and denitrification system

The simultaneous anammox and denitrification (SAD) system has received growing interest for the enhanced nitrogen removal, while the ecological traits of microbial community including spatial distribution characteristics, assembly processes and interspecies interactions have not been fully unraveled. The present study applied metagenomics and ecological analysis methods to gain the ecological traits of microbial communities in the SAD system across different organic substrate loadings. Results showed that organic matter significantly affected the bioreactor performance, and the optimal total nitrogen removal efficiency reached 93.4 ± 0.7% under the COD concentrations of 180 ± 18.2 mg/L. Functional organisms including Candidatus Brocadia (3.9%), Denitratisoma (1.6%), Dokdonella (4.4%) and Thauera (4.6%) obviously enriched under the optimal organic loading conditions. Moreover, microbial communities were significantly governed by deterministic process under high organic concentrations, and the denitrifying organisms displayed important ecological roles in the communities. Although anammox bacteria obviously enriched at the middle of bioreactor, it possessed the highest expression activities at both bottom and middle sites. Denitrifying bacteria that enriched at the bottom sites strongly achieved nitrate reduction and provided nitrite for anammox bacteria, while these organisms trended to compete nitrite with anammox bacteria at the middle site. These findings highlight the importance of microbial ecology in the SAD systems, which may expand our understanding of the synergistic patterns between anammox and denitrifying bacteria.

Soil microbial co-occurrence networks become less connected with soil development in a high Arctic glacier foreland succession

Classically, ecologists have considered that biota becomes more integrated and interdependent with ecosystem development in primary successional environments. However, recent work on soil microbial communities suggests that there may in fact be no change in network integration over successional time series. Here, we performed a test of this principle by identifying network-level topological features of the soil microbial co-occurrence networks in the primary successional foreland environment of the retreating high-Arctic glacier of Midtre Lovénbreen, Svalbard. Soil was sampled at sites along the foreland of inferred ages 10-90 years since deglaciation. DNA was extracted and amplicon sequenced for 16 s rRNA genes for bacteria and ITS1 region for fungi. Despite the chronologically-related soil pH decline and organic C/N accumulation, analysis on network-level topological features showed network integration did not change with inferred chronological ages, whereas network integration declined with decreasing pH and increasing total organic carbon (TOC) - both factors that can be viewed as an indicator of soil development. We also found that bacteria played a greater role in the network structure than fungi, with all keystone species in the microbial co-occurrence network being bacteria species. Both number and relative abundance of the keystone species were significantly higher when soil pH increased or TOC decreased. It appears that in the more extreme and less productive conditions of early primary succession, integration between members of soil biota into consortia may play a greater role in niche adaptation and survival. Our finding also emphasizes that ecosystem development is not simply a product of time but is influenced by locally heterogeneous factors.

Pseudomonas sp. TCd-1 significantly alters the rhizosphere bacterial community of rice in Cd contaminated paddy field

Cadmium (Cd) pollution of paddy soils is one of the main concerns causing food security and environmental problems. Microbial bioremediation is an effective and eco-friendly measure that uses microbes to reduce Cd accumulation in crops. Additionally, rhizosphere bacterial communities also act essential roles in crop tolerance of heavy metals. However, the effects of inoculations with Cd resistant bacteria on crop rhizosphere bacterial communities under Cd exposure are largely unknown. In this study, we used high-throughput 16S rRNA gene sequencing technologies to explore the community structure and co-occurrence network of the rhizosphere bacterial communities associated with the rice crop under different Cd treatments and the application of Cd-tolerant strain Pseudomonas sp. TCd-1. We found that the strain TCd-1 both significantly reduced the rhizobacterial alpha diversity and changed the beta diversity. PERMANOVA and NMDS analysis showed that Cd stress and TCd-1 strain could act as strong environmental filters resulting in observable differentiation of rhizobacterial community composition among different groups. In addition, RDA results indicated that the rhizosphere pH, root Cd content, catalase (CAT), urease (URE), gibberellic acid (GA3) exert significant association with rhizosphere bacterial assembly. PICRUSt analysis revealed that the TCd-1 strain improved the metabolic capacity of rhizosphere bacteria under Cd stress. Furthermore, co-occurrence network topological features and keystone taxa also varied among different groups. This study could provide necessary insights into developing an efficient bioremediation and safe production of rice crops in Cd contaminated paddy fields with the application of Pseudomonas sp. TCd-1 strain, as well as advance our understanding of the principles of rhizosphere bacterial community assembly under Cd stress.

Different Assembly Patterns of Planktonic and Sedimentary Bacterial Community in a Few Connected Eutrophic Lakes

The mechanism of bacterial community assembly has been the hot spot in the field of microbial ecology and it is difficult to quantitatively estimate the influences of different ecological processes. Here, a total of 23 pairs of planktonic and sedimentary samples were collected from five lakes in Wuhan, China. significant higher α-diversity (p < 0.001) and β-diversity (p < 0.001) of bacterial communities were observed in sediment than those in water. Some phylum had linear relationships with the comprehensive TSI (TSIc) by regression analysis. Non-metric multidimensional scaling (NMDS) and redundancy analysis (RDA) revealed that the depth of water, NO3−-N, NH4+-N, PO43−, and CODcr were the key environmental variables in planktonic bacterial communities, whereas in sediment they were the depth, NO3−-N, and NH4+-N. Furthermore, variation partitioning analysis (VPA) showed that spatial and environmental factors could only explain 40.2% and 27.9% of the variation in planktonic and sedimentary bacterial communities, respectively. More importantly, null model analysis suggested that different assembly mechanisms were found between in water and in sediment with the fact that planktonic bacterial community assembly was mainly driven by dispersal limitation process whereas variable selection process played a vital role in that of sediment.

Conceptual strategies for characterizing interactions in microbial communities

Understanding the sets of inter- and intraspecies interactions in microbial communities is a fundamental goal of microbial ecology. However, the study and quantification of microbial interactions pose several challenges owing to their complexity, dynamic nature, and the sheer number of unique interactions within a typical community. To overcome such challenges, microbial ecologists must rely on various approaches to distill the system of study to a functional and conceptualizable level, allowing for a practical understanding of microbial interactions in both simplified and complex systems. This review broadly addresses the role of several conceptual approaches available for the microbial ecologist’s arsenal, examines specific tools used to accomplish such approaches, and describes how the assumptions, expectations, and philosophies underlying these tools change across scales of complexity.

Dynamic changes of rhizosphere soil bacterial community and nutrients in cadmium polluted soils with soybean-corn intercropping

BACKGROUND

Soybean-corn intercropping is widely practised by farmers in Southwest China. Although rhizosphere microorganisms are important in nutrient cycling processes, the differences in rhizosphere microbial communities between intercropped soybean and corn and their monoculture are poorly known. Additionally, the effects of cadmium (Cd) pollution on these differences have not been examined. Therefore, a field experiment was conducted in Cd-polluted soil to determine the effects of monocultures and soybean-corn intercropping systems on Cd concentrations in plants, on rhizosphere bacterial communities, soil nutrients and Cd availability. Plants and soils were examined five times in the growing season, and Illumina sequencing of 16S rRNA genes was used to analyze the rhizosphere bacterial communities.

RESULTS

Intercropping did not alter Cd concentrations in corn and soybean, but changed soil available Cd (ACd) concentrations and caused different effects in the rhizosphere soils of the two crop species. However, there was little difference in bacterial community diversity for the same crop species under the two planting modes. Proteobacteria, Chloroflexi, Acidobacteria, Actinobacteria and Firmicutes were the dominant phyla in the soybean and corn rhizospheres. In ecological networks of bacterial communities, intercropping soybean (IS) had more module hubs and connectors, whereas intercropped corn (IC) had fewer module hubs and connectors than those of corresponding monoculture crops. Soil organic matter (SOM) was the key factor affecting soybean rhizosphere bacterial communities, whereas available nutrients (N, P, K) were the key factors affecting those in corn rhizosphere. During the cropping season, the concentration of soil available phosphorus (AP) in the intercropped soybean-corn was significantly higher than that in corresponding monocultures. In addition, the soil available potassium (AK) concentration was higher in intercropped soybean than that in monocropped soybean.

CONCLUSIONS

Intercropped soybean-corn lead to an increase in the AP concentration during the growing season, and although crop absorption of Cd was not affected in the Cd-contaminated soil, soil ACd concentration was affected. Intercropped soybean-corn also affected the soil physicochemical properties and rhizosphere bacterial community structure. Thus, intercropped soybean-corn was a key factor in determining changes in microbial community composition and networks. These results provide a basic ecological framework for soil microbial function in Cd-contaminated soil.

Rhizosphere bacteriome structure and functions

Microbial composition and functions in the rhizosphere—an important microbial hotspot—are among the most fascinating yet elusive topics in microbial ecology. We used 557 pairs of published 16S rDNA amplicon sequences from the bulk soils and rhizosphere in different ecosystems around the world to generalize bacterial characteristics with respect to community diversity, composition, and functions. The rhizosphere selects microorganisms from bulk soil to function as a seed bank, reducing microbial diversity. The rhizosphere is enriched in Bacteroidetes, Proteobacteria, and other copiotrophs. Highly modular but unstable bacterial networks in the rhizosphere (common for r-strategists) reflect the interactions and adaptations of microorganisms to dynamic conditions. Dormancy strategies in the rhizosphere are dominated by toxin–antitoxin systems, while sporulation is common in bulk soils. Functional predictions showed that genes involved in organic compound conversion, nitrogen fixation, and denitrification were strongly enriched in the rhizosphere (11–182%), while genes involved in nitrification were strongly depleted.

Understanding soil microbiota dynamics is key the development of soil-based sustainable agriculture and conservation strategies. This meta-analysis shows that bulk soil functions as a seed bank for the rhizosphere, which encompasses a rich microbiota adapted to dynamic conditions in hotpots.

Storage tank as a pretreatment unit for rainwater cleaner production: Role of biofilm bacterial communities and functional genera in water quality improvement

This study investigated the water purification function and mechanism of biofilm in storage tank, with a view to using it as a pretreatment unit for rainwater cleaner production. Shortening the hydraulic retention time (HRT) of storage tank from 12 to 4 h improved the pollutants removal performance and reduced the suspended bacteria counts. The accumulation of abundant taxa and succession of rare taxa were observed with biofilm development. Positive correlations within and across different bacterial taxa were dominant in the network, and some rare genera (Ralstonia and Micropruina) were identified as hub bacteria. Candidatus Nitrospira nitrosa and Nitrospira sp. ENR4 were two identified complete ammonia oxidizers. Denitrifying bacteria tended to enrich and formed more complex interactions over time. The main nitrogen metabolism pathways may be ammonia assimilatory, complete denitrification and dissimilatory/assimilatory nitrate reduction. HRT was negatively correlated with most dominant genera, and contributed 20.35% to the variation of functional taxa. This study highlights the self-purification function and micro-ecology of storage tank, and provides a new insight for its role in rainwater cleaner production process.

Effect of sludge retention time on microbial succession and assembly in thermal hydrolysis pretreated sludge digesters: Deterministic versus stochastic processes

Thermal hydrolysis process (THP) assisted anaerobic digestion (AD) has been demonstrated to be an efficient approach to improve biogas production and solids reduction. Given the faster reaction kinetics in the THP-AD system, reduction of sludge retention time (SRT) is possible. However, a comprehensive understanding of the effects of sludge retention time (SRT) on microbial dynamics and community assemblages is still lacking in THP-AD systems. Thus, twelve THP-AD reactors were operated at different SRTs (10-30 d) to fulfill the knowledge gap. Results showed that, although all the bioreactors displayed good performance, shorter SRT reactors (SRT 10 d) took a longer time to reach the stable state. The total biogas production at SRT of 10 d was lower than that at other longer SRTs, attributing to the limited hydrolytic/fermentative capacities of AD microbiomes. Different SRTs resulted in distinct succession patterns of AD microbiomes. THP sludge reduced the microbial diversity in all the bioreactors over time, but longer SRTs maintained higher biodiversity. Null model analysis suggested that THP-AD microbial community assembly was predominately driven by deterministic selection at the tested SRT range, but stochasticity increased with elevated SRTs, likely attributing to the immigrants from the feedstock. Phylogenetic molecular ecological networks (pMENs) analysis revealed more stable network structures at longer SRTs, evidenced by the lower modularity, shorter harmonic geodesic distance, and higher connectivity. The potential keystone taxa under varied SRTs were identified, some of which were hydrolytic/fermentative bacteria (e.g., Peptostreptococcus, Lutispora, Synergistaceae), suggesting that these species related to organic hydrolysis/fermentation even with low-abundance could still play pivotal ecological roles in maintaining the THP-AD microbial community structure and functions. Collectively, this study provides comprehensive and in-depth insights into the mechanisms underlying community assembly in THP-AD reactors, which could aid in diagnosing system stability.

Understanding Interaction Patterns within Deep-Sea Microbial Communities and Their Potential Applications

Environmental microbes living in communities engage in complex interspecies interactions that are challenging to decipher. Nevertheless, the interactions provide the basis for shaping community structure and functioning, which is crucial for ecosystem service. In addition, microbial interactions facilitate specific adaptation and ecological evolution processes particularly essential for microbial communities dwelling in resource-limiting habitats, such as the deep oceans. Recent technological and knowledge advancements provide an opportunity for the study of interactions within complex microbial communities, such as those inhabiting deep-sea waters and sediments. The microbial interaction studies provide insights into developing new strategies for biotechnical applications. For example, cooperative microbial interactions drive the degradation of complex organic matter such as chitins and celluloses. Such microbiologically-driven biogeochemical processes stimulate creative designs in many applied sciences. Understanding the interaction processes and mechanisms provides the basis for the development of synthetic communities and consequently the achievement of specific community functions. Microbial community engineering has many application potentials, including the production of novel antibiotics, biofuels, and other valuable chemicals and biomaterials. It can also be developed into biotechniques for waste processing and environmental contaminant bioremediation. This review summarizes our current understanding of the microbial interaction mechanisms and emerging techniques for inferring interactions in deep-sea microbial communities, aiding in future biotechnological and therapeutic applications.

An Exploratory Study on the Microbiome of Northern and Southern Populations of Ixodes scapularis Ticks Predicts Changes and Unique Bacterial Interactions

The black-legged tick (Ixodes scapularis) is the primary vector of Borrelia burgdorferi, the causative agent of Lyme disease in North America. However, the prevalence of Lyme borreliosis is clustered around the Northern States of the United States of America. This study utilized a metagenomic sequencing approach to compare the microbial communities residing within Ix. scapularis populations from northern and southern geographic locations in the USA. Using a SparCC network construction model, we performed potential interactions between members of the microbial communities from Borrelia burgdorferi-infected tissues of unfed and blood-fed ticks. A significant difference in bacterial composition and diversity was found between northern and southern tick populations. The network analysis predicted a potential antagonistic interaction between endosymbiont Rickettsia buchneri and Borrelia burgdorferi sensu lato. The network analysis, as expected, predicted significant positive and negative microbial interactions in ticks from these geographic regions, with the genus Rickettsia, Francisella, and Borreliella playing an essential role in the identified clusters. Interactions between Rickettsia buchneri and Borrelia burgdorferi sensu lato need more validation and understanding. Understanding the interplay between the microbiome and tick-borne pathogens within tick vectors may pave the way for new strategies to prevent tick-borne infections.

New insights into metabolomics profile generation in fermented tea: the relevance of bacteria and metabolites in Fuzhuan brick tea

BACKGROUND

The contribution of bacteria to fermented tea is not clear and the associated research is relatively limited. To reveal the role of microorganisms in fermented tea processing, the microbial community and metabolites of Fuzhuan brick tea (FBT), a Chinese traditional fermented tea, were revealed via high-throughput sequencing and liquid chromatography-mass spectrometry (LC-MS).

RESULTS

In FBT, bacterial communities had a higher abundance and diversity, Lactococcus and Bacillus were the main bacteria, and Eurotium was the predominant fungus. The predictive metabolic function indicated the pathways of cellular growth, environmental information, genetics and material metabolism of bacterial communities were abundant, whereas the fungal community predictive metabolic function was almost saprotroph. Using LC-MS, 1143 and 536 metabolites were defined in positive and negative ion mode, respectively. There were essential correlations between bacterial populations and metabolites, such that Bacillus was correlated significantly with 44 metabolites (P < 0.05) and Enterococcus was significantly associated with 15 metabolites (P < 0.05). Some of the main active components were significantly correlated with the bacteria, such as Enterococcus, Lactococcus and Carnobacterium.

CONCLUSION

Not only Eurotium, but also the bacteria were involved in the changes of metabolomics profile in fermented FBT. The present study assists in providing new insights into metabolomics profile generation in fermented tea. The present research lays a foundation for controlling the FBT fermentation by artificial inoculation to improve quality. © 2021 Society of Chemical Industry.

Soil Bacterial Community Shifts Are Driven by Soil Nutrient Availability along a Teak Plantation Chronosequence in Tropical Forests in China

Simple Summary

Tropical forests play an important role in the global carbon cycle, especially in the context of global climate change. Soil microorganisms are essential to the functions, services, and productivity of terrestrial ecosystems as a link to maintain the connections and interactions between the aboveground and belowground ecosystems. The interactions between plants and the soil microbiome are crucial for plant growth, health, and resistance to stressors. However, information on the response of soil microbial communities to a chronosequence of woody plants is lacking, especially in tropical forests. This study compares the soil properties, diversity, composition, and co-occurrence patterns of bacterial communities in the rhizosphere and bulk soils along a teak plantation chronosequence. The results show that the composition and co-occurrence patterns of the bacterial communities are statistically different among the plantations, while stand age has no significant impact on soil bacterial alpha diversity. The results further show that soil nutrients play a key role in shaping the soil bacterial community. The study also provides information about the dynamics and characteristics of these soil bacterial communities and adds valuable information that may underpin new strategies for the management of teak plantations.

Abstract

Soil bacterial communities play crucial roles in ecosystem functions and biogeochemical cycles of fundamental elements and are sensitive to environmental changes. However, the response of soil bacterial communities to chronosequence in tropical ecosystems is still poorly understood. This study characterized the structures and co-occurrence patterns of soil bacterial communities in rhizosphere and bulk soils along a chronosequence of teak plantations and adjacent native grassland as control. Stand ages significantly shifted the structure of soil bacterial communities but had no significant impact on bacterial community diversity. Bacterial community diversity in bulk soils was significantly higher than that in rhizosphere soils. The number of nodes and edges in the bacterial co-occurrence network first increased and then decreased with the chronosequence. The number of strongly positive correlations per network was much higher than negative correlations. Available potassium, total potassium, and available phosphorus were significant factors influencing the structure of the bacterial community in bulk soils. In contrast, urease, total potassium, pH, and total phosphorus were significant factors affecting the structure of the bacterial community in the rhizosphere soils. These results indicate that available nutrients in the soil are the main drivers regulating soil bacterial community variation along a teak plantation chronosequence.

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

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

Solid-phase denitrification in high salinity and low-temperature wastewater treatment

Nitrogen removal from wastewater is often deteriorated under high salinity and low temperature. Solid-phase denitrification (SPD) might improve total nitrogen removal efficiency (TNRE) by stably supplying carbon resources under adverse conditions. In this study, an SPD biofilm reactor was successfully established by inoculating halophilic sludge and filling poly (butanediol succinate) (PBS) granules, and achieved over 96% TNRE at low temperature. More extracellular polysaccharides were produced at low temperature. Microbial network analysis evidenced dominant heterotrophic denitrifiers (Marinicella, Fusibacter, Saccharicrinis and Vitellibacter) at 25 °C were replaced by genera Melioribacter, Marinobacter, Desulfatitalea and Thiomicrospira at 15 °C. At low temperature, genes nirS and narG might be mainly responsible for denitrification. Fluorescence spectrum coupled with fluorescence regional integration and parallel factor analysis revealed low temperature increased the proportion of proteins of soluble microbial products. This study provides guidance for the practical application of SPD in the treatment of high salinity and low-temperature wastewater.

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.

RNA-Seq Reveals that Light and Darkness Are Different Stimuli in Freshwater Heterotrophic Actinobacteria

Light is a ubiquitous source of both energy and information in surface environments, and regulates gene expression not only in photosynthetic microorganisms, but in a broad range of photoheterotrophic and heterotrophic microbes as well. Actinobacteria are keystone species in surface freshwater environments, where the ability to sense light could allow them to coordinate periods of nutrient uptake and metabolic activity with primary production. The model freshwater Actinobacteria Rhodoluna (R.) lacicola strain MWH-Ta8 and Aurantimicrobium (A.) photophilum strain MWH-Mo1 grow faster in the light than in the dark, but do not use light energy to support growth. Here, we characterize transcription throughout a light-dark cycle in R. lacicola and A. photophilum. In both species, some genes encoding carbohydrate metabolism and storage are upregulated in the light. However, expression of genes of the TCA cycle is only coordinated with light availability in R. lacicola. In fact, the majority of genes that respond to light and darkness in these two species are different, even though their light-responsive phenotypes are similar. The ability to respond to light and darkness may be widespread in freshwater Actinobacteria, but the genetic networks controlled by these two stimuli may vary significantly.

Impact of different types of anthropogenic pollution on bacterial community and metabolic genes in urban river sediments

Sediment bacterial communities play a crucial role in the biogeochemical cycle of nutrient elements in urban river. However, the distribution of nitrogen cycle genes on bacterial communities in urban rivers sediments is largely unknown. Here, 16S rRNA amplicon sequencing was used to analyze the composition, co-occurrence patterns and nitrogen cycle process of bacterial communities in urban river sediments under the influence of different exogenous pollution. The results revealed that bacterial communities had significant spatial heterogeneity in river sediments of different polluted areas, and the input of different exogenous pollutants shaped the abundance and distribution of nitrogen cycle-related genes in the sediments. In addition, denitrification process played a leading role in the nitrogen cycle of river sediments, and the genes associated with the nitrification process were rarely observed in all samples. The important bacterial taxonomic biomarkers of nitrogen cycling-related genes screened by random forest algorithm were Synergistia, WS6_Dojkabacteria and Caldisericia. Meanwhile, different co-occurrence patterns observed in different types of polluted areas clarified the impact of environmental filtration and niche differentiation on bacterial communities. In conclusion, this study reveals the nitrogen cycle process and the distribution of related genes mediated by bacterial communities under the impact of different anthropogenic contamination, and provides novel insights for the assembly of bacterial communities in urban river sediments.

Characterization of tetracycline-resistant microbiome in soil-plant systems by combination of H218O-based DNA-Stable isotope probing and metagenomics

The emergence and spread of antibiotic resistance have been considered as a global health threat. However, effective methods to identify antibiotic-resistant bacteria (ARB) in complex microbial community are lacking, and the potential transmission pathways of ARB and antibiotic resistance genes (ARGs) in the soil-plant system remain scarce. Here in this study, tetracycline was chosen as the target antibiotic due to its globally wide usage and clinical significance. DNA-based stable isotope probing with H₂¹⁸O was applied to identify the tetracycline-resistant bacteria from soil-plant systems. Eighteen-year organic fertilization significantly changed the composition of the tetracycline-resistant microbiome in the soil-wheat system and resulted in a higher relative abundance of ARGs in the wheat endophyte. Rhizosphere harboring the most diverse ARGs and mobile genetic elements was identified as a hot spot for horizontal gene transfer and an important bridge between bulk soil and wheat endophyte. Micrococcaceae and Sphingomonadaceae carrying ARGs associated with abundant mobile genetic elements, were identified as the core bacterial taxa in long-term manure-amended and untreated soil-wheat systems, respectively. This method contributes to a more precise track of ARB in the environment, and our work depicts the high potential of ARG transfer in the rhizosphere mediated by the core species.

Dynamic carbon flux network of a diverse marine microbial community

The functioning of microbial ecosystems has important consequences from global climate to human health, but quantitative mechanistic understanding remains elusive. The components of microbial ecosystems can now be observed at high resolution, but interactions still have to be inferred e.g., a time-series may show a bloom of bacteria X followed by virus Y suggesting they interact. Existing inference approaches are mostly empirical, like correlation networks, which are not mechanistically constrained and do not provide quantitative mass fluxes, and thus have limited utility. We developed an inference method, where a mechanistic model with hundreds of species and thousands of parameters is calibrated to time series data. The large scale, nonlinearity and feedbacks pose a challenging optimization problem, which is overcome using a novel procedure that mimics natural speciation or diversification e.g., stepwise increase of bacteria species. The method allows for curation using species-level information from e.g., physiological experiments or genome sequences. The product is a mass-balancing, mechanistically-constrained, quantitative representation of the ecosystem. We apply the method to characterize phytoplankton-heterotrophic bacteria interactions via dissolved organic matter in a marine system. The resulting model predicts quantitative fluxes for each interaction and time point (e.g., 0.16 µmolC/L/d of chrysolaminarin to Polaribacter on April 16, 2009). At the system level, the flux network shows a strong correlation between the abundance of bacteria species and their carbon flux during blooms, with copiotrophs being relatively more important than oligotrophs. However, oligotrophs, like SAR11, are unexpectedly high carbon processors for weeks into blooms, due to their higher biomass. The fraction of exudates (vs. grazing/death products) in the DOM pool decreases during blooms, and they are preferentially consumed by oligotrophs. In addition, functional similarity of phytoplankton i.e., what they produce, decouples their association with heterotrophs. The methodology is applicable to other microbial ecosystems, like human microbiome or wastewater treatment plants.

Bacterial community dynamics on bats and the implications for pathogen resistance

The bats skin microbiota plays an important role in reducing pathogen infection, including the deadly fungal pathogen Pseudogymnoascus destructans, the causative agent of white-nose syndrome. However, the dynamic of skin bacterial communities response to environmental perturbations remains poorly described. We characterized skin bacterial community over time and space in Rhinolophus ferrumequinum, a species with high resistance to the infection with P. destructans. We collected environmental covariate data to determine what factors influenced changes in community structure. We observed significant temporal and spatial shifts in the skin bacterial community, which was mainly associated with variation in operational taxonomic units. The skin bacterial community differed by the environmental microbial reservoirs and was most influenced by host body condition, bat roosting temperature and geographic distance between sites, but was not influenced by pathogen infection. Furthermore, the skin microbiota was enriched in particular taxa with antifungal abilities, such as Enterococcus, Burkholderia, Flavobacterium, Pseudomonas, Corynebacterium and Rhodococcus. And specific strains of Pseudomonas, Corynebacterium and Rhodococcus even inhibited P. destructans growth. Our findings provide new insights in characterizing the variation in bacterial communities can inform us about the processes of driving community assembly and predict the host's ability to resist or survive pathogen infection.

Lake microbial communities are not resistant or resilient to repeated large-scale natural pulse disturbances

Opportunities to study community level responses to extreme natural pulse disturbances in unaltered ecosystems are rare. Lake sediment records that span thousands of years can contain well resolved sediment pulses, triggered by earthquakes. These paleo-records provide a means to study repeated pulse disturbance and processes of resistance (insensitivity to disturbance) and ecological resilience (capacity to regain structure, function and process). In this study, sedimentary DNA was extracted from a sediment core from Lake Paringa (New Zealand) that is situated in a near natural catchment. Metabarcoding and inferred functions were used to assess the lake microbial community over the past 1,100 years - a period that included four major earthquakes. Microbial community composition and function differed significantly between highly perturbed (postseismic, c. 50 yrs) phases directly after the earthquakes and more stable (interseismic, c. 250 yr) phases, indicating a lack of community resistance. Although community structure differed significantly in successive postseismic phases, function did not, suggesting potential functional redundancy. Significant differences in composition and function in successive interseismic phases demonstrates communities are not resilient to large-scale natural pulse disturbances. The clear difference in structure and function, and high number of indicator taxa (responsible for driving differences in communities between phases) in the fourth interseismic phase likely represents a regime shift, possibly due to the two-fold increase in sediment and terrestrial biospheric organic carbon fluxes recorded following the fourth earthquake. Large pulse disturbances that enhance sediment inputs into lake systems may produce an underappreciated mechanism that destabilises lake ecosystem processes.