Characterizing ecological relationships between viruses, bacteria and phytoplankton in the ocean is critical to understanding the ecosystem; however, these relationships are infrequently investigated together. To understand the dynamics of microbial communities and environmental factors in harmful algal blooms (HABs), we examined the environmental factors and microbial communities during Akashiwo sanguinea HABs in the Jangmok coastal waters of South Korea by metagenomics. Specific bacterial species showed complex synergistic and antagonistic relationships with the A. sanguinea bloom. The endoparasitic dinoflagellate Amoebophrya sp. 1 controlled the bloom dynamics and correlated with HAB decline. Among nucleocytoplasmic large DNA viruses (NCLDVs), two Pandoraviruses and six Phycodnaviruses were strongly and positively correlated with the HABs. Operational taxonomic units of microbial communities and environmental factors associated with A. sanguinea were visualized by network analysis: A. sanguinea–Amoebophrya sp. 1 (r = .59, time lag: 2 days) and A. sanguinea–Ectocarpus siliculosus virus 1 in Phycodnaviridae (0.50, 4 days) relationships showed close associations. The relationship between A. sanguinea and dissolved inorganic phosphorus relationship also showed a very close correlation (0.74, 0 day). Microbial communities and the environment changed dynamically during the A. sanguinea bloom, and the rapid turnover of microorganisms responded to ecological interactions. A. sanguinea bloom dramatically changes the environments by exuding dissolved carbohydrates via autotrophic processes, followed by changes in microbial communities involving host‐specific viruses, bacteria and parasitoids. Thus, the microbial communities in HAB are composed of various organisms that interact in a complex manner.

The total dinoflagellate cyst community and the cysts of Cochlodinium polykrikoides in the surface sediments of South Sea (Tongyeong coast), South Korea, were analysed using next-generation sequencing (NGS) and morphological approaches. Dinoflagellate cysts can be highly abundant (111–4,087 cysts g⁻¹ dry weight) and have diverse species composition. A total of 35 taxa of dinoflagellate cysts representing 16 genera, 21 species (including four unconfirmed species), and 14 complex species were identified by NGS analysis. Cysts of Scrippsiella spp (mostly Scrippsiella trochoidea) were the most dominant and Polykrikos schwartzii, Pentapharsodinium dalei, Ensiculifera carinata, and Alexandrium catenella/tamarense were common. Thus, a combination of NGS and morphological analysis is effective for studying the cyst communities present in a given environment. Although C. polykrikoides developed massive blooms during 2013–2014, microscopy revealed low density of their cysts, whereas no cysts were detected by NGS. However, the vegetative C. polykrikoides not appeared during 2015–2017 in spite of the observation of C. polykrikoides cysts. This suggests that the C. polykrikoides blooms were not due to development of their cysts but to other factors such as currents transporting them to a marine environment suitable for their growth.

We investigated the dynamics of the bacterial composition and metabolic function within Akashiwo sanguinea bloom using a 100-L indoor microcosm and metagenomic next-generation sequencing. We found that the bacterial community was classified into three groups at 54% similarity. Group I was associated with “during the A. sanguinea bloom stage” and mainly consisted of Alphaproteobacteria, Flavobacteriia and Gammaproteobacteria. Meanwhile, groups II and III were associated with the “late bloom/decline stage to post-bloom stage” with decreased Flavobacteriia and Gammaproteobacteria in these stages. Upon the termination of the A. sanguinea bloom, the concentrations of inorganic nutrients (particularly PO₄³⁻, NH₄⁺ and dissolved organic carbon) increased rapidly and then decreased. From the network analysis, we found that the A. sanguinea node is associated with certain bacteria. After the bloom, the specific increases in NH₄⁺ and PO₄³⁻ nodes are associated with other bacterial taxa. The changes in the functional groups of the bacterial community from chemoheterotrophy to nitrogen association metabolisms were consistent with the environmental impacts during and after A. sanguinea bloom. Consequently, certain bacterial communities and the environments dynamically changed during and after harmful algal blooms and a rapid turnover within the bacterial community and their function can respond to ecological interactions.