
==== Front
iScience
iScience
iScience
2589-0042
Elsevier

S2589-0042(24)02070-4
10.1016/j.isci.2024.110845
110845
Article
The genus Limnospira contains only two species both unable to produce microcystins: L. maxima and L. platensis comb. nov.
Pinchart Pierre-Etienne 12
Marter Pia 34
Brinkmann Henner 3
Quilichini Yann 1
Mysara Mohamed 56
Petersen Jörn joern.petersen@dsmz.de
3∗
Pasqualini Vanina 1
Mastroleo Felice felice.mastroleo@sckcen.be
57∗∗
1 UMR SPE, UAR Stella Mare, Université de Corse, 20250 Corte, France
2 Fédération des Spiruliniers de France (FSF), 34800 Clermont-l'Hérault, France
3 Leibniz Institute DSMZ – German Collection of Microorganisms and Cell Cultures, 38124 Braunschweig, Germany
4 Julius Kühn-Institut (JKI), 06484 Quedlinburg, Germany
5 Microbiology Unit, Nuclear Medical Applications, Belgian Nuclear Research Centre, SCK CEN, 2400 Mol, Belgium
6 Bioinformatics Group, Center for Informatics Science, Nile University, Giza 12677, Egypt
∗ Corresponding author joern.petersen@dsmz.de
∗∗ Corresponding author felice.mastroleo@sckcen.be
7 Lead contact

30 8 2024
20 9 2024
30 8 2024
27 9 1108456 5 2024
28 6 2024
27 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by-nc/4.0/ This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).
Summary

Spirulina is the commercial name for edible cyanobacteria of the genus Limnospira. The taxonomy of this genus is confusing with four species distributed in two lineages. Furthermore, the species Limnospira fusiformis has been cited as toxic by potentially producing microcystins. Taxonomic ambiguity combined with suspected health concerns constitute a major issue for spirulina producers. In a collection of six cultivars and one ecotype, we identified strains of the two lineages through metagenetic and morphological analyses. We demonstrated that the genus Limnospira only comprises two distinct species according to genomic comparisons of three genomes obtained in this study and 19 reference genomes. We showed that the V3-V4 region of the 16S rRNA gene is sufficient to identify the genus Limnospira and to distinguish the two species. Toxinogenesis investigations on eleven genomes from each Limnospira species revealed no genes involved in cyanotoxin synthesis, reflecting the inability of this genus to produce microcystins.

Graphical abstract

Highlights

• The genus Limnospira comprises two species: L. maxima and L. platensis comb. nov.

• The V3-V4 region of the 16S rRNA gene distinguishes these two Limnospira species

• The genus Limnospira is unable to produce microcystins

Molecular biology; Toxicology; Microbiology

Subject areas

Molecular biology
Toxicology
Microbiology
Published: August 30, 2024
==== Body
pmcIntroduction

Spirulina is the traditional and commercial name for edible filamentous cyanobacteria, which was first assigned to the genera Spirulina and Arthrospira,1 and more recently to the genus Limnospira.2 Spirulina is naturally present in alkaline fresh water, rich in carbonates, bicarbonates, nitrates, phosphates, and iron. It grows mainly in the intertropical regions with a wide geographical distribution (Africa, Asia, and Americas) and also in Southern Europe.3,4,5,6 Spirulina has long been used as human food, notably by the Kanembou (Lake Chad) and the Aztecs (Lake Texcoco).7 Its exceptional nutritional value has led the United Nations Food and Agriculture Organization (FAO) to consider spirulina a food of major interest against malnutrition for the food security and the response to emergency food situations.8 Spirulina also has strong potential for health, especially in the fields of cardiovascular diseases, viral infections, cancer prevention and therapy, immune response, as well as diabetes and cholesterol control.9,10

Many non-edible cyanobacteria are known to produce various cyanotoxins, and “blue-green algae” dietary supplements have long been suspected of containing toxic cyanobacteria and cyanotoxins.11 In 2018, after a routine check to identify a new strain of spirulina, three French producers growing the Arthrospira fusiformis Camargue ecotype were forced to empty their ponds following a warning about this strain’s suspected toxicity. Indeed, Limnospira (Arthrospira) fusiformis is cited in inventories of toxic cyanobacteria,12,13 which refer to the publications of Ballot et al. reporting production of microcystins by Kenyan strains of spirulina classified as A. fusiformis.14,15 Since then, numerous genomic studies on several spirulina strains have consistently shown the absence of toxicity of Limnospira sp.6,16,17,18 However, none of these strains was identified as L. fusiformis in these publications. To reassure the producers and food safety authorities, one or more strains of L. fusiformis should be identified to confirm that this species cannot produce cyanotoxins.19

Another question for spirulina producers concerns the identification and naming of cultivated spirulina strains. Arthrospira platensis was for a long time considered as the only species of spirulina.20,21,22 The United States and European authorities consider two species, Arthrospira maxima and A. platensis.23,24 The species A. fusiformis is sometimes cited as a synonym of A. platensis or as a synonym of A. maxima.24,25 When assessing commercially grown taxa in the light of the new genus Limnospira, the authors included three species (L. fusiformis, Limnospira indica, and Limnospira maxima) without ruling on the species A. platensis.2 This assignment is compatible with the current “List of Prokaryotic names with Standing in Nomenclature” (LPSN).26 Several studies have shown the existence of two or three lineages,3,6,27,28,29 but two recent studies proposed that the genus Limnospira is monospecific.30,31 However, the “paracas” strain is perceived by producers as a species different from other cultivated Limnospira strains due to its morphotype and its higher productivity. This strain also appears to have a specific sensitivity to a cyanophage that causes “sudden and massive deaths.”32 The phylogenetic markers used in different studies—internally transcribed spacer (ITS), gvp gene cluster, 16S rRNA gene, phycocyanin intergenic spacer (cpcBA-IGS), and 16S rRNA-ITS sequences—did not allow to determine with certainty the number of species belonging to the genus Limnospira. DNA-DNA hybridization (DDH) is considered “the gold standard” for species delineation in prokaryotic taxonomy.33,34 To our knowledge, no taxonomic study based on this method, either in vitro or in silico, has been conducted on Limnospira and Arthrospira genomes. Accordingly, the current study was aimed to determine the number of different bacterial species designated as “spirulina” and to clarify their nomenclature.

In this study, we present the results of 16S rRNA gene amplicon sequencing and metagenomic analysis, complemented with microscopic observations, carried out on seven spirulina culture samples of different origin, including the French “camarguaise” strain. The results of comparisons between 22 genomes, or metagenome-assembled genomes, carried out with the “average orthologous nucleotide identity” (OrthoANI) and “digital DNA-DNA hybridization” (dDDH) methods, allowed us to propose a taxonomic revision of the genus Limnospira. We also present an advanced genomic analysis of the identified L. fusiformis and other strains, targeting their potential ability to produce toxins. It is discussed if it is still advisable to mention L. (A.) fusiformis as a microcystin-producing species.

In the following, we use the genus names Limnospira and Arthrospira to refer to the spirulina strains, depending on the context. The strains mentioned in this study are named according to the names used in the GenBank database and other cited studies. Taxonomic changes will be proposed in the discussion.

Results

16S rRNA gene amplicon sequencing

The 16S rRNA gene sequences (400–405 bp) resulting from amplicon sequencing of the V3 and V4 hypervariable regions allowed the identification of two operational taxonomic units (OTUs) corresponding to Limnospira (Arthrospira) and five OTUs corresponding to other cyanobacteria (Table 1).Table 1 Cyanobacterial OTUs

OTU	Taxa1	pid2	Samples	
CAM_GP	ETH_PS	ETH_VR	IND_PP	PAR_GP	PAR_RB	SEN_CD	
OTU1_TG (clade I)	Limnospira fusiformis, L. maxima	100	19582 (95.28%)	9199 (100%)	25572 (88.76%)	18618 (99.63%)		280 (11.02%)	14954 (84.73%)	
OTU1_AA (clade II)	L. fusiformis, L. maxima, Arthrospira platensis	99.26			27 (0.09%)	5 (0.03%)	31784 (100%)	2 (0.08%)	19 (0.11%)	
OTU6_CA	Sodalinema sp.	100	2 (0.01%)		3210 (11.14%)	65 (0.35%)		2079 (81.82%)	2677 (15.17%)	
OTU26	Cyanobium sp.	97.50	667 (3.25%)							
OTU48	Cyanocohniella sp.	96.54	295 (1.44%)							
OTU65	Haloleptolyngbya sp.	99.26						180 (0.07%)		
OTU310	Nodosilinea sp.	97.77	5 (0.02%)							
		Totals	20551 (100%)	9199 (100%)	28809 (100%)	18688 (100%)	31784 (100%)	2541 (100%)	17650 (100%)	
Reads and proportion of each sample.

1 “Best hit” taxa from Blastn queries with 100% recovery query and higher percent identity.

2 Percentage identity between OTU and “best Hit” taxa sequences.

Concerning OTU1_AA, the alignments obtained with L. fusiformis and L. maxima on the one hand, and with A. platensis on the other hand have the same percentage identity (402/405, 99.26%), but the three nucleotide differences do not concern the same alignment positions.

Examination of the alignment of Limnospira OTUs sequences with other sequences from spirulina strains—notably the reference strains of the type species L. maxima UTEX 2720, L. fusiformis SAG 85.79, L. indica PCC 8005, and A. platensis PCC 7345—revealed three diagnostic nucleotide sites—namely positions 211, 237, and 241—that allowed us to differentiate between the two clades (Figure 1). OTU1_TG was found in L. maxima, L. fusiformis, and L. indica strains, while OTU1_AA is characteristic for A. platensis.Figure 1 Alignment position 200 to 250 of the V3-V4 region of the 16S rRNA gene

Differences between OTU1_TG (clade I) and OTU1_AA (clade II) are located in positions 211, 237, and 241. Complete sequences in Table S1.

Regardless of its comparably short sequence length, the V3-V4 region of the 16S rRNA gene is sufficient to distinguish between the two clades of Limnospira sp., which is in accordance with the phylogenomic trees and comparative genome analyses conducted in this study (see the following).

Morphology and ultrastructure

The seven spirulina culture samples exhibit a great diversity of Limnospira morphotypes, including spiral or helical (Figures 2A, 2B, 3A, and 3B), tightly coiled (Figures 2C and 3C), wavy or weakly spiral (Figures 2D and 3G), and straight trichomes (Figures 3D–3F).Figure 2 Scanning electronic microscopy (SEM) pictures of Limnospira strains

CAM_GP (A), ETH_PS (B), IND_PP (C), and PAR_GP (D). Scale bar: 50 μm.

Figure 3 Light microscopy (LM) views (f-m) of spirulina strains

CAM_GP (A), ETH_PS (B), IND_PP (C and H), ETH_VR (D), SEN_CD (E), PAR_RB (F), and PAR_GP (G). Scale bar: 50 μm.

The CAM_GP and ETH_PS samples are composed of Limnospira spiral trichomes (Figures 2A, 2B, 3A, and 3B) that match the “fusiformis” species description.35

The IND_PP sample is composed of Limnospira tightly coiled trichomes corresponding to the “Lonar” strain description (Figure 3C).21 Under laboratory conditions, the original monomorphotype (Figure 3C) evolved toward several morphotypes with helical trichomes (Figure 3H), which correspond to a former report of “Arthrospira maxima” strains.5

The ETH_VR and SEN_CD samples are mainly composed of Limnospira straight trichomes and are heavily contaminated by a filamentous cyanobacterium, a putative representative of the genus Sodalinema (Figures 3D and 3E). In these samples, some specimens are spiraled and match the “fusiformis” species description.35

The PAR_RB sample, a mixture of spirulina strains from both clades, was according to the 16S rRNA gene amplicon sequencing analysis also contaminated with a cyanobacterium of the genus Sodalinema (Table 1). It is mainly composed of straight Limnospira trichomes and some weakly spiraled specimens like the “paracas” strain (Figure 3F).21,31,36

The PAR_GP sample is composed of weakly spiraled Limnospira trichomes corresponding to the “paracas” strain trichome shape (Figure 3G).

Clade I samples (CAM_GP, ETH_PS, IND_PP, ETH_VR, SEN_CD, and PAR_RB; Table 1) present irregular helical shapes, with a pitch generally smaller than the trichome diameter, and with gradually attenuated trichome and cells diameters at the ends. We have observed a development from these screw-like coiled trichomes toward straight shapes in all our clade I samples in laboratory cultures. The clade II sample, the “paracas” strain (PAR-GP), presents a regular helical shape with a pitch approximately double the diameter. This shape appears wavy in light microscopy (LM) observations (Figure 3G). The morphology of this strain remains stable, without the appearance of straight shapes, in laboratory culture for more than three years.

The presence of a calyptra, a characteristic of the genus Limnospira,2 is observed on most trichome ends in both lineages. However, the calyptra and the attenuation of the diameter of the trichome appear to be more pronounced for the “fusiformis”-like strains (Figures 4A–4D).Figure 4 SEM views (scale bar: 5 μm) of Limnospira

CAM_GP (A), ETH_PS (B), PAR_GP (C and D) showing ends of trichomes and TEM views (scale bar: 1μ m) of Limnospira - PAR_GP (E and F) showing the ultrastructure including thylakoid membranes (Th), ribosomes (Ri), and carboxysomes (Ca). Blue triangles indicate the calyptra at the end of the trichomes.

Transmission electronic microscopy (TEM) observations did not reveal any remarkable differences between the two clades of Limnospira (Figure S1). In particular, we observed a thickened cell wall or calyptra at the end of trichomes and thylakoids with an irregular arrangement in PAR_GP samples (clade II, Figures 4E and 4F), as described by Nowicka-Krawczyk et al. for L. fusiformis strains belonging to clade I.2

Metagenome assembly and binning of spirulina

Metagenome sequencing, assembly and binning of two spirulina cultures from clade I (ETH_PS and CAM_GP) and one culture from clade II (PAR_GP) resulted in a single cyanobacterial metagenome-assembled genome (MAG) of the genus Limnospira per sample (Table 2). All three MAGs have a good quality with a calculated completeness of at least 96% and a contamination rate of less than 1% irrespective of the comparably large number of contigs. Mapping of the Illumina reads on the MAGs resulted in 1600-fold to 2900-fold genome coverage, which is nearly two orders of magnitude greater than required for reliable binning. For future studies, the sequencing depth, which was comparable to a former metagenome study of non-axenic limnic cyanobacteria,37 could be significantly reduced, thus paving the way for cost-effective metagenomic monitoring of spirulina cultivation. The analysis of the Limnospira MAGs showed a comparable G + C content (44.4%–44.7%), coding density (0.83–0.84) and a draft genome size ranging between 4.9 Mbp and 5.4 Mbp (Table 2). Accordingly, metagenome sequencing revealed similar genomic characteristics for Limnospira strains of both clades.Table 2 Limnospira MAGs and their characteristics

Sample-MAG	Clade	Taxon1	Comple-teness (%)	Contamination (%)	Mean contig length	N50 (contigs)	Coding density	# contigs	Genome size	G + C (%)	Cove-rage	
ETH_PS-MAG01	I	Limnospira fusiformis	97.09	00.44	6052	16808	0.83	881	5,332,277	44.67	1613	
CAM_GP-MAG01	I	Limnospira fusiformis	96.00	00.66	13773	20213	0.84	356	4,903,374	44.51	2912	
PAR_GP-MAG01	II	« Limnospira platensis »	98.14	00.00	7946	33073	0.83	678	5,387,508	44.43	1958	
1 According GTDB-Tk classification.

Phylogenomic analyses

Our three newly established MAGs were analyzed with 19 reference genomes of the genus Limnospira and ten genomes from closely related cyanobacteria all belonging to clade A of the cyanobacterial species tree.37,38 The phylogenomic protein tree based on 18,746 variable amino acid positions showed a common branching of all 22 analyzed Limnospira genomes in a distinct subtree that is supported by 100 bootstrap proportion (BP; Figure S2). The Limnospira subtree is divided in two distinct clades both supported by 100% BP. Clade I comprise two newly established metagenomes (CAM_GP and ETH_PS), several genomes of strains that are designated L. maxima, L. indica, and L. fusiformis, and it contains also putative “Arthrospira” genomes that were apparently misclassified (strain C1) or unclassified (strains O9.13F, TJSD091, and TJSD092). Clade II comprises genomes of different L. platensis strains including the new MAG from the PAR_GP culture. However, the internal branching pattern especially of clade I remained unresolved and the position of the comparably fast evolving Arthrospira sp. O9.13F genome might be influenced by long-branch attraction.39

With the aim of a better phylogenetic resolution, we reanalyzed the relationships of the 22 Limnospira strains with a maximum likelihood (ML) tree based on 32,918 variable nucleotide positions (Figure 5). The subanalysis showed a distinct branching pattern of the eleven A. platensis strains in clade II, comprising a common branch (98% BP) with three subtrees ([1] PCC 7345 – NIES-39, [2] Paraca, PAR_GP, and [3] YZ – FACHB-971) and A. platensis FACHB-439 in a basal position. In contrast, the close relationship of the strains in clade 1 with a largely unresolved topology was confirmed by the nucleotide tree. A notable finding was the common branching of Arthrospira sp. O9.13F with L. indica PCC 8005 and our L. maxima ETH_PS MAG, which is supported by 100% BP regardless of its accelerated evolutionary rate. Our phylogenomic analyses, as well as a previous phylogenetic analyze based on 16S-ITS (Roussel et al., 2023), showed that the supposed Arthrospira strains C1, O9.13F, TJSD091, and TJSD092 are located in clade I. However, our phylogenetic analyses showed a very close relationship of all strains of clade I, but provided no evidence for the presence of distinct lineages corresponding to the three previously proposed species (L. maxima, L. fusiformis, and L. indica).Figure 5 Phylogenomic maximum likelihood tree of three newly established Limnospira Metagenome-Assembled Genomes (MAGs) from spirulina cultures (CAM_GP, ETH_PS, and PAR_GP) and 19 cyanobacterial reference genomes from the genus Limnospira

The RaxML tree was constructed from 32,918 variable nucleotide positions based on an amino acid alignment of 92 housekeeping genes under the GTR4Γ model using 100 bootstrap replicates. Bootstrap values >50% are shown at the nodes. MAGs from the current study are highlighted in blue and bold. Genome assembly accession numbers are listed in key resources table. The bar represents relative number of changes per amino acid position.

Average Nucleotide Identity

The 231 pairwise orthologous average nucleotide identity (OrthoANI) comparisons among 22 Limnospira and Arthrospira strains are presented in Figure 6.Figure 6 Pairwise OrthoANI among Limnospira and Arthrospira strains

The Limnospira OrthoANI analysis is consistent with the results of the phylogenomic analyses (Figures 5 and S2) and allows to distinguish between the two clades (Figure 6). Clade I include all strains designated as L. maxima, L. fusiformis, and L. indica. Clade II includes only strains of the A. platensis lineage. The intra-clade OrthoANI values are greater than 98.7% with an average of 99.3% for clade I and 99.4% for clade II. The OrthoANI values between strains of different clades are ranging between 93.5% and 94.7% with an average of 93.9%. All values are below the ANI threshold of 95% (or 96%), which is generally accepted as the intra-species boundary for prokaryotes.40 Accordingly, all taxa of clade I belong to the same bacterial species, while the taxa of clade II represent a closely related second species of the genus Limnospira.

Digital DNA-DNA hybridization

The 231 pairwise genome comparisons using the dDDH among 22 Limnospira and Arthrospira strains are presented in Figure 7.Figure 7 Pairwise dDDH among Limnospira and Arthrospira strains

In agreement with the phylogenomic trees (Figures 5 and S2) and the OrthoANI analysis (Figure 6), the dDDH analysis allowed to distinguish between clade I and II (Figure 7). The intra-clade dDDH values are ranging between 85.4% and 99.4%. The inter-clade dDDH values are much lower at 49.6%–56.6%. They are clearly below the threshold of 70% for the delineation of two bacterial species,34,41 which confirms the conclusion of our OrthoANI analysis that clade I and clade II represent two different species of the genus Limnospira. In addition, the high dDDH values for all strains except O9.13F from clade I (≥89.9%) reflect their very close relationship. Arthrospira O9.13F has the most divergent genome of clade I (dDDH: 85.4%–89.3%), but nevertheless definitely belongs to the same species. This conclusion is independently supported by the placement of strain O9.13F in the phylogenomic tree (Figure 5). The contradictory position within a 16S-ITS tree in a recent study reflects the rather limited phylogenetic resolution of the 16S rRNA gene.31

Toxinogenesis

Blastn searches with the microcystin gene clusters of Microcystis aeruginosa K-139, Planktothrix agardhii NIVA-CY 126/8, and Anabaena sp. 90 (see method details section) revealed no hits in the three Limnospira MAGs established in the current study and 19 reference genomes of Limnospira and Arthrospira (see key resources table), which documents the absence of very closely related sequences. Analogous blastp and tblastn analyses with Myc proteins for microcystin biosynthesis (see method details section) did not result in the detection of any enzymes related to non-ribosomal peptide synthases or polyketide synthases. Accordingly, our homology searches revealed no indication for the presence of genes related to the biosynthesis of this cyanotoxin. In order to confirm this conclusion, our three newly established MAGs and 19 reference genomes of the genus Limnospira (Arthrospira) were analyzed with the antiSMASH 7.1.0 software. However, neither non-ribosomal peptide synthase (NRPS) nor polyketide synthase (PKS) clusters could be detected in any of the analyzed genomes. This result allowed us to conclude that none of the analyzed 22 strains of the genus Limnospira (Arthrospira) has the genetic capacity to produce microcystins or other cyanotoxins whose biosynthesis gene clusters are known.

Discussion

Taxonomy

There is a consensus on the existence of two distinct lineages in the genus Limnospira—using chronologically ITS region, gvp gene cluster, 16S rRNA gene, cpcBA-IGS, ANI, isDDH, and 16S-ITS3,6,27,29,30,31—which should be distinguished according to the criteria of prokaryotic taxonomy to solve the dilemma of the four species described under the “Botanical Code” (ICN, L. maxima, L. fusiformis, L. indica, “L. [A.] platensis”). However, a recent publication, based on the analysis of the 16S rRNA gene, 16S-ITS sequences, ANI and morphological characteristics, concluded that the genus Limnospira is monospecific with two distinct clades.31 The inter-clade ANI value reported in the aforementioned publication (93.4%) should lead to the conclusion that the genus Limnospira contains two species when using the ANI threshold of 95 or 96% for the differentiation between species.40,42 Values generated by the original ANI and OrthoANI are comparable, and the cut-off for species delineation is 95%∼96% for both methods.43 The OrthoANI values for Limnospira genomes compared in this study clearly support the presence of two species based on inter-clade OrthoANI values ranging between 93.5% and 94.6% and an intra-clade OrthoANI above 98.7% for clade I and above 99.0% for clade II. These results are comparable to those obtained for other genera of cyanobacteria. Species delineation in the cyanobacterial genera Raphidiopsis and Microcystis could be assessed with species boundaries occurring at about 96% ANI at the genomic level.44 ANI analysis on three Microcoleus species revealed inter-species ANI values between 86% and 93% and intra-species ANI values of more than 98%.45

DDH has been used by taxonomists since the 1960s as the most important criterion for the delineation of prokaryotic species.46 However, the determination of wet-lab DDH values is laborious and depends on DNA isolated from purified bacterial isolates, which formerly excluded non-axenic cyanobacteria from this crucial tool in bacterial taxonomy. With the aid of established algorithms, it is possible to calculate in silico digital DDH (dDDH) values that are very close to the wet-lab DDH values and have the same delineation threshold between species.34,47 A value below 70% DDH (or dDDH) is the generally accepted threshold for delineating species41,48 and a value of 79–80% has been proposed for subspecies delineation.49 Furthermore, progress in metagenome sequencing and binning paved the way to analyze the MAGs of non-axenic cyanobacteria with dDDH. The dDDH values, established in this study, are ranging between 85.4% and 95.2% for intra-clade comparisons and between 49.6% and 56.6% for inter-clade comparisons. These values clearly documented that the two clades of Limnospira correspond to two different species. Furthermore, the intra-species comparison within the same clade showed that none of the two species comprises further subspecies.

The separation of spirulina into two clades is confirmed by the different methods used in this study, which give consistent results on the affiliation of the different strains to one or the other clade. In particular, the sequences of the hypervariable regions V3 and V4 of the 16S rRNA gene allow an unequivocally identification of the genus Limnospira and are sufficient to distinguish its two species in agreement with genome-based methods (Figure 1).

The phylogenomic tree (Figure 5) documented that clade I includes the type species of the genus Limnospira Nowicka-Krawczyk et al. 2019, which is designated L. fusiformis (Woronichin 1934) Nowicka-Krawczyk et al. 2019 (reference strain: SAG 85.79) and L. indica (Desikachary and Jeeji Bai 1992) Nowicka-Krawczyk et al. 2019 (reference strain: PCC 8005).2 In addition, three diagnostic nucleotide positions in the V3-V4 region of the 16S rRNA gene showed that L. maxima (Setchell and Gardner 1917) Nowicka-Krawczyk et al. 2019 (reference strain: UTEX LB 2720) also belongs to clade 1 (Figure 1), and a comparison of the complete 16S rRNA genes revealed an identity of 100% with L. fusiformis SAG 85.79 and L. indica PCC 8005.31 Our OrthoANI comparison (Figure 6) and dDDH analysis (Figure 7) provided clear evidence that all strains of clade I belong to the same species. Accordingly, the current study showed that the species names Limnospira fusiformis (Woronichin 1934) Nowicka-Krawczyk et al. 2019 and Limnospira indica (Desikachary and Jeeji Bai 1992) Nowicka-Krawczyk et al. 2019 are heterotypic synonyms of the type species Limnospira maxima (Setchell and Gardner 1917) Nowicka-Krawczyk et al. 2019. Taxa of clade II includes the strain PCC 7345 of the species Arthrospira platensis Gomont 1892, which was recently proposed to be considered as the reference strain,31 as well as other strains assigned to the same species (Figure 5). This species was not transferred to the genus Limnospira given the benthic character in the description of A. platensis by Gomont in 1892.2 A recent study showed that A. platensis can present benthic and planktonic forms and suggests that this species should therefore be placed in the genus Limnospira.50 The 16S rRNA gene sequences published in the aforementioned study fall into our Limnospira clade II. Therefore, we propose the second species of the genus Limnospira to be named L. platensis. The type strain of this species is Limnospira platensis PCC 7345. A list of more than 130 spirulina strains and their correct species identification is provided in the supplementary material (Table S2).

Spirulina exhibits a diversity and morphological plasticity that does not allow a formal identification of one or the other species.3,5,21,31 Nevertheless, we observed that the shape of the trichomes of L. platensis strain paracas (PAR_GP) is distinguished by the regularity of the pitch and the diameter of the helix. This regularity of coiling is also visible in the images of L. platensis PCC 7345 available in the photo library of the Pasteur Institute collection (phototheque.pasteur.fr), as well as in the images of other strains of L. platensis.50 However, this distinguishing feature must be confirmed by further studies on strains of the two Limnospira species, whose identification will ideally be based on (meta-)genome sequencing.

Taxonomic considerations on the genus Limnospira (Nowicka-Krawczyk et al. 2019)

Biological material or DNA was neither available for the first described Arthrospira maxima species (Setchell and Gardner 1917) nor for the first described Arthrospira platensis species (Gomont 1892).

The original description of the genus Limnospira2 is given with the following amendments including data that describes the type specie (L. maxima) and his heterotypic synonyms (L. fusiformis and L. indica).

Class: Cyanophyceae.

Order: Oscillatoriales

Family: Microcoleaceae

Type species: Limnospira maxima (Setchell and Gardner 1917) Nowicka-Krawczyk et al. 2019.

Morphology: Trichomes solitary, isopolar, cylindrical, unbranched, more or less tapering toward the ends, not or slightly constricted at cross-walls, blue green or dark green. Trichomes are usually spirally coiled and can present different shapes—screw-like coiled, tightly coiled, weakly spiraled—with a tendency to loosening the coils until straight shape. Spiral diameter and pitch are usually irregular. Cells are always shorter (2–7 μm) than wide (6–12(16) μm). Terminal cells rounded or subcapitate with thickened outer cell wall or calyptra.

Ultrastructure: Thylakoids have an irregular arrangement. Aerotopes are facultatively present.

Ecology: Planktic in freshwater, brackish and alkaline-saline waters. Usually photoautotrophic and able of mixotrophy.51,52

Etymology: Límni- [Gr., λίμνη], lake, lagoon; -speíra [Gr., σπείρα], spiral; the generic name refers to the morphology of trichomes and typical habitat, i.e., lakes.

Limnospira maxima (Setchell and Gardner 1917) Nowicka-Krawczyk et al. 2019. Basonym: Arthrospira maxima (Setchell and Gardner 1917)

Heterotypic synonym: Limnospira fusiformis (Woronichin 1934) Nowicka-Krawczyk et al. 2019; Limnospira indica (Desikachary and Jeeji Bai 1992) Nowicka-Krawczyk et al. 2019.

Type strain: PCC 8005. Deposited in Pasteur Cyanobacteria Collection (France, as “Arthrospira sp.”), in DSMZ-German Collection of Microorganisms and Cell Cultures (Germany) and in SCK CEN (Belgium).

Limnospira platensis (Gomont 1892) comb. nov.

Basonyms: Arthrospira platensis (Gomont 1892), Spirulina platensis (Geitler 1925).

Type strain: PCC 7345. Deposited in Pasteur Cyanobacteria Collection (France, as “Arthrospira sp.”) and in The University of Texas Culture Collection of Algae (USA, = UTEX LB 1926 as “Spirulina platensis”).

The description of Limonospira platensis is identical to that of Limnospira genus except for the following modifications.

Morphology: Regular helical shape with constant diameter and pitch, the latter greater than the diameter. The shape is stable without tendency to loosening the coils.

Ecology: Planktic, and rarely benthic.

Limnospira toxinogenesis

Examination of the three Limnospira sp. MAGs from this study and 19 reference genomes showed that the genus Limnospira do not contain any genes that are involved in the biosynthesis of microcystin. The results of Ballot et al. who reported about microcystins in Kenyan strains of L. fusiformis have never been reproduced by other studies focusing on the detection of toxins in Limnospira cultures.14,15,53,54 The hypothesis that the “Arthrospira fusiformis” cultures examined by Ballot et al. were contaminated by toxin-producing cyanobacteria of other genera cannot be ruled out. The protocol for isolating A. fusiformis has not been described, the isolates were not deposited in a public culture collection, and Lakes Sonachi and Bogoria, from which the studied “A. fusiformis” isolates originated, contain the microcystin-producing picocyanobacteria Syncechococcus sp. and Synechocystis sp.14,15 Moreover, ten genes (mcyA to mcyJ), which are located in a 55 kb gene cluster in the genome of Microcystis aeruginosa PCC 7806SL and Anabaena sp. 90, are required for the biosynthesis of microcystin. None of these genes or corresponding proteins could be found in the 22 genomes studied, including four strains previously identified as L. (A.) fusiformis (SAG 85.79, KN01, CAM_GP, and ETH_PS). Therefore, the assumption that certain strains of Limnospira can produce microcystins does not appear to be a credible hypothesis. Consequently, we propose that Limnospira sp., in particular L. maxima (previously L. fusiformis), should no longer be included in the lists of toxic cyanobacteria.

Limitations of the study

The Kenyan strains of “Arthrospira fusiformis” studied by Ballot et al. in 2005 are no longer available. We would have liked to include them in our study to check their supposed toxicity.

The observation of a morphological difference between the two Limnospira species was carried out with a single strain of Limnospira platensis. Although this observation is confirmed by images available on the web and in a publication,50 it is necessary to complete the observations on a larger number of strains.

Resource availability

Lead contact

Further information and requests can be directed to Dr. Felice Mastroleo (felice.mastroleo@sckcen.be).

Materials availability statement

This study did not generate new unique reagents.

Data and code availability

• The genetic data used in this article are available in NCBI GenBank. Metagenome data are deposited in the following BioProjects: PRJNA1085564 (Limnospira maxima CAM_GP), PRJNA1085565 (Limnospira maxima ETH_PS), and PRJNA1085629 (Limnospira platensis PAR_GP). Fastq files of the 16S rRNA amplicon were submitted to NCBI’s Short Read Archive and are publicly available under accession number PRJNA1095467.

• This paper does not report original code.

• Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

Acknowledgments

We would like to thank the spirulina producers—Gilles Planchon, Philippe Schoeffert, Camille Dameron, Vincent Rioux, and Rémy Bosc—who provided us with the spirulina strains used in this study. We would also like to thank Markus Göker for helpful discussion on taxonomy. This study is financially supported by the French Government and the Corsican Regional Council (CPER Gerhyco project) in the framework of the work of P.-E.P. toward his PhD. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; and in the writing of the manuscript and the decision to publish the results.

Author contributions

Conceptualization, P.-E.P., J.P., V.P., and F.M.; methodology, P.-E.P., P.M., and H.B.; validation, V.P., J.P., and F.M.; formal analysis, P.-E.P., P.M., and H.B.; investigation, P.-E.P.; resources, V.P.; data curation, P.-E.P. and P.M.; writing – original draft preparation, P.-E.P.; writing – review and editing, F.M., J.P., and P.-E.P.; visualization, P.-E.P. and J.P.; supervision, V.P., J.P., and F.M.; project administration, P.-E.P., V.P., and F.M.; funding acquisition, V.P. and F.M. All authors have read and agreed to the published version of the manuscript.

Declaration of interests

The authors declare no competing interests.

STAR★Methods

Key resources table

REAGENT or RESOURCE	SOURCE	IDENTIFIER	
Biological samples	
	
Spirulina cultures	This paper	Table 1	
	
Deposited data	
	
Genome sequences of 19 Limnospira and Arthrospira reference strains	NCBI GenBank	Accession numbers in Table 2	
Meta Assembled Genomes (MAGs)	This paper	NCBI Bioproject: PRJNA1085564, PRJNA1085565, PRJNA1085629	
V3-V4 region of the 16S rRNA gene sequences	This paper	NCBI Bioproject: PRJNA1095467	
	
Software and algorithms	
	
OCToPUS pipeline v.1.0	Mysara et al.55	https://github.com/M-Mysara/OCToPUS	
Ribosomal Database Project (RDP v.18)	Wang et al.56	https://rdp.cme.msu.edu/	
MEGA 11	Tamura et al.57	https://www.megasoftware.net/	
MEGAHIT	Li et al.58	https://github.com/voutcn/megahit	
MaxBin 2.0	Wu et al.59	https://sourceforge.net/projects/maxbin2/	
MetaBAT	Kang et al.60	https://bitbucket.org/berkeleylab/metabat/src/master/	
Concoct	Alneberg et al.59	https://github.com/BinPro/CONCOCT	
DAS Tool	Sieber et al.61	https://github.com/cmks/DAS_Tool	
CheckM	Parks et al.62	https://ecogenomics.github.io/CheckM/	
DFAST	Tanizawa et al.63	https://github.com/nigyta/dfast_core/	
InterProScan	Jones et al.64	https://github.com/ebi-pf-team/interproscan	
SignalP	Almagro Armenteros et al.65	https://services.healthtech.dtu.dk/services/SignalP-5.0/	
Genome Taxonomy Database Toolkit GTDB-Tk v2.1.0	Parks et al.66	https://github.com/ecogenomics/gtdbtk.	
RAxML v8.2.10	Stamataki.67	https://github.com/stamatak/standard-RAxML.	
Orthologous Average Nucleotide Identity Tool v0.93.1.	Lee et al.43	https://www.ezbiocloud.net/tools/orthoani	
Genome-to-Genome Distance Calculator v3.0	Meier-Kolthoff et al.47	https://ggdc.dsmz.de/ggdc.php/	
BLAST+	Camacho et al.68	https://blast.ncbi.nlm.nih.gov/Blast.cgi	
antiSMASH 7.1.0	Blin et al.69	https://antismash.secondarymetabolites.org	

Experimental model and study participant details

Spirulina strains

Sixteen samples of different spirulina cultures, from producers of the “Fédération des Spiruliniers de France (FSF, www.spiruliniersdefrance.fr)”, constituted the initial collection of spirulina strains. A survey of producers made it possible to identify the initial suppliers of the strains used, and to reduce our sampling to seven strains originating from France, Ethiopia, India, Peru and Senegal. These spirulina strains and their uses in this study are listed below.Strain ID	Putative origin	Provider	16S rRNA (V3-V4)	Microscopy	Genome assembly	
CAM_GP	France (Camargue)	spirulinasolution.fr	•	•	•	
ETH_PS	Ethiopia	spirulinepourquoipas.com	•	•	•	
ETH_VR	Ethiopia	spirulinedugarlaban.com	•	•	–	
IND_PP	India (Lonar)	laspiraleverte.com	•	•	–	
PAR_GP	Peru (Paracas)	spirulinasolution.fr	•	•	•	
PAR_RB	Peru (Paracas)	algosud.com	•	•	–	
SEN_CD	Senegal	spirulinedeslandes.com	•	•	–	

Optical microscopy observations, unless otherwise stated, and DNA extractions for amplicon sequencing were carried out on the samples received from the producers; electron microscopy observations and metagenomics were carried out on samples sub-cultivated in the laboratory. Spirulina samples were cultivated in cell culture flask with vented cap (Greiner Bio-One, Belgium), in modified Zarrouk medium70 at room temperature and following day/night regime illumination. Flasks were continuously mixed using an orbital shaker at 120 rpm (Edmund Bülher, Germany).

Method details

Microscopy

Optical LM observations were carried out on living and cultivated samples using a Zeiss Axiolab microscope, an Axiocam 305 camera and Zeiss Zen 3.6 software.

For scanning electronic microscopy (SEM) and TEM, the specimens were fixed in cold (4°C) and saline (10 g NaCl per liter) solution of 2.5% glutaraldehyde in 0.1 M sodium cacodylate buffer at pH 7.2.

For SEM observations, one milliliter drop of this cell suspension were deposed on poly-L-lysine coated slides. After settling for 2 h, the slides were rinsed in Milli-Q water and dehydrated in a graded ethanol series: 30%, 50%, 70%, 90% and 100% (30 min for each immersion step). The last immersion in ethanol 100% was repeated twice. The samples were critical point dried in an Emitech K850 instrument, mounted on aluminum stubs with the use of double-sided carbon adhesive and coated with gold/palladium in a Quorum Technologies SC7640 sputter coater. The samples were examined under a Hitachi S-3400-N scanning electron microscope operated at an accelerating voltage of 10 kV.

For TEM observations, specimens’ pellets obtained after centrifugation, were post-fixed in 2% osmium tetroxide in 0.1 M cacodylate buffer (pH 7.2) for 1 h at room temperature. The postfixed specimens were processed in graded ethyl alcohol, propylene oxide, Spurr resin and polymerized for 24 h at 60°C. Ultrathin sections (60–90 nm) were cut on an ultramicrotome Power tome PC-RMC Boeckeler, placed on 300 and 200 mesh copper grids and stained with uranyl acetate followed by lead citrate. The samples were examined under a transmission electron microscope Hitachi H-7650 at an accelerating voltage of 80 kV.

DNA isolation

DNA from 5 mL spirulina culture samples was isolated using the protocol adapted for filamentous cyanobacteria of the genus Limnospira (Arthrospira) developed by Morin et al.71 Briefly, cells were collected by centrifugation (10 min, 3500 rpm) and resuspended into 0.5 mL of 0.15 M NaCl and 0.1 M EDTA. Following three freeze-thaw cycles, the cells were collected by centrifugation (10 min, 8000 rpm), resuspended in 0.5 mL TE Buffer (10 mM Tris, 1 mM EDTA, pH 8.0) and transferred into fresh tubes for enzymatic cell wall lysis using 100 μL of 50 mg/mL lysozyme, for 30 min at 37°C. Proteins were digested with 5 μL of 50 mg/mL proteinase K in 2% SDS final concentration, for 1 h at 37°C. Polysaccharides, proteins and cell wall debris were thereafter removed by selective precipitation: 150 μL of 5 M NaCl was added to the tubes, followed by 0.1 volume of a 10% cetyltrimethylammonium bromide (CTAB) stock solution and incubated at 65°C for 10 min. Nucleic acids purification was achieved by extraction in 1 volume of chloroform:isoamyl alcohol. The tubes were placed on ice for 30 min to allow precipitation of CTAB complexes, before being centrifuged (10 min, 8000 rpm). The supernatant was transferred into a fresh tube, gently mixed with 0.6 volume of isopropanol until DNA precipitated. The DNA pellets were recovered by centrifugation, (10 min, 4°C, 15000 rpm) and washed with 1 mL of ice-cold 70% ethanol to remove any residual salt. After a final centrifugation (5 min, 4°C, 15000 rpm), the supernatant was discarded, and the pellets were dried before being resuspended in 100 μL TE Buffer. The samples were treated with 1 μL of RNAse (10 mg/mL) for 1 h at 37°C and then kept at −80°C until further use.

16S rRNA gene amplicon sequencing

High-throughput paired-end amplicon sequencing of the V3-V4 hypervariable region of the 16S rRNA gene was conducted on BaseClear’s Illumina MiSeq (V3 chemistry) platform according to the manufacturer’s guidelines. Sequence data of the 16S rRNA gene amplicons were processed using OCToPUS pipeline v.1.0,55 and analyzed as previously described,72 generating the OTUs. The OTUs were classified using Ribosomal Database Project (RDP v.18) dataset56 and their purity were assessed using the oligotyping approach73 and classified using Blastn algorithm choosing the rRNA/ITS database. The sequence alignments and their examination were carried out with MEGA 11.57

Metagenome sequencing and binning

High-throughput short-read metagenome sequencing was conducted as previously described.37 The NEBNext Ultra II FS DNA Library Prep Kit (New England Biolabs) was used for the construction of the Illumina libraries. Sequencing was performed on the Illumina NovaSeq 6000 system using the v3 chemistry with 600 cycles. The fastq-mcf tool of ea-utils v1.04.803 was used for quality control and adapter clipping.

We used our previously established metagenome assembly and binning pipeline for the analysis of short-read Illumina sequences.37 In short, processed sequence data were assembled with MEGAHIT,58 metagenomic binning was conducted with three methods, namely MaxBin 2.0,59 MetaBAT,60 Concoct,74 and subsequently analyzed with DAS Tool61 to retrieve the most reliable MAGs. The quality of the MAGs was investigated with CheckM62 and coverage values were calculated per bbmap. Finally, genomes were annotated by DFAST,63 InterProScan64 and SignalP.65 The initial taxonomic classification of the cyanobacterial MAGs was conducted with the Genome Taxonomy Database Toolkit GTDB-Tk v2.1.0 on the GTDB reference data version r207.66

Phylogenomic analyses

Genome-based phylogenomic analyses of cyanobacteria were conducted as previously reported.75 Briefly, the phylogenomic trees were reconstructed from concatenated amino acid or nucleotide alignments of 92 housekeeping genes (Up-to-Date-Bacterial Core Gene [UBCG] dataset.76 The alignments were manually refined with the edit-option of the MUST package77 and the application of G-blocks;78 The ML trees were calculated with RAxML v8.2.10 and inferred under a GTR+4Γ model using 100 bootstrap replicates.67

Genome comparison

ANI was calculated using the OrthoANI algorithm with OAT (Orthologous ANI Tool) version 0.93.1. OrthoANI is highly correlated with ANI and show approximately 0.1% higher values than the original ANI.43

DDDH was calculated using Genome-to-Genome Distance Calculator (GGDC 3.0, formula 2) with the platform available under https://tygs.dsmz.de.34,47 The Limnospira and Arthrospira genomes used for the OrthoANI and the dDDH analysis are those published in GenBank as available on 09/30/2023, complemented with the later released assembly of the reference strain PCC 7345 of A. platensis (see key resources table for details).

Toxinogenesis analysis

The search for nucleotide and protein sequences for the synthesis of microcystins, was carried out based on the current knowledge of cyanotoxins and their metabolic synthesis pathways.79,80 In summary, microcystins are secondary metabolites, cyclic heptapeptides whose biosynthesis depends on ten genes (mcyA-J), grouped in a cluster of approximately 55kbp, which encode a set of multifunctional enzymes: non-ribosomal peptide synthetases (NRPSs), polyketide synthetases (PKSs), tailoring enzymes and a transporter enzyme. All sequences were retrieved from GenBank.

Genes involved in microcystin biosynthesis

Sequence description	Organism	Length (bp)	Accession	
mcyD, mcyE, mcyF, mcyG, mcyH, mcyI, mcyJ, dnaN genes	Microcystis aeruginosa K-139	36,468	AB032549	
complete gene cluster: mcyA to mcyJ	Planktothrix agardhii NIVA-CYA 126/8	55,575	AJ441056	
complete gene cluster: mcyA to mcyJ	Anabaena sp. 90	57,817	AJ536156	

Proteins involved in microcystin biosynthesis by Microcystis aeruginosa PCC 7806SL

Protein	Name	Length (aa)	Accession	
McyA	non-ribosomal peptide synthetase	2,795	ARI82521	
McyB	non-ribosomal peptide synthetase	2,126	ARI82522	
McyC	non-ribosomal peptide synthetase	1,291	ARI82523	
McyD	type I polyketide synthase	3,906	ARI82520	
McyE	hybrid non-ribosomal peptide synthetase/type I polyketide synthase	3,487	ARI82519	
McyF	aspartate/glutamate racemase family protein	251	ARI82518	
McyG	type I polyketide synthase	2,631	ARI82517	
McyH	ATP-binding cassette domain-containing protein	585	ARI82516	
McyI	phosphoglycerate dehydrogenase	337	ARI82515	
McyJ	methyltransferase domain-containing protein	311	ARI82514	

Nucleotide sequence was searched using the blastn algorithm and protein sequence was investigated with tblastn and blastp.68 We also used the software ‘antibiotics and secondary metabolite analysis shell - antiSMASH 7.1.0’69 to detect possible NRPS and PKS in the Limnospira genomes.

Quantification and statistical analysis

Phylogenomic analysis under ML (Figures 5 and S2) was conducted with RAxML (Randomized Axelerated ML) version 8.2.10.67 Number of bootstrap replicates (n = 100) are mentioned in figure legend. Bootstrap values are show at the phylogenomic tree nodes.

Additional resources

No additional resources.

Supplemental information

Document S1. Figures S1 and S2, Tables S1 and S2

Supplemental information can be found online at https://doi.org/10.1016/j.isci.2024.110845.
==== Refs
References

1 Castenholz R.W. Rippka R. Herdman M. Wilmotte A. Form-Arthrospira Bergeys Man. Syst. Archaea Bact. 1 2015 1 3 10.1002/9781118960608.gbm00433
2 Nowicka-Krawczyk P. Mühlsteinová R. Hauer T. Detailed characterization of the Arthrospira type species separating commercially grown taxa into the new genus Limnospira (Cyanobacteria) Sci. Rep. 9 2019 694 10.1038/s41598-018-36831-0 30679537
3 Baurain D. Renquin L. Grubisic S. Scheldeman P. Belay A. Wilmotte A. Remarkable conservation of Internally Transcribed Spacer sequences of Arthrospira (Spirulina) (Cyanophyceae, Cyanobacteria) strains from four continents and of recent and 30-year-old dried samples from Africa J. Phycol. 38 2002 384 393 10.1046/j.1529-8817.2002.01010.x
4 Fužinato S. Fodora A. Subakov-Simić G. Arthrospira fusiformis (Voronichin) Komarek et Lund (Cyanoprokaryota) A new species for Europe Algol. Stud. 134 2010 17 24 10.1127/1864-1318/2010/0134-0017
5 Sili C. Torzillo G. Vonshak A. Arthrospira (Spirulina) Ecology of cyanobacteria II 2012 Springer 677 705
6 Misztak A.E. Waleron M. Furmaniak M. Waleron M.M. Bazhenova O. Daroch M. Waleron K.F. Comparative genomics and physiological investigation of a new Arthrospira/Limnospira strain O9. 13F isolated from an alkaline, winter freezing, siberian lake Cells 10 2021 3411 10.3390/cells10123411 34943919
7 Ciferri O. Spirulina, the Edible Microorganism Microbiol. Rev. 47 1983 551 578 10.1128/mr.47.4.551-578.1983 6420655
8 Capelli B. Cysewski G.R. Potential health benefits of spirulina microalgae Nutrafoods 9 2010 19 26 10.1007/BF03223332
9 Wu Q. Liu L. Miron A. Klímová B. Wan D. Kuča K. The antioxidant, immunomodulatory, and anti-inflammatory activities of Spirulina: an overview Arch. Toxicol. 90 2016 1817 1840 10.1007/s00204-016-1744-5 27259333
10 Furmaniak M.A. Misztak A.E. Franczuk M.D. Wilmotte A. Waleron M. Waleron K.F. Edible cyanobacterial genus Arthrospira: actual state of the art in cultivation methods, genetics, and application in medicine Front. Microbiol. 8 2017 2541 10.3389/fmicb.2017.02541 29326676
11 Carmichael W.W. The Toxins of Cyanobacteria Sci. Am. 270 1994 78 86 10.1038/scientificamerican0194-78 8284661
12 Bernard C. Ballot A. Thomazeau S. Maloufi S. Furey A. Mankiewicz Boczek J. Pawlik Skowronska B. Capelli C. Salmaso N. Cyanobacteria Associated with the Production of Cyanotoxins Handbook on cyanobacterial monitoring and cyanotoxin analysis 2017 Wiley 501 525
13 Chorus I. Welker M. Toxic Cyanobacteria in Water: A Guide to Their Public Health Consequences, Monitoring and Management 2021 Taylor & Francis
14 Ballot A. Krienitz L. Kotut K. Wiegand C. Metcalf J.S. Codd G.A. Pflugmacher S. Cyanobacteria and cyanobacterial toxins in three alkaline Rift Valley lakes of Kenya—Lakes Bogoria, Nakuru and Elmenteita J. Plankton Res. 26 2004 925 935 10.1093/plankt/fbh084
15 Ballot A. Krienitz L. Kotut K. Wiegand C. Pflugmacher S. Cyanobacteria and cyanobacterial toxins in the alkaline crater lakes Sonachi and Simbi, Kenya Harmful Algae 4 2005 139 150 10.1016/j.hal.2004.01.001
16 Janssen P.J. Morin N. Mergeay M. Leroy B. Wattiez R. Vallaeys T. Waleron K. Waleron M. Wilmotte A. Quillardet P. Genome Sequence of the Edible Cyanobacterium Arthrospira sp. PCC 8005 J. Bacteriol. 192 2010 2465 2466 10.1128/JB.00116-10 20233937
17 Fujisawa T. Narikawa R. Okamoto S. Ehira S. Yoshimura H. Suzuki I. Masuda T. Mochimaru M. Takaichi S. Awai K. Genomic structure of an economically important cyanobacterium, Arthrospira (Spirulina) platensis NIES-39 DNA Res. 17 2010 85 103 10.1093/dnares/dsq004 20203057
18 Cheevadhanarak S. Paithoonrangsarid K. Prommeenate P. Kaewngam W. Musigkain A. Tragoonrung S. Tabata S. Kaneko T. Chaijaruwanich J. Sangsrakru D. Draft genome sequence of Arthrospira platensis C1 (PCC9438) Stand. Genomic Sci. 6 2012 43 53 10.4056/sigs.2525955 22675597
19 Pinchart P.-E. Leruste A. Pasqualini V. Mastroleo F. Microcystins and Cyanobacterial Contaminants in the French Small-Scale Productions of Spirulina (Limnospira sp.) Toxins 15 2023 354 10.3390/toxins15060354 37368655
20 Falquet J. Hurni J.P. Spiruline: Aspects Nutritionnels (Flamant Vert France) 1986
21 Fox R.D. Spirulina : Production & Potential (Edisud) 1996
22 Jourdan J.-P. Grow Your Own Spirulina 2001 Geneva Switz
23 Marles R.J. Barrett M.L. Barnes J. Chavez M.L. Gardiner P. Ko R. Mahady G.B. Low Dog T. Sarma N.D. Giancaspro G.I. United States Pharmacopeia Safety Evaluation of Spirulina Crit. Rev. Food Sci. Nutr. 51 2011 593 604 10.1080/10408391003721719 21793723
24 ANSES Avis de l’Agence nationale de sécurité sanitaire de l’alimentation, de l’environnement et du travail relatif aux « risques liés à la consommation de compléments alimentaires contenant de la spiruline » 2017
25 Li R. Debella H.J. Carmichael W.W. Isolates identifiable as Arthrospira maxima and Arthrospira fusiformis (Oscillatoriales, Cyanobacteria) appear identical on the basis of a morphological study in culture and 16S rRNA gene sequences Phycologia 40 2001 367 371 10.2216/i0031-8884-40-4-367.1
26 Parte A.C. Sardà Carbasse J. Meier-Kolthoff J.P. Reimer L.C. Göker M. List of Prokaryotic names with Standing in Nomenclature (LPSN) moves to the DSMZ Int. J. Syst. Evol. Microbiol. 70 2020 5607 5612 10.1099/ijsem.0.004332 32701423
27 Scheldeman P. Baurain D. Bouhy R. Scott M. Mühling M. Whitton B.A. Belay A. Wilmotte A. Arthrospira (‘Spirulina’) strains from four continents are resolved into only two clusters, based on amplified ribosomal DNA restriction analysis of the internally transcribed spacer FEMS Microbiol. Lett. 172 1999 213 222 10.1111/j.1574-6968.1999.tb13471.x 10188250
28 Choi G.-G. Ahn C.-Y. Oh H.-M. Phylogenetic relationships of Arthrospira strains inferred from 16S rRNA gene and cpcBA-IGS sequences ALGAE 27 2012 75 82 10.4490/algae.2012.27.2.075
29 Miklaszewska M. Waleron M. Morin N. Calusinska M. Wilmotte A. Tandeau De Marsac N. Rippka R. Waleron K. Elucidation of the gas vesicle gene clusters in cyanobacteria of the genus Arthrospira (Oscillatoriales, Cyanophyta) and correlation with ITS phylogeny Eur. J. Phycol. 47 2012 233 244 10.1080/09670262.2012.692817
30 Papapanagiotou G. Gkelis S. Taxonomic revision of commercially used Arthrospira (Cyanobacteria) strains: a polyphasic approach Eur. J. Phycol. 54 2019 595 608 10.1080/09670262.2019.1624832
31 Roussel T. Halary S. Duval C. Piquet B. Cadoret J.-P. Vernès L. Bernard C. Marie B. Monospecific renaming within the cyanobacterial genus Limnospira (Spirulina) and consequences for food authorization J. Appl. Microbiol. 134 2023 lxad159 10.1093/jambio/lxad159
32 Jacquet S. Zhong X. Parvathi A. Ram A.S.P. First description of a cyanophage infecting the cyanobacterium Arthrospira platensis (Spirulina) J. Appl. Phycol. 25 2013 195 203 10.1007/s10811-012-9853-x
33 Stackebrandt E. Goebel B.M. Taxonomic Note: A Place for DNA-DNA Reassociation and 16S rRNA Sequence Analysis in the Present Species Definition in Bacteriology Int. J. Syst. Evol. Microbiol. 44 1994 846 849 10.1099/00207713-44-4-846
34 Meier-Kolthoff J.P. Auch A.F. Klenk H.-P. Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions BMC Bioinf. 14 2013 60 10.1186/1471-2105-14-60
35 Komarek J. Anagnostidis K. Süßwasserflora von Mitteleuropa, Bd. 19/2: Cyanoprokaryota 2005 Springer
36 Jourdan J.-P. Cultivez Votre Spiruline Edt Antenna Technol. 2006 1 146
37 Marter P. Huang S. Brinkmann H. Pradella S. Jarek M. Rohde M. Bunk B. Petersen J. Filling the gaps in the cyanobacterial tree of life—metagenome analysis of Stigonema ocellatum DSM 106950, Chlorogloea purpurea SAG 13.99 and Gomphosphaeria aponina DSM 107014 Genes 12 2021 389 10.3390/genes12030389 33803228
38 Shih P.M. Wu D. Latifi A. Axen S.D. Fewer D.P. Talla E. Calteau A. Cai F. Tandeau de Marsac N. Rippka R. Improving the coverage of the cyanobacterial phylum using diversity-driven genome sequencing Proc. Natl. Acad. Sci. 110 2013 1053 1058 10.1073/pnas.1217107110 23277585
39 Brinkmann H. Van der Giezen M. Zhou Y. Poncelin de Raucourt G. Philippe H. An empirical assessment of long-branch attraction artefacts in deep eukaryotic phylogenomics Syst. Biol. 54 2005 743 757 10.1080/10635150500234609 16243762
40 Richter M. Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition Proc. Natl. Acad. Sci. 106 2009 19126 19131 10.1073/pnas.0906412106 19855009
41 Goris J. Konstantinidis K.T. Klappenbach J.A. Coenye T. Vandamme P. Tiedje J.M. DNA–DNA hybridization values and their relationship to whole-genome sequence similarities Int. J. Syst. Evol. Microbiol. 57 2007 81 91 10.1099/ijs.0.64483-0 17220447
42 Jain C. Rodriguez-R L.M. Phillippy A.M. Konstantinidis K.T. Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries Nat. Commun. 9 2018 5114 10.1038/s41467-018-07641-9 30504855
43 Lee I. Ouk Kim Y. Park S.-C. Chun J. OrthoANI: an improved algorithm and software for calculating average nucleotide identity Int. J. Syst. Evol. Microbiol. 66 2016 1100 1103 10.1099/ijsem.0.000760 26585518
44 Willis A. Woodhouse J.N. Defining Cyanobacterial Species: Diversity and Description Through Genomics Crit. Rev. Plant Sci. 39 2020 101 124 10.1080/07352689.2020.1763541
45 Bouma-Gregson K. Crits-Christoph A. Olm M.R. Power M.E. Banfield J.F. Microcoleus (Cyanobacteria) form watershed-wide populations without strong gradients in population structure Mol. Ecol. 31 2022 86 103 10.1111/mec.16208 34608694
46 Wayne L.G. Moore W.E.C. Stackebrandt E. Kandler O. Colwell R.R. Krichevsky M.I. Truper H.G. Murray R.G.E. Grimont P.A.D. Brenner D.J. Report of the Ad Hoc Committee on Reconciliation of Approaches to Bacterial Systematics Int. J. Syst. Evol. Microbiol. 37 1987 463 464 10.1099/00207713-37-4-463
47 Meier-Kolthoff J.P. Carbasse J.S. Peinado-Olarte R.L. Göker M. TYGS and LPSN: a database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes Nucleic Acids Res. 50 2022 D801 D807 10.1093/nar/gkab902 34634793
48 Chun J. Oren A. Ventosa A. Christensen H. Arahal D.R. da Costa M.S. Rooney A.P. Yi H. Xu X.-W. De Meyer S. Trujillo M.E. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes Int. J. Syst. Evol. Microbiol. 68 2018 461 466 10.1099/ijsem.0.002516 29292687
49 Meier-Kolthoff J.P. Hahnke R.L. Petersen J. Scheuner C. Michael V. Fiebig A. Rohde C. Rohde M. Fartmann B. Goodwin L.A. Complete genome sequence of DSM 30083T, the type strain (U5/41T) of Escherichia coli, and a proposal for delineating subspecies in microbial taxonomy Stand. Genomic Sci. 9 2014 2 10.1186/1944-3277-9-2 25780495
50 Santos K.R.d.S. Hentschke G.S. Ferrari G. Andreote A.P.D. Fiore M.d.F. Vasconcelos V. Sant’Anna C.L. Molecular, morphological and ecological studies of Limnospira platensis (Cyanobacteria), from saline and alkaline lakes, Pantanal Biome, Brazil Front. Environ. Sci. 11 2023 1204787 10.3389/fenvs.2023.1204787
51 Mühling M. Belay A. Whitton B.A. Screening Arthrospira (Spirulina) strains for heterotrophy J. Appl. Phycol. 17 2005 129 135 10.1007/s10811-005-7214-8
52 Cavallini A. Torre S. Usai L. Casula M. Fais G. Nieri P. Concas A. Lutzu G.A. Effect of cheese whey on phycobiliproteins production and FAME profile by Arthrospira platensis (Spirulina): Promoting the concept of a circular bio-economy Sustain. Chem. Pharm. 40 2024 101625 10.1016/j.scp.2024.101625
53 Mussagy A. Annadotter H. Cronberg G. An experimental study of toxin production in Arthrospira fusiformis (Cyanophyceae) isolated from African waters Toxicon 48 2006 1027 1034 10.1016/j.toxicon.2006.08.013 17049366
54 Yang Y. Park Y. Cassada D.A. Snow D.D. Rogers D.G. Lee J. In vitro and in vivo safety assessment of edible blue-green algae, Nostoc commune var. sphaeroides Kützing and Spirulina plantensis Food Chem. Toxicol. 49 2011 1560 1564 10.1016/j.fct.2011.03.052 21473896
55 Mysara M. Njima M. Leys N. Raes J. Monsieurs P. From reads to operational taxonomic units: an ensemble processing pipeline for MiSeq amplicon sequencing data GigaScience 6 2017 1 10 10.1093/gigascience/giw017
56 Wang Q. Garrity G.M. Tiedje J.M. Cole J.R. Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy Appl. Environ. Microbiol. 73 2007 5261 5267 10.1128/AEM.00062-07 17586664
57 Tamura K. Stecher G. Kumar S. MEGA11: molecular evolutionary genetics analysis version 11 Mol. Biol. Evol. 38 2021 3022 3027 10.1093/molbev/msab120 33892491
58 Li D. Liu C.-M. Luo R. Sadakane K. Lam T.-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph Bioinformatics 31 2015 1674 1676 10.1093/bioinformatics/btv033 25609793
59 Wu Y.-W. Simmons B.A. Singer S.W. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets Bioinformatics 32 2016 605 607 10.1093/bioinformatics/btv638 26515820
60 Kang D.D. Froula J. Egan R. Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities PeerJ 3 2015 e1165 10.7717/peerj.1165
61 Sieber C.M.K. Probst A.J. Sharrar A. Thomas B.C. Hess M. Tringe S.G. Banfield J.F. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy Nat. Microbiol. 3 2018 836 843 10.1038/s41564-018-0171-1 29807988
62 Parks D.H. Imelfort M. Skennerton C.T. Hugenholtz P. Tyson G.W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes Genome Res. 25 2015 1043 1055 10.1101/gr.186072.114 25977477
63 Tanizawa Y. Fujisawa T. Nakamura Y. DFAST: a flexible prokaryotic genome annotation pipeline for faster genome publication Bioinformatics 34 2018 1037 1039 29106469
64 Jones P. Binns D. Chang H.-Y. Fraser M. Li W. McAnulla C. McWilliam H. Maslen J. Mitchell A. Nuka G. InterProScan 5: genome-scale protein function classification Bioinformatics 30 2014 1236 1240 10.1093/bioinformatics/btu031 24451626
65 Almagro Armenteros J.J. Tsirigos K.D. Sønderby C.K. Petersen T.N. Winther O. Brunak S. von Heijne G. Nielsen H. SignalP 5.0 improves signal peptide predictions using deep neural networks Nat. Biotechnol. 37 2019 420 423 10.1038/s41587-019-0036-z 30778233
66 Parks D.H. Chuvochina M. Rinke C. Mussig A.J. Chaumeil P.-A. Hugenholtz P. GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy Nucleic Acids Res. 50 2022 D785 D794 10.1093/nar/gkab776 34520557
67 Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies Bioinformatics 30 2014 1312 1313 10.1093/bioinformatics/btu033 24451623
68 Camacho C. Coulouris G. Avagyan V. Ma N. Papadopoulos J. Bealer K. Madden T.L. BLAST+: architecture and applications BMC Bioinf. 10 2009 421 10.1186/1471-2105-10-421
69 Blin K. Shaw S. Augustijn H.E. Reitz Z.L. Biermann F. Alanjary M. Fetter A. Terlouw B.R. Metcalf W.W. Helfrich E.J.N. antiSMASH 7.0: new and improved predictions for detection, regulation, chemical structures and visualisation Nucleic Acids Res. 51 2023 W46 W50 10.1093/nar/gkad344 37140036
70 Cogne G. Lehmann B. Dussap C.-G. Gros J.-B. Uptake of macrominerals and trace elements by the cyanobacterium Spirulina platensis (Arthrospira platensis PCC 8005) under photoautotrophic conditions: Culture medium optimization Biotechnol. Bioeng. 81 2003 588 593 10.1002/bit.10504 12514808
71 Morin N. Vallaeys T. Hendrickx L. Natalie L. Wilmotte A. An efficient DNA isolation protocol for filamentous cyanobacteria of the genus Arthrospira J. Microbiol. Methods 80 2010 148 154 10.1016/j.mimet.2009.11.012 20004220
72 Segers C. Mysara M. Claesen J. Baatout S. Leys N. Lebeer S. Verslegers M. Mastroleo F. Intestinal mucositis precedes dysbiosis in a mouse model for pelvic irradiation ISME Commun. 1 2021 24 10.1038/s43705-021-00024-0 36737646
73 Eren A.M. Maignien L. Sul W.J. Murphy L.G. Grim S.L. Morrison H.G. Sogin M.L. Freckleton R. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data Methods Ecol. Evol. 4 2013 1111 1119 10.1111/2041-210X.12114 24358444
74 Alneberg J. Bjarnason B.S. De Bruijn I. Schirmer M. Quick J. Ijaz U.Z. Lahti L. Loman N.J. Andersson A.F. Quince C. Binning metagenomic contigs by coverage and composition Nat. Methods 11 2014 1144 1146 10.1038/NMETH.3103 25218180
75 Methner A. Kuzyk S.B. Petersen J. Bauer S. Brinkmann H. Sichau K. Wanner G. Wolf J. Neumann-Schaal M. Henke P. Thiorhodovibrio frisius and Trv. litoralis spp. nov., Two Novel Members from a Clade of Fastidious Purple Sulfur Bacteria That Exhibit Unique Red-Shifted Light-Harvesting Capabilities Microorganisms 11 2023 2394 10.3390/microorganisms11102394 37894052
76 Na S.-I. Kim Y.O. Yoon S.-H. Ha S.M. Baek I. Chun J. UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction J. Microbiol. 56 2018 280 285 10.1007/s12275-018-8014-6 29492869
77 Philippe H. MUST, a computer package of Management Utilities for Sequences and Trees Nucleic Acids Res. 21 1993 5264 5272 10.1093/nar/21.22.5264 8255784
78 Talavera G. Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments Syst. Biol. 56 2007 564 577 10.1080/10635150701472164 17654362
79 Bashir F. Bashir A. Bouaïcha N. Chen L. Codd G.A. Neilan B. Xu W.-L. Ziko L. Rajput V.D. Minkina T. Cyanotoxins, biosynthetic gene clusters, and factors modulating cyanotoxin biosynthesis World J. Microbiol. Biotechnol. 39 2023 241 10.1007/s11274-023-03652-x 37394567
80 Gugger M. Boullié A. Laurent T. Cyanotoxins and Other Bioactive Compounds from the Pasteur Cultures of Cyanobacteria (PCC) Toxins 15 2023 388 10.3390/toxins15060388 37368689
