Characterization of Taxonomic and Functional Dynamics Associated with Harmful Algal Bloom Formation in Recreational Water Ecosystems

Harmful algal bloom (HAB) formation leads to the eutrophication of water ecosystems and may render recreational lakes unsuitable for human use. We evaluated the applicability and comparison of metabarcoding, metagenomics, qPCR, and ELISA-based methods for cyanobacteria/cyanotoxin detection in bloom and non-bloom sites for the Great Lakes region. DNA sequencing-based methods robustly identified differences between bloom and non-bloom samples (e.g., the relative prominence of Anabaena and Planktothrix). Shotgun sequencing strategies also identified the enrichment of metabolic genes typical of cyanobacteria in bloom samples, though toxin genes were not detected, suggesting deeper sequencing or PCR methods may be needed to detect low-abundance toxin genes. PCR and ELISA indicated microcystin levels and microcystin gene copies were significantly more abundant in bloom sites. However, not all bloom samples were positive for microcystin, possibly due to bloom development by non-toxin-producing species. Additionally, microcystin levels were significantly correlated (positively) with microcystin gene copy number but not with total cyanobacterial 16S gene copies. In summary, next-generation sequencing-based methods can identify specific taxonomic and functional targets, which can be used for absolute quantification methods (qPCR and ELISA) to augment conventional water monitoring strategies.