Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2007, p. 6543-6550, Vol. 73, No. 20
0099-2240/07/$08.00+0 doi:10.1128/AEM.01377-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Applied Chemistry and Microbiology, University of Helsinki, Viikki Biocenter, P.O. Box 56, 00014 University of Helsinki, Finland
Received 21 June 2007/ Accepted 20 August 2007
|
|
|---|
|
|
|---|
Recently, microcystin-LR was detected in low concentrations at the entrance to the Gulf of Finland (17). The toxin producer could not be identified, but it was suspected to be Anabaena (17). However, to date, no clear evidence has been presented that links microcystin production in the Baltic Sea directly to the genus Anabaena. It is currently believed that the toxicity of cyanobacterial blooms in the Baltic Sea is attributed exclusively to the genus Nodularia spumigena, which produces the hepatotoxic cyclic pentapeptide nodularin (40). In freshwater environments, Aphanizomenon flos-aquae produces neurotoxins (43) but, for the Baltic Sea, it has been reported to be nontoxic (37, 40).
Microcystin-producing Anabaena strains have been identified in freshwater environments in Canada, Denmark, Egypt, Finland, France, and Norway (39). Microcystins are specific inhibitors of eukaryotic serine/threonine protein phosphatases 1 and 2A (30) and have been linked to liver cancer in humans (20). Although the binding mechanisms of microcystins to protein phosphatases are documented in detail (11), the ecological function of these toxins is largely unknown. Microcystins are small cyclic peptides comprised of seven amino acids (cyclo-D-Ala1-X2-D-MeAsp3-Z4-Adda5-D-Glu6-Mdha7). Structural variation has been encountered in all seven positions, but the most variable amino acids are found at the X and Z positions. The X and Z positions contain L-amino acids, the most common of which are leucine and arginine, respectively, but other amino acids are also found (5). Position seven, not normally referred to as a variable residue, usually contains N-methyldehydroalanine, but seven other amino acids have been reported at this position (5). Position five is occupied by 3-amino-9-methoxy-2,6,8-trimethyl-10-phenyldeca-4,6-dienoic acid (Adda), an unusual beta-amino acid found only in microcystins and nodularins, and is believed to play a critical role in the toxicity of these compounds (38). Over 70 different microcystin variants have been described previously, while the number of new microcystins reported continues to grow (4). Microcystin variants differ in toxicity with acute 50% lethal dose values ranging from 50 to 1,200 µg kg–1 (mouse, intraperitoneally), although nontoxic variants also exist (39).
We isolated a range of planktonic Anabaena strains from the Gulf of Finland and screened the cultures for microcystins by using liquid chromatography-mass spectrometry (LC-MS). We identified a number of Anabaena strains producing microcystins. The results presented here unequivocally identify Anabaena as a microcystin producer in the Baltic Sea. Furthermore, our results demonstrate unexpected genetic heterogeneity of Baltic Sea Anabaena strains.
|
|
|---|
![]() View larger version (27K): [in a new window] |
FIG. 1. Locations of sampling stations in the Gulf of Finland. Open-sea Anabaena strains were isolated from stations 1 to 15 in 2004. Anabaena strains 315 and 318 were isolated from the coastal waters of Helsinki (the locations are indicated with white arrows). (Modified from reference 19a.)
|
|
View this table: [in a new window] |
TABLE 1. Number of Anabaena strains isolated from each sampling station and data for environmental parametersa
|
Water samples were collected by a rosette sampler from depths of 0, 3, and 7 m (1 liter in volume) to make a direct toxin analysis. Water samples were filtered through a 1.0-µm filter (47-mm polycarbonate filters; Osmonics, Minnetonka, MN) and subsequently frozen at 20°C. Water sample handling was performed as sample handling for the Anabaena strain was, but the water sample was additionally bath sonicated (Sonorex Super 10P; Bandelin Electronic, Berlin, Germany) for 15 min and mixing and sonication were repeated three times. Prior to microcystin analysis, the supernatant was concentrated by vacuum centrifuging (Heto vacuum centrifuge; Heto-Holten A/S, Allerød, Denmark) and reconstituted to 50 µl.
Microcystin analysis.
The analysis of microcystins was carried out using a high-performance liquid chromatograph, combined with a diode array detector (Agilent 1100; Agilent Technologies, Santa Clara, CA) and a mass spectrophotometer (Agilent XCT Plus ion trap). Microcystins from the methanolic extracts were separated with a Zorbax Eclipse XDB-C8 column (4.6 by 150 mm; particle diameter, 5 µm; Agilent Technologies, Santa Clara, CA). For liquid chromatography, the mobile phase consisted of a gradient of 0.1% aqueous (water for high-pressure liquid chromatography, CHROMASOLV Plus; Sigma-Aldrich, Steinheim, Germany) formic acid (50% solution, Fluka; Sigma-Aldrich, Steinheim, Germany) (solvent A) and 0.1% formic acid in acetonitrile (Sigma-Aldrich, Steinheim, Germany) (solvent B). The linear gradient was as follows: 20% solvent A at 0 min, 80% solvent A at 50 min, and 100% solvent A at 65 min. A flow rate of 0.6 ml min–1 was used with the column temperature set to 40°C. Electrospray ionization was performed in positive ion mode. Nebulizer gas (N2) pressure was 50 lb/in2, and drying gas flow and temperature were 10 liters min–1 and 350°C, respectively. The capillary voltage was set at 3,270 V, with a capillary exit offset value of 317.4 V, a skimmer 1 potential of 41.5 V, and a trap drive value of 82.8. Spectra were recorded as averages of four using ultra scan mode and a scan range from 50 to 1,200 m/z. MS2 spectra were recorded in an auto-MS mode by using the following parameters: 5 to 10 precursor ions from 800 to 1,200 m/z, an isolation width of 4.0 m/z, and a fragmentation amplitude value of 0.50 V.
Microcystin analysis of filtered water sample was carried out as described earlier, with the following exceptions: the separation column was a Luna C18(2) column (150 by 2.0 mm; particle diameter, 5 µm; Phenomenex, Torrance, CA), the solvent A concentration was 60% at 60 min, the gradient flow rate was 0.15 ml min–1, the nebulizer gas (N2) pressure was 30 lb/in2, and the drying flow was 8 liters min–1.
Anabaena strains 90 (9, 42) and 66A (32, 42) as well as Microcystis viridis strain NIES102 (18) were used as reference strains. The identification of the microcystins was based on characteristic MH+ values corresponding with the range of published microcystins and the loss of neutral fragment 134 in the ion source (an event not so clearly observed from the MS2 spectrum). In addition, the occurrence of high intensity ion m/z 599 [(MeAsp)-Arg-Adda-(Glu) with H+] and less intense ion m/z 375 or 361 [Adda-134-Glu-(M)dha with H+; not present in all MS2 spectra] in the MS2 spectrum as well as comparison of LC-MS properties of microcystins produced by reference strains were used for identification. Microcystins have a characteristic UV spectrum with a maximum absorbance of 238 nm, and in the case of major variants, the UV spectrum attested also to the presence of microcystins. Peak areas from the chromatographic MH+ ion signals of the methanolic extracts and the MCYST-LR standard (a gift from Zbigniew Grzonka, Faculty of Chemistry, University of Gda
sk, Poland) were used for microcystin quantitation.
Collecting environmental data from the sampling stations.
Water samples for the determination of phosphate, total phosphorus, silicate, nitrate, nitrite, ammonium, and total nitrogen concentrations were obtained with a rosette sampler, and measurements were carried out at the time of sampling by using an autoanalyzer and following the guidelines of the Baltic Sea Monitoring Programme (13). A conductivity-temperature-depth probe provided continuous measurements of temperature, oxygen, and salinity. However, conductivity-temperature-depth sampling was not conducted at station 2 (Table 1).
Statistical analysis.
Since all microcystin-producing Anabaena strains were detected at one sampling station, station 15 (see Results), we used principal component analysis (PCA) to relate the presence of the microcystin-producing Anabaena strains to environmental conditions. PCA analysis was performed with Canoco for Windows, version 4.52 (48), and the PCA plot was constructed in CanoDraw 4.5.
16S rRNA gene PCR amplification and sequencing.
DNA from the strains was extracted using a modified cetyltrimethylammonium bromide (CTAB)-based extraction method (19). The quality and quantity of extracted DNA were determined by UV spectrophotometric measurements. The 16S rRNA gene was amplified from Anabaena strains with the primer pair pA (6) and B23S (27) as previously described (10). The 16S rRNA gene was sequenced using the internal sequencing primers 16S544R, 16S1092R, and 16S979F (35). The sizes of the PCR amplification products were checked in agarose gels, and PCR products were purified using Montage PCR centrifugal filter devices (Millipore, Billerica, MA). Sequencing was performed in the DNA sequencing laboratory of the Institute of Biotechnology, University of Helsinki, Helsinki, Finland.
Phylogenetic analysis of 16S rRNA gene sequences.
Phylogenetic analysis was performed for microcystin-producing Anabaena strains in order to resolve their relationship to Anabaena strains originating from freshwater environments. Representative 16S rRNA gene sequences from toxic and nontoxic Anabaena spp. from a range of environments were aligned using ClustalW program, version 1.4, as implemented in the BioEdit sequence alignment editor (version 7.0.1). The alignment was refined manually as ambiguous sites were excluded. Altogether, 1,354 bp of the 16S rRNA gene was used for phylogenetic analyses. Trees were constructed using neighbor-joining (NJ), maximum parsimony (MP), and maximum likelihood (ML) algorithms in the PAUP, version 10b, program (47). The GTR+I+G evolutionary model of substitution was used for NJ and ML using parameters (base frequencies, rate matrix of substitution types, and shape of gamma distribution) estimated from the data. In MP and ML analyses, we used heuristic searches, tree bisection, and reconnection branch rearrangement with a rearrangement limit of 1,000. To assess the reliability of the tree constructions, 1,000 bootstrap resamplings were performed for each tree with the same heuristic search. Nodularia sp. strain PCC73104/1 (accession number AJ133184), Nodularia sp. strain PCC7804 (AJ133181), Cyanospira rippkae (AY038036), Anabaenopsis sp. strain PCC9215 (AY038033), and Anabaena cylindrica (AF091150) were used as outgroups in all analyses. Uncorrected P distances of 16S rRNA gene sequences were calculated in PAUP, version 10b.
Nucleotide sequence accession numbers.
The nucleotide sequences of 16S rRNA genes have been deposited in the GenBank database under the accession numbers EF547190 to EF547196.
|
|
|---|
Strains BIR246, BIR250A, BIR257, BIR258, and BIR260 as well as strain 318 produced the same set of microcystin variants (Fig. 2). Four dominant variants were detected, and they represented 80 to 90% of the total microcystin content (Fig. 2). The MS2 spectrum of variant A strongly indicated that the structure for A was [D-Asp3]MCYST-HtyR (Fig. 2 and 3). We identified variant B as MCYST-HtyR (Fig. 2 and 3), which has been documented from Anabaena flos-aquae NRC 525-17 (12) and is the only variant having an m/z value of 1,059. Variants C and D were classified as [D-Asp3]MCYST-LR and MCYST-LR, respectively (Fig. 2 and 3), because their fragmentation of the MH+ ions and retention times were in agreement with the same microcystins produced by reference strain Anabaena 90 (9, 42). In addition to these four main microcystin variants, 23 other minor variants were described for these six Anabaena strains (data not shown).
![]() View larger version (23K): [in a new window] |
FIG. 2. Dominant microcystin variants produced by seven Anabaena strains. The y axis represents the relative proportion (%) of each microcystin variant out of all variants of each strain produced. Strains BIR246, BIR250A, BIR257, BIR258, BIR260, and 318 produced the same dominant variants, whereas the microcystin pattern produced by strain 315 was different. Variant A, [D-Asp3]MCYST-HtyR; variant B, MCYST-HtyR, variant C, [D-Asp3]MCYST-LR; variant D, MCYST-LR; variant E, [D-Asp3, Dha7]MCYST-HtyR; and variant F, [Dha7]MCYST-HtyR.
|
![]() View larger version (33K): [in a new window] |
FIG. 3. Microcystin structures of the main variants produced by seven Anabaena strains.
|
Toxin analysis of the filtered water sample.
Nodularin and two variants of microcystins were detected by an LC-MS analysis of a filtered water sample from sampling station 15 (Fig. 4). The detected microcystin variants were MCYST-HtyR (m/z, 1,059) and MCYST-LR (m/z, 995). The total microcystin concentration in the water sample was 0.02 µg liter–1.
![]() View larger version (13K): [in a new window] |
FIG. 4. LC-MS chromatogram of the filtered water sample. Peaks representing nodularin, MCYST-HtyR, and MCYST-LR are indicated with black arrows. The y axis represents the intensity of each peak, and the x axis represents the retention time (in minutes).
|
![]() View larger version (16K): [in a new window] |
FIG. 5. PCA plot based on the environmental parameters measured from sampling stations (values listed in Table 1). Numbers correlate to the station numbering in Fig. 1 as well as in Table 1 (station 1 was sampled twice [indicated in the PCA plot as 1a and 1b]). Low-salinity and low-silicate concentrations were the variables to explain the divergence of station 15 (station indicated with a square) from the other stations. All open-sea microcystin-producing Anabaena strains were isolated from the station 15. The first and second principal components accounted for 94.4% of the total variance.
|
![]() View larger version (56K): [in a new window] |
FIG. 6. Neighbor-joining tree based on 16S rRNA gene sequences (1,354 bp). Four branches, including the microcystin-producing Anabaena strain isolated from the Gulf of Finland, are boxed. Numbers at the nodes indicate bootstrap values of more than 50% for NJ, MP, and ML analysis. The taxa Nodularia sp. strain PCC73104/1 (accession number AJ133184), Nodularia sp. strain PCC7804 (AJ133181), Cyanospira rippkae (AY038036), Anabaenopsis sp. strain PCC9215 (AY038033), and Anabaena cylindrica (AF091150) were used as outgroups in all analyses. , hepatotoxic strains studied by LC-MS.
|
|
View this table: [in a new window] |
TABLE 2. Uncorrected P similarities of 16S rRNA gene sequences (1,354 bp) of the Baltic Anabaena strains
|
|
|
|---|
All seven toxic Anabaena strains produced two to four dominant microcystin variants. It is common for Anabaena to produce a number of main toxin variants simultaneously (41). The 50% lethal dose (LD50) toxicity values for [D-Asp3]MCYST-HtyR (variant A) (Fig. 2 and 3) and [D-Asp3]MCYST-LR (variant C) are 160 to 300 µg kg–1 (39). For MCYST-HtyR (variant B), the LD50 value ranges between 80 and 100 µg kg–1, and for MCYST-LR (variant D), the LD50 value is 50 µg kg–1 (39). Also, variants [D-Asp3,Dha7]MCYST-HtyR (variant E) and [Dha7]MCYST-HtyR (variant F) have proved to be toxic (39). Variants [D-Asp3]MCYST-LR and MCYST-LR are prevalent in Anabaena strains isolated from Finnish freshwater habitats (42), and MCYST-LR was detected also at the entrance to the Gulf of Finland (17) as well as at the coastal waters of the southern Baltic proper (28).
Two microcystin variants, MCYST-LR and MCYST-HtyR, together with the nodularin, were detected in the filtered water sample from station 15. The concentration of MCYST-LR was about twice as much as MCYST-HtyR, whereas MCYST-HtyR was produced in greater amounts than MCYST-LR was in most of the isolated strains. Anabaena strains also produced demethylated microcystin variants, which were not detected in the filtered water sample. Temperature might be one factor contributing to this difference because a rise in temperature has been reported to increase the demethylated microcystin variants produced in freshwater Anabaena strains (36). At the time of sampling, the temperature of the surface water column at station 15 was 15.8°C, whereas isolated strains were grown at 22 to 24°C in the laboratory.
All of the microcystin-producing Anabaena strains were isolated from the easternmost sampling station (station 15). Thus, it is possible that some environmental factor(s) at sampling station 15 favored microcystin-producing strains. It is generally accepted that cyanobacteria produce toxins in conditions that are high in the nutrients that usually also favor the growth of these organisms (39). The concentration of PO4 phosphorus, which is considered to promote microcystin production in freshwater Anabaena strains (36), was nevertheless lower at station 15 than at stations at the entrance to the Gulf of Finland (Table 1). The PCA demonstrated that low concentrations of silicate and salinity distinguished station 15 from other stations. However, cyanobacteria are not known to require silicate and the connection between the occurrence of microcystin-producing Anabaena and low silicate concentrations remains speculative. Based on the PCA analysis, salinity appears to be the best explanation for the confined presence of microcystin-producing Anabaena strains to the eastern part of the Gulf of Finland. Furthermore, the coastal waters of Helsinki City, where two toxic Anabaena strains (315 and 318) were isolated, had salinities lower than 5.7
in 1997 and 1998 (according to monitoring data of the Helsinki City Environment Centre). The Gulf of Finland has a strong salinity gradient from the western and more saline part to the less-saline eastern part, with salinity being affected by freshwater inflow (21). Salinity has been identified as a possible controlling factor for diazotrophic cyanobacterial blooms, e.g., for distribution (26) and genetic diversity (23) of Aphanizomenon in the Baltic Sea. Thus, we might speculate that low-salinity waters favor toxin-producing Anabaena. However, future studies are needed to assess the effect of different salinities on the fitness of the Baltic Sea Anabaena strains and their microcystin production.
Microcystin-producing Baltic Sea Anabaena strains were genetically heterogeneous, since studied strains dispersed to three different branches in the complex 16S rRNA gene tree. However, since the 16S rRNA gene similarities of the Baltic microcystin-producing Anabaena strains were at least 97.9% (Table 2), they all belonged to the same species according to the "species" description based on 16S rRNA gene similarity of >97.5% (44). Conversely, in the Baltic Sea plankton, hepatotoxic Nodularia spumigena (3, 24) as well as nontoxic Aphanizomenon flos-aquae (2, 23) populations have shown to be genetically homogeneous. Thus, the Anabaena population appears to be genetically more diverse than Aphanizomenon or Nodularia populations in the Baltic Sea.
In previous studies, microcystin-producing freshwater Anabaena strains have been genetically heterogeneous but have grouped together (including also a few nontoxic Anabaena strains) when distinguished by 16S rRNA genes (29, 34), internal transcribed spacer sequence 1-S (ITS1-S) (the spacer region of the ribosomal operon) (11), rpoB (RNA polymerase beta-subunit) (34), and rbcLX (RubisCO) (34) genes. Although the phylogeny of the Anabaena strains from freshwater environments has been studied intensively, the paucity of strains has impeded studies of the phylogeny of the brackish Anabaena. Based on our results, microcystin-producing Baltic Sea Anabaena strains are genetically more diverse than are freshwater hepatotoxic Anabaena strains in the 16S rRNA tree. It has been suggested that toxic and nontoxic cyanobacteria might be separated using proxies, such as 16S rRNA, to indirectly determine toxin production. For example, to identify Australian neurotoxic Anabaena circinalis strains, specific clustering of A. circinalis in rpoC1 (7) and 16S rRNA (1) gene trees has been used for detection. Similarly, hepatotoxic and nontoxic Nodularia strains have been distinguished by 16S rRNA genes as well as by hetR, ITS1-S, and PC-IGS loci (24, 25). Even recently, Janse et al. (15) suggested that the toxic Microcystis strains could be separated from the nontoxic ones by using rRNA ITS sequences (15). However, the unsuitability of molecular taxonomy to separate toxic and nontoxic strains has been stated previously especially in the cases of Microcystis (29, 33) and Planktothrix (29) and our results proved that unsuitability in case of the Anabaena strains.
To relate the toxin production unequivocally to a certain genus or species requires strain isolation or single-filament or aggregate picking. Altogether, only a few Anabaena strains have been isolated from the Baltic Sea (K. Rantasärkkä, unpublished data) and Anabaena has not been included in toxin screenings (e.g., see references 37 and 45), probably due to its role (considered minor) as a component of cyanobacterial blooms in the Baltic Sea. However, in the Gulf of Finland, Anabaena species may exist even as dominating components in phytoplankton, together with Aphanizomenon (17). Our results demonstrate that the assumption that Anabaena strains are nontoxic in the Baltic Sea is incorrect, and future toxin screenings should take into account the presence of microcystin variants produced by Anabaena, especially in the low-salinity coastal waters, where the recreational value is the greatest. Further work will be required to determine the distribution of microcystin-producing Anabaena in the Baltic Sea.
We thank Kerttu Koskenniemi and Kaisa Rantasärkkä for assistance on the R/V Aranda in collecting the strains. Kerttu Koskenniemi is also acknowledged for kindly preparing the map. We are indebted to Lyudmila Saari and Eeva Eronen, who purified and maintained the cultures. The critical comments of Hermanni Kaartokallio, Harri Kuosa, and Christina Lyra improved the final manuscript. The crew of R/V Aranda (vessel of the Finnish Institute of Marine Research) is greatly appreciated for their kind assistance.
Published ahead of print on 31 August 2007. ![]()
|
|
|---|
ski. 2003. Nodularia spumigena blooms and the occurrence of hepatotoxin in the Gulf of Gda
sk. Oceanologia 45:305-316.
tovská, R. Willame, A. Rantala, L. Hoffmann, J. Komárek, and K. Sivonen. 2005. Phylogenetic and morphological evaluation of the genera Anabaena, Aphanizomenon, Trichormus and Nostoc (Nostocales, Cyanobacteria). Int. J. Syst. Evol. Microbiol. 55:11-26.
milauer. 2002. Canoco reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5). Microcomputer Power, Ithaca, NY.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»