Contrasting Network Features between Free-Living and Particle-Attached Bacterial Communities in Taihu Lake

Response of particle-attached and free-living bacterial communities to Microcystis blooms

The massive proliferation of Microcystis threatens freshwater ecosystems and degrades water quality globally. Understanding the mechanisms that contribute to Microcystis growth is crucial for managing Microcystis blooms. The lifestyles of bacteria can be classified generally into two groups: particle-attached (PA; > 3 µm) and free-living (FL; 0.2-3.0 µm). However, little is known about the response of PA and FL bacteria to Microcystis blooms. Using 16S rRNA gene high-throughput sequencing, we investigated the stability, assembly process, and co-occurrence patterns of PA and FL bacterial communities during distinct bloom stages. PA bacteria were phylogenetically different from their FL counterparts. Microcystis blooms substantially influenced bacterial communities. The time decay relationship model revealed that Microcystis blooms might increase the stability of both PA and FL bacterial communities. A contrasting community assembly mechanism was observed between the PA and FL bacterial communities. Throughout Microcystis blooms, homogeneous selection was the major assembly process that impacted the PA bacterial community, whereas drift explained much of the turnover of the FL bacterial community. Both PA and FL bacterial communities could be separated into modules related to different phases of Microcystis blooms. Microcystis blooms altered the assembly process of PA and FL bacterial communities. PA bacterial community appeared to be more responsive to Microcystis blooms than FL bacteria. Decomposition of Microcystis blooms may enhance cooperation among bacteria. Our findings highlight the importance of studying bacterial lifestyles to understand their functions in regulating Microcystis blooms. KEY POINTS: • Microcystis blooms alter the assembly process of PA and FL bacterial communities • Microcystis blooms increase the stability of both PA and FL bacterial communities • PA bacteria seem to be more responsive to Microcystis blooms than FL bacteria.

Community structure and carbon metabolism functions of bacterioplankton in the Guangdong coastal zone

Coastal ecosystems are an important region for biogeochemical cycling, are a hotspot of anthropogenic disturbance and play a crucial role in global carbon cycling through the metabolic activities of bacterioplankton. Bacterioplankton can be broadly classified into two lifestyles: free-living (FL) and particle-attached (PA). However, how coastal bacterioplankton the community structure, co-occurrence networks and carbon metabolic functions with different lifestyles are differentiated is still largely unknown. Understanding these processes is necessary to better determine the contributions of coastal bacterioplankton to carbon cycling. Here, the characteristics of community structure and carbon metabolism function of bacterioplankton with two lifestyles in the coastal areas of Guangdong Province were investigated using amplicon sequencing, metagenomic, and metatranscriptomic techniques. The results show that the main bacterioplankton responsible for carbon metabolism were the Pseudomonadota, Bacteroidota, and Actinomycetota. The microbial community structure, carbon metabolic function, and environmental preferences differ between different lifestyles. FL and PA bacteria exhibited higher carbon fixation and degradation potentials, respectively. A range of environmental factors, such as dissolved oxygen, pH, and temperature, were associated with the community structure and carbon metabolic functions of the bacterioplankton. Human activities, such as nutrient discharge, may affect the distribution of functional genes and enhance the carbon degradation functions of bacterioplankton. In conclusion, this study increased the understanding of the role of microorganisms in regulating carbon export in coastal ecosystems with intense human activity.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00245-x.

Insights into cyanobacterial blooms through the lens of omics

Cyanobacteria are oxygen-producing photosynthetic bacteria that convert carbon dioxide into biomass upon exposure to sunlight. However, favorable conditions cause harmful cyanobacterial blooms (HCBs), which are the dense accumulation of biomass at the water surface or subsurface, posing threats to freshwater ecosystems and human health. Understanding the mechanisms underlying cyanobacterial bloom formation is crucial for effective management. In this regard, recent advancements in omics technologies have provided valuable insights into HCBs, which have raised expectations to develop more effective control methods in the near future. This literature review aims to present the genomic architecture, adaptive mechanisms, microbial interactions, and ecological impacts of HCBs through the lens of omics. Genomic analysis indicates that the genome plasticity of cyanobacteria has enabled their resilience and effective adaptation to environmental changes. Transcriptomic investigations have revealed that cyanobacteria use various strategies for adapting to environmental stress. Additionally, metagenomic and metatranscriptomic analyses have emphasized the significant role of the microbial community in regulating HCBs. Finally, we offer perspectives on potential opportunities for further research in this field.

Aerobic anoxygenic phototrophs play important roles in nutrient cycling within cyanobacterial Microcystis bloom microbiomes

BACKGROUND

During the bloom season, the colonial cyanobacterium Microcystis forms complex aggregates which include a diverse microbiome within an exopolymer matrix. Early research postulated a simple mutualism existing with bacteria benefitting from the rich source of fixed carbon and Microcystis receiving recycled nutrients. Researchers have since hypothesized that Microcystis aggregates represent a community of synergistic and interacting species, an interactome, each with unique metabolic capabilities that are critical to the growth, maintenance, and demise of Microcystis blooms. Research has also shown that aggregate-associated bacteria are taxonomically different from free-living bacteria in the surrounding water. Moreover, research has identified little overlap in functional potential between Microcystis and members of its microbiome, further supporting the interactome concept. However, we still lack verification of general interaction and know little about the taxa and metabolic pathways supporting nutrient and metabolite cycling within Microcystis aggregates.

RESULTS

During a 7-month study of bacterial communities comparing free-living and aggregate-associated bacteria in Lake Taihu, China, we found that aerobic anoxygenic phototrophic (AAP) bacteria were significantly more abundant within Microcystis aggregates than in free-living samples, suggesting a possible functional role for AAP bacteria in overall aggregate community function. We then analyzed gene composition in 102 high-quality metagenome-assembled genomes (MAGs) of bloom-microbiome bacteria from 10 lakes spanning four continents, compared with 12 complete Microcystis genomes which revealed that microbiome bacteria and Microcystis possessed complementary biochemical pathways that could serve in C, N, S, and P cycling. Mapping published transcripts from Microcystis blooms onto a comprehensive AAP and non-AAP bacteria MAG database (226 MAGs) indicated that observed high levels of expression of genes involved in nutrient cycling pathways were in AAP bacteria.

CONCLUSIONS

Our results provide strong corroboration of the hypothesized Microcystis interactome and the first evidence that AAP bacteria may play an important role in nutrient cycling within Microcystis aggregate microbiomes. Video Abstract.

Microbial interaction patterns and nitrogen cycling regularities in lake sediments under different trophic conditions

Exploring how nitrogen (N) cycling microbes interact in eutrophic lake sediments and how biogenic elements influence the nitrogen cycle is crucial for understanding biogeochemical cycles and nitrogen accumulation mechanisms. In this study, sediment samples were collected from various areas of Taihu Lake with different trophic conditions in all four seasons from 2015 to 2017. Using high-throughput sequencing and molecular ecological network analysis, we investigated the microbial interaction patterns and the role of nitrogen cycling in sediments from lakes with different trophic conditions. The results showed distinct structures of sediment microbial networks between lake areas with different trophic conditions. In the more eutrophic region, network indices indicate higher transfer efficiency of energy, material, and information, more significant competition, and weaker niche differentiation of the microbial community. The sedimentary environment in the moderately eutrophic area exhibited greater potential for denitrification, nitrification, and anammox compared to the mesotrophic area, but the inhibition between N functional microbes and limitations in N removal processes were also more likely to occur. The topological structure of the networks showed that the carbon (C), sulfur (S), and iron (Fe) cycles had a strong influence on the nitrogen cycle in both lake areas. In the moderately eutrophic lake area, C- and S-cycling functional bacteria facilitated a closed cycle of the coupled N fixation-nitrification-DNRA (dissimilatory nitrate reduction to ammonium) process and reduced N removal. In the mesotrophic lake area, C- and S-cycling functional bacteria promoted both N fixation and mineralization, and Fe-cycling functional bacteria coupled with denitrifiers enhanced the nitrogen removal process of products from nitrogen fixation and mineralization. This study improved the understanding of the nitrogen cycling mechanism in lake sediments under different trophic conditions.

Variation in resource competition traits among Microcystis strains is affected by their microbiomes

Freshwater harmful algal blooms are often dominated by Microcystis, a phylogenetically cohesive group of cyanobacteria marked by extensive genetic and physiological diversity. We have previously shown that this genetic diversity and the presence of a microbiome of heterotrophic bacteria influences competitive interactions with eukaryotic phytoplankton. In this study, we sought to explain these observations by characterizing Monod equation parameters for resource usage (maximum growth rate μ max, half-saturation value for growth K s, and quota) as a function of N and P levels for four strains (NIES-843, PCC 9701, PCC 7806 [WT], and PCC 7806 ΔmcyB) in presence and absence of a microbiome derived from Microcystis isolated from Lake Erie. Results indicated limited differences in maximum growth rates but more pronounced differences in half-saturation values among Microcystis strains. The largest impact of the microbiome was reducing the minimal nitrogen concentration sustaining growth and reducing half saturation values, with variable results depending on the Microcystis strain. Microcystis strains also differed from each other in their N and P quotas and the extent to which microbiome presence affected them. Our data highlight the importance of the microbiome in altering Microcystis-intrinsic traits, strain competitive hierarchies, and thus bloom dynamics. As quota, μ max, and K s are commonly used in models for harmful algal blooms, our data suggest that model improvement may be possible by incorporating genotype dependencies of resource-use parameters.

The Assembly Process of Free-Living and Particle-Attached Bacterial Communities in Shrimp-Rearing Waters: The Overwhelming Influence of Nutrient Factors Relative to Microalgal Inoculation

The ecological functions of bacterial communities vary between particle-attached (PA) lifestyles and free-living (FL) lifestyles, and separately exploring their community assembly helps to elucidate the microecological mechanisms of shrimp rearing. Microalgal inoculation and nutrient enrichment during shrimp rearing are two important driving factors that affect rearing-water bacterial communities, but their relative contributions to the bacterial community assembly have not been evaluated. Here, we inoculated two microalgae, Nannochloropsis oculata and Thalassiosira weissflogii, into shrimp-rearing waters to investigate the distinct effects of various environmental factors on PA and FL bacterial communities. Our study showed that the composition and representative bacteria of different microalgal treatments were significantly different between the PA and FL bacterial communities. Regression analyses and Mantel tests revealed that nutrients were vital factors that constrained the diversity, structure, and co-occurrence patterns of both the PA and FL bacterial communities. Partial least squares path modeling (PLS-PM) analysis indicated that microalgae could directly or indirectly affect the PA bacterial community through nutrient interactions. Moreover, a significant interaction was detected between PA and FL bacterial communities. Our study reveals the unequal effects of microalgae and nutrients on bacterial community assembly and helps explore microbial community assembly in shrimp-rearing ecosystems.

Distinct community assembly processes and habitat specialization driving the biogeographic patterns of abundant and rare bacterioplankton in a brackish coastal lagoon

The ecological diversity patterns and community assembly processes along spatio-temporal scales are least studied in the bacterioplankton sub-communities of brackish coastal lagoons. We examined the biogeographic patterns and relative influences of different assembly processes in structuring the abundant and rare bacterioplankton sub-communities of Chilika, the largest brackish water coastal lagoon of India. Rare taxa demonstrated significantly higher α- and β-diversity and biogeochemical functions than abundant taxa in the high-throughput 16S rRNA gene sequence dataset. The majority of the abundant taxa (91.4 %) were habitat generalists with a wider niche breadth (niche breadth index, B = 11.5), whereas most of the rare taxa (95.2 %) were habitat specialists with a narrow niche breadth (B = 8.9). Abundant taxa exhibited a stronger distance-decay relationship and higher spatial turnover rate than rare taxa. β-diversity partitioning revealed that the contribution of species turnover (72.2-97.8 %) was greater than nestedness (2.2-27.8 %) in causing the spatial variation in both abundant and rare taxa. Null model analyses revealed that the distribution of abundant taxa was mostly structured by stochastic processes (62.8 %), whereas deterministic processes (54.1 %) played a greater role in the rare taxa. However, the balance of these two processes varied across spatio-temporal scales in the lagoon. Salinity was the key deterministic factor controlling the variation of both abundant and rare taxa. Potential interaction networks showed a higher influence of negative interactions, indicating that species exclusion and top-down processes played a greater role in the community assembly. Notably, abundant taxa emerged as keystone taxa across spatio-temporal scales, suggesting their greater influences on other bacterial co-occurrences and network stability. Overall, this study provided detailed mechanistic insights into biogeographic patterns and underlying community assembly processes of the abundant and rare bacterioplankton over spatio-temporal scales in a brackish lagoon.

First report on the bacterial community composition, diversity, and functions in Ramsar site of Central Himalayas, Nepal

Wetland bacterial communities are highly sensitive to altered hydrology and the associated change in water physicochemical and biological properties leading to shifts in community composition and diversity, hence affecting the ecological roles. However, relevant studies are lacking in the wetlands of central Himalayas Nepal. Thus, we aimed to explore the variation of bacterial communities, diversity, and ecologic functions in the wet and dry periods of a wetland (designed as Ramsar site, Ramsar no 2257) by using 16S rRNA gene-based Illumina MiSeq sequencing. We reported a pronounced variation in water physicochemical and biological properties (temperature, pH, Chla, DOC, and TN), bacterial diversity, and community composition. Bacterial communities in the dry season harbored significantly higher alpha diversity, while significantly higher richness and abundance were reflected in the wet season. Our results uncovered the effect of nutrients on bacterial abundance, richness, and community composition. Fourteen percent of the total OTUs were shared in two hydrological periods, and the largest portion of unique OTUs (58%) was observed in the dry season. Planctomycetes and Bacteroidetes dominated the wet season exclusive OTUs; meanwhile, Actinobacteria dominated the dry season exclusive OTUs. Bacteria in these wetlands exhibited divergent ecological functions during the dry and wet seasons. By disclosing the variation of water bacterial communities in different hydrologic periods and their relationship with environmental factors, this first-hand work in the Ramsar site of Nepal will develop a baseline dataset for the scientific community that will assist in understanding the wetland's microbial ecology and biogeography.

Habitats modulate influencing factors shaping the spatial distribution of bacterial communities along a Tibetan Plateau riverine wetland

The Tibetan Plateau riverine wetland is very sensitive to global climate change. Understanding the mechanisms that maintain the spatial patterns of bacterial communities provides insight into the dominant biogeochemical processes within the plateau riverine wetlands. Nonetheless, the spatial distribution of bacterial communities along these wetlands has rarely been explored. We investigated the spatial patterns of bacterial community within rhizosphere soil, bulk soil, and sediment samples collected along the Yarlung Tsangpo riverine wetland (YTRW), the longest plateau riverine wetland in China. Our results indicated that the diversity of bacterial communities in all three habitats increased significantly along the YTRW. The slope of the linear relationship between distance and bacterial community diversity in sediment was steeper than those for bulk and rhizosphere soils. Furthermore, bacterial communities in all three habitats showed significant distance-decay relationships. A combination of historical factors (geographical distance and climatic factors) and contemporary environmental heterogeneity (edaphic properties) controlled spatial distributions of bacterial communities in all three habitats, although climatic factors were predominant. Climatic factors affected rhizosphere bacterial communities more than those in bulk soil and sediment. Co-occurrence network analysis revealed that the potential interactions between bacterial taxa may decrease along the YTRW. This field investigation highlighted that the climatic factors strongly influenced the spatial distribution of bacterial communities along the YTRW; however, habitat differences among rhizosphere soil, bulk soil, and sediment samples affected the relative importance of climatic factors on spatial distributions of the associated bacterial communities. These findings would improve the understanding of biogeochemical processes in these typical habitats and potential alterations provoked by climate change.

Spatiotemporal variations, assembly processes, and co-occurrence patterns of particle-attached and free-living bacteria in a large drinking water reservoir in China

Particle-attached (PA) and free-living (FL) bacterial communities are sensitive to pollutant concentrations and play an essential role in biogeochemical processes and water quality maintenance in aquatic ecosystems. However, the spatiotemporal variations, assembly processes, co-occurrence patterns, and environmental interactions of PA and FL bacteria in drinking water reservoirs remain as yet unexplored. To bridge this gap, we collected samples from 10 sites across four seasons in Lake Tianmu, a large drinking water reservoir in China. Analysis of 16S rRNA gene libraries demonstrated spatiotemporal variations in bacterial diversity and identified differences in bacterial community composition (BCC) between PA and FL lifestyles. Capacity for nitrogen respiration, nitrogen fixation, and nitrate denitrification was enriched in the PA lifestyle, while photosynthesis, methylotrophy, and methanol oxidation were enriched in the FL lifestyle. Deterministic processes, including interspecies interactions and environmental filtration, dominated the assembly of both PA and FL bacterial communities. The influence of environmental filtration on the FL community was stronger than that on the PA community, indicating that bacteria in the FL lifestyle were more sensitive to environmental variation. Co-occurrence patterns and keystone taxa differed between PA and FL lifestyles. The ecological functions of keystone taxa in the PA lifestyle were associated with the supply and recycling of nutrients, while those in FL were associated with the degradation of complex pollutants. PA communities were more stable than FL communities in the face of changing environmental conditions. Nutrients (e.g., TDN and NO₃ ^-) and abiotic and biotic factors (e.g., WT and Chl-a) exerted positive and negative effects, respectively, on the co-occurrence networks of both lifestyles. These results improve our understanding of assembly processes, co-occurrence patterns, and environmental interactions within PA and FL communities in a drinking water reservoir.

Microcystis pangenome reveals cryptic diversity within and across morphospecies

Microcystis, a common harmful algal bloom (HAB) taxon, threatens water supplies and human health, yet species delimitation is contentious in this taxon, leading to challenges in research and management of this threat. Historical and common morphology-based classifications recognize multiple morphospecies, most with variable and diverse ecologies, while DNA sequence-based classifications indicate a single species with multiple ecotypes. To better delimit Microcystis species, we conducted a pangenome analysis of 122 genomes. Core- and non-core gene phylogenetic analyses placed 113 genomes into 23 monophyletic clusters containing at least two genomes. Overall, genome-related indices revealed that Microcystis contains at least 16 putative genospecies. Fifteen genospecies included at least one Microcystis aeruginosa morphospecies, and 10 genospecies included two or more morphospecies. This classification system will enable consistent taxonomic identification of Microcystis and thereby aid in resolving some of the complexities and controversies that have long characterized eco-evolutionary research and management of this important HAB taxon.

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

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

Temporal and functional interrelationships between bacterioplankton communities and the development of a toxigenic Microcystis bloom in a lowland European reservoir

The cyanobacteria-associated microbiome is constantly reshaped by bloom development. However, the synergistic-antagonistic nature of the relationships between Microcystis and its microbiome still remains unclear. Therefore, temporal changes of bacterioplankton communities and their functional potential through different developing stages of a Microcystis toxigenic bloom were investigated, considering bacterioplankton assemblages as particle-attached (PAB) and free-living (FLB) bacteria. 16S rRNA sequencing revealed that PAB were represented by Proteobacteria and Cyanobacteria, while FLB by Proteobacteria and Actinobacteria. Network and ordination analyses indicated that PAB inter-relationships were more complex-numerous connections between taxa with stronger correlations, than FLB-rather influenced by physico-chemical parameters. PAB in pre-summer was diverse with Proteobacteria containing potential taxa involved in nitrogen-transforming processes. In mid-summer, PAB presented a mix-bloom dominated by Snowella, Aphanizomenon, and Microcystis, which were succeeded by toxigenic Microcystis in post-summer. Both periods were associated to potential taxa with parasitic/predatory lifestyles against cyanobacteria. In post-summer, Sutterellaceae were recognized as poor water quality indicators, and their strong association with Microcystis could have represented an increased threat for that period. Microcystis was a major factor significantly reducing PAB diversity and evenness, suggesting that it negatively influenced bacterioplankton assemblages, probably also altering the overall community functional potential.

Salinity and seasonality shaping free-living and particle-associated bacterioplankton community assembly in lakeshores of the northeastern Qinghai-Tibet Plateau

Microorganisms in lakeshore zones are essential for pollution interception and biodiversity maintenance. However, the biogeographic patterns of bacterioplankton communities in lakeshore zones and the mechanisms that driving them are poorly understood. We analyzed the 16 S rRNA gene sequences of particle-associated (PA) and free-living (FL) bacterioplankton communities in the lakeshore zones of 14 alpine lakes in two seasons on Qinghai-Tibet Plateau to investigate the bacterial diversity, composition and assembly processes. Our results revealed that PA and FL bacterioplankton communities were driven by both seasonality and salinity in the lakeshores on Qinghai-Tibet Plateau. Compared to FL bacterioplankton, PA bacterioplankton communities were more susceptible to seasonality than spatial salinity. FL bacterioplankton communities were more salinity constrained than the PA counterpart. Besides, the Stegen null model analyses have validated a quantitative bias on stochastic processes at different spatial scales. At a regional scale, stochasticity was the predominant assembly process in both PA and FL bacterioplankton. While at a subregional scale, dispersal limitation was the main contributor of stochastic processes for PA bacterioplankton in summer and heterogeneous selection was the dominant deterministic processes in winter, whereas the community assembly of FL bacterioplankton was more stochastic processes (i.e., dispersal limitation) dominated in the freshwater type but deterministic process (i.e., heterogeneous selection) increased with increasing salinity. Our study provides new insights into both significant spatiotemporal patterns and distinct assembly processes of PA and FL bacterioplankton in alpine lakeshores on the northeastern Qinghai-Tibet Plateau.

Comparison of assembly process and co-occurrence pattern between planktonic and benthic microbial communities in the Bohai Sea

Unraveling the mechanisms structuring microbial community is a central goal in microbial ecology, but a detailed understanding of how community assembly processes relate to living habitats is still lacking. Here, via 16S rRNA gene amplicon sequencing, we investigated the assembly process of microbial communities in different habitats [water verse sediment, free-living (FL) verse particle-associated (PA)] and their impacts on the inter-taxa association patterns in the coastal Bohai Sea, China. The results showed clear differences in the composition and diversity of microbial communities among habitats, with greater dissimilarities between water column and sediment than between FL and PA communities. The microbial community assembly was dominated by dispersal limitation, ecological drift, and homogeneous selection, but their relative importance varied in different habitats. The planktonic communities were mainly shaped by dispersal limitation and ecological drift, whereas homogeneous selection played a more important role in structuring the benthic communities. Furthermore, the assembly mechanisms differed between FL and PA communities, especially in the bottom water with a greater effect of ecological drift and dispersal limitation on the FL and PA fractions, respectively. Linking assembly process to co-occurrence pattern showed that the relative contribution of deterministic processes (mainly homogeneous selection) increased under closer co-occurrence relationships. By contrast, stochastic processes exerted a higher effect when there were less inter-taxa connections. Overall, our findings demonstrate contrasting ecological processes underpinning microbial community distribution in different habitats including different lifestyles, which indicate complex microbial dynamic patterns in coastal systems with high anthropogenic perturbations.

Habitat-dependent prokaryotic microbial community, potential keystone species, and network complexity in a subtropical estuary

Microbes (e.g., bacteria and archaea) are indispensable components for the key biological processes of estuarine ecosystems and three main habitats (sediment, particle, and water) are harboring diverse estuarine microbes. However, we still know little about how the microbial community structures, potential keystone species, and network properties change among these three habitats in estuarine ecosystems. In this study, we collected size-fractioned water and sediment samples from the Pearl River Estuary to reveal their microbial diversity, community structures, network properties, and potential keystone taxa. We found that the sediment microbial community was remarkably more diverse than particle-attached (PA) and free-living (FL) communities, whereas its ecological network was less complex in terms of node distance and connectivity. TOC was determined as the main driver of sediment community, while the PA and FL communities were predominantly shaped by NO₂^-, non-ionic ammonia (NH) and pH. Among the bulk water, there were no significant differences between PA and FL communities in diversity, community structure, and network complexity. However, the PA community was more susceptible to metal elements, suggesting their higher level of involvement in physiological metabolism. Potential keystone taxa among community networks were taxonomically divergent in three habitats. Specifically, Synechococcales (Cyanobacteria) and Actinomarinales (Actinobacteria) exclusively served as the module-hubs in FL network, while members from phylum Proteobacteria and Bacteroidetes were the module-hubs and connectors in PA network. Potential keystone taxa in sediment network were more diverse and covered 9 phyla, including the only archaeal lineage Bathyarchaeia (Crenarchaeota). Overall, our study provided more detailed information about estuarine microbial communities in three habitats, especially the potential keystone species, which provided new perspectives on evaluating further effects of anthropogenic disturbances on estuarine microbes and facilitated the environment monitoring based on microbial community.

Cyanopeptides restriction and degradation co-mediate microbiota assembly during a freshwater cyanobacterial harmful algal bloom (CyanoHAB)

Cyanobacterial harmful algal blooms (CyanoHABs) are globally intensifying and exacerbated by climate change and eutrophication. However, microbiota assembly mechanisms underlying CyanoHABs remain elusive. Especially, cyanopeptides, as a group of bioactive secondary metabolites of cyanobacteria, could affect microbiota assembly and ecosystem function. Here, the trajectory of cyanopeptides was followed and linked to microbiota during Microcystis-dominated CyanoHABs in Lake Taihu, China. The most abundant cyanopeptide classes detected included microginin, spumigin, microcystin, nodularin and cyanopeptolin with total MC-LR-equivalent concentrations between 0.23 and 2051.54 ppb, of which cyanotoxins beyond microcystins (e.g., cyanostatin B and nodularin_R) far exceeded reported organismal IC₅₀ and negatively correlated with microbiota diversity, exerting potential collective eco-toxicities stronger than microcystins alone. The microbial communities were differentiated by size fraction and sampling date throughout CyanoHABs, and surprisingly, their variances were better explained by cyanopeptides (19-38%) than nutrients (0-16%). Cyanopeptides restriction (e.g., inhibition) and degradation were first quantitatively verified as the deterministic drivers governing community assembly, with stochastic processes being associated with the interplay between cyanopeptide dynamics and lake microbiota. This study presents an emerging paradigm in which cyanopeptides restriction and degradation co-mediate lake water microbiota assembly, unveiling new insights into the ecotoxicological significance of CyanoHABs to freshwater ecosystems.

Microbial Community Structure and Ecological Networks during Simulation of Diatom Sinking

Microbial-mediated utilization of particulate organic matter (POM) during its downward transport from the surface to the deep ocean constitutes a critical component of the global ocean carbon cycle. However, it remains unclear as to how high hydrostatic pressure (HHP) and low temperature (LT) with the sinking particles affects community structure and network interactions of the particle-attached microorganisms (PAM) and those free-living microorganisms (FLM) in the surrounding water. In this study, we investigated microbial succession and network interactions in experiments simulating POM sinking in the ocean. Diatom-derived ¹³C- and ¹²C-labeled POM were used to incubate surface water microbial communities from the East China Sea (ECS) under pressure (temperature) of 0.1 (25 °C), 20 (4 °C), and 40 (4 °C) MPa (megapascal). Our results show that the diversity and species richness of the PAM and FLM communities decreased significantly with HHP and LT. Microbial community analysis indicated an increase in the relative abundance of Bacteroidetes at high pressure (40 MPa), mostly at the expense of Gammaproteobacteria, Alphaproteobacteria, and Gracilibacteria at atmospheric pressure. Hydrostatic pressure and temperature affected lifestyle preferences between particle-attached (PA) and free-living (FL) microbes. Ecological network analysis showed that HHP and LT enhanced microbial network interactions and resulted in higher vulnerability to networks of the PAM communities and more resilience of those of the FLM communities. Most interestingly, the PAM communities occupied most of the module hubs of the networks, whereas the FLM communities mainly served as connectors of the modules, suggesting their different ecological roles of the two groups of microbes. These results provided novel insights into how HHP and LT affected microbial community dynamics, ecological networks during POM sinking, and the implications for carbon cycling in the ocean.

Pollution gradients shape the co-occurrence networks and interactions of sedimentary bacterial communities in Taihu Lake, a shallow eutrophic lake

The co-occurrence networks and interactions of bacterial communities in sediments are highly variable with environmental factors, which are vital to the nutrient biogeochemical cycle, pollutants biodegradation, and microbial community stability in lake ecosystems. Although pollution gradients reflect environmental variation comprehensively, few studies have characterized the changes in co-occurrence networks and interactions of bacterial communities along sediment pollution gradients. In order to investigate the impact of pollution gradients on compositions, co-occurrence networks, and interactions of sedimentary microbial communities, we studied the bacterial communities in the sediments of a typical shallow eutrophic lake, Taihu Lake, along pollution gradients using 16S rRNA gene high-throughput sequencing technology. All the sediment sampling sites were classified into mild, moderate, and severe pollution groups according to the sediments' physicochemical properties. Our results showed that the taxon richness was lowest in the severe pollution group, and the diversity of species decreased with the level of pollution. The complexity of the co-occurrence network decreased as the level of pollution increased, and the severe pollution group was characterized by a small-world network. The relative abundance of Proteobacteria, Bacteroidetes, and Chlorobi increased significantly as the level of pollution increased (P < 0.05). Strong inter-phyla co-occurrence or co-exclusion patterns demonstrated that the strength of interactions was enhanced in the severe pollution group, indicating stronger cooperative or competitive relationships. Chloroflexales and Chlorobiales were unique keystone taxa in the severe pollution group. The results of this study indicate that severe pollution reduces microbial diversity and network complexity, which may lead to community instability. The competition for nutrients of some copiotrophic bacteria may be enhanced as the level of pollution increased. The unique keystone taxa may contribute to photosynthesis and pollutant degradation in the severe pollution group. These findings expand our understanding of variation in bacterial co-occurrence networks and interactions along sediment pollution gradients.

Genetic detection of freshwater harmful algal blooms: A review focused on the use of environmental DNA (eDNA) in Microcystis aeruginosa and Prymnesium parvum

Recurrence and severity of harmful algal blooms (HABs) are increasing due to a number of factors, including human practices and climate change. Sensitive and robust methods that allow for early and expedited HAB detection across large landscape scales are needed. Among the suite of HAB detection tools available, a powerful option exists in genetics-based approaches utilizing environmental sampling, also termed environmental DNA (eDNA). Here we provide a detailed methodological review of three HAB eDNA approaches (quantitative PCR, high throughput sequencing, and isothermal amplification). We then summarize and synthesize recently published eDNA applications covering a variety of HAB surveillance and research objectives, all with a specific emphasis in the detection of two widely problematic freshwater species, Microcystis aeruginosa and Prymnesium parvum. In our summary and conclusion we build on this literature by discussing ways in which eDNA methods could be advanced to improve HAB detection. We also discuss ways in which eDNA data could be used to potentially provide novel insight into the ecology, mitigation, and prediction of HABs.

Diversity Distribution, Driving Factors and Assembly Mechanisms of Free-Living and Particle-Associated Bacterial Communities at a Subtropical Marginal Sea

Free-living (FL) and particle-associated (PA) bacterioplankton communities play critical roles in biogeochemical cycles in the ocean. However, their community composition, assembly process and functions in the continental shelf and slope regions are poorly understood. Based on 16S rRNA gene amplicon sequencing, we investigated bacterial communities’ driving factors, assembly processes and functional potentials at a subtropical marginal sea. The bacterioplankton community showed specific distribution patterns with respect to lifestyle (free living vs. particle associated), habitat (slope vs. shelf) and depth (surface vs. DCM and Bottom). Salinity and water temperature were the key factors modulating turnover in the FL community, whereas nitrite, silicate and phosphate were the key factors for the PA community. Model analyses revealed that stochastic processes outweighed deterministic processes and had stronger influences on PA than FL. Homogeneous selection (Hos) was more responsible for the assembly and turnover of FL, while drift and dispersal limitation contributed more to the assembly of PA. Importantly, the primary contributor to Hos in PA was Gammaproteobacteria:Others, whereas that in FL was Cyanobacteria:Bin6. Finally, the PICRUSt2 analysis indicated that the potential metabolisms of carbohydrates, cofactors, amino acids, terpenoids, polyketides, lipids and antibiotic resistance were markedly enriched in PA than FL.

Effects of shading levels on the composition and co-occurrence patterns of bacterioplankton and epibiotic bacterial communities of Cabomba caroliniana

Epibiotic bacterial community colonized on the plant leaf plays important roles in promoting plant growth and nutrient absorption, but is sensitive to environmental changes. As one of the most important environmental factors affecting the growth of plants and photosynthetic microorganisms, light may affect the diversity, composition, and interactions of the epibiotic bacterial community. Submerged plants in the aquatic ecosystem may be more sensitive to light intensity variations compared to the terrestrial plants since they usually receive less light. However, the effects of light on the interactions between the submerged plants and their epibiotic microbial communities remain uncertain. Here we used the 16S rRNA gene high-throughput sequencing to investigate the diversity and composition of the bacterioplankton and epibiotic bacterial communities of the Cabomba caroliniana under four different shading levels. A total of 24 water and leaf samples were collected from the experimental microcosms near Lake Taihu. We found the epibiotic bacterial community possessed a higher diversity than that of the bacterioplankton community, although the alpha diversity of the bacterioplankton community was more susceptible to different levels of shading. SourceTracker analysis revealed that with the increase of shading, the colonization of bacterioplankton to epibiotic bacteria decreased. Network analysis showed that the bacterial community network at 50% shading level had the lowest modularity and highest clustering coefficient compared to the bacterial community networks of other shading levels. Our findings provided new understandings of the effects of different light intensities on the epibiotic bacterial communities of submerged macrophytes.

The Lifestyle-Dependent Microbial Interactions Vary Between Upstream and Downstream of the Three Gorges Dam

Dams represent the most significant anthropogenic disturbance to global rivers. Previous studies have shown that free-living and particle-attached microbes exhibited differentially in river and reservoir ecosystems. However, little is known about the dam’s effect on their co-occurrence patterns. Here, a random matrix theory (RMT)-based network approach was used to construct microbial ecological networks for free-living and particle-attached communities in the immediate vicinity of the Three Gorges Dam (TGD), based on a high-throughput sequencing of 16S rRNA gene. Microbial distribution pattern showed that differences caused by lifestyle (free-living vs. particle-attached) were greater than those caused by geographic position (upstream vs. downstream of the TGD). Network analysis revealed higher connectivity and a lower number of modules in the overall downstream networks. Furthermore, considering the lifestyle, the network structures and properties for free-living and particle-attached microbes were different between upstream and downstream of the dam. Specifically, free-living communities located upstream of the dam exhibited a more complex co-occurrence pattern than the particle-attached communities, whereas the opposite was true for those located downstream of the dam. This variation indicated a strong impact of the dam on microbial interactions for microbes with similar lifestyle in the vicinity of the dam. We identified 112 persistent operational taxonomic unit (OTU)-level species that stably coexisted regardless of lifestyle and geographic positions. These persistent species occupied 21.33–25.57% of the total nodes in each network, and together with their first neighbors, they contributed more than 50% of the nodes and edges belonging to each network. Furthermore, we found that taxonomic affiliations for central nodes (with high degree) varied in these persistent species sub-networks. Collectively, our findings expand the current understanding of the dam’s effect on species interaction variation patterns for free-living and particle-attached communities in the vicinity of the dam, which are more complex than traditional alpha and beta microbial diversity.

Disentangling the seasonal co-occurrence patterns and ecological stochasticity of planktonic and benthic bacterial communities within multiple lakes

Both the planktonic bacterial community (PBC) and benthic bacterial community (BBC) are important for biogeochemical processes in freshwater lakes. Despite their ecological significance, little is known about their seasonal co-occurrence patterns and the ecological processes that drive them. In this study, we aimed to investigate the ecological associations among bacterial taxa and assembly processes of PBC and BBC in different seasons. We used 16S rRNA gene high-throughput sequencing of a total of 150 water and sediment samples collected from multiple lakes distributed in an urban region of China during winter and summer. Our results revealed that PBC showed stronger seasonal variations in co-occurrence patterns than BBC, suggesting that BBC had greater temporal stability than PBC. Winter PBC network was characterized by higher connectivity and complexity, and thereby the formation of a highly stable community structure; whereas lower connectivity arising from the presence of fewer predicted keystone taxa (hubs and connectors in a network) was destabilizing to summer PBC network. In addition, the phylum Firmicutes identified as a putative keystone taxon of PBC in both seasons played a non-negligible role in maintaining network structure which may result from strong functional associations with other bacterioplankton. Temperature and pH were the best explanatory factors predicting the seasonal co-occurrence patterns of PBC and BBC, respectively. Normalized stochasticity ratio based on null-model analysis indicated that deterministic processes overwhelmed stochastic processes in governing the assembly of PBC and BBC in both seasons. However, we observed a greater influence of ecological stochasticity on BBC assembly than PBC assembly in both seasons. Taken together, these findings provide insights into understanding the impacts of habitat heterogeneity and seasonal variability on microbial assemblage patterns in lake ecosystems.

The Relationships Between the Free-Living and Particle-Attached Bacterial Communities in Response to Elevated Eutrophication

Exploring the relationships between free-living (FL) and particle-attached (PA) bacterial communities can provide insight into their connectivity and the partitioning of biogeochemical processes, which is crucial to understanding the elemental cycles and metabolic pathways in aquatic ecosystems. However, there is still intense debate about that whether FL and PA fractions have the same assemblage. To address this issue, we investigated the extent of similarity between FL and PA bacterial communities along the environmental gradients in Lake Wuli, China. Our results revealed that the west Lake Wuli was slightly eutrophic and the east lake was moderately and highly eutrophic. The alpha-diversity of the FL bacterial communities was significantly lower than that of the PA fraction in the west lake, whereas the alpha-diversity of the two fractions was comparable in the east lake. The beta-diversity of both communities significantly differed in the west lake, whereas it resembled that in the east lake. Moreover, functional prediction analysis highlighted the significantly larger differences of metabolic functions between the FL and PA fractions in the west lake than in the east lake. Suspended particles and carbon resource promote the similarity between the FL and PA fractions. Collectively, our result reveals a convergent succession of aquatic communities along the eutrophic gradient, highlighting that the connectivity between FL and PA bacterial communities is nutrient related.

Illumina MiSeq sequencing and network analysis the distribution and co-occurrence of bacterioplankton in Danjiangkou Reservoir, China

Network analysis has contributed to studies of the interactions of microorganisms and the identification of key populations. However, such analysis has rarely been conducted in the study of reservoir bacterioplankton communities. This study investigated the bacterioplankton community composition in the surface water of the Danjiangkou Reservoir using the Illumina MiSeq sequencing platform. We observed that the bacterioplankton community primarily consisted of 27 phyla and 336 genera, including Actinobacteria, Proteobacteria, and Bacteroidetes, demonstrating the richness of the community composition. Redundancy analysis of the bacterioplankton communities and environmental variables showed that the total nitrogen (TN), pH, chemical oxygen demand (COD), and permanganate index (COD_Mn) were important factors affecting the bacterioplankton distribution. Network analysis was performed using the relative abundances of bacterioplankton based on the phylogenetic molecular ecological network (pMEN) method. The connectivity of node i within modules (Zi), the connectivity of node i among modules (Pi), and the number of key bacteria were high at the Taizishan and Heijizui sites, which were associated with higher TN contents than at the other sites. Among the physicochemical properties of water, TN, ammonia nitrogen (NH₄-N), pH, COD, and dissolved oxygen (DO) might have great influences on the functional units of the bacterial communities in bacterioplankton molecular networks. This study improves the understanding of the structure and function of bacterioplankton communities in the Danjiangkou Reservoir.

The global Microcystis interactome

Bacteria play key roles in the function and diversity of aquatic systems, but aside from study of specific bloom systems, little is known about the diversity or biogeography of bacteria associated with harmful cyanobacterial blooms (cyanoHABs). CyanoHAB species are known to shape bacterial community composition and to rely on functions provided by the associated bacteria, leading to the hypothesized cyanoHAB interactome, a coevolved community of synergistic and interacting bacteria species, each necessary for the success of the others. Here, we surveyed the microbiome associated with Microcystis aeruginosa during blooms in 12 lakes spanning four continents as an initial test of the hypothesized Microcystis interactome. We predicted that microbiome composition and functional potential would be similar across blooms globally. Our results, as revealed by 16S rRNA sequence similarity, indicate that M. aeruginosa is cosmopolitan in lakes across a 280° longitudinal and 90° latitudinal gradient. The microbiome communities were represented by a wide range of operational taxonomic units and relative abundances. Highly abundant taxa were more related and shared across most sites and did not vary with geographic distance, thus, like Microcystis, revealing no evidence for dispersal limitation. High phylogenetic relatedness, both within and across lakes, indicates that microbiome bacteria with similar functional potential were associated with all blooms. While Microcystis and the microbiome bacteria shared many genes, whole-community metagenomic analysis revealed a suite of biochemical pathways that could be considered complementary. Our results demonstrate a high degree of similarity across global Microcystis blooms, thereby providing initial support for the hypothesized Microcystis interactome.

Network Analysis Reveals Seasonal Patterns of Bacterial Community Networks in Lake Taihu under Aquaculture Conditions

Bacterial communities play essential roles in multiple ecological processes, such as primary production and nutrient recycling in aquatic systems. However, although the composition, diversity and function of bacterial communities have been well studied, little is known about the interactions and co-occurrence characteristics of these communities, let alone their seasonal patterns. To investigate the seasonal variations of bacterial community interactions, we collected water samples from four seasons in Lake Taihu and applied network analysis to reveal bacterial community interactions. Bacterial community networks were non-random in structure, and interactions among bacterial taxa in the networks varied markedly in different seasons. The autumnal bacterial network was the largest and most complex among obtained networks, whereas the spring correlation network was the simplest, having no module hubs or connectors. The important species of the networks were the dominant bacterial phyla/classes (e.g., Alphaproteobacteria and Bacteroidetes), although their relative abundance varied among seasons. The relationships between species and measured environmental variables changed over seasons; fewer environmental factors were correlated with bacterial species in the spring bacterial network, while we observed a greater number of species-environment correlations in the winter network. Our study highlights the seasonal differences in bacterial community interactions and expands our understanding of freshwater microbial ecology in systems affected by aquaculture.

Heterogeneity of interactions of microbial communities in regions of Taihu Lake with different nutrient loadings: A network analysis

To investigate the differences in the interactions of microbial communities in two regions in Taihu Lake with different nutrient loadings [Meiliang Bay (MLB) and Xukou Bay (XKB)], water samples were collected and both intra- and inter-kingdom microbial community interactions were examined with network analysis. It is demonstrated that all of the bacterioplankton, microeukaryotes and inter-kingdom communities networks in Taihu Lake were non-random. For the networks of bacterioplankton and inter-kingdom community in XKB, higher clustering coefficient and average degree but lower average path length indexes were observed, indicating the nodes in XKB were more clustered and closely connected with plenty edges than those of MLB. The bacterioplankton and inter-kingdom networks were considerably larger and more complex with more module hubs and connectors in XKB compared with those of MLB, whereas the microeukaryotes networks were comparable and had no module hubs or connectors in the two lake zones. The phyla of Acidobacteria, Cyanobacteria and Planctomycetes maintained greater cooperation with other phyla in XKB, rather than competition. The relationships between microbial communities and environmental factors in MLB were weaker. Compared with the microbial community networks of XKB, less modules in networks of MLB were significantly correlated with total phosphorous and total nitrogen.