Spatial and Temporal Diversity of Cyanometabolites in the Eutrophic Curonian Lagoon (SE Baltic Sea)

: This work aims to determine the proﬁles of cyanopeptides and anatoxin synthetized by cyanobacteria in the Lithuanian part of the Curonian Lagoon (SE Baltic Sea) and to characterize their spatial and temporal patterns in this ecosystem. Cyanometabolites were analysed by a LC-MS/MS system and were coupled to a hybrid triple quadrupole/linear ion trap mass spectrometer. During the investigation period (2013–2017), 10 microcystins, nodularin, anatoxin-a, 16 anabaenopeptins, including 1 oscillamide, 12 aeruginosins, 1 aeruginosamide, 3 cyanopeptolins and 4 microginins were detected. The most frequently detected metabolites were found at all investigated sites. Demethylated microcystin variants and anabaenopeptins had the strongest relationship with Planktothrix agardhii , while non-demethylated microcystin variants and anatoxin had the strongest relationship with Microcystis spp. Low concentrations of some microcystins: [Asp 3 ]MC-RR, MC-RR, MC-LR, as well as a few other cyanopeptides: AP-A and AEG-A were found during the cold period (December– March). Over the study period, Aphanizomenon , Planktothrix and Microcystis were the main dominant cyanobacteria species, while Planktothrix , Microcystis , and Dolichospermum were potentially producers of cyanopeptides and anatoxin detected in samples from the Curonian Lagoon.


Introduction
Cyanobacteria are prokaryotes found in a wide variety of aquatic environments. They grow at high densities, forming 'blooms' that increase in extension and frequency, and follow anthropogenic activities and climate changes [1][2][3]. In temperate regions, cyanobacteria 'blooms' usually occur during the summer and autumn periods, while some cyanobacteria can also be found in winter and spring [4][5][6][7]. Cyanobacteria blooms are composed of toxic and nontoxic organisms representing different taxonomic units. Cyanobacteria species from at least 16 different genera, belonging to the orders Oscillatoriales, Nostocales, Chroococcales, and Synechococcales are considered as the main bloom forming species, and more than 100 different cyanobacteria taxa (species or genera) produce toxins [8][9][10]. The mixtures of various secondary metabolites are detected during naturally occurring cyanobacteria 'blooms'. These include low molecular weight toxins, such as anatoxins, saxitoxins, β-methylamino-L-alanine (BMAA) and cylindrospermopsin, and a wide array of oligopeptides, termed cyanopeptides, including microcystins [11]. To date, 2031 cyanobacterial metabolites have been described, among which approximately 65% are peptides [10]. The majority of these peptides range from 400 to 1900 Da, but most of them are present in the 1000-1100 Da range. Cyanopeptides consist of cyclic and linear compounds [10]. Amongst the bioactive cyanopeptides, microcystins are the most geographically widespread and are not restricted to any climatic zone or geographic range [12]. Over 286 microcystin variants have been identified and structurally characterized so far [13,14]. Other peptide classes

Study Area
The Curonian Lagoon is a large (surface area 1584 km 2 ) and shallow (mean depth 3.8 m) estuarine lagoon in the south-eastern (SE) Baltic Sea. The lagoon is mainly a freshwater body connected with the sea by a narrow Klaipėda Strait; brackish water intrusions enter the lagoon irregularly and the salinity in the northern part fluctuates from 0 to 7 psu [35]. The highest annual average salinities of 3-5.5 psu are determined in the Klaipeda Strait. It gradually decreases towards the south and reaches a freshwater salinity level (0.5 psu) at a distance of approximately 35 km from the sea entrance [35]. The highest water temperature is observed in July-August with a monthly average of 20.2 • C [36].
In the Curonian Lagoon, five stations were selected, representing different hydrological conditions. Station 1 (Nida) is situated in the sheltered zone of the central part of the lagoon (

Sampling and Measurement of Environmental Parameters
Sub-surface samples were collected every 1-3 weeks in the Lithuan Curonian Lagoon. In 2013, samples were collected from August to Octob from July to November (63 samples in total) during the cyanobacteria veg (stations 1-4), while from 2015 (starting from October) to 2017 samples we

Sampling and Measurement of Environmental Parameters
Sub-surface samples were collected every 1-3 weeks in the Lithuanian part of the Curonian Lagoon. In 2013, samples were collected from August to October, and in 2014 from July to November (63 samples in total) during the cyanobacteria vegetation period  (Figure 1). Phytoplankton samples were fixed with acidic Lugol's iodine solution and analysed using the Utermöhl method up to species or a higher taxonomic level; phytoplankton abundance and biomass (wet weight) according to the closest geometrical shape were calculated based on HELCOM recommendations [37]. Physical water parameters: water temperature (in situ using Eco-Sense ® ODO200) and salinity (after transporting samples into the laboratory, using Mettler Toledo, FiveEasy Cond meter F30) were measured during each sampling. Ice presence/absence was recorded as well. Chlorophyll a (Chl a), including all phytoplankton taxa, was measured in situ using a Fluorometer (FluoroProbeII, bbe Moldaenke).

Analysis of Cyanometabolites
For the analyses of cyanometabolites, water samples of 100-1000 mL (depending on algae density) were filtered on GF/F filters (ø 47 mm). The material on the filter was extracted with 75% methanol (3 mL) and sonicated with an ultrasonic disruptor for 1 min (HD 2070 Sonopuls, Bandeline, Berlin, Germany; 20 kz, 25% duty cycle) [38], followed by a 10 min sonication in a water bath (Sonorex, Bandelin, Berlin, Germany) [39]. After centrifugation at 12,000 g for 15 min the supernatant was separated and analysed with a Agilent HPLC system (Agilent technologies, Waldboronn, Germany) coupled to a hybrid triple quadrupole/linear ion trap mass spectrometer QTRAP LC-MS/MS (QTRAP5500, Applied Biosystems, Sciex; Canada) as described in [40]. Peptides were separated on a Zorbax Eclipse XDB-C18 column (4.6 × 150 mm; 5 µm) (Agilent Technologies, Santa Clara, California, USA) by gradient elution with a mixture of 5% acetonitrile in water and 100% acetonitrile, both with 0.1% formic acid. To determine the profile of peptides in the samples, an information dependent acquisition mode (IDA) was used. The structure of the detected peptides was suggested based on collected fragmentation spectra and the spectra published by other authors (e.g., [20,41,42], (Table S1)), while their relative quantities were estimated based on peak areas in the LC-MS chromatogram. For the quantitative analysis of microcystins, nodularin (NOD), and anatoxin-a, for whom standards were available, the multiple reaction monitoring mode (MRM) with the following transitions was used: 1038 → 135q (q-quantifier), 1038 → 103 (

Data Analyses
For phytoplankton data analysis, potentially toxic cyanobacteria species were pooled into main genera: Aphanizomenon The redundancy analysis (RDA) was applied using Brodgar (2.7.5) and R (3.3.3) packages to test the relationships between the abiotic factors (water temperature, salinity), the biomass of the different cyanobacteria species (used as explanatory variables) and cyanometabolites produced by cyanobacteria as response variables. All response and biological explanatory variable data were square root-transformed. The Brodgar generated RDA biplots that were interpreted were based on explanatory factor line directions and lengths and response variable lines [43]. Only the most abundant cyanometabolites (found in more than 12% of the samples and NOD) were applied for the statistical analyses.

ER REVIEW 5 of 22
The Baltic Sea water intrusions into the Curonian Lagoon were defined by salinity changes and periodically observed at the stations near the Klaipeda strait. Salinity fluctuations at St. Two and St. Three were dynamic during the summer-autumn period (ranged from 0.1 to 6.6), while at St. Five the fluctuations occurred very often and periodically throughout all the periods of investigation (from 0.2 to 6.6 psu) ( Figure 2).

Diversity of Cyanometabolites at the Study Sites
In total, 48 variants of cyanobacteria metabolites were detected in the Curonian Lagoon during the study period. The compounds belonged to six different peptide classes: microcystins (MC) (including nodularin (NOD), anabaenopeptins (AP), aeruginosines (AER), cyanopeptolins (CP), microginins (Mg), aeruginosamide (AEG), and one alkaloid class-anatoxin (ANTX) (Table S1). Microcystins, anabaenopeptins, and aeruginosins accounted for the highest variety of cyanopeptides, while the lowest variety was observed for the remaining ones. The lowest numbers of peptides were found in the samples collected from the stations influenced by the Nemunas River (St. Four; 21 compounds) and the Baltic Sea (St. Five; 23 compounds), while the highest numbers were detected at the stations situated along the western coast of the lagoon (St. One and St. Two; 35 and 31 compounds, respectively) ( Table S1). The highest diversity of metabolites (40 compounds) was detected in 2014. The Baltic Sea water intrusions into the Curonian Lagoon were defined by salinity changes and periodically observed at the stations near the Klaipeda strait. Salinity fluctuations at St. Two and St. Three were dynamic during the summer-autumn period (ranged from 0.1 to 6.6), while at St. Five the fluctuations occurred very often and periodically throughout all the periods of investigation (from 0.2 to 6.6 psu) ( Figure 2).

Diversity of Cyanometabolites at the Study Sites
In total, 48 variants of cyanobacteria metabolites were detected in the Curonian Lagoon during the study period. The compounds belonged to six different peptide classes: microcystins (MC) (including nodularin (NOD), anabaenopeptins (AP), aeruginosines (AER), cyanopeptolins (CP), microginins (Mg), aeruginosamide (AEG), and one alkaloid class-anatoxin (ANTX) (Table S1). Microcystins, anabaenopeptins, and aeruginosins accounted for the highest variety of cyanopeptides, while the lowest variety was observed for the remaining ones. The lowest numbers of peptides were found in the samples collected from the stations influenced by the Nemunas River (St. Four; 21 compounds) and the Baltic Sea (St. Five; 23 compounds), while the highest numbers were detected at the stations situated along the western coast of the lagoon (St. One and St. Two; 35 and 31 compounds, respectively) ( Table S1). The highest diversity of metabolites (40 compounds) was detected in 2014.

Station One (Nida)
In total, 34 metabolites were found at St. One with a maximum of 19 different compounds per sample (Figure 3a, Table S1). The highest diversity of cyanometabolites, including the highest number of aeruginosins (up to 10 variants), was observed in 2014. The most common metabolites throughout the study were AEG-A and AP-A (70% of all samples), Osc-Y (58%), AP-B, AP-F, and AP662 in 43% all samples (Table S1).
Nine microcystins and anatoxin-a were found at Station One during the whole study period (Table S1) The highest cyanobacteria biomass (166 mg L −1 ) and Chl a (1065 µg L −1 ) was recorded in September 2014 ( Figure 3c). Microcystis spp. and Planktothrix agardhii were found in all of the samples, including the winter period. P. agardhii constituted about a quarter of the total cyanobacteria biomass on average. The highest biomass of P. agardhii was recorded in October 2014 (8.3 mg L −1 ), while the Microcystis spp. biomass was only 3.2 mg L −1 (July 2014). These cyanobacteria accounted for 26% and 12%, respectively of the total cyanobacteria biomass on average. Aphanizomenon spp. was found in 80% of the samples and accounted for more than 90% of the total cyanobacteria biomass in September-October. Aphanizomenon was more abundant in 2013-2014 as compared to all the investigated periods. The highest biomass was recorded in September 2014 (160.4 mg L −1 ). Woronichinia was found in 97.9% of the total samples; the largest part of the total cyanobacteria biomass occurred in October-March. The highest biomass of Woronichinia was observed in October 2015 (6.9 mg L −1 ).

Station Two (Juodkrantė)
In total, 31 cyanometabolites were found at St. Two with a maximum of 22 different compounds per sample in 2014 ( Figure 4a, Table S1). In terms of toxic cyanometabolites, ten different variants were detected: eight microcystins, anatoxin, and nodularin. Microcystins  (Table S1). Nodularin was detected in September-November 2014; the concentrations varied from 0.06 to 0.12 µg L −1 . ANTX-a was detected in half of the samples, with the highest concentration (0.51 µg L −1 ) measured in August 2013. Nine microcystins and anatoxin-a were found at Station One during the whole study period (Table S1) Table S1). Microcystin MC-RR was de- In terms of other cyanopeptides, 10 anabaenopeptin variants were detected (Table S1). The most common were Osc-Y (found in all of the same samples, except November 2014) and AP-A (found in 81.3% of the samples). Other metabolites were much less common: AP-F, AP-B, and AER 126A were found in 56.3% of the samples, AEG-A and AP676 were found in 50% of the samples (Table S1). Mg 591B was found twice in August and October 2013, while CP978 was found once in October 2014.
The highest phytoplankton biomass and Chl a concentrations were measured in September 2014 (65.3 mg L −1 and 402.5 µg L −1 , respectively), when Aphanizomenon dominated (47.7 mg L −1 ) (Figure 4c). Compared to other cyanobacteria, Aphanizomenon constituted the highest contribution to the total cyanobacteria biomass (up to 97% of the total cyanobacteria biomass). Microcystis spp. and P. agardhii were found in all of the samples. Microcystis constituted from 0.2% to 87.4% of the total cyanobacteria biomass, and the highest biomass was measured in September 2014 (10.9 mg L −1 ). P. agardhii constituted up to 63.8% of the total biomass and reached the maximum of 6.3 mg L −1 in September 2014 as well. Woronichinia was found mainly in 2013, and its highest biomass was observed in October 2013 (0.55 mg L −1 ).

Station Three (Dreverna)
In total, 27 cyanometabolites were found at St. Three with a maximum diversity of 16 different compounds per sample ( Figure 5a, Table S1).  (Table S1) In terms of other cyanopeptides, nine anabaenopeptins and six aeruginosins were found. The most common anabaenopeptins were AP-A, Osc-Y, and AP-F found in 42.3%, 42.3% and 35.7% of the samples, respectively. The Mg 591B was found in September 2013 and from September to October 2014, while CP978 was found only in August 2014.
Station Three is characterised by the lowest phytoplankton biomass and Chl a concentration ( Figure 5c). The highest Chl a concentration was recorded in October 2014 (76 µg L −1 ), while the highest cyanobacteria biomass (10.5 mg L −1 ), dominated by Aphanizomenon (9.2 mg L −1 ), was recorded in early September 2014. During the study period, the total biomass of P. agardhii did not exceed 1.8 mg L −1 , except September 2013 (4.5 mg L −1 ).

Station Four (Ventė)
Out of all of the investigated stations, the lowest variety of cyanometabolites was found in samples from Station Four, where only 21 compounds in total and a maximum of 14 different metabolites per sample were detected (Figure 6a). The MC-RR variant constituted the largest contribution to the overall bioactive cyanometabolite concentration (Figure 6b In terms of other cyanopeptides, nine anabaenopeptins and six aeruginosins were found. The most common anabaenopeptins were AP-A, Osc-Y, and AP-F found in 42.3%, 42.3% and 35.7% of the samples, respectively. The Mg 591B was found in September 2013 and from September to October 2014, while CP978 was found only in August 2014. Station Three is characterised by the lowest phytoplankton biomass and Chl a con- ]MC-HtyR, ANTX-a and nodularin were found rarely over the investigated period; the concentrations of ANTX-a and nodularin did not exceed 0.08 µg L −1 (Figure 6b).  Among other cyanopeptides, anabaenopeptins were the most dominant compoundsin total, seven anabaenopeptin variants were detected (Figure 6a). The most found anabaenopeptins were AP662, AP676 and Osc-Y (47.1%, 35.3% and 35.3% of the samples, respectively). The aeruginosin and cyanopeptolin variants were found only in 2014.
Aphanizomenon and Microcystis were found in 88.2% of the samples (Figure 6c). The highest biomass of Aphanizomenon was measured in October 2014 (17.7 mg L −1 ) and this coincided with the highest Chl a concentration (101 µg L −1 ). The biomass of Microcystis was much lower compared to Aphanizomenon and did not exceed 2.0 mg L −1 . Dolichospermum was found in 58.8% of the samples. The highest biomass of this genus was measured in August-September 2013 and early August 2014 (about 7.7 mg L −1 ). P. agardhii was found in 76.5% of the samples, and the highest biomass was measured in September 2013 (4.0 mg L −1 ).

Station Five (Smiltynė)
A relatively low diversity of cyanometabolites (23) was also characteristic for Station Five; in total six different microcystin variants were detected (Table S1). In terms of other cyanopeptides, AEG-A, AP-A and Osc-Y were the most common compounds (53.1%, 43.8%, and 40.6% of all samples, respectively). No metabolites were found from late October 2015 to May 2016 except from AEG-A which was detected in December ( Figure 7a, Table S1).
Planktothrix and Woronichinia were the most common cyanobacteria species at Station Five (100% and 96.9% of the samples, respectively) (Figure 7c). The highest biomasses of Woronichinia and Planktothrix were recorded accordingly in late December 2015 (5.2 mg L −1 ) and June 2017 (4.4 mg L −1 ). Aphanizomenon was found in 75% of samples and was the dominant species during the August-October months. The highest biomass was measured in October 2015 (11.6 mg L −1 ) and in September 2017 (9.7 mg L −1 ). The highest Chl a concentration (92.3 µg L −1 ) was recorded in October 2015.

Statistical Analysis
According to the RDA analysis, the environmental factors in combination with the dominant cyanobacteria biomass accounted for 52% of the variance in the concentration of cyanometabolites from the samples taken from the Curonian Lagoon, while two axes explained 83% of the variation (   Bioactive cyanopeptides (MC-RR, MC-LR, MC-YR) and anatoxin ANTX-a were closely associated with the occurrence of Microcystis spp., Dolichospermum spp., and water temperature. All statistical relationships were statistically significant (p < 0.05), however, Microcystis spp. explains 27.73% (F = 14), Dolichospermum spp.-10.02% (F = 3) and water temperature-8.67% (F = 7) of the variation in peptide concentrations. AEG-A had the closest relationship with the genus Aphanizomenon, while the contribution of Aphanizomenon to the variation was 10.76% (F = 5, p < 0.05).
The presence of nodularin in the Curonian Lagoon was very closely connected to the presence of brackish water, which explains 7.17% (F = 6, p = 0.001) of the data distribution; and both of them were negatively related with the water temperature. The impact of different years and Woronichinia had small though statistically significant contributions to the variation of data (contribution 10.45%, F = 6, p = 0.001 and contribution 3.57%, F = 3, p < 0.05, accordingly).

Discussion
In the present study, we provide the first detailed description of the occurrence and the diversity of cyanopeptides and anatoxin in the Curonian Lagoon (Lithuanian part). Cyanobacteria usually produce toxins and other secondary cyanometabolites as a potential defence mechanism against surrounding aquatic organisms. Cyanometabolites could inhibit the growth of some bacteria and can be toxic to zooplankton ( [11,23] and references therein). As many of the compounds are biologically active and frequently toxic, they can pose a risk to the aquatic ecosystem and human health. The problem is global, therefore The presence of nodularin in the Curonian Lagoon was very closely connected to the presence of brackish water, which explains 7.17% (F = 6, p = 0.001) of the data distribution; and both of them were negatively related with the water temperature. The impact of different years and Woronichinia had small though statistically significant contributions to the variation of data (contribution 10.45%, F = 6, p = 0.001 and contribution 3.57%, F = 3, p < 0.05, accordingly).

Discussion
In the present study, we provide the first detailed description of the occurrence and the diversity of cyanopeptides and anatoxin in the Curonian Lagoon (Lithuanian part). Cyanobacteria usually produce toxins and other secondary cyanometabolites as a potential defence mechanism against surrounding aquatic organisms. Cyanometabolites could inhibit the growth of some bacteria and can be toxic to zooplankton ( [11,23] and references therein). As many of the compounds are biologically active and frequently toxic, they can pose a risk to the aquatic ecosystem and human health. The problem is global, therefore extended knowledge about cyanobacteria occurrence and metabolic diversity is needed [44]. Studies on cyanobacteria have been carried out in the Curonian Lagoon for many years, but they have been limited to the structure of the cyanobacterial community and the occurrence of microcystins, nodularin and anatoxin [28,30]. In terms of water quality, the concentrations of microcystins are the parameters for which, according to the WHO and EPA, the recreational values of water bodies are assessed [17,45]. However, other metabolites such as cyanopeptolins, microginins, aeruginosines sometimes have equivalent or even stronger bioactive effects compared to microcystins [11]. Moreover, the water quality in coastal areas can be affected by the water outflows from the eutrophic estuaries/lagoons [32,46]; therefore, an assessment of the diversity and seasonal variation of cyanometabolites can help to improve the estimation of associated risks.
In our study, a shift from the dominance of MC-RR in 2013 and 2014 to the dominance of [Asp 3 ]MC-RR in 2015-2017 was observed (Figures 3-6). However, the highest microcystin concentration was of [Asp 3 ]MC-RR in 2014 (St. One). Microcystin [Asp 3 ]MC-RR is almost by 3-fold more toxic than its methylated variant MC-RR [17]; therefore, this shift can affect the water quality. Microcystin MC-RR was strongly linked with the presence of Microcystis spp. and Dolichopermum spp., while [Asp 3 ]MC-RR is strongly linked with the presence of Planktothrix (Figure 8). Our results confirm the previous findings that Microcystis and Dolichospermum species are the main producers of non-demethylated microcystins, while Planktothrix is usually linked with demethylated forms of microcystins [40,62]. Moreover, previous studies have shown that Microcystis from the Curonian Lagoon contain the mcyE gene responsible for microcystins synthesis [63]. Moreover, concentrations of [Asp 3 ]MC-RR were high in the Curonian Lagoon throughout the study period and are also the highest in other European lakes [60]. The total microcystin concentration recorded in this work exceeded the recommended guideline values for recreational waters provided by US Environmental Protection Agency (US EPA) (8 µg L −1 ) in 12.6 % of cases. The guideline levels were exceeded at all stations except Station Three. However, these concentrations were within the limits provided by the World Health Organization (24 µg L −1 ); the highest total microcystin concentration recorded in the lagoon was 22.8 µg L −1 . Climatic conditions play an important role in shaping the composition of cyanobacteria species in the water body [64]. Calm and hot weather favour A. flosaquae blooms (this study in 2014, [24,64]). On the other hand, during rainy and windy summers, water mixing, and turbidity increases, resulting in A. flosaquae replacement by other cyanobacteria species (this study, y. 2017 [24,64]). Despite the low biomass of cyanobacteria in 2017 (summer maximum 5.1 mg L −1 ), the number of detected metabolites was similar to other years (Table S1). Therefore, it is likely that A. flosaquae has no significant contribution to the number of metabolites detected in the samples. Despite the dominance of A. flosaquae, other species, particularly P. agardhii, Dolichospermum spp., and Microcystis spp. were considered as the main producers of cyanometabolites.
Cyanobacteria of the main genera recorded in the Curonian Lagoon, Aphanizomenon, Microcystis, Planktothrix, and Woronichinia, were also observed during the cold period, and even under the ice cover (this study Figures 3 and 7, [5,6,[65][66][67]. Low concentrations of [Asp 3 ]MC-RR, MC-RR, and MC-LR, as well as AP-A and AEG-A, indicating the presence of cyanobacteria, were found in the Curonian Lagoon during the December-March period at water temperature 2.9 • C, on average. It is known that microcystins accumulate in the tissues of a wide range of aquatic consumers [68,69]. Therefore, it is not surprising that microcystins have been found in sediments, fish, and mussel tissues in the spring [29,33]. The presence of nodularin at concentrations in the range of 0.01 to 0.28 µg L −1 , in the Curonian Lagoon was associated with higher water salinity values. It is known that N. spumigena, the only nodularin producer in the Baltic Sea [70,71], could live both in brackish and fresh waters [41,60]. The presence of N. spumigena and nodularin in the Curonian Lagoon is associated with Baltic Sea water intrusions to the lagoon [28,32]. According to Zemlys et al. [35], saline water intrusions are more frequent in the northern part of the lagoon, while in extreme cases they could reach the central part. During our study, the farthest station where nodularin was detected was near the Nemunas river avandelta (St. Four). Although Baltic Sea strains of Nodularia can grow in culture at 0 g L −1 NaCl at least for 70 h [72], it seems that this species has not established a permanent population in the Curonian Lagoon. The detection of the cell-bound nodularin did not always coincide with the detection of N. spumigena. Therefore, we can assume that the microscopy method may not reveal the presence of rare species of cyanobacteria.
ANTX-a is a relatively less common bioactive cyanometabolite present in waterbodies worldwide. According to Mantzouki et al. [60], it was found in 39% of the investigated 137 European lakes. In the Baltic Sea region, ANTX-a was found in the Gulf of Finland [49], the Gulf of Gdansk [73], and the south-eastern Baltic coast [32], while it was absent in the Baltic Proper [74,75], the Danish straits [76], and the Vistula Lagoon [77]. In the Curonian Lagoon, anatoxin was not found during the snapshot surveys, even during the blooms of different cyanobacteria (e.g., M. aeruginosa, A. flosaquae, P. agardhii) [30,31], while during long-term surveys, the toxin was reported to be present in 31% (this study) or 47% [32] of the investigated samples. During our study, ANTX-a was found in the lagoon from July to October, but not every year. According to the statistical analysis, the closest relationship of this cyanomethabolite was established with the Microcystis genus. However, only few Microcystis strains from Japan and Portugal have been found to produce ANTX-a so far [78][79][80]. It is known that ANTX-a can be also produced by Aphanizomenon, Dolichospermum (previously known as Anabaena), Gomphosphaeria, Limnothrix, Lyngbya, Planktothrix, Planktolyngbia, Pseudanabaena, and Synechocystis ( [81,82] and references therein). In the Lithuanian lakes, the presence of ANTX-a significantly correlated only with the C. issatschenkoi; however, none of the isolated strains of the species showed ANTX-a production [83]. The concentrations of ANTX-a (0-0.515 µg L −1 ) detected in the Curonian Lagoon are similar to the concentrations determined in Lithuanian lakes [83] and other European lakes [60].
There are numerous references on the composition of cyanopeptides in lakes, while estuaries are poorly studied. Estuaries could be more diverse in species composition, because of the presence of both fresh and fresh-brackish water organisms [6]; therefore, we could expect more potential cyanopeptide producers. In many estuaries and coastal waters, Microcystis spp. and Dolichospermum (previously known as Anabaena) are considered as the main potential microcystin producers ( [46] and references therein). Species representing Aphanizomenon and Planktothrix genera are less common in estuaries, despite the fact that Aphanizomenon thrive in brackish waters (e.g., Baltic Sea, [84]), meanwhile Planktothrix can grow and produce toxins in waters up to 7.5 g L −1 of NaCl [85]. In the Curonian lagoon, depending on the year and season, Aphanizomenon, Planktothrix, Microcystis, Dolichospermum or Woronichinia can have the largest contribution to the total cyanobacterial biomass.
Until recently, P. agardhii was observed only occasionally in the Curonian Lagoon and reached less than 5% of the total phytoplankton abundance ( [86] and references therein). Based on the available literature, from the late 1990s, P. agardhii are found all year round (this study, [6,65]) and could account for 13-57% of the total cyanobacteria biomass [25,33,87]. Planktothrix tolerates fluctuations in fresh-brackish waters, turbid and or shady conditions [5,64], which are characteristic to the Curonian Lagoon too. Demethylated microcystin variants ([Asp 3 ]MC-RR and [Asp 3 ]MC-LR) are usually associated with the presence of P. agardhii [40,62]. In previous studies [28,30], only [Asp 3 ]MC-RR was found in the Curonian Lagoon during the dominance of Planktothrix. Microcystins [Asp 3 ]MC-RR and [Asp 3 ]MC-LR were the most common toxins in our study. However, we cannot unequivocally attribute the presence of the compounds to a particular cyanobacterium. It has been shown that P. agardhii from the Curonian Lagoon has the mcyE gene, responsible for toxin production [88].
Usually, more attention is given to widely known bioactive cyanopeptides (e.g., microcystins) and anatoxins, their composition, quantity, and production regulation, etc. Nevertheless, other cyanopeptides could be more harmful to humans and aquatic organisms [11]. In the Curonian Lagoon, anabaenopeptins showed the closest statistical coincidence with P. agardhii and had the highest number of variants, of which the most frequent were AP-A, AP-B, AP-F, and Osc-Y. These peptides are also commonly detected in other freshwater bodies dominated by P. agardhii [40,89], Microcystis spp. [89][90][91], Coelosphaeriaceae (e.g., Woronichinia compacta) [92], in the strain of Anabaena90 isolated from the fresh environment [93]. The same metabolites AP-A, Osc-Y, AP915, AP807, and AP-D were found in the brackish environment (Gulf of Gdansk, south Baltic Sea), dominated by N. spumigena (50% of the cyanobacterial biomass), A. flosaquae (40%), and Dolichospermum spp. (10%) [94]. According to Janssen ([11] and references therein), in European lakes the concentration of anabaenopeptins usually is higher as compared to microcystins.
A. flosaquae is responsible for the largest and longest cyanobacteria blooms in the Curonian Lagoon (e.g., this study, [6,24,95]). So far, there are no data on its ability to produce microcystins either in the Baltic Sea, or in the Curonian Lagoon [12,30], meanwhile, it is considered as the producer of other classes of bioactive metabolites [12]. Despite the weak relationship between Aphanizomenon abundance and concentrations of cyanopeptides, the highest correlation based on the RDA analysis was observed with the aeruginosamide AEG-A. This confirms the findings indicating that A. flosaquae and C. issatschenkoi are correlated with aeruginosamide [96]. During previous studies in the Curonian Lagoon, AEG-A was found in both samples from A. flosaquae and in M. aeruginosa dominated communities [31]. However, in that study, AEG-A was found in 52% of the samples collected all year round and in samples dominated by different cyanobacteria: Aphanizomenon, Dolichospermum, Microcystis, P. agardhii, Woronichinia. Aeruginosamide is also commonly found in waterbodies dominated by Microcystis spp. [56,89]. In Poland aeruginosamide was found in several freshwater samples dominated by P. agardhii [40]. Recently, 18 different aeruginosamide variants were found in a Limnoraphis CCNP1324 culture isolated from the brackish Baltic Sea [97].
Aeruginosins also have a number of structural variants. To date, about 80 variants of the compound from Microcystis, Nodularia, and Planktothrix have been described [10,98]. In phytoplankton samples collected from the Curonian Lagoon, 12 different aeruginosin variants were detected in the summer-autumn period, but they were not common. Typically, the highest number of aeruginosin variants were recorded when Aphanizomenon had the highest biomass. Aeruginosins with m/z 618 and m/z 568 were found in previous research conducted in the Curonian Lagoon during the dominance of A. flosaquae and M. aeruginosa [31]. Recently, it was found that A. flosaquae strains from USA and Finland have aeruginosin biosynthetic gene clusters [99]; therefore, it is possible that in the Curonian Lagoon, A. flosaquae is one of the aeruginosin producers as well. The other two species, P. agardhii and Microcystis spp., which are known to produce aeruginosins [10,98], usually coexist as subdominant species in the Curonian Lagoon. The most common, AER 126A and AER 126B, were previously isolated from P. agardhii [100] and were found in a P. agardhii dominated freshwater pond as well (Poland, [40]).
Cyanopeptolins and microginins are two other groups with approximately 210 and 72 variants, respectively [12,101]. According to the literature, there are several species in the Curonian Lagoon that can synthesize cyanopeptolins, namely Microcystis, Planktothrix, and Dolichospermum [12], while microginins are synthesized mainly by Microcystis, Planktothrix, Synechococcus, and Woronichinia naegeliana [12,101,102]. The amount of data on the occurrence of these two groups of peptides is rapidly increasing in recent years (e.g., [99] and references therein, [40,56,89,103]), however, they are still under-researched compounds [11,101]. During the whole period of our study, we found only three cyanopeptolins of different variants and four microginin variants. Moreover, both cyanopeptolins and microginins were found in the Curonian Lagoon, only occasionally in the summerautumn period and in less than 4% of all samples, except for microginin Mg 591B, which was identified in 12% of the samples. M. aeruginosa strains isolated from Portuguese freshwaters were the source/producers of 11 cyanopeptolins and three microginin variants [104]; from Greek freshwaters, 51 different microginin variants were detected [101]. In the Curonian Lagoon, cyanopeptolins were detected during the dominance of Microcystis and P. agardhii, while microginins were present during the predominance of P. agardhii. The low number of cyanopeptolins and microginins found in the Curonian Lagoon might be attributed to the narrowly specific cyanobacteria species / strain's ability to produce these two compounds, or to unfavourable estuarine conditions for the production of the compounds.
The eutrophic, dynamic, and heterogeneous environment of the Curonian Lagoon is favourable for the growth, diversity, and high abundance of different cyanobacteria species and, consequently, supports the diversity and high concentration of cyanopeptides. Over the past decades, a shorter or absent ice period and increasing water temperature have been observed during the winter in the Curonian Lagoon [105]. Due to this reason, cyanobacteria (especially overwintering species such as Planktothrix or Limnothrix) and their produced cyanometabolites may be present in water bodies during the winter period [5]. This, in turn, led to the detection of a greater variety and higher concentrations of cyanopeptides, also in winter. Moreover, aeruginosins, anabaenopeptins, cyanopeptolins, and microginins can be as toxic as microcystins [11]. These cyanopeptides may have synergistic effects on aquatic organisms when found together [11,89]. In general, the effects of cyanopeptides on human health have been more frequently investigated than the ecological function of these metabolites in the water ecosystems [106,107].
Despite the fact that Aphanizomenon was the dominant cyanobacterium species in our study of the Curonian Lagoon, it was Planktothrix, Microcystis, and Dolichospermum that were the main producers of cyanopeptides and anatoxin in the water body. The cyanometabolites richness shown in our work can form a good starting point for further studies on their role in ecosystem functioning, including their indirect impact on water quality at bathing sites. Knowledge of the effects of environmental conditions on the production and fate of these compounds are important to the general understanding of water ecosystem dynamics.