A cohesive strain cluster causes benthic cyanotoxic blooms in rivers worldwide.

Over the last two decades, proliferations of benthic cyanobacteria producing derivatives of anatoxin-a have been reported in rivers worldwide. Here, we follow up on such a toxigenic event happening in the Areuse river in Switzerland and investigate the diversity and genomics of major bloom-forming riverine benthic cyanobacteria. We show, using 16S rRNA-based community profiling, that benthic communities are dominated by Oscillatoriales. We correlate the detection of one sequence variant matching the species with the presence of anatoxin-a derivatives and use long-read metagenomics to assemble complete circular genomes of the strain. The main dihydro-anatoxin-a-producing strain in the Areuse is distinct from strains isolated in New Zealand, the USA, and Canada, but forms a monophyletic strain cluster with them with average nucleotide identity values close to the species threshold. Compared to the rest of the genus, the toxin-producing strains encode a 15 % smaller genome, lacking genes for the synthesis of some essential vitamins. Toxigenic mats harbor a distinct microbiome dominated by proteobacteria and bacteroidetes, which may support cyanobacterial growth by providing them with essential nutrients. We recommend that strains closely related to be monitored internationally in order to help predict and mitigate similar cyanotoxic events.