Succession of diversity, assembly mechanisms, and activities of the microeukaryotic community throughout Scrippsiella acuminata (Dinophyceae) bloom phases

Harmful algal bloom (HAB) is a rapidly expanding marine ecological hazard. Although numerous studies have been carried out about the ecological impact and the ecological mechanism of HAB outbreaks, few studies have comprehensively addressed the shifts of species composition, metabolic activity level, driving factors and community assembly mechanisms of microeukaryotic plankton in the course of the bloom event. To fill the gap of research, we conducted 18S ribosomal DNA and RNA sequencing during the initiation, development, sustenance and decline stages of a Scrippsiella acuminata (S. acuminata) bloom at the coastal sea of Fujian Province, China. We found that the bloom event caused a decrease in microeukaryotic plankton species diversity and increase in community homogeneity. Our results revealed that the RNA- and DNA-inferred communities were similar, but α-diversity was more dynamic in RNA- than in DNA-inferred communities. The main taxa with high projected metabolic activity (with RNA:DNA ratio as the proxy) during the bloom included dinoflagellates, Cercozoa, Chlorophyta, Protalveolata, and diatoms. The role of deterministic processes in microeukaryotic plankton community assembly increased during the bloom, but stochastic processes were always the dominant assembly mechanism throughout the bloom process. Our findings improve the understanding of temporal patterns, driving factors and assembly mechanisms underlying the microeukarytic plankton community in a dinoflagellate bloom.

A comprehensive insight into the abundance and community of anammox bacteria in sediments of Hangzhou Bay, China

A total of 13 surface sediments were collected from Hangzhou Bay (HZB) for an investigation into the distribution and influencing factors of anammox bacterial community. The anammox bacterial 16S rRNA and hzo genes ranged between 2.34 × 105 to 9.22 × 105 copies/g and 3.68 × 105 to 1.70 × 106 copies/g, respectively. The results of high throughput sequencing (HTS) revealed that the obtained OTUs were affiliated with five known genera, named Ca. Scalindua, Ca. Jettenia, Ca. Brocadia, Ca. Kuenenia and Ca. Anammoxoglobus. RDA analysis indicated that salinity, pH, and water depth influenced the anammox bacterial community. Furthermore, network analysis identified Ca. Scalindua as a key genus. Neutral community model (NCM) and modified stochasticity ratio (MST) indicated that the deterministic process dominated the anammox bacterial community assembly. Overall, this study offers a more comprehensive understanding of the abundance and community of anammox bacteria in the sediments of HZB.

Microbial one‑carbon and nitrogen metabolisms are beneficial to the reservoir recovery after cyanobacterial bloom

Cyanobacterial blooms have profound effects on the structure and function of plankton communities in inland waters, but few studies have focused on the effects of microbial-based processes in one‑carbon and nitrogen cycling on water quality improvement following the bloom. Here, we compared the structure and function of the bacterial community, focusing on microbial one‑carbon and nitrogen metabolisms during and after a cyanobacterial Microcystis bloom in a deep subtropical reservoir. Our data showed that microbial one‑carbon and nitrogen cycles were closely related to different periods of the bloom, and the changes of functional genes in microbial carbon and nitrogen cycling showed the same consistent trend as that of Methylomonas sp. With the receding of the bloom, the abundance of Methylomonas as well as the functional genes of microbial one‑carbon and nitrogen cycling reached the peak and then recovered. Our results indicate that microbial one‑carbon and nitrogen metabolisms were beneficial to the recovery of water quality from the cyanobacterial bloom. This study lays a foundation for a deep understanding of the cyanobacterial decomposition mediated by microbes in one‑carbon and nitrogen cycles in inland freshwaters.

Rare biosphere regulates the planktonic and sedimentary bacteria by disparate ecological processes in a large source water reservoir

The bacteria in the water column and surface sediments are inherently intertwined and inseparable in aquatic ecosystems, yet little is known about the integrated spatiotemporal dynamics and driving mechanisms of both planktonic and sedimentary bacterial communities in reservoirs. By investigating the planktonic and sedimentary bacteria during four seasons from 88 samples of 11 representative sites across the Danjiangkou reservoir, we depicted an integrated biogeographic pattern of bacterial communities in the water source of the world's largest water diversion project. Our study revealed both planktonic (mantel r = 0.502, P = 0.001) and sedimentary (mantel r = 0.131, P = 0.009) bacterial communities were significantly correlated with environmental heterogeneity, but a weak disparity along spatial heterogeneity, and the significant seasonal dynamics of planktonic (mantel r = 0.499, P = 0.001) rather than sedimentary bacteria. Particularly, rare biosphere played a main role in determining the community succession in the reservoir. It not only exhibited a more striking environmental separation than abundant taxa but also was an essential part in mediating spatiotemporal shifts of planktonic bacteria and maintaining the stability of bacterial community. These rare bacteria were respectively mediated by stochastic (62.68%) and selective (79.60%) processes in water and sediments despite abundant taxa being largely determined by stochastic processes (86.88-93.96%). Overall, our study not only fills a gap in understanding the bacterial community dynamics and underlying drivers in source water reservoirs, but also highlights the particular importance of rare bacteria in mediating biogeochemical cycles in world's large reservoir ecosystems.

Antibiotic-accelerated cyanobacterial growth and aquatic community succession towards the formation of cyanobacterial bloom in eutrophic lake water

Antibiotics can stimulate the growth of model cyanobacterial species under pure culture conditions, but their influence on cyanobacterial blooms in natural aquatic ecosystems remains unclear. In this study, three commonly detected antibiotics (sulfamethoxazole, tetracycline, and ciprofloxacin) and their ternary mixture were proved to selectively stimulate (p < 0.05) the growth and photosynthetic activity of cyanobacteria in an aquatic microcosm at an environmentally relevant exposure dose of 300 ng/L under both oligotrophic and eutrophic conditions. Under the eutrophic condition, cyanobacteria reached a bloom density of 1.61 × 10⁶ cells/mL in 15 days without antibiotics, while the cyanobacteria exposed to tetracycline, sulfamethoxazole, ciprofloxacin, and their ternary mixture exceeded this bloom density within only 10, 8, 7, and 6 days, respectively. Principal coordinate analysis indicated that the antibiotic contaminants accelerated the prokaryotic community succession towards the formation of a cyanobacterial bloom by promoting the dominance of Microcystis, Synechococcus, and Oscillatoria under the eutrophic condition. After 15 days of culture, the antibiotic exposure increased the density of cyanobacteria by 1.38-2.31-fold and 2.28-3.94-fold under eutrophic and oligotrophic conditions, respectively. Antibiotic exposure generated higher stimulatory effects on cyanobacterial growth under the oligotrophic condition, but the antibiotic(s)-treated cyanobacteria did not form a bloom due to nutrient limitation. Redundancy analysis indicated that the three target antibiotics and their ternary mixture affected the prokaryotic community structure in a similar manner, while tetracycline showed some differences compared to sulfamethoxazole, ciprofloxacin, and the ternary antibiotic mixture with regard to the regulation of the eukaryotic community structure. This study demonstrates that antibiotic contaminants accelerate the formation of cyanobacterial blooms in eutrophic lake water and provides insights into the ecological effects of antibiotics on aquatic microbial communities.

Nitrogen cycling processes and the role of multi-trophic microbiota in dam-induced river-reservoir systems

The nitrogen (N) cycle is one of the most important nutrient cycles in river systems, and it plays an important role in maintaining biogeochemical balance and global climate stability. One of the main ways that humans have altered riverine ecosystems is through the construction of hydropower dams, which have major effects on biogeochemical cycles. Most previous studies examining the effects of damming on N cycling have focused on the whole budget or flux along rivers, and the role of river as N sources or sinks at the global or catchment scale. However, so far there is still lack of comprehensive and systematic summarize on N cycling and the controlling mechanisms in reservoirs affected by dam impoundment. In this review, we firstly summarize N cycling processes along the longitudinal riverine-transition-lacustrine gradient and the vertically stratified epilimnion-thermocline-hypolimnion gradient. Specifically, we highlight the direct and indirect roles of multi-trophic microbiota and their interactions in N cycling and discuss the main factors controlling these biotic processes. In addition, future research directions and challenges in incorporating multi-trophic levels in bioassessment, environmental flow design, as well as reservoir regulation and restoration are summarized. This review will aid future studies of N fluxes along dammed rivers and provide an essential reference for reservoir management to meet ecological needs.

Transformations from specialists to generalists cause bacterial communities are more stable than micro-eukaryotic communities under anthropogenic activity disturbance

Different microbial components have different responses to environmental disturbances. Here, we found that the planktonic bacterial and micro-eukaryotic communities had different responses to anthropogenic activity disturbance in a subtropical river, because they had different survival strategies (generalist and specialist). We used nutrients (nitrogen and phosphorus) as indicators of anthropogenic activities. We found that river stretch 1 showed low nutrient concentrations from October 2018 to September 2019. However, a nutrient disturbance was observed in river stretch 2. The nutrient concentrations increased largely in December and January but recovered to low values in June. Bacterial communities had higher resilience under this disturbance than micro-eukaryotic communities in river stretch 2. The bacterial community composition were quite different between the two river stretches in December and January but were similar in June and July. However, the differences of micro-eukaryotic community composition between the two river stretches were always high during the study period. The bacterial communities in river stretch 2 contained more generalists and nutrient tolerant specialists. The bacterial nutrient tolerant specialists rapidly decreased in the low nutrient months and were replaced by the generalists. Bacteria which were involved in this shifts accounted for 29.3% of the total abundance. However, the micro-eukaryotic communities in river stretch 2 contained more moderate generalists. These moderate generalists were insensitive to the variation of nutrients and only 19.56% of the micro-eukaryotes had significant responses to the disturbance. The survival strategies caused bacterial communities had higher adaptability than eukaryotes to environmental fluctuation.

Underestimated methane production triggered by phytoplankton succession in river-reservoir systems: Evidence from a microcosm study

The impoundment of dammed rivers accelerates phytoplankton succession from river-dominated to lake-dominated species. Little is known about the role of phytoplankton succession in methane (CH₄) production. In this study, we performed a 61-day microcosm investigation to simulate the collapse processes of Cyclotella meneghiniana (river-dominated algae) and Chlorella pyrenoidosa and Microcystis aeruginosa (lake-dominated algae). The results suggested that different methanogenic conditions were induced by the collapse of river-and lake-dominated algae. The rapid settlement of C. meneghiniana induced aerobic conditions in the water that inhibited anaerobic CH₄ production and intensified CH₄ oxidation as a result of an increase in pmoA. However, the decomposition of C. pyrenoidosa and M. aeruginosa depleted dissolved oxygen and provided abundant labile organic matter, which jointly elevated mcrA and the mcrA/pmoA ratio. Under this condition, anaerobic CH₄ production was the dominant pathway for the mineralization of algae-derived carbon. Finally, the CH₄ produced per unit of particulate total carbon (identified as the carbon content of the algal biomass) by C. pyrenoidosa and M. aeruginosa was 16.29-fold and 8.56-fold higher, respectively, than that produced by C. meneghiniana. These observations provided evidence that lake-dominated algae played a more vital role in CH₄ production than river-dominated algae when algal succession occurred. This discovery might be a new and vital, yet largely underestimated CH₄ emission pathway in river-reservoir systems, that should be considered when evaluating the effect of hydraulic projects on greenhouse gas emissions.

Reservoir water stratification and mixing affects microbial community structure and functional community composition in a stratified drinking reservoir

To investigate how the aquatic bacterial community of a stratified reservoir drives the evolution of water parameters, the microbial community structure and network characteristics of bacteria in a stratified reservoir were investigated using Illumina MiSeq sequencing technology. A total of 42 phyla and 689 distinct genera were identified, which showed significant seasonal variation. Additionally, stratified variations in the bacterial community strongly reflected the vertical gradient and seasonal changes in water temperature, dissolved oxygen, and nutrition concentration. Furthermore, principal coordinate analysis indicated that most microorganisms were likely influenced by changes in water stratification conditions, exhibiting significant differences during the stratification period and mixing period based on Adonis, MRPP, and Anosim. Compared to the stratification period, 123 enhanced operational taxonomic units (OTUs; 29%) and 226 depleted OTUs (52%) were identified during the mixing period. Linear discriminant analysis effect size results showed that 15 major genera were enriched in the mixing period and 10 major genera were enriched in the stratification period. Importantly, network analysis revealed that the keystone species belonged to hgcI_clade, CL500-29, Acidibacter, Paucimonas, Flavobacterium, Prochlorothrix, Xanthomonadales, Chloroflexia, Burkholderiales, OPB56, KI89A_clade, Synechococcus, Caulobacter or were unclassified. Redundancy analysis showed that temperature, dissolved oxygen, pH, chlorophyll-α, total phosphorus, nitrate, and ammonia were important factors influencing the water bacterial community and function composition, which were consistent with the results of the Mantel test analysis. Furthermore, random forest analysis showed that temperature, dissolved oxygen, ammonia, and total dissolved phosphorous were the most important variables predicting water bacterial community and function community α- and β-diversity (P < 0.05). Overall, these results provide insight into the interactions between the microbial community and water quality evolution mechanism in Zhoucun reservoir.

In Vitro Toxicological Screening of Stable and Senescing Cultures of Aphanizomenon, Planktothrix, and Raphidiopsis

Toxicity of cyanobacteria is the subject of ongoing research, and a number of toxic metabolites have been described, their biosynthesis pathways have been elucidated, and the mechanism of their action has been established. However, several knowledge gaps still exist, e.g., some strains produce hitherto unknown toxic compounds, while the exact dynamics of exerted toxicity during cyanobacterial growth still requires further exploration. Therefore, the present study investigated the toxicity of extracts of nine freshwater strains of Aphanizomenon gracile, an Aphanizomenon sp. strain isolated from the Baltic Sea, a freshwater strain of Planktothrix agardhii, and two strains of Raphidiopsis raciborskii obtained from 25- and 70-day-old cultures. An in vitro experimental model based on Cyprinus carpio hepatocytes (oxidative stress markers, DNA fragmentation, and serine/threonine protein activity) and brain homogenate (cholinesterase activity) was employed. The studied extracts demonstrated toxicity to fish cells, and in general, all examined extracts altered at least one or more of considered parameters, indicating that they possess, to some degree, toxic potency. Although the time from which the extracts were obtained had a significant importance for the response of fish cells, we observed strong variability between the different strains and species. In some strains, extracts that originated from 25-day-old cultures triggered more harmful effects on fish cells compared to those obtained from 70-day-old cultures, whereas in other strains, we observed the opposite effect or a lack of a significant change. Our study revealed that there was no clear or common pattern regarding the degree of cyanobacterial bloom toxicity at a given stage of development. This means that young cyanobacterial blooms that are just forming can pose an equally toxic threat to aquatic vertebrates and ecosystem functioning as those that are stable or old with a tendency to collapse. This might be largely due to a high variability of strains in the bloom.

Dispersal limitation drives biogeographical patterns of anammox bacterial communities across the Yangtze River

Over the past few decades, anaerobic ammonium oxidation (anammox) has been extensively documented at different scales in natural ecosystems. Previous studies have stated that the community composition of anammox bacteria is shaped mainly by environmental factors, whereas spatial factors have been largely overlooked. This study investigated biogeographical patterns of anammox bacterial communities using 42 sediment samples along a 4300-km stretch of the Yangtze River, the longest river in Asia. A significant distance-decay relationship was observed for anammox bacterial community similarity, which was significantly influenced by mean dendritic distance rather than environmental factors. This implied that dispersal limitation plays an important role in shaping biogeographical pattern of anammox community. Furthermore, our results revealed that neutral processes played vital role in shaping community assembly of anammox bacteria, and their communities were seriously dispersal limited. These findings contrast with previous studies on community similarities between broad taxonomic groups, which are mainly determined by niche-based selection owing to greater niche distances within broad taxonomic groups than in anammox bacteria. Importantly, the slope of the distance-decay curve was much steeper than previously reported for whole bacteria, which indicating the species turnover rate of anammox bacterial community was significantly higher than that of the whole bacterial community. Anammox bacteria harbor strong adhesion ability and low dispersal potential, and ultimately exhibited a high species turnover rate. Together, in the context of biogeography, our results highlight the importance of dispersal limitation in shaping the biogeographical pattern of anammox bacterial community.

The impacts of different high-throughput profiling approaches on the understanding of bacterial antibiotic resistance genes in a freshwater reservoir

Bacterial antibiotic resistance genes (ARGs), a kind of emerging environmental pollutants, greatly threat human health through pathogenic bacteria. High-throughput quantitative PCR (HT-qPCR) and metagenomic approaches are two popular tools applied in aquatic environmental ARGs monitoring. However, current poor knowledge of different ARG profiling approaches' impacts on the understanding of the ARGs in aquatic environments greatly limit the further field application of these two approaches. For filling such knowledge gap, this study simultaneously employed these two approaches to examine and compare the ARGs in a freshwater reservoir across space and time. We found metagenomic approach detected more ARG subtypes and much higher bacitracin resistance genes' abundances than HT-qPCR. In general, HT-qPCR and metagenomics analyses both revealed similar ARG dynamic patterns and co-occurrence patterns between ARGs and bacterial taxa as well as the relationships between ARGs and environmental factors. Our results indicated the impacts of different ARG profiling approaches on the understanding of bacterial ARGs might be minor or negligible. HT-qPCR approach has the superiorities of time-saving, absolute quantification, low requirement for bioinformatics skills but also has some drawbacks including higher PCR amplification & primer bias, higher primer dependency and relative lower ARG subtype quantification capability compared to metagenomic approach. We suggest HT-qPCR approach can be employed for routine aquatic environmental monitoring, and metagenomic approach could be applied in comprehensive surveys for getting more ARG subtype information. Our data can be a useful reference for choosing right ARG profiling approaches for bacterial ARGs monitoring and risk assessment in aquatic environments.

Community dynamics of free-living and particle-attached bacteria following a reservoir Microcystis bloom

The composition of microbial communities can vary at the microspatial scale between free-living (FL) and particle-attached (PA) niches. However, it remains unclear how FL and PA bacterial communities respond to cyanobacterial blooms across water depths. Here, we examined the community dynamics of the FL (0.2-3 μm) and PA (>3 μm) bacterioplankton based on 16S rRNA gene high-throughput sequencing in a subtropical stratified reservoir under Microcystis aeruginosa bloom and non-bloom conditions. Both FL and PA bacterioplankton communities showed different responses in alpha- and beta-diversities to the bloom, suggesting the idea that the responses of bacterial community could depend on lifestyle. Specifically, abundant PA subcommunities showed a greater variation between bloom and non-bloom groups than abundant FL ones. In contrast, rare FL subcommunities exhibited a stronger response to water depth than rare PA ones. Furthermore, the rare taxa exhibited a preference for PA status, shaped and stimulated by the M. aeruginosa bloom. Our analyses also showed that PA bacterial communities were generally more diverse and appeared to be more responsive to routinely measured environmental variables than FL bacteria. Microcystis blooms had a facilitative influence on specific bacteria by mediating the transitions from free-living to particle-attached lifestyles. Altogether, these findings highlight the importance of bacterial lifestyle and abundance in understanding the dynamics of microbial community in cyanobacterial bloom aquatic ecosystem.

Metagenomics Response of Anaerobic Ammonium Oxidation (anammox) Bacteria to Bio-Refractory Humic Substances in Wastewater

Anammox-based processes have been widely applied for the treatment of wastewater (e.g., wastewater irrigation systems and constructed wetland) which consists of bio-refractory humic substances. Nonetheless, the impacts of bio-refractory humic substances on anammox consortia are rarely reported. In the present study, three identical lab-scale anammox reactors (i.e., HS0, HS1 and HS10), two of which were dosed with humic substances at 1 and 10 mg·L−1, respectively, were operated for nearly one year. The long-term operation of the reactors showed that the presence of humic substances in influent had no significant influence on nitrogen removal rates. Despite this, comparative metagenomics showed changes in anammox microbiota structure during the exposure to humic substance; e.g., the relative abundance of Candidatus Kuenenia was lower in HS10 (18.5%) than that in HS0 (22.8%) and HS1 (21.7%). More specifically, a lower level of humic substances (1 mg·L−1) in influent led to an increase of genes responsible for signal transduction, likely due to the role of humic substances as electron shuttles. In contrast, a high level of humic substances (10 mg·L−1) resulted in a slight decrease of functional genes associated with anammox metabolism. This may partially be due to the biodegradation of the humic substances. In addition, the lower dosage of humic substances (1 mg·L−1) also stimulated the abundance of hzs and hdh, which encode two important enzymes in anammox reaction. Overall, this study indicated that the anammox system could work stably over a long period under humic substances, and that the process was feasible for leachate treatment.

Characterization of bacterial and microbial eukaryotic communities associated with an ephemeral hypoxia event in Taihu Lake, a shallow eutrophic Chinese lake

While the important roles of microbial communities in oceanic hypoxic zones were beginning to be understood, little is known about microbial community associated with this phenomenon in shallow lakes. To address this deficit, both the bacterial and microbial eukaryotic communities of an ephemeral hypoxic area of Taihu Lake were characterized. The hypoxia provided nutritional niches for various bacteria, which results in high abundance and diversity. Specific bacterial groups, such as vadinBC27 subgroup of Bacteroidetes, Burkholderiales, Rhodocyclales, Pseudomonas, and Parcubacteria, were dominated in hypoxic sites and relevant to the fermentation, denitrification, nitrification, and sulfur metabolism. Conversely, most of microbial eukaryotes disappeared along with the decline of DO. An unexpected dominance of fungi was observed during hypoxia, which partly explained by the accumulation of toxic algae. Mucor was the single dominant genus in the hypoxic zone. We proposed that this group might cooperate with bacterial communities in the anaerobic degradation of algal biomass and woody materials. Generally, the hypoxic microbiome in shallow lakes is mainly involved in fermentative metabolism depending on phytodetritus and is potentially influenced by terrestrial sources. This study provided new insights into the unique microbiome in short-term hypoxia in shallow lakes and lays the foundation for studies that will enhance our understanding of the microbial players associated with hypoxia and their adaption strategy on the global scale.

The antibiotic resistome of free-living and particle-attached bacteria under a reservoir cyanobacterial bloom

In freshwater systems, both antibiotic resistance genes (ARGs) and cyanobacterial blooms attract global public health concern. Cyanobacterial blooms can greatly impact bacterial taxonomic communities, but very little is known about the influence of the blooms on antibiotic resistance functional community. In this study, the ARGs in both free-living (FL) and particle-attached (PA) bacteria under bloom and non-bloom conditions were simultaneously investigated in a subtropical reservoir using high-throughput approaches. In total, 145 ARGs and 9 mobile genetic elements (MGEs) were detected. The most diverse and dominant of which (68.93%) were multidrug resistance genes and efflux pump mechanism. The richness of ARGs in both FL and PA bacteria was significantly lower during the bloom period compared with non-bloom period. The abundance of ARGs in FL bacteria was significantly lower under bloom condition than in the non-bloom period, but the abundance of ARGs in PA bacteria stayed constant. More importantly, the resistant functional community in PA bacteria was more strongly influenced by the cyanobacterial bloom than in the FL bacteria, although >96% ARGs were shared in both FL and PA bacteria or both bloom and non-bloom periods. We also compared the community compositions between taxonomy and function, and found antibiotic resistant communities were highly variable and exhibited lower similarity between bloom and non-bloom periods than seen in the taxonomic composition, with an exception of FL bacteria. Altogether, cyanobacterial blooms appear to have stronger inhibitory effect on ARG abundance in FL bacteria, and stronger influence on antibiotic resistant community composition in PA bacteria. Our results further suggested that both neutral and selective processes interactively affected the ARG composition dynamics of the FL and PA bacteria. However, the antibiotic resistant community of FL bacteria exhibited a higher level of temporal stochasticity following the bloom event than PA bacteria. Therefore, we emphasized the bacterial lifestyles as an important mechanism, giving rise to different responses of antibiotic resistant community to the cyanobacterial bloom.

Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom

Plankton communities normally consist of few abundant and many rare species, yet little is known about the ecological role of rare planktonic eukaryotes. Here we used a 18S ribosomal DNA sequencing approach to investigate the dynamics of rare planktonic eukaryotes, and to explore the co-occurrence patterns of abundant and rare eukaryotic plankton in a subtropical reservoir following a cyanobacterial bloom event. Our results showed that the bloom event significantly altered the eukaryotic plankton community composition and rare plankton diversity without affecting the diversity of abundant plankton. The similarities of both abundant and rare eukaryotic plankton subcommunities significantly declined with the increase in time-lag, but stronger temporal turnover was observed in rare taxa. Further, species turnover of both subcommunities explained a higher percentage of the community variation than species richness. Both deterministic and stochastic processes significantly influenced eukaryotic plankton community assembly, and the stochastic pattern (e.g., ecological drift) was particularly pronounced for rare taxa. Co-occurrence network analysis revealed that keystone taxa mainly belonged to rare species, which may play fundamental roles in network persistence. Importantly, covariations between rare and non-rare taxa were predominantly positive, implying multispecies cooperation might contribute to the stability and resilience of the microbial community. Overall, these findings expand current understanding of the ecological mechanisms and microbial interactions underlying plankton dynamics in changing aquatic ecosystems.

Dynamics of Bacterial and Fungal Communities during the Outbreak and Decline of an Algal Bloom in a Drinking Water Reservoir

The microbial communities associated with algal blooms play a pivotal role in organic carbon, nitrogen and phosphorus cycling in freshwater ecosystems. However, there have been few studies focused on unveiling the dynamics of bacterial and fungal communities during the outbreak and decline of algal blooms in drinking water reservoirs. To address this issue, the compositions of bacterial and fungal communities were assessed in the Zhoucun drinking water reservoir using 16S rRNA and internal transcribed spacer (ITS) gene Illumina MiSeq sequencing techniques. The results showed the algal bloom was dominated by Synechococcus, Microcystis, and Prochlorothrix. The bloom was characterized by a steady decrease of total phosphorus (TP) from the outbreak to the decline period (p < 0.05) while Fe concentration increased sharply during the decline period (p < 0.05). The highest algal biomass and cell concentrations observed during the bloom were 51.7 mg/L and 1.9×10⁸ cell/L, respectively. The cell concentration was positively correlated with COD_Mn (r = 0.89, p = 0.02). Illumina Miseq sequencing showed that algal bloom altered the water bacterial and fungal community structure. During the bloom, the dominant bacterial genus were Acinetobacter sp., Limnobacter sp., Synechococcus sp., and Roseomonas sp. The relative size of the fungal community also changed with algal bloom and its composition mainly contained Ascomycota, Basidiomycota and Chytridiomycota. Heat map profiling indicated that algal bloom had a more consistent effect upon fungal communities at genus level. Redundancy analysis (RDA) also demonstrated that the structure of water bacterial communities was significantly correlated to conductivity and ammonia nitrogen. Meanwhile, water temperature, Fe and ammonia nitrogen drive the dynamics of water fungal communities. The results from this work suggested that water bacterial and fungal communities changed significantly during the outbreak and decline of algal bloom in Zhoucun drinking water reservoir. Our study highlights the potential role of microbial diversity as a driving force for the algal bloom and biogeochemical cycling of reservoir ecology.