Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed

The biological functions of microcystins

Microcystins are potent hepatotoxins predominantly produced by bloom-forming freshwater cyanobacteria (e.g., Microcystis, Planktothrix, Dolichospermum). Microcystin biosynthesis involves large multienzyme complexes and tailoring enzymes encoded by the mcy gene cluster. Mutation, recombination, and deletion events have shaped the mcy gene cluster in the course of evolution, resulting in a large diversity of microcystin congeners and the natural coexistence of toxic and non-toxic strains. The biological functions of microcystins and their association with algal bloom formation have been extensively investigated over the past decades. This review synthesizes recent advances in decoding the biological role of microcystins in carbon/nitrogen metabolism, antioxidation, colony formation, and cell-to-cell communication. Microcystins appear to adopt multifunctional roles in cyanobacteria that reflect the adaptive plasticity of toxic cyanobacteria to changing environments.

Metabologenomics reveals strain-level genetic and chemical diversity of Microcystis secondary metabolism

Microcystis spp. are renowned for producing the hepatotoxin microcystin in freshwater cyanobacterial harmful algal blooms around the world, threatening drinking water supplies and public and environmental health. However, Microcystis genomes also harbor numerous biosynthetic gene clusters (BGCs) encoding the biosynthesis of other secondary metabolites, including many with toxic properties. Most of these BGCs are uncharacterized and currently lack links to biosynthesis products. However, recent field studies show that many of these BGCs are abundant and transcriptionally active in natural communities, suggesting potentially important yet unknown roles in bloom ecology and water quality. Here, we analyzed 21 xenic Microcystis cultures isolated from western Lake Erie to investigate the diversity of the biosynthetic potential of this genus. Through metabologenomic and in silico approaches, we show that these Microcystis strains contain variable BGCs, previously observed in natural populations, and encode distinct metabolomes across cultures. Additionally, we find that the majority of metabolites and gene clusters are uncharacterized, highlighting our limited understanding of the chemical repertoire of Microcystis spp. Due to the complex metabolomes observed in culture, which contain a wealth of diverse congeners as well as unknown metabolites, these results underscore the need to deeply explore and identify secondary metabolites produced by Microcystis beyond microcystins to assess their impacts on human and environmental health.IMPORTANCEThe genus Microcystis forms dense cyanobacterial harmful algal blooms (cyanoHABs) and can produce the toxin microcystin, which has been responsible for drinking water crises around the world. While microcystins are of great concern, Microcystis also produces an abundance of other secondary metabolites that may be of interest due to their potential for toxicity, ecological importance, or pharmaceutical applications. In this study, we combine genomic and metabolomic approaches to study the genes responsible for the biosynthesis of secondary metabolites as well as the chemical diversity of produced metabolites in Microcystis strains from the Western Lake Erie Culture Collection. This unique collection comprises Microcystis strains that were directly isolated from western Lake Erie, which experiences substantial cyanoHAB events annually and has had negative impacts on drinking water, tourism, and industry.

Spatio-temporal connectivity of a toxic cyanobacterial community and its associated microbiome along a freshwater-marine continuum

Due to climate changes and eutrophication, blooms of predominantly toxic freshwater cyanobacteria are intensifying and are likely to colonize estuaries, thus impacting benthic organisms and shellfish farming representing a major ecological, health and economic risk. In the natural environment, Microcystis form large mucilaginous colonies that influence the development of both cyanobacterial and embedded bacterial communities. However, little is known about the fate of natural colonies of Microcystis by salinity increase. In this study, we monitored the fate of a Microcystis dominated bloom and its microbiome along a French freshwater-marine gradient at different phases of a bloom. We demonstrated changes in the cyanobacterial genotypic composition, in the production of specific metabolites (toxins and compatible solutes) and in the heterotrophic bacteria structure in response to the salinity increase. In particular M. aeruginosa and M. wesenbergii survived salinities up to 20. Based on microcystin gene abundance, the cyanobacteria became more toxic during their estuarine transfer but with no selection of specific microcystin variants. An increase in compatible solutes occurred along the continuum with extensive trehalose and betaine accumulations. Salinity structured most the heterotrophic bacteria community, with an increased in the richness and diversity along the continuum. A core microbiome in the mucilage-associated attached fraction was highly abundant suggesting a strong interaction between Microcystis and its microbiome and a likely protecting role of the mucilage against an osmotic shock. These results underline the need to better determine the interactions between the Microcystis colonies and their microbiome as a likely key to their widespread success and adaptation to various environmental conditions.

Comparative genomic analysis of Microcystis strain diversity using conserved marker genes

Microcystis-dominated cyanobacterial harmful algal blooms (cyanoHABs) have a global impact on freshwater environments, affecting both wildlife and human health. Microcystis diversity and function in field samples and laboratory cultures can be determined by sequencing whole genomes of cultured isolates or natural populations, but these methods remain computationally and financially expensive. Amplicon sequencing of marker genes is a lower cost and higher throughput alternative to characterize strain composition and diversity in mixed samples. However, the selection of appropriate marker gene region(s) and primers requires prior understanding of the relationship between single gene genotype, whole genome content, and phenotype. To identify phylogenetic markers of Microcystis strain diversity, we compared phylogenetic trees built from each of 2,351 individual core genes to an established phylogeny and assessed the ability of these core genes to predict whole genome content and bioactive compound genotypes. We identified single-copy core genes better able to resolve Microcystis phylogenies than previously identified marker genes. We developed primers suitable for current Illumina-based amplicon sequencing with near-complete coverage of available Microcystis genomes and demonstrate that they outperform existing options for assessing Microcystis strain composition. Results showed that genetic markers can be used to infer Microcystis gene content and phenotypes such as potential production of bioactive compounds , although marker performance varies by bioactive compound gene and sequence similarity. Finally, we demonstrate that these markers can be used to characterize the Microcystis strain composition of laboratory or field samples like those collected for surveillance and modeling of Microcystis-dominated cyanobacterial harmful algal blooms.

Two duckweed species exhibit variable tolerance to microcystin-LR exposure across genotypic lineages

Cyanotoxins produced by harmful cyanobacteria blooms can damage freshwater ecosystems and threaten human health. Floating macrophytes may be used as a means of biocontrol by limiting light and resources available to cyanobacteria. However, genetic variation in macrophyte sensitivity to cyanotoxins could influence their suitability as biocontrol agents. We investigated the influence of such intraspecific variation on the response of two rapidly growing duckweed species, Lemna minor and Spirodela polyrhiza, often used in nutrient and metal bioremediation. We assessed two biomarkers related to productivity (biomass and chlorophyll A production) and two related to fitness measures (population size and growth rate). Fifteen genetic lineages of each species were grown in media containing common cyanotoxin microcystin-LR at ecologically relevant concentrations or control media for a period of twelve days. Genotype identity had a strong impact on all biomarker responses. Microcystin concentration slightly increased the final population sizes of both macrophyte species with a marginal effect on growth rate of L. minor and the chlorophyll A production of S. polyrhiza, but overall these species were very tolerant of microcystin. The strong tolerance supports the potential use of these plants as bioremediators of cyanobacterial blooms. However, differential impact of microcystin exposure discovered in single lineage models among genotypes indicates a potential for cyanotoxins to act as selective forces, necessitating attention to genotype selection for bioremediation.

GST-Mu of Cristaria plicata is regulated by Nrf2/Keap1 pathway in detoxification microcystin and has antioxidant function

Glutathione S-transferase is a crucial phase II metabolic enzyme involved in detoxification and metabolism in aquatic organisms. This study aimed to investigate the regulation of Nrf2/Keap1 pathway on microcystin-induced CpGST-Mu expression and CpGST-Mu resistance to hydrogen peroxide. A mu class GST from Cristaria plicata (CpGST-Mu) was identified. The full-length cDNA was 1026 bp, with an open reading frame of 558 bp. Subcellular localization revealed that CpGST-Mu was localized in cytoplasm. The optimum pH and temperature for the catalytic activity of CpGST-Mu protein was pH 6 and 40 °C, respectively. The results of Real-time quantitative PCR showed that CpGST-Mu mRNA was constitutively expressed in tissues, with the highest expression level in hepatopancreas and the lowest expression level in gill. The mRNA level of CpGST-Mu was significantly increased under the stress of microcystins and hydrogen peroxide. CpGST-Mu had an antagonistic effect on hydrogen peroxide. In the knockdown experiments, the mRNA levels of CpGST-Mu exhibited corresponding changes while Nrf2 and Keap1 genes were individually knocked down. These findings indicated that GST-Mu exhibited antioxidant properties and its expression was regulated by Nrf2/Keap1 signaling pathway. The study provided new information on the function of GST-Mu and could contribute to future studies on how to excrete microcystins in molluscs.

Microcystin congeners in Lake Erie follow the seasonal pattern of nitrogen availability

Cyanobacteria harmful algal blooms produce many toxic secondary metabolites called cyanotoxins. The most studied group of cyanotoxins are microcystins (MC), with over 300 congeners reported. MC-LR is the most studied congener because of its abundance and toxicity. Recent toxicology studies suggest that more hydrophobic MC congeners such as MC-LA, MC-LF, and MC-LW may be less abundant but up to seven times more toxic than MC-LR, whereas, MC-RR's toxicity is only one-fifth that of MC-LR. Hence, understanding the environmental stressors that change the MC congener profile is critical to assessing the negative impact on environmental and human health. A two-year field and experimental study investigated seasonal and spatial changes of MC congener profiles in the western basin of Lake Erie. Both studies showed that nitrogen enrichment favored the production of nitrogen-rich MC-RR (C49H75N13O12). The field study showed that nitrogen depletion favored the low-nitrogen MC-LA (C46H67N7O12). MC-LR (a medium N level, C49H75N10O12) accounted for ∼30% to 50% of the total MC concentration and was stable across nitrogen concentrations. Using the relative toxicity and concentrations of each MC congener, both LC-MS/MS and ELISA overestimated the toxicity early bloom (July) and underestimated it late bloom (September). On 24 July 2019, highly toxic MC-LW and MC-LF were detected at nearshore stations with relative toxicity exceeding drinking water standards. This study demonstrated that the less toxic, high nitrogen MC-RR dominated under nitrogen-replete conditions in the early season, whereas the more toxic, less nitrogen MC-LA dominated under nitrogen-limited conditions later in the season.

The Western Lake Erie culture collection: A promising resource for evaluating the physiological and genetic diversity of Microcystis and its associated microbiome

Cyanobacterial harmful algal blooms (cyanoHABs) dominated by Microcystis spp. have significant public health and economic implications in freshwater bodies around the world. These blooms are capable of producing a variety of cyanotoxins, including microcystins, that affect fishing and tourism industries, human and environmental health, and access to drinking water. In this study, we isolated and sequenced the genomes of 21 primarily unialgal Microcystis cultures collected from western Lake Erie between 2017 and 2019. While some cultures isolated in different years have a high degree of genetic similarity (genomic Average Nucleotide Identity >99%), genomic data show that these cultures also represent much of the breadth of known Microcystis diversity in natural populations. Only five isolates contained all the genes required for microcystin biosynthesis while two isolates contained a previously described partial mcy operon. Microcystin production within cultures was also assessed using Enzyme-Linked Immunosorbent Assay (ELISA) and supported genomic results with high concentrations (up to 900 μg L⁻¹) in cultures with complete mcy operons and no or low toxin detected otherwise. These xenic cultures also contained a substantial diversity of bacteria associated with Microcystis, which has become increasingly recognized as an essential component of cyanoHAB community dynamics. These results highlight the genomic diversity among Microcystis strains and associated bacteria in Lake Erie, and their potential impacts on bloom development, toxin production, and toxin degradation. This culture collection significantly increases the availability of environmentally relevant Microcystis strains from temperate North America.

Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

Abundance trade-offs and dominant taxa maintain the stability of the bacterioplankton community underlying Microcystis blooms

Microcystis blooms are an intractable global environmental problem that pollute water and compromise ecosystem functioning. Closed-lake management practices keep lakes free of sewage and harmful algae invasions and have succeeded in controlling local Microcystis blooms; however, there is little understanding of how the bacterioplankton communities associated with Microcystis have changed. Here, based on metagenomic sequencing, the phyla, genera, functional genes and metabolic functions of the bacterioplankton communities were compared between open lakes (underlying Microcystis blooms) and closed lakes (no Microcystis blooms). Water properties and zooplankton density were investigated and measured as factors influencing blooms. The results showed that (1) the water quality of closed lakes was improved, and the nitrogen and phosphorus concentrations were significantly reduced. (2) The stability of open vs. closed-managed lakes differed notably at the species and genus levels (p < 0.01), but no significant variations were identified at the phylum and functional genes levels (p > 0.05). (3) The relative abundance of Microcystis (Cyanobacteria) increased dramatically in the open lakes (proportions from 1.44 to 41.76%), whereas the relative abundance of several other dominant genera of Cyanobacteria experienced a trade-off and decreased with increasing Microcystis relative abundance. (4) The main functions of the bacterioplankton communities were primarily related to dominant genera of Proteobacteria and had no significant relationship with Microcystis. Overall, the closed-lake management practices significantly reduced nutrients and prevented Microcystis blooms, but the taxonomic and functional structures of bacterioplankton communities remained stable overall.

Transcriptionally active nitrogen fixation and biosynthesis of diverse secondary metabolites by Dolichospermum and Aphanizomenon-like Cyanobacteria in western Lake Erie Microcystis blooms

Cyanobacterial harmful algal blooms (cyanoHABs) in the western basin of Lake Erie are dominated by microcystin producing Microcystis spp., but other cyanobacterial taxa that coexist in these communities may play important roles in production of toxins and shaping bloom dynamics and community function. In this study, we used metagenomic and metatranscriptomic data from the 2014 western Lake Erie cyanoHAB to explore the genetic diversity and biosynthetic potential of cyanobacteria belonging to the Anabaena, Dolichospermum, Aphanizomenon (ADA) clade. We reconstructed two near-complete metagenome-assembled genomes from two distinct ADA clade species, each containing biosynthetic gene clusters that encode novel and known secondary metabolites, including those with toxic and/or known taste and odor properties, that were transcriptionally active. However, neither ADA metagenome-assembled genome contained genes encoding guanitoxins, anatoxins, or saxitoxins, which are known to be produced by ADA. The ADA cyanobacteria accounted for most of the metagenomic and metatranscriptomic reads from nitrogen fixation genes, suggesting they were the dominant N-fixers at the times and stations sampled. Despite their relatively low abundance, our results highlight the possibility that ADA taxa could influence the water quality and ecology of Microcystis blooms, although the extent of these impacts remains to be quantified.

Widespread formation of intracellular calcium carbonates by the bloom-forming cyanobacterium Microcystis

The formation of intracellular amorphous calcium carbonates (iACC) has been recently observed in a few cultured strains of Microcystis, a potentially toxic bloom-forming cyanobacterium found worldwide in freshwater ecosystems. If iACC-forming Microcystis are abundant within blooms, they may represent a significant amount of particulate Ca. Here, we investigate the significance of iACC biomineralization by Microcystis. First, the presence of iACC-forming Microcystis cells has been detected in several eutrophic lakes, indicating that this phenomenon occurs under environmental conditions. Second, some genotypic (presence/absence of ccyA, a marker gene of iACC biomineralization) and phenotypic (presence/absence of iACC) diversity have been detected within a collection of strains isolated from one single lake. This illustrates that this trait is frequent but also variable within Microcystis even at a single locality. Finally, one-third of publicly available genomes of Microcystis were shown to contain the ccyA gene, revealing a wide geographic and phylogenetic distribution within the genus. Overall, the present work shows that the formation of iACC by Microcystis is common under environmental conditions. While its biological function remains undetermined, this process should be further considered regarding the biology of Microcystis and implications on the Ca geochemical cycle in freshwater environments.

Effect of hydrological modification on the potential toxicity of Microcystis aeruginosa complex in Salto Grande reservoir, Uruguay

It is widely known that the environmental conditions caused by the construction of reservoirs favor the proliferation of toxic cyanobacteria and the formation of blooms due to the high residence time of the water, low turbidity, temperature regimes, among others. Microcystin-producing cyanobacteria such as those from the Microcystis aeruginosa complex (MAC) are the most frequently found organisms in reservoirs worldwide, being the role of the environment on microcystin production poorly understood. Here, we addressed the community dynamics and potential toxicity of MAC cyanobacteria in a subtropical reservoir (Salto Grande) located in the low Uruguay river. Samples were taken from five different sites (upstream, inside the reservoir and downstream) during contrasting seasons (summer and winter) to analyze: (i) the MAC community structure by amplicon sequencing of the phycocyanin gene spacer, (ii) the genotype diversity of microcystin-producing MAC by high resolution melting analysis of the mcyJ gene, and (iii) the abundance and mcy transcription activity of the microcystin-producing (toxic) fraction. We found that MAC diversity decreased from summer to winter but, despite the observed changes in MAC community structure, the abundance of toxic organisms and the transcription of mcy genes were always higher inside the reservoir, regardless of the season. Two different genotypes of toxic MAC were detected inside the reservoir, one associated with low water temperature (15 °C) and one thriving at high water temperature (31 °C). These findings indicate that the environmental conditions inside the reservoir reduce community diversity while promoting the proliferation of toxic genotypes that actively transcribe mcy genes, whose relative abundance will depend on the water temperature.

Trait biases in microbial reference genomes

Common culturing techniques and priorities bias our discovery towards specific traits that may not be representative of microbial diversity in nature. So far, these biases have not been systematically examined. To address this gap, here we use 116,884 publicly available metagenome-assembled genomes (MAGs, completeness ≥80%) from 203 surveys worldwide as a culture-independent sample of bacterial and archaeal diversity, and compare these MAGs to the popular RefSeq genome database, which heavily relies on cultures. We compare the distribution of 12,454 KEGG gene orthologs (used as trait proxies) in the MAGs and RefSeq genomes, while controlling for environment type (ocean, soil, lake, bioreactor, human, and other animals). Using statistical modeling, we then determine the conditional probabilities that a species is represented in RefSeq depending on its genetic repertoire. We find that the majority of examined genes are significantly biased for or against in RefSeq. Our systematic estimates of gene prevalences across bacteria and archaea in nature and gene-specific biases in reference genomes constitutes a resource for addressing these issues in the future.

Cryptic Diversity of Black Band Disease Cyanobacteria in Siderastrea siderea Corals Revealed by Chemical Ecology and Comparative Genome-Resolved Metagenomics

Black band disease is a globally distributed and easily recognizable coral disease. Despite years of study, the etiology of this coral disease, which impacts dozens of stony coral species, is not completely understood. Although black band disease mats are predominantly composed of the cyanobacterial species Roseofilum reptotaenium, other filamentous cyanobacterial strains and bacterial heterotrophs are readily detected. Through chemical ecology and metagenomic sequencing, we uncovered cryptic strains of Roseofilum species from Siderastrea siderea corals that differ from those on other corals in the Caribbean and Pacific. Isolation of metabolites from Siderastrea-derived Roseofilum revealed the prevalence of unique forms of looekeyolides, distinct from previously characterized Roseofilum reptotaenium strains. In addition, comparative genomics of Roseofilum strains showed that only Siderastrea-based Roseofilum strains have the genetic capacity to produce lasso peptides, a family of compounds with diverse biological activity. All nine Roseofilum strains examined here shared the genetic capacity to produce looekeyolides and malyngamides, suggesting these compounds support the ecology of this genus. Similar biosynthetic gene clusters are not found in other cyanobacterial genera associated with black band disease, which may suggest that looekeyolides and malyngamides contribute to disease etiology through yet unknown mechanisms.

Microcystis pangenome reveals cryptic diversity within and across morphospecies

Microcystis, a common harmful algal bloom (HAB) taxon, threatens water supplies and human health, yet species delimitation is contentious in this taxon, leading to challenges in research and management of this threat. Historical and common morphology-based classifications recognize multiple morphospecies, most with variable and diverse ecologies, while DNA sequence-based classifications indicate a single species with multiple ecotypes. To better delimit Microcystis species, we conducted a pangenome analysis of 122 genomes. Core- and non-core gene phylogenetic analyses placed 113 genomes into 23 monophyletic clusters containing at least two genomes. Overall, genome-related indices revealed that Microcystis contains at least 16 putative genospecies. Fifteen genospecies included at least one Microcystis aeruginosa morphospecies, and 10 genospecies included two or more morphospecies. This classification system will enable consistent taxonomic identification of Microcystis and thereby aid in resolving some of the complexities and controversies that have long characterized eco-evolutionary research and management of this important HAB taxon.

Not all bad: Gyromitrin has a limited distribution in the false morels as determined by a new ultra high-performance liquid chromatography method

Gyromitrin (acetaldehyde N-methyl-N-formylhydrazone) and its homologs are deadly mycotoxins produced most infamously by the lorchel (also known as false morel) Gyromitra esculenta, which is paradoxically consumed as a delicacy in some parts of the world. There is much speculation about the presence of gyromitrin in other species of the lorchel family (Discinaceae), but no studies have broadly assessed its distribution. Given the history of poisonings associated with the consumption of G. esculenta and G. ambigua, we hypothesized that gyromitrin evolved in the last common ancestor of these taxa and would be present in their descendants with adaptive loss of function in the nested truffle clade, Hydnotrya. To test this hypothesis, we developed a sensitive analytical derivatization method for the detection of gyromitrin using 2,4-dinitrobenzaldehyde as the derivatization reagent. In total, we analyzed 66 specimens for the presence of gyromitrin over 105 tests. Moreover, we sequenced the nuc rDNA internal transcribed spacer region ITS1-5.8S-ITS2 (ITS barcode) and nuc 28S rDNA to assist in species identification and to infer a supporting phylogenetic tree. We detected gyromitrin in all tested specimens from the G. esculenta group as well as G. leucoxantha. This distribution is consistent with a model of rapid evolution coupled with horizontal transfer, which is typical for secondary metabolites. We clarified that gyromitrin production in Discinaceae is both discontinuous and more limited than previously thought. Further research is required to elucidate the gyromitrin biosynthesis gene cluster and its evolutionary history in lorchels.