Applicability of massively parallel sequencing on monitoring harmful algae at Varna Bay in the Black Sea

In this study the plankton diversity in 13 environmental samples from Varna Bay (in the western Black Sea) was analyzed using massively parallel sequencing (MPS). This preliminary study was undertaken to assess the potential of this technology for future implementation in monitoring programs in the Black Sea. Amplicon sequences of the 18S rRNA gene (V4-5 regions) were obtained using the Illumina MiSeq 250PE platform. A total of 1137 operational taxonomic units (OTUs) were obtained among which 242 OTUs with >0.990 BLAST top hit similarity (21.3% of all detected OTUs) closely related to sequences belonging to −protists. A large portion (175 OTUs=72.3%) was identified at the species levels, including species typical for the Bulgarian Black Sea plankton community, as well as many that haven’t been reported earlier in the Bulgarian Black Sea coast (124 OTUs=51.2%). Dinoflagellates were represented by the highest species number (77 OTUs comprising 31.8% of protist species), with dominant genera Gyrodinium and Heterocapsa. The present survey revealed the presence of 12 species listed as harmful, some of which have been previously overlooked, such as Cochlodinium polykrikoides, Karenia bicuneiformis, and Karlodinium veneficum. Species identification was possible for 10.3–36.0% of the detected OTUs in the six major supergroups. The frequency in Rhizaria was significantly lower than that in other major groups (p