Responses of microeukaryotic community structure to a Phaeocystis globosa bloom in a semi-enclosed subtropical bay.

The occurrence of Phaeocystis globosa, a harmful algal bloom species in Chinese coastal waters, has significant impacts on marine organisms and poses a threat to the safety of coastal nuclear power plants. Although previous studies have established a close association between P. globosa blooms and the bacterial community, the relationship between the microeukaryotic community and P. globosa blooms remains poorly understood. In this study, the variations in the microeukaryotic community resulting from a P. globosa bloom were analyzed using 18S rRNA gene amplicon sequencing. The results indicated that the diversity of the microeukaryotic community during the bloom phase was significantly higher than that during the dissipation phase. The microeukaryotic community compositions varied significantly between the two phases of the P. globosa bloom. During the bloom phase, the dominant microeukaryotic was Viridiplantae, which was then replaced by Dinoflagellata during the dissipation phase. Co-occurrence network analysis showed that the relationship among the microeukaryotic community during the bloom phase was more complex than that during the dissipation phase, and the keystone taxa varied as the bloom progressed. Additionally, microeukaryotic community assembly was primarily driven by stochastic processes during the bloom phase based on the β-nearest taxon distance, whereas it was driven by both deterministic processes and stochastic processes during the dissipation phase. Overall, our findings provide novel insight into the mechanisms and interactions involved in microeukaryotic community dynamics in environments disturbed by P. globosa blooms.