Co-occurrence patterns and assembly processes of microeukaryotic communities in an early-spring diatom bloom

The interaction and assembly processes of microeukaryotic community compositions (MECs) are rarely elucidated in environment with strong disturbance such as harmful algal blooms. To fill this gap, we analyzed changes of MECs induced by a diatom bloom using 18S rRNA gene amplicon sequencing. The MECs were mainly dominated by Cercozoa (average relative abundance, 49.2%), Diatom (25.5%) and Dinoflagellata (15.6%). MECs changed significantly (ANOSIM P  < 0.01) in four-bloom stages. Environmental factors including pH, DO, nitrate and phosphate, together with bacterial communities could significantly influence the variation of MECs. Co-occurrence network analysis revealed a complex interaction between microeukaryotic and bacterial communities. Most OTUs in modules of the co-occurrence network were specific to one particular bloom stage. Phylogenetic based β-nearest taxon distance analyses revealed that stochastic processes mainly dominated microeukaryotic community assembly in the initial and after-bloom stage. However, microeukaryotic community assembly in middle and late stage of the bloom were driven by deterministic processes. In conclusion, both stochastic and deterministic processes play important roles in distinct bloom stages. These findings may expand current understandings of assembly mechanisms and microbial interactions underlying microeukaryotic dynamics in eutrophic aquatic ecosystems where harmful algal blooms occurred frequently.