Next Article in Journal
Detection of Extremely Low Level Ciguatoxins through Monitoring of Lithium Adduct Ions by Liquid Chromatography-Triple Quadrupole Tandem Mass Spectrometry
Next Article in Special Issue
Microbial Community Response to Granular Peroxide-Based Algaecide Treatment of a Cyanobacterial Harmful Algal Bloom in Lake Okeechobee, Florida (USA)
Previous Article in Journal
Effects of a Curcumin/Silymarin/Yeast-Based Mycotoxin Detoxifier on Redox Status and Growth Performance of Weaned Piglets under Field Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microbial Diversity Impacts Non-Protein Amino Acid Production in Cyanobacterial Bloom Cultures Collected from Lake Winnipeg

by
Stephanie L. Bishop
1,2,*,†,
Julia T. Solonenka
1,†,
Ryland T. Giebelhaus
1,3,4,
David T. R. Bakker
1,
Isaac T. S. Li
1 and
Susan J. Murch
1
1
Department of Chemistry, University of British Columbia, Syilx Okanagan Nation Territory, Kelowna, BC V1V 1V7, Canada
2
Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
3
Department of Chemistry, University of Alberta, Edmonton, AB T6G 2N4, Canada
4
The Metabolomics Innovation Centre, Edmonton, AB T6G 2N4, Canada
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxins 2024, 16(4), 169; https://doi.org/10.3390/toxins16040169
Submission received: 14 February 2024 / Revised: 10 March 2024 / Accepted: 20 March 2024 / Published: 26 March 2024

Abstract

:
Lake Winnipeg in Manitoba, Canada is heavily impacted by harmful algal blooms that contain non-protein amino acids (NPAAs) produced by cyanobacteria: N-(2-aminoethyl)glycine (AEG), β-aminomethyl-L-alanine (BAMA), β-N-methylamino-L-alanine (BMAA), and 2,4-diaminobutyric acid (DAB). Our objective was to investigate the impact of microbial diversity on NPAA production by cyanobacteria using semi-purified crude cyanobacterial cultures established from field samples collected by the Lake Winnipeg Research Consortium between 2016 and 2021. NPAAs were detected and quantified by ultra-performance liquid chromatography–tandem mass spectrometry (UPLC-MS/MS) using validated analytical methods, while Shannon and Simpson alpha diversity scores were determined from 16S rRNA metagenomic sequences. Alpha diversity in isolate cultures was significantly decreased compared to crude cyanobacterial cultures (p < 0.001), indicating successful semi-purification. BMAA and AEG concentrations were higher in crude compared to isolate cultures (p < 0.0001), and AEG concentrations were correlated to the alpha diversity in cultures (r = 0.554; p < 0.0001). BAMA concentrations were increased in isolate cultures (p < 0.05), while DAB concentrations were similar in crude and isolate cultures. These results demonstrate that microbial community complexity impacts NPAA production by cyanobacteria and related organisms.
Key Contribution: Recent studies show that inter-species interactions may be associated with the production of neurotoxic non-protein amino acids (NPAAs) by cyanobacteria and related organisms. Our study highlights that the differences in NPAA concentrations in cyanobacterial bloom cultures are linked to the complexity of the microbial communities present in mixed-species cultures.

Graphical Abstract

1. Introduction

Harmful cyanobacterial blooms are a global public health concern and consistently affect Lake Winnipeg in Manitoba, Canada [1,2]. Lake Winnipeg is the tenth largest freshwater lake in the world with two distinct basins (north and south) separated by a narrow channel [3]. In the mid-1990s, an increased input flow from the Red River watershed incited an ecological state change of the lake [4,5,6]. Lake Winnipeg has been subsequently impacted by heightened agricultural, industrial, and urban activities, in addition to greater spring nutrient loads and rising global temperatures that have accelerated eutrophication [7,8,9]. These factors have resulted in increased cyanobacterial bloom frequency and severity [2,6,10,11,12,13]. The Lake Winnipeg Research Consortium Inc. (Gimli, MB, Canada) was founded in August 1998 to coordinate scientific research on the lake following the 1997 Red River flood [14].
One of the main concerns arising from Lake Winnipeg cyanobacterial blooms is the production of cyanobacterial toxins, including non-protein amino acids (NPAAs) which are associated with adverse health effects [1,15,16,17]. β-N-methylamino-L-alanine (BMAA) is a NPAA first identified in cycad seeds from Guam collected by the ethnobotanist Marjorie Grant Whiting [18,19,20]. Further investigation showed that the BMAA found in cycads was attributable to the presence of symbiotic cyanobacteria growing on cycad roots [21,22]. This resulted in the hypothesis that BMAA could bioaccumulate in ecosystems and food webs, thereby leading to a higher consumption by humans than previously thought [23]. Since then, this NPAA and its isomers 2,4-diaminobutyric acid (DAB), N-(2-aminoethyl)glycine (AEG), and β-aminomethyl-L-alanine (BAMA) have been detected in all five cyanobacterial sections from fresh water and marine ecosystems around the world [17,22,24]. More than 50 publications have documented the presence of BMAA in over 30 cyanobacterial genera since 2005; however, discrepancies in the literature have highlighted the importance of using validated and fit-for-purpose methods to analyze NPAAs in cyanobacterial matrices [25,26,27].
The presence of NPAAs in Lake Winnipeg is concerning due to the wide variety of neurotoxic effects they display, including misincorporation into proteins [21,28,29,30,31] and inducing oxidative stress in several neuronal cell lines through mechanisms including glutamate receptor excitotoxicity, cystine/glutamate antiporter inhibition, and modulation of the canonical Wnt signaling pathway [32,33,34,35,36,37]. Since these compounds are generally found in relatively low concentration in the environment, they may accumulate in tissues over time to produce neurotoxic effects [21]. The neurotoxicity of BMAA and DAB has been demonstrated with in vitro and in vivo models including rodents, birds, primates, and zebrafish [20,38,39,40]. In one study, vervets fed a low dose of BMAA (210 mg/kg/day) over 140 days developed neuropathy—including neurofibrillary tangles and β-amyloid deposits—consistent with neurodegeneration [40]. Therefore, understanding the factors that influence NPAA accumulations in areas with dense human populations that experience cyanobacterial blooms is necessary to mitigate health concerns due to the risk of long-term exposure.
Presently, the biosynthetic route to BMAA and other NPAAs in cyanobacteria is under investigation. NPAAs are likely ancient molecules and have been found in early earth geobiological samples such as meteorites and microbialites [41,42]. DAB is particularly abundant in nature and can be produced by several species of plants and microbes [41,43,44]. Recently, a bioinformatics approach was used to explore mechanisms for the 2,3-diaminopropanoic acid-dependent biosynthesis of BMAA in plants, bacteria, and cyanobacteria [44]. This approach was complicated by variations in analytical methods and cyanobacterial culture state of the data available; of the 28 cyanobacterial studies used, only 25% explicitly used all-axenic cultures, while others employed mono-cyanobacterial cultures, environmental samples, or did not specify. However, achieving axenicity in algal cultures is markedly difficult due to the presence of obligate symbionts and heterotrophic “helper” bacteria, and even with flow cytometry and 16S rRNA amplicon sequencing it is possible to miss contaminants [45]. As such, the impact of contaminating microbial species on NPAA biosynthesis in cyanobacterial cultures is unknown.
Our previous collaboration with the Lake Winnipeg Research Consortium characterized the concentration profiles of NPAAs by cyanobacteria and related organisms in field samples collected in 2016 from diverse sites around Lake Winnipeg [17]. The objective of the current study was to understand how microbial diversity affects NPAA production in non-axenic cyanobacterial cultures. We hypothesized that cyanobacterial cultures from Lake Winnipeg would produce different amounts of NPAAs pre- and post- attempted purification of bacterial contaminants. To investigate our hypothesis, we tracked culture quality during the course of semi-purification and used 16S rRNA amplicon sequencing to determine the metagenomic composition of each culture. We then used a validated analytical method to determine NPAA concentrations in crude and semi-purified cultures, and correlated these concentrations to the alpha diversity in the cultures [17,46]. The results of this study offer an insight into how microbial community composition impacts NPAA profiles in cyanobacteria.

2. Results

2.1. Establishment and Semi-Purification of Cyanobacterial Cultures from Lake Winnipeg

In total, 27 crude cyanobacterial cultures were established from the Lake Winnipeg collections from 2016–2019. Semi-purification of these cultures and additional samples collected in 2021 yielded 16 additional isolate cyanobacterial cultures. Visual observation of cyanobacterial cultures using an inverted microscope revealed mixed various cyanobacteria from the orders Synechococcales (individual spherical unicellular cells), Chroococcales (stacked colonies consisting of 2–4 oblong cells), and various plant and algae-like organisms distributed across the lake (Figure S1). Two of the samples were identified as Prochlorococcus marinus. The records of qualitative culture tracking (Figure S2) and metagenomic compositions (Figure S3) are provided in the Supplementary Materials; classifications were identical at the bootstrap levels of 50, 80, and 95. Although complete axenicity was not achieved in all isolates, a quarter of our isolated cultures contained >95% cyanobacterial DNA by 16S rRNA amplicon sequencing, and the proportion of cyanobacterial operational taxonomic units (OTUs) was relatively greater in the isolated cyanobacteria compared to the crude cultures.

2.2. Distribution of NPAAs in Crude and Isolate Cyanobacterial Cultures

2.2.1. Concentration Profiles of NPAAs across Lake Winnipeg

The NPAA concentration profiles across Lake Winnipeg varied by isomer, as well as the relative purity of the culture (Figure 1; Table S1). All established cultures produced detectable levels of total AEG. The highest levels of AEG in crude cultures (681 ± 50 ng/g) were produced by cyanobacteria grown from collections at Station 9 and in isolate cyanobacterial cultures from Station 5 NS (248 ± 53 ng/g), with both locations in the south basin of the lake. BAMA was detected in 24 of 27 crude cultures (89%) and in 14 of 16 isolate cultures (88%). The highest levels of BAMA were found in crude cultures established from Station 7NS (83 ± 27 ng/g) in the south basin, and in the isolate cultures from Station W4 (188 ± 97 ng/g) in the north basin of the lake. BMAA was detected in 52% (14 of 27) of the crude cyanobacterial cultures and in 63% (10 of 16) of the isolate cultures. The highest levels of BMAA in crude cultures (484 ± 53 ng/g) were established from Station 3C in the south basin of the lake, and in isolate cultures from Station W1 (73 ± 35 ng/g) in the north basin of the lake. DAB was measured in all cyanobacterial cultures, with the highest levels in crude cultures (8.4 ± 0.55 µg/g) found in a collection from Station 20 in the north basin of the lake and isolated cultures from Station 5 NS in the south basin (18.7 ± 2.34 µg/g).

2.2.2. Occurrence of NPAAs in Paired Crude–Isolate Cyanobacterial Cultures

Four of the crude and isolate cultures were derived from the same parent field sample, corresponding to a cyanobacterial bloom sample collected in 2019, and samples from Station W9, Station 20, and Station 19 collected in 2018 (shown on the maps in Figure 1). The phylogenetic composition of the cultures pre- and post-semi-purification are shown in Figure 2. As expected, in all four crude–isolate pairs, the percentage of reads corresponding to cyanobacteria was relatively higher in the isolate cultures compared to the crude cultures. The four paired crude–isolate cultures were all putatively identified as Chroococcales via morphological inspection; samples Bloom 2019, Station W9 2018, and Station 20 2018 contained oblong cells forming stacks of 2–4, while sample Station 19 2018 formed three-dimensional clusters of spherical cells. Samples Bloom 2019 and Station W9 2018 appeared identical in size (approximately 5–6 microns), while the cells in sample Station 20 2018 were similar, albeit slightly larger in size (approximately 6–8 microns) (Figure S1). Furthermore, samples Bloom 2019 and Station W9 2018 were putatively identified as containing the same cyanobacterial species based on morphological inspection (Figure S1).
A comparison of NPAA concentration profiles between crude and isolate cultures showed that statistically higher levels of AEG were produced in crude cultures (average 138 ± 16 ng/g) compared to isolate cultures (average 43 ± 7.8 ng/g; p < 0.0001; Figure 3a). AEG levels were statistically higher in the crude cultures of all four crude–isolate pairs. BMAA concentrations were statistically higher in crude cultures (average 136 ± 18 ng/g) compared to isolate cultures (average 25 ± 5.3 ng/g; p < 0.0001; Figure 3c). BMAA was detected in three of the isolate cultures and two of the paired crude cultures, with a significant decrease in the isolate culture established from Station 20 in the north basin. Conversely, BAMA concentrations were statistically higher in isolate cultures (average 53 ± 9.7 ng/g) compared to crude cultures (average 31 ± 3.0 ng/g; p < 0.05; Figure 3b). The levels of BAMA were generally higher in the paired crude cultures compared to isolate cultures, but not statistically different. DAB was found at similar levels in isolate cultures (average 2.5 ± 0.7 µg/g) compared to crude cultures (average 2.7 ± 0.2 µg/g; Figure 3d). The levels of DAB were significantly lower for paired isolate cultures established from the bloom site in 2019, as well as Stations 19 and 20.

2.3. Analysis of NPAA Profiles and Alpha Diversity in Cyanobacterial Cultures

At a general level, the alpha diversity of the isolate cultures was significantly decreased compared to the crude cultures (p < 0.0001; Figure 4a). We then compared the influence of Shannon diversity scores on NPAA concentrations in crude and isolate cultures via the Pearson correlation test and a two-way analysis of covariance (ANCOVA; α = 0.05). The Pearson correlation test revealed that AEG was moderately correlated with Shannon diversity (r = 0.554; p < 0.0001) and all other NPAAs were weakly correlated with Shannon diversity (r < 0.250). The two-way ANCOVA revealed a significant influence of both culture type (p < 0.0001) and Shannon diversity (p = 0.014) on AEG concentrations. Culture type also had a significant influence on BMAA concentrations (p = 0.026) and DAB concentrations (p = 0.004), but none of BAMA, BMAA, or DAB were significantly influenced by Shannon diversity scores. The ordination plot showing the beta diversity (Bray–Curtis distances) revealed clustering of the crude cultures, indicating similarity in species composition (Figure 4b). Additionally, the principal component analysis biplot of NPAA concentrations showed that AEG and BMAA were associated with crude cultures and a higher Shannon diversity, while BAMA and DAB were associated with isolate cultures and a lower Shannon diversity (Figure 4c). None of the plots showed specific patterns related to the lake section of the initial collection site (statistical ellipses not shown); rather, data clustering was associated with the culture type (Figure S4).

3. Discussion

Recently, there has been growing interest in exploring the role of environmental factors in the production of BMAA and other NPAAs by cyanobacteria and organisms present in complex ecosystems [47]. The cultivation of axenic algal and cyanobacterial cultures presents many challenges [45], but our data indicate that our strategy of the qualitative visual tracking of culture purity corresponded appropriately with our metagenomic findings. This showed that relatively successful cyanobacterial purification was accomplished without the use of additional chemical or instrumental strategies. Previous studies have shown that NPAA production in aquatic organisms is species-specific [43,48]. While we were not able to identify the primary species associated with the production of NPAAs, the metagenomic analysis of paired crude–isolate cultures presented here suggests that the relative complexity of a mixed-species culture influences NPAA levels. Additionally, our beta diversity analysis showed clustering of crude cultures compared to isolate cultures, highlighting the basal microbial community present within Lake Winnipeg [49]. The paired crude–isolate cultures from the same collection location contained residual bacteria from the genera Dietzia spp. and Alkalibacterium spp., both of which are chemoorganotrophic species [50,51]. AEG, BAMA, and BMAA are not well-characterized bacteria, but some microbial species are known producers of DAB, which could explain the higher concentration of DAB detected in most samples compared to its isomers [43,44].
Additionally, our data provide insights into the effects of microbial community composition on the production of NPAAs. Comparison of these data to our previous work showed that concentrations of NPAAs in uncultured field samples were orders of magnitude higher than when these field samples were cultured in vitro [17]. However, the concentrations for all in vitro cultures in this study as well as our previous study were much more similar, suggesting that growing cyanobacteria in optimized laboratory conditions in a nutritionally replete growth medium may alter the ability of some cyanobacterial strains to produce detectable BMAA [52,53]. The concentrations of NPAAs in environmental and lab-cultured samples vary greatly between studies and these results are further complicated by the wide range of analytical methods used to detect these compounds in biological matrices [25,27]. A recent study indicated that BMAA concentrations in cyanobacterial strains isolated from Northern Polish water bodies ranged from not detected (ND) to 0.8 µg/g, while cultured cyanobacterial strains ranged from ND to 0.6 µg/g (dry weight (DW); total BMAA) [54]. Our findings fall within the range of these values, as total BMAA concentrations ranged from ND to 0.48 µg/g DW and furthermore indicate that BMAA production in cyanobacteria is not strongly correlated to the amount of contaminant bacteria in culture, as BMAA concentrations were not significantly influenced by the alpha diversity of the culture. This is in line with previously published in vitro diatom studies showing that the lack of culture axenicity did not significantly affect concentrations of BMAA and other NPAAs in these species [43].
Researchers are using both in vitro and field monitoring approaches to uncover the dynamics between cyanobacteria and other BMAA-producing organisms in complex microbial communities [55,56,57,58,59]. BMAA metabolism in cyanobacteria is hypothesized to be closely associated with environmental nitrogen levels and nitrogen fixation [24,60,61]. Many recent studies have examined the effects of nitrogen availability on NPAA production by cyanobacteria and related organisms using in vitro techniques or by tracking bloom samples over a season [24,60,61,62]. During non-favorable environmental conditions, such as a bloom outbreak or collapse, BMAA may accumulate and act as an antagonistic agent to provide a competitive advantage for non-nitrogen fixing cyanobacteria [15,63,64,65,66,67]. Additionally, certain microbial phyla may have evolved to co-exist with NPAA-producing organisms using specific mechanisms of protection. A new line of research shows that the gut microbiota of cycad-feeding insects, consisting primarily of Proteobacteria, produce iron-chelating metabolites that allow these insects to tolerate toxin-rich diets including BMAA [68].
In summary, the metabolic functions of NPAAs including BMAA have not been fully elucidated, and more research is needed to understand the role of these compounds in harmful algal blooms. Researchers are beginning to develop new in vitro approaches, such as co-culture systems and isotope tracing methods, to analyze metabolite exchange in multi-species microbial communities [69,70]. These emerging strategies could be used to trace the metabolic production pathways and exchange of NPAAs in controlled multi-species communities that include cyanobacteria and related organisms. These approaches would complement existing NPAA field monitoring strategies [61]. Elucidating the specific factors—including nutrient availability and antagonistic or mutualistic microbial interactions—that affect cyanotoxin production in multi-species communities where numerous cyanobacteria species co-exist in varying abundances will lead to the development of better mitigation strategies to reduce cyanotoxin exposures in areas experiencing harmful algal blooms.

4. Materials and Methods

4.1. Sample Collection and Cyanobacterial Culture Establishment

4.1.1. Field Sample Collection

Samples were collected by the Lake Winnipeg Research Consortium in 50 mL conical centrifuge tubes (Corning CentriStar; Corning, New York, NY, USA) from the surface (0–0.5 m depth) of standard stations and from cyanobacterial bloom sites during the voyage of the research ship Namao in the summer and fall seasons of 2016–2019, and 2021. Two samples were collected from each location during the summer and fall voyages of the Namao. All samples were identified by the station or GPS location of collection. The samples were shipped at an ambient temperature for the establishment of in vitro cultures of the cyanobacteria. The samples were couriered to the University of British Columbia Okanagan within 1 week of the voyage.

4.1.2. Crude Cyanobacterial Culture Establishment

Upon arrival to the University of British Columbia Okanagan, field samples were briefly vortexed to homogenize cells. Using sterile inoculation loops (Nunc™ Disposable Loops and Needles; Thermo Scientific, Waltham, MA, USA), 100 µL of cells from each sample were spread onto Petri dishes (VWR, Mississauga, ON, Canada) containing complete BG11 media (Sigma Aldrich, Oakville, ON, Canada) solidified with agar (Fisher Scientific, Ottawa, ON, Canada; 4 g agar/500 mL media) and autoclaved (Steris, Mentor, OH, USA; 121 °C at 18 psi for 20 min), as described previously [17]. Crude cultures were incubated in a controlled-environment growth room at 27 °C with a 16 h photoperiod (Sunblaster™ full spectrum daylight fluorescent bulbs; 40 μmol/m2/s) and subcultured monthly by streaking onto fresh BG11 agar plates using sterile inoculation loops. After five months of sequential culture, crude cyanobacterial cultures were transferred into bioreactor boxes (Liquid Lab Vessel; Caisson Laboratories, Smithfield, VA, USA), to which 50 mL liquid BG11 media were added. These boxes were incubated in the same controlled-environment room with regular replenishment of media to mitigate volume loss due to evaporation.

4.1.3. Isolate Cyanobacterial Culture Establishment

To generate isolate cyanobacterial cultures, crude cultures were screened at 3× magnification (Nikon SMZ745 stereo microscope, Brighton, MI, USA), and cyanobacterial colonies were subsequently isolated and cultured onto individual Petri dishes with sterile inoculation loops. Isolate cultures were incubated in the previously described controlled-environment room and were routinely monitored at 3× magnification for the presence of contaminating bacteria. The culture quality was ranked qualitatively on a scale of 1–5, with 1 indicating heavy contamination (i.e., no regions of uncontaminated cyanobacterial growth) and 5 meaning no visible contamination (Figure S2). When the culture growth appeared to plateau, indicated by no further increases in colony size, they were checked under the stereo microscope and clean regions were marked for subculturing to gradually purify the samples of contaminants.

4.2. Detection and Quantification of Non-Protein Amino Acids

4.2.1. Sample Preparation

The total NPAA acid content in crude cultures was measured by pipetting approximately 1.5 mL of cells from each bioreactor box or Petri dish into a pre-weighed microcentrifuge tube (1.5 mL; Eppendorf, Hamburg, Germany), which was frozen at −20 °C. Frozen cells, including their extracellular matrices, were lyophilized to complete dryness (FreeZone 4.5; Labconco, Kansas City, MI, USA). Sample preparation, detection, and quantification followed the previously published validated analytical methods, with the data collection time extended to include BAMA [17,46]. In brief, approximately 10–20 mg of lyophilized cells was accurately weighed using an analytical balance (Ohaus; VWR, West Chester, PA, USA) and transferred into a glass hydrolysis vial (15 mm × 45 mm; Fisher Scientific). 1000 µL of 6 N HCl were added to the hydrolysis vial, purged with nitrogen gas for 30 s, then hydrolyzed at 110 °C for 18 h on a heating block (VWR Standard Dry Block Heater). After cooling, 400 µL of sample was pipetted into a centrifuge filter tube (0.2 µm nylon membrane; VWR) and centrifuged (VWR Galaxy 16DH centrifuge) for 5 min at 13,000 rpm (16,000× g). A 10 µL aliquot of the filtrate was placed into a polypropylene centrifuge tube (1.5 mL; Fisher Scientific) and dried to complete dryness in a Speedvac (Labconco Centrivap, VWR). The dried-down filtrate was then reconstituted with a borate buffer (pH 10) and derivatized (see Section 4.2.2).
For separate measurement of protein-bound and free NPAA content in isolate cultures, approximately 10–20 mg of fresh cells were collected directly from the Petri dish into pre-weighed microcentrifuge tubes and analyzed without lyophilization. The free amino acid portion was extracted by pipetting 200 μL of 0.1 N TCA to each sample and vortex-mixing at a maximum speed for 30 s. The samples were then centrifuged for 5 min at 13,000 rpm and the supernatant was pipetted into a centrifuge filter tube and centrifuged for 5 min at 13,000 rpm to yield the free amino acid extract. The pellet was resuspended in 100 μL 0.1× phosphate-buffered saline and transferred to a glass hydrolysis vial; 900 μL of 6 N HCl was added to the vial and the rest of the procedure followed as previously stated to collect the protein-bound amino acid portion.

4.2.2. Derivatization Reaction and Quantification of NPAAs

The sample filtrates were reconstituted with 80 µL a 0.2 M borate buffer and derivatized with 20 µL AccQ-Fluor reagent (Waters Corp., Mississauga, ON, Canada; reconstituted as per the kit instructions). Complete reaction was ensured by vortex-mixing (Vortex Genie 2; Scientific Industries, Bohemia, NY, USA) and incubation for 10 min at 55 °C prior to analysis. Lysine was also monitored during analysis to ensure complete derivatization [71]. Samples were quantified by comparison to authentic standards (β-N-methylamino-L-alanine: Sigma Aldrich, CAS No. 16012-55-8; N-(2-aminoethyl)glycine: Sigma Aldrich, CAS No. 24123-14-6; 2,4-diaminobutyric acid: TCI America, CAS No. 1883-09-06), and a custom synthesis of β-aminomethyl-L-alanine described in the supplemental section to [17], derivatized in the same manner for each batch of samples, as described previously. The 100 µL aliquot of derivatized sample was then transferred to an autosampler vial (2 mL amber glass with pre-slit Teflon-coated caps; Waters Corp.) fitted with a conical bottom spring insert (250 µL glass; Canadian Life Science, Peterborough, ON, Canada) for separation and detection.
All samples were analyzed by UPLC-MS/MS (Waters Acquity I-Class BSM, SM-FTN, column heater/cooler; Waters Xevo TQ-S with ESCi probe, MassLynx v4.1) with electrospray ionization (ESI+) and multiple reaction monitoring (MRM) on a Waters Acquity UPLC BEH C18 (2.1 mm × 100 mm; 1.7 µm particle size) column, as described in detail previously. A gradient elution of 20 mM ammonium formate at pH of 5.0 adjusted with glacial acetic acid (Fisher) (Solvent A) was used, and methanol (Fisher, Optima™ LC/MS grade) (Solvent B) with initial conditions of 90:10, reaching 50:50 at 7.00 min, 25:75 at 7.50 min and re-equilibrating to 90:10 at 9.60 min to 12 min before the next injection. The needle and seal wash solvent consisted of 90% acetonitrile (Fisher, Optima™ LC/MS grade) and 10% water (Millipore Ultrapure, 18.2 MΩ cm). The method developed by Glover et al. (2015) was fully validated by a Single Laboratory Validation protocol and there was no significant difference in the performance characteristics of the current analysis. Multiple reaction monitoring (MRM) transitions and other mass spectrometry parameters are presented in Table S2. The complete baseline separation of all analytes was achieved, ensuring accurate quantification of each analyte. All standard curves had excellent linearity, with R2 > 0.99 in the linear dynamic range (0.37–1530 pg on column for each analyte). The method detection limit (MDL) was calculated by the U.S. Environmental Protection Agency guidelines, defined as the minimum measured concentration of a substance that can be reported with 99% confidence that the measured concentration is distinguishable from the method blank results [72]. The MDLs using the spiked blank method for AEG, BAMA, BMAA, and DAB were determined to be 0.34, 0.67, 0.48, and 1.19 pg on the column, respectively. The limits of quantification (LOQ) for AEG, BAMA, BMAA, and DAB were 1.17, 2.33, 1.67, and 4.11 pg on the column, respectively. The MDL for BMAA in a cyanobacterial matrix (Nostoc sp.) was 5.82 ng/g (dry weight). Due to endogenously high levels of AEG, BAMA, and DAB in this matrix, producing a signal-to-noise ratio of >20 in the samples, we could not calculate the matrix MDL for those analytes.

4.3. Taxonomic Analysis of Cyanobacterial Cultures and Data Visualization

4.3.1. Microscopic Images and Visual Identification of Cyanobacterial Orders

Images were acquired using an inverted microscope (Olympus IX83; Olympus Canada Inc., Toronto, ON, Canada) at 40× magnification. The transmission images of the crude and isolate cultures were obtained using full color DIC (Differential Interference Contrast), while fluorescence images of the crude cultures were obtained by a non-staining epifluorescence method with an mCherry dichroic filter cube (excitation: 562 ± 20 nm, emission: 641 ± 37.5 nm). Microscopic images were edited using open-source ImageJ software (v1.51). The dominant cyanobacterial order/s present within each isolate culture were putatively identified using the morphological characteristics described in [73], the Cyanobacterial Key developed by Robin A. Matthews and Geoffrey B. Matthews at Western Washington University in 2018 (http://snoringcat.net/cyanobacteria_key/index.html, accessed on 25 April 2023), the cyanobacterial gallery maintained by Dr. Jason Oyadomari at Finlandia University (https://www.keweenawalgae.mtu.edu/index.htm, accessed on 25 April 2023), and the taxonomic key in the 2018 “A Guide to Cyanobacteria; Identification and Impact” [74].

4.3.2. Next-Generation Sequencing (NGS) of Cyanobacterial Cultures

DNA Extraction

Cyanobacterial genomic DNA was extracted from the cultures using the Soil DNA Isolation Kit (Norgen Biotek #64000, Thorold, Schmon Pkwy, Canada) with some additional modifications to increase disruption of the cyanobacterial membranes. The frozen cell pellets were resuspended in 750 μL kit-provided Lysis Buffer G and vortexed for several seconds; the resuspensions were then transferred to kit-provided bead tubes. Of the kit-provided Lysis Additive A, 200 µL and 50 μL of 0.2 μg/μL lysozyme were then added, and the samples were vortexed at a maximum speed for 10 min, followed by sonication for 5 min. The samples were centrifuged (16,000× g, 2 min) and the supernatant transferred to a fresh tube. Of the Binding Buffer I, 100 µL was added, and the tubes were inverted 8× and incubated at 4 °C for 5 min. They were centrifuged (16,000× g, 4 min) and the supernatant again transferred. Of the OSR solution, 50 µL was added and the tubes were inverted 8× and incubated at 4 °C for 5 min. They were centrifuged (16,000× g, 4 min), and 700 μL of the supernatant was transferred to fresh tubes. Of the Lysis Buffer QP, 400 µL and 550 μL of ethanol were added, and the samples were vortexed and then passed through kit-provided spin columns (3 × 11,000× g, 30 s). Of the Binding Buffer B, 550 µL was passed through the spin columns (11,000× g, 1 min), then they were washed with 500 μL Wash Solution A (2 × 11,000× g, 1 min). The columns were dried by centrifugation (16,000× g, 4 min) and the columns transferred to the elution tubes, to which 50 μL Elution Buffer B was added. The elution tubes were incubated for 1 min at room temperature, and finally the genomic DNA was collected by elution (11,000× g, 1 min). Eluted gDNA was stored at −20 °C until it was sent for next-generation sequencing (NGS) analysis.

16S rRNA Sequencing

Illumina MiSeq 16S rRNA sequencing of the samples was performed at the Sequencing and Bioinformatics Consortium at the University of British Columbia (Vancouver, BC, Canada). A 25 μL aliquot of each gDNA sample was submitted. The V3 and V4 regions of the 16S rRNA genes were PCR-amplified using literature-derived bacterial primers and then modified with Illumina sequencing adaptors and dual-index barcodes with the Nextera XT DNA Index kit. Libraries were then pooled and sequenced on a MiSeq Nano flow cell (Illumina, San Diego, CA, USA) to generate paired-end 250 bp reads. Illumina MiSeq Reporter (v2.6.2.3) was used to generate bioinformatic reports with taxonomic information for each sample.

DADA2 Taxonomic Identification

DADA2 (v1.16) was used to process the Illumina paired-end fastq files [75]. In brief, the read qualities of the samples were inspected and filtered; the first 20 and final 10 nucleotides from each sequence were trimmed to remove primer sequences and poor-quality regions. The forward and reverse reads were then merged, and chimeras removed. Taxonomy was assigned using CyanoSeq (v1.1.2), consisting of a curated cyanobacterial 16S rRNA database containing 5411 cyanobacterial sequences [76] and SILVA (v138.1), and a database containing more than 2 million 16S/18S rRNA bacterial, archaiac, and eukaryotic sequences [77]. Taxonomic assignments were performed at the bootstrap levels of 50, 80, and 95. Shannon and Simpson indexes of alpha diversity, and Bray–Curtis distances for beta diversity were calculated using phyloseq (v1.44) [78].

4.3.3. Statistical Analysis and Visualization

Quantification, statistical analysis, and visualization of NPAA data was performed using Microsoft Excel®, R (v4.3.1) in RStudio®, GraphPad Prism (v9.5.1), and Adobe®, Illustrator (v27.5). Data from the chromatographic analysis of NPAAs were collated in Excel® to calculate the concentrations of each isomer in each sample. All averaged concentration data are presented with the standard error of the mean. For statistical purposes, undetected (ND) values were replaced with a value of 1/10th of the method detection limit for each compound and back-calculated to units of ng/g based on the starting sample weight and sample dilution scheme. As none of the raw concentration values for NPAAs followed a normal distribution by the D’Agostino and Pearson test, we log-transformed all values before performing the Pearson correlation test and a two-way analysis of covariance (ANCOVA; α = 0.05). After the log-transformation, all concentration values followed a normal distribution except for BMAA due to the high number of samples with an ND value. The statistical tests used for each analysis are described in the figure captions. A site map and GPS coordinates of the collection sites used for this study were provided by the Lake Winnipeg Research Consortium. NPAA concentrations were site-mapped using the GeoJSON package (v0.3.4) from R [79] and the geomap provided by the University of Manitoba Centre for Earth Observation Science [80]. All other data were plotted using GraphPad Prism or the ggplot2 package [81] in R.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins16040169/s1, Figure S1: Microscopic images of cyanobacterial cultures; Figure S2: Culture quality tracking spanning a period of 445 days; Figure S3: Metagenomic profiles of crude and isolate cultures; Table S1: Total β-N-methylamino-L-alanine (BMAA), 2,4-diaminobutyric acid (DAB), N-(2-aminoethyl)glycine (AEG), and β-aminomethyl-L-alanine (BAMA) concentrations for all bloom and isolate cyanobacterial cultures analyzed in this study, including collection sites and years; Figure S4: Correlations of non-protein amino acid (NPAA) concentrations to culture type; Table S2: Multiple reaction monitoring (MRM) transitions for validated ultra-performance liquid chromatography–tandem mass spectrometry (UPLC-MS/MS) method with the AccQ-Fluor reagent. Supplemental data file: Lake Winnipeg Metagenomic Analysis.xlxs.

Author Contributions

Conceptualization, S.L.B., J.T.S. and S.J.M.; methodology, all authors; validation, S.L.B. and J.T.S.; formal analysis, S.L.B. and J.T.S.; investigation, S.L.B., J.T.S., R.T.G. and D.T.R.B.; resources, I.T.S.L. and S.J.M.; data curation, S.L.B. and J.T.S.; writing—original draft preparation, S.L.B. and J.T.S.; writing—review and editing, all authors; visualization, S.L.B. and J.T.S.; supervision, I.T.S.L. and S.J.M.; project administration, I.T.S.L. and S.J.M.; funding acquisition, S.J.M. All authors have read and agreed to the published version of the manuscript.

Funding

Funding from the Natural Sciences and Engineering Research Council of Canada (grant number #RGPIN-2018-03883) and the Canadian Foundation for Innovation is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Summarized data are available in the Supporting Information or as a supplemental data file. Any additional data are available from the authors by request.

Acknowledgments

We would like to acknowledge Karen Scott, the Lake Winnipeg Research Consortium, and the crew of the Namao for sharing sample collections.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bashir, F.; Bashir, A.; Bouaïcha, N.; Chen, L.; Codd, G.A.; Neilan, B.; Xu, W.L.; Ziko, L.; Rajput, V.D.; Minkina, T.; et al. Cyanotoxins, Biosynthetic Gene Clusters, and Factors Modulating Cyanotoxin Biosynthesis. World J. Microbiol. Biotechnol. 2023, 39, 241. [Google Scholar] [CrossRef] [PubMed]
  2. Carmichael, W.W.; Boyer, G.L. Health Impacts from Cyanobacteria Harmful Algae Blooms: Implications for the North American Great Lakes. Harmful Algae 2016, 54, 194–212. [Google Scholar] [CrossRef] [PubMed]
  3. Brunskill, G.J.; Elliott, S.E.M.; Campbell, P. Morphometry, Hydrology, and Watershed Data Pertinent to the Limnology of Lake Winnipeg; Department of Fisheries and Oceans: Winnipeg, MB, Canada, 1980. [Google Scholar]
  4. Scott, K.J.; Mccullough, G.; Stainton, M.; Ayles, B.; Hann, B. Lake Winnipeg-State of the Science What Is the Scientific Basis for Understanding and Protecting Lake Winnipeg? Lake Winnipeg Research Consortium Workshop Report: Winnipeg, MB, Canada, 2011. [Google Scholar]
  5. Wassenaar, L.I.; Rao, Y.R. Lake Winnipeg: The Forgotten Great Lake. J. Great Lakes Res. 2012, 38, 1–5. [Google Scholar] [CrossRef]
  6. McCullough, G.K.; Page, S.J.; Hesslein, R.H.; Stainton, M.P.; Kling, H.J.; Salki, A.G.; Barber, D.G. Hydrological Forcing of a Recent Trophic Surge in Lake Winnipeg. J. Great Lakes Res. 2012, 38, 95–105. [Google Scholar] [CrossRef]
  7. Bhateria, R.; Jain, D. Water Quality Assessment of Lake Water: A Review. Sustain. Water Resour. Manag. 2016, 2, 161–173. [Google Scholar] [CrossRef]
  8. Cottingham, K.L.; Ewing, H.A.; Greer, M.L.; Carey, C.C.; Weathers, K.C. Cyanobacteria as Biological Drivers of Lake Nitrogen and Phosphorus Cycling. Ecosphere 2015, 6, 1–19. [Google Scholar] [CrossRef]
  9. Spears, B.M.; Carvalho, L.; Perkins, R.; Paterson, D.M. Effects of Light on Sediment Nutrient Flux and Water Column Nutrient Stoichiometry in a Shallow Lake. Water Res. 2008, 42, 977–986. [Google Scholar] [CrossRef] [PubMed]
  10. Zastepa, A.; Pick, F.R.; Blais, J.M. Fate and Persistence of Particulate and Dissolved Microcystin-LA from Microcystis Blooms. Hum. Ecol. Risk Assess. 2014, 20, 1670–1686. [Google Scholar] [CrossRef]
  11. Pick, F.R. Blooming Algae: A Canadian Perspective on the Rise of Toxic Cyanobacteria. Can. J. Fish. Aquat. Sci. 2016, 73, 1149–1158. [Google Scholar] [CrossRef]
  12. Cirés, S.; Ballot, A. A Review of the Phylogeny, Ecology and Toxin Production of Bloom-Forming Aphanizomenon Spp. and Related Species within the Nostocales (Cyanobacteria). Harmful Algae 2016, 54, 21–43. [Google Scholar] [CrossRef]
  13. Paerl, H.W.; Otten, T.G. Harmful Cyanobacterial Blooms: Causes, Consequences, and Controls. Microb. Ecol. 2013, 65, 995–1010. [Google Scholar] [CrossRef]
  14. Scott, K.J. Lake Winnipeg Research Consortium Inc. 2016/2017 Annual Report, Science and Education Programs; Lake Winnipeg Research Consortium: Winnipeg, MB, Canada, 2017. [Google Scholar]
  15. Zhang, Y.; Whalen, J.K. Production of the Neurotoxin Beta-N-Methylamino-l-Alanine May Be Triggered by Agricultural Nutrients: An Emerging Public Health Issue. Water Res. 2020, 170, 115335. [Google Scholar] [CrossRef] [PubMed]
  16. McKindles, K.M.; Zimba, P.V.; Chiu, A.S.; Watson, S.B.; Gutierrez, D.B.; Westrick, J.; Kling, H.; Davis, T.W. A Multiplex Analysis of Potentially Toxic Cyanobacteria in Lake Winnipeg during the 2013 Bloom Season. Toxins 2019, 11, 587. [Google Scholar] [CrossRef]
  17. Bishop, S.L.; Kerkovius, J.K.; Menard, F.; Murch, S.J. N-β-Methylamino-L-Alanine and Its Naturally Occurring Isomers in Cyanobacterial Blooms in Lake Winnipeg. Neurotox. Res. 2018, 33, 133–142. [Google Scholar] [CrossRef]
  18. Whiting, M.G. Toxicity of Cycads. Econ. Bot. 1963, 17, 270–302. [Google Scholar] [CrossRef]
  19. Vega, A.; Bell, E.A. α-Amino-β-Methylaminopropionic Acid, a New Amino Acid from Seeds of Cycas Circinalis. Phytochemistry 1967, 6, 759–762. [Google Scholar] [CrossRef]
  20. Spencer, P.S.; Nunn, P.B.; Hugon, J.; Ludolph, A.C.; Stephen, M.; Roy, D.N.; Robertson, R.C.; Spencer, P.S.; Nunn, P.B.; Hugon, J.; et al. Guam Amyotrophic Lateral Sclerosis--Parkinsonism--Dementia Linked to a Plant Excitant Neurotoxin. Science 1987, 237, 517–522. [Google Scholar] [CrossRef] [PubMed]
  21. Murch, S.J.; Cox, P.A.; Banack, S.A. A Mechanism for Slow Release of Biomagnified Cyanobacterial Neurotoxins and Neurodegenerative Disease in Guam. Proc. Natl. Acad. Sci. USA 2004, 101, 12228–12231. [Google Scholar] [CrossRef] [PubMed]
  22. Cox, P.A.; Banack, S.A.; Murch, S.J.; Rasmussen, U.; Tien, G.; Bidigare, R.R.; Metcalf, J.S.; Morrison, L.F.; Codd, G.A.; Bergman, B. Diverse Taxa of Cyanobacteria Produce Beta-N-Methylamino-L-Alanine, a Neurotoxic Amino Acid. Proc. Natl. Acad. Sci. USA 2005, 102, 5074–5078. [Google Scholar] [CrossRef]
  23. Cox, P.A.; Sacks, O.W. Cycad Neurotoxins, Consumption of Flying Foxes, and ALS-PDC Disease in Guam. Neurology 2002, 58, 956–959. [Google Scholar] [CrossRef]
  24. Yan, B.Y.; Liu, Z.Q.; Huang, R.; Xu, Y.P.; Liu, D.M.; Wang, W.; Zhao, Z.W.; Cui, F.Y.; Shi, W.X. Impact Factors on the Production of Beta-Methylamino-L-Alanine (BMAA) by Cyanobacteria. Chemosphere 2020, 243, 125355. [Google Scholar] [CrossRef] [PubMed]
  25. Dunlop, R.A.; Banack, S.A.; Bishop, S.L.; Metcalf, J.S.; Murch, S.J.; Davis, D.A.; Stommel, E.W.; Karlsson, O.; Brittebo, E.B.; Chatziefthimiou, A.D.; et al. Is Exposure to BMAA a Risk Factor for Neurodegenerative Diseases? A Response to a Critical Review of the BMAA Hypothesis. Neurotox. Res. 2021, 39, 81–106. [Google Scholar] [CrossRef]
  26. Wang, Z.Q.; Wang, S.; Zhang, J.Y.; Lin, G.M.; Gan, N.; Song, L.; Zeng, X.; Zhang, C.C. Investigation on Cyanobacterial Production of the Proposed Neurotoxin β-N-Methylamino-L-Alanine (BMAA). Water Biol. Secur. 2023, 2, 100208. [Google Scholar] [CrossRef]
  27. Bishop, S.L.; Murch, S.J. A Systematic Review of Analytical Methods for the Detection and Quantification of β-: N -Methylamino-l-Alanine (BMAA). Analyst 2020, 145, 13–28. [Google Scholar] [CrossRef] [PubMed]
  28. Glover, W.B.; Mash, D.C.; Murch, S.J. The Natural Non-Protein Amino Acid N-β-Methylamino-L-Alanine (BMAA) Is Incorporated into Protein during Synthesis. Amino Acids 2014, 46, 2553–2559. [Google Scholar] [CrossRef] [PubMed]
  29. Dunlop, R.A.; Cox, P.A.; Banack, S.A.; Rodgers, K.J. The Non-Protein Amino Acid BMAA Is Misincorporated into Human Proteins in Place of L-Serine Causing Protein Misfolding and Aggregation. PLoS ONE 2013, 8, e75376. [Google Scholar] [CrossRef] [PubMed]
  30. Han, N.C.; Bullwinkle, T.J.; Loeb, K.F.; Faull, K.F.; Mohler, K.; Rinehart, J.; Ibba, M. The Mechanism of β-N-Methylamino-L-Alanine Inhibition of TRNA Aminoacylation and Its Impact on Misincorporation. J. Biol. Chem. 2020, 295, 1402–1410. [Google Scholar] [CrossRef]
  31. Korn, A.; Höfling, C.; Zeitschel, U.; Krueger, M.; Roßner, S.; Huster, D. Incorporation of the Nonproteinogenic Amino Acid β-Methylamino-Alanine Affects Amyloid β Fibril Properties and Toxicity. ACS Chem. Neurosci. 2020, 11, 1038–1047. [Google Scholar] [CrossRef] [PubMed]
  32. Liu, X.Q.; Rush, T.; Zapata, J.; Lobner, D. β-N-Methylamino-l-Alanine Induces Oxidative Stress and Glutamate Release through Action on System Xc-. Exp. Neurol. 2009, 217, 429–433. [Google Scholar] [CrossRef]
  33. Albano, R.; Lobner, D. Transport of BMAA into Neurons and Astrocytes by System Xc-. Neurotox. Res. 2018, 33, 1–5. [Google Scholar] [CrossRef]
  34. Schneider, T.; Simpson, C.; Desai, P.; Tucker, M.; Lobner, D. Neurotoxicity of Isomers of the Environmental Toxin L-BMAA. Toxicon 2020, 184, 175–179. [Google Scholar] [CrossRef] [PubMed]
  35. Lobner, D.; Piana, P.M.T.; Salous, A.K.; Peoples, R.W. β-N-Methylamino-l-Alanine Enhances Neurotoxicity through Multiple Mechanisms. Neurobiol. Dis. 2007, 25, 360–366. [Google Scholar] [CrossRef] [PubMed]
  36. Burton, B.; Collins, K.; Brooks, J.; Marx, K.; Renner, A.; Wilcox, K.; Moore, E.; Osowski, K.; Riley, J.; Rowe, J.; et al. The Biotoxin BMAA Promotes Dysfunction via Distinct Mechanisms in Neuroblastoma and Glioblastoma Cells. PLoS ONE 2023, 18, e0278793. [Google Scholar] [CrossRef] [PubMed]
  37. Martin, R.M.; Stallrich, J.; Bereman, M.S. Mixture Designs to Investigate Adverse Effects upon Co-Exposure to Environmental Cyanotoxins. Toxicology 2019, 421, 74–83. [Google Scholar] [CrossRef] [PubMed]
  38. Vega, A.; Bell, E.A.; Nunn, P.B. The Preparation of L- and d-α-Amino-β-Methylaminopropionic Acids and the Identification of the Compound Isolated from Cycas Circinalis as the l-Isomer. Phytochemistry 1968, 7, 1885–1887. [Google Scholar] [CrossRef]
  39. Davis, D.A.; Cox, P.A.; Banack, S.A.; Lecusay, P.D.; Garamszegi, S.P.; Hagan, M.J.; Powell, J.T.; Metcalf, J.S.; Palmour, R.M.; Beierschmitt, A.; et al. L-Serine Reduces Spinal Cord Pathology in a Vervet Model of Preclinical ALS/MND. J. Neuropathol. Exp. Neurol. 2020, 79, 393–406. [Google Scholar] [CrossRef] [PubMed]
  40. Cox, P.A.; Davis, D.A.; Mash, D.C.; Metcalf, J.S.; Banack, S.A. Dietary Exposure to an Environmental Toxin Triggers Neurofibrillary Tangles and Amyloid Deposits in the Brain. Proc. R. Soc. B Biol. Sci. 2016, 283, 20152397. [Google Scholar] [CrossRef]
  41. Nunn, P.B.; Codd, G.A. Metabolic Solutions to the Biosynthesis of Some Diaminomonocarboxylic Acids in Nature: Formation in Cyanobacteria of the Neurotoxins 3-N-Methyl-2,3-Diaminopropanoic Acid (BMAA) and 2,4-Diaminobutanoic Acid (2,4-DAB). Phytochemistry 2017, 144, 253–270. [Google Scholar] [CrossRef]
  42. Bishop, S.L.; Tymm, F.J.M.; Perry, K.; Kerkovius, J.K.; Menard, F.; Brady, A.; Slater, G.; Lim, D.S.S.; Metcalf, J.S.; Banack, S.A.; et al. Early-Earth Nonprotein Amino Acid Metabolites in Modern Cyanobacterial Microbialites. Environ. Chem. Lett. 2020, 18, 467–473. [Google Scholar] [CrossRef]
  43. Reveillon, D.; Sechet, V.; Hess, P.; Amzil, Z. Production of BMAA and DAB by Diatoms (Phaeodactylum tricornutum, Chaetoceros Sp., Chaetoceros calcitrans and, Thalassiosira pseudonana) and Bacteria Isolated from a Diatom Culture. Harmful Algae 2016, 58, 45–50. [Google Scholar] [CrossRef]
  44. Mantas, M.J.Q.; Nunn, P.B.; Codd, G.A.; Barker, D. Genomic Insights into the Biosynthesis and Physiology of the Cyanobacterial Neurotoxin 3-N-Methyl-2,3-Diaminopropanoic Acid (BMAA). Phytochemistry 2022, 200, 113198. [Google Scholar] [CrossRef] [PubMed]
  45. Pokorny, L.; Hausmann, B.; Pjevac, P.; Schagerl, M. How to Verify Non-Presence—The Challenge of Axenic Algae Cultivation. Cells 2022, 11, 2594. [Google Scholar] [CrossRef] [PubMed]
  46. Glover, W.B.; Baker, T.C.; Murch, S.J.; Bnown, P.N. Determination of Beta-N-Methylamino-L-Alanine, N-(2-Aminoethyl)Glycine, and 2,4-Diaminobutyric Acid in Food Products Containing Cyanobacteria by Ultra-Performance Liquid Chromatography and Tandem Mass Spectrometry: Single-Laboratory Validation. J. AOAC Int. 2015, 98, 1559–1565. [Google Scholar] [CrossRef]
  47. Koksharova, O.A.; Safronova, N.A. Non-Proteinogenic Amino Acid β-N-Methylamino-L-Alanine (BMAA): Bioactivity and Ecological Significance. Toxins 2022, 14, 539. [Google Scholar] [CrossRef]
  48. Lage, S.; Strom, L.; Godhe, A.; Rydberg, S. Kinetics of Beta-N-Methylamino-L-Alanine (BMAA) and 2, 4-Diaminobutyric Acid (DAB) Production by Diatoms: The Effect of Nitrogen. Eur. J. Phycol. 2019, 54, 115–125. [Google Scholar] [CrossRef]
  49. Kuzyk, S.B.; Ma, X.; Yurkov, V. Seasonal Dynamics of Lake Winnipeg’s Microbial Communities Reveal Aerobic Anoxygenic Phototrophic Populations Coincide with Sunlight Availability. Microorganisms 2022, 10, 1690. [Google Scholar] [CrossRef] [PubMed]
  50. Ishikawa, M.; Tanasupawat, S.; Nakajima, K.; Kanamori, H.; Ishizaki, S.; Kodama, K.; Okamoto-Kainuma, A.; Koizumi, Y.; Yamamoto, Y.; Yamasato, K. Alkalibacterium thalassium Sp. Nov., Alkalibacterium pelagium Sp. Nov., Alkalibacterium putridalgicola Sp. Nov. and Alkalibacterium kapii Sp. Nov., Slightly Halophilic and Alkaliphilic Marine Lactic Acid Bacteria Isolated from Marine Organisms and Salted Foods Collected in Japan and Thailand. Int. J. Syst. Evol. Microbiol. 2009, 59, 1215–1226. [Google Scholar] [CrossRef] [PubMed]
  51. Koerner, R.J.; Goodfellow, M.; Jones, A.L. The Genus Dietzia: A New Home for Some Known and Emerging Opportunist Pathogens. FEMS Immunol. Med. Microbiol. 2009, 55, 296–305. [Google Scholar] [CrossRef] [PubMed]
  52. Li, A.F.; Tian, Z.J.; Li, J.; Yu, R.C.; Banack, S.A.; Wang, Z.Y. Detection of the Neurotoxin BMAA within Cyanobacteria Isolated from Freshwater in China. Toxicon 2010, 55, 947–953. [Google Scholar] [CrossRef]
  53. Esterhuizen, M.; Downing, T.G. Beta-N-Methylamino-L-Alanine (BMAA) in Novel South African Cyanobacterial Isolates. Ecotoxicol. Environ. Saf. 2008, 71, 309–313. [Google Scholar] [CrossRef]
  54. Błaszczyk, A.; Siedlecka-Kroplewska, K.; Woźniak, M.; Mazur-Marzec, H. Presence of SS-N-Methylamino-L-Alanine in Cyanobacteria and Aquatic Organisms from Waters of Northern Poland; BMAA Toxicity Studies. Toxicon 2021, 194, 90–97. [Google Scholar] [CrossRef] [PubMed]
  55. Gross, M. The Mysteries of the Diatoms. Curr. Biol. 2012, 22, R581–R585. [Google Scholar] [CrossRef] [PubMed]
  56. Lage, S.; Burian, A.; Rasmussen, U.; Costa, P.R.; Annadotter, H.; Godhe, A.; Rydberg, S. BMAA Extraction of Cyanobacteria Samples: Which Method to Choose? Environ. Sci. Pollut. Res. 2016, 23, 338–350. [Google Scholar] [CrossRef] [PubMed]
  57. Jiang, L.; Eriksson, J.; Lage, S.; Jonasson, S.; Shams, S.; Mehine, M.; Ilag, L.L.; Rasmussen, U. Diatoms: A Novel Source for the Neurotoxin BMAA in Aquatic Environments. PLoS ONE 2014, 9, e84578. [Google Scholar] [CrossRef] [PubMed]
  58. Violi, J.P.; Facey, J.A.; Mitrovic, S.M.; Colville, A.; Rodgers, K.J. Production of β-Methylamino-L-Alanine (BMAA) and Its Isomers by Freshwater Diatoms. Toxins 2019, 11, 512. [Google Scholar] [CrossRef] [PubMed]
  59. Li, A.F.; Yan, Y.J.; Qiu, J.B.; Yan, G.W.; Zhao, P.; Li, M.; Ji, Y.; Wang, G.X.; Meng, F.P.; Li, Y.; et al. Putative Biosynthesis Mechanism of the Neurotoxin Beta-N-Methylamino-L-Alanine in Marine Diatoms Based on a Transcriptomics Approach. J. Hazard. Mater. 2023, 441, 129953. [Google Scholar] [CrossRef]
  60. Downing, S.; Banack, S.A.; Metcalf, J.S.; Cox, P.A.; Downing, T.G. Nitrogen Starvation of Cyanobacteria Results in the Production of Beta-N-Methylamino-L-Alanine. Toxicon 2011, 58, 187–194. [Google Scholar] [CrossRef]
  61. Esterhuizen-Londt, M.; Downing, S.; Downing, T.G. Improved Sensitivity Using Liquid Chromatography Mass Spectrometry (LC-MS) for Detection of Propyl Chloroformate Derivatised Beta-N-Methylamino-L-Alanine (BMAA) in Cyanobacteria. Water SA 2011, 37, 133–138. [Google Scholar] [CrossRef]
  62. Peters, S.J.; Rodgers, K.J.; Mitrovic, S.M.; Bishop, D.P. The Changes in Cyanobacterial Concentration of β-Methylamino-L-Alanine during a Bloom Event. Molecules 2022, 27, 7382. [Google Scholar] [CrossRef]
  63. Berntzon, L.; Erasmie, S.; Celepli, N.; Eriksson, J.; Rasmussen, U.; Bergman, B. BMAA Inhibits Nitrogen Fixation in the Cyanobacterium Nostoc Sp. PCC 7120. Mar. Drugs 2013, 11, 3091–3108. [Google Scholar] [CrossRef]
  64. Popova, A.A.; Rasmussen, U.; Semashko, T.A.; Govorun, V.M.; Koksharova, O.A. Stress Effects of Cyanotoxin Beta-Methylamino-L-Alanine (BMAA) on Cyanobacterial Heterocyst Formation and Functionality. Environ. Microbiol. Rep. 2018, 10, 369–377. [Google Scholar] [CrossRef]
  65. Popova, A.A.; Semashko, T.A.; Kostina, N.V.; Rasmussen, U.; Govorun, V.M.; Koksharova, O.A. The Cyanotoxin BMAA Induces Heterocyst Specific Gene Expression in Anabaena Sp. PCC 7120 under Repressive Conditions. Toxins 2018, 10, 478. [Google Scholar] [CrossRef]
  66. Jiao, Y.Y.; Chen, Q.K.; Chen, X.; Wang, X.; Liao, X.W.; Jiang, L.J.; Wu, J.; Yang, L.Y. Occurrence and Transfer of a Cyanobacterial Neurotoxin Beta-Methylamino-L-Alanine within the Aquatic Food Webs of Gonghu Bay (Lake Taihu, China) to Evaluate the Potential Human Health Risk. Sci. Total Environ. 2014, 468, 457–463. [Google Scholar] [CrossRef]
  67. Pip, E.; Munford, K.E.; Bowman, L. Seasonal Nearshore Occurrence of the Neurotoxin β-N-Methylamino-L-Alanine (BMAA) in Lake Winnipeg, Canada. Environ. Pollut. 2016, 5, 110. [Google Scholar] [CrossRef]
  68. Gutiérrez-García, K.; Whitaker, M.R.L.; Bustos-Díaz, E.D.; Salzman, S.; Ramos-Aboites, H.E.; Reitz, Z.L.; Pierce, N.E.; Cibrián-Jaramillo, A.; Barona-Gómez, F. Gut Microbiomes of Cycad-Feeding Insects Tolerant to β-Methylamino-L-Alanine (BMAA) Are Rich in Siderophore Biosynthesis. ISME Commun. 2023, 3, 122. [Google Scholar] [CrossRef]
  69. Kost, C.; Patil, K.R.; Friedman, J.; Garcia, S.L.; Ralser, M. Metabolic Exchanges Are Ubiquitous in Natural Microbial Communities. Nat. Microbiol. 2023, 8, 2244–2252. [Google Scholar] [CrossRef] [PubMed]
  70. Bishop, S.L.; Drikic, M.; Wacker, S.; Chen, Y.Y.; Kozyrskyj, A.L.; Lewis, I.A. Moving beyond Descriptive Studies: Harnessing Metabolomics to Elucidate the Molecular Mechanisms Underpinning Host-Microbiome Phenotypes. Mucosal Immunol. 2022, 15, 1071–1084. [Google Scholar] [CrossRef] [PubMed]
  71. Glover, W.B.; Murch, S.J. Analysis of B-Methylamino-Alanine in Environmental and Food Samples. In Phytochemicals-Biosynthesis, Function and Application; Jetter, R., Ed.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 149–174. [Google Scholar]
  72. EPA Definition and Procedure for the Determination of the Method Detection Limit, Revision 2; United States Environmental Protection Agency Office of Science and Technology: Washington, DC, USA, 2016.
  73. Komárek, J.; Katovský, J.; Mare, J.; Johansen, J.R. Taxonomic Classification of Cyanoprokaryotes (Cyanobacterial Genera) 2014, Using a Polyphasic Approach. Preslia 2014, 86, 295–335. [Google Scholar]
  74. Nienaber, M.A.; Steinitz-Kannan, M. Identification and Impact; University Press of Kentucky: Lexington, KY, USA, 2018; ISBN 9780813175591. [Google Scholar]
  75. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-Resolution Sample Inference from Illumina Amplicon Data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef]
  76. Lefler, F.W.; Berthold, D.E.; Dail Laughinghouse IV, H. CyanoSeq: A Database of Cyanobacterial 16S RRNA Sequences with Curated Taxonomy. J. Phycol. 2023, 59, 470–480. [Google Scholar] [CrossRef]
  77. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef] [PubMed]
  78. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [PubMed]
  79. Chamberlain, S.; Ooms, J. Classes for “GeoJSON”. Available online: https://github.com/ropensci/geojson, https://docs.ropensci.org/geojson/ (accessed on 18 April 2023).
  80. Agriculture and Agri Food Canada. Lake Winnipeg Outline; Canadian Watershed Information Network, V2; University of Manitoba Centre for Earth Observation Science: Winnipeg, Canada, 2017. [Google Scholar]
  81. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: Basel, Switzerland, 2016. [Google Scholar]
Figure 1. Concentration profiles and distribution of non-protein amino acids in Lake Winnipeg. (a) Lake Winnipeg (highlighted in the red box) is located in Manitoba, Canada. (b) N-(2-aminoethyl)glycine (AEG); (c) β-aminomethyl-L-alanine (BAMA); (d) β-N-methylamino-L-alanine (BMAA); and (e) 2,4-diaminobutyric acid (DAB) concentrations in crude (brown; n = 27 samples with 3 technical replicates) and isolate (blue; n = 16 samples with 4 technical replicates) cultures. Original field samples were collected between 2016 and 2021. Stations containing paired crude and isolate cultures from the same original field sample are shown on each map. Canada map outline adapted from the Atlas of Canada, Natural Resources Canada, and Lake Winnipeg outline adapted from the University of Manitoba Centre for Earth Observation Science.
Figure 1. Concentration profiles and distribution of non-protein amino acids in Lake Winnipeg. (a) Lake Winnipeg (highlighted in the red box) is located in Manitoba, Canada. (b) N-(2-aminoethyl)glycine (AEG); (c) β-aminomethyl-L-alanine (BAMA); (d) β-N-methylamino-L-alanine (BMAA); and (e) 2,4-diaminobutyric acid (DAB) concentrations in crude (brown; n = 27 samples with 3 technical replicates) and isolate (blue; n = 16 samples with 4 technical replicates) cultures. Original field samples were collected between 2016 and 2021. Stations containing paired crude and isolate cultures from the same original field sample are shown on each map. Canada map outline adapted from the Atlas of Canada, Natural Resources Canada, and Lake Winnipeg outline adapted from the University of Manitoba Centre for Earth Observation Science.
Toxins 16 00169 g001
Figure 2. Phylum breakdown of paired crude (brown) and isolate (blue) samples from the same parent field samples. Cultures were established from locations and years (a) Bloom 2019, (b) Station W9 2018, (c) Station 20 2018, and (d) Station 19 2018. Cultures shown in (a) and (b) were putatively identified as containing identical cyanobacterial species.
Figure 2. Phylum breakdown of paired crude (brown) and isolate (blue) samples from the same parent field samples. Cultures were established from locations and years (a) Bloom 2019, (b) Station W9 2018, (c) Station 20 2018, and (d) Station 19 2018. Cultures shown in (a) and (b) were putatively identified as containing identical cyanobacterial species.
Toxins 16 00169 g002
Figure 3. Comparison of non-protein amino acid concentrations in crude and isolate cultures. (a) N-(2-aminoethyl)glycine (AEG); (b) β-aminomethyl-L-alanine (BAMA); (c) β-N-methylamino-L-alanine (BMAA); and (d) 2,4-diaminobutyric acid (DAB) concentrations in all crude (brown; n = 27 samples with three technical replicates each, and isolate (blue; n = 16 samples with four technical replicates each) cultures and comparisons between paired crude and isolate cultures from the same original field sample. Error bars indicate standard error of the mean and not detected values were estimated as 1/10th of the method detection limit for the calculation of statistical differences via an unpaired t-test (α = 0.05) where * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001. Samples Bloom 2019 and Station W9 2018, highlighted in grey, were putatively identified as containing identical cyanobacterial species.
Figure 3. Comparison of non-protein amino acid concentrations in crude and isolate cultures. (a) N-(2-aminoethyl)glycine (AEG); (b) β-aminomethyl-L-alanine (BAMA); (c) β-N-methylamino-L-alanine (BMAA); and (d) 2,4-diaminobutyric acid (DAB) concentrations in all crude (brown; n = 27 samples with three technical replicates each, and isolate (blue; n = 16 samples with four technical replicates each) cultures and comparisons between paired crude and isolate cultures from the same original field sample. Error bars indicate standard error of the mean and not detected values were estimated as 1/10th of the method detection limit for the calculation of statistical differences via an unpaired t-test (α = 0.05) where * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001. Samples Bloom 2019 and Station W9 2018, highlighted in grey, were putatively identified as containing identical cyanobacterial species.
Toxins 16 00169 g003
Figure 4. Correlations between microbial diversity and non-protein amino acid (NPAA) production in cultures. (a) Alpha diversity measures (Shannon, Simpson) of crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site. Statistically significant differences were calculated by an unpaired t-test (α = 0.05). (b) Ordination (non-metric multi-dimensional scaling; NMDS) plot showing the beta diversity measure (Bray–Curtis distances) of crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site. (c) Principal component analysis biplot of auto-scaled NPAA concentrations from crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site and with alpha values (shading) showing the Shannon diversity scores. Undetected values were estimated as 1/10th of the method detection limit for statistical purposes.
Figure 4. Correlations between microbial diversity and non-protein amino acid (NPAA) production in cultures. (a) Alpha diversity measures (Shannon, Simpson) of crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site. Statistically significant differences were calculated by an unpaired t-test (α = 0.05). (b) Ordination (non-metric multi-dimensional scaling; NMDS) plot showing the beta diversity measure (Bray–Curtis distances) of crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site. (c) Principal component analysis biplot of auto-scaled NPAA concentrations from crude cultures (circles) compared to isolate cultures (triangles), colored according to the location of initial field collection site and with alpha values (shading) showing the Shannon diversity scores. Undetected values were estimated as 1/10th of the method detection limit for statistical purposes.
Toxins 16 00169 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bishop, S.L.; Solonenka, J.T.; Giebelhaus, R.T.; Bakker, D.T.R.; Li, I.T.S.; Murch, S.J. Microbial Diversity Impacts Non-Protein Amino Acid Production in Cyanobacterial Bloom Cultures Collected from Lake Winnipeg. Toxins 2024, 16, 169. https://doi.org/10.3390/toxins16040169

AMA Style

Bishop SL, Solonenka JT, Giebelhaus RT, Bakker DTR, Li ITS, Murch SJ. Microbial Diversity Impacts Non-Protein Amino Acid Production in Cyanobacterial Bloom Cultures Collected from Lake Winnipeg. Toxins. 2024; 16(4):169. https://doi.org/10.3390/toxins16040169

Chicago/Turabian Style

Bishop, Stephanie L., Julia T. Solonenka, Ryland T. Giebelhaus, David T. R. Bakker, Isaac T. S. Li, and Susan J. Murch. 2024. "Microbial Diversity Impacts Non-Protein Amino Acid Production in Cyanobacterial Bloom Cultures Collected from Lake Winnipeg" Toxins 16, no. 4: 169. https://doi.org/10.3390/toxins16040169

APA Style

Bishop, S. L., Solonenka, J. T., Giebelhaus, R. T., Bakker, D. T. R., Li, I. T. S., & Murch, S. J. (2024). Microbial Diversity Impacts Non-Protein Amino Acid Production in Cyanobacterial Bloom Cultures Collected from Lake Winnipeg. Toxins, 16(4), 169. https://doi.org/10.3390/toxins16040169

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