A cohesive Microcoleus 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 <i>Microcoleus</i> sequence variant matching the <i>Microcoleus anatoxicus</i> 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 <i>Microcoleus</i> 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 <i>M. anatoxicus</i> be monitored internationally in order to help predict and mitigate similar cyanotoxic events.