Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale

Cyanobacteria are oxygenic photosynthetic bacteria that perform a substantial part of the global primary production. Some species are responsible for catastrophic environmental events, called blooms, which have become increasingly common in lakes and freshwater bodies as a consequence of global changes. Genotypic diversity is considered essential for marine cyanobacterial population, allowing it to cope with spatio-temporal environmental variations and to adapt to specific micro-niches in the ecosystem. This aspect is underestimated in the study of bloom development, however, and given little notice in studies of the ecology of harmful cyanobacteria. Here we compared the genomes of four strains of Aphanizomenon gracile, a species of filamentous toxinogenic cyanobacteria (Nostocales) found worldwide in fresh and brackish water. Millimeter-sized fascicles were isolated from a single water sample and have been maintained in culture since 2010. A comparative study revealed extensive heterogeneity in gene contents, despite similar genome size and high similarity indices. These variations were mainly associated with mobile genetic elements and biosynthetic gene clusters. For some of the latter, metabolomic analysis confirmed the production of related secondary metabolites, such as cyanotoxins and carotenoids, which are thought to play a fundamental role in the cyanobacterial fitness. Altogether, these results demonstrated that an A. gracile bloom could be a highly diverse population at low spatial scale and raised questions about potential exchanges of essential metabolites between individuals.