ORIGINAL RESEARCH article
Insight Into the Molecular Mechanisms for Microcystin Biodegradation in Lake Erie and Lake Taihu
- 1Department of Microbiology, The University of Tennessee, Knoxville, Knoxville, TN, United States
- 2Department of Biology, James Madison University, Harrisonburg, VA, United States
- 3Great Lakes Institute for Environmental Research, University of Windsor, Windsor, ON, Canada
- 4Department of Biological Sciences, Bowling Green State University, Bowling Green, OH, United States
- 5Department of Chemistry, College of Environmental Science and Forestry, State University of New York, Syracuse, NY, United States
Microcystins are potent hepatotoxins that are frequently detected in fresh water lakes plagued by toxic cyanobacteria. Microbial biodegradation has been referred to as the most important avenue for removal of microcystin from aquatic environments. The biochemical pathway most commonly associated with the degradation of microcystin is encoded by the mlrABCD (mlr) cassette. The ecological significance of this pathway remains unclear as no studies have examined the expression of these genes in natural environments. Six metatranscriptomes were generated from microcystin-producing Microcystis blooms and analyzed to assess the activity of this pathway in environmental samples. Seventy-eight samples were collected from Lake Erie, United States/Canada and Lake Tai (Taihu), China, and screened for the presence of mlr gene transcripts. Read mapping to the mlr cassette indicated transcripts for these genes were absent, with only 77 of the collective 3.7 billion reads mapping to any part of the mlr cassette. Analysis of the assembled metatranscriptomes supported this, with only distantly related sequences identified as mlrABC-like. These observations were made despite the presence of microcystin and over 500,000 reads mapping to the mcy cassette for microcystin production. Glutathione S-transferases and alkaline proteases have been previously hypothesized to be alternative pathways for microcystin biodegradation, and expression of these genes was detected across space and time in both lakes. While the activity of these alternative pathways needs to be experimentally confirmed, they may be individually or collectively more important than mlr genes in the natural environment. Importantly, the lack of mlr expression could indicate microcystin biodegradation was not occurring in the analyzed samples. This study raises interesting questions about the ubiquity, specificity and locality of microcystin biodegradation, and highlights the need for the characterization of relevant mechanisms in natural communities to understand the fate of microcystin in the environment and risk to public health.
Microcystins (MCs) are the one of the most frequently detected cyanotoxins within freshwater harmful cyanobacterial blooms (cyanoHABs). With cyanoHABS on the rise around the world, understanding the fate of MCs in the environment has assumed greater importance as they pose both ecological and public health risks. MCs can be potent protein phosphatase 1 and 2A inhibitors as well as potential tumor promoters. Consumption of MCs can result in acute hepatocytosis, cancer and various gastrointestinal problems (Carmichael, 2001; Zurawell et al., 2005; Carmichael and Boyer, 2016). Conventional municipal water treatments are effective at removing or inactivating cyanobacterial cells, but the methods risk the release of MCs from the cell. MCs are not removed by flocculation processes and require secondary treatments that include the use of activated carbon, which is effective but costly (He et al., 2016). Indeed, liver cancer and colorectal cancer has been linked to nearby bodies of water plagued with toxic cyanobacteria in the United States, Serbia and China (Carmichael, 2001; Hernández et al., 2009; Svircev et al., 2009). Cases of mortality and morbidity due to cyanotoxins have been reported since the 1800s involving birds, livestock, dogs, fish, and even humans (Chorus and Bartram, 1999; Wood, 2016). The number of cases has risen gradually over time and it is predicted that the number of individuals affected is greatly underestimated (Wood, 2016).
Microcystins are chemically stable and resistant to many abiotic factors, which has led to numerous studies questioning their fate in the environment (Harada and Tsuji, 1998). Whereas their disappearance has been attributed to some abiotic effects including dilution, adsorption, and photodegradation (Harada and Tsuji, 1998; Corbel et al., 2014), they are generally considered to be resistant to abiotic factors. MCs are also resistant to the activity of common peptidases (Tsuji et al., 1994; Mazur and Plinski, 2001; Smith et al., 2010) due to their cyclic structure and the alternating incorporation of the non-protein “R” stereoisomers of the amino acid, coupled with “iso” peptide bonds formed through the side chains of glutamate and aspartate. Despite these challenges, microorganisms that span multiple phyla and even domains are reported to degrade MCs (Supplementary Table S1). Further, experiments assessing the biodegradation of MCs which employ natural microbial communities from lakes (Christoffersen et al., 2002; Edwards et al., 2008; Eleuterio and Batista, 2010; Mou et al., 2013; Dziga et al., 2017; Lezcano et al., 2018), drinking water reservoirs (Cousins et al., 1996; Ho et al., 2012), estuary and sea water (Lemes et al., 2008), water treatment facilities (Lam et al., 1995; Saito et al., 2003b; Ma et al., 2014), soil (Miller and Fallowfield, 2001; Bibo et al., 2008; Cao et al., 2018; Redouane et al., 2019), lake sediments (Rapala et al., 1994; Cousins et al., 1996; Chen X. et al., 2010; Song et al., 2014; Li et al., 2016; Zhu et al., 2019) and biofilters (Grützmacher et al., 2002; Ho et al., 2006; Lee et al., 2006; Ho et al., 2007; Hoefel et al., 2009; Ho et al., 2012; Kumar et al., 2019) have demonstrated that biodegradation capacity is seemingly ubiquitous across a wide range of environments and occurs under a variety of conditions (Li et al., 2017). For these reasons, the current literature often refers to biodegradation as the most important route for the disappearance of MCs in nature (Christoffersen et al., 2002; Holst et al., 2003; Chen et al., 2008; Grützmacher et al., 2009; Bukowska et al., 2018).
The reported genetic pathway involved in MC biodegradation involves four enzymes encoded by mlrA, mlrB, mlrC, and mlrD (Bourne et al., 1996; Bourne et al., 2001). This pathway has been primarily associated with the biodegradation of MC-LR, a prevalent congener in nature (Carmichael and Boyer, 2016). MlrA, termed microcystinase, is a metalloprotease linearizing the MC-LR at the Arg-ADDA bond. It is commonly used as marker gene and considered to catalyze the most important step in biodegradation of MCs (Dziga et al., 2012; Dexter et al., 2018; Lezcano et al., 2018). MlrB is a serine protease that drives the hydrolysis of the linearized MC-LR at the Ala-Leu bond, forming a tetrapeptide. The tetrapeptide is then subject to further degradation by a second metalloprotease, MlrC, and it was demonstrated that MlrC can also cleave the linearized MC at the -Adda group (Dziga et al., 2012). MlrD is a putative transporter, but whether MC biodegradation occurs intracellularly or extracellularly remains inconclusive (Bourne et al., 2006).
This pathway has been extensively studied in culture (Supplementary Table S1, Bourne et al., 2006; Dziga et al., 2012; Maseda et al., 2012; Dziga et al., 2016; Dexter et al., 2018; as examples), although not all MC degrading isolates have tested positive for these genes. Further, the environmental relevance and distribution of the mlr cassette in the environment remains unclear. The mlrA gene has been detected in DNA isolated from microbial communities with MC degradation capabilities from freshwater samples (Mou et al., 2013; Kansole and Lin, 2016; Lezcano et al., 2018), but more commonly environments where MC biodegradation has been tested are not screened for use of this pathway. More often, the mlrA gene has been detected in enriched samples such as biofilms, constructed wetlands, filtration sand and water samples from treatment plants (Grützmacher et al., 2002; Bourne et al., 2006; Ho et al., 2006; Hoefel et al., 2009; Ho et al., 2012). Other genes in the cassette are not as widely studied in the environment and sequences of proteins with confirmed function are lacking, particularly for mlrC. By use of the same PCR primers (Saito et al., 2003a) and narrowly targeted degenerate qPCR primers (Hoefel et al., 2009) used amongst many studies, the understanding of the diversity of this pathway remains uncertain and makes a negative PCR result difficult to interpret. Further, since the methods used to screen for mlrA in native microbial communities have been based on DNA sequences (Hoefel et al., 2009; Mou et al., 2013; Zhu et al., 2014; Lezcano et al., 2018), only the potential of the pathway to be a route for MC biodegradation has thus far been considered. The activity of all the genes in the mlr pathway, either in the form of transcribed RNA or expressed proteins, has not been measured in natural samples.
Our goal in this study was to survey the natural environment for evidence of microcystin biodegradation during toxic blooms by examining the relative expression of the genes in the mlr pathway. We also screened for transcripts of other hypothesized mechanisms for MC biodegradation, including glutathione S-transferases (GSTs), CAAX proteases and alkaline proteases (Takenaka and Watanabe, 1997; Mou et al., 2013). Knowledge of MC biodegradation in the natural environment on a molecular level has potential use in predicting the capacity of the native microbial community to remove MCs from that environment and identifying environmentally relevant mechanisms that are responsible. This was done using shotgun metatranscriptomes from six surveys in eutrophic freshwater lakes with a combined total of 3.7 billion reads. The use of metatranscriptomes relative to DNA-based assays allows for a more accurate representation of the active function of the microbial community. Metatranscriptomes can also provide insight on the proportional representation of a gene’s activity within a community and provide the ability to detect more divergent sequences of the genes of interest.
Materials and Methods
Sample Collection and Sequencing
A total of six surveys, encompassing seventy-eight metatranscriptomes with a wide range of environmental conditions and MC concentrations, were used to screen for the expression of genes involved in MC biodegradation (Figure 1). Thirty-seven of the metatranscriptomes were generated from surface waters in the western basin of Lake Erie (United States/Canada) with sixteen from a bloom event in early August 2014 across eight stations (Steffen et al., 2017), seven tracking a bloom over a diurnal cycle in late August 2014 (Davenport et al., 2019), five from a bloom event in July 2013 from two stations (Steffen et al., 2017) and nine from three stations in August 2012 (Steffen et al., 2015). The Lake Erie dataset from early August of 2014 by Steffen et al. (2017) was of special interest because samples were taken at and around the City of Toledo’s water intake where dissolved MC concentrations in the source water were above 10 μg/L and exceeded 2.5 μg/L in the finished water, resulting in a 3 day water ban affecting 400,000 people (Carmichael and Boyer, 2016; Steffen et al., 2017).
Figure 1. Map of Lake Erie (A) and Taihu (B). Stations sampled for metatranscriptomes are depicted as green (Lake Erie) or blue (Taihu) circles. Coordinates for each sample location can be found in Supplementary Tables S3, S4. Particulate microcystin and chl a dynamics depict the range of bloom severity and toxicity amongst samples used in the metatranscriptomic analysis (C). The particulate microcystin and chl a concentrations for each sample can be found in Supplementary Tables S3, S4.
Forty-one metatranscriptomes were generated from surface water collected in the northern part of Lake Tai (Taihu), China in 2013 and 2014 during bloom events. The Taihu 2014 metatranscriptomes were composed of thirty-five samples collected across nine locations at monthly intervals from June to October, capturing an entire bloom event (Stough et al., 2017; Tang et al., 2018). Six metatranscriptomes were generated from samples collected in August 2013 from six stations. Sample collection, RNA extraction and sequencing for the Taihu 2013 metatranscriptomes were performed using the same approaches as previously described for published studies on Lake Erie (Steffen et al., 2015). The Taihu 2013 metatranscriptomes are publicly available on the MG-RAST (Meyer et al., 2008) server (ID: mgm4768721, mgm4768720, mgm4768724, mgm4768726, mgm4768728, mgm47687130) and the Lake Erie 2014 diel metatranscriptomes are available from the NCBI SRA database (SRP128942, SRP128954, SRP128945, SRP117911, SRP117914, SRA117922, SRA117915). More information on data availability, location, dates of sampling and environmental parameters for all metatranscriptomes can be found in Supplementary Tables S2–S4 or in the citations above.
Read Recruitment to the mlr Cassette
Expression of mlr pathway was first evaluated by recruiting reads from Lake Erie and Taihu metatranscriptomes to the mlrA-D from Sphingomonas sp. ACM3962 in CLC Genomics Workbench (Qiagen, Hilden, Germany). Reads were trimmed with a quality limit of 0.05 and an allowance of up to two ambiguous base pairs using CLC Genomics Workbench. Recruitments were performed using only reads >100 bp for all metatranscriptomes and a similarity fraction of 0.8 and a length fraction of 0.5 to encompass known and potential unknown diversity of these genes. Since the 2012 Lake Erie dataset generated reads of only 50 bp, slightly more stringent parameters were used to reduce non-specific mapping (length fraction 0.97 and similarity fraction 0.8). Reads were also recruited using the same metrics to mlrA-D from Sphingopyxis sp. strain C-1, Rhizobium sp. TH, and Novosphingobium sp. THN1. For comparisons, these reads were recruited to the mcy operon (mcyA-J) from the model species Microcystis aeruginosa NIES 843 with the more stringent similarity fraction and length fraction of 0.9 that still allowed for the inclusion of other toxin producing species (Steffen et al., 2017).
Screening for mlrABC, GSTs and Alkaline Proteases in Assembled Data
Reads were assembled from individual samples into contiguous sequences (“contigs”) in CLC Genomics Workbench using default parameters (bubble size = 50). All contigs >200 bp from the metatranscriptomes were screened using BLASTx against a locally curated database containing the currently available MlrA, MlrB, and MlrC amino acid sequences in NCBI’s GenBank. Contigs with positive hits to MlrA, MlrB, and MlrC (contigs with an E-value of <0.00001) were collected and translated to amino acid sequences using the appropriate reading frame identified by BLASTx. To compare active sites, sequences for MlrA, MlrB, and MlrC were aligned to their respective candidate contigs in CLC Genomics Workbench. This method of mlr detection was also applied to a metagenome from a study on Lake Erie where MC biodegradation and mlrA-like genes have previously been reported (Mou et al., 2013) as well as a metatranscriptome generated from samples from a winter diatom bloom (Edgar et al., 2016). All contigs detected using these methods can be found in the Supplementary Data Sheet S1.
Due to the phylogenetic diversity of GSTs and alkaline proteases, these sequences were identified in the metatranscriptomes by first uploading all assembled data to MGRAST. Contigs were annotated using the best representative hit and an E-value of <0.00001 against the SEED database for functional classification and Refseq for taxonomic assignment. It should be noted that no contigs were annotated as mlr genes. Contigs were first filtered by phyla, only including those from the phyla Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes. These phyla were chosen because in most of the metatranscriptomes, these were classified as the most highly active non-cyanobacterial prokaryotic transcripts. In addition, all currently known MC degrading bacteria are within these phylogenetic groups. Contigs or portions of contigs annotated as GSTs (E.C. 22.214.171.124) filtered by representative hit were downloaded and reads were mapped with high stringency (0.97 similarity fraction and length fraction) and normalized by contig length and library size to determine expression of these genes in each sample. This stringent parameter was used so as to only recruit very specific reads to each contig. The only annotated alkaline protease contigs present in any metatranscriptomes from the phyla of interest were annotated as secreted alkaline metalloproteinase PrtA/B/C/G (E.C. 3.3.24.-). Expression of these genes were determined as previously described for GSTs. Shadeplots were constructed in GraphPad Prism v7.03 after being log10 transformed for visualization, and the shading of each box represents levels of total expression in a sample from a particular gene; the darker the shading (the more positive the number) indicates higher expression. All assembled data can be found on MGRAST and sample IDs are available in Supplementary Tables S3, S4.
Phylogenetic Analysis of MlrA, MlrB and MlrC Candidate Contigs
Amino acid reference sequences for CAAX protease and bacteriocin processing (CPBP) intramembrane metalloprotease family proteins were downloaded from UNIPROT as well as the top BLASTx results from the NCBI Refseq database for each mlrA-like contig were used to build a reference tree for phylogenetic analysis of MlrA candidate contigs. The 100 most closely related sequences in the Refseq database to MlrB and MlrC according to BLASTx were used as reference sequences for MlrB and MlrC candidate contigs. Bacterial beta-lactamase amino acid sequences from UNIPROT with confirmed function were also incorporated into MlrB reference tree. Alignments of reference sequences were performed using MUSCLE in Mega v7.0 with eight iterations (Tamura et al., 2007). Maximum likelihood trees were made using the PhyML server according to the LG model, and likelihood ratios and branch support were calculated using a Shimodaira-Hasegawa (SH)-like approximate likelihood ratio test (aLRT-SH-like, Guindon et al., 2010). Contigs were incorporated into the tree using pplacer1 as in Krausfeldt et al. (2017).
Measured particulate concentrations of MC verified that the majority of the samples were toxic and in conjunction with chl a concentrations indicated these samples ranged in bloom severity and toxicity levels (Figure 1C and Supplementary Tables S3, S4). Measurements of dissolved MC, particulate MC and chl a from 2015–2017 in Lake Erie indicated the presence of particulate MC and chl a corresponded to the presence of dissolved MC throughout blooms events (Supplementary Figure S1).
Read Mapping to the mlr Cassette
Reads from each metatranscriptome were first recruited to all four genes in the mlr pathway from Sphingomonas sp. ACM3962 with parameters set to account for diversity of known mlr sequences (Supplementary Figure S2). Only 77 reads out of the total of over 3.5 billion reads surveyed mapped to genes in the mlr pathway; 50 from Lake Erie, and 27 from Taihu (Figures 2A,B). Notably, mapped reads did not provide full coverage of any of the individual genes in the mlrABCD cassette. Reads were also mapped to the mlr cassette identified in other MC degraders, and similarly, few reads mapped with no full coverage of any gene in the cassette (Supplementary Table S5). Using the same parameters for the mcy operon, over 3,000,000 reads mapped to the mcy genes. When using more stringent parameters, which increased the similarity necessary for mapping, almost 500,000 reads mapped to the mcy operon, the genes for MC production (Figures 2A,B).
Figure 2. Number of reads mapped on a log scale from metatranscriptomes from Lake Erie (A) and Taihu (B) to the mlr cassette (red bars) for Sphingomonas ACM3962 for microcystin biodegradation (at a similarity fraction of 0.8 and length fraction of 0.5) and mcy cassette (black bars) from Microcystis aeruginosa NIES843 for microcystin production (at a similarity and length fraction of 0.9).
Screening for mlrA, mlrB, and mlrC
Reads from each metatranscriptome were assembled into contigs and screened for the marker gene mlrA using BLASTx. Only ten contigs were identified by this analysis to be mlrA-like. After screening these contigs against the Refseq database using BLASTx, four contigs were annotated as part of the same protein family (CPBP) as mlrA (Pei et al., 2011) due to the presence of an abi protein domain, and these were considered potential mlrA candidate transcripts. Three contigs were from the 2014 diel study on Lake Erie and one was from Taihu in 2013. BLASTx indicated all sequences shared ∼37 to ∼57% identity and ∼54 to ∼84% positives with at least one MlrA sequence with known function, but the E-values were all >1 × 10–14 (Supplementary Data Sheet S1). Further phylogenetic analysis using translated protein sequences was performed, and all of these contigs were very distantly related to sequences within the MlrA clade (Figure 3). Over 200 contigs were identified as mlrB candidates while only nine were identified as mlrC candidates (Supplementary Figure S3). Phylogenetic analysis indicated there were seven clades of MlrB-like contigs, and these contigs were detected in both lakes. Both MlrB and MlrC candidate sequences were distantly related to the tightly clustered clades of MlrB and MlrC sequences with confirmed function. Alignments of reference MlrA and MlrC sequences with respective contigs indicated the MlrA and MlrC candidate transcripts did not contain the appropriate active site residues needed for functionality against the MC molecule or target degradation products (Supplementary Figures S4, S5), although some of the contigs were not long enough to be aligned to that region. Many MlrB candidate sequences did contain the residues within the active site needed to function properly against linearized MC (Supplementary Figure S6).
Figure 3. Maximum likelihood phylogenetic tree with candidate mlrA sequences aligned using pplacer indicated by large closed circles (light blue = Lake Erie Diel 2014; green = Lake Taihu 2013). The letters correspond to which sample they were found: A: Taihu 2013 at the dock; B: Lake Erie (Diel) 2014 1600 h; C: Lake Erie (Diel) 2014 2200 h; D: Lake Erie (Diel) 2014 2200 h. The mlrA clade is collapsed for visualization purposes in a closed red wedge and other collapsed branches are denoted by open wedges. Smaller closed circles at nodes represent likelihood ratios. ∙∙∙ > 90, ∙∙ > 70, ∙ > 50.
For comparison, this technique was also used on 454-generated metagenomic data from Mou et al. (2013) from mesocosm experiments using water from Lake Erie to screen for mlrA. Four metagenomes were generated: two that followed treatment with MCs and two which were controls to examine changes in community structure and function upon MC additions. One 500 bp sequence from one MC treatment shared 95.63% identity with mlrA from Sphingopyxis sp. C1 (a sequence used for screening in the analysis above) with an E-value of 0, demonstrating the ability of the technique used here to identify genes of interest. Other sequences that shared identity with mlrA had an E-value > 0.00001, indicating there was a very low probability that these shared the same function as mlrA. Analysis of a metatranscriptome from a non-cyanobacterial bloom event in the winter yielded no mlrA or mlrA-like sequences.
Determining Expression of mlrA-Like Genes and Genes in Alternative Pathways
To determine if the putative CPBP genes (mlrA-like genes), identified by BLASTx were potential candidates for MC removal, reads from the samples in which the contigs were found were mapped to these respective contigs. Reads mapped accounted for a low percentage of the total number of reads in each individual sample, with only 72 reads mapping in total (Figure 4 and Supplementary Tables S3, S4). Annotations assigned from the SEED database were used to try to identify more contigs that were transcripts for proteins with the CPBP family. All of them were annotated as abortive infection proteins and of cyanobacterial origin, so they were excluded from further analysis.
Figure 4. Shadeplots representing a qualitative comparison of reads that mapped to candidate mlrA sequences (mlrA-like), GSTs and alkaline protease (AP) contigs from non-cyanobacterial prokaryotes, and mcyD from M. aeruginosa NIES843, across all samples collected from Lake Erie (A) and Lake Taihu (B). The reads that mapped to any putative genes of interest were totaled within each sample and normalized by contig length and total library size. Shading in each box represents total expression levels for each gene; the darker the color (more positive the number) the higher the expression. For visualization purposes, these numbers were log transformed. A dashed line within a box depicts a value of 0. Library sizes and values used for the shadeplot prior to transformation can be found in Supplementary Tables S3, S4. Samples in the shadeplot also appear in the same order as they do in the tables for reference.
Contigs were identified across all studies in Lake Erie and Taihu as putative GSTs and putative alkaline proteases, with the exception of Lake Erie in 2012 (Figure 4). Reads from samples with putative GSTs or alkaline proteases were mapped to the respective contigs to compare relative abundances of reads from these genes to mlrA-like genes (Figure 4 and Supplementary Tables S3, S4). Expression of mlrA-like genes was generally lower than GSTs or alkaline proteases in samples where mlrA-like genes were detectable. Trends in temporal and spatial distribution of GST and alkaline protease transcripts were more similar to those observed for mcyD transcripts.
With the prevalence and increased severity of cyanoHABs worldwide, it is important to understand the fate of the toxins produced during these blooms. The capacity for MC biodegradation has been demonstrated in microbial communities from a variety of environments all over the world suggesting biodegradation of MCs is ubiquitous. The mlr pathway is currently the only characterized mechanism for MC biodegradation, however, there are many gaps regarding the diversity and ubiquity of this pathway especially in the natural environment. Many studies only use PCR approaches to determine if mlr is present in either isolates, sand filters or environmental samples (as examples, Chen J. et al., 2010; Yang et al., 2014a; Lezcano et al., 2016; Dziga et al., 2017; Thees et al., 2018). One limitation to PCR-based approaches for detecting mlr genes is that more diverse sequences may be missed due to mismatch of primers. Currently, all of the publicly available mlrA sequences share 87–100 percent identity on both the nucleotide and amino acid level (Supplementary Figure S2). Repeated use of the primers based on highly identical sequences may explain why there are still few representative mlrA sequences available outside of Alphaproteobacterial lineage and consequently, could explain why many MC degraders remain “negative” for a microcystinase homolog. Even degenerate primers that can capture more diverse sequences can be difficult to construct when reference sequences are highly identical. Genome sequencing of MC degraders, shotgun metagenomes and metatranscriptomes provide the ability to expand the knowledge about the diversity of genes not yet discovered, but few studies have employed these mechanisms to study MC biodegradation specifically (Mou et al., 2013; Maghsoudi et al., 2016; Wang et al., 2018; Zhang et al., 2018).
In this study, 78 metatranscriptomes encompassing both daily and seasonal bloom dynamics across several years were analyzed to investigate mechanisms of MC biodegradation during toxic cyanobacterial blooms. Despite the notion that MC biodegradation is ubiquitous across different environmental samples worldwide and the presence of active toxin-producing cyanobacterial species in each sample from Lake Erie and Taihu, there was little evidence to suggest the mlr pathway was an active mechanism for MC removal from these lakes. Very few reads mapped to the mlr cassette from four MC degrading bacteria and read mapping did not provide full coverage across any of the genes to indicate transcripts were present. These results suggested that none of the genes in the mlr pathway were expressed concomitantly or alone in these samples. This was supported by phylogenetic analysis of MlrA candidate sequences identified by BLASTx: all were distantly related to MlrA suggesting transcripts for this marker gene for MC biodegradation was absent. However, these sequences were part of the same protein family as MlrA, the CPBP family. While the function of these proteins in prokaryotes remains largely unclear, they are typically annotated as CAAX proteases or abortive infection proteins, and the literature suggests they may confer immunity to antimicrobial peptides secreted by the self or other bacteria (Pei et al., 2011; Barrett et al., 2012). A role for these proteins in MC biodegradation has been hypothesized previously (Kansole and Lin, 2016), but this was based on presence/absence and remains to be experimentally confirmed. The analyses presented here suggested CAAX protease transcripts were not prevalent across or within samples, and the transcripts detectable were not highly expressed and did not contain the appropriate active site for functionality (Dziga et al., 2012).
In agreement with the read recruitments, BLASTx screening confirmed that MlrB and MlrC were likely not being expressed in any of the metatranscriptomes either. All MlrB and MlrC candidate transcripts were distantly related to MlrB and MlrC sequences with known function against MCs. Similarly to MlrA candidates, few MlrC candidate sequences were identified. The diversity and prevalence of MlrB candidate sequences is interesting as MlrB has been called a silent member of the mlr gene cluster in some MC degraders (Jiang et al., 2011), and yet also a better biomarker for MC biodegradation than mlrA (Tsao et al., 2017). Many of the MlrB candidates contained the residues within the active site required to be functional, but this active site is conserved within the Beta-Lactamase protein family, and may not be a good predictor of function against linearized MC. Therefore, phylogenetic analyses alone were not sufficient to conclude these were indeed distantly related proteins with the same function. Still, this may have suggested there more diverse mlr genes that would not be captured by currently published primer sets or that there are more promiscuous enzymes in the environment with the ability to degrade MC. Ultimately these data suggest that transcripts from the mlr pathway were absent.
Microcystin biodegradation has been observed by the natural microbial community residing in Lake Erie (Mou et al., 2013) and Taihu (Chen et al., 2008; Li et al., 2016) previously, and several bacterial strains reported to degrade MCs have been isolated from both lakes as well (Hu et al., 2009; Chen J. et al., 2010; Jiang et al., 2011; Yang et al., 2014a, b; Zhu et al., 2016; Krishnan et al., 2018; Thees et al., 2018). The potential for MC biodegradation via the mlr pathway was also evident in both lakes; several mlrA+ MC degraders have been isolated from samples collected from Taihu (Supplementary Table S1), and the present study confirmed mlrA was detectable in the metagenomes from Lake Erie. Interestingly, only one mlrA+ isolate has been collected from the water column in Taihu (Jiang et al., 2011) while several others were isolated from sludge and sediment (Yang et al., 2018a, b; Supplementary Table S1). Indeed, MC biodegradation by the mlr pathway might be more prevalent in lake sediments (Wu et al., 2015; Bukowska et al., 2018) which were not specifically examined here. This is potentially supported by previous observations in Lake Erie: MC degrading bacteria cultured were all mlrA– and isolated from surface water or a visible cyanobacterial bloom. It is possible that mlr harboring microorganisms more commonly occupy niches besides the water column, like the sediments, in which case the timing and location of MC biodegradation by the mlr pathway in the water column will be highly dependent on lake geomorphology and climatic events. This may explain why no mlr transcripts were detectable in the metatranscriptomes in the current study despite the presence of MC; they were all generated from water or scum samples collected at the surface.
Assuming the genetic potential for the mlr pathway did exist in the water column, the absence of an active mlr pathway could to due to a number of factors. It is possible mlr+ microbes existed in low abundance during the blooms due to competition. This was proposed in Lezcano et al. (2018) where increased abundances of mlrA gene copies were observed during bloom decline in the San Juan Reservoir in Spain, and the work by Zhu et al. (2014) on Lake Dianchi in China generally supported these observations. However, since the general activity (by proxy of transcripts) of other non-cyanobacterial prokaryotes is readily detectable even during dense blooms (Steffen et al., 2012; Chen et al., 2018), this suggested that competition with cyanobacteria was not solely responsible for the low abundance of mlr+ microorganisms or absence of mlr gene expression. Although the sample collection encompassed a wide range of environmental conditions (temperature, pH, nutrient concentrations, and oxygen levels), rates of MC biodegradation and the effect of environmental factors vary across different habitats (Li et al., 2017). Activity of the mlr pathway could have required specific conditions not captured in this study. However, rates of biodegradation are usually higher in lakes with bloom history (Li et al., 2017), which may mean transcripts were actually fleeting, and they were missed in this analysis. It is also possible these genes were transcribed at such a low level that they were not detectable or that they were simply missed during the sampling time despite the spatial and temporal sampling coverage of areas that experience blooms in these lakes.
Microcystin does induce mlr gene expression in culture (Shimizu et al., 2011), and many studies have enriched for MC degraders by one or more additions of MCs (i.e., Krishnan et al., 2018; Thees et al., 2018), suggesting activity of the mlr pathway could have been substrate limited. It is reasonable to assume that MC concentrations in the water column would be higher during bloom collapse, as observed in Lezcano et al. (2018), since MC remains intracellular (in the particulate phase) until released after cell lysis (into the dissolved phase). This could increase the abundance of mlr+ microbes and expression or mlr genes. However, dissolved MCs are not only present during bloom collapse; particulate MC in Lake Erie corresponded to the presence of dissolved MC (Supplementary Figure S1), and several studies that have detected dissolved MC in the water column during blooms in Taihu (Song et al., 2007; Su et al., 2018). This is likely due to cell turnover or lytic phage infection, which is evident in many of the samples analyzed here (Steffen et al., 2015; Steffen et al., 2017; Stough et al., 2017) and corresponded with very high particulate MC levels (several samples reaching >100 μg/L in Taihu). This highlights that MC was not only an available substrate at the very end of a bloom season in Lake Erie and Taihu. If the activity of the mlr pathway is dependent on higher concentrations of MC, it may not be a dominant pathway for MC biodegradation under the most environmentally relevant conditions in these lakes. Indeed, meta-analyses capturing “non-bloom” events or community responses before, during and after a bloom collapse would be useful to address the potential seasonal relevance of the mlr pathway in Lake Erie and Taihu. Currently, there are very few or no “late bloom or no bloom” metatranscriptomes to date, nor are there any metatranscriptomes generated from samples where MC biodegradation activity has been confirmed. These types community analysis will be essential in comprehensively determining the relevance of the mlr pathway in Lake Erie and Taihu.
Alternatively, the existence of other pathways for MC biodegradation has been suggested numerous times in the literature when mlr genes have not been detectable. No other mechanisms have been fully elucidated to date, but a few have been hypothesized. GSTs have been implicated in the detoxification of MCs in some higher plants, invertebrates and vertebrates by first conjugating to the MCs to be targeted for removal by further processes (Pflugmacher et al., 1998; Campos and Vasconcelos, 2010; Antas et al., 2018). While the role of GSTs in bacterial degradation of MCs has not been confirmed, they can play a role in biodegradation of other toxic compounds (Allocati et al., 2009). Further, Mou et al. (2013) observed higher abundances of GSTs where MCs were degraded in Lake Erie as well as other genes involved in xenobiotic metabolism. Takenaka and Watanabe (1997) provided experimental evidence for the role of alkaline proteases from Pseudomonas aeruginosa in the biodegradation of MCs, though the specific alkaline protease was not well described. With the exception of Lake Erie 2012, GSTs and alkaline proteases were highly expressed in samples collected from Lake Erie and Taihu compared to mlrA-like sequences in this study. These genes are inherently multifunctional and linking them to this process would require experimental confirmation. However, if their activity against MC is confirmed, the prevalent expression of these genes across time and space would suggest these pathways may be a more important mechanism of MC removal than the mlr pathway. The samples collected for the Lake Erie 2012 study were collected on 20 μm mesh and may indicate microbes expressing these genes are more prevalent in the free-living community rather than within cyanobacterial colonies or large particles. Transcripts for serine proteases (E.C. 3.4.21-) were sporadically observed as well (data not shown), demonstrated to be relevant in MC degrading Lactobacilli (Nybom et al., 2012). This is a broad categorization for enzymes that include common peptidases like trypsin, proteinase K, chymotrypsin and pepsin, and previous studies have shown that these enzymes in isolation can degrade MCs (Kankaanpää et al., 2005). Others have noted this would be a slow process compared to other peptides and proteins (Smith et al., 2010). Collectively, the observations in this study raise interesting questions about the specificity, diversity and ubiquity of MC biodegradation pathways in the environment.
The characterization of both mlr and non-mlr mechanisms for MC biodegradation in the natural environment is needed, as well as further investigation into the locality, abundance and activity of native MC degrading microorganisms during blooms to understand of the fate of cyanotoxins like MCs and the risk to human health. In the present study, an extensive and quantitative survey of six studies yielding 78 metatranscriptomes from two eutrophic freshwater lakes suggested that mlr genes were not being expressed in either of these systems. The samples were taken during MC-producing bloom events where cyanobacteria were actively transcribing mcy genes for MC production and MC was measurable throughout, yet mlr transcripts were not detectable. Despite the detection of a small number of potential mlrA-like genes sequences, these do not appear to be highly expressed or widely distributed across samples examined in this multi-survey study. Primarily, these samples were collected from Microcystis blooms, but samples collected at station 1163 in Sandusky Bay of Lake Erie suggest this pathway was not active during Planktothrix blooms either (Hampel et al., 2019). This is in accordance with previous work, especially from Lake Erie, and supports the on-going hypothesis that other pathways are involved in this process. While there is a need for experimental confirmation of activity against MCs, expression of GSTs and evidence for other protease activity warrants further investigation of multiple pathways for MC biodegradation in environmental samples. The complete lack of mlr detection suggested either that other pathways for MC biodegradation are more active during blooms or there are specific niches where this process occurs that may not include the water column. Importantly, this could also indicate that MC biodegradation was not occurring. If the latter is true, this highlights the crucial need for routine detection and implementation of appropriate water treatments as biodegradation of MCs may not be as universally important as previously thought.
Data Availability Statement
LK, GLB, and SW designed the study. MS extracted RNA from Taihu in 2013. GSB and RM provided the data from the Lake Erie diel study. LK analyzed the metatranscriptomes. All authors drafted the manuscript.
This work was made possible due to grant NA11NOS4780021 from NOS National Center for Coastal Ocean Science (NCCOS) under their Harmful Algal Bloom Prevention Control and Mitigation program to SW and GB, and National Science Foundation awards (DEB-1240870; IOS-1451528) to SW. This work was partially supported by funding from the NIH (1P01ES028939-01) and NSF (OCE-1840715) to the Bowling Green State University Lake Erie Center for Fresh Waters and Human Health. Sequencing work conducted by the United States Department of Energy Joint Genome Institute (RMLM), a Department of Energy Office of Science User Facility, was supported by the Office of Science of the United States Department of Energy under Contract No. DE-AC02-05CH11231.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Dr. Xiaozhen Mou for providing access to metagenomes. We also acknowledge the generous support of the Kenneth and Blaire Mossman endowment to the The University of Tennessee, Knoxville.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02741/full#supplementary-material
Antas, P., Carneiro, M., Reis, B., Castelo-Branco, R., Azevedo, J., Urbatzka, R., et al. (2018). GST transcriptional changes induced by a toxic Microcystis aeruginosa strain in two bivalve species during exposure and recovery phases. Ecotoxicology 27, 1272–1280. doi: 10.1007/s10646-018-1980-y
Bibo, L., Yan, G., Bangding, X., Jiantong, L., and Yongding, L. (2008). A laboratory study on risk assessment of microcystin-RR in cropland. J. Environ. Manag. 86, 566–574. doi: 10.1016/j.jenvman.2006.12.040
Bourne, D. G., Blakeley, R. L., Riddles, P., and Jones, G. J. (2006). Biodegradation of the cyanobacterial toxin microcystin LR in natural water and biologically active slow sand filters. Water Res. 40, 1294–1302. doi: 10.1016/j.watres.2006.01.022
Bourne, D. G., Jones, G. J., Blakeley, R. L., Jones, A., Negri, A. P., and Riddles, P. (1996). Enzymatic pathway for the bacterial degradation of the cyanobacterial cyclic peptide toxin microcystin LR. Appl. Envion. Micro. 62, 4086–4094.
Bourne, D. G., Riddles, P., Jones, G. J., Smith, W., and Blakeley, R. L. (2001). Characterisation of a gene cluster involved in bacterial degradation of the cyanobacterial toxin microcystin LR. Environ. Toxicol. 16, 523–534. doi: 10.1002/tox.10013
Bukowska, A., Kaliński, T., and Chróst, R. J. (2018). Degradation of microcystins by water and bottom sediment bacterial communities from a eutrophic freshwater lake. Aquat. Microbl. Ecol. 82, 129–144. doi: 10.3354/ame01887
Cao, Q., Steinman, A. D., Yao, L., and Xie, L. (2018). Effects of light, microorganisms, farming chemicals and water content on the degradation of microcystin-LR in agricultural soils. Ecotoxicol. Environ. Safa 156, 141–147. doi: 10.1016/j.ecoenv.2018.03.030
Carmichael, W. W., and Boyer, G. L. (2016). Health impacts from cyanobacteria harmful algae blooms: Implications for the North American Great Lakes. Harmful Algae. 54, 194–212. doi: 10.1016/j.hal.2016.02.002
Chen, J., Hu, L. B., Zhou, W., Yan, S. H., Yang, J. D., Xue, Y. F., et al. (2010). Degradation of microcystin-LR and RR by a Stenotrophomonas sp. strain EMS isolated from Lake Taihu, China. Int. J. Mol Sci. 11, 896–911. doi: 10.3390/ijms11030896
Chen, X., Yang, X., Yang, L., Xiao, B., Wu, X., Wang, J., et al. (2010). An effective pathway for the removal of microcystin LR via anoxic biodegradation in lake sediments. Water Res. 44, 1884–1892. doi: 10.1016/j.watres.2009.11.025
Chen, W., Song, L., Peng, L., Wan, N., Zhang, X., and Gan, N. (2008). Reduction in microcystin concentrations in large and shallow lakes: water and sediment-interface contributions. Water Res. 42, 763–773. doi: 10.1016/j.watres.2007.08.007
Chen, Z., Zhang, J., Li, R., Tian, F., Shen, Y., Xie, X., et al. (2018). Metatranscriptomics analysis of cyanobacterial aggregates during cyanobacterial bloom period in Lake Taihu, China. Environ. Sci. Poll. Res. 25, 4811–4825. doi: 10.1007/s11356-017-0733-4
Christoffersen, K., Lyck, S., and Winding, A. (2002). Microbial activity and bacterial community structure during degradation of microcystins. Aquat. Microb. Ecolog. 27, 125–136. doi: 10.1016/j.envpol.2017.08.013
Corbel, S., Mougin, C., and Bouaïcha, N. (2014). Cyanobacterial toxins: modes of actions, fate in aquatic and soil ecosystems, phytotoxicity and bioaccumulation in agricultural crops. Chemosphere 96, 1–15. doi: 10.1016/j.chemosphere.2013.07.056
Davenport, E. J., Neudeck, M. J., Matson, P. G., Bullerjahn, G. S., Davis, T. W., and Wilhelm, S. W. (2019). Metatranscriptomic analyses of diel metabolic functions during a Microcystis bloom in western Lake Erie (USA). Front. Microbiol. 10:2081. doi: 10.3389/fmicb.2019.02081
Dexter, J., Dziga, D., Lv, J., Zhu, J., Strzalka, W., Maksylewicz, A., et al. (2018). Heterologous expression of mlrA in a photoautotrophic host–engineering cyanobacteria to degrade microcystins. Environ. Pol. 237, 926–935. doi: 10.1016/j.envpol.2018.01.071
Dziga, D., Maksylewicz, A., Maroszek, M., Budzyńska, A., Napiorkowska-Krzebietke, A., Toporowska, M., et al. (2017). The biodegradation of microcystins in temperate freshwater bodies with previous cyanobacterial history. Ecotoxicol. Environ. Saf. 145, 420–430. doi: 10.1016/j.ecoenv.2017.07.046
Dziga, D., Wladyka, B., Zielińska, G., Meriluoto, J., and Wasylewski, M. (2012). Heterologous expression and characterisation of microcystinase. Toxicon 59, 578–586. doi: 10.1016/j.toxicon.2012.01.001
Dziga, D., Zielinska, G., Wladyka, B., Bochenska, O., Maksylewicz, A., Strzalka, W., et al. (2016). Characterization of enzymatic activity of MlrB and MlrC proteins involved in bacterial degradation of cyanotoxins microcystins. Toxins 8:76. doi: 10.3390/toxins8030076
Edgar, R., Morris, P., Rozmarynowycz, M., D’souza, N., Moniruzzaman, M., Bourbonniere, R., et al. (2016). Adaptations to photoautotrophy associated with seasonal ice cover in a large lake revealed by metatranscriptome analysis of a winter diatom bloom. J. Great Lakes Res. 42, 1007–1015. doi: 10.1016/j.jglr.2016.07.025
Eleuterio, L., and Batista, J. R. (2010). Biodegradation studies and sequencing of microcystin-LR degrading bacteria isolated from a drinking water biofilter and a fresh water lake. Toxicon 55, 1434–1442. doi: 10.1016/j.toxicon.2010.02.020
Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010
Hampel, J. J., McCarthy, M. J., Neudeck, M., Bullerjahn, G. S., McKay, R. M. L., and Newell, S. E. (2019). Ammonium recycling supports toxic Planktothrix blooms in sandusky bay, lake erie: evidence from stable isotope and metatranscriptome data. Harmful Algae. 81, 42–52. doi: 10.1016/j.hal.2018.11.011
Harada, K.-I., and Tsuji, K. (1998). Persistence and decomposition of hepatotoxic microcystins produced by cyanobacteria in natural environment. J. Toxicol. 17, 385–403. doi: 10.3109/15569549809040400
He, X., Liu, Y. L., Conklin, A., Westrick, J., Weavers, L. K., Dionysiou, D. D., et al. (2016). Toxic cyanobacteria and drinking water: impacts, detection, and treatment. Harmful Algae. 54, 174–193. doi: 10.1016/j.hal.2016.01.001
Hernández, J. M., Lopez-Rodas, V., and Costas, E. (2009). Microcystins from tap water could be a risk factor for liver and colorectal cancer: a risk intensified by global change. Med. Hypotheses 72, 539–540. doi: 10.1016/j.mehy.2008.11.041
Ho, L., Dreyfus, J., Boyer, J., Lowe, T., Bustamante, H., Duker, P., et al. (2012). Fate of cyanobacteria and their metabolites during water treatment sludge management processes. Sci. Total Environ. 424, 232–238. doi: 10.1016/j.scitotenv.2012.02.025
Ho, L., Gaudieux, A. L., Fanok, S., Newcombe, G., and Humpage, A. R. (2007). Bacterial degradation of microcystin toxins in drinking water eliminates their toxicity. Toxicon 50, 438–441. doi: 10.1016/j.toxicon.2007.03.021
Ho, L., Meyn, T., Keegan, A., Hoefel, D., Brookes, J., Saint, C. P., et al. (2006). Bacterial degradation of microcystin toxins within a biologically active sand filter. Water Res. 40, 768–774. doi: 10.1016/j.watres.2005.12.009
Hoefel, D., Adriansen, C. M., Bouyssou, M. A., Saint, C. P., Newcombe, G., and Ho, L. (2009). Development of an mlrA gene-directed TaqMan PCR assay for quantitative assessment of microcystin-degrading bacteria within water treatment plant sand filter biofilms. Appl. Environ. Micro. 75, 5167–5169. doi: 10.1128/AEM.00036-09
Holst, T., Jørgensen, N. O., Jørgensen, C., and Johansen, A. (2003). Degradation of microcystin in sediments at oxic and anoxic, denitrifying conditions. Water Res. 37, 4748–4760. doi: 10.1016/s0043-1354(03)00413-5
Hu, L. B., Yang, J. D., Zhou, W., Yin, Y. F., Chen, J., and Shi, Z. Q. (2009). Isolation of a Methylobacillus sp. that degrades microcystin toxins associated with cyanobacteria. New Biotech. 26, 205–211. doi: 10.1016/j.nbt.2009.09.001
Jiang, Y., Shao, J., Wu, X., Xu, Y., and Li, R. (2011). Active and silent members in the mlr gene cluster of a microcystin-degrading bacterium isolated from Lake Taihu. China. FEMS Microbiol. Lett. 322, 108–114. doi: 10.1111/j.1574-6968.2011.02337.x
Kankaanpää, H. T., Holliday, J., Schröder, H., Goddard, T. J., von Fister, R., and Carmichael, W. W. (2005). Cyanobacteria and prawn farming in northern New South Wales, Australia—a case study on cyanobacteria diversity and hepatotoxin bioaccumulation. Toxicol. Appl. Pharm. 203, 243–256. doi: 10.1016/j.taap.2004.04.012
Krausfeldt, L. E., Tang, X., van de Kamp, J., Gao, G., Bodrossy, L., Boyer, G. L., et al. (2017). Spatial and temporal variability in the nitrogen cyclers of hypereutrophic Lake Taihu. FEMS Micro. Ecol. 93:fix204. doi: 10.1093/femsec/fix024
Krishnan, A., Zhang, Y.-Q., and Mou, X. (2018). Isolation and characterization of microcystin-degrading bacteria from Lake Erie. Bull. Environ. Contam Toxicol. 101, 617–623. doi: 10.1007/s00128-018-2468-4
Kumar, P., Rubio, H. D. P., Hegde, K., Brar, S. K., Cledon, M., Kermanshahi-pour, A., et al. (2019). Agro-industrial residues as a unique support in a sand filter to enhance the bioactivity to remove microcystin-Leucine aRginine and organics. Sci. Total Environ. 670, 971–981. doi: 10.1016/j.scitotenv.2019.03.260
Lam, A. K.-Y., Fedorak, P. M., and Prepas, E. E. (1995). Biotransformation of the cyanobacterial hepatotoxin microcystin-LR, as determined by HPLC and protein phosphatase bioassay. Environ. Sci. Tech. 29, 242–246. doi: 10.1021/es00001a030
Lee, Y.-J., Jung, J.-M., Jang, M.-H., Ha, K., and Joo, G.-J. (2006). Degradation of microcystins by adsorbed bacteria on a granular active carbon(GAC) filter during the water treatment process. J. Environ. Biol. 37, 317–322.
Lemes, G. A., Kersanach, R., Pinto, L. D. S., Dellagostin, O. A., Yunes, J. S., and Matthiensen, A. (2008). Biodegradation of microcystins by aquatic Burkholderia sp. from a South Brazilian coastal lagoon. Ecotoxicol. Environ. Saf. 69, 358–365. doi: 10.1016/j.ecoenv.2007.03.013
Lezcano, M. Á, Morón-López, J., Agha, R., López-Heras, I., Nozal, L., Quesada, A., et al. (2016). Presence or absence of mlr genes and nutrient concentrations co-determine the microcystin biodegradation efficiency of a natural bacterial community. Toxins 8:318. doi: 10.3390/toxins8110318
Lezcano, M. Á, Quesada, A., and El-Shehawy, R. (2018). Seasonal dynamics of microcystin-degrading bacteria and toxic cyanobacterial blooms: interaction and influence of abiotic factors. Harmful Algae 71, 19–28. doi: 10.1016/j.hal.2017.11.002
Li, J., Li, J., Shi, G., Mei, Z., Wang, R., and Li, D. (2016). Discerning biodegradation and adsorption of microcystin-LR in a shallow semi-enclosed bay and bacterial community shifts in response to associated process. Ecotoxicol. Environ. Saf. 132, 123–131. doi: 10.1016/j.ecoenv.2016.05.033
Li, J., Li, R., and Li, J. (2017). Current research scenario for microcystins biodegradation–A review on fundamental knowledge, application prospects and challenges. Sci. Total Environ. 595, 615–632. doi: 10.1016/j.scitotenv.2017.03.285
Ma, G., Pei, H., Hu, W., Xu, X., Ma, C., and Li, X. (2014). The removal of cyanobacteria and their metabolites through anoxic biodegradation in drinking water sludge. Bioresour. Tech. 165, 191–198. doi: 10.1016/j.biortech.2014.03.068
Maghsoudi, E., Fortin, N., Greer, C., Maynard, C., Pagé, A., Duy, S. V., et al. (2016). Cyanotoxin degradation activity and mlr gene expression profiles of a Sphingopyxis sp. isolated from Lake Champlain. Canada. Environ. Sci. 18, 1417–1426. doi: 10.1039/c6em00001k
Maseda, H., Shimizu, K., Doi, Y., Inamori, Y., Utsumi, M., Sugiura, N., et al. (2012). MlrA located in the inner membrane is essential for initial degradation of microcystin in Sphingopyxis sp. C-1. JPN J. Water Treat Biol. 48, 99–107. doi: 10.2521/jswtb.48.99
Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E. M., Kubal, M., et al. (2008). The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC bioinformatics 9:386. doi: 10.1186/1471-2105-9-386
Mou, X., Lu, X., Jacob, J., Sun, S., and Heath, R. (2013). Metagenomic identification of bacterioplankton taxa and pathways involved in microcystin degradation in Lake Erie. PLoS One 8:e61890. doi: 10.1371/journal.pone.0061890
Nybom, S., Dziga, D., Heikkilä, J., Kull, T., Salminen, S., and Meriluoto, J. (2012). Characterization of microcystin-LR removal process in the presence of probiotic bacteria. Toxicon 59, 171–181. doi: 10.1016/j.toxicon.2011.11.008
Pei, J., Mitchell, D. A., Dixon, J. E., and Grishin, N. V. (2011). Expansion of type II CAAX proteases reveals evolutionary origin of γ-secretase subunit APH-1. J. Mol. Bio. 410, 18–26. doi: 10.1016/j.jmb.2011.04.066
Pflugmacher, S., Wiegand, C., Oberemm, A., Beattie, K. A., Krause, E., Codd, G. A., et al. (1998). Identification of an enzymatically formed glutathione conjugate of the cyanobacterial hepatotoxin microcystin-LR: the first step of detoxication. Biochim. Biophys. Acta Gen.Sub. 1425, 527–533. doi: 10.1016/s0304-4165(98)00107-x
Rapala, J., Lahti, K., Sivonen, K., and Niemelä, S. (1994). Biodegradability and adsorption on lake sediments of cyanobacterial hepatotoxins and anatoxin-a. Lett. Appl. Micro. 19, 423–428. doi: 10.1111/j.1472-765x.1994.tb00972.x
Redouane, E. M., El, S. A. Z., El, F. K., Oufdou, K., Oudra, B., Lahrouni, M., et al. (2019). Mode of action and fate of microcystins in the complex soil-plant ecosystems. Chemosphere 225, 270–281. doi: 10.1016/j.chemosphere.2019.03.008
Saito, T., Okano, K., Park, H.-D., Itayama, T., Inamori, Y., Neilan, B. A., et al. (2003a). Detection and sequencing of the microcystin LR-degrading gene, mlrA, from new bacteria isolated from Japanese lakes. FEMS Microbiol. Lett. 229, 271–276. doi: 10.1016/s0378-1097(03)00847-4
Saito, T., Sugiura, N., Itayama, T., Inamori, Y., and Matsumura, M. (2003b). Biodegradation of microcystis and microcystins by indigenous nanoflagellates on biofilm in a practical treatment facility. Environ. Tech. 24, 143–151. doi: 10.1080/09593330309385545
Shimizu, K., Maseda, H., Okano, K., Itayama, T., Kawauchi, Y., Chen, R., et al. (2011). How microcystin-degrading bacteria express microcystin degradation activity. Lakes Reserv. 16, 169–178. doi: 10.1111/j.1440-1770.2011.00480.x
Smith, J. L., Schulz, K. L., Zimba, P. V., and Boyer, G. L. (2010). Possible mechanism for the foodweb transfer of covalently bound microcystins. Ecotoxicol. Environ. Saf. 73, 757–761. doi: 10.1016/j.ecoenv.2009.12.003
Song, L., Chen, W., Peng, L., Wan, N., Gan, N., and Zhang, X. (2007). Distribution and bioaccumulation of microcystins in water columns: a systematic investigation into the environmental fate and the risks associated with microcystins in Meiliang Bay, Lake Taihu. Water Res. 41, 2853–2864. doi: 10.1016/j.watres.2007.02.013
Steffen, M. M., Belisle, B. S., Watson, S. B., Boyer, G. L., Bourbonniere, R. A., and Wilhelm, S. W. (2015). Metatranscriptomic evidence for co-occurring top-down and bottom-up controls on toxic cyanobacterial communities. Appl. Environ. Micro. 81, 3268–3276. doi: 10.1128/AEM.04101-14
Steffen, M. M., Davis, T. W., McKay, R. M., Bullerjahn, G. S., Krausfeldt, L. E., Stough, J. M., et al. (2017). Ecophysiological examination of the Lake Erie Microcystis bloom in 2014: linkages between biology and the water supply shutdown of Toledo. Ohio. Environ. Sci. Technol. 51, 6745–6755. doi: 10.1021/acs.est.7b00856
Steffen, M. M., Li, Z., Effler, T. C., Hauser, L. J., Boyer, G. L., and Wilhelm, S. W. (2012). Comparative metagenomics of toxic freshwater cyanobacteria bloom communities on two continents. PLoS One 7:e44002. doi: 10.1371/journal.pone.0044002
Stough, J. M., Tang, X., Krausfeldt, L. E., Steffen, M. M., Gao, G., Boyer, G. L., et al. (2017). Molecular prediction of lytic vs lysogenic states for microcystis phage: metatranscriptomic evidence of lysogeny during large bloom events. PLoS One 12:e0184146. doi: 10.1371/journal.pone.0184146
Su, X., Steinman, A. D., Xue, Q., Zhao, Y., and Xie, L. (2018). Evaluating the contamination of microcystins in Lake Taihu, China: the application of equivalent total MC-LR concentration. Ecol. Indicators. 89, 445–454. doi: 10.1016/j.ecolind.2017.11.042
Svircev, Z., Krstic, S., Miladinov-Mikov, M., Baltic, V., and Vidovic, M. (2009). Freshwater cyanobacterial blooms and primary liver cancer epidemiological studies in Serbia. J. Environ. Sci. Health Part C. 27, 36–55. doi: 10.1080/10590500802668016
Tang, X., Krausfeldt, L. E., Shao, K., LeCleir, G. R., Stough, J. M. A., Gao, G., et al. (2018). Seasonal gene expression and the ecophysiological implications of toxic Microcystis aeruginosa blooms in Lake Taihu. Environ. Sci. Tech. 52, 11049–11059. doi: 10.1021/acs.est.8b01066
Thees, A., Atari, E., Birbeck, J., Westrick, J. A., and Huntley, J. F. (2018). Isolation and characterization of Lake Erie bacteria that degrade the cyanobacterial microcystin toxin MC-LR. J. Great Lakes Res. 45, 138–149. doi: 10.1016/j.jglr.2018.10.013
Tsao, S., Wei, D. J., Chang, Y. T., and Lee, J. F. (2017). Aerobic biodegradation of microcystin-LR by an indigenous bacterial mixed culture isolated in Taiwan. Int. Biodeterior. Biodegrad. 124, 101–108. doi: 10.1016/j.ibiod.2017.04.011
Tsuji, K., Naito, S., Kondo, F., Ishikawa, N., Watanabe, M. F., Suzuki, M., et al. (1994). Stability of microcystins from cyanobacteria: effect of light on decomposition and isomerization. Environ. Sci. Tech. 28, 173–177. doi: 10.1021/es00050a024
Wang, J., Wang, C., Li, J., Bai, P., Li, Q., Yuan, S. M., et al. (2018). Comparative genomics of degradative Novosphingobium strains with special reference to the microcystin-degrading Novosphingobium sp. THN1. Front. Micro. 9:2238. doi: 10.3389/fmicb.2018.02238
Wu, X., Wang, C., Tian, C., Xiao, B., and Song, L. (2015). Evaluation of the potential of anoxic biodegradation of intracellular and dissolved microcystins in lake sediments. J. Haz Mat. 286, 395–401. doi: 10.1016/j.jhazmat.2015.01.015
Yang, F., Guo, J., Huang, F., Massey, I., Huang, R., Li, Y., et al. (2018a). Removal of microcystin-LR by a novel native effective bacterial community designated as YFMCD4 isolated from Lake Taihu. Toxins 10:363. doi: 10.3390/toxins10090363
Yang, F., Massey, I. Y., Guo, J., Yang, S., Pu, Y., Zeng, W., et al. (2018b). Microcystin-LR degradation utilizing a novel effective indigenous bacterial community YFMCD1 from Lake Taihu. J. Toxicol. Environ. Health Part A 81, 184–193. doi: 10.1080/15287394.2018.1423803
Yang, F., Zhou, Y., Sun, R., Wei, H., Li, Y., Yin, L., et al. (2014a). Biodegradation of microcystin-LR and-RR by a novel microcystin-degrading bacterium isolated from Lake Taihu. Biodeg 25, 447–457. doi: 10.1007/s10532-013-9673-y
Yang, F., Zhou, Y., Yin, L., Zhu, G., Liang, G., and Pu, Y. (2014b). Microcystin-degrading activity of an indigenous bacterial strain Stenotrophomonas acidaminiphila MC-LTH2 isolated from Lake Taihu. PLoS One 9:e86216. doi: 10.1371/journal.pone.0086216
Zhang, L., Lyu, K., Wang, N., Gu, L., Sun, Y., Zhu, X., et al. (2018). Transcriptomic analysis reveals the pathways associated with resisting and degrading microcystin in Ochromonas. Environ. Sci. Technol. 52, 11102–11113. doi: 10.1021/acs.est.8b03106
Zhu, F.-P., Han, Z. L., Duan, J. L., Shi, X. S., Wang, T. T., Sheng, G. P., et al. (2019). A novel pathway for the anaerobic biotransformation of microcystin-LR using enrichment cultures. Environ. Pol. 247, 1064–1070. doi: 10.1016/j.envpol.2019.02.013
Zhu, L., Wu, Y., Song, L., and Gan, N. (2014). Ecological dynamics of toxic Microcystis spp. and microcystin-degrading bacteria in Dianchi Lake. China. Appl. Env. Micro. 80, 1874–1881. doi: 10.1128/AEM.02972-13
Zhu, X., Shen, Y., Chen, X., Hu, Y. O. O., Xiang, H., Tao, J., et al. (2016). Biodegradation mechanism of microcystin-LR by a novel isolate of Rhizobium sp. TH and the evolutionary origin of the mlrA gene. Int. Biodeter. Biodegrad. 115, 17–25. doi: 10.1016/j.ibiod.2016.07.011
Zurawell, R. W., Chen, H., Burke, J. M., and Prepas, E. E. (2005). Hepatotoxic cyanobacteria: a review of the biological importance of microcystins in freshwater environments. J. Toxicol. Environ. Health Part B. 8, 1–37. doi: 10.1080/10937400590889412
Keywords: biodegradation, cyanotoxins, cyanobacteria, harmful algal blooms, RNA-sequencing
Citation: Krausfeldt LE, Steffen MM, McKay RM, Bullerjahn GS, Boyer GL and Wilhelm SW (2019) Insight Into the Molecular Mechanisms for Microcystin Biodegradation in Lake Erie and Lake Taihu. Front. Microbiol. 10:2741. doi: 10.3389/fmicb.2019.02741
Received: 02 July 2019; Accepted: 11 November 2019;
Published: 10 December 2019.
Edited by:Pia H. Moisander, University of Massachusetts Dartmouth, United States
Reviewed by:Fei Yang, Central South University, China
Rachel Susan Poretsky, University of Illinois at Chicago, United States
Copyright © 2019 Krausfeldt, Steffen, McKay, Bullerjahn, Boyer and Wilhelm. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Steven W. Wilhelm, firstname.lastname@example.org