Next Article in Journal
Assessment of Fungal Contamination in Fish Feed from the Lake Victoria Basin, Uganda
Next Article in Special Issue
Insight into Unprecedented Diversity of Cyanopeptides in Eutrophic Ponds Using an MS/MS Networking Approach
Previous Article in Journal
Leukotoxins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cylindrospermopsin- and Deoxycylindrospermopsin-Producing Raphidiopsis raciborskii and Microcystin-Producing Microcystis spp. in Meiktila Lake, Myanmar

1
Norwegian Institute for Water Research, Gaustadalléen 21, N-0349 Oslo, Norway
2
Forest Research Institute, 15013 Yezin, Myanmar
3
Department of Natural Sciences and Environmental Health, University of South- Eastern Norway, Gullbringvegen 36, N-3800 Bø, Norway
4
Department of Sustainable Agro-ecosystem and Bioresources, Research and Innovation Centre, Fondazione Edmund Mach, Via E. Mach 1, 38010 San Michele all’Adige, Italy
5
National Research Council, 1411 Oxford Street, Halifax, NS B3H 3Z1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 12 March 2020 / Revised: 2 April 2020 / Accepted: 3 April 2020 / Published: 7 April 2020

Abstract

:
Meiktila Lake is a shallow reservoir located close to Meiktila city in central Myanmar. Its water is used for irrigation, domestic purposes and drinking water. No detailed study of the presence of cyanobacteria and their potential toxin production has been conducted so far. To ascertain the cyanobacterial composition and presence of cyanobacterial toxins in Meiktila Lake, water samples were collected in March and November 2017 and investigated for physico-chemical and biological parameters. Phytoplankton composition and biomass determination revealed that most of the samples were dominated by the cyanobacterium Raphidiopsis raciborskii. In a polyphasic approach, seven isolated cyanobacterial strains were classified morphologically and phylogenetically as R. raciborskii, and Microcystis spp. and tested for microcystins (MCs), cylindrospermopsins (CYNs), saxitoxins and anatoxins by enzyme-linked immunosorbent assay (ELISA) and liquid chromatography–mass spectrometry (LC–MS). ELISA and LC–MS analyses confirmed CYNs in three of the five Raphidiopsis strains between 1.8 and 9.8 μg mg−1 fresh weight. Both Microcystis strains produced MCs, one strain 52 congeners and the other strain 20 congeners, including 22 previously unreported variants. Due to the presence of CYN- and MC-producing cyanobacteria, harmful effects on humans, domestic and wild animals cannot be excluded in Meiktila Lake.
Key Contribution: This study confirmed the production of CYN and deoxyCYN by Raphidiopsis raciborskii strains and numerous MCs by Microcystis strains isolated from Meiktila Lake in Myanmar. The MCs included many novel congeners demonstrated by LC–MS and chemical derivatization methods. Among these were the rarely reported L-Glu and L-dihydrotyrosine-containing congeners. This is the first finding of toxin-producing cyanobacteria in a Myanmar waterbody.

1. Introduction

Many lakes and reservoirs worldwide are affected by periodic cyanobacterial dominance or even cyanobacterial blooms. Such mass developments of cyanobacteria are typical for eutrophic conditions and are often induced by nutrient enrichment caused by increased agricultural, urban and industrial activities and are also expected to increase due to regional and global climate change [1]. Various cyanobacterial species forming such blooms are potential producers of hepatotoxic or neurotoxic compounds and their presence is often associated with animal poisonings and a threat to human health [2].
Myanmar is characterized by the presence of several natural lakes and numerous man-made reservoirs. Meiktila Lake is one of the numerous reservoirs in Myanmar and was built in ancient times, dating from an unknown period [3] but most likely in the reign of King Narapathisithu (1173–1210) [4]. Today the lake is divided by a dam into a northern and a southern part (Figure 1) [5]. Meiktila Lake is exposed to sedimentation due to deforestation in the catchment and especially the northern part has been partially filled with sediment over a period of more than 100 years [6]. The priority use of water from Meiktila Lake is drinking water, water for domestic purposes and for irrigation, although the lake is also polluted with domestic waste water, street runoff and solid waste [6,7].
Only limited information is available about the limnological characteristics of Meiktila Lake and other freshwater habitats in Myanmar. In 1995, a study of the algal flora of Meiktila Lake was reported [4]. A recent study described the investigation of physical parameters, macrophyte and phytoplankton composition in the period 2011–2014 in Meiktila Lake [5]. Twenty taxa of aquatic macrophytes including helophytes have been documented in Meiktila Lake [5]. Several heterocytous cyanobacterial taxa, e.g., Anabaena, Anabaenopsis and Calothrix and a few nonheterocytous cyanobacterial taxa e.g., Aphanocapsa, Chroococcus, Microcystis, Arthrospira and Oscillatoria have been reported but not further investigated [4,5]. Neither study mentioned the presence of the cyanobacterium Raphidiopsis raciborskii (formerly Cylindrospermopsis raciborskii) (Woloszynska) Aguilera, Berrendero Gómez, Kastovsky, Echenique & Salerno nor the presence of the microcystin (MC) and cylindropsermopsin (CYN) groups of cyanobacterial toxins, which are documented from many lakes in Asia [8,9,10].
We suspect that the number of cyanobacterial species documented to date in Meiktila Lake was underestimated and that various toxin-producing cyanobacteria were present in the cyanobacterial community. There is clearly a lack of information about cyanobacteria and the production of cyanobacterial toxins in Meiktila Lake and the recently described potentially toxic cyanobacterium Microcystis is most likely not the only potential toxin-producing cyanobacterium. The ongoing pollution of the lake suggests the potential for occurrence of more frequent severe cyanobacterial blooms, which would have a negative impact on the use of the lake for drinking water and domestic purposes by the residents. This study aimed therefore to investigate the presence of cyanobacteria and their potential toxins in Meiktila Lake, applying modern analytical methods in a polyphasic approach to elucidate in detail the cyanobacterial composition, phylogeny and toxin production and toxin profiles.

2. Results

2.1. Physico-Chemical Parameters

At both sampling dates in March and November 2017, all sampling points in Meiktila Lake were characterized by water temperatures of 26.2–28.0 °C, pH of 8.5–9.3 and conductivities of 580–729 μS cm−1. Secchi depth was between 0.8 m (at sampling point MK1) and 1.8 m (at MK3). Total phosphorus and total nitrogen concentrations were 12–23 and 360–570 μg L−1, respectively.

2.2. Phytoplankton Community

At all sampling stations, and at both sampling dates, cyanobacteria were the dominant group in the phytoplankton in both parts of Meiktila lake together with diatoms (Bacillariophyceae) Cryptophyceae, Chlorophyceae and Euglenophyceae (Table 1). The most dominant cyanobacterium was R. raciborskii, which comprised biomasses between 0.2 and 1.9 mg L−1 fresh weight (FW), or 27%–91% of the cyanobacterial biomass at the sampling points MK1, MK2, MK4 and MK5. Other cyanobacteria present in the samples belonged to the genera Aphanocapsa, Aphanothece, Chroococcus, Merismopedia, Limnothrix, Microcystis, Planktolyngbya, Planktothrix, Sphaerospermopsis and Synechococcus. They together comprised biomasses between 0.06 and 1.2 mg L−1 cyanobacterial wet weight. At MK3, the biomass of R. raciborskii at both sampling dates (0.02–0.08 mg L−1) was lower than at the other sampling points MK1, MK2, MK4 and MK5 (0.20–1.94 mg L−1) (data not shown). The Microcystis biomass was lower than the Raphidiopsis biomass at all sampling points and sampling dates and ranged from 0.003 to 0.16 mg L−1 (data not shown).

2.3. Morphological and Phylogenetic Characterization

Seven potentially toxin-producing cyanobacterial strains were isolated from Meiktila Lake (Table 2). Based on morphological features, e.g., presence and form of colonies or filaments, vegetative cells and heterocytes, five of the isolated cyanobacterial strains were identified as R. raciborskii and two strains as M. aeruginosa and M. novacekii, respectively (Figure 2). However, Harke et al. [11] suggested all Microcystis warrant placement into the same species complex. Therefore, we use “Microcystis” instead of species names in the following parts of the manuscript.
The Raphidiopsis strains were mostly characterized by straight tapered filaments. The filament length and width varied between 8.8–90 × 1.9–5.8 µm. Heterocytes were observed in some filaments of all isolated strains. Akinetes were not observed in any of the investigated strains. As in the cultured strains, only a few of the filaments possessed heterocytes in the environmental samples. The two Microcystis strains were characterized by cell diameters ranging from 4.2 to 6.6 µm (strain AB2017/15) and from 3.7 to 5.8 µm (strain AB2017/14) (data not shown).
The morphological determination of the isolated strains was supported by phylogenetic analyses (Figure 3; Figure 4). Phylogenetic relationships of the investigated strains are presented in the maximum-likelihood (ML) tree of the 16S rRNA gene of Cylindrospermopsis/Raphidiopsis (Figure 3) and a separate ML tree of the Microcystis 16S rRNA gene (Figure 4). In the ML tree in Figure 3, the Raphidiopsis strains from Meiktila Lake grouped together with 16S rRNA gene sequences derived from Cylindrospermopsis and Raphidiopsis strains from Asia, Europe, Africa, Australia and North America (cluster I). The CYN-producing and nonCYN-producing Raphidiopsis strains from Meiktila Lake could not be distinguished phylogenetically using 16S rRNA gene and had similar 16 rRNA gene sequences (Figure 3). In cluster II, strains from North and South America, (USA, Mexico, Brazil), North Africa (Tunisia), Southwest Europe (Spain) and New Zealand, grouped together. Both Microcystis strains from Meiktila Lake possessed similar 16S rRNA gene sequences and clustered together with 16S rRNA gene sequences of Microcystis from Europe, Asia, Africa and South America (Figure 4).

2.4. Identification of Cyanobacterial Toxins and Toxin-Producing Strains

Three of the five investigated Raphidiopsis strains produced CYNs in variable amounts by either enzyme-linked immunosorbent assay (ELISA) or liquid chromatography with tandem mass spectrometry (LC-MS/MS) (Table 3). Concentrations of CYNs were 1.8–4.3 µg mg−1 FW by ELISA. Using LC–MS/MS, CYN concentrations of 1.7–2.5 µg mg−1 FW and deoxyCYN from 1.3 to 7.3 µg mg−1 FW were detected. In the three CYN-producing strains, deoxyCYN comprised 43%–75% of the total CYNs.
All investigated Raphidiopsis strains tested negative for saxitoxins (STXs), anatoxins (ATXs) and MCs by ELISA. Both Microcystis strains tested negative for CYNs, STXs and ATXs by ELISA but were identified as MC-producers by ELISA and their MC profiles were therefore investigated by high resolution LC–MS/MS (LC–HRMS/MS).
Underivatized samples were analysed by LC-HRMS/MS in positive and negative ionisation modes as previously described [12,13] and then after reaction with mercaptoethanol (targeting Mdha7/Dha7 moieties in MCs) [14] and Oxone/DMSO (targets sulfide groups in methionine and Cys/GSH conjugates of MCs) [15]. Results of these analyses are summarised in Figure 5, Table 4 and Table S1. Peaks from putative MCs were identified by their reaction with mercaptoethanol, production of characteristic product ions in data-dependent and/or data-independent acquisition (DDA and/or DIA) LC-MS/MS screens, and possessing plausible potential elemental formulae based on both positive and negative mode full scan HRMS. These peaks were then targeted by LC-HRMS/MS at suitable collision energies to obtain structurally informative HRMS/MS spectra to assist with identification, compared by LC–HRMS with samples containing some of the putative MCs, and subjected to selective oxidation to detect the presence of sulfide moieties that could be present in some of the MCs.
Peaks were only considered to be MCs if they: 1, showed apparent pseudo-molecular ions appropriate for a MC in both positive and negative ionisation modes; 2, showed one or more of the characteristic MC fragments shown in Figure 5; 3, displayed appropriate chemical reactivity for the putative structure, and; 4, displayed retention times (tR) and charge states (z) appropriate to the putative structure (e.g., based on the apparent number of polar and charged residues, such as Arg). MCs were considered “confirmed” (1, 3, 13, 14, 17, 18, 21, 25, 26 and 41) if they behaved identically in all respects to the standards (Table 4). Structures were considered “probable” if they behaved identically in all respects to a compound already identified with high probability in an available sample (2, 4, 8, 11, 12, 20, 28, 29, 36, 45 and 50). For compounds for which standards or appropriate samples were not available, these were regarded as “probable” if, in addition to displaying the appropriate physical and chemical characteristics (Table 4 and Table S1), they also displayed interpretable MS/MS spectra that were clearly consistent with the proposed structure by comparison with related compounds (57, 9, 10, 15, 16, 24, 27, 3032, 3740, 4244, 4649, 51, 53 and 54). Compounds were considered tentative if there was limited MS/MS spectral evidence (23) or if the evidence was ambiguous (e.g., several isomers were present that showed indistinguishable MS/MS spectra, i.e., 52, 55 and 56). Compounds designated “unidentified” were definitively identified as MCs, but the spectral data was insufficient to identify them (19, 22 and 3335). All compounds listed in Table 4 as containing Mdha7 and which gave adequate signal-to-noise in their MS/MS spectra in positive mode, showed product ions at m/z 135.0804, 375.1914 and 446.2286, indicative of the presence of Adda5D-Glu6–Mdha–D-Ala1, while in those listed as containing Mser7 or Dha7 the latter two product ions were heavier, or lighter, by a mass corresponding to H2O or CH2, respectively (all with Δm < 5 ppm), and the presence of these units is implicit in the discussion of the structural elucidation in Section 3.
In culture AB2017/14, 52 microcystin congeners were detected by LC–HRMS, with a total concentration of 1100 µg g−1 FW. Twenty-one of these were unidentified or previously unreported variants. In culture AB2017/15, 20 microcystin variants (of which six were previously unreported) were detected, with a total concentration of 14000 µg g−1 FW. The microcystin variants and the concentrations found in each strain are shown in Table 4 and Table S1.

3. Discussion

This study clearly demonstrates for the first time the presence of CYN- and deoxyCYN-producing R. raciborskii and MC-producing Microcystis in the phytoplankton community of Meiktila Lake in Myanmar. The relatively high biomass of R. raciborskii, up to 1.9 mg L−1 in the phytoplankton community of Meiktila Lake, is expected to cause elevated concentrations of CYNs in the lake water. The results suggest a higher risk for humans and animals to be affected by CYNs than by MCs, although this could be affected by variations in the biomass of CYN-producing R. raciborskii versus MC-producing Microcystis. Variations in cyanobacterial bloom composition and toxin production are influenced by abiotic factors such as nutrients, temperature and light and by biotic factors such as grazing, parasitism and predation [16,17]. The distribution of CYN/deoxyCYN-producing and nonproducing Raphidiopsis strains in Meiktila Lake is likely to vary over time, and dominance by a Raphidiopsis strain such as AB2017/13 would lead to CYN and deoxyCYN concentrations up to 20 µg L−1 for the highest Raphidiopsis biomasses measured in this study. It is therefore expected that CYN/ deoxyCYN concentrations in the lake water will at times exceed the guideline value for CYN in drinking water of 1 µg L−1 [18]. The tolerable daily intake value of 0.03 µg kg−1 for a person of 70 kg body weight would be exceeded after the intake of slightly more than 100 mL of lake water if CYN and deoxyCYN are similarly toxic. However, the toxicity of deoxyCYN to humans is not yet clear. According to Norris et al. [19], deoxyCYN does not contribute significantly to the toxicity of R. raciborskii. In contrast, cell viability assays showed that deoxyCYN was only slightly less toxic than CYN and most likely operates by similar toxicological mechanisms [20]. The potential risk of deoxyCYN for humans needs therefore to be clarified [20]. The use of Meiktila Lake water for drinking water, irrigation, domestic purposes or animal consumption is complicated by the fact that an unknown proportion of CYNs can be extracellular and is therefore not eliminated by filtration. Lake water contaminated by CYN and other toxins like MC-LR can lead to morphological and physiological changes and potential loss of productivity by agricultural plants, and bioaccumulation of cyanotoxins in the tissues of edible terrestrial plants in a concentration-dependent manner has been reported [21].
Griffith and Saker [22] have shown that in stationary phase of cultures, more than 50% of CYN can be extracellular. In environmental samples, the same authors found that extracellular CYN could exceed 90%. Boiling in water does not significantly degrade CYN within 15 min [23]. The removal of extracellular CYN/deoxyCYN therefore needs other methods, like the use of activated carbon, membrane filtration or chemical inactivation (Ultraviolet (UV), or oxidants) [24]. The presence of Raphidiopsis strains in Meiktila Lake that do not produce CYN makes it likely that CYN concentrations in the lake will vary considerably depending on the ratio of the two chemotypes in the phytoplankton community. As Meiktila Lake water is used for domestic purposes (drinking water, irrigation, washing of clothes, and personal hygiene), regular monitoring of cyanobacterial biomass and CYNs is recommended.
R. raciborskii has not been described from Myanmar water bodies and was not observed in a phytoplankton community study conducted in Meiktila Lake from 2011 to 2012 [5]. Raphidiopsis (and Cylindrospermopsis) spp., however, have been described from various other Southeast Asian freshwater habitats [25,26,27] and other water habitats worldwide [8]. R. raciborskii is only known to produce CYNs in Australia and the Asian countries China, Japan, Vietnam and Thailand and to produce STXs in Brazil [26,28,29,30,31,32]. The prime radiation centre of R. raciborskii is thought to be in Africa, with a second radiation centre in Australia [8]. Our 16S rRNA gene analysis confirms the close relationship of the Raphidiopsis strains from Meiktila Lake to other Raphidiopsis strains from Asia, Europe and Australia. Our 16S rRNA gene tree also clearly supports the suggested movement of Raphidiopsis from the American continents to Southwest Europe and North Africa and probably further to Greece and China, as has been described [33,34]. The close relationship of Raphidiopsis strains from Australia, Asia and Europe does not, however, explain why CYN- and deoxyCYN-producing strains have only been found in Asia and Australia but not in Europe, Africa or the Americas. Parts of, or the whole, CYN gene cluster could have been lost during the spread from Asia westwards, or only nonCYN-producing strains may have spread to Europe. Our finding of nontoxic Raphidiopsis strains in Meiktila Lake supports the latter hypothesis.
Both Microcystis strains AB 2017/14 and AB2017/15 isolated from Meiktila Lake are confirmed microcystin producers and are closely related to Microcystis strains from Africa, Europe and Asia based on 16S rRNA gene phylogeny. Both strains had identical 16S rRNA gene sequences but were clearly distinguished chemically by their MC congener profiles. Fifty-six microcystin variants were found in Microcystis strains AB2017/14 and AB2017/15 isolated from Meiktila Lake (Figure 6).
In order to reliably estimate the quantities of the MCs in the extracts by LC-HRMS, it was necessary to characterize, and if possible, identify all of them. The reason for this is that the response in LCMS can be expected to vary from congener-to-congener, primarily due to variations in the number of easily ionisable amino acid residues, especially Arg, present in the MC’s structure. Only the identification of previously unreported MCs (see Bouaïcha et al. [35]) in the cultures is discussed further (i.e., 7, 9, 19, 22, 24, 27, 3035, 37, 38, 40, 42, 46, 5153, 55 and 56) but spectra of all compounds for which adequate MS/MS spectra were obtained are available in the Supporting Information.
Three of the compounds were sulfide-containing variants (5, 15, and 31) which reacted when the extract was oxidised with Oxone/DMSO (Table 4). The first two have been reported in cultures and blooms [15,36], and their characteristics were fully consistent with those reported for 5 and 15 here, and in the case of 15 its oxidation product (MC-M(O)R) showed characteristic product ions including neutral loss of CH4OS and displayed an MS/MS spectrum (Figures S8, S11 and S14) identical to that reported previously for MC-M(O)R [15]. The third of sulfide-containing MC was identified as MC-RM (31) based on its physical and chemical properties (Table 4), which were essentially identical to those of 15 except for its longer tR and that its MS/MS spectrum closely paralleled that of MC-RA (28) and displayed product ions characteristic of an MC with one Arg at position-2 rather than at the more common position-4 (Figures S8, S14 and S15). For example, fragments at m/z 440.2263 (C18H30O6N7+, Δm = 2.2 ppm, from Mdha7D-Ala1–Arg2D-Masp3) and 731.3716 (C34H51O10N8+, Δm = −0.9 ppm, from Adda5D-Glu6–Mdha7D-Ala1–Arg2D-Masp3), together with the complete absence of a product ion at m/z 599.3552 (from Arg4–Adda5D-Glu6), confirmed Arg at position-2 and Met at position-4 (see Okello et al. [37] for assigned product ions from MC-YR) of 31.
Eight of the compounds (4, 9, 24, 27, 3638 and 40) showed characteristics of MCs containing one or more Glu residues at position-2 or -4. Two of these (9 and 27) had formulae consistent with MC-RE or MC-ER (Table 4). Compound 9 gave product ions (Figures S2 and S5–S7) typical of a MC with Arg at position-4, including m/z 599.3536 (C31H47O6N6+, Δm = −2.6 ppm, from Arg4–Adda5D-Glu6), 284.1238 (C12H18O5N3+, Δm = −1.1 ppm, from Mdha7D-Ala1–Glu2), and 286.1497 (C11H20O4N5+, Δm = −4.3 ppm, from D-Masp3–Arg4), indicating 9 to be MC-ER. The MS/MS spectrum of 27 (Figures S2 and S15) included product ions at m/z 440.2248 (C18H30O6N7+, Δm = −0.9 ppm, from Mdha7D-Ala1–Arg2D-Masp3), and 731.3678 (C34H51O10N8+, Δm = −6.1 ppm, from Adda5D-Glu6–Mdha7D-Ala1–Arg2D-Masp3) which, together with the complete absence of a product ion at m/z 599.3552, confirmed Arg at position-2 and Glu at position-4 (see Okello et al. [37]) for assigned product ions from MC-RY), showing 27 to be MC-RE. The characteristics of 37 (Table 4) were consistent with MC-EE. In addition, 37 gave product ions (Figures S27–S33) including m/z 276.1189 (C10H18O6N3+, Δm = −0.4 ppm, from D-Masp3–Glu4), 405.1605 (C15H25O9N4+, Δm = −2.9. ppm, from Glu2D-Masp3–Glu4) and 575.2703 (C28H39O9N4+, Δm = −1.5 ppm, from Adda5D-Glu6–Mdha7D-Ala1–Glu2) that confirmed 37 as MC-EE. Compound 40 was identified as MC-LE based on the characteristics presented in Table 4, as well as product ions (Figures S27–S33) observed in its MS/MS spectra, including m/z 460.2397 (C19H34N5O8+, Δm = −1.2, from D-Ala1–Leu2-D-Masp3–Glu4) and 397.2073 (C18H29N4O6+, Δm = −2.2, from Mdha7D-Ala1–Leu2-D-Masp3), confirmed its identity as MC-LE (40). Earlier-eluting desmethylated D-variants of 9, 27, 37 and 40 were similarly identified as the corresponding D-Asp3-congeners [D-Asp3]MC-ER (4), [D-Asp3]MC-RE (24), [D-Asp3]MC-EE (36) and [D-Asp3]MC-LE (38) based on analysis of their MS/MS spectra (Figures S21–S27) and characteristics presented in Table 4. Furthermore, the 4 and 36 in this sample coeluted with, and gave identical product ion spectra to, [D-Asp3]MC-ER (4) and [D-Asp3]MC-EE (36) identified [12] in an extract of a culture of Planktothrix prolifica NIVA-CYA544.
Twenty-one conventional late-eluting nonArg-containing MCs (30, 32 and 3656) were detected. Of these, the identities of four that contained Glu2 or Glu4 (3638 and 40) were discussed above. The remaining previously unreported nonArg MCs were 30, 32, 42, 46, 5153, 55 and 56. Compound 42 had the same characteristics as MC-LAbu (45) (Table 4), however, its MS/MS spectrum (Figures S34–S36) was consistent with MC-HilA. In particular, product ions at m/z 573.3270 (C30H45N4O7+, Δm = −2.3 ppm, from Adda5D-Glu6–Mdha7D-Ala1–Hil2 minus C9H10O (cf m/z 559.3126 for 41 and 45)) and 411.2231 (C19H31N4O6+, Δm = −1.6 ppm, from Mdha7D-Ala1–Hil2-D-Masp3 (cf m/z 397.2082 for 41 and 45)) as well as a range of other ions indicated the identity as MC-HilA (42), although the actual connectivity of the carbons in the amino acid side-chain at position-2 cannot be determined by mass spectrometry. A related compound (46) had characteristics (Table 4) and gave product ions (Figures S40–S42) that were consistent with MC-HilAbu. Product ions included m/z 573.3260 (C30H45N4O7+, Δm = −3.9, from Adda5D-Glu6–Mdha7D-Ala1–Hil2 minus C9H10O), 232.1291 (C9H18N3O4+, Δm = −0.4, from D-Masp3–Abu4) and 430.2651 (C19H36N5O6+, Δm = −2.1, from D-Ala1–Hil2-D-Masp3–Abu4 (cf m/z 402.2347 for 41 and 416.2504 42). This data unambiguously shows the presence of an extra CH2 group in both amino acid-2 and -4 in 46, relative to MC-LA (41) and is consistent with MC-HilAbu (46). Compound 51 displayed characteristics consistent with MC-FV (Table 4), as well as product ions (Figures S49–S51) at m/z 246.1456 (C10H20N3O4+, Δm = −0.4, from D-Masp3–Val4), 593.2960 (C32H41N4O7+, Δm = −1.6, from Adda5D-Glu6–Mdha7D-Ala1–Phe2 minus C9H10O (cf. 559.3126 for 41)) and 464.2494 (C22H34N5O6+, Δm = −3.1, from D-Ala1–Phe2D-Masp3–Val4 (cf. 402.2347 for 41)). This establishes an extra C5H2 and 4 RDBE in amino acid-2 and C2H4 in amino acid-4, relative to MC-LA (41), consistent with MC-FV (51). Compound 53 had characteristics consistent with MC-WV (Table 4). This was supported by its MS/MS spectra (Figures S49–S51), which included product ions at m/z 246.1456 (C10H20N3O4+, Δm = 2.9, from D-Masp3–Val4), 632.3063 (C34H42N5O7+, Δm = −2.5, from Adda5D-Glu6–Mdha7D-Ala1–Trp2 minus C9H10O (cf. 559.3126 for 41)) and 503.2605 (C24H35N6O6+, Δm = −2.6, from D-Ala1–Trp2D-Masp3–Val4 (cf. 402.2347 for 41)). This indicates the presence of an extra C5HN and 6 RDBE in amino acid-2 and C2H4 in amino acid-4, relative to MC-LA (41), consistent with MC-WV (53). Later-eluting isomers of 50, 51 and 54 were also present (i.e., 52, 55 and 56), with identical characteristics (Table 5) and product ion spectra (Figures S52–S56). These compounds all contain branching amino acids at the variable position-2 (nominally Leu for 52 and 56) or -4 (Val for 55), and most likely the isomers present result from changes to this branching (e.g., Ile or 2-aminohexanoic acid at position-2, and 2-aminopentanoic acid or isovaline at position-4).
In addition, 30 and 32 differed from each other by CH2 and had characteristics consistent with MC-(H2)YA (32) and MC-(H2)YG (30), respectively (Table 4). These fragmented somewhat differently from typical Arg-free MCs such as MC-LA (41) (Figure S17). Both compounds showed weak product ions at m/z 155.0815 and 580.3017, indicating that 30 and 32 both contained Adda5D-Glu6–Mdha7D-Ala1. However, both 30 and 32 also gave product ions at m/z 320.1605, 611.3075 and 745.3807 (cf. 268.1650, 559.3117 and 693.3854 in MC-LA (41)), indicating the presence of an extra C3O and 3 RDBE at amino acid-2 relative to 41, consistent with the presence of the unusual amino acid L-dihydrotyrosine ((H2)Y) at position-2. Compounds 30 and 32 gave product ions as m/z 449.2033 and 449.2015 (C21H29N4O7+, Δ 0.5 and −3.6 ppm, respectively, from Mdha7D-Ala1–(H2)Tyr2D-Masp3; cf. m/z 397.2082 for 41, from Mdha7D-Ala1–Leu2D-Masp3). Thus, the difference in mass (14.0157, i.e., CH2) between 30 and 32 lies not in residue-3 (D-Masp3 vs D-Asp3) as might be expected but in residue-4. Thus, 32 is identified as MC-(H2)YA, and 30 as MC-(H2)YG, which appears to be the first MC so far reported [35] with Gly at position-4.
Compounds 6 and 7 coeluted but differed by a mass corresponding to CH2, with 6 having the same accurate mass, MS/MS spectrum (Figure S3) and retention time as MC-(H2)YR (6) identified in a sample from a recent study [13]. The MS of 7 was consistent with [D-Asp3]MC-(H2)YR, and the MS/MS spectra of 6 and 7 were very similar. The MS/MS spectrum of 7 included product ions (Figure S4) at m/z 155.0813 (Δm = −1.6 ppm, from Mdha7D-Ala1) and 599.3525 (C31H47O6N6+, Δm = −4.5 ppm, from Arg4–Adda5D-Glu6), 120.0806 (C8H10N+, Δm = −1.6 ppm, from (H2)Tyr), and 320.1611 (C16H22O4N3+, Δm = 2.0 ppm, from Mdha7D-Ala1–(H2)Tyr2), 272.1353 (C10H18O4N5+, Δm = 2.0 ppm, from D-Asp3–Arg4) and 714.3802 (C35H52O9N7+, Δm = −2.7 ppm, from D-Asp3–Arg4–Adda5D-Glu6), showing that 7 is [D-Asp3]MC-(H2)YR.
Five MCs were present (19, 22 and 3335) whose structures could not be identified from their characteristics (Table 4 and Table S1) or product ion spectra (Figures S12 and S18–S20).Compounds 33 and 34 gave apparent m/z values that did not correspond to known or plausible MC variants and differed from each other by a mass corresponding to CH2. They did not contain Arg but contained one extra nitrogen atom, 5 extra RDBE and, more surprisingly, one less oxygen atom than MC-LA (41). Compound 33 gave product ions at m/z 375.1911 (C20H27O5N2+, Δm = −0.9 ppm, from Adda5D-Glu6–Mdha7 minus C9H10O), 446.2284 (C23H32N3O6+, Δm = −0.4 ppm, from Adda5D-Glu6–Mdha7D-Ala1 minus C9H10O) and 580.2986 (C32H42N3O7+, Δm = −5.4 ppm, from Adda5D-Glu6–Mdha7D-Ala1), however, fragments containing amino acids 2–4 were either shifted or absent. Although the data appear to be consistent with analogues containing Orn and Phe at positions 2 and 4 with an amide linkage to the neighbouring D-Asp3/Masp3 residue, further data are required for even a tentative structural assignment. The apparent pseudomolecular ion isotope envelopes of 33 and 34 (Figure S60) displayed unusual patterns that suggest that these compounds may be reacting during ionisation, possibly including dehydration, further complicating mass spectral analysis.
As with 30 and 32, compound 35 showed product ions at m/z 213.0866, 320.1583, 375.1909, 446.2279, 449.2013, 509.2680 and 611.3040, consistent with the presence of Adda5D-Glu6–Mdha7D-Ala1–(H2)Tyr2D-Masp3, which would require an MC with an amino acid side-chain at position-4 possessing an unprecedented C3H3 and 2 RDBE. This could possibly be due to the presence of a larger fragile amino acid at this position that undergoes ready elimination during MS, and the full structure of 35 remains undetermined.
Compounds 19 and 22 had identical product ion spectra including ions at m/z 135.0803, 213.0866, 269.1235, 375.1907, 446.2281, 599.3533, 640.2807 and 710.3886, identical to those from MC-LR, indicating the presence of D-Masp3–Arg4–Adda5D-Glu6–Mdha7D-Ala1 and that these compounds differed from MC-LR (17) and MC-WR (26) only in the amino acid at position-2. Product ions at 469.1815 and 486.2111 (both from X2D-Masp3–Arg4) and 173.0709 and 144.0455 (from X2 and 173.0709−CH3N) were consistent with this and indicated the presence of a side chain at position-2 containing C9H6NO and 7 RDBE (cf. C9H8N and 6 RDBE for the side chain of Trp in 26), suggesting the presence of an unidentified oxidised variant of Trp present in the two isomeric MCs, 19 and 22.
Between them, the two cultures contained 34 known and 22 previously unreported MCs of which 14 were assigned probable structures based on their chemical and mass spectrometric properties. In all cases, the number of Arg residues present in the MCs were reliably determined from their elemental compositions and charge-state, even for congeners for which definitive structures could not be established. Quantitation of the individual congeners was then performed from negative (Table 4) and positive mode full scan chromatograms, which gave essentially identical results, relative to 3-point calibration curves of the appropriate MC reference materials (RMs) containing no, one or two Arg residues ([D-Leu1]MC-LY [38], MC-LR (17) and MC-RR (3), respectively). AB2017/14 contained 52 MCs, of which only 31 had been previously reported, and the newly reported MCs constituted nearly 20% of the total identified MC content. AB2017/15 contained fewer MCs (20) of which only six had not been reported previously, but these constituted nearly 34% of the total MC content in this culture. Notable amongst the newly reported MCs are those containing Glu at positions 2 and 4 (4, 9, 24, 27, 3638, and 40), which constituted 1.7% of the MCs in AB2017/14 and 37.5% of the MCs in AB2017/15. MCs of this type have only been reported previously from two sources [12,39]. AB2017/14 also contained a number of congeners that appeared to be derived from MC-WR (19, 20, 22, 23 and 26) that together constituted 20.2% of its MC content, and an unusually high number (25) of nonArg-containing MCs (30, 3235 and 3756) together constituting 20.5% of the total MC content. These results underscore both the diversity of MCs that may be present in a single sample and the potential difficulty of reliably quantitating the total MCs using traditional targeted LC–MS/MS methods. These factors may contribute to the reported apparent overestimation of MC levels when using less-targeted methods such as ELISA and PP2A inhibition, relative to highly congener-targeted LC–MS/MS approaches [40].
Although both strains produce a variety of MC variants, the risk of harmful effects caused by microcystins is likely to be low due to the low Microcystis biomasses observed in Meiktila Lake during this study. An increase in Microcystis biomass cannot, however, be excluded due to the pollution from various sources, e.g., waste water and street runoff [6]. This will most likely lead to an increase in microcystin concentrations in the lake, with potential harmful effects on humans, domestic and wild animals using the untreated lake water. The toxicity of microcystins found in the strains isolated from Meiktila Lake varies from highly toxic variants like MC-LR or MC-LA to less toxic variants such as MC-RR [41]. The toxicity of most of the MC variants found in this study has not yet been described, which makes a full risk assessment difficult. The production of a high number of MC variants (up to 47) has also been shown for two other Microcystis strains isolated from South African Haartbeestpoort Dam and Japanese Lake Kasumigaura [42,43].
Shallow lakes like Meiktila Lake are often characterized by competition between macrophytes and phytoplankton. High nutrient loading and phytoplankton growth lead to turbid water conditions and prevent the growth of macrophytes [44]. In Meiktila Lake the present turbid conditions are explained by deforestation and erosion in the catchment area, pollution with domestic waste water, street runoff and solid waste [6]. The relatively strong growth of R. raciborskii most likely is an additional reason for the observed turbidity at sampling points MK1, MK2, MK4 and MK5 in Meiktila Lake. The Potamogeton belt which separates MK2 and MK3 seems to act as a kind of filter or barrier because the Secchi depth at MK3 (1.8 m) was considerably higher than at MK1 (0.8 m). According to Van Donk and Van de Bund [45], macrophytes significantly modify the composition of the phytoplankton community and lead to a decrease in its abundance and biomass. It has been shown that certain macrophyte species exhibit allelopathic activity against certain phytoplankton species [46]. The obvious decrease of Raphidiopsis biomass from MK2 to MK3 may be therefore attributable to allelopathic effects of the macrophyte species in the belt. Allelopathy has been suggested by several authors to be responsible for observed phytoplankton patterns in whole-lake studies of vegetated, shallow lakes, but evidence for or against allelopathy has not been provided [47,48,49].

4. Conclusions

In conclusion, this is the first report of R. raciborskii and Microcystis in Meiktila Lake, Myanmar. Three of five Raphidiopsis isolates produced CYN and deoxyCYN, like other Raphidiopsis/Cylindrospermopsis isolates described from other Asian countries and Australia. Both Microcystis strains isolated from Meiktila Lake produced at least 56 MC variants (52 for AB2017/14 and 20 for AB2017/15), including 22 previously undescribed congeners. Harmful effects on humans and animals using Meiktila Lake as a water source cannot therefore be excluded.

5. Materials and Methods

5.1. Study Area, Measurements and Sampling

Meiktila Lake is a shallow reservoir, located close to Meiktila city in central Myanmar in the Mandalay region at an altitude of 230 m (Figure 1). During the rainy season, from April/May to October/November, it receives water from Mondaing Dam, located ca. 15 km west of Meiktila at an altitude of 245 m. Depending on the season, Meiktila Lake covers an area of around 54 km2 with a maximum water depth of 10 m. It is divided by a dam into a northern and a southern lake [5]. Five sampling points were selected, three in the northern part of Meiktila Lake and two in the southern part (Figure 1, Table 5).
Sampling points 2 and 3 in the Northern Part of Meiktila Lake were separated by a broad Potamogeton belt. At all five sampling points in March and November 2017, in situ measurements of water temperature, pH, conductivity and dissolved oxygen were conducted, and integrated water samples were taken (1 m steps up to max. 3 m water depth) for the analysis of chemical parameters (ammonium, nitrate, total nitrogen, soluble reactive phosphorous, total phosphorous, Ca, phytoplankton composition and biomass and for the isolation of cyanobacterial strains. For quantitative phytoplankton analysis, a 50 mL subsample was removed from a sample taken from integrated samples and preserved with Lugol’s solution and a concentrated net sample (mesh size 20 µm) was taken and preserved by addition of formaldehyde (4% final concentration) for qualitative analysis. A 50 mL water sample for isolation of cyanobacteria was taken at each sampling point and kept in a cool shady place and gently shaken twice per day before further treatment in Norway.

5.2. Phytoplankton Analysis

The Lugol-fixed phytoplankton samples were counted in sedimentation chambers (Hydro-Bios Apparatebau GmbH Kiel, Germany) using an inverted microscope (Leica DMi8; Ortomedic, Oslo, Norway) according to Utermöhl [50]. Phytoplankton biomass was calculated by geometrical approximations using the computerized counting programme Opticount (SequentiX - Digital DNA Processing, Klein Raden, Germany). The specific density of phytoplankton cells was calculated as 1 g cm−3.

5.3. Isolation of Strains and Morphological Characterization

Using a microcapillary, single colonies of Microcystis and filaments of Raphidiopsis were isolated. They were washed five times and placed in wells on microtiter plates containing 300 µL Z8 medium [51]. After successful growth, the samples were placed in 50 mL Erlenmeyer flasks containing 20 mL Z8 medium and maintained at 22 °C. Strains were classified based on morphological traits [52,53]. Morphological examination was conducted using a Leica DM2500 light microscope, Leica DFC450 camera and Leica Application Suite software (LAS) (Leica, Oslo, Norway). The morphological identification was based on the following criteria: (i) size of vegetative cells and heterocytes and (ii) nature and shape of filaments or colonies. Length and width of 50–250 vegetative cells or filaments and of 20–50 heterocytes were measured. Akinetes were not detected in the samples. All strains used in this study (Table 2) are maintained at the Norwegian Institute for Water Research, Oslo, Norway.

5.4. Genomic DNA Extraction, PCR Amplification and Sequencing

Genomic DNA was extracted according to Ballot et al. [54]. All PCRs were performed on a Bio-Rad CFX96 Real-Time PCR Detection System (Bio-Rad Laboratories, Oslo, Norway) using the iProof High-Fidelity PCR Kit (Bio-Rad Laboratories, Oslo, Norway). The 16S rRNA gene of the isolated strains from Meiktila Lake was amplified using the primers as described by Ballot et al. [54]. PCR products were visualized by 1% agarose gel electrophoresis with GelRed staining (GelRed Nucleic Acid Gel Stain, Biotium, Fremont, CA, USA) and UV illumination.
Amplified 16S rRNA gene products were purified through Qiaquick PCR purification columns (Qiagen, Hilden, Germany). Sequencing of the purified 16S rRNA gene products was performed using the same primers as for PCR and intermediate primers as described in Ballot et al. [54]. For each PCR product, both strands were sequenced on an ABI 3730 Avant genetic analyser using the BigDye terminator V.3.1 cycle sequencing kit (Applied Biosystems, Thermo Fisher Scientific Oslo, Norway) according to the manufacturer’s instructions.

5.5. Phylogenetic Analysis

Sequences of the 16S rRNA gene of the cyanobacterial strains were analysed using the Seqassem software package (version 07/2008) and the Align MS Windows-based manual sequence alignment editor (version 03/2007) (SequentiX - Digital DNA Processing, Klein Raden, Germany). Segments with highly variable and ambiguous regions and gaps making proper alignment impossible were excluded from the analyses.
A 16S rRNA gene set containing 1135 positions was used in the phylogenetic tree for Cylindrospermopsis/Raphidiopsis. Sphaerospermopsis aphanizomenoides (LN846954) was employed as the outgroup, five Raphidiopsis strains from Meiktila Lake and 35 additional Cylindrospermopsis/Raphidiopsis sequences derived from GenBank were included in the analyses. A set containing 1426 positions was used for the Microcystis 16S rRNA gene analysis. Chroococcus subviolaceus (MF072353) was employed as the outgroup, two strains from Meiktila Lake and 38 additional Microcystis sequences derived from GenBank were included in the analysis. Phylogenetic trees for 16S rRNA genes were constructed using the ML algorithm in Mega v. 7 [55]. In the ML analyses, evolutionary substitution models were evaluated using Mega v. 7 [55]. The HKY+G+I evolutionary model was found to be the best-fitting evolutionary model for the Nostocales 16S rRNA gene tree and T92+G+I for the Microcystis 16S rRNA gene tree. ML analyses of both trees were performed with 1000 bootstrap replicates using Mega v. 7 [55]. The sequence data were submitted to the European Nucleotide Archive (ENA) under the accession numbers listed in Table 2.

5.6. Toxin Analysis

5.6.1. Material

LC–MS/MS utilised standards of ATX-a (Tocris Bioscience, Bristol, UK), homoATX-a (Novakits, Nantes, France) and CYN (Vinci Biochem, Vinci, Italy) and certified reference materials (CRMs) of STX, dcSTX, NeoSTX, GTX1, GTX4, GTX5 and C1 and C2 toxins (National Research Council of Canada, Halifax, NS, Canada (NRC)). LC–HRMS utilised CRMs of MC-RR (3), MC-LR (17) and [Dha7]MC-LR (18) (NRC) and an RM of [D-Leu1]MC-LY [38]. Additional RMs of [D-Asp3]MC-RR (1), D-Asp3]MC-LR (13), MC-YR (14), MC-HilR (21) containing traces of MC-FR (25), MC-WR (26) and MC-LA (41) were prepared at NRC from commercial samples (Enzo Life Sciences, Farmingdale, NY, USA), and extracts containing an array of other identified MCs were available from previous work [12,13,15]. Standards for the Adda-ELISA and for the CYN, ATX and STX ELISAs were as provided with the kits (Abraxis LLC, Warminister, PA, USA).

5.6.2. ELISA for MCs, CYNs, ATXs and STXs

Fresh culture material of two Microcystis and five Raphidiopsis strains was frozen and thawed three times. The R. raciborskii strains were tested for CYNs using the Abraxis Cylindrospermopsin ELISA kit (Abraxis LLC, Warminister, PA, USA) following the manufacturer’s instructions. The test is a direct competitive ELISA that detects cylindrospermopsin but also recognizes deoxycylindrospermopsin and 7-epi-cylindrospermopsin. The ELISA results do not distinguish between dissolved and cell-bound toxins. Both Microcystis strains were tested for microcystins using the Abraxis Microcystins/Nodularins (ADDA) ELISA kits (Abraxis LLC, Warminister, PA, USA). The test is an indirect competitive ELISA designed to detect Adda, (3-amino-9-methoxy-2,6,8-trimethyl-10-phenyldeca-4,6-dienoic acid), based on specific recognition of the Adda moiety [56]. ADDA is a nonprotein amino acid and is the most common side chain at position-5 in microcystins (Figure 5).
All strains were also tested for saxitoxins and anatoxin-a using the Abraxis Saxitoxins (PSP) and Abraxis Anatoxin (VFDF) ELISA kits (Abraxis LLC, Warminister, PA, USA). The saxitoxin ELISA is a direct competitive ELISA that detects saxitoxin based on specific antibody recognition but also recognizes other saxitoxins (e.g., dcSTX, GTXs, lyngbyatoxin, NeoSTX) to varying degrees according to the manufacturer’s instructions. The test for anatoxin-a is a direct competitive ELISA that detects anatoxin-a based on specific antibody recognition but also recognizes homoanatoxin according to the manufacturer’s instructions The colour reaction of all ELISA tests was evaluated at 450 nm on a Perkin Elmer1420 Multilabel counter Victor3 (Perkin Elmer, Waltham, MA, USA), and concentrations were evaluated by manual analysis of the absorbance data as recommended by the vendor.

5.6.3. Microcystin Analysis by LC–HRMS

Fresh culture material of both Microcystis strains was prepared for LC–HRMS by freeze-thawing (3 times), diluting with an equal volume of MeOH and filtering (0.22 μm) [57]). LC–HRMS/MS analysis was performed on a Q Exactive-HF Orbitrap mass spectrometer equipped with a HESI-II heated electrospray ionization interface (ThermoFisher Scientific, Waltham, MA, USA) using an Agilent 1200 LC system including a binary pump, autosampler and column oven (Agilent, Santa Clara, CA, USA). Analyses were performed with SymmetryShield 3.5 µm C18 column (100 × 2.1 mm; Waters, Milford, MA, USA) held at 40 °C with mobile phases A and B of H2O and CH3CN, respectively, each of which contained formic acid (0.1% v/v). Gradient elution (0.3 mL min−1) was from 20% to 90% B over 18 min, then to 100% B over 0.1 min and a hold at 100% B (2.9 min), then returned to 20% B over 0.1 min with a hold at 20% B (3.9 min) to equilibrate the column (total run time 25 min). Injection volume was typically 1–5 µL.
The MS was operated in positive ion mode and calibrated from m/z 74 to 1622. The spray voltage was 3.7 kV, the capillary temperature was 350 °C and the sheath and auxiliary gas flow rates were 25 and 8 units, respectively, with MS data acquired from 2 to 20 min. Mass spectral data were collected using a combined full scan (FS) and data independent acquisition (DIA) method. FS data was collected from m/z 500 to 1400 using the 60,000-resolution setting, an AGC target of 1 × 106 and a max IT of 100 ms. DIA data was collected using the 15,000 resolution setting, an AGC target of 2 × 105, maxIT set to “auto” and a stepped collision energy of 30, 60 and 80 eV. Precursor isolation windows were 62 m/z wide with centering at m/z 530, 590, 650, 710, 770, 830, 890, 950, 1010, 1070, 1130, 1190, 1250, 1310 and 1370. Mass spectral data were also collected using a combined full scan (FS) and top-10 data-dependent acquisition (DDA) method. Data was acquired as for DIA but with an exclusion list generated from a blank injection and an inclusion list (both at ± 5 ppm) from a publicly available database of MC m/z values [58], except that maxIT was set to 100 ms and dynamic exclusion 5.0 s and “if idle pick others” were selected. Putative MCs detected using the above FS/DIA method were further probed in a targeted manner using the parallel reaction monitoring scan (PRM) mode with a 0.7 m/z precursor isolation window, typically using the 30,000-resolution setting, an AGC target of 5 × 105 and a max IT of 400 ms. Typical collision energies were: stepped CE at 30 and 35 eV for MCs with no Arg; stepped CE at 60, 65 and 70 eV for MCs with one Arg; and CE at 65 eV for [M+H]+ and stepped CE at 20, 25 and 30 eV for [M+2H]2+ of MCs with two Arg groups. Full scan chromatograms were obtained in MS-SIM mode as for DIA but with resolution 120,000 and max IT 300 ms.
In negative mode, the mass spectrometer was calibrated from m/z 69 to 1780 and the spray voltage was −3.7 kV, while the capillary temperature, sheath and auxiliary gas flow rates were the same as for positive mode. Mass spectrometry data were collected in FS/DIA scan mode as above using a scan range of m/z 750–1400, a resolution setting of 60,000, AGC target of 1 × 106 and a max IT of 100 ms. For DIA, MS/MS data was collected from m/z 93 to 1400 using a resolution setting of 15,000, AGC target of 2 × 105, max IT set to “auto” and stepped collision energy 65 and 100 eV. Isolation windows were 45 m/z wide and centered at m/z 772, 815, 858, 902, 945, 988, 1032, 1075, 1118, 1162, 1205, 1248, 1294, 1335 and 1378. Mass spectral data were also collected using a combined full scan (FS) and top-10 data dependent acquisition (DDA) method. Data were acquired as for DIA but with an exclusion list generated from a blank injection and an inclusion list from a publicly available list (both at ± 5 ppm) of MC m/z values [58], except that maxIT was set to 100 ms and dynamic exclusion 5.0 s and “if idle pick others” were selected. Full scan chromatograms were obtained over a scan range m/z 750–1400 at a resolution setting of 120,000 using an AGC target of 1 × 106 and a max IT of 300 ms.
Thiol derivatizations were performed by addition of (NH4)2CO3 (0.1 M, 200 µL) to the filtered extract (200 µL), with 200 µL transferred to two LC-MS vials. To one vial was added 1 µL of a 1:1 mixture of mercaptoethanol and d4-mercaptoethanol (Sigma–Aldrich, St. Louis, MO, USA), while 1 µL of water was added to the other vial as a control. Oxidations were performed by addition of DMSO (5 µL) and Oxone (10 mg/mL in water; 25 µL) to 50 µL of extract [15]. Samples and reactions were placed in the sample tray (held at 15 °C) for analysis, and the reactions were monitored periodically until completion and then analysed.

5.6.4. CYN, deoxyCYN, ATXs and STXs Analysis by LC-MS/MS

Extraction was performed on freeze-dried cultures (40 mL), according to the protocol in [59]. In brief, dry material was treated with 6 mL of 50% methanol and sonicated (Omniruptor4000 probe sonicator, Omni-Inc., Kennesaw, MA, USA) for 10 min in pulsed mode (50%) using 160 W power. An aliquot of the solution was then filtered on Phenex RC syringe filters (0.2 μm; Phenomenex, Castel Maggiore, Italy) and analyzed by LC–MS/MS.
LC-MS/MS analysis were performed using a Waters Acquity UPLC system (Waters, Milford, MA, USA) coupled to a SCIEX 4000 QTRAP mass spectrometer (AB Sciex Pte. Ltd., Singapore). Chromatographic separation of analytes was performed using a HILIC column (Ascentis Express OH5, 2.7 μm, 50 × 2.1 mm; Merck Life Science S.r.l., Milan, Italy), while MS detection was performed using positive electrospray ionization using scheduled Multiple Reaction Monitoring. Details of the experimental set up are as described by Cerasino et al. [60], and the method was suitable for the detection and quantification of the following toxins: ATX-a, homoATX-a, CYN, STX, dcSTX, NeoSTX, GTX1, GTX4, GTX5 and C1 and C2 [60]. Quantification limits were 0.2–200 µg L–1. Other toxic alkaloids not available as pure standards were also screened but only for tentative analysis (hydroxy-, epoxy- and homo-ATXs, deoxyCYN, dcNeoSTX, GTX2/3, dcGTX2/3 and C3 and C4 toxins) using equivalent detection settings to their most similar analogs.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/4/232/s1, Figures S1–S59: LC-HRMS/MS spectra of 1, 3, 4, 6, 7, 9, 1115, 17, 1925, and 2756, Figure S60: LC-HRMS spectra of unidentified microcystins 3335, Table S1: Identities of microcystins detected by LC-HRMS/MS analysis in M. aeruginosa strains AB2017/14 and /15, their retention times (tR), concentrations, elemental compositions, observed m/z values in positive and negative ionisation modes, whether they matched retention times in reference samples and whether characteristic microcystin product ions were observed.

Author Contributions

Conceptualization, A.B.; methodology, A.B., C.O.M., L.C., V.H.; validation, A.B., C.O.M., L.C.; investigation, A.B., C.O.M., L.C., M.M., T.S.; data curation, A.B., C.O.M., L.C.; writing—original draft preparation, A.B., T.S., C.O.M.; writing—review and editing, A.B., C.O.M., M.M., L.C., T.S., V.H.; supervision, A.B.; project administration A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Project “Integrated Water Resources Management (IWRM) – Institutional Building and Training in Myanmar (funded by the Norwegian Ministry of Foreign Affairs and the Royal Norwegian Embassy in Myanmar).

Acknowledgments

We thank K. Thomas for preparation of RMs of [D-Asp3]MC-RR (1), D-Asp3]MC-LR (13), MC-YR (14), MC-HilR (21) containing traces of MC-FR (25), MC-WR (26) and MC-LA (41) at National Research Council, Halifax, Canada in partnership with the Norwegian Veterinary Institute (Oslo, Norway).

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Paerl, H.W.; Huisman, J. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environ. Microbiol. Rep. 2009, 1, 27–37. [Google Scholar] [CrossRef] [PubMed]
  2. Sivonen, K. Cyanobacterial toxins. In Encyclopedia of Microbiology; Schaechter, M., Ed.; Elsevier: Oxford, UK, 2009; pp. 290–307. [Google Scholar]
  3. Harvey, G.E. History of Burma from the Earliest Times to 10 March 1824: The Beginning of the English Conquest, 1; Longmans, Green and Company: London, UK, 1925. [Google Scholar]
  4. Nyunt, K.K. Algal flora of Meiktila Lake and its surrounding area. Master´s Thesis, Mandalay University, Mandalay, Myanmar, 1995. [Google Scholar]
  5. Hlaing, N.K. Limnological study of Meiktila Lake, Mandalay Division. Ph.D. Thesis, University of Mandalay, Mandalay, Myanmar, 2014. [Google Scholar]
  6. Klink, T.; de Jongh, I.; Leushuis, M.; Bijlmakers, L.; Speets, R.; Brinkman, L. Scoping Study Meiktila Lake Area IWRM; Royal HaskoningDHV, Arcadis, Rebel Group, Dutch Water Authorities: Den Haag, The Netherlands, 2015; pp. 1–50. [Google Scholar]
  7. Aung, T.T.; Zin, T. Zooplankton community in Meiktila Lake, Meiktila, central Myanmar; Universities Research Journal; The Government of The Republic of the Union of Myanmar Ministry of Education: Yangon, Myanmar, 2015; Volume 7, pp. 1–11. [Google Scholar]
  8. Antunes, J.T.; Leao, P.N.; Vasconcelos, V.M. Cylindrospermopsis raciborskii: Review of the distribution, phylogeography, and ecophysiology of a global invasive species. Front. Microbiol. 2015, 6, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Burford, M.A.; Beardall, J.; Willis, A.; Orr, P.T.; Magalhaes, V.F.; Rangel, L.M.; Azevedo, S.; Neilan, B.A. Understanding the winning strategies used by the bloom-forming cyanobacterium Cylindrospermopsis raciborskii. Harmful Algae 2016, 54, 44–53. [Google Scholar] [CrossRef] [PubMed]
  10. Mowe, M.A.D.; Mitrovic, S.M.; Lim, R.P.; Furey, A.; Yeo, D.C.J. Tropical cyanobacterial blooms: A review of prevalence, problem taxa, toxins and influencing environmental factors. J. Limnol. 2015, 74, 205–224. [Google Scholar] [CrossRef] [Green Version]
  11. Harke, M.J.; Steffen, M.M.; Gobler, C.J.; Otten, T.G.; Wilhelm, S.W.; Wood, S.A.; Paerl, H.W. A review of the global ecology, genomics, and biogeography of the toxic cyanobacterium, Microcystis spp. Harmful Algae 2016, 54, 4–20. [Google Scholar] [CrossRef] [Green Version]
  12. Mallia, V.; Uhlig, S.; Rafuse, C.; Meija, J.; Miles, C.O. Novel microcystins from Planktothrix prolifica NIVA-CYA 544 identified by LC-MS/MS, functional group derivatization and 15N-labeling. Mar. Drugs 2019, 17, 643. [Google Scholar] [CrossRef] [Green Version]
  13. Yilmaz, M.; Foss, A.J.; Miles, C.O.; Özen, M.; Demir, N.; Balci, M.; Beach, D.G. Comprehensive multi-technique approach reveals the high diversity of microcystins in field collections and an associated isolate of Microcystis aeruginosa from a Turkish lake. Toxicon 2019, 167, 87–100. [Google Scholar] [CrossRef]
  14. Miles, C.O.; Sandvik, M.; Haande, S.; Nonga, H.; Ballot, A. LC-MS analysis with thiol derivatization to differentiate [Dhb7]- from [Mdha7]-microcystins: analysis of cyanobacterial blooms, Planktothrix cultures and European crayfish from Lake Steinsfjorden, Norway. Environ. Sci. Technol. 2013, 47, 4080–4087. [Google Scholar] [CrossRef]
  15. Miles, C.O.; Melanson, J.E.; Ballot, A. Sulfide oxidations for LC-MS analysis of methionine-containing microcystins in Dolichospermum flos-aquae NIVA-CYA 656. Environ. Sci. Technol. 2014, 48, 13307–13315. [Google Scholar] [CrossRef]
  16. Lürling, M.; Van Oosterhout, F.; Faassen, E. Eutrophication and warming boost cyanobacterial biomass and microcystins. Toxins 2017, 9, 64. [Google Scholar] [CrossRef] [Green Version]
  17. Guedes, I.A.; Rachid, C.T.C.C.; Rangel, L.M.; Silva, L.H.S.; Bisch, P.M.; Azevedo, S.M.F.O.; Pacheco, A.B.F. Close link between harmful cyanobacterial dominance and associated bacterioplankton in a tropical eutrophic reservoir. Front. Microbiol. 2018, 9, 13. [Google Scholar] [CrossRef] [PubMed]
  18. Humpage, A.R.; Falconer, I.R. Oral toxicity of the cyanobacterial toxin cylindrospermopsin in male Swiss albino mice: determination of no observed adverse effect level for deriving a drinking water guideline value. Environ. Toxicol. 2003, 18, 94–103. [Google Scholar] [CrossRef] [PubMed]
  19. Norris, R.L.; Eaglesham, G.K.; Pierens, G.; Shaw, G.R.; Smith, M.J.; Chiswell, R.K.; Seawright, A.A.; Moore, M.R. Deoxycylindrospermopsin, an analog of cylindrospermopsin from Cylindrospermopsis raciborskii. Environ. Toxicol. 1999, 14, 163–165. [Google Scholar] [CrossRef]
  20. Neumann, C.; Bain, P.; Shaw, G. Studies of the comparative in vitro toxicology of the cyanobacterial metabolite deoxycylindrospermopsin. J. Toxicol. Environ. Health, Part A 2007, 70, 1679–1686. [Google Scholar] [CrossRef]
  21. Machado, J.; Campos, A.; Vasconcelos, V.; Freitas, M. Effects of microcystin-LR and cylindrospermopsin on plant–soil systems: a review of their relevance for agricultural plant quality and public health. Environ. Res. 2017, 153, 191–204. [Google Scholar] [CrossRef] [Green Version]
  22. Griffiths, D.J.; Saker, M.L. The Palm Island mystery disease 20 years on: A review of research on the cyanotoxin cylindrospermopsin. Environ. Toxicol. 2003, 18, 78–93. [Google Scholar] [CrossRef]
  23. Chiswell, R.K.; Shaw, G.R.; Eaglesham, G.; Smith, M.J.; Norris, R.L.; Seawright, A.A.; Moore, M.R. Stability of cylindrospermopsin, the toxin from the cyanobacterium, Cylindrospermopsis raciborskii: Effect of pH, temperature, and sunlight on decomposition. Environ. Toxicol. 1999, 14, 155–161. [Google Scholar] [CrossRef]
  24. Environmental Protection Agency. Cyanobacteria and Cyanotoxins: Information for Drinking Water Systems; United States Environmental Protection Agency: Washington, DC, USA, 2014; pp. 1–11. [Google Scholar]
  25. Li, R.H.; Carmichael, W.W.; Brittain, S.; Eaglesham, G.K.; Shaw, G.R.; Liu, Y.D.; Watanabe, M.M. First report of the cyanotoxins cylindrospermopsin and deoxycylindrospermopsin from Raphidiopsis curvata (cyanobacteria). J. Phycol. 2001, 37, 1121–1126. [Google Scholar] [CrossRef]
  26. Nguyen, T.T.L.; Hoang, T.H.; Nguyen, T.K.; Duong, T.T. The occurrence of toxic cyanobacterium Cylindrospermopsis raciborskii and its toxin cylindrospermopsin in the Huong River, Thua Thien Hue province, Vietnam. Environ. Monit. Assess. 2017, 189, 11. [Google Scholar] [CrossRef]
  27. Dao, T.S.; Cronberg, G.; Nimptsch, J.; Do-Hong, L.-C.; Wiegand, C. Toxic cyanobacteria from Tri An Reservoir, Vietnam. Nowa Hedwigia 2010, 90, 433–448. [Google Scholar] [CrossRef]
  28. Hawkins, P.R.; Chandrasena, N.R.; Jones, G.J.; Humpage, A.R.; Falconer, I.R. Isolation and toxicity of Cylindrospermopsis raciborskii from an ornamental lake. Toxicon 1997, 35, 341–346. [Google Scholar] [CrossRef]
  29. Saker, M.L.; Neilan, B.A. Varied diazotrophies, morphologies, and toxicities of genetically similar isolates of Cylindrospermopsis raciborskii (Nostocales, Cyanophyceae) from northern Australia. Appl. Environ. Microbiol. 2001, 67, 1839–1845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Li, R.; Carmichael, W.W.; Brittain, S.; Eaglesham, G.K.; Shaw, G.R.; Mahakhant, A.; Noparatnaraporn, N.; Yongmanitchai, W.; Kaya, K.; Watanabe, M.M. Isolation and identification of the cyanotoxin cylindrospermopsin and deoxy-cylindrospermopsin from a Thailand strain of Cylindrospermopsis raciborskii (Cyanobacteria). Toxicon 2001, 39, 973–980. [Google Scholar] [CrossRef]
  31. Chonudomkul, D.; Yongmanitchai, W.; Theeragool, G.; Kawachi, M.; Kasai, F.; Kaya, K.; Watanabe, M.M. Morphology, genetic diversity, temperature tolerance and toxicity of Cylindrospermopsis raciborskii (Nostocales, Cyanobacteria) strains from Thailand and Japan. FEMS Microbiol. Ecol. 2004, 48, 345–355. [Google Scholar] [CrossRef] [PubMed]
  32. Zarenezhad, S.; Sano, T.; Watanabe, M.M.; Kawachi, M. Evidence of the existence of a toxic form of Cylindrospermopsis raciborskii (Nostocales, Cyanobacteria) in Japan. Phycol. Res. 2012, 60, 98–104. [Google Scholar] [CrossRef]
  33. Cirés, S.; Wörmer, L.; Ballot, A.; Agha, R.; Wiedner, C.; Velazquez, D.; Casero, M.C.; Quesada, A. Phylogeography of cylindrospermopsin and paralytic shellfish toxin-producing Nostocales cyanobacteria from Mediterranean Europe (Spain). Appl. Environ. Microbiol. 2014, 80, 1359–1370. [Google Scholar] [CrossRef] [Green Version]
  34. Panou, M.; Zervou, S.K.; Kaloudis, T.; Hiskia, A.; Gkelis, S. A Greek Cylindrospermopsis raciborskii strain: missing link in tropic invader’s phylogeography tale. Harmful Algae 2018, 80, 96–106. [Google Scholar] [CrossRef]
  35. Bouaïcha, N.; Miles, C.O.; Beach, D.G.; Labidi, Z.; Djabri, A.; Benayache, N.Y.; Nguyen-Quang, T. Structural diversity, characterization and toxicology of microcystins. Toxins 2019, 11, 714. [Google Scholar] [CrossRef] [Green Version]
  36. Foss, A.J.; Miles, C.O.; Samdal, I.A.; Lovberg, K.E.; Wilkins, A.L.; Rise, F.; Jaabaek, J.A.H.; McGowan, P.C.; Aubel, M.T. Analysis of free and metabolized microcystins in samples following a bird mortality event. Harmful Algae 2018, 80, 117–129. [Google Scholar] [CrossRef]
  37. Okello, W.; Portmann, C.; Erhard, M.; Gademann, K.; Kurmayer, R. Occurrence of microcystin-producing cyanobacteria in Ugandan freshwater habitats. Environ. Toxicol. 2010, 25, 367–380. [Google Scholar] [CrossRef] [Green Version]
  38. LeBlanc, P.; Merkley, N.; Thomas, K.; Lewis, N.I.; Békri, K.; Renaud, S.L.; Pick, F.R.; McCarron, P.; Miles, C.O.; Quilliam, M.A. Isolation and characterization of [D-Leu1]microcystin-LY from Microcystis aeruginosa CPCC-464. Toxins 2020, 12, 77. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Namikoshi, M.; Yuan, M.; Sivonen, K.; Carmichael, W.W.; Rinehart, K.L.; Rouhiainen, L.; Sun, F.; Brittain, S.; Otsuki, A. Seven new microcystins possessing two L-glutamic acid units, isolated from Anabaena sp. strain 186. Chem. Res. Toxicol. 1998, 11, 143–149. [Google Scholar] [CrossRef] [PubMed]
  40. Foss, A.J.; Aubel, M.T. Using the MMPB technique to confirm microcystin concentrations in water measured by ELISA and HPLC (UV, MS, MS/MS). Toxicon 2015, 104, 91–101. [Google Scholar] [CrossRef] [PubMed]
  41. Rinehart, K.L.; Namikoshi, M.; Choi, B.W. Structure and biosynthesis of toxins from blue–green algae (cyanobacteria). J. Appl. Phycol. 1994, 6, 159–176. [Google Scholar] [CrossRef]
  42. Ballot, A.; Sandvik, M.; Rundberget, T.; Botha, C.J.; Miles, C.O. Diversity of cyanobacteria and cyanotoxins in Hartbeespoort Dam, South Africa. Mar. Freshwat. Res. 2014, 65, 175–189. [Google Scholar] [CrossRef] [Green Version]
  43. Fewer, D.P.; Rouhiainen, L.; Jokela, J.; Wahlsten, M.; Laakso, K.; Wang, H.; Sivonen, K. Recurrent adenylation domain replacement in the microcystin synthetase gene cluster. BMC Evol. Biol. 2007, 7, 183. [Google Scholar] [CrossRef] [Green Version]
  44. Scheffer, M.; Carpenter, S.; Foley, J.A.; Folke, C.; Walker, B. Catastrophic shifts in ecosystems. Nature 2001, 413, 591–596. [Google Scholar] [CrossRef]
  45. van Donk, E.; van de Bund, W.J. Impact of submerged macrophytes including charophytes on phyto- and zooplankton communities: allelopathy versus other mechanisms. Aquat. Bot. 2002, 72, 261–274. [Google Scholar] [CrossRef]
  46. Gao, Y.N.; Dong, J.; Fu, Q.Q.; Wang, Y.P.; Chen, C.; Li, J.H.; Li, R.; Zhou, C.J. Allelopathic effects of submerged macrophytes on phytoplankton. Allelopathy J. 2017, 40, 1–22. [Google Scholar]
  47. Mjelde, M.; Faafeng, B.A. Ceratophyllum demersum hampers phytoplankton development in some small Norwegian lakes over a wide range of phosphorus concentrations and geographical latitude. Freshwat. Biol. 1997, 37, 355–365. [Google Scholar] [CrossRef]
  48. Blindow, I.; Hargeby, A.; Andersson, G. Seasonal changes of mechanisms maintaining clear water in a shallow lake with abundant Chara vegetation. Aquat. Bot. 2002, 72, 315–334. [Google Scholar] [CrossRef]
  49. Lombardo, P. Applicability of littoral food-web biomanipulation for lake management purposes: snails, macrophytes, and water transparency in northeast Ohio shallow lakes. Lake Reserv. Manage. 2005, 21, 186–202. [Google Scholar] [CrossRef]
  50. Utermöhl, H. Zur Vervollkommnung der quantitativen Phytoplankton-Methodik. Mitt. Int. Ver. Theor. Angew. Limnol. 1958, 9, 1–38. [Google Scholar] [CrossRef]
  51. Kotai, J. Instructions for Preparation of Modified Nutrient Solution Z8 for Algae.—Publication B-11/69; Norwegian Institute for Water Research: Oslo, Norway, 1972. [Google Scholar]
  52. Komárek, J.; Anagnostidis, K. Cyanoprokaryota 1. Chroococcales.—19/1. In Süsswasserflora von Mitteleuropa; Ettl, H., Gärtner, G., Heynig, H., Mollenhauer, D., Eds.; Gustav Fischer: Jena-Stuttgart-Lübeck-Ulm, Germany, 1999; p. 548. [Google Scholar]
  53. Komárek, J. Cyanoprokaryota. 3. Heterocytous genera. In Süsswasserflora von Mitteleuropa/Freshwater flora of Central Europe; Büdel, B., Krienitz, L., Gärtner, G., Schagerl, M., Eds.; Springer Spektrum: Berlin-Heidelberg, Germany, 2013; p. 1130. [Google Scholar]
  54. Ballot, A.; Cerasino, L.; Hostyeva, V.; Cirés, S. Variability in the sxt gene clusters of PSP toxin producing Aphanizomenon gracile strains from Norway, Spain, Germany and North America. PLoS One 2016, 11, e0167552. [Google Scholar] [CrossRef] [PubMed]
  55. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Fischer, W.J.; Garthwaite, I.; Miles, C.O.; Ross, K.M.; Aggen, J.B.; Chamberlin, A.R.; Towers, N.R.; Dietrich, D.R. Congener-independent immunoassay for microcystins and nodularins. Environ. Sci. Technol. 2001, 35, 4849–4856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Miles, C.O.; Sandvik, M.; Nonga, H.E.; Rundberget, T.; Wilkins, A.L.; Rise, F.; Ballot, A. Thiol derivatization for LC-MS identification of microcystins in complex matrices. Environ. Sci. Technol. 2012, 46, 8937–8944. [Google Scholar] [CrossRef] [Green Version]
  58. Miles, C.O.; Stirling, D. Toxin mass list, version 15b. Available online: https://www.researchgate.net/publication/324039408_Toxinmasslist_COM_v15b (accessed on 11 October 2018).
  59. Cerasino, L.; Meriluoto, J.; Bláha, L.; Carmeli, S.; Kaloudis, T.; Mazur-Marzec, A. SOP 6. Extraction of cyanotoxins from cyanobacterial biomass. In Handbook of Cyanobacterial Monitoring and Cyanotoxin Analysis; Meriluoto, J., Spoof, L., Codd, G.A., Eds.; Wiley: Chichester, UK, 2017; pp. 350–353. [Google Scholar]
  60. Cerasino, L.; Capelli, C.; Salmaso, N. A comparative study of the metabolic profiles of common nuisance cyanobacteria in southern perialpine lakes. Adv. Oceanogr. Limnol. 2017, 8, 22–32. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of Meiktila Lake. The map shows the locations of water sampling (Stations MK1-MK5). The location of Meiktila Lake in Myanmar is shown in the inset.
Figure 1. Map of Meiktila Lake. The map shows the locations of water sampling (Stations MK1-MK5). The location of Meiktila Lake in Myanmar is shown in the inset.
Toxins 12 00232 g001
Figure 2. Micrographs of cyanobacteria investigated in this study. (a) Microcystis novacekii (AB2017/14); (b) Microcystis aeruginosa (AB2017/15); (c) Raphidiopsis raciborskii (AB2017/05); (d) Raphidiopsis raciborskii (AB2017/09); (e) Raphidiopsis raciborskii (AB2017/12); (f) Raphidiopsis raciborskii (AB2017/13); (g) Raphidiopsis raciborskii (AB2017/16). Scale bars indicate 50 µm.
Figure 2. Micrographs of cyanobacteria investigated in this study. (a) Microcystis novacekii (AB2017/14); (b) Microcystis aeruginosa (AB2017/15); (c) Raphidiopsis raciborskii (AB2017/05); (d) Raphidiopsis raciborskii (AB2017/09); (e) Raphidiopsis raciborskii (AB2017/12); (f) Raphidiopsis raciborskii (AB2017/13); (g) Raphidiopsis raciborskii (AB2017/16). Scale bars indicate 50 µm.
Toxins 12 00232 g002
Figure 3. ML tree based on partial 16S rRNA gene sequences of 40 Raphidiopsis/Cylindrospermospis strains. Outgroup = Sphaerospermopsis aphanizomenoides (LN846954). Cluster I includes Cylindrospermopsis and Raphidiopsis strains from Asia, Europe, Africa, Australia and North America, cluster II includes strains from North and South America (USA, Mexico, Brazil), North Africa (Tunisia), Southwest Europe (Spain) and New Zealand. Strains from this study are marked in bold. Bootstrap values above 50 are included. The scale bar indicates 0.5% sequence divergence.
Figure 3. ML tree based on partial 16S rRNA gene sequences of 40 Raphidiopsis/Cylindrospermospis strains. Outgroup = Sphaerospermopsis aphanizomenoides (LN846954). Cluster I includes Cylindrospermopsis and Raphidiopsis strains from Asia, Europe, Africa, Australia and North America, cluster II includes strains from North and South America (USA, Mexico, Brazil), North Africa (Tunisia), Southwest Europe (Spain) and New Zealand. Strains from this study are marked in bold. Bootstrap values above 50 are included. The scale bar indicates 0.5% sequence divergence.
Toxins 12 00232 g003
Figure 4. ML tree based on partial 16S rRNA gene sequences of 40 Microcystis strains. Outgroup = Chroococcus subviolaceus (MF072353). Strains from this study are marked in bold. Bootstrap values above 50 are included. The scale bar indicates 2% sequence divergence.
Figure 4. ML tree based on partial 16S rRNA gene sequences of 40 Microcystis strains. Outgroup = Chroococcus subviolaceus (MF072353). Strains from this study are marked in bold. Bootstrap values above 50 are included. The scale bar indicates 2% sequence divergence.
Toxins 12 00232 g004
Figure 5. Structures and exact negative ionisation m/z of microcystins (MCs) identified in cultures AB2017/14 and AB2017/15 in this study, showing characteristic product ions at m/z 135.0804 (positive) and 128.0353 (negative) LC–MS/MS spectra (see Table 4 and Table S1). The origins of additional positive mode product ions containing Mdha7 (R3 = Me) at m/z 375.1914 and 446.2286 are also shown. Note that the corresponding product ions containing Dha7 (R3 = H) have m/z 361.1758 and 432.2129, and m/z 393.2020 and 464.2391 for Mser7 (R4 = OH). A full version of this table including positive ionisation data is shown in the Supporting Information (Table S1). Abbreviations: Abu, aminobutyric acid; Dha, dehydroalanine; (H2)Tyr, dihydrotyrosine; (H4)Tyr, 4,5,6,7-tetrahydrotyrosine; Kyn, kynurenine; Mdha, N-methyldehydroalanine; Mlan, N-methyllanthionine; Mser, N-methylserine; Oia, oxindolyalanine; Tyr(OMe), methoxytyrosine.
Figure 5. Structures and exact negative ionisation m/z of microcystins (MCs) identified in cultures AB2017/14 and AB2017/15 in this study, showing characteristic product ions at m/z 135.0804 (positive) and 128.0353 (negative) LC–MS/MS spectra (see Table 4 and Table S1). The origins of additional positive mode product ions containing Mdha7 (R3 = Me) at m/z 375.1914 and 446.2286 are also shown. Note that the corresponding product ions containing Dha7 (R3 = H) have m/z 361.1758 and 432.2129, and m/z 393.2020 and 464.2391 for Mser7 (R4 = OH). A full version of this table including positive ionisation data is shown in the Supporting Information (Table S1). Abbreviations: Abu, aminobutyric acid; Dha, dehydroalanine; (H2)Tyr, dihydrotyrosine; (H4)Tyr, 4,5,6,7-tetrahydrotyrosine; Kyn, kynurenine; Mdha, N-methyldehydroalanine; Mlan, N-methyllanthionine; Mser, N-methylserine; Oia, oxindolyalanine; Tyr(OMe), methoxytyrosine.
Toxins 12 00232 g005
Figure 6. LC–HRMS full scan extracted ion chromatograms (3.8–19.4 min, positive ionization mode) of extracts from (A) Microcystis culture AB2017/14 and (B) Microcystis culture AB2017/15. Chromatograms were produced by extracting at m/z (± 5 ppm) for all MCs listed in Table 4 (see Table S1 for positive ionisation m/z values). Note that some of the smaller peaks are not labelled on the chromatograms, and the peak marked with an asterisk is not from a MC.
Figure 6. LC–HRMS full scan extracted ion chromatograms (3.8–19.4 min, positive ionization mode) of extracts from (A) Microcystis culture AB2017/14 and (B) Microcystis culture AB2017/15. Chromatograms were produced by extracting at m/z (± 5 ppm) for all MCs listed in Table 4 (see Table S1 for positive ionisation m/z values). Note that some of the smaller peaks are not labelled on the chromatograms, and the peak marked with an asterisk is not from a MC.
Toxins 12 00232 g006
Table 1. Biomass (mg L−1 FW) of phytoplankton groups at sampling points MK1–MK5 in Meiktila Lake in March and November of 2017.
Table 1. Biomass (mg L−1 FW) of phytoplankton groups at sampling points MK1–MK5 in Meiktila Lake in March and November of 2017.
Sampling PointMK1MK 1MK 2MK 2MK 3MK 3MK 4MK 4MK 5MK 5
Sampling DateMarNovMarNovMarNovMarNovMarNov
Phytoplankton Group
Bacillariophyceae0.5250.1491.3160.1010.0350.3010.0830.0640.0880.059
Chlorophyceae0.1430.0550.1820.1460.1020.0450.1010.1020.1070.165
Chrysophyceae0.01100.00200.0450.0210.0060.00700.004
Conjugatophyceae 0.0990.0060.14100.01300.0080.0170.0450
Cryptophyceae0.0230.1320.0070.1700.1550.2990.1270.2750.1060.138
Cyanobacteria1.0631.9651.6401.0500.6060.0781.0641.9900.7512.397
Dinophyceae0.1860.0050.22200.28500.0840.0060.0780
Euglenophyceae0.2550.0290.1590.0050.0140.0500.0090.0020.0530.021
Eustigmatophyceae00000000.00200
Klebsormidiophyceae000000.023000.0020
Prymnesiophyceae00.0020.0000.0020.0030.0040.0020.0030.0080.007
Trebouxiophyceae000.007000.0050.005000.001
Xanthophyceae000.00800.00500.0030.00100
Total2.3062.3443.6841.4741.2640.8271.4932.4671.2382.790
Table 2. Strains isolated from Meiktila Lake, strain codes and European Nucleotide Archive (ENA) accession numbers.
Table 2. Strains isolated from Meiktila Lake, strain codes and European Nucleotide Archive (ENA) accession numbers.
SpeciesStrainAccession nr.
16S rRNA Gene
Raphidiopsis
R. raciborskiiAB2017/05LR590626
R. raciborskiiAB2017/09LR590627
R. raciborskiiAB2017/12LR590628
R. raciborskiiAB2017/13LR590629
R. raciborskiiAB2017/16LR746263
Microcystis
MicrocystisAB2017/14LR590630
MicrocystisAB2017/15LR590631
Table 3. Concentrations (µg mg-1 FW) of CYNs by ELISA and of CYN and deoxyCYN by LC–MS/MS in cultured R. raciborskii strains isolated from Meiktila Lake*.
Table 3. Concentrations (µg mg-1 FW) of CYNs by ELISA and of CYN and deoxyCYN by LC–MS/MS in cultured R. raciborskii strains isolated from Meiktila Lake*.
StrainELISALC-MS/MS
CYNs CYNdeoxyCYNCYN (%)deoxyCYN (%)
AB2017/092.18--3862
AB2017/05n.d.n.d.n.d.n.d.n.d.
AB2017/161.841.651.255743
AB2017/12n.d.n.d.n.d.n.d.n.d.
AB2017/134.312.467.292575
*- = biomass not determined; n.d. = not detected; FW = fresh weight; percentages are of total CYNs by LC–MS/MS.
Table 4. Identities of microcystins detected by LC-HRMS/MS analysis in Microcystis strains AB2017/14 and /15 isolated from Meiktila Lake, their retention times (tR), concentrations, relative abundances (%) and observed m/z values in negative ionisation modea.
Table 4. Identities of microcystins detected by LC-HRMS/MS analysis in Microcystis strains AB2017/14 and /15 isolated from Meiktila Lake, their retention times (tR), concentrations, relative abundances (%) and observed m/z values in negative ionisation modea.
m/zCompound NameConfidencetR (min)Concentrationb
AB2017/14AB2017/15
μg g−1%μg g−1%
11022.5443[D-Asp3]MC-RRcConfirmed4.182.40.211002.47.21
21022.5454[Dha7]MC-RRcProbable4.530.10.0131.80.23
31036.5597MC-RRcConfirmed4.5519.41.741458.010.49
4995.4858[D-Asp3]MC-ERcProbable6.14ND-119.60.86
51114.5657MC-LR–CysdProbable6.160.10.01ND-
61045.5378MC-(H2)YRcProbable6.1832.42.90ND-
71031.5224[D-Asp3]MC-(H2)YRcProbable6.2032.92.95ND-
81033.5381[D-Asp3]MC-(H4)YRcProbable6.270.60.05148.11.07
91009.5015MC-ERcProbable6.550.60.0537.80.27
101029.5072[D-Asp3]MC-YRcProbable6.59ND-144.11.04
111047.5540MC-(H4)YRcProbable6.646.40.57310.82.24
121011.5530[Mser7]MC-LRProbable6.800.40.0413.30.10
13979.5273[D-Asp3]MC-LRcConfirmed6.8940.83.651870.313.45
141043.5224MC-YRcConfirmed6.999.10.81317.42.28
151011.4999MC-MRc,d,eProbable6.993.60.32ND-
161073.5329MC-Y(OMe)RcProbable7.077.20.64ND-
17993.5435MC-LRcConfirmed7.12183.016.383332.523.97
18979.5289[Dha7]MC-LRcConfirmed7.130.80.0737.40.27
191080.5170Unidentified MCcUnidentified7.1559.05.28ND-
201070.5332MC-KynRcProbable7.354.50.40ND-
211007.5582MC-HilRcConfirmed7.3841.83.7417.90.13
221080.5169UnidentifiedcUnidentified7.3859.65.34ND-
231082.5311MC-OiaRcTentative7.40102.49.17ND-
24995.4857[D-Asp3]MC-REcProbable7.43ND-110.10.79
251027.5269MC-FRcConfirmed7.48255.222.85ND-
261066.5380MC-WRcConfirmed7.640.70.06ND-
271009.5015MC-REcProbable7.793.80.341281.09.21
28951.4955MC-RAcProbable8.24131.16ND-
29965.5116MC-RAbucProbable8.837.10.64ND-
30946.4577MC-(H2)YGcProbable9.065.40.48ND-
311011.4966MC-RMc,d,eProbable9.230.20.02ND-
32960.4738MC-(H2)YAcProbable9.3410.10.90ND-
33953.4783Unidentified MCcUnidentified9.50121.07ND-
34967.4942Unidentified MCcUnidentified9.782.30.21ND-
35984.4735Unidentified MCcUnidentified9.963.70.33ND-
36968.4272[D-Asp3]MC-EEcProbable11.43ND-403.32.90
37982.4432MC-EEcProbable12.580.90.081533.511.03
38952.4690[D-Asp3]MC-LEcProbable12.821.20.11358.82.58
39894.4631[D-Asp3]MC-LAcProbable13.326.50.58ND-
40966.4849MC-LEcProbable13.6912.91.151376.19.90
41908.4788MC-LAcConfirmed14.8744.84.01ND-
42922.4942MC-HilAcProbable15.295.90.53ND-
43981.4737MC-WAcProbable15.884.40.39ND-
44942.4629MC-FAcProbable15.9542.33.79ND-
45922.4943MC-LAbucProbable16.0916.61.49ND-
46936.5110MC-HilAbucProbable16.552.20.20ND-
47995.4898MC-WAbucProbable17.2013.41.20ND-
48956.4784MC-FAbucProbable17.3926.42.36ND-
49936.5106MC-LVcProbable17.623.60.32ND-
50950.5254MC-LLcProbable17.933.90.35ND-
51970.4944MC-FVcProbable18.400.70.06ND-
52950.5262iso-MC-LLcTentative18.491.20.11ND-
531009.5041MC-WVcProbable18.542.90.26ND-
541023.5210MC-WLcProbable18.742.80.25ND-
55970.4942iso-MC-FVcTentative18.792.10.19ND-
561023.5216iso-MC-WLcTentative19.161.30.12ND-
a A comprehensive version of this table, including positive and negative ionisation MS data, reactivity towards thiols and mild oxidising agents, number of rings plus double-bond equivalents (RDBE) and presence of characteristic ions observed in positive and negative ionisation MS/MS spectra, is in the Supporting Information (Table S1) together with LC–HRMS/MS spectra (Figures S1–S59). b Concentration expressed per weight of biomass (FW) and as a percentage of total microcystins detected in each culture); ND = not detected; c Reacted with mercaptoethanol; d Oxidised by Oxone/DMSO; e HRMS/MS spectrum of oxidation product obtained.
Table 5. Sampling points and sampling depth in Meiktila Lake for chemical and biological measurements.
Table 5. Sampling points and sampling depth in Meiktila Lake for chemical and biological measurements.
Sampling Point Water Depth (m)Depth of Integrated Sample (m)Geographical Position
MK1 3.30–1N 20° 52’ 59.196, E 95° 51’ 12.204
MK220–1N 20° 53’ 21.48, E 95° 51’ 2.124
MK32.50–1N 20° 54’ 16.38, E 95° 50’ 40.092
MK44.40–2N 20° 52’ 21.468, E 95° 51’ 12.528
MK57.10–3N 20° 51’ 58.752, E 95° 51’ 18.936

Share and Cite

MDPI and ACS Style

Ballot, A.; Swe, T.; Mjelde, M.; Cerasino, L.; Hostyeva, V.; Miles, C.O. Cylindrospermopsin- and Deoxycylindrospermopsin-Producing Raphidiopsis raciborskii and Microcystin-Producing Microcystis spp. in Meiktila Lake, Myanmar. Toxins 2020, 12, 232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12040232

AMA Style

Ballot A, Swe T, Mjelde M, Cerasino L, Hostyeva V, Miles CO. Cylindrospermopsin- and Deoxycylindrospermopsin-Producing Raphidiopsis raciborskii and Microcystin-Producing Microcystis spp. in Meiktila Lake, Myanmar. Toxins. 2020; 12(4):232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12040232

Chicago/Turabian Style

Ballot, Andreas, Thida Swe, Marit Mjelde, Leonardo Cerasino, Vladyslava Hostyeva, and Christopher O. Miles. 2020. "Cylindrospermopsin- and Deoxycylindrospermopsin-Producing Raphidiopsis raciborskii and Microcystin-Producing Microcystis spp. in Meiktila Lake, Myanmar" Toxins 12, no. 4: 232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12040232

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop