Jun R. Yang

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.

Although it is widely recognized that cyanobacterial blooms have substantial influence on the plankton community in general, their correlations with the whole community of eukaryotic plankton at longer time scales remain largely unknown. Here, we investigated the temporal dynamics of eukaryotic plankton communities in two subtropical reservoirs over a 6-year period (2010-2015) following one cyanobacterial biomass cycle-the cyanobacterial bloom (middle 2010), cyanobacteria decrease (late 2010-early 2011), non-bloom (2011-2014), cyanobacteria increase, and second bloom (late 2014-2015). The eukaryotic community succession that strongly correlated with this cyanobacterial biomass cycle was divided into four periods, and each period had distinct characteristics in cyanobacterial biomass and environments in both reservoirs. Integrated co-occurrence networks of eukaryotic plankton based on the whole study period revealed that the cyanobacterial biomass had remarkably high network centralities, and the eukaryotic OTUs that had stronger correlations with the cyanobacterial biomass exhibited higher centralities. The integrated networks were also modularly responded to different eukaryotic succession periods, and therefore correlated with the cyanobacterial biomass cycle. Moreover, sub-networks based on the different eukaryotic succession periods indicated that the eukaryotic co-occurrence patterns were not constant but varied largely associating with the cyanobacterial biomass. Based on these long-term observations, our results reveal that the cyanobacterial biomass cycle created distinct niches between persistent bloom, non-bloom, decrease and increase of cyanobacteria, and therefore associated with distinct eukaryotic plankton patterns. Our results have important implications for understanding how complex aquatic plankton communities respond to cyanobacterial blooms under the changing environments.