Next Article in Journal
Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell
Next Article in Special Issue
Toxic or Otherwise Harmful Algae and the Built Environment
Previous Article in Journal
Botulinum Toxin Injection for Medically Refractory Neurogenic Bladder in Children: A Systematic Review
Previous Article in Special Issue
Algal Toxic Compounds and Their Aeroterrestrial, Airborne and other Extremophilic Producers with Attention to Soil and Plant Contamination: A Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 13
Issue 7
10.3390/toxins13070448
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Report on Microcystis as a Potential Microviridin Producer in Bulgarian Waterbodies

by
Blagoy Uzunov
1,*,
Katerina Stefanova
2,
Mariana Radkova
2,
Jean-Pierre Descy
3,
Georg Gärtner
4 and
Maya Stoyneva-Gärtner
1,*
1
Department of Botany, Faculty of Biology, Sofia University, 8 blvd. Dragan Zankov, 1164 Sofia, Bulgaria
2
AgroBioInstitute, Bulgarian Agricultural Academy, 8 blvd. Dragan Zankov, 1164 Sofia, Bulgaria
3
Unité d’Océanographie Chimique, Université de Liège, Sart Tilman, 4000 Liège, Belgium
4
Institut für Botanik der Universität Innsbruck, Sternwartestrasse 15, 6020 Innsbruck, Austria
*
Authors to whom correspondence should be addressed.
Toxins 2021, 13(7), 448; https://doi.org/10.3390/toxins13070448
Submission received: 30 May 2021 / Revised: 24 June 2021 / Accepted: 25 June 2021 / Published: 28 June 2021
(This article belongs to the Special Issue Ecology of Aerophytic and Freshwater Cyanobacteria and Algae with Special Reference to Toxic Species)
Browse Figures
Versions Notes

Abstract

:
Bulgaria, situated on the Balkan Peninsula, is rich in small and shallow, natural and man-made non-lotic waterbodies, which are threatened by blooms of Cyanoprokaryota/Cyanobacteria. Although cyanotoxins in Bulgarian surface waters are receiving increased attention, there is no information on microviridins and their producers. This paper presents results from a phytoplankton study, conducted in August 2019 in three lakes (Durankulak, Vaya, Uzungeren) and five reservoirs (Duvanli, Mandra, Poroy, Sinyata Reka, Zhrebchevo) in which a molecular-genetic analysis (PCR based on the precursor mdnA gene and subsequent translation to amino acid alignments), combined with conventional light microscopy and an HPLC analysis of marker pigments, were applied for the identification of potential microviridin producers. The results provide evidence that ten strains of the genus Microcystis, and of its most widespread species M. aeruginosa in particular, are potentially toxigenic in respect to microviridins. The mdnA sequences were obtained from all studied waterbodies and their translation to amino-acid alignments revealed the presence of five microviridin variants (types B/C, Izancya, CBJ55500.1 (Microcystis 199), and MC19, as well as a variant, which was very close to type A). This study adds to the general understanding of the microviridin occurrence, producers, and sequence diversity.
Keywords:
coastal lake; Cyanobacteria; Cyanoprokaryota; cyanotoxins; harmful algal blooms; Microcystis aeruginosa; reservoir
Key Contribution: The paper presents the first documentation of potential microviridin producers in Bulgarian non-lotic waterbodies from the genus Microcystis, their genetic diversity, and the possibility for release of five microviridin variants.

1. Introduction

Currently, the problems caused by harmful algal blooms (HABs) of Cyanoprokaryota/Cyanobacteria (hereafter abbreviated as CyanoHABs) and their toxins (cyanotoxins) are recognized worldwide [1,2,3]. However, since the beginning of the 21st century, a series of studies have demonstrated that the assessment of cyanobacterial toxicity cannot solely rely on the commonly known cyanotoxins, implying the activity of other compounds, the structure and function of which need elucidation (for details see [4]). In this regard, beyond the best studied and routinely monitored cyanotoxins from the specific group of cyclic non-ribosomally produced heptapeptides, named microcystins [1,2,3], more than six hundred other peptides or peptidic metabolites, commonly unified as cyanopeptins (CNPs), were isolated from the vast metabolite repertory of cyanoprokaryotes [5,6,7,8]. Despite some variations in composition and amounts, these CNPs were regularly found during CyanoHABs [7]. Nowadays, CNPs are gaining strong attention because it has been reported that metabolomic profiles consisting of different CNPs may affect the invertebrates and fish populations [7]. The CNPs include the unique 16-membered family of the toxins microviridins (MVs), named after their first discovered producer—Microcystis viridis (A. Braun) Lemmermann (strain NIES-102) [9]. MVs are peculiar ribosomally synthesized and post-translationally modified peptides (depsipeptides), most of which inhibit serine proteases and have been found in different species of aquatic Microcystis Lemmermann, Anabaena Bory ex Bornet et Flahault s.l., Planktothrix Anagnostidis et Komárek (Syn. Oscillatoria Vaucher ex Gomont p.p.), Nodularia Mertens ex Bornet et Flahault and Radiocystis Skuja, and in some aeroterrestrial/freshwater strains of Nostoc Vaucher ex Bornet et Flahault (e.g., [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38]). Some features given to MVs were suggested to be related to their allelochemical properties and to the affecting of proteolysis, in particular [8,39]. However, further research in this field is needed to clarify the function and potential ecological role of MVs (e.g., [8,35]). The biosynthesis of MVs was elucidated and it was demonstrated that the microviridin (MV) biosynthetic clusters have different organizations [8,18]. Although the subsequent application of the PCR approach revealed a far more widespread occurrence of MV genes than expected, and their global distribution was suggested, the natural diversity of MV precursor sequences remained almost unknown and MVs have been detected infrequently in cyanoprokaryote peptide screenings [8,18].
Bulgaria, situated on the Balkan Peninsula, has numerous (ca. 10,000) wetlands; however, they occupy less than 1% of its territory and, being mostly shallow, small, and located in lowlands or plains, are vulnerable to threats of CyanoHABs [40]. Although the first records on algal flora dated from the end of the 18th century and targeted investigations on cyanotoxins started after 2000 [41], there are no studies on MVs or their producers in the country. This paper provides first data on the presence and diversity of potential MV producers, and on their contribution to the phytoplankton in eight selected Bulgarian waterbodies, obtained after combined studies by conventional light microscopy (LM), HPLC pigment marker analysis, and molecular-genetic studies. The latter include PCR amplification of the entire MV precursor gene A (mdnA) [8,18], the translation of obtained nucleotide sequences to amino acids, and the determination of their alignment with the intention to identify the MV types. The results provide evidence for ten different Microcystis strains as potential producers of five MV variants.

2. Results

2.1. Phytoplankton Species Composition and Abundance, Obtained by Light Microscopy (LM)

In total, 171 species were identified using LM in the phytoplankton of the eight studied waterbodies. They belonged to seven algal phyla, namely Cyanoprokaryota, Chlorophyta, Streptophyta, Ochrophyta, Cryptophyta, Pyrrhophyta, and Euglenophyta (Figure 1). Cyanoprokaryotes, represented by 49 species, had an important contribution to the total phytoplankton biodiversity (29%), occupying second position after Chlorophyta (49%).
The total number of phytoplankton species ranged between six (reservoir Sinyata Reka) and 65 (reservoir Duvanli), and the number of cyanoprokaryotes was between five (reservoir Poroy) and 25 (reservoir Mandra, site East). They contributed from 14% (lake Uzungeren) to 83% (reservoir Sinyata Reka) of the total biodiversity at each site (Figure 2).
Cyanoprokaryotes comprised from 11% (Lake Durankulak, Eastern site) to 99% (reservoir Sinyata Reka) of the total phytoplankton biomass (Figure 3). On average they contributed 62% to the total phytoplankton biomass (Figure 3). In most sites where the cyanoprokaryote biomass exceeded this average value, the phytoplankton was dominated by filamentous species: Aphanizomenon cf. klebahnii Elenkin ex Pechar (dominant in all sites of the coastal reservoirs Mandra and Poroy), Sphaerospermopsis torques-reginae (Komárek) Werner, Laughinghouse IV, Fior et Sant’Anna (dominant in the inland reservoir Sinyata Reka), Planktothrix isothrix (Skuja) Komárek et Komárková and P. suspensa (Pringsheim) Anagnostidis & Komárek (co-dominants in the coastal Lake Vaya), Pseudanabaena limnetica (Lemmermann) Komárek and Spirulina cf. laxissima G. S. West (co-dominants in the inland reservoir Duvanli).
Considering the results on potential MV-producers, obtained by the PCR analysis, below we present in more detail the LM data on the genus Microcystis, from which four species (morphospecies) were identified: M. aeruginosa (Kützing) Kützing, M. cf. comperei Komárek, M. natans Lemmermann ex Skuja, and M. wesenbergii (Komárek) Komárek (Table 1). In some samples we also found separate cells, disintegrated colonies, or initial colonies, as well as colonies with transitional morphology, for which species identification was not possible (Table 1). Microcystis was observed in almost all waterbodies (except Vaya and Uzungeren) with a low number of certainly identified species, ranging in separate sites between one (in the reservoir Poroy) and three (in the reservoir Duvanli) (Figure 2, Table 1). Its contribution to the total phytoplankton biomass was different, but always extremely low—between <0.05% and <1% (Table 1). The only exception was the higher amount of M. aeruginosa (<5%) in the eastern part of the reservoir Mandra (Table 1).

2.2. Results on General Phytoplankton Composition from HPLC Analysis of Marker Pigments

According to HPLC determination of marker pigments concentrations, processed with CHEMTAX for estimating phytoplankton class abundance [42,43,44,45], cyanoprokaryotes dominated the phytoplankton of most waterbodies (Figure 4).
The values of chlorophyll a indicated the meso- to hypertrophic status of the studied waterbodies. The chlorophyll a concentration ranged between 6 (coastal lake Durankulak, site East) to 316 µg L−1 (reservoir Sinyata Reka) and 83 µg L−1 (coastal lake Vaya). The contribution of cyanoprokaryotes to the phytoplankton biomass ranged in a similar way, being the lowest in the eastern part of Durankulak (8%) and the highest in Vaya (97%) and Sinyata Reka (81%) (Figure 4).

2.3. Results from PCR Analysis for Microcystin-Producing Strains

The precursor gene mdnA from the MV synthetized gene cluster [18] was successfully amplified in all ten investigated metagenomic DNA samples from the studied waterbodies and 22 sequences were obtained. Checking these in the National Centre for Biotechnology Information (NCBI) genetic database [46] by the standard Basic Local Alignment Search Tool (BLAST) [47], revealed ten different sequences, which represent ten strains. They showed high homology (97–100%) with different mdnA-containing Microcystis strains, published in NCBI [46]. All obtained sequences with their corresponding highly homologous NCBI strains are divided in two clusters in the constructed phylogenetic tree (Figure 5).
The first cluster contains two clearly defined subclusters. Subcluster I is formed mainly by the sequences, isolated from the reservoirs Poroy and Sinyata Reka, which show high homology (98.7–99.1%) with two strains of M. aeruginosa—NIES-298 and NIES-2481. Subcluster II contains most of the sequences (eight), isolated from the reservoir Mandra, which show 97.3–100% homology with two Microcystis strains unidentified to species level (Microcystis sp. FN 668693.1 and Microcystis sp. 199), the first of which was uncultured.
The second cluster contains sequences, which show high homology (99.1–99.6%) with the strain Microcystis sp. MC19. They are isolated from five different waterbodies, including two sequences from both sides of the reservoir Mandra, and single sequences from the reservoirs Duvanli and Zhrebchevo, and from the coastal lakes Durankulak and Uzungeren as well.
Interestingly, the isolate from the coastal Lake Vaya differs from all other obtained sequences and the two single sequences obtained from both sites of Lake Durankulak are quite different: the sequence from the western part of the lake is situated in cluster I, whereas the sequence from the eastern part of the same lake is positioned in cluster II (Figure 5).
We intended to identify the MV variants through the translation of the obtained mdnA sequences and their comparison with published leader peptide sequences and amino acid alignments [8,18], processed by the Vector NTI Advance 11.5 software package (Version 11.5, Invitrogen Corporation, Carisbad, CA 92008, USA, 2010) and with part of the consensus sequence designed by WebLogo [48].
The determined MV sequences (Figure 6) show that most of our strains should be capable of producing MV with the highly conserved PFFARFL motif from the α-helix of the MV leader peptide structure [8] (Figure 6). The exceptions are in: (1) the partially sequenced strain Blu from the small inland reservoir Sinyata Reka, which contained less amino acids in comparison with all other strains and in which only GRFL was detected in the α-helix; (2) the strain Man(W) 5, isolated from the Mandra reservoir, in which PLFARFL was found in the α-helix (Figure 6). The KYPSD sequence from the standard MV consensus sequence TXK(Y/W/F)PS DW(E/G)(E/D), firstly named as the MV core motif [18], and later represented mostly as the larger canonical TxKxPSD motif [8,25], was found in all mdnA sequences isolated in this study (Figure 6). The order of the neighboring amino acids before and after this core motif differed in the obtained strains (Figure 6) and was compared with the standard published parts of MV sequences [8,18] (Table 2).
The comparison of data demonstrated that:
(1) The peculiar MV sequence with the shorter α-helix, obtained from the strain Blu, isolated from the small inland reservoir Sinyata Reka, was close but not identical with the MV sequence from M. aeruginosa NIES-298 (Figure 6). According to the alignment of amino acids neighboring to KYPSD, both sequences (Blu and NIES298) coincided completely with MV variant B/C (Table 2);
(2) The MV sequence from the single isolate Por 2 from the small coastal reservoir Poroy, was very close but not identical to the MV variants of M. viridis NIES 102 and M. aeruginosa 843 (Figure 6). According to the sequences neighboring to KYPSD, it was very close to MV variant A (Table 2);
(3) The MV sequence of the strain Por 1 (and of its identical strains Por 3,4,5) from the same Poroy reservoir completely coincided with the MV variant, isolated from Microcystis sp. strain Izancya (FN667620.1) and named MV Izancya [18] (Figure 6; Table 2). This MV has not been referred to as one of the main MV variants, published by do Amaral et al. [8];
(4) Most MV sequences obtained from Mandra, as well as the sequences from Vaya and Durankulak (western part) are identical with the alignment of Microcystis sp. FN668693 (except the above-mentioned Man (W) 5). They all are very close to the MV of Microcystis sp. 199, which, due to a partially known sequence [46] has a shorter published α-helix (Figure 6). The MV of Microcystis sp. 199, named MV 199 [18], later has been referred to as a separate MV variant, labelled as CBJ55500.1 (Microcystis 199) [8]. According to the alignment of amino acids neighboring to KYPSD, all MV sequences of the strains Dur (W), Man (E) 1,2,4, Man (W) 5, and Vai coincide completely with the MV variant CBJ55500.1 (Microcystis 199) (Table 2).
(5) The five MV sequences, obtained from the strains isolated from Durankulak (eastern part), Duvanli, Mandra (eastern part), Uzungeren and Zhrebchevo, are similar to the MV variant of the strain Microcystis sp. MC19, despite the last two having a shorter α-helix (Figure 6). According to the alignment of amino acids neighboring to KYPSD, all MV sequences isolated from the strains Man (E) 3, Dur (E), Duv, Uz, and Zh are completely similar with one of the MV variants, identified from field samples by Ziemert et al. [18] (Table 2).

3. Discussion

The results obtained during this study provide the first evidence for the presence of potential MV producers in Bulgarian waterbodies. According to the distribution of the precursor gene mdnA, all ten new potential MV-producing strains belonged to the genus Microcystis and to its species M. aeruginosa in particular (Figure 5). This morphospecies was the most frequent in the studied waterbodies, especially when considering that some of the separate cells and disintegrated or transitional colonies may also belong to it (Figure 1, Table 1). This consideration can explain the PCR result for the mdnA sequence in all waterbodies, including those for which certain Microcystis morphospecies have not been indicated (Figure 5). Our finding of M. aeruginosa as the most common potential MV producer is in complete agreement with the conclusion in the summary of Amaral et al. [8] (p.14), who stated that “the genus Microcystis and the species M. aeruginosa are the largest producers of microviridins—currently, of the 25 isolated microviridins, 11 belong to the genus Microcystis, and eight of these belong to the species M. aeruginosa”. At present, it is impossible to refer to the three other Microcystis species detected by LM (Microcystis cf. comperei, Microcystis natans, and M. wesenbergii) as potential MV producers. However, it should be noted that some of the obtained sequences belong to three unidentified Microcystis strains, available in NCBI [46] (Figure 5). Interestingly, the single strain Por 2, isolated from the small coastal reservoir Poroy and its MV sequence are close, but not identical to the first known MV producer and its MV—M. viridis NIES-102. By LM M. viridis was not identified during this study nor in our samples from the same coastal wetlands Durankulak, Poroy, Uzungeren, Mandra and Vaya, as well as in the inland reservoir Sinyata Reka, collecte d in the previous year, 2018 [49]. This species was, however, present but rare in the samples from the coastal Lake Vaya collected in August 2004 [50,51].
We did not obtain any PCR signal in relation with the mdnA gene and other cyanoprokaryote species, found as dominants during the blooms (Figure 3), namely Aphanizomenon cf. klebahnii, Sphaerospermopsis torques-reginae, Planktothrix isothrix and P. suspensa, Pseudanabaena limnetica, and Spirulina cf. laxissima. This result agrees with the fact that these algae have not been pointed out as potential MV producers [8,9,11,12,13,15,16,19,23,26,27,28]. However, Ziemert et al. [18] (p. 3570) stressed that “The fact that no PCR product was obtained thus does not necessarily imply that the strains lack the capacity to produce microviridin-like peptides. Rather, it is possible that the orthologous genes in these strains are in a different order or in an independent position in the genome”. Therefore, future genetic and biochemical studies are needed to identify all possible MV producers.
The translation, applied to the mdnA sequences, revealed four variants of MVs, which could be related with certainty to known MV types: B/C, MV from M. isancya, CBJ55500.1 (Microcystis 199), and MV from Microcystis MC19. One more MV variant was very close, but not identical to MV type A and to MVs, isolated from M. viridis NIES 102 and M. aeruginosa 843 [8,18] (Figure 6, Table 2).
All results on the geographic diversity of Microcystis strains, obtained both by LM and molecular-genetic methods, as well as on the potential release of five MV variants, occurred in a context of high cyanoprokaryote contribution to phytoplankton biodiversity and biomass, as shown by LM and HPLC analyses (Figure 2, Figure 3 and Figure 4). All observations point to the meso- to hypertrophic character of the eight studied waterbodies and once more stress the strong vulnerability to CyanoHABs of the shallow lowland and plain wetlands of Bulgaria [40,41], in which four general types of cyanotoxins (microcystins, anatoxins, saxitoxins, and cylindrospermopsins) have been found (e.g., [41,52]). Although Microcystis was found in small amounts, likely due to strong rains preceding the sampling [53], this genus was commonly reported among the most widespread causative agents of CyanoHABs in Bulgaria [41,51] and its mass development in dry summer periods is to be expected. Therefore, the finding of its MV-producing strains in the country stresses the need for future larger scale studies of these peculiar cyanotoxins.

4. Materials and Methods

4.1. Sites and Sampling

The sampling was carried out in the period 14–21 August 2019 in eight waterbodies situated in Central and Eastern Bulgaria (Figure 7, Table 3). Detailed data on the morphometry, historical development, usage and conservation value of the studied three coastal lakes (Durankulak, Vaya and Uzungeren) and five reservoirs (two of which are coastal—Poroy and Mandra, and three are inland—Duvanli, Sinyata Reka, and Zhrebchevo) can be found in the Inventory of Bulgarian wetlands [40]. Here we only note that: (1) coastal lake Durankulak is a protected area and is included in the Red List of Bulgarian Wetlands [40] in the category Critically Endangered; (2) the reservoir Zhrebchevo is included in the Appendix 1 “List of complex and significant reservoirs” of the Water Act (State Gazette 67/1999) [40]; (3) the reservoirs Mandra and Zhrebchevo are large (3366 and 1851 ha, respectively), while the reservoirs Duvanli, Poroy and Sinyata Reka are small (27, 223 and 6 ha, respectively); (4) The reservoirs Zhrebchevo, Duvanli and Sinyata Reka are plain (200–600 m a.s.l.) waterbodies, while Durankulak, Poroy, Vaya, Mandra and Uzungeren are lowland lakes and reservoirs (<0–200 m a.s.l.) [40]—Table 3.
These water bodies were chosen according to the results obtained from a larger study of cyanotoxins and their producers in 28 Bulgarian waterbodies in relation to human health risks and national security problems [53,54]. In the context of this study, sampling sites were chosen after observation of each waterbody by a drone DJI Mavic 2 Enterprise Dual Pro (DJI Technology Co., Ltd. Shenzhen, China), supplied by a photo camera and capable of measuring the surface water temperature. Therefore, from two of the selected waterbodies we sampled at two sites (Table 3). The phytoplankton sampling was performed from inflatable boats. Data on geographical coordinates, altitude, water temperature, pH, total dissolved solids, oxygen concentration, and conductivity were taken by the water monitoring instruments Aquameter AM-200 and Aquaprobe AP-2000 (Aquaread Ltd., Broadstairs, UK). Water transparency was measured with a Secchi disk. The Aqualytic AL410 Photometer from AQUALYTIC® (Dortmund, Germany) was used for the ex situ measurement of total nitrogen (TN) and total phosphorus (TP) [53,54]. All obtained data are shown in Table 3.
The phytoplankton samples were collected in equal amounts of 0.5 L from the surface water layer (0–20 cm) for all types of further studies. For phytoplankton identification and counts they were fixed immediately to a 2% final formalin concentration and transported to the lab, where they were further concentrated by sedimentation to a volume of 50 mL [52,53,54]. The water samples for pigment analyses and PCR studies, within a few hours after collection, were filtered under a mild vacuum through Macherey-Nagel GF5 filters (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany) with porosity 0.4 µm and Whatman 0.45 µm cellulose filters Whatman NC45 ST/Sterile EO (Merck KGaA, Darmstadt, Germany), respectively. The obtained filters were immediately placed in 15 mL sterile plastic tubes (Falcon) and preserved in a dry ice for transportation to the labs and further treatment [53,54].

4.2. Phytoplankton Species Composition and Abundance Assessment by Conventional Light Microscopy (LM)

The phytoplankton species composition was identified following standard taxonomic sources, with updates from Algaebase [55] on non-permanent slides under magnification 100× and immersion using a Motic B1 microscope, supplied by a Moticam 2.0 mp camera with the Motic Images 2 Plus software program (Version 2.0, Moticam, Hong Kong, China, 2017). The subsequent algal counts were done on the same microscope using the Thoma blood-counting chamber, conducting eight counts for each site. The abundance of each species was estimated in both cell numbers and biomass, the last by using the stereometrical approximations and subsequent weight recalculation [53,54,56,57]. Since it was supposed that for the estimation of toxin concentrations the biomass is a better parameter than cell numbers [58], only biomass data were presented here.

4.3. Phytoplankton Composition Assessment by HPLC Marker Pigment Analysis

A pigment analysis of phytoplankton in the water samples was done by HPLC with subsequent application of CHEMTAX [43,44,45] following the standard operational procedure (SOP5), described by Descy [42] in relation to the basic guidelines for detection and monitoring of toxic cyanoprokaryotes. Pigment extraction was carried in 90% HPLC-grade acetone and all processing steps were identical to those of our previous studies in the region [52,53,54].
The use of CHEMTAX allowed estimating the contribution to chlorophyll a of the main algal phyla: Cyanoprokaryota (identified with two pigment types), Chlorophyta and Streptophyta (presented commonly as green algae), Ochrophyta (identified mainly as classes Chrysophyceae and Bacillariophyceae), Cryptophyta, Pyrrhophyta, and Euglenophyta.

4.4. Molecular-Genetic Studies

The metagenomic DNA was isolated from the filters obtained after sampling in the field. Amplification was conducted with a specific pair of primers mult fw (TCACTCGAAATTACCAGAGGAA) and mdn mult rv (GGTGTAATCAAGAAAAGTGCT), designed to specifically amplify the conserved flanking regions of the precursor gene mdnA from the MV cluster [18]. This gene was chosen because it has been proven that it occurs in the majority of gene clusters responsible for the biosynthesis of MVs [8]. This precursor gene is needed for the production of prepeptide (immature peptide, leader peptide), required as a first step in the production of a completely tricyclic N-acetylated MV [8,59].
For the amplification procedure, the steps described in the protocol of the manufacturer of Taq polymerase (My Taq HS mix, BIOLINE, USA Inc., Taunton, MA, USA) were followed. Firstly, a 95 °C denaturation step was applied in the duration of 3 min. It was followed by 35 cycles of denaturation at 95 °C for 15 s, annealing at 59 °C for 30 s, and synthesis at 72 °C for 25 s. The final synthesis was at 72 °C in a duration of 3 min. The PCR mix contained 10µL MyTaqHS Mix (BIOLINE), and 40 nM of each primer. The amplification reaction ran on a thermal cycler (QB-96 apparatus, Quanta Biotech, Byfleet, Surrey, UK).
The amplified fragments were purified from the gel and sent for sequencing by Macrogene Inc. (Seoul, Korea). The obtained results were compared with sequences available in the National Centre for Biotechnology Information (NCBI) genetic database [46] using a search by the standard Basic Local Alignment Search Tool (BLAST) [47]. Then, only the sequences with the highest homology were used for the construction of the phylogenetic tree by application of the program Mega 6.0 and the neighbor-joining method, with 1000 bootstrap replications [60].
Additionally, three mdnA libraries were developed after the cloning of fragments, isolated from the coastal reservoirs Mandra (East and West) and Poroy. Five individual clones from each library were sent for sequencing by Macrogene Inc. (Seoul, Korea) and included in the constructed phylogenetic tree. There, the accession numbers, obtained from NCBI [46] for the new mdnA-based strains (MZ274329-MZ274337) are shown in brackets. The partial short sequence Blu (131 bp) has no NCBI number. It is generally similar to Por1, but has three different nucleotides, which occupy the 75 (G), 81 (T), and 175 (G) positions.

Author Contributions

Conceptualization, M.S.-G., B.U., K.S., and M.R.; sample collection, B.U., M.S.-G., J.-P.D., G.G., and M.R.; investigation, LM, M.S.-G., G.G., and B.U.; analysis of phytoplankton by HPLC, J.-P.D.; molecular-genetic analyses, K.S. and M.R.; resources, B.U.; writing—original draft preparation, M.S.-G.; writing—review and editing, M.S.-G., J.-P.D., G.G., B.U., M.R., and K.S.; project administration, B.U.; funding acquisition, B.U., M.S.-G., M.R., and G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Scientific Research Fund of the Ministry of Education and Science of Bulgaria (SRF-MESB), grant numbers KP-06-OPR06/2/18.12.2018 and KP-06-OPR03/18/19.12.2018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in [NCBI] at [https://www.ncbi.nlm.nih.gov/ (assessed on 26 May 2021)], reference number [MZ274329-MZ274337].

Acknowledgments

The authors would like once more to acknowledge SRF-MESB for the project grants and the European Cooperation in Science and Technology, COST Action ES 1105 ‘CYANOCOST—Cyanobacterial blooms and toxins in water resources: Occurrence, impacts and management’ for adding value to this paper through networking and knowledge sharing with European experts and researchers in the field.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Svirčev, Z.; Lalić, D.; Savić, G.B.; Tokodi, N.; Backović, D.D.; Chen, L.; Meriluoto, J.; Codd, G.A. Global geographical and historical overview of cyanotoxin distribution and cyanobacterial poisonings. Arch. Toxicol. 2019, 93, 2429–2481. [Google Scholar] [CrossRef]
  2. Buratti, F.M.; Manganelli, M.; Vichi, S.; Stefanelli, M.; Scardala, S.; Testai, E.; Funari, E. Cyanotoxins: Producing organisms, occurrence, toxicity, mechanism of action and human health toxicological risk evaluation. Arch. Toxicol. 2017, 91, 1049–1130. [Google Scholar] [CrossRef]
  3. Meriluoto, J.; Spoof, L.; Codd, J. (Eds.) Handbook of Cyanobacterial Monitoring and Cyanotoxin Analysis; John Wiley & Sons, Ltd.: Chichester, UK, 2017. [Google Scholar] [CrossRef]
  4. Falfushynska, H.; Horyn, O.; Osypenko, I.; Rzymski, P.; Wejnerowski, Ł.; Dziuba, M.K.; Sokolova, I.M. Multibiomarker-based assessment of toxicity of central European strains of filamentous cyanobacteria Aphanizomenon gracile and Raphidiopsis raciborskii to zebrafish Danio rerio. Water Res. 2021, 194, 116923. [Google Scholar] [CrossRef]
  5. Welker, M.; Von Döhren, H. Cyanobacterial peptides–Nature’s own combinatorial biosynthesis. FEMS Microbiol. Rev. 2006, 30, 530–563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Janssen, E.M.L. Cyanobacterial peptides beyond microcystins–A review on co-occurrence, toxicity, and challenges for risk assessment. Water Res. 2019, 151, 488–499. [Google Scholar] [CrossRef]
  7. Kust, A.; Řeháková, K.; Vrba, J.; Maicher, V.; Mareš, J.; Hrouzek, P.; Chiriac, M.-C.; Benedová, Z.; Tesařová, B.; Saurav, K. Insight into unprecedented diversity of cyanopeptides in eutrophic ponds using an ms/ms networking approach. Toxins 2020, 12, 561. [Google Scholar] [CrossRef]
  8. do Amaral, S.C.; Monteiro, P.R.; Neto, J.D.S.P.; Serra, G.M.; Gonçalves, E.C.; Xavier, L.P.; Santos, A.V. Current knowledge on microviridin from cyanobacteria. Mar. Drugs 2021, 19, 17. [Google Scholar] [CrossRef] [PubMed]
  9. Ishitsuka, M.O.; Kusumi, T.; Kakisawa, H.; Kunimitsu, K.; Watanabe, M.M. Microviridin. A novel tricyclic depsipeptide from the toxic cyanobacterium Microsystis viridis. J. Am. Chem. Soc. 1990, 112, 8180–8182. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Li, K.; Yang, G.; McBride, J.L.; Bruner, S.D.; Ding, Y. A distributive peptide cyclase processes multiple microviridin core peptides within a single polypeptide substrate. Nat. Commun. 2018, 9, 1780. [Google Scholar] [CrossRef]
  11. Okino, T.; Matsuda, H.; Murakami, M.; Yamaguchi, K. New microviridins, elastase inhibitors from the blue-green alga Microcystis aeruginosa. Tetrahedron 1995, 51, 10679–10686. [Google Scholar] [CrossRef]
  12. Shin, H.J.; Murakami, M.; Matsuda, H.; Yamaguchi, K. Microviridins D-F, serine protease inhibitors from the cyanobacterium Oscillatoria agardhii (NIES-204). Tetrahedron 1996, 52, 8159–8168. [Google Scholar] [CrossRef]
  13. Rohrlack, T.; Christoffersen, K.; Hansen, P.E.; Zhang, W.; Czarnecki, O.; Henning, M.; Fastner, J.; Erhard, M.; Neilan, B.A.; Kaebernick, M. Isolation, characterization, and quantitative analysis of microviridin J, a new Microcystis metabolite toxic to Daphnia. J. Chem. Ecol. 2003, 29, 1757–1770. [Google Scholar] [CrossRef] [PubMed]
  14. Rohrlack, T.; Christoffersen, K.; Kaebernick, M.; Neilan, B.A. Cyanobacterial protease inhibitor microviridin J causes a lethal molting disruption in Daphnia pulicaria. Appl. Environ. Microbiol. 2004, 70, 5047–5050. [Google Scholar] [CrossRef] [Green Version]
  15. Vegman, M.; Carmeli, S. Three aeruginosins and a microviridin from a bloom assembly of Microcystis spp. collected from a fishpond near Kibbutz Lehavot HaBashan, Israel. Tetrahedron 2014, 70, 6817–6824. [Google Scholar] [CrossRef]
  16. Rounge, T.B.; Rohrlack, T.; Nederbragt, A.J.; Kristensen, T.; Jakobsen, K.S. A genome-wide analysis of nonribosomal peptide synthetase gene clusters and their peptides in a Planktothrix rubescens strain. BMC Genom. 2009, 10, 396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Ziemert, N.; Ishida, K.; Liaimer, A.; Hertweck, C.; Dittmann, E. Ribosomal synthesis of tricyclic depsipeptides in bloom-forming cyanobacteria. Angew. Chem. Int. Ed. Engl. 2008, 47, 7756–7759. [Google Scholar] [CrossRef] [PubMed]
  18. Ziemert, N.; Ishida, K.; Weiz, A.; Hertweck, C.; Dittmann, E. Exploiting the natural diversity of microviridin gene clusters for discovery of novel tricyclic depsipeptides. Appl. Environ. Microbiol. 2010, 76, 3568–3574. [Google Scholar] [CrossRef] [Green Version]
  19. Philmus, B.; Christiansen, G.; Yoshida, W.Y.; Hemscheidt, T.K. Post-translational modification in microviridin biosynthesis. ChemBioChem 2008, 9, 3066–3073. [Google Scholar] [CrossRef]
  20. Philmus, B.; Guerrette, J.P.; Hemscheidt, T.K. Substrate specificity and scope of MvdD, a GRASP-like ligase from the microviridin biosynthetic gene cluster. ACS Chem. Biol. 2009, 4, 429–434. [Google Scholar] [CrossRef]
  21. Reshef, V.; Carmeli, S. New microviridins from a water bloom of the cyanobacterium Microcystis aeruginosa. Tetrahedron 2006, 62, 7361–7369. [Google Scholar] [CrossRef]
  22. Funk, M.; Van Der Donk, W. Ribosomal natural products, tailored to fit. Acc. Chem. Res. 2017, 50, 1577–1586. [Google Scholar] [CrossRef] [PubMed]
  23. Fujii, K.; Sivonen, K.; Naganawa, E.; Harada, K. Non-toxic peptides from toxic cyanobacteria, Oscillatoria agardhii. Tetrahedron 2000, 56, 725–733. [Google Scholar] [CrossRef]
  24. Gatte-Picchi, D.; Weiz, A.; Ishida, K.; Hertweck, C.; Dittmann, E. Functional analysis of environmental DNA-derived microviridins provides new insights into the diversity of the tricyclic peptide family. Appl. Environ. Microbiol. 2014, 80, 1380–1387. [Google Scholar] [CrossRef] [Green Version]
  25. Sieber, S.; Grendelmeier, S.M.; Harris, L.A.; Mitchell, D.A.; Gademann, K. Microviridin 1777: A toxic chymotrypsin inhibitor discovered by a metabologenomic approach. J. Nat. Prod. 2020, 83, 438–446. [Google Scholar] [CrossRef]
  26. Welker, M.; Brunke, M.; Preussel, K.; Lippert, I.; von Döhren, H. Diversity and distribution of Microcystis (cyanobacteria) oligopeptide chemotypes from natural communities studies by single-colony mass spectrometry. Microbiology 2004, 150, 1785–1796. [Google Scholar] [CrossRef] [Green Version]
  27. Welker, M.; Christiansen, G.; von Döhren, H. Diversity of coexisting Planktothrix (cyanobacteria) chemotypes deduced by mass spectral analysis of microcystins and other oligopeptides. Arch. Microbiol. 2005, 182, 288–298. [Google Scholar] [CrossRef]
  28. Martins, J.; Saker, M.; Moreira, C.; Welker, M.; Fastner, J.; Vasconcelos, V. Peptide diversity in strains of the cyanobacterium Microcystis aeruginosa isolated from Portuguese water supplies. Appl. Environ. Microbiol. 2009, 82, 951–961. [Google Scholar] [CrossRef]
  29. Martins, J.; Vasconcelos, V. Cyanobactins from cyanobacteria: Current genetic and chemical state of knowledge. Mar. Drugs 2015, 13, 6910–6946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Arnison, P.G.; Bibb, M.J.; Bierbaum, G.; Bowers, A.A.; Bugni, T.S.; Bulaj, G.; Camarero, J.A.; Campopiano, D.J.; Challis, G.L.; Clardy, J.; et al. Ribosomally synthesized and post-translationally modified peptide natural products: Overview and recommendations for a universal nomenclature. Nat. Prod. Rep. 2014, 30, 108–160. [Google Scholar] [CrossRef]
  31. Sandonato, B.B.; Santos, V.G.; Luizete, M.F.; Bronzel, J.L., Jr.; Eberlin, M.N.; Milagre, H.M.S. MALDI imaging mass spectrometry of fresh water cyanobacteria: Spatial distribution of toxins and other metabolites. J. Braz. Chem. Soc. 2017, 28, 521–528. [Google Scholar] [CrossRef]
  32. Pereira, D.; Pimenta, A.; Giani, A. Profiles of toxic and non-toxic oligopeptides of Radiocystis fernandoii (Cyanobacteria) exposed to three different light intensities. Microbiol. Res. 2012, 167, 413–421. [Google Scholar] [CrossRef] [PubMed]
  33. Murakami, M.; Sun, Q.; Ishida, K.; Matsuda, H.; Okino, T.; Yamaguchi, K. Microvirdins, elastase inhibitors from the cyanobacterium Nostoc minutum (NIES-26). Phytochemistry 1997, 45, 1197–1202. [Google Scholar] [CrossRef]
  34. Dehm, D.; Krumbholz, J.; Baunach, M.; Wiebach, V.; Hinrichs, K.; Guljamow, A.; Tabuchi, T.; Jenke-Kodama, H.; Süssmuth, R.; Dittmann, E. Unlocking the spatial control of secondary metabolism uncovers hidden natural product diversity in Nostoc punctiforme. ACS Chem. Biol. 2019, 14, 1271–1279. [Google Scholar] [CrossRef] [PubMed]
  35. Gärtner, G.; Stoyneva-Gärtner, M.; Uzunov, B. Algal toxic compounds and their aeroterrestrial, airborne and other extremophilic producers with attention to soil and plant contamination: A review. Toxins 2021, 13, 322. [Google Scholar] [CrossRef]
  36. Rohrlack, T.; Christiansen, G.; Kurmayer, R. Putative antiparasite defensive system involving ribosomal and nonribosomal oligopeptides in Cyanobacteria of the genus Planktothrix. Appl. Environ. Microbiol. 2013, 79, 2642–2647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kaebernick, M.; Rohrlack, T.; Christoffersen, K.; Neilan, B.A. A spontaneous mutant of microcystin biosynthesis: Genetic characterization and effect on Daphnia. Environ. Microbiol. 2001, 3, 669–679. [Google Scholar] [CrossRef]
  38. Czarnecki, O.; Henning, M.; Lippert, I.; Welker, M. Identification of peptide metabolites of Microcystis (Cyanobacteria) that inhibit trypsin-like activity in planktonic herbivorous Daphnia (Cladocera). Environ. Microbiol. 2006, 8, 77–87. [Google Scholar] [CrossRef]
  39. Ghosh, S.K.; Bagchi, D.; Bagchi, S.N. Proteolytic activity in Microcystis aeruginosa PCC7806 is inhibited by a trypsin-inhibitory cyanobacterial peptide with a partial structure of microviridin. J. Appl. Phycol. 2008, 20, 1045–1052. [Google Scholar] [CrossRef]
  40. Michev, T.; Stoyneva, M. (Eds.) Inventory of Bulgarian Wetlands and Their Biodiversity; Elsi-M: Sofia, Bulgaria, 2007. [Google Scholar]
  41. Stoyneva-Gärtner, M.P.; Descy, J.-P.; Latli, A.; Uzunov, B.; Pavlova, V.; Bratanova, Z.l.; Babica, P.; Maršálek, B.; Meriluoto, J.; Spoof, L. Assessment of cyanoprokaryote blooms and of cyanotoxins in Bulgaria in a 15-years period (2000–2015). Adv. Oceanogr. Limnol. 2017, 8, 131–152. [Google Scholar] [CrossRef] [Green Version]
  42. Descy, J.P. SOP5: Estimation of cyanobacteria biomass by marker pigment analysis. In Handbook of Cyanobacterial Monitoring and Cyanotoxin Analysis; Meriluoto, J., Spoof, L., Codd, J., Eds.; John Wiley & Sons, Ltd.: Chichester, UK, 2017; pp. 343–349. [Google Scholar]
  43. Sarmento, H.; Descy, J.-P. Use of marker pigments and functional groups for assessing the status of phytoplankton assemblages in lakes. J. Appl. Phycol. 2008, 20, 1001–1011. [Google Scholar] [CrossRef]
  44. Mackey, M.D.; Mackey, D.J.; Higgins, H.W.; Wright, S.W. CHEMTAX—A program for estimating class abundances from chemical markers: Application to HPLC measurements of phytoplankton. Mar. Ecol. Prog. Ser. 1996, 144, 265–283. [Google Scholar] [CrossRef] [Green Version]
  45. Wright, S.W.; Jeffrey, S.W. Pigment markers for phytoplankton production. In Marine Organic Matter: Biomarkers, Isotopes and DNA. The Handbook of Environmental Chemistry; Volkman, J.K., Ed.; Springer: Berlin/Heidelberg, Germany, 2006; Volume 2N, pp. 71–104. [Google Scholar]
  46. NCBI: National Centre for Biotechnology Information (NCBI). Available online: https://www.ncbi.nlm.nih.gov/ (accessed on 12 May 2021).
  47. BLAST: Basic Local Alignment Search Tool (BLAST). Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi (accessed on 12 May 2021).
  48. WebLogo. Available online: http://weblogo.berkeley.edu/ (accessed on 27 May 2021).
  49. Radkova, M.; Stefanova, K.; Uzunov, B.; Gärtner, G.; Stoyneva-Gärtner, M. Morphological and molecular identification of microcystin-producing cyanobacteria in nine shallow Bulgarian water bodies. Toxins 2020, 12, 39. [Google Scholar] [CrossRef] [Green Version]
  50. Dimitrova, R.; Nenova, E.; Uzunov, B.; Shishiniova, M.; Stoyneva, M. Phytoplankton composition of Vaya Lake (2004–2006). Bulg. J. Agric. Sci. 2014, 20 (Suppl. 1), 165–172. [Google Scholar]
  51. Descy, J.-P.; Stoyneva-Gärtner, M.P.; Uzunov, B.A.; Dimitrova, P.H.; Pavlova, V.T.; Gärtner, G. Studies on cyanoprokaryotes of the water bodies along the Bulgarian Black Sea Coast (1890–2017): A review, with special reference to new, rare and harmful taxa. Acta Zool. Bulgar. 2018, 43–52. [Google Scholar]
  52. Stoyneva-Gärtner, M.P.; Uzunov, B.A.; Descy, J.-P.; Gärtner, G.; Draganova, P.H.; Borisova, C.I.; Pavlova, V.; Mitreva, M. Pilot application of drone-observations and pigment marker detection by HPLC in the studies of CyanoHABs in Bulgarian inland waters. Mar. Freshw. Res. 2019, 71, 606–616. [Google Scholar] [CrossRef]
  53. Stoyneva-Gärtner, M.; Stefanova, K.; Descy, J.-P.; Uzunov, B.; Radkova, M.; Pavlova, V.; Mitreva, M.; Gärtner, G. Microcystis aeruginosa and M. wesenbergii were the primary planktonic microcystin producers in several Bulgarian waterbodies (August 2019). Appl. Sci. 2021, 11, 357. [Google Scholar] [CrossRef]
  54. Uzunov, B.; Stefanova, K.; Radkova, M.; Descy, J.-P.; Gärtner, G.; Stoyneva-Gärtner, M. Microcystis species and their toxigenic strains in phytoplankton of ten Bulgarian wetlands (August 2019). Botanica 2021, 27, 77–94. [Google Scholar]
  55. AlgaeBase. Available online: http://www.algaebase.org/ (accessed on 26 May 2021).
  56. Rott, E. Some results from phytoplankton counting intercalibration. Schweiz. Z. Hydrol. 1981, 43, 34–62. [Google Scholar] [CrossRef]
  57. Stoyneva, M.P.; Descy, J.-P.; Vyverman, W. Green algae in Lake Tanganyika: Is morphological variation a response to seasonal changes? Hydrobiologia 2007, 578, 7–16. [Google Scholar] [CrossRef]
  58. Padisák, J.; Chorus, I.; Welker, M.; Maršálek, B.; Kurmayer, R. Laboratory analyses of cyanobacteria and water chemistry. In Toxic Cyanobacteria in Water. A Guide to Their Public Health Consequences, Monitoring and Management, 2nd ed.; Chorus, I., Welker, M., Eds.; CRC Press: London, UK, 2021; pp. 689–743. [Google Scholar]
  59. Li, K.; Condurso, H.; Li, G.; Ding, Y.; Bruner, S. Structural basis for precursor protein-directed ribosomal peptide macrocyclization. Nat. Chem. Biol. 2016, 12, 973–979. [Google Scholar] [CrossRef] [PubMed]
  60. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Toxins 13 00448 g001 550
Figure 1. The general phytoplankton species composition in eight waterbodies in Central and Eastern Bulgaria sampled in August 2019. The number of identified taxa and their percentage contribution to the total phytoplankton biodiversity in each phylum is indicated.
Figure 1. The general phytoplankton species composition in eight waterbodies in Central and Eastern Bulgaria sampled in August 2019. The number of identified taxa and their percentage contribution to the total phytoplankton biodiversity in each phylum is indicated.
Toxins 13 00448 g001
Toxins 13 00448 g002 550
Figure 2. The relative contribution of Cyanoprokaryota to the phytoplankton biodiversity (total number of species) by sites in the studied Bulgarian non-lotic waterbodies (right axis) and number of certainly identified by light microscopy morphospecies of Microcystis (left axis). Abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Figure 2. The relative contribution of Cyanoprokaryota to the phytoplankton biodiversity (total number of species) by sites in the studied Bulgarian non-lotic waterbodies (right axis) and number of certainly identified by light microscopy morphospecies of Microcystis (left axis). Abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Toxins 13 00448 g002
Toxins 13 00448 g003 550
Figure 3. The relative contribution of Cyanoprokaryota to the total phytoplankton biomass, estimated from microscope counts, in the phytoplankton of the studied Bulgarian non-lotic waterbodies (August 2019). Legend abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Figure 3. The relative contribution of Cyanoprokaryota to the total phytoplankton biomass, estimated from microscope counts, in the phytoplankton of the studied Bulgarian non-lotic waterbodies (August 2019). Legend abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Toxins 13 00448 g003
Toxins 13 00448 g004 550
Figure 4. The cyanoprokaryote contribution to the total phytoplankton biomass according to the pigment markers, analyzed by HPLC (expressed as percentage contribution to chlorophyll a, calculated using CHEMTAX) and chlorophyll a concentration in the studied non-lotic Bulgarian waterbodies (August 2019). Abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Figure 4. The cyanoprokaryote contribution to the total phytoplankton biomass according to the pigment markers, analyzed by HPLC (expressed as percentage contribution to chlorophyll a, calculated using CHEMTAX) and chlorophyll a concentration in the studied non-lotic Bulgarian waterbodies (August 2019). Abbreviations: DRN E—Lake Durankulak (site East), DRN W—Lake Durankulak (site West), MND E—reservoir Mandra (site East), MND W—reservoir Mandra (site West), UZNG—Lake Uzungeren, ZHRB—reservoir Zhrebchevo, DVNL—reservoir Duvanli, SNRK—reservoir Sinyata Reka.
Toxins 13 00448 g004
Toxins 13 00448 g005 550
Figure 5. The neighbor-joining phylogenic tree constructed with sequences, obtained from the three libraries from the coastal reservoirs Mandra (sites East and West) and Poroy, combined with the sequences from the PCR fragments, amplified from all studied Bulgarian waterbodies. The bootstrap value is a percentage of 1000 resamplings. The new obtained sequences are supplied with the NCBI accession numbers MZ274329-MZ274337. For the identical sequences (IS), obtained during this study, only one accession number is provided in each cluster or subcluster as follows: (1) The IS from Poroy (Por 1, 3-5)—MZ274329; (2) The IS from Durankulak (Dur W) and Mandra (Man (E) 4)—MZ274331; (3) The IS from Mandra (Man (W) 1,2,4 and Man (E) 1,2,5)—MZ274332; (4) The IS from Duvanlii (Duv) and identical sequences from Uzungeren (Uz) and Zhrebchevo (Zh)—MZ274336. The partial sequence Blu (131 bp) from the reservoir Sinyata Reka (=Blue River) has no NCBI number, but is generally similar to Por1 (explanations are provided in the text).
Figure 5. The neighbor-joining phylogenic tree constructed with sequences, obtained from the three libraries from the coastal reservoirs Mandra (sites East and West) and Poroy, combined with the sequences from the PCR fragments, amplified from all studied Bulgarian waterbodies. The bootstrap value is a percentage of 1000 resamplings. The new obtained sequences are supplied with the NCBI accession numbers MZ274329-MZ274337. For the identical sequences (IS), obtained during this study, only one accession number is provided in each cluster or subcluster as follows: (1) The IS from Poroy (Por 1, 3-5)—MZ274329; (2) The IS from Durankulak (Dur W) and Mandra (Man (E) 4)—MZ274331; (3) The IS from Mandra (Man (W) 1,2,4 and Man (E) 1,2,5)—MZ274332; (4) The IS from Duvanlii (Duv) and identical sequences from Uzungeren (Uz) and Zhrebchevo (Zh)—MZ274336. The partial sequence Blu (131 bp) from the reservoir Sinyata Reka (=Blue River) has no NCBI number, but is generally similar to Por1 (explanations are provided in the text).
Toxins 13 00448 g005
Toxins 13 00448 g006 550
Figure 6. Comparison of the amino acid alignment of the translated mdnA sequences, isolated from eight Bulgarian waterbodies (red font), with the alignments of the translated mdnA sequences in NCBI strains of M. aeruginosa and M. viridis (black font), processed by the Vector NTI v.10 software package. Below is a part of the consensus sequence generated by WebLogo [48]. The specific motif PFFARFL from the α-helix of the leader peptide and the core motif KYPSD are outlined. Legend: Blu—reservoir Sinyata Reka (=Blue River); Dur (E)—Lake Durankulak, site East; Dur (W)—Lake Durankulak, site West; Duv—reservoir Duvanli; ManE—reservoir Mandra, site East; Man (W)—Mandra West part; Por—reservoir Poroy; Vai—Lake Vaya; Uz—Lake Uzungeren; Zh—reservoir Zhrebchevo.
Figure 6. Comparison of the amino acid alignment of the translated mdnA sequences, isolated from eight Bulgarian waterbodies (red font), with the alignments of the translated mdnA sequences in NCBI strains of M. aeruginosa and M. viridis (black font), processed by the Vector NTI v.10 software package. Below is a part of the consensus sequence generated by WebLogo [48]. The specific motif PFFARFL from the α-helix of the leader peptide and the core motif KYPSD are outlined. Legend: Blu—reservoir Sinyata Reka (=Blue River); Dur (E)—Lake Durankulak, site East; Dur (W)—Lake Durankulak, site West; Duv—reservoir Duvanli; ManE—reservoir Mandra, site East; Man (W)—Mandra West part; Por—reservoir Poroy; Vai—Lake Vaya; Uz—Lake Uzungeren; Zh—reservoir Zhrebchevo.
Toxins 13 00448 g006
Toxins 13 00448 g007 550
Figure 7. Map of R. Bulgaria with locations of the studied waterbodies (in case of different sampling sites, they are labelled as E and W) (modified after Ginkgo Maps (http://www.ginkgomaps.com (assessed on 26 May 2021)) and Google Maps (https://www.google.com/maps (assessed on 26 May 2021)).
Figure 7. Map of R. Bulgaria with locations of the studied waterbodies (in case of different sampling sites, they are labelled as E and W) (modified after Ginkgo Maps (http://www.ginkgomaps.com (assessed on 26 May 2021)) and Google Maps (https://www.google.com/maps (assessed on 26 May 2021)).
Toxins 13 00448 g007
Table
Table 1. Distribution of Microcystis taxa identified by LM and their contribution to the total phytoplankton biomass in the studied Bulgarian waterbodies (WBs) in August 2019: Res.—reservoir; MA—Microcystis aeruginosa; MC—Microcystis cf. comperei; MN—Microcystis natans; MW—Microcystis wesenbergii; SS/DC—separate cells or disintegrated colonies; TTs—colonies with transitional morphology; n.d.—not detected.
Table 1. Distribution of Microcystis taxa identified by LM and their contribution to the total phytoplankton biomass in the studied Bulgarian waterbodies (WBs) in August 2019: Res.—reservoir; MA—Microcystis aeruginosa; MC—Microcystis cf. comperei; MN—Microcystis natans; MW—Microcystis wesenbergii; SS/DC—separate cells or disintegrated colonies; TTs—colonies with transitional morphology; n.d.—not detected.
WaterbodyMAMCMNMWSS/DCTTs
Lake Durankulak East<0.5%n.d.n.d.n.d.<0.05%.n.d.
Lake Durankulak Westn.d.n.d.n.d.<1%<0.05%.<0.05%
Lake Vayan.d.n.dn.dn.d<0.05%n.d.
Res. Poroyn.d.n.d.n.d.<0.5%.n.d.n.d.
Res. Mandra East<5%n.d.n.d.<0.5%n.d.<0.05%
Res. Mandra West<0.1%n.d.n.d.<0.5%<0.05%.<0.05%
Lake Uzungerenn.d.n.dn.dn.d<0.05%n.d.
Res. Zhrebchevo<0.5%n.d.<0.05%.n.d.<0.05%.n.d.
Res. Duvanlin.d.<0.01%<0.1%<1%<0.1%<0.05%
Res. Sinyata Reka<0.5%n.d.n.d.<1%<0.05%.<0.05%
Table
Table 2. Comparison of standard published microviridin sequences [8,18] with the microviridin sequences, obtained from eight lowland Bulgarian waterbodies (represented in red font) according to the alignment of the amino acids before (Z1–Z3) and after (EE(D)Z4) the TX KYPSD W motif [8,18]. Those marked with an asterisk (*) indicate the difference in letters, published by different authors (Y and D in [8], and F and E in [18]). The dark yellow color marks the core motif sequence KYPSD, by yellow color are shown the not changing parts of the sequences. Identical sequences are colored in dark green, and similar but not completely identical sequences are shown in light green. The grey color marks differences in the similar sequences. Legend: Blu—reservoir Sinyata Reka (=Blue River); Dur (E)—Lake Durankulak, site East; Dur (W)—Lake Durankulak, site West; Duv—reservoir Duvanli; ManE—reservoir Mandra, site East; Man (W)—Mandra West part; Por—reservoir Poroy; Vai—Lake Vaya; Uz—Lake Uzungeren; Zh—reservoir Zhrebchevo.
Table 2. Comparison of standard published microviridin sequences [8,18] with the microviridin sequences, obtained from eight lowland Bulgarian waterbodies (represented in red font) according to the alignment of the amino acids before (Z1–Z3) and after (EE(D)Z4) the TX KYPSD W motif [8,18]. Those marked with an asterisk (*) indicate the difference in letters, published by different authors (Y and D in [8], and F and E in [18]). The dark yellow color marks the core motif sequence KYPSD, by yellow color are shown the not changing parts of the sequences. Identical sequences are colored in dark green, and similar but not completely identical sequences are shown in light green. The grey color marks differences in the similar sequences. Legend: Blu—reservoir Sinyata Reka (=Blue River); Dur (E)—Lake Durankulak, site East; Dur (W)—Lake Durankulak, site West; Duv—reservoir Duvanli; ManE—reservoir Mandra, site East; Man (W)—Mandra West part; Por—reservoir Poroy; Vai—Lake Vaya; Uz—Lake Uzungeren; Zh—reservoir Zhrebchevo.
Microviridin Variants/Amino AcidsZ1Z2Z3TXKYPSDWEE/DZ4
Microviridin AYGGTFKYPSDWEEY
Por 2YGGTLKYPSDWEEY
MV Microcystis aeruginosa NIES 843YGGTFKYPSDWEDY
MV Microcystis viridis NIES 102YGGTFKYPSDWEEY
Microviridin B/CFGTTLKYPSDWEEY
Blu (and M. aeruginosa NIES 298)FGTTLKYPSDWEEY
Microviridin D/KY*(F)GNTMKYPSDWED*(E)Y
Microviridin E/F FSTYKYPSDWEDF
Microviridin G/HYPQTLKYPSDWEEY
Microviridin IYPTTLKYPSDWEDY
Microviridin J ISTRKYPSDWEEW
Microviridin LYGGTFKYPSDWEDY
Microviridin SD1684, SD1634, SD1652 TATRKYPSDWEDY
Microviridin LH1667 YSTFKYPSDWEDY
Microviridin 1777YNVTFKYPSDWEDY
CBJ55500.1 (Microcystis 199)YGVTLKYPSDWEEF
Dur (W), Man (E) 1, 2, 4, Man (W) 5, VaiYGVTLKYPSDWEEF
MV from a field sample MV/MC19YNVTLKYPSDWEEF
Dur (E), Man (E) 3, Uz, ZhYNVTLKYPSDWEEF
MV IzancyaFGTTFKYPSDWEDY
Por 1FGTTFKYPSDWEDY
MV Microcystis NIES 100FGTTFKYPSDWEDF
MV Microcystis PCC 9805 TSTRKYPSDWEEF
MV Microcystis NIES 103YGGTFKYPSDWEEY
MV Microcystis PCC 7005GRGTLKYPSDWEES
MV from a field sample YSTRKYPSDWEEF
MV from a field sampleANVTLKYPSDWGEF
MV from a field sampleYGGTLKYPSDWEDY
MV from a field sampleYGSTFKYPSDWEDF
MV from a field sampleYEVTLKYPSDWEEF
Table
Table 3. Sampling sites (organized by order of sampling, 14–20 August 2019) and their main environmental parameters. Legend: IBWXXXX—number of the waterbody in the Inventory of Bulgarian non-lotic wetlands [40]); Alt—altitude (m a.s.l.); Area—in ha; WT—water temperature (°C); SD—Secchi depth (m); EC—electric conductivity (S m−1); TD—total dissolved solids (µg L−1); DO—oxygen concentration (mg L−1); TP—total phosphorus (mg L−1); TN—total nitrogen (mg L−1). For details see the text of the paper.
Table 3. Sampling sites (organized by order of sampling, 14–20 August 2019) and their main environmental parameters. Legend: IBWXXXX—number of the waterbody in the Inventory of Bulgarian non-lotic wetlands [40]); Alt—altitude (m a.s.l.); Area—in ha; WT—water temperature (°C); SD—Secchi depth (m); EC—electric conductivity (S m−1); TD—total dissolved solids (µg L−1); DO—oxygen concentration (mg L−1); TP—total phosphorus (mg L−1); TN—total nitrogen (mg L−1). For details see the text of the paper.
WaterbodyAltLatitudeLongitudeAreaWTpHSDECTDDOTPTN
Lake Durankulak (IBW0216)
West		243°40.0006′29°32.6166′ 26.58.890.60.9746317.860.300.66
East443°40.5355′28°33.0806′36926.78.910.60.1056806.040.330.63
Res. Poroy (IBW3038)4342°43.3403′27°37.5255′22327.58.050.40.6444167.60.100.31
Res. Mandra (IBW1720)
West		742°24.0295′27°19.1194′ 25.887.90.450.6764367.930.660.46
East842°25.9303′27°26.7652′336627.28.460.450.5783757.871.51.8
Lake Vaya (IBW0191)−242°30.5940′27°22.075′246327.99.220.150.490177.690.500.26
Lake Uzungeren (IBW0710)−342°26.1551′27°27.2235′17727.68.450.450.17511329.70.400.28
Res. Zhrebchevo (IBW2545)25342°36.6024′25°51.2345′185127.67.700.70.3582338.010.100.18
Res. Duvanli (IBW1483)25042°23.1851′24°43.1000′2726.38.760.44.0502917.090.100.25
Res. Sinyata Reka (IBW1890)31742°28.1480′24°42.217′627.49.720.50.4703059.36254.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Uzunov, B.; Stefanova, K.; Radkova, M.; Descy, J.-P.; Gärtner, G.; Stoyneva-Gärtner, M. First Report on Microcystis as a Potential Microviridin Producer in Bulgarian Waterbodies. Toxins 2021, 13, 448. https://doi.org/10.3390/toxins13070448

AMA Style

Uzunov B, Stefanova K, Radkova M, Descy J-P, Gärtner G, Stoyneva-Gärtner M. First Report on Microcystis as a Potential Microviridin Producer in Bulgarian Waterbodies. Toxins. 2021; 13(7):448. https://doi.org/10.3390/toxins13070448

Chicago/Turabian Style

Uzunov, Blagoy, Katerina Stefanova, Mariana Radkova, Jean-Pierre Descy, Georg Gärtner, and Maya Stoyneva-Gärtner. 2021. "First Report on Microcystis as a Potential Microviridin Producer in Bulgarian Waterbodies" Toxins 13, no. 7: 448. https://doi.org/10.3390/toxins13070448

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