Identification of promoter elements in the Dolichospermum circinale AWQC131C saxitoxin gene cluster and the experimental analysis of their use for heterologous expression.

BACKGROUND: Dolichospermum circinale is a filamentous bloom-forming cyanobacterium responsible for biosynthesis of the paralytic shellfish toxins (PST), including saxitoxin. PSTs are neurotoxins and in their purified form are important analytical standards for monitoring the quality of water and seafood and biomedical research tools for studying neuronal sodium channels. More recently, PSTs have been recognised for their utility as local anaesthetics. Characterisation of the transcriptional elements within the saxitoxin (sxt) biosynthetic gene cluster (BGC) is a first step towards accessing these molecules for biotechnology. RESULTS: In D. circinale AWQC131C the sxt BGC is transcribed from two bidirectional promoter regions encoding five individual promoters. These promoters were identified experimentally using 5' RACE and their activity assessed via coupling to a lux reporter system in E. coli and Synechocystis sp. PCC 6803. Transcription of the predicted drug/metabolite transporter (DMT) encoded by sxtPER was found to initiate from two promoters, PsxtPER1 and PsxtPER2. In E. coli, strong expression of lux from PsxtP, PsxtD and PsxtPER1 was observed while expression from Porf24 and PsxtPER2 was remarkably weaker. In contrast, heterologous expression in Synechocystis sp. PCC 6803 showed that expression of lux from PsxtP, PsxtPER1, and Porf24 promoters was statistically higher compared to the non-promoter control, while PsxtD showed poor activity under the described conditions. CONCLUSIONS: Both of the heterologous hosts investigated in this study exhibited high expression levels from three of the five sxt promoters. These results indicate that the majority of the native sxt promoters appear active in different heterologous hosts, simplifying initial cloning efforts. Therefore, heterologous expression of the sxt BGC in either E. coli or Synechocystis could be a viable first option for producing PSTs for industrial or biomedical purposes.

Proteogenomics of a saxitoxin-producing and non-toxic strain of Anabaena circinalis (cyanobacteria) in response to extracellular NaCl and phosphate depletion.

In Australia, saxitoxin production is strain dependent within the bloom-forming freshwater cyanobacterium Anabaena circinalis. Freshwater cyanobacteria are exposed to rapid fluctuations in environmental nutrient concentrations, and their adaption is vital for competition, succession and dominance. Two elements of environmental significance, phosphorus and sodium chloride, are proposed to play a role in bloom development and saxitoxin biosynthesis respectively. The aim of our study was to comparatively analyse the model saxitoxin-producing A. circinalis AWQC131C and non-toxic A. circinalis AWQC310F at the genomic level and proteomic level, in response to phosphate depletion and increased extracellular NaCl. When challenged, photosynthesis, carbon/nitrogen metabolisms, transcription/translation, oxidative stress and nutrient transport functional categories demonstrated the largest changes in protein abundance. In response to increased NaCl, SxtC, a protein conserved in all known saxitoxin biosynthetic pathways, was downregulated. Additionally, toxin quantification revealed a decrease in total saxitoxin and decarbomoyl-gonyautoxin2/3 content in response to the NaCl treatment. In response to phosphate depletion, the toxic and non-toxic strain displayed similar proteomic profiles, although the toxic strain did not alter the abundance of as many proteins as the non-toxic strain. These findings have important implications for the future, since response and adaption mechanisms are directly related to in situ dominance of cyanobacteria.

Genome mining for natural product biosynthetic gene clusters in the Subsection V cyanobacteria.

BACKGROUND: Cyanobacteria are well known for the production of a range of secondary metabolites. Whilst recent genome sequencing projects has led to an increase in the number of publically available cyanobacterial genomes, the secondary metabolite potential of many of these organisms remains elusive. Our study focused on the 11 publically available Subsection V cyanobacterial genomes, together with the draft genomes of Westiella intricata UH strain HT-29-1 and Hapalosiphon welwitschii UH strain IC-52-3, for their genetic potential to produce secondary metabolites. The Subsection V cyanobacterial genomes analysed in this study are reported to produce a diverse range of natural products, including the hapalindole-family of compounds, microcystin, hapalosin, mycosporine-like amino acids and hydrocarbons. RESULTS: A putative gene cluster for the cyclic depsipeptide hapalosin, known to reverse P-glycoprotein multiple drug resistance, was identified within three Subsection V cyanobacterial genomes, including the producing cyanobacterium H. welwitschii UH strain IC-52-3. A number of orphan NRPS/PKS gene clusters and ribosomally-synthesised and post translationally-modified peptide gene clusters (including cyanobactin, microviridin and bacteriocin gene clusters) were identified. Furthermore, gene clusters encoding the biosynthesis of mycosporine-like amino acids, scytonemin, hydrocarbons and terpenes were also identified and compared. CONCLUSIONS: Genome mining has revealed the diversity, abundance and complex nature of the secondary metabolite potential of the Subsection V cyanobacteria. This bioinformatic study has identified novel biosynthetic enzymes which have not been associated with gene clusters of known classes of natural products, suggesting that these cyanobacteria potentially produce structurally novel secondary metabolites.

Comparative proteomics reveals that a saxitoxin-producing and a nontoxic strain of Anabaena circinalis are two different ecotypes.

In Australia, saxitoxin production is restricted to the cyanobacterial species Anabaena circinalis and is strain-dependent. We aimed to characterize a saxitoxin-producing and nontoxic strain of A. circinalis at the proteomic level using iTRAQ. Seven proteins putatively involved in saxitoxin biosynthesis were identified within our iTRAQ experiment for the first time. The proteomic profile of the toxic A. circinalis was significantly different from the nontoxic strain, indicating that each is likely to inhabit a unique ecological niche. Under control growth conditions, the saxitoxin-producing A. circinalis displayed a higher abundance of photosynthetic, carbon fixation and nitrogen metabolic proteins. Differential abundance of these proteins suggests a higher intracellular C:N ratio and a higher concentration of intracellular 2-oxoglutarate in our toxic strain compared with the nontoxic strain. This may be a novel site for posttranslational regulation because saxitoxin biosynthesis putatively requires a 2-oxoglutarate-dependent dioxygenase. The nontoxic A. circinalis was more abundant in proteins, indicating cellular stress. Overall, our study has provided the first insight into fundamental differences between a toxic and nontoxic strain of A. circinalis, indicating that they are distinct ecotypes.

Comparative analysis of hapalindole, ambiguine and welwitindolinone gene clusters and reconstitution of indole-isonitrile biosynthesis from cyanobacteria.

BACKGROUND: The hapalindole-type family of natural products is a group of hybrid isoprenoid-indole alkaloids, produced solely by members of the Subsection V cyanobacterial strains. This family broadly includes the hapalindoles, welwitindolinones, fisherindoles and ambiguines amongst others, all of which have an isonitrile- or isothiocyanate-containing indole alkaloid skeleton, with a cyclized isoprene unit. The hapalindoles are diversified into the welwitindolinones, fischerindoles and ambiguines through the employment of tailoring oxygenase, methyltransferase and prenyltransferase enzymes. We compare the genetic basis for the biosynthesis of this diverse group of natural products and identify key early biosynthetic intermediates. RESULTS: Whole genome sequencing of freshwater and terrestrial cyanobacteria Westiella intricata UH strain HT-29-1, Hapalosiphon welwitschii UH strain IC-52-3, Fischerella ambigua UTEX 1903 and Fischerella sp. ATCC 43239 led to the identification of a candidate hapalindole-type gene cluster in each strain. These were compared with the recently published ambiguine and welwitindolinone gene clusters and four unpublished clusters identified within publicly available genomes. We present detailed comparative bioinformatic analysis of the gene clusters and the biosynthesis of a pivotal indole-isonitrile intermediate resulting in both cis and trans geometrical isomers. Enzyme analyses and metabolite extractions from two hapalindole-producing Fischerella strains indicate the presence of cis and trans indole-isonitriles as biosynthetic intermediates in the early steps of the pathway. CONCLUSIONS: Interestingly, the organization of the welwitindolinone gene cluster is conserved in all producing strains but distinct from the hapalindole and ambiguine clusters. Enzymatic assays using WelI1 and WelI3 from Westiella intricata UH strain HT-29-1 demonstrated the ability to catalyze the formation of both cis and trans geometrical isomers when using a cell lysate. The enzymatic and metabolic characterization of both cis and trans indole-isonitrile intermediates implies conservation of their stereochemical integrity towards members of the ambiguine and welwitindolinone products. In summary, we present data that supports a unified biosynthetic pathway towards hapalindoles in nine individual species of cyanobacteria. Diversification of the pathway occurs later through the employment of specialized enzymatic steps towards fischerindoles, ambiguines and welwitindolinones.

Environmental conditions that influence toxin biosynthesis in cyanobacteria.

Over the past 15 years, the genetic basis for production of many cyanobacterial bioactive compounds has been described. This knowledge has enabled investigations into the environmental factors that regulate the production of these toxins at the molecular level. Such molecular or systems level studies are also likely to reveal the physiological role of the toxin and contribute to effective water resource management. This review focuses on the environmental regulation of some of the most relevant cyanotoxins, namely the microcystins, nodularin, cylindrospermopsin, saxitoxins, anatoxins and jamaicamides.

Both Constitutive and Infection-Responsive Secondary Metabolites Linked to Resistance against Austropuccinia psidii (Myrtle Rust) in Melaleuca quinquenervia.

Austropuccinia psidii is a fungal plant pathogen that infects species within the Myrtaceae, causing the disease myrtle rust. Myrtle rust is causing declines in populations within natural and managed ecosystems and is expected to result in species extinctions. Despite this, variation in response to A. psidii exist within some species, from complete susceptibility to resistance that prevents or limits infection by the pathogen. Untargeted metabolomics using Ultra Performance Liquid Chromatography with Ion Mobility followed by analysis using MetaboAnalyst 3.0, was used to explore the chemical defence profiles of resistant, hypersensitive and susceptible phenotypes within Melaleuca quinquenervia during the early stages of A. psidii infection. We were able to identify three separate pools of secondary metabolites: (i) metabolites classified structurally as flavonoids that were naturally higher in the leaves of resistant individuals prior to infection, (ii) organoheterocyclic and carbohydrate-related metabolites that varied with the level of host resistance post-infection, and (iii) metabolites from the terpenoid pathways that were responsive to disease progression regardless of resistance phenotype suggesting that these play a minimal role in disease resistance during the early stages of colonization of this species. Based on the classes of these secondary metabolites, our results provide an improved understanding of key pathways that could be linked more generally to rust resistance with particular application within Melaleuca.

Neurotoxic alkaloids: saxitoxin and its analogs.

Saxitoxin (STX) and its 57 analogs are a broad group of natural neurotoxic alkaloids, commonly known as the paralytic shellfish toxins (PSTs). PSTs are the causative agents of paralytic shellfish poisoning (PSP) and are mostly associated with marine dinoflagellates (eukaryotes) and freshwater cyanobacteria (prokaryotes), which form extensive blooms around the world. PST producing dinoflagellates belong to the genera Alexandrium, Gymnodinium and Pyrodinium whilst production has been identified in several cyanobacterial genera including Anabaena, Cylindrospermopsis, Aphanizomenon Planktothrix and Lyngbya. STX and its analogs can be structurally classified into several classes such as non-sulfated, mono-sulfated, di-sulfated, decarbamoylated and the recently discovered hydrophobic analogs--each with varying levels of toxicity. Biotransformation of the PSTs into other PST analogs has been identified within marine invertebrates, humans and bacteria. An improved understanding of PST transformation into less toxic analogs and degradation, both chemically or enzymatically, will be important for the development of methods for the detoxification of contaminated water supplies and of shellfish destined for consumption. Some PSTs also have demonstrated pharmaceutical potential as a long-term anesthetic in the treatment of anal fissures and for chronic tension-type headache. The recent elucidation of the saxitoxin biosynthetic gene cluster in cyanobacteria and the identification of new PST analogs will present opportunities to further explore the pharmaceutical potential of these intriguing alkaloids.

Draft Genome Sequence of the Fungus Lecanicillium psalliotae Strain HWLR35, Isolated from a Wheat Leaf Infected with Leaf Rust (Caused by Puccinia triticina).

Lecanicillium psalliotae is an entomopathogenic, mycoparasitical, and nematophagous fungus known to produce antibiotic and antifungal compounds. Here, we report the first 36-Mb draft genome sequence of L. psalliotae strain HWLR35. The draft genome contains 197 scaffolds and is predicted to have 11,009 protein-coding genes.

Structural determinants and modulation of substrate specificity in phenylalanine-tyrosine ammonia-lyases.

Aromatic amino acid ammonia-lyases catalyze the deamination of L-His, L-Phe, and L-Tyr, yielding ammonia plus aryl acids bearing an alpha,beta-unsaturated propenoic acid. We report crystallographic analyses of unliganded Rhodobacter sphaeroides tyrosine ammonia-lyase (RsTAL) and RsTAL bound to p-coumarate and caffeate. His 89 of RsTAL forms a hydrogen bond with the p-hydroxyl moieties of coumarate and caffeate. His 89 is conserved in TALs but replaced in phenylalanine ammonia-lyases (PALs) and histidine ammonia-lyases (HALs). Substitution of His 89 by Phe, a characteristic residue of PALs, yields a mutant with a switch in kinetic preference from L-Tyr to L-Phe. Structures of the H89F mutant in complex with the PAL product, cinnamate, or the PAL-specific inhibitor, 2-aminoindan-2-phosphonate (AIP), support the role of position 89 as a specificity determinant in the family of aromatic amino acid ammonia-lyases and aminomutases responsible for beta-amino acid biosynthesis.

Advances in genomics, transcriptomics and proteomics of toxin-producing cyanobacteria.

A common misconception persists that the genomes of toxic and non-toxic cyanobacterial strains are largely conserved with the exception of the presence or absence of the genes responsible for toxin production. Implementation of -omics era technologies has challenged this paradigm, with comparative analyses providing increased insight into the differences between strains of the same species. The implementation of genomic, transcriptomic and proteomic approaches has revealed distinct profiles between toxin-producing and non-toxic strains. Further, metagenomics and metaproteomics highlight the genomic potential and functional state of toxic bloom events over time. In this review, we highlight how these technologies have shaped our understanding of the complex relationship between these molecules, their producers and the environment at large within which they persist.

Evidence for differences in the metabolism of saxitoxin and C1+2 toxins in the freshwater cyanobacterium Cylindrospermopsis raciborskii T3.

The activity of paralytic shellfish poisoning (PSP) toxins biosynthetic enzymes was assayed in the cyanobacterium Cylindrospermopsis raciborskii T3 after inhibiting protein synthesis with chloramphenicol (CAM). The production of C1+2 and saxitoxin (STX) was sensitive to CAM with STX levels decreasing by 70% after 24-h exposure to the antibiotic. PSP toxin production was strongly promoted by arginine supplementation, with a maximum 476% increase in intracellular STX concentrations after 24-h exposure to 10 mM of the amino acid. However, arginine had no stimulating effect on PSP toxin levels if supplemented in combination with CAM at 10 microg l(-1). Addition of agmatine and proline to C. raciborskii T3 cultures in the presence of 10 microg l(-1) CAM increased C1+2 toxins levels, while having a negative or no effect on STX accumulation. In vitro, PSP toxin levels increased naturally in cyanobacterial extracts, with CAM and arginine having no influence on either C1+2 or STX synthesis. The evidence presented in this study suggests a possible difference between the metabolism of STX and the C1+2 toxins and indicated a high turnover rate of STX biosynthetic enzymes in C. raciborskii T3.

Exploring cyanobacterial genomes for natural product biosynthesis pathways.

Cyanobacteria produce a vast array of natural products, some of which are toxic to human health, while others possess potential pharmaceutical activities. Genome mining enables the identification and characterisation of natural product gene clusters; however, the current number of cyanobacterial genomes remains low compared to other phyla. There has been a recent effort to rectify this issue by increasing the number of sequenced cyanobacterial genomes. This has enabled the identification of biosynthetic gene clusters for structurally diverse metabolites, including non-ribosomal peptides, polyketides, ribosomal peptides, UV-absorbing compounds, alkaloids, terpenes and fatty acids. While some of the identified biosynthetic gene clusters correlate with known metabolites, genome mining also highlights the number and diversity of clusters for which the product is unknown (referred to as orphan gene clusters). A number of bioinformatic tools have recently been developed in order to predict the products of orphan gene clusters; however, in some cases the complexity of the cyanobacterial pathways makes the prediction problematic. This can be overcome by the use of mass spectrometry-guided natural product genome mining, or heterologous expression. Application of these techniques to cyanobacterial natural product gene clusters will be explored.

Absence of detectable levels of the cyanobacterial toxin (microcystin-LR) carry-over into milk.

The potential for the carry-over of the cyanobacterial toxin, microcystin-LR, from feed to milk was assessed using four Holstein-Friesian cows in a 4 week feeding trial. Two cows were used as control and the other two dosed daily at increasing weekly concentrations of microcystins from zero to a maximum dosage of 13 microg toxin kg x (-1) d x (-1) (or 7.4 mg toxin day(-1)). The absence of any deviation from the control in terms of physiological response and plasma indicators (total bilirubin, gamma-glutamyl transpeptidase and alkaline phosphatase) suggests that the microcystin-LR dosage did not have a detrimental effect on cattle liver function or milk yield during the course of the study. While the milk production did decrease over the period of the trial, no difference was observed between control and dosed cattle. Protein phosphatase inhibition assays were successfully used to determine the presence of microcystin-LR in prepared milk samples with an average recovery of 88% for samples spiked with 0.6 microg x l(-1) microcystin-LR. The level of microcystin-LR in all milk samples during the trial was less than 0.2 microg x l(-1). This suggests that after digestion, microcystin--LR is either not present in milk or sufficiently modified to render it non-toxic.

Molecular detection of genes responsible for cyanobacterial toxin production in the genera Microcystis, Nodularia, and Cylindrospermopsis.

Cyanobacteria are ubiquitous in the freshwater environment. Their success as a group in a wide range of aquatic habitats has been attributed to their unique physiological characteristics and their high adaptive ability over a wide range of environmental conditions. They are capable of reaching very high biomass levels, often dominating the other aquatic biota, and under some circumstances can accumulate near the water surface, producing scums. Such cyanobacterial "blooms" are of particular concern in reservoirs used to supply potable water. Dense aggregations of cyanobacterial cells may block water filters, and many species produce compounds that affect the taste and odor of water supplies. Of greatest concern, however, is the potential of many bloom-forming cyanobacteria to produce a wide range of toxic substances. These natural compounds, known as cyanotoxins, are chemically diverse and are usually either neuro- or hepatotoxic in pathology.

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.

Discovery of two cyanobacterial phenylalanine ammonia lyases: kinetic and structural characterization.

Phenylalanine ammonia lyase (PAL) catalyzes the deamination of phenylalanine to cinnamate and ammonia. While PALs are common in terrestrial plants where they catalyze the first committed step in the formation of phenylpropanoids, only a few prokaryotic PALs have been identified to date. Here we describe for the first time PALs from cyanobacteria, in particular, Anabaena variabilis ATCC 29413 and Nostoc punctiforme ATCC 29133, identified by screening the genome sequences of these organisms for members of the aromatic amino acid ammonia lyase family. Both PAL genes associate with secondary metabolite biosynthetic gene clusters as observed for other eubacterial PAL genes. In comparison to eukaryotic homologues, the cyanobacterial PALs are 20% smaller in size but share similar substrate selectivity and kinetic activity toward L-phenylalanine over L-tyrosine. Structure elucidation by protein X-ray crystallography confirmed that the two cyanobacterial PALs are similar in tertiary and quatenary structure to plant and yeast PALs as well as the mechanistically related histidine ammonia lyases.

Molecular basis for chloronium-mediated meroterpene cyclization: cloning, sequencing, and heterologous expression of the napyradiomycin biosynthetic gene cluster.

Structural inspection of the bacterial meroterpenoid antibiotics belonging to the napyradiomycin family of chlorinated dihydroquinones suggests that the biosynthetic cyclization of their terpenoid subunits is initiated via a chloronium ion. The vanadium-dependent haloperoxidases that catalyze such reactions are distributed in fungi and marine algae and have yet to be characterized from bacteria. The cloning and sequence analysis of the 43-kb napyradiomycin biosynthetic cluster (nap) from Streptomyces aculeolatus NRRL 18422 and from the undescribed marine sediment-derived Streptomyces sp. CNQ-525 revealed 33 open reading frames, three of which putatively encode vanadium-dependent chloroperoxidases. Heterologous expression of the CNQ-525-based nap biosynthetic cluster in Streptomyces albus produced at least seven napyradiomycins, including the new analog 2-deschloro-2-hydroxy-A80915C. These data not only revealed the molecular basis behind the biosynthesis of these novel meroterpenoid natural products but also resulted in the first in vivo verification of vanadium-dependent haloperoxidases.

Structural characterization of in vitro and in vivo intermediates on the loading module of microcystin synthetase.

The microcystin family of toxins is the most common cause of hepatotoxicity associated with water blooms of cyanobacterial genera. The biosynthetic assembly line producing the toxic cyclic peptide, microcystin, contains an adenylation-peptidyl carrier protein didomain (A-PCP) at the N-terminus of the initiator module McyG (295 kDa) that has been postulated to activate and load the starter unit phenylacetate for formation of the unusual aromatic beta-amino acid residue, Adda, before subsequent extension. Characterization of the McyG A-PCP didomain (78 kDa) using ATP-PP i exchange assays and mass spectrometry revealed that assorted phenylpropanoids are preferentially activated and loaded onto the PCP carrier domain rather than phenylacetate itself. For the first time, thioesters formed in vivo were detected directly using large molecule mass spectrometry. Additionally substrates were cleaved using a type II thioesterase for structural elucidation by small molecule mass spectrometry. Unprecedented features of the McyG A-PCP didomain include the in vivo acylation of the holo PCP with exogenous and endogenous substrates, along with the ability of the apo protein to retain the acyl-AMP intermediate during affinity purification. These results imply that phenylpropanoids are preferentially loaded onto the McyG PCP; however one carbon must be excised following extension of the starter unit with malonyl-CoA in order to generate the expected polyketide chain which leads us to ponder the novel biochemistry by which this occurs.

Characterization of the nodularin synthetase gene cluster and proposed theory of the evolution of cyanobacterial hepatotoxins.

Nodularia spumigena is a bloom-forming cyanobacterium which produces the hepatotoxin nodularin. The complete gene cluster encoding the enzymatic machinery required for the biosynthesis of nodularin in N. spumigena strain NSOR10 was sequenced and characterized. The 48-kb gene cluster consists of nine open reading frames (ORFs), ndaA to ndaI, which are transcribed from a bidirectional regulatory promoter region and encode nonribosomal peptide synthetase modules, polyketide synthase modules, and tailoring enzymes. The ORFs flanking the nda gene cluster in the genome of N. spumigena strain NSOR10 were identified, and one of them was found to encode a protein with homology to previously characterized transposases. Putative transposases are also associated with the structurally related microcystin synthetase (mcy) gene clusters derived from three cyanobacterial strains, indicating a possible mechanism for the distribution of these biosynthetic gene clusters between various cyanobacterial genera. We propose an alternative hypothesis for hepatotoxin evolution in cyanobacteria based on the results of comparative and phylogenetic analyses of the nda and mcy gene clusters. These analyses suggested that nodularin synthetase evolved from a microcystin synthetase progenitor. The identification of the nodularin biosynthetic gene cluster and evolution of hepatotoxicity in cyanobacteria reported in this study may be valuable for future studies on toxic cyanobacterial bloom formation. In addition, an appreciation of the natural evolution of nonribosomal biosynthetic pathways will be vital for future combinatorial engineering and rational design of novel metabolites and pharmaceuticals.