Evidence for positive selection acting on microcystin synthetase adenylation domains in three cyanobacterial genera

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.

Identification of a homoarginine biosynthetic gene from a microcystin biosynthetic pathway in Fischerella sp. PCC 9339

Microcystins (MCs) are a family of chemically diverse toxins produced by numerous distantly related cyanobacteria. They are potent inhibitors of eukaryotic protein phosphatases 1 and 2A and are responsible for the toxicosis and death of wild and domestic animals around the world. Microcystins are synthesized on large enzyme complexes comprised of peptide synthetases, polyketide synthases, and additional modifying enzymes. Bioinformatic analysis identified the presence of an additional uncharacterized enzyme in the microcystin (mcy) biosynthetic gene cluster in Fischerella sp. PCC 9339, which we named McyK, that lacked a clearly defined role in the biosynthesis of microcystin. Further bioinformatic analysis suggested that McyK belongs to the inosamine-phosphate amidinotransferase family and could be involved in synthesizing homo amino acids. Quadrupole time-of-flight tandem mass spectrometry (Q-TOFMS/MS) analysis confirmed that Fischerella sp. PCC 9339 produces MC-Leucine2-Homoarginine4(MC-LHar) and [Aspartic acid3]MC-Leucine2-Homoarginine4 ([Asp3]MC-LHar) as the dominant chemical variants. We hypothesized that the McyK enzyme might be involved in the production of microcystin variants containing homoarginine (Har) in the strain. Heterologous expression of a codon-optimized mcyK gene in Escherichia coli confirmed that McyK is responsible for the synthesis of L-Har. These results confirm the production of MC-LHar, a novel microcystin chemical variant [Asp3]MC-LHar, and a new microcystin biosynthetic enzyme involved in supply of the rare homo-amino acid Har to the microcystin biosynthetic pathway in Fischerella sp. PCC 9339. This study provides new insights into the logic underpinning the biosynthesis of microcystin chemical variants and broadens our knowledge of structural diversity of the microcystin family of toxins.

Cyanotoxin-encoding genes as powerful predictors of cyanotoxin production during harmful cyanobacterial blooms in an inland freshwater lake: Evaluating a novel early-warning system

Freshwater harmful cyanobacterial blooms (HCBs) potentially produce excessive cyanotoxins, mainly microcystins (MCs), significantly threatening aquatic ecosystems and public health. Accurately predicting HCBs is thus essential to developing effective HCB mitigation and prevention strategies. We previously developed a novel early-warning system that uses cyanotoxin-encoding genes to predict cyanotoxin production in Harsha Lake, Ohio, USA, in 2015. In this study, we evaluated the efficacy of the early-warning system in forecasting the 2016 HCB in the same lake. We also examined potential HCB drivers and cyanobacterial community composition. Our results revealed that the cyanobacterial community was stable at the phylum level but changed dynamically at the genus level over time. Microcystis and Planktothrix were the major MC-producing genera that thrived in June and July and produced high concentrations of MCs (peak level 10.22 μg·L^-1). The abundances of the MC-encoding gene cluster mcy and its transcript levels significantly correlated with total MC concentrations (before the MC concentrations peaked) and accurately predicted MC production as revealed by logistic equations. When the Microcystis-specific gene mcyG reached approximately 1.5 × 10³ copies·mL^-1 or when its transcript level reached approximately 2.4 copies·mL^-1, total MC level exceeded 0.3 μg L^-1 (a health advisory limit) approximately one week later (weekly sampling scheme). This study suggested that cyanotoxin-encoding genes are promising predictors of MC production in inland freshwater lakes, such as Harsha Lake. The evaluated early-warning system can be a useful tool to assist lake managers in predicting, mitigating, and/or preventing HCBs.

Mining of Cyanobacterial Genomes Indicates Natural Product Biosynthetic Gene Clusters Located in Conjugative Plasmids

Microbial natural products are compounds with unique chemical structures and diverse biological activities. Cyanobacteria commonly possess a wide range of biosynthetic gene clusters (BGCs) to produce natural products. Although natural product BGCs have been found in almost all cyanobacterial genomes, little attention has been given in cyanobacterial research to the partitioning of these biosynthetic pathways in chromosomes and plasmids. Cyanobacterial plasmids are believed to disperse several natural product BGCs, such as toxins, by plasmids through horizontal gene transfer. Therefore, plasmids may confer the ability to produce toxins and may play a role in the evolution of diverse natural product BGCs from cyanobacteria. Here, we performed an analysis of the distribution of natural product BGCs in 185 genomes and mapped the presence of genes involved in the conjugation in plasmids. The 185 analyzed genomes revealed 1817 natural products BGCs. Individual genomes contained 1–42 biosynthetic pathways (mean 8), 95% of which were present in chromosomes and the remaining 5% in plasmids. Of the 424 analyzed cyanobacterial plasmids, 12% contained homologs of genes involved in conjugation and natural product biosynthetic pathways. Among the biosynthetic pathways in plasmids, manual curation identified those to produce aeruginosin, anabaenopeptin, ambiguine, cryptophycin, hassallidin, geosmin, and microcystin. These compounds are known toxins, protease inhibitors, odorous compounds, antimicrobials, and antitumorals. The present study provides in silico evidence using genome mining that plasmids may be involved in the distribution of natural product BGCs in cyanobacteria. Consequently, cyanobacterial plasmids have importance in the context of biotechnology, water management, and public health risk assessment. Future research should explore in vivo conjugation and the end products of natural product BGCs in plasmids via chemical analyses.

The Landscape of Recombination Events That Create Nonribosomal Peptide Diversity

Nonribosomal peptides (NRP) are crucial molecular mediators in microbial ecology and provide indispensable drugs. Nevertheless, the evolution of the flexible biosynthetic machineries that correlates with the stunning structural diversity of NRPs is poorly understood. Here, we show that recombination is a key driver in the evolution of bacterial NRP synthetase (NRPS) genes across distant bacterial phyla, which has guided structural diversification in a plethora of NRP families by extensive mixing and matching of biosynthesis genes. The systematic dissection of a large number of individual recombination events did not only unveil a striking plurality in the nature and origin of the exchange units but allowed the deduction of overarching principles that enable the efficient exchange of adenylation (A) domain substrates while keeping the functionality of the dynamic multienzyme complexes. In the majority of cases, recombination events have targeted variable portions of the Acore domains, yet domain interfaces and the flexible Asub domain remained untapped. Our results strongly contradict the widespread assumption that adenylation and condensation (C) domains coevolve and significantly challenge the attributed role of C domains as stringent selectivity filter during NRP synthesis. Moreover, they teach valuable lessons on the choice of natural exchange units in the evolution of NRPS diversity, which may guide future engineering approaches.

Algal Toxic Compounds and Their Aeroterrestrial, Airborne and other Extremophilic Producers with Attention to Soil and Plant Contamination: A Review

The review summarizes the available knowledge on toxins and their producers from rather disparate algal assemblages of aeroterrestrial, airborne and other versatile extreme environments (hot springs, deserts, ice, snow, caves, etc.) and on phycotoxins as contaminants of emergent concern in soil and plants. There is a growing body of evidence that algal toxins and their producers occur in all general types of extreme habitats, and cyanobacteria/cyanoprokaryotes dominate in most of them. Altogether, 55 toxigenic algal genera (47 cyanoprokaryotes) were enlisted, and our analysis showed that besides the “standard” toxins, routinely known from different waterbodies (microcystins, nodularins, anatoxins, saxitoxins, cylindrospermopsins, BMAA, etc.), they can produce some specific toxic compounds. Whether the toxic biomolecules are related with the harsh conditions on which algae have to thrive and what is their functional role may be answered by future studies. Therefore, we outline the gaps in knowledge and provide ideas for further research, considering, from one side, the health risk from phycotoxins on the background of the global warming and eutrophication and, from the other side, the current surge of interest which phycotoxins provoke due to their potential as novel compounds in medicine, pharmacy, cosmetics, bioremediation, agriculture and all aspects of biotechnological implications in human life.

Coherence of Microcystis species revealed through population genomics

Microcystis is a genus of freshwater cyanobacteria, which causes harmful blooms in ecosystems worldwide. Some Microcystis strains produce harmful toxins such as microcystin, impacting drinking water quality. Microcystis colony morphology, rather than genetic similarity, is often used to classify Microcystis into morphospecies. Yet colony morphology is a plastic trait, which can change depending on environmental and laboratory culture conditions, and is thus an inadequate criterion for species delineation. Furthermore, Microcystis populations are thought to disperse globally and constitute a homogeneous gene pool. However, this assertion is based on relatively incomplete characterization of Microcystis genomic diversity. To better understand these issues, we performed a population genomic analysis of 33 newly sequenced genomes mainly from Canada and Brazil. We identified 17 Microcystis clusters of genomic similarity, five of which correspond to monophyletic clades containing at least three newly sequenced genomes. Four out of these five clades match to named morphospecies. Notably, M. aeruginosa is paraphyletic, distributed across 12 genomic clusters, suggesting it is not a coherent species. A few clades of closely related isolates are specific to a unique geographic location, suggesting biogeographic structure over relatively short evolutionary time scales. Higher homologous recombination rates within than between clades further suggest that monophyletic groups might adhere to a Biological Species-like concept, in which barriers to gene flow maintain species distinctness. However, certain genes-including some involved in microcystin and micropeptin biosynthesis-are recombined between monophyletic groups in the same geographic location, suggesting local adaptation.

The Biosynthesis of Rare Homo-Amino Acid Containing Variants of Microcystin by a Benthic Cyanobacterium

Microcystins are a family of chemically diverse hepatotoxins produced by distantly related cyanobacteria and are potent inhibitors of eukaryotic protein phosphatases 1 and 2A. Here we provide evidence for the biosynthesis of rare variants of microcystin that contain a selection of homo-amino acids by the benthic cyanobacterium Phormidium sp. LP904c. This strain produces at least 16 microcystin chemical variants many of which contain homophenylalanine or homotyrosine. We retrieved the complete 54.2 kb microcystin (mcy) gene cluster from a draft genome assembly. Analysis of the substrate specificity of McyB₁ and McyC adenylation domain binding pockets revealed divergent substrate specificity sequences, which could explain the activation of homo-amino acids which were present in 31% of the microcystins detected and included variants such as MC-LHty, MC-HphHty, MC-LHph and MC-HphHph. The mcy gene cluster did not encode enzymes for the synthesis of homo-amino acids but may instead activate homo-amino acids produced during the synthesis of anabaenopeptins. We observed the loss of microcystin during cultivation of a closely related strain, Phormidium sp. DVL1003c. This study increases the knowledge of benthic cyanobacterial strains that produce microcystin variants and broadens the structural diversity of known microcystins.

Interactions between Microcystis aeruginosa and coexisting bisphenol A at different nitrogen levels

Microcystis aeruginosa is known as a main contributor of cyanobacterial bloom. However, factors that drive its formation and dispersion remain poorly understood. The cellular-level responses to nutrient drivers of eutrophication were investigated. The results showed that growth rate of M. aeruginosa was significantly enhanced with the increasing bisphenol A (BPA) and nitrogen (N) level. Stress of BPA significantly inhibited cellular density, chlorophyll-a content across all the nutrient conditions, while Fv/Fm and rETRmax value were promoted by BPA. Responses of reactive oxygen species (ROS) value, superoxide dismutase (SOD) activity and malodialdehyde (MDA) content indicated that nitrogen deficiency and BPA caused oxidative stress to M. aeruginosa. Besides, nitrogen and BPA regulated the production and release of microcystins (MCs). M. aeruginosa exposed to BPA caused 95 up-regulated proteins, which was primarily associated with photosynthesis, nitrogen metabolism, glycolysis/glyconeogenesis and carbon fixation in photosynthetic organisms. The 91 down-regulated proteins were related to quorum sensing, longevity regulating and cell cycle-caulobacter, confirming that the driving force of regulating the change of cellular density and genes expression weakened. These findings provide important clues to elucidate the combined regulatory mechanisms of cyanobacterial blooms triggered by endocrine-disrupting compounds and environmental factors and help to effectively prevent and reduce cyanobacterial blooms.

Spatial-temporal survey of Microcystis oligopeptide chemotypes in reservoirs with dissimilar waterbody features and their relation to genetic variation

The ability of cyanobacteria to produce toxins and other secondary metabolites is patchily distributed in natural populations, enabling the use of cellular oligopeptide compositions as markers to classify strains into ecologically-relevant chemotypical subpopulations. The composition and spatiotemporal distribution of Microcystis chemotypes within and among waterbodies was studied at different time scales by analyzing (i) Microcystis strains isolated between 1998 and 2007 from different Spanish reservoirs and (ii) individual Microcystis aeruginosa colonies collected from pelagic and littoral habitats in Valmayor reservoir (Spain) during a bloom. No agreement between chemotypes and both morphotypes and genotypes (based on cpcBA-IGS, 16S-23S rRNA ITS and mcyB genes) was found, suggesting that oligopeptide profiles in individual strains evolve independently across morphospecies and phylogenetic genotypes, and that the diversity of microcystin variants produced cannot be explained by mcyB gene variations alone. The presence of identical chemotypes in spatially-distant reservoirs with dissimilar trophic state, lithology or depth indicate that waterbody characteristics and geographical boundaries weakly affect chemotype composition and distribution. At smaller spatiotemporal scales (i.e. during bloom), M. aeruginosa populations showed high number of chemotypes, as well as marked differences in chemotype composition and relative abundance among the littoral and pelagic habitats. This indicates that the factors influencing chemotype composition, relative abundance and dynamics operate at short spatial and temporal scales, and supports emerging hypotheses about interactions with antagonistic microorganisms as possible drivers for widespread chemical polymorphisms in cyanobacteria.

Consensus assessment of the contamination level of publicly available cyanobacterial genomes

Publicly available genomes are crucial for phylogenetic and metagenomic studies, in which contaminating sequences can be the cause of major problems. This issue is expected to be especially important for Cyanobacteria because axenic strains are notoriously difficult to obtain and keep in culture. Yet, despite their great scientific interest, no data are currently available concerning the quality of publicly available cyanobacterial genomes. As reliably detecting contaminants is a complex task, we designed a pipeline combining six methods in a consensus strategy to assess the contamination level of 440 genome assemblies of Cyanobacteria. Two methods are based on published reference databases of ribosomal genes (SSU rRNA 16S and ribosomal proteins), one is indirectly based on a reference database of marker genes (CheckM), and three are based on complete genome analysis. Among those genome-wide methods, Kraken and DIAMOND blastx share the same reference database that we derived from Ensembl Bacteria, whereas CONCOCT does not require any reference database, instead relying on differences in DNA tetramer frequencies. Given that all the six methods appear to have their own strengths and limitations, we used the consensus of their rankings to infer that >5% of cyanobacterial genome assemblies are highly contaminated by foreign DNA (i.e., contaminants were detected by 5 or 6 methods). Our results will help researchers to check the quality of publicly available genomic data before use in their own analyses. Moreover, we argue that journals should make mandatory the submission of raw read data along with genome assemblies in order to facilitate the detection of contaminants in sequence databases.

Adaptation of a Chytrid Parasite to Its Cyanobacterial Host Is Hampered by Host Intraspecific Diversity

Experimental evolution can be used to test for and characterize parasite and pathogen adaptation. We undertook a serial-passage experiment in which a single parasite population of the obligate fungal (chytrid) parasite Rhizophydium megarrhizum was maintained over a period of 200 days under different mono- and multiclonal compositions of its phytoplankton host, the bloom-forming cyanobacterium Planktothrix. Despite initially inferior performance, parasite populations under sustained exposure to novel monoclonal hosts experienced rapid fitness increases evidenced by increased transmission rates. This demonstrates rapid adaptation of chytrids to novel hosts and highlights their high evolutionary potential. In contrast, increased fitness was not detected in parasites exposed to multiclonal host mixtures, indicating that cyanobacterial intraspecific diversity hampers parasites adaptation. Significant increases in intensity of infection were observed in monoclonal and multiclonal treatments, suggesting high evolvability of traits involved in parasite attachment onto hosts (i.e., encystment). A comparison of the performance of evolved and unevolved (control) parasite populations against their common ancestral host did not reveal parasite attenuation. Our results exemplify the ability of chytrid parasites to adapt rapidly to new hosts, while providing experimental evidence that genetic diversity in host populations grants increased resistance to disease by hindering parasite adaptation.

Monitoring and research of microcystins and environmental factors in a typical artificial freshwater aquaculture pond

Freshwater aquaculture ponds are important artificially regulated aquatic ecosystems which provide a large number of freshwater fish products in China. The cyanobacteria bloom and microcystin (MC) pollution caused by anthropogenic eutrophication have attracted much attention due to their toxic effects. To provide an insight into the cyanobacterial problem in the ponds, the environmental parameters and MCs of a typical artificial pond in the Yangtze River Delta region of China were monitored and studied from May to December 2015. During the monitoring period, the ponds were in serious eutrophication with total phosphorus (TP) concentrations between 0.95 and 1.80 μg/L, and total nitrogen (TN) concentrations between 1.1 and 4.86 μg/L. High feed coefficient and high fish stock were the main reasons for the eutrophication. The results showed that the water temperature was the key factor that affected the cyanobacteria blooming in the pond. The chlorophyll a concentration was significantly positively correlated with the cyanobacteria density during the blooming season. MC-LR and MC-RR existed simultaneously and showed a significant positive correlation. The peak concentrations of dissolved MC-LR and MC-RR in the pond water were 40.6 and 4.7 μg/L, respectively, which is considered highly toxic. Free MC-LR and MC-RR were also found in the aquaculture products. MC-LR concentrations in the bighead carp (Aristichthys nobilis) liver and shrimp (Macrobrachium nipponense) muscle were up to 2.64 and 4.17 μg/kg, respectively. MC-RR concentration was up to 1.89 μg/kg in the black carp (Mylopharyngodon piceus) liver. The results implied the potential health risks for citizens and pets caused by current artificial freshwater aquaculture pond systems.

Biochemical Dissection of the Natural Diversification of Microcystin Provides Lessons for Synthetic Biology of NRPS

The cyanobacterial hepatotoxin microcystin is assembled at a non-ribosomal peptide synthetase (NRPS) complex. The enormous structural diversity of this peptide, which is also found in closely related strains, is the result of frequent recombination events and point mutations. Here, we have compared the in vitro activation profiles of related monospecific and multispecific modules that either strictly incorporate leucine or arginine or incorporate chemically diverse amino acids in parallel into microcystin. By analyzing di- and tri-domain proteins we have dissected the role of adenylation and condensation domains for substrate specificity. We have further analyzed the role of subdomains and provide evidence for an extended gatekeeping function for the condensation domains of multispecific modules. By reproducing natural point mutations, we could convert a monospecific module into a multispecific module. Our findings may inspire novel synthetic biology approaches and demonstrate how recombination platforms of NRPSs have developed in nature.

Chemotype diversity in Planktothrix rubescens (cyanobacteria) populations is correlated to lake depth

The cyanobacterial species Planktothrix rubescens is known to preferably inhabit deep, stratified, oligo- to mesotrophic lakes. It is also known for the production of diverse bioactive peptides, including the hepatotoxic microcystins. A number of studies showed that P. rubescens populations generally consist of multiple distinct genotypes or chemotypes, respectively. In the present study, variability of chemotype diversity was analysed. Filaments of P. rubescens were isolated from water samples originating from 10 European lakes and analysed by MALDI-TOF mass spectrometry. In most of the analysed filaments multiple peptides belonging to multiple peptide classes could be detected. A resulting data matrix of 964 filaments and 37 individual peptides was subjected to correspondence analysis and K-means clustering. From the latter analysis the distribution of chemotypes among the lakes was established and diversity estimated by computing Shannon-Indices. Diversity varied strongly among lakes with the lowest diversity found in non-alpine lakes. Further, chemotype diversity was strongly correlated to the maximum depth of the sampled lakes in alpine and non-alpine lakes. The possible influence of both factors, geographic isolation and water column depth, on the observed patterns of chemotype diversity of P. rubescens populations is discussed.

Chemical and Genetic Diversity of Nodularia spumigena from the Baltic Sea

Nodularia spumigena is a toxic, filamentous cyanobacterium occurring in brackish waters worldwide, yet forms extensive recurrent blooms in the Baltic Sea. N. spumigena produces several classes of non-ribosomal peptides (NRPs) that are active against several key metabolic enzymes. Previously, strains from geographically distant regions showed distinct NRP metabolic profiles. In this work, conspecific diversity in N. spumigena was studied using chemical and genetic approaches. NRP profiles were determined in 25 N. spumigena strains isolated in different years and from different locations in the Baltic Sea using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Genetic diversity was assessed by targeting the phycocyanin intergenic spacer and flanking regions (cpcBA-IGS). Overall, 14 spumigins, 5 aeruginosins, 2 pseudaeruginosins, 2 nodularins, 36 anabaenopeptins, and one new cyanopeptolin-like peptide were identified among the strains. Seven anabaenopeptins were new structures; one cyanopeptolin-like peptide was discovered in N. spumigena for the first time. Based on NRP profiles and cpcBA-IGS sequences, the strains were grouped into two main clusters without apparent influence of year and location, indicating persistent presence of these two subpopulations in the Baltic Sea. This study is a major step in using chemical profiling to explore conspecific diversity with a higher resolution than with a sole genetic approach.

Is qPCR a Reliable Indicator of Cyanotoxin Risk in Freshwater?

The wide distribution of cyanobacteria in aquatic environments leads to the risk of water contamination by cyanotoxins, which generate environmental and public health issues. Measurements of cell densities or pigment contents allow both the early detection of cellular growth and bloom monitoring, but these methods are not sufficiently accurate to predict actual cyanobacterial risk. To quantify cyanotoxins, analytical methods are considered the gold standards, but they are laborious, expensive, time-consuming and available in a limited number of laboratories. In cyanobacterial species with toxic potential, cyanotoxin production is restricted to some strains, and blooms can contain varying proportions of both toxic and non-toxic cells, which are morphologically indistinguishable. The sequencing of cyanobacterial genomes led to the description of gene clusters responsible for cyanotoxin production, which paved the way for the use of these genes as targets for PCR and then quantitative PCR (qPCR). Thus, the quantification of cyanotoxin genes appeared as a new method for estimating the potential toxicity of blooms. This raises a question concerning whether qPCR-based methods would be a reliable indicator of toxin concentration in the environment. Here, we review studies that report the parallel detection of microcystin genes and microcystin concentrations in natural populations and also a smaller number of studies dedicated to cylindrospermopsin and saxitoxin. We discuss the possible issues associated with the contradictory findings reported to date, present methodological limitations and consider the use of qPCR as an indicator of cyanotoxin risk.

Combined LC-MS/MS and Molecular Networking Approach Reveals New Cyanotoxins from the 2014 Cyanobacterial Bloom in Green Lake, Seattle

Cyanotoxins obtained from a freshwater cyanobacterial collection at Green Lake, Seattle during a cyanobacterial harmful algal bloom in the summer of 2014 were studied using a new approach based on molecular networking analysis of liquid chromatography tandem mass spectrometry (LC-MS/MS) data. This MS networking approach is particularly well-suited for the detection of new cyanotoxin variants and resulted in the discovery of three new cyclic peptides, namely microcystin-MhtyR (6), which comprised about half of the total microcystin content in the bloom, and ferintoic acids C (12) and D (13). Structure elucidation of 6 was aided by a new microscale methylation procedure. Metagenomic analysis of the bloom using the 16S-ITS rRNA region identified Microcystis aeruginosa as the predominant cyanobacterium in the sample. Fragments of the putative biosynthetic genes for the new cyanotoxins were also identified, and their sequences correlated to the structure of the isolated cyanotoxins.

Cyanobacterial Toxic and Bioactive Peptides in Freshwater Bodies of Greece: Concentrations, Occurrence Patterns, and Implications for Human Health

Cyanobacterial harmful algal blooms represent one of the most conspicuous waterborne microbial hazards in aquatic environments mostly due to the production of toxic secondary metabolites, mainly microcystins (MCs). Other bioactive peptides are frequently found in cyanobacterial blooms, yet their concentration and ecological relevance is still unknown. In this paper we studied the presence and concentration of cyanobacterial peptides (microcystins, anabaenopeptins, anabaenopeptilides) in 36 Greek freshwater bodies, using HPLC-DAD, ELISA, and PP1IA. Microcystins were found in more than 90% of the samples investigated, indicating that microcystin-producing strains seem to also occur in lakes without blooms. Microcystins MC-RR, MC-LR, and MC-YR were the main toxin constituents of the bloom samples. Anabaenopeptin A and B were predominant in some samples, whereas anabaenopeptolide 90A was the only peptide found in Lake Mikri Prespa. The intracellular concentrations of anabaenopeptins produced by cyanobacterial bloom populations are determined for the first time in this study; the high (>1000 µg·L(-1)) anabaenopeptin concentration found indicates there may be some impacts, at least on the ecology and the food web structure of the aquatic ecosystems. The maximum intracellular MC values measured in Lakes Kastoria and Pamvotis, exceeding 10,000 µg·L(-1), are among the highest reported.

The genetic basis for O-acetylation of the microcystin toxin in cyanobacteria

Microcystins are a family of cyclic peptide toxins produced by cyanobacteria. They are responsible for the toxicosis and death of wild and domestic animals throughout the world. They display extensive variation in amino acid composition and functional group chemistry. O-acetylated microcystins are frequently produced by free-living and symbiotic strains of the genus Nostoc. Here, we show that the production of acetylated microcystins is catalyzed by an acetyl-coenzyme A-dependent O-acetyltransferase (McyL) encoded in the 57 kb microcystin synthetase gene cluster of Nostoc sp. 152. Phylogenetic analysis demonstrates that McyL belongs to a family of enzymes that inactivate antibiotics through O-acetylation. The McyL enzyme has a relaxed substrate specificity, allowing the preparation of semisynthetic microcystins. This study sheds light on the evolutionary origins and genetic diversity of an important class of enzymes involved in antibiotic resistance.

Oligopeptides as biomarkers of cyanobacterial subpopulations. Toward an understanding of their biological role

Cyanobacterial oligopeptides comprise a wide range of bioactive and/or toxic compounds. While current research is strongly focused on exploring new oligopeptide variants and their bioactive properties, the biological role of these compounds remains elusive. Oligopeptides production abilities show a remarkably patchy distribution among conspecific strains. This observation has prompted alternative approaches to unveil their adaptive value, based on the use of cellular oligopeptide compositions as biomarkers of intraspecific subpopulations or chemotypes in freshwater cyanobacteria. Studies addressing the diversity, distribution, and dynamics of chemotypes in natural systems have provided important insights into the structure and ecology of cyanobacterial populations and the adaptive value of oligopeptides. This review presents an overview of the fundamentals of this emerging approach and its most relevant findings, and discusses our current understanding of the role of oligopeptides in the ecology of cyanobacteria.

Congruence between mcy based genetic type and microcystin composition within the populations of toxic Microcystis in a Plateau Lake, China

Genetic diversity and differential microcystin contributions within the populations of toxic Microcystis in freshwater ecosystems have not been fully discussed. To illustrate this issue, we sequenced clone libraries of the first adenylation domain of mcyB (mcyBA1) from Microcystis populations in a plateau lake. Phylogenetic analysis revealed two divergent groups of McyBA1 with two closely related subgroups within Group I. While neutral evolution was indicated on the whole McyBA1 domain, 1-5 recombination breakpoints and several codons under positive or negative selections were found. Significant seasonal changes of McyBA1 diversity were observed. Group I composed the major proportion of the McyBA1 pool throughout the growing season in the whole lake. Group IA and Group II denominated Microcystis strains isolated from this lake were characterized by preference production of microcystin-RR (62-85%) and microcystin-LR (> 98%) respectively. We detected the intracellular microcystins in lake water and microcystin-RR was a main variant (mostly > 50%). In summary, McyBA1 subgroups were dominant within the population of toxic Microcystis and contributed the predominance of microcystin-RR in the lake. The differences of substrate preference in microcystin biosynthesis among groups were caused by neutral evolution and homologous recombination.

From green to red: horizontal gene transfer of the phycoerythrin gene cluster between Planktothrix strains

Horizontal gene transfer is common in cyanobacteria, and transfer of large gene clusters may lead to acquisition of new functions and conceivably niche adaption. In the present study, we demonstrate that horizontal gene transfer between closely related Planktothrix strains can explain the production of the same oligopeptide isoforms by strains of different colors. Comparison of the genomes of eight Planktothrix strains revealed that strains producing the same oligopeptide isoforms are closely related, regardless of color. We have investigated genes involved in the synthesis of the photosynthetic pigments phycocyanin and phycoerythrin, which are responsible for green and red appearance, respectively. Sequence comparisons suggest the transfer of a functional phycoerythrin gene cluster generating a red phenotype in a strain that is otherwise more closely related to green strains. Our data show that the insertion of a DNA fragment containing the 19.7-kb phycoerythrin gene cluster has been facilitated by homologous recombination, also replacing a region of the phycocyanin operon. These findings demonstrate that large DNA fragments spanning entire functional gene clusters can be effectively transferred between closely related cyanobacterial strains and result in a changed phenotype. Further, the results shed new light on the discussion of the role of horizontal gene transfer in the sporadic distribution of large gene clusters in cyanobacteria, as well as the appearance of red and green strains.

Harmful cyanobacterial blooms: causes, consequences, and controls

Cyanobacteria are the Earth's oldest oxygenic photoautotrophs and have had major impacts on shaping its biosphere. Their long evolutionary history (≈ 3.5 by) has enabled them to adapt to geochemical and climatic changes, and more recently anthropogenic modifications of aquatic environments, including nutrient over-enrichment (eutrophication), water diversions, withdrawals, and salinization. Many cyanobacterial genera exhibit optimal growth rates and bloom potentials at relatively high water temperatures; hence global warming plays a key role in their expansion and persistence. Bloom-forming cyanobacterial taxa can be harmful from environmental, organismal, and human health perspectives by outcompeting beneficial phytoplankton, depleting oxygen upon bloom senescence, and producing a variety of toxic secondary metabolites (e.g., cyanotoxins). How environmental factors impact cyanotoxin production is the subject of ongoing research, but nutrient (N, P and trace metals) supply rates, light, temperature, oxidative stressors, interactions with other biota (bacteria, viruses and animal grazers), and most likely, the combined effects of these factors are all involved. Accordingly, strategies aimed at controlling and mitigating harmful blooms have focused on manipulating these dynamic factors. The applicability and feasibility of various controls and management approaches is discussed for natural waters and drinking water supplies. Strategies based on physical, chemical, and biological manipulations of specific factors show promise; however, a key underlying approach that should be considered in almost all instances is nutrient (both N and P) input reductions; which have been shown to effectively reduce cyanobacterial biomass, and therefore limit health risks and frequencies of hypoxic events.

Convergent evolution of [D-Leucine(1)] microcystin-LR in taxonomically disparate cyanobacteria

Background

Many important toxins and antibiotics are produced by non-ribosomal biosynthetic pathways. Microcystins are a chemically diverse family of potent peptide toxins and the end-products of a hybrid NRPS and PKS secondary metabolic pathway. They are produced by a variety of cyanobacteria and are responsible for the poisoning of humans as well as the deaths of wild and domestic animals around the world. The chemical diversity of the microcystin family is attributed to a number of genetic events that have resulted in the diversification of the pathway for microcystin assembly.

Results

Here, we show that independent evolutionary events affecting the substrate specificity of the microcystin biosynthetic pathway have resulted in convergence on a rare [D-Leu¹] microcystin-LR chemical variant. We detected this rare microcystin variant from strains of the distantly related genera Microcystis, Nostoc, and Phormidium. Phylogenetic analysis performed using sequences of the catalytic domains within the mcy gene cluster demonstrated a clear recombination pattern in the adenylation domain phylogenetic tree. We found evidence for conversion of the gene encoding the McyA₂ adenylation domain in strains of the genera Nostoc and Phormidium. However, point mutations affecting the substrate-binding sequence motifs of the McyA₂ adenylation domain were associated with the change in substrate specificity in two strains of Microcystis. In addition to the main [D-Leu¹] microcystin-LR variant, these two strains produced a new microcystin that was identified as [Met¹] microcystin-LR.

Conclusions

Phylogenetic analysis demonstrated that both point mutations and gene conversion result in functional mcy gene clusters that produce the same rare [D-Leu¹] variant of microcystin in strains of the genera Microcystis, Nostoc, and Phormidium. Engineering pathways to produce recombinant non-ribosomal peptides could provide new natural products or increase the activity of known compounds. Our results suggest that the replacement of entire adenylation domains could be a more successful strategy to obtain higher specificity in the modification of the non-ribosomal peptides than point mutations.

Gene flow, recombination, and selection in cyanobacteria: population structure of geographically related Planktothrix freshwater strains

Several Planktothrix strains, each producing a distinct oligopeptide profile, have been shown to coexist within Lake Steinsfjorden (Norway). Using nonribosomal peptide synthetase (NRPS) genes as markers, it has been shown that the Planktothrix community comprises distinct genetic variants displaying differences in bloom dynamics, suggesting a Planktothrix subpopulation structure. Here, we investigate the Planktothrix variants inhabiting four lakes in southeast of Norway utilizing both NRPS and non-NRPS genes. Phylogenetic analyses showed similar topologies for both NRPS and non-NRPS genes, and the lakes appear to have similar structuring of Planktothrix genetic variants. The structure of distinct variants was also supported by very low genetic diversity within variants compared to the between-variant diversity. Incongruent topologies and split decomposition revealed recombination events between Planktothrix variants. In several strains the gene variants seem to be a result of recombination. Both NRPS and non-NRPS genes are dominated by purifying selection; however, sites subjected to positive selection were also detected. The presence of similar and well-separated Planktothrix variants with low internal genetic diversity indicates gene flow within Planktothrix populations. Further, the low genetic diversity found between lakes (similar range as within lakes) indicates gene flow also between Planktothrix populations and suggests recent, or recurrent, dispersals. Our data also indicate that recombination has resulted in new genetic variants. Stability within variants and the development of new variants are likely to be influenced by selection patterns and within-variant homologous recombination.

Nodularin, a cyanobacterial toxin, is synthesized in planta by symbiotic Nostoc sp

The nitrogen-fixing bacterium, Nostoc, is a commonly occurring cyanobacterium often found in symbiotic associations. We investigated the potential of cycad cyanobacterial endosymbionts to synthesize microcystin/nodularin. Endosymbiont DNA was screened for the aminotransferase domain of the toxin biosynthesis gene clusters. Five endosymbionts carrying the gene were screened for bioactivity. Extracts of two isolates inhibited protein phosphatase 2A and were further analyzed using electrospray ionization mass spectrometry (ESI-MS)/MS. Nostoc sp. 'Macrozamia riedlei 65.1' and Nostoc sp. 'Macrozamia serpentina 73.1' both contained nodularin. High performance liquid chromatography (HPLC) HESI-MS/MS analysis confirmed the presence of nodularin at 9.55±2.4 ng μg-1 chlorophyll a in Nostoc sp. 'Macrozamia riedlei 65.1' and 12.5±8.4 ng μg-1 Chl a in Nostoc sp. 'Macrozamia serpentina 73.1' extracts. Further scans indicated the presence of the rare isoform [L-Har(2)] nodularin, which contains L-homoarginine instead of L-arginine. Nodularin was also present at 1.34±0.74 ng ml(-1) (approximately 3 pmol per g plant ww) in the methanol root extracts of M. riedlei MZ65, while the presence of [L-Har(2)] nodularin in the roots of M. serpentina MZ73 was suggested by HPLC HESI-MS/MS analysis. The ndaA-B and ndaF genomic regions were sequenced to confirm the presence of the hybrid polyketide/non-ribosomal gene cluster. A seven amino-acid insertion into the NdaA-C1 domain of N. spumigena NSOR10 protein was observed in all endosymbiont-derived sequences, suggesting the transfer of the nda cluster from N. spumigena to terrestrial Nostoc species. This study demonstrates the synthesis of nodularin and [L-Har(2)] nodularin in a non-Nodularia species and the production of cyanobacterial hepatotoxin by a symbiont in planta.

The chemical ecology of cyanobacteria

This review covers the literature on the chemically mediated ecology of cyanobacteria, including ultraviolet radiation protection, feeding-deterrence, allelopathy, resource competition, and signalling. To highlight the chemical and biological diversity of this group of organisms, evolutionary and chemotaxonomical studies are presented. Several technologically relevant aspects of cyanobacterial chemical ecology are also discussed.

Phylogenetic inference of colony isolates comprising seasonal Microcystis blooms in Lake Taihu, China

Blooms of the toxin-producing cyanobacterium, Microcystis spp., are an increasingly prevalent water quality problem and health hazard worldwide. China's third largest lake, Lake Taihu, has been experiencing progressively more severe Microcystis blooms over the past three decades. In 2009 and 2010, individual Microcystis colonies, consisting of four different morphospecies, were isolated and genotyped using a whole-cell multiplex PCR assay. The 16S-23S rDNA-ITS sequences were aligned based on Bayesian inference and indicated that one morphospecies was genetically unique (Microcystis wesenbergii) and three were indistinguishable (Microcystis aeruginosa, Microcystis flos-aquae, and Microcystis ichthyoblabe). Microcystin (mcyB) genes were detected intermittently in two of the morphospecies while the other two morphospecies lacked the mcyB gene in all samples. Water temperature was found to influence bloom formation and morphotype prevalence, and chlorophyll a and temperature were positively and significantly correlated with microcystin concentration. Cooler water temperatures promoted toxigenic strains of Microcystis. Wind appeared to influence the distribution of morphotypes across the lake, with M. aeruginosa and M. ichthyoblabe being more susceptible to wind stress than M. wesenbergii and M. flos-aquae. The results of this study indicated that the blooms were composed of a variety of Microcystis morphospecies, with more genotypes observed than can be attributed to individual morphotypes. We conclude that morphology is not a reliable indicator of toxigenicity in Lake Taihu, and caution should be exercised when the M. aeruginosa morphotype is present because it is capable of producing MC-LR, the most toxic microcystin isoform.

Genetic variation of adenylation domains of the anabaenopeptin synthesis operon and evolution of substrate promiscuity

Anabaenopeptins (AP) are bioactive cyclic hexapeptides synthesized nonribosomally in cyanobacteria. APs are characterized by several conserved motifs, including the ureido bond, N-methylation in position 5, and d-Lys in position 2. All other positions of the AP molecule are variable, resulting in numerous structural variants. We have identified a nonribosomal peptide synthetase (NRPS) operon from Planktothrix agardhii strain CYA126/8 consisting of five genes (apnA to apnE) encoding six NRPS modules and have confirmed its role in AP synthesis by the generation of a mutant via insertional inactivation of apnC. In order to correlate the genetic diversity among adenylation domains (A domains) with AP structure variation, we sequenced the A domains of all six NRPS modules from seven Planktothrix strains differing in the production of AP congeners. It is remarkable that single strains coproduce APs bearing either of the chemically divergent amino acids Arg and Tyr in exocyclic position 1. Since the A domain of the initiation module (the ApnA A₁ domain) has been proposed to activate the amino acid incorporated into exocyclic position 1, we decided to analyze this domain both biochemically and phylogenetically. Only ApnA A₁ enzymes from strains producing AP molecules containing Arg or Tyr in position 1 were found to activate these two chemically divergent amino acids in vitro. Phylogenetic analysis of apn A domain sequences revealed that strains with a promiscuous ApnA A₁ domain are derived from an ancestor that activates only Arg. Surprisingly, positive selection appears to affect only three codons within the apnA A₁ gene, suggesting that this remarkable promiscuity has evolved from point mutations only.

Diversity of and selection acting on cylindrospermopsin cyrB gene adenylation domain sequences in Florida

Aphanizomenon ovalisporum is the only confirmed cylindrospermopsin producer identified in the United States to date. On the other hand, Cylindrospermopsis raciborskii is a prominent feature of many lakes in Florida and other regions of the United States. To see the variation in cylindrospermopsin cyrB gene adenylation domain sequences and possibly discover new cylindrospermopsin producers, we collected water samples for a 3-year period from 17 different systems in Florida. Positive amplicons were cloned and sequenced, revealing that approximately 92% of sequences were A. ovalisporum-like (>99% identity). Interestingly, 6% of sequences were very similar (>99% identity) to cyrB sequences of C. raciborskii from Australia and of Aphanizomenon sp. from Germany. Neutrality tests suggest that A. ovalisporum-like cyrB adenylation domain sequences are under purifying selection, with abundant low-frequency polymorphisms within the population. On the other hand, when compared between species by codon-based methods, amino acids of CyrB also seem to be under purifying selection, in accordance with the one proposed amino acid thought to be activated by the CyrB adenylation domain.

Phylogenetic methods in natural product research

Natural products researchers are increasingly employing evolutionary analyses of genes and gene products that rely on phylogenetic trees. The field of phylogenetic inference and of evolutionary analyses based on phylogenies is growing at an amazing rate, making it difficult to keep up with the latest methodologies. Here, we summarize phylogenetic applications in natural products research, and review methods and software useful for carrying out analyses inferring or using phylogenetic trees. We include an updated overview of available alignment methods and programs, as well as a selection of some useful phylogenetic analysis tools. This review covers primarily the period 2000-2009 for applications of phylogenetic methods in natural product research, and 1990-2009 for phylogenetic methods, with some references going back to the 1960s.

Recombination, cryptic clades and neutral molecular divergence of the microcystin synthetase (mcy) genes of toxic cyanobacterium Microcystis aeruginosa

Background

The water-bloom-forming cyanobacterium Microcystis aeruginosa is a known producer of various kinds of toxic and bioactive chemicals. Of these, hepatotoxic cyclic heptapeptides microcystins have been studied most intensively due to increasing concerns for human health risks and environmental damage. More than 70 variants of microcystins are known, and a single microcystin synthetase (mcy) gene cluster consisting of 10 genes (mcyA to mcyJ) has been identified to be responsible for the production of all known variants of microcystins. Our previous multilocus sequence typing (MLST) analysis of the seven housekeeping genes indicated that microcystin-producing strains of M. aeruginosa are classified into two phylogenetic groups.

Results

To investigate whether the mcy genes are genetically structured similarly as in MLST analysis of the housekeeping genes and to identify the evolutionary forces responsible for the genetic divergence of these genes, we used 118 mcy-positive isolates to perform phylogenetic and population genetic analyses of mcy genes based on three mcy loci within the mcy gene cluster (mcyD, mcyG, and mcyJ), none of which is involved in the production of different microcystin variants. Both individual phylogenetic analysis and multilocus genealogical analysis of the mcy genes divided our isolates into two clades, consistent with the MLST phylogeny based on seven housekeeping loci. No shared characteristics within each clade are known, and microcystin analyses did not identify any compositional trend specific to each clade. Statistical analyses for recombination indicated that recombination among the mcy genes is much more frequent within clades than between, suggesting that recombination has been an important force maintaining the cryptic divergence of mcy genes. On the other hand, a series of statistical tests provided no strong evidence for selection to explain the deep divergence of the mcy genes. Furthermore, analysis of molecular variance (AMOVA) indicated a low level of geographic structuring in the genetic diversity of mcy.

Conclusion

Our phylogenetic analyses suggest that the mcy genes of M. aeruginosa are subdivided into two cryptic clades, consistent with the phylogeny determined by MLST. Population genetic analyses suggest that these two clades have primarily been maintained as a result of homology-dependent recombination and neutral genetic drift.