Multilocus sequence typing (MLST) reveals high genetic diversity and clonal population structure of the toxic cyanobacterium Microcystis aeruginosa

Microcystis aeruginosa is one of the most prevalent bloom-forming cyanobacteria and has been the cause of increasing public health concern due to the production of hepatotoxins (microcystins). To investigate the genetic diversity, clonality and evolutionary genetic background with regard to the toxicity of M. aeruginosa, a multilocus sequence typing (MLST) scheme was developed, based on seven selected housekeeping loci (ftsZ, glnA, gltX, gyrB, pgi, recA and tpi). Analysis of a collection of 164 isolates from Japan and other countries identified 79 unique sequence types (STs), revealing a high level of genetic diversity (H=0.951). Although recombination between loci was indicated to be substantial by Shimodaira-Hasegawa (SH) tests, multilocus linkage disequilibrium analyses indicated that recombination between strains probably occurs at some frequency but not to the extent at which alleles are associated randomly, suggesting that the population structure of M. aeruginosa is clonal. Analysis of subsets of strains also indicated that the clonal population structure is maintained even in a local population. Phylogenetic analysis based on the concatenated sequences of seven MLST loci demonstrated that microcystin-producing genotypes are not monophyletic, providing further evidence for the gain and loss of toxicity during the intraspecific diversification of M. aeruginosa. However, toxic strains are genetically distinct from non-toxic strains in MLST allelic profiles, and it was also shown that non-toxic strains harbouring toxin genes fall into a single monophyletic clade, except for one case. These results suggest that the toxicity of M. aeruginosa is relatively stable in the short term, and therefore can be unequivocally characterized by MLST. The MLST scheme established here will be of great help for future detailed population genetic studies of M. aeruginosa.