Temporal heterogeneity of bacterial communities and their responses to Raphidiopsis raciborskii blooms

To elucidate the interspecies connectivity between cyanobacteria and other bacteria (noncyanobacteria), microbial diversity and composition were investigated through high-throughput sequencing (HTS) in a drinking water reservoir in Chongqing city, Southwest China, during Raphidiopsis raciborskii blooms. Significant temporal changes were observed in microbial community composition during the sampling period, primarily reflected by variations in relative bacterial abundance. The modularity analysis of the network demonstrated that the bacterial community forms co-occurrence/exclusion patterns in response to variations in environmental factors. Moreover, five modules involved in the dynamic phases of the R. raciborskii bloom were categorized into the Pre-Bloom, Bloom, Post-Bloom, and Non-Bloom Groups. The reservoir was eutrophic (i.e., the average concentrations of total nitrogen (TN) and total phosphorus (TP) were 2.32 and 0.07 mg L^-1, respectively) during the investigation; however, Pearson's correlation coefficient showed that R. raciborskii was not significantly correlated with nitrogen and phosphorus. However, other environmental factors, such as water temperature, pH, and the permanganate index, were positively correlated with R. raciborskii. Importantly, Proteobacteria (α-, γ-Proteobacteria), Acidobacteria, Chloroflexi, and Firmicutes were preferentially associated with increased R. raciborskii blooms. These results suggested that the transition of R. raciborskii bloom-related microbial modules and their keystone species could be crucial in the development and collapse of R. raciborskii blooms and could provide a fundamental basis for understanding the linkage between the structure and function of the microbial community during bloom dynamics.