Intra

ISME Communications volume 3, Article number: 57 (2023) Cite this article

1 Altmetric

Metrics details

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.

The ecological success of cyanobacteria is partly a result of the intra-specific diversity of ecotypes [1,2,3], groups of individuals sharing close growth optima for environmental factors such as temperature and light intensity, which allow the population to cope with spatiotemporal variations in their habitat. Phylogenetic inferences based on single (e.g., Internal Transcribed Spacer, ITS) or multiple genes allow to group these ecotype members within a clade, suggesting the selection for adaptative capabilities. In terms of gene content, however, the ecotypes can also show significant divergence from each other. For instance the gene sets involved in the phosphorus assimilation, an essential abiotic factor for cyanobacterial growth, diverge between ecotypes from P-rich or P-limited aquatic ecosystems [4,5,6]. For the well-studied cyanobacterial species Prochlorochoccus marinus (Synechococcales) and Microcystis aeruginosa (Chroococcales), a single ecotype can also contain many strains with distinct genotypes. This so-called microdiversity relies primarily in these species on a high genomic plasticity [7], driven in part by mobile genetic elements [7,8,9] and is assumed to be involved in individuals adaptation to specific micro-niches [10].

Aphanizomenon gracile is a nitrogen-fixing filamentous cyanobacterium able to form fascicles and is responsible for toxic blooms in freshwater and brackish aquatic ecosystems. It belongs to the order, Nostocales, about which no data are available on genetic intra-population diversity. In this study, we sequenced the genomes of four strains of A. gracile isolated from a single sample of water to question the intra-population diversity at low-spatial scale.

We sequenced and assembled the genomes of four A. gracile strains (PMC627.10, PMC638.10, PMC644.10 and PMC649.10) firstly isolated from a single sample of water and maintained in monoculture (see Supplementary Methods). The resulting genomes displayed almost-perfect completeness (from 99.18% for PMC649.10 to 100% for PMC638.10), but were composed of numerous sequences, between 65 (PMC627.10) and 238 (PMC644.10), because of a large number of repeated sequences. The four genomes were very similar in size (5.40 ± 0.03 Mb), GC content (38.35 ± 0.05%) and number of tRNAs (from 40 to 42) (Table 1). Each strain shared with the others a minimal sequence similarity of 99.86% with the 16S rRNA coding gene, while the ITS sequences were all identical (Fig. S1). At the whole genome level, the proximity among strains was supported by pairwise ANI (Average Nucleotide Identity, Table 1) with values ≥99.12%. Such a level of similarity usually leads to the conclusion that related organisms belong to the same operational taxonomic unit and by extension to the same ecotype as the phylogenetically distant marine unicellular pico-cyanobacterium, Prochlorococcus marinus, in which microdiversity has been thoroughly studied [11]. In the absence of physiological data confirming the ecotypes of the A. gracile strains, these findings allow us to interpret the internal variability of the A. gracile strains in a microdiversity context.

Despite the absence of self-evident differences between genetic markers, comparative genomics revealed divergences among the strains. A striking representation of this is given by the low synteny index of 77.01 ± 1.31%, similar to those reported among 12 Microcystis aeruginosa strains sampled from distant geographical locations (76 ± 4%) [7]. As exemplified for PMC627.10, the variations of the synteny index along the genome were significantly associated with the COGs frequency distribution (Pearson's chi-squared test, p < 2.2e−16) meaning that the low synteny was due to the presence or absence of strain-specific genes rather than intra-genomic rearrangements (Fig. 1a). In total, among the 5911 COGs predicted in all genomes (4972 ± 28 per strain; Table 1), the core genome was represented by only 67.97% of them (4018) (Fig. 1b). The accessory gene fraction consisted of almost half (48.54%) singletons, i.e., strain-specific genes (between 192 and 305 by genome). The latter fraction was significantly enriched in genes encoding prokaryotic defense systems and transposases. These subcategories of "Replication, Recombination and Repair" (COG category L) as well as Orphans (t test H0: no difference, p < 0.01) are usually considered as markers of mobile genetic elements, here potentially involved in genomic plasticity in Aphanizomenon. Remarkably, COGs related to the biosynthesis of secondary metabolites (COG category Q) were also enriched in accessory fractions (Fig. 1c).

a Representation of the PMC627.10 A. gracile genome displaying from the inner to the outer circle: (i) the coding DNA sequence (CDS) on both strands with biosynthetic gene clusters (BGCs) in black, (ii) frequency (%) of CDS-associated clusters of orthologous genes (COGs) in the four A. gracile genomes, (iii) the synteny index, and (iv) the largest (>10kb) assembled scaffolds representing >99% of the genome. b Venn diagram of the distribution of COGs and singletons with white labels. c Percentage of COG functional categories with significant differences (Student's t test with p < 0.01) between core (white) and flexible (red) gene sets, including a focus on the "Replication, Recombination and Repair" COG category (L). COG categories: B-"Chromatin structure/dynamics", C-"Energy production/conversion", E-"Amino acid metabolism/transport", F-"Nucleotide metabolism/transport", G-"Carbohydrate metabolism/transport", H-"Coenzyme metabolism", I-"Lipid metabolism", J-"Translation", O-"Post-translational modification/protein turnover/chaperone functions", P-"Inorganic ion transport/metabolism", Q-"Secondary metabolism", and S-"Function unknown". d Presence or absence of BGCs (COG category Q) among A. gracile strains and closed Nostocales taxa gathered by BGC types (‘-‘, uncharacterized product). The solid circles indicate BGCs whose products have been formally identified by mass spectrometry. Left: Phylogenomic tree (see Supplementary Methods and Table S1 for the list of genes used). Right: histogram displaying the number of BGCs by strain colored by BGC type. NRPS/PKS non-ribosomal peptide synthetase/polyketide synthetase, RiPPs ribosomally synthesized, post-translationally modified peptides, terpenes and others (uncharacterized). The lists of BGC and analytes detected in each strain can be found in Supplementary Tables S2 and S3, respectively.

The study of A. gracile secondary metabolite biosynthetic gene clusters (BGC) revealed that their repertoires were unique to each A. gracile strain, with only 8 BGCs common to all strains out of a total of 22 distinct BGC characterized in these genomes (Fig. 1d). Although variable, the composition of this BGC set appeared to be clade-associated in the Nostocales tree, as the number of BGCs in common with A. gracile genomes decreased with increasing distance from their clade in the phylogenetic tree. Accessory BGCs were notably involved in the production of carotenoids and shinorines, both implicated in several adaptative processes including photoprotection [12, 13], and cyanotoxins (e.g. puwainaphycins and saxitoxins), whose production in the associated cultures was confirmed by high-resolution mass spectrometry (Fig. S2, Table S3). Although not-well characterized, cyanotoxins are thought to play key roles in adaptation, given the high cost of BGC maintenance in genomes and resource requirements for their production. For example, saxitoxins, which were only synthesized in our study by PMC 627.10 and PMC 638.10 strains from synthesized from their 26.5kb long BGC, are presumed to regulate Na+ homeostasis and have been proved to be actively extracellularly exported in salinity stress conditions [14].

Lastly, the group of accessory genes with resolutive functional annotations at the gene name level, revealed potential discrepancies in important metabolic pathways, such as sulfur metabolism and regulation of photosynthesis, although they did not belong to the COG categories enriched in the flexible fraction (Supplementary Results and Table S4).

After ten years in culture, genetic drift could be responsible to some extent for the observed genomic diversity among strains (e.g., [15]). However, because the number of proteins in the pangenome is much higher than that of individual A. gracile strains, which is constrained by their genome size, this diversity probably reflects true genomic diversity in the natural population, as has already been observed at the level of individual marker genes, such as sxtA [16]. Our results lead to the conclusion that the intra-population micro-diversity could be a feature largely shared within the cyanobacterial phylum. On the one hand, the presence of accessory genes related to the regulation of photosynthesis and sulfur metabolism implies that many strains with diverse adaptive capabilities can thrive until bloom development and co-occur at low spatial scale. On the other hand, the inter-strain variability in the production of essential secondary metabolites for which biosynthetic genes are especially enriched in accessory gene fraction, suggests that some of these molecules could be shared to the surrounding cyanobacteria by the producing strains, allowing the population to benefit from their adaptive functions while limiting the resource requirements for individuals. Thus, evaluation of the extent of biodiversity within a blooming cyanobacteria population could be a cornerstone for the establishment of a holistic theoretical framework of the dynamics of toxinogenic cyanobacteria blooms.

PMC strains are available upon request (https://mcam.mnhn.fr/en/cyanobacteria-and-live-microalgae-1281). Raw reads were deposited into the GENBANK Sequence Read Archive (SRA) database under the BioProject PRJNA693796 as well as A. gracile genome assemblies (JAQSYX000000000, JAQSYY000000000, JAQSYZ000000000, JAQSZA000000000).

Yepremian C, Gugger MF, Briand E, Catherine A, Berger C, Quiblier C, et al. Microcystin ecotypes in a perennial Planktothrix agardhii bloom. Water Res. 2007;41:4446–56.

Article CAS PubMed Google Scholar

Kardinaal WEA, Tonk L, Janse I, Hol S, Slot P, Huisman J, et al. Competition for light between toxic and nontoxic strains of the harmful cyanobacterium Microcystis. Appl Environ Microbiol. 2007;73:2939–46.

Article PubMed PubMed Central Google Scholar

Sabart M, Misson B, Jobard M, Bronner G, Donnadieu-Bernard F, Duffaud E, et al. Genetic diversity along the life cycle of the cyanobacterium Microcystis: highlight on the complexity of benthic and planktonic interactions. Environ Microbiol. 2015;17:901–11.

Article PubMed Google Scholar

Martiny AC, Coleman ML, Chisholm SW. Phosphate acquisition genes in Prochlorococcus ecotypes: evidence for genome-wide adaptation. Proc Natl Acad Sci USA. 2006;103:12552–7.

Article CAS PubMed PubMed Central Google Scholar

Martínez A, Osburne MS, Sharma AK, DeLong EF, Chisholm SW. Phosphite utilization by the marine picocyanobacterium Prochlorococcus MIT9301. Environ Microbiol. 2012;14:1363–77.

Article PubMed Google Scholar

Jackrel SL, White JD, Evans JT, Buffin K, Hayden K, Sarnelle O, et al. Genome evolution and host-microbiome shifts correspond with intraspecific niche divergence within harmful algal bloom-forming Microcystis aeruginosa. Mol Ecol. 2019;28:3994–4011.

Article CAS PubMed Google Scholar

Humbert J-F, Barbe V, Latifi A, Gugger M, Calteau A, Coursin T, et al. A tribute to disorder in the genome of the bloom-forming freshwater cyanobacterium Microcystis aeruginosa. PloS One. 2013;8:e70747.

Article CAS PubMed PubMed Central Google Scholar

Pérez-Carrascal OM, Terrat Y, Giani A, Fortin N, Greer CW, Tromas N, et al. Coherence of Microcystis species revealed through population genomics. ISME J. 2019;13:2887–2900.

Article PubMed PubMed Central Google Scholar

Coleman ML, Sullivan MB, Martiny AC, Steglich C, Barry K, Delong EF, et al. Genomic islands and the ecology and evolution of Prochlorococcus. Science. 2006;311:1768–70.

Article CAS PubMed Google Scholar

Biller SJ, Berube PM, Lindell D, Chisholm SW. Prochlorococcus: the structure and function of collective diversity. Nat Rev Microbiol. 2015;13:13–27.

Article CAS PubMed Google Scholar

Kashtan N, Roggensack SE, Rodrigue S, Thompson JW, Biller SJ, Coe A, et al. Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus. Science. 2014;344:416–20.

Article CAS PubMed Google Scholar

Takaichi S. Carotenoids in Algae: distributions, biosyntheses and functions. Mar Drugs. 2011;9:1101–18.

Article CAS PubMed PubMed Central Google Scholar

Singh SP, Klisch M, Sinha RP, Häder D-P. Effects of abiotic stressors on synthesis of the mycosporine-like amino acid shinorine in the cyanobacterium anabaena variabilis PCC 7937. Photochem Photobiol. 2008;84:1500–5.

Article CAS PubMed Google Scholar

Soto-Liebe K, Méndez MA, Fuenzalida L, Krock B, Cembella A, Vásquez M. PSP toxin release from the cyanobacterium Raphidiopsis brookii D9 (Nostocales) can be induced by sodium and potassium ions. Toxicon Off J Int Soc Toxinology. 2012;60:1324–34.

Article CAS Google Scholar

Willis A, Bent SJ, Jameson ID. Morphological changes and genome evolution in Raphidiopsis raciborskii CS-506 after 23 years in living culture. Appl Phycol. 2022;3:189–98.

Article Google Scholar

Ledreux A, Thomazeau S, Catherine A, Duval C, Yéprémian C, Marie A, et al. Evidence for saxitoxins production by the cyanobacterium Aphanizomenon gracile in a French recreational water body. Harmful Algae. 2010;10:88–97.

Article CAS Google Scholar

Download references

SH would like to thank Claude Yéprémian, for sharing his extensive knowledge about cyanobacteria. The authors thank the "Cyanobacteria and Eukaryotic Microalgae Collections" and the "Plateau technique de Spectrométrie de Masse Bio-organique" (UMR 7245 MCAM) of the MNHN, Paris, France.

This work was supported by the ATM grant "ECOL-Cyanomique" from the AVIV department of the MNHN, Paris. SKT was funded by an ATER grant from the MNHN.

Muséum National d’Histoire Naturelle, CNRS, UMR7245 Mécanismes de Communication et Adaptation des Micro-organismes, 12 rue Buffon, 75005, Paris, France

Sébastien Halary, Sébastien Duperron, Sandra Kim Tiam, Charlotte Duval, Cécile Bernard & Benjamin Marie

UMR5557 Laboratoire d’Ecologie Microbienne, Université de Lyon, 43 bd du 11 novembre 1918, Villeurbanne, F-69622, Lyon, France

Sandra Kim Tiam

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

BM and SH conceived this study. CD conducted biomass culture and DNA extraction. SH performed bioinformatics. SKT and BM performed metabolomics. SH, BM, SD and CB analyzed the results. SH wrote the paper with the contributions of all co-authors.

Correspondence to Sébastien Halary.

The authors declare no competing interests.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Halary, S., Duperron, S., Kim Tiam, S. et al. Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale. ISME COMMUN. 3, 57 (2023). https://doi.org/10.1038/s43705-023-00263-3

Download citation

Received: 23 November 2022

Revised: 09 May 2023

Accepted: 24 May 2023

Published: 07 June 2023

DOI: https://doi.org/10.1038/s43705-023-00263-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative