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.     

Naturally mosaic operons for secondary metabolite biosynthesis: variability and putative horizontal transfer of discrete catalytic domains of the epothilone polyketide synthase locus

A putative instance of horizontal gene transfer (HGT) involving adjacent, discrete beta-ketoacyl synthase (KS), acyl carrier protein (ACP) and nonribosomal peptide synthase (NRPS) domains of the epothilone Type I polyketide biosynthetic gene cluster from the myxobacterium Sorangium cellulosom was identified using molecular phylogenetics and sequence analyses. The specific KS domain of the module EPO B fails to cluster phylogenetically with other epothilone KS sequences present at this locus, in contrast to what is typically observed in many other Type I polyketide synthase (PKS) biosynthetic loci. Furthermore, the GC content of the epoB KS, epoA ACP and NRPS domains differs significantly from the base composition of other epothilone domain sequences. In addition, the putatively transferred epothilone loci are located near previously identified transposon-like sequences. Lastly, comparison with other KS loci revealed another possible case of horizontal transfer of secondary metabolite genes in the genus Pseudomonas. This study emphasizes the use of several lines of concordant evidence (phylogenetics, base composition, transposon sequences) to infer the evolutionary history of particular gene and enzyme sequences, and the results support the idea that genes coding for adaptive traits, e.g. defensive natural products, may be prone to transposition between divergent prokaryotic taxa and genomes.     

Evolution of the aflatoxin gene cluster

WhyAspergillus species produce aflatoxin remains an unsolved question. In this report we suggest that evolution of the aflatoxin biosynthesis gene cluster has been a multistep process. More than 300 million years ago a primordial cluster of genes allowed production of anthraquinones that may have served as insect attractants to facilitate spore dispersal. Later adaptive evolutionary steps introduced genes into the cluster that encoded enzymes associated with fungal virulence. These genes may have allowed the otherwise saprophytic fungi to be better able to colonize living plants. Later, genes for production of aflatoxins B1 and G1 were added to the basal cluster. Loss of the ability to produce aflatoxin G1 occurred with the divergence ofA. flavus, a species that, perhaps, was more successful than its ancestors at colonizing plants. This logical progression in evolutionary development of the aflatoxin biosynthetic cluster fits the phylogenetic data as well as known chemical reactivity of the initially formed anthraquinone polyketide metabolites.     

Identification and analysis of the chivosazol biosynthetic gene cluster from the myxobacterial model strain Sorangium cellulosum So ce56

Myxobacteria belonging to the genus Sorangium are known to produce a variety of biologically active secondary metabolites. Chivosazol is a macrocyclic antibiotic active against yeast, filamentous fungi and especially against mammalian cells. The compound specifically destroys the actin skeleton of eucaryotic cells and does not show activity against bacteria. Chivosazol contains an oxazole ring and a glycosidically bound 6-deoxyglucose (except for chivosazol F). In this paper we describe the biosynthetic gene cluster that directs chivosazol biosynthesis in the model strain Sorangium cellulosum So ce56. This biosynthetic gene cluster spans 92 kbp on the chromosome and contains four polyketide synthase genes and one hybrid polyketide synthase/nonribosomal peptide synthetase gene. An additional gene encoding a protein with similarity to different methyltransferases and presumably involved in post-polyketide modification was identified downstream of the core biosynthetic gene cluster. The chivosazol biosynthetic gene locus belongs to the recently identified and rapidly growing class of trans-acyltransferase polyketide synthases, which do not contain acyltransferase domains integrated into the multimodular megasynthetases.     

Natural biocombinatorics in the polyketide synthase genes of the actinobacterium Streptomyces avermitilis

Modular polyketide synthases (PKSs) of bacteria provide an enormous reservoir of natural chemical diversity. Studying natural biocombinatorics may aid in the development of concepts for experimental design of genes for the biosynthesis of new bioactive compounds. Here we address the question of how the modularity of biosynthetic enzymes and the prevalence of multiple gene clusters in Streptomyces drive the evolution of metabolic diversity. The phylogeny of ketosynthase (KS) domains of Streptomyces PKSs revealed that the majority of modules involved in the biosynthesis of a single compound evolved by duplication of a single ancestor module. Using Streptomyces avermitilis as a model organism, we have reconstructed the evolutionary relationships of different domain types. This analysis suggests that 65% of the modules were altered by recombinational replacements that occurred within and between biosynthetic gene clusters. The natural reprogramming of the biosynthetic pathways was unambiguously confined to domains that account for the structural diversity of the polyketide products and never observed for the KS domains. We provide examples for natural acyltransferase (AT), ketoreductase (KR), and dehydratase (DH)–KR domain replacements. Potential sites of homologous recombination could be identified in interdomain regions and within domains. Our results indicate that homologous recombination facilitated by the modularity of PKS architecture is the most important mechanism underlying polyketide diversity in bacteria.     

Naturally mosaic operons for secondary metabolite biosynthesis: variability and putative horizontal transfer of discrete catalytic domains of the epothilone polyketide synthase locus

A putative instance of horizontal gene transfer (HGT) involving adjacent, discrete -ketoacyl synthase (KS), acyl carrier protein (ACP) and nonribosomal peptide synthase (NRPS) domains of the epothilone Type I polyketide biosynthetic gene cluster from the myxobacterium Sorangium cellulosom was identified using molecular phylogenetics and sequence analyses. The specific KS domain of the module EPO B fails to cluster phylogenetically with other epothilone KS sequences present at this locus, in contrast to what is typically observed in many other Type I polyketide synthase (PKS) biosynthetic loci. Furthermore, the GC content of the epoB KS, epoA ACP and NRPS domains differs significantly from the base composition of other epothilone domain sequences. In addition, the putatively transferred epothilone loci are located near previously identified transposon-like sequences. Lastly, comparison with other KS loci revealed another possible case of horizontal transfer of secondary metabolite genes in the genus Pseudomonas. This study emphasizes the use of several lines of concordant evidence (phylogenetics, base composition, transposon sequences) to infer the evolutionary history of particular gene and enzyme sequences, and the results support the idea that genes coding for adaptive traits, e.g. defensive natural products, may be prone to transposition between divergent prokaryotic taxa and genomes.Communicated by W. Arber     

Microcystin biosynthesis in planktothrix: genes,evolution, and manipulation

Microcystins represent an extraordinarily large family of cyclic heptapeptide toxins that are nonribosomally synthesized by various cyanobacteria. Microcystins specifically inhibit the eukaryotic protein phosphatases 1 and 2A. Their outstanding variability makes them particularly useful for studies on the evolution of structure-function relationships in peptide synthetases and their genes. Analyses of microcystin synthetase genes provide valuable clues for the potential and limits of combinatorial biosynthesis. We have sequenced and analyzed 55.6 kb of the potential microcystin synthetase gene (mcy) cluster from the filamentous cyanobacterium Planktothrix agardhii CYA 126. The cluster contains genes for peptide synthetases (mcyABC), polyketide synthases (PKSs; mcyD), chimeric enzymes composed of peptide synthetase and PKS modules (mcyEG), a putative thioesterase (mcyT), a putative ABC transporter (mcyH), and a putative peptide-modifying enzyme (mcyJ). The gene content and arrangement and the sequence of specific domains in the gene products differ from those of the mcy cluster in Microcystis, a unicellular cyanobacterium. The data suggest an evolution of mcy clusters from, rather than to, genes for nodularin (a related pentapeptide) biosynthesis. Our data do not support the idea of horizontal gene transfer of complete mcy gene clusters between the genera. We have established a protocol for stable genetic transformation of Planktothrix, a genus that is characterized by multicellular filaments exhibiting continuous motility. Targeted mutation of mcyJ revealed its function as a gene coding for a O-methyltransferase. The mutant cells produce a novel microcystin variant exhibiting reduced inhibitory activity toward protein phosphatases.     

Ancient gene duplications and the root(s) of the tree of life

Summary. Tracing organismal histories on the timescale of the tree of life remains one of the challenging tasks in evolutionary biology. The hotly debated questions include the evolutionary relationship between the three domains of life (e.g., which of the three domains are sister domains, are the domains para-, poly-, or monophyletic) and the location of the root within the universal tree of life. For the latter, many different points of view have been considered but so far no consensus has been reached. The only widely accepted rationale to root the universal tree of life is to use anciently duplicated paralogous genes that are present in all three domains of life. To date only few anciently duplicated gene families useful for phylogenetic reconstruction have been identified. Here we present results from a systematic search for ancient gene duplications using twelve representative, completely sequenced, archaeal and bacterial genomes. Phylogenetic analyses of identified cases show that the majority of datasets support a root between Archaea and Bacteria; however, some datasets support alternative hypotheses, and all of them suffer from a lack of strong phylogenetic signal. The results are discussed with respect to the impact of horizontal gene transfer on the ability to reconstruct organismal evolution. The exchange of genetic information between divergent organisms gives rise to mosaic genomes, where different genes in a genome have different histories. Simulations show that even low rates of horizontal gene transfer dramatically complicate the reconstruction of organismal evolution, and that the different most recent common molecular ancestors likely existed at different times and in different lineages. Correspondence and reprints: Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT 06269-3125, U.S.A. Present address: Genome Atlantic, Department of Biochemistry and Molecular Biology, Dalhousie University, Halifax, Nova Scotia, Canada.     

From genetic diversity to metabolic unity: studies on the biosynthesis of aurafurones and aurafuron-like structures in myxobacteria and streptomycetes

The myxobacterial polyketide secondary metabolites aurafuron A and B were identified by genome mining in the myxobacterial strain Stigmatella aurantiaca DW4/3-1. The compounds contain an unusual furanone moiety and resemble metabolites isolated from soil-dwelling and marine actinobacteria, a fungus and mollusks. We describe here the cloning and functional analysis of the aurafuron biosynthetic gene cluster, including site-directed mutagenesis and feeding studies using labeled precursors. The polyketide core of the aurafurones is assembled by a modular polyketide synthase (PKS). As with many such systems described from myxobacteria, the aurafuron PKS exhibits a number of unusual features, including the apparent iterative use of a module, redundant modules and domains, a trans acting dehydratase and the absence of a terminal thioesterase domain. Four oxidoreductases are encoded within the gene locus, some of which likely participate in formation of the furanone moiety via a Baeyer-Villiger type oxidation. Indeed, inactivation of a gene encoding a cytochrome P₄₅₀ monooxygenase completely abolished production of both compounds. We also compare the complete gene locus to biosynthetic gene clusters from two Streptomyces sp., which produce close structural analogues of the aurafurones. A portion of the post-PKS biosynthetic machinery is strikingly similar in all three cases, in contrast to the PKS genes, which are highly divergent. Phylogenetic analysis of the ketosynthase domains further indicates that the PKSs have developed independently (polyphyletically) during evolution. These findings point to a currently unknown but important biological function of aurafuron-like compounds for the producing organisms.     

Phylogenetic Analysis of Antibiotic Glycosyltransferases

Catalyzed by a family of enzymes called glycosyltransferases, glycosylation reactions are essential for the bioactivities of secondary metabolites such as antibiotics. Due to the special characters of antibiotic glycosyltransferases (AGts), antibiotics can function by attaching some unusual deoxy-sugars to their aglycons. Comprehensive similarity searches on the amino acid sequences of AGts have been performed. We reconstructed the molecular phylogeny of AGts with neighbor-joining, maximum-likelihood, and Bayesian methods of phylogenetic inference. The phylogenetic trees show a distinct separation of polyene macrolide (PEM) AGts and other polyketide AGts. The former are more like eukaryotic glycosyltransferases and were deduced to be the results of horizontal gene transfer from eukaryotes. Protein tertiary structural comparison also indicated that some glycopeptide AGts (Gtf-proteins) have a close evolutionary relationship with MurGs, essential glycosyltransferases involved in maturation of bacterial cell walls. The evolutionary relationship of glycopeptide antibiotic biosynthetic gene clusters was speculated according to the phylogenetic analysis of Gtf-proteins. Considering the fact that polyketide AGts and Gtf-proteins are all GT Family 1 members and their aglycon acceptor biosynthetic patterns are very similar, we deduced that AGts and the synthases of their aglycon acceptors have some evolutionary relevance. Finally, the evolutionary origins of AGts that do not fall into GT Family 1 are discussed, suggesting that their ancestral proteins appear to be derived from various proteins responsible for primary metabolism. [Reviewing Editor: Dr. Niles Lehman]     

Ancient Horizontal Gene Transfer from Bacteria Enhances Biosynthetic Capabilities of Fungi

Background

Polyketides are natural products with a wide range of biological functions and pharmaceutical applications. Discovery and utilization of polyketides can be facilitated by understanding the evolutionary processes that gave rise to the biosynthetic machinery and the natural product potential of extant organisms. Gene duplication and subfunctionalization, as well as horizontal gene transfer are proposed mechanisms in the evolution of biosynthetic gene clusters. To explain the amount of homology in some polyketide synthases in unrelated organisms such as bacteria and fungi, interkingdom horizontal gene transfer has been evoked as the most likely evolutionary scenario. However, the origin of the genes and the direction of the transfer remained elusive.

Methodology/Principal Findings

We used comparative phylogenetics to infer the ancestor of a group of polyketide synthase genes involved in antibiotic and mycotoxin production. We aligned keto synthase domain sequences of all available fungal 6-methylsalicylic acid (6-MSA)-type PKSs and their closest bacterial relatives. To assess the role of symbiotic fungi in the evolution of this gene we generated 24 6-MSA synthase sequence tags from lichen-forming fungi. Our results support an ancient horizontal gene transfer event from an actinobacterial source into ascomycete fungi, followed by gene duplication.

Conclusions/Significance

Given that actinobacteria are unrivaled producers of biologically active compounds, such as antibiotics, it appears particularly promising to study biosynthetic genes of actinobacterial origin in fungi. The large number of 6-MSA-type PKS sequences found in lichen-forming fungi leads us hypothesize that the evolution of typical lichen compounds, such as orsellinic acid derivatives, was facilitated by the gain of this bacterial polyketide synthase.     

Wide distribution of O157-antigen biosynthesis gene clusters in Escherichia coli

Most Escherichia coli O157-serogroup strains are classified as enterohemorrhagic E. coli (EHEC), which is known as an important food-borne pathogen for humans. They usually produce Shiga toxin (Stx) 1 and/or Stx2, and express H7-flagella antigen (or nonmotile). However, O157 strains that do not produce Stxs and express H antigens different from H7 are sometimes isolated from clinical and other sources. Multilocus sequence analysis revealed that these 21 O157:non-H7 strains tested in this study belong to multiple evolutionary lineages different from that of EHEC O157:H7 strains, suggesting a wide distribution of the gene set encoding the O157-antigen biosynthesis in multiple lineages. To gain insight into the gene organization and the sequence similarity of the O157-antigen biosynthesis gene clusters, we conducted genomic comparisons of the chromosomal regions (about 59 kb in each strain) covering the O-antigen gene cluster and its flanking regions between six O157:H7/non-H7 strains. Gene organization of the O157-antigen gene cluster was identical among O157:H7/non-H7 strains, but was divided into two distinct types at the nucleotide sequence level. Interestingly, distribution of the two types did not clearly follow the evolutionary lineages of the strains, suggesting that horizontal gene transfer of both types of O157-antigen gene clusters has occurred independently among E. coli strains. Additionally, detailed sequence comparison revealed that some positions of the repetitive extragenic palindromic (REP) sequences in the regions flanking the O-antigen gene clusters were coincident with possible recombination points. From these results, we conclude that the horizontal transfer of the O157-antigen gene clusters induced the emergence of multiple O157 lineages within E. coli and speculate that REP sequences may involve one of the driving forces for exchange and evolution of O-antigen loci.     

Non-ribosomal peptide synthetases: Identifying the cryptic gene clusters and decoding the natural product

Non-ribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) present in bacteria and fungi are the major multi-modular enzyme complexes which synthesize secondary metabolites like the pharmacologically important antibiotics and siderophores. Each of the multiple modules of an NRPS activates a different amino or aryl acid, followed by their condensation to synthesize a linear or cyclic natural product. The studies on NRPS domains, the knowledge of their gene cluster architecture and tailoring enzymes have helped in the in silico genetic screening of the ever-expanding sequenced microbial genomic data for the identification of novel NRPS/PKS clusters and thus deciphering novel non-ribosomal peptides (NRPs). Adenylation domain is an integral part of the NRPSs and is the substrate selecting unit for the final assembled NRP. In some cases, it also requires a small protein, the MbtH homolog, for its optimum activity. The presence of putative adenylation domain and MbtH homologs in a sequenced genome can help identify the novel secondary metabolite producers. The role of the adenylation domain in the NRPS gene clusters and its characterization as a tool for the discovery of novel cryptic NRPS gene clusters are discussed.     

Comparison of cyanopeptolin genes in Planktothrix, Microcystis, and Anabaena strains: evidence for independent evolution within each genus

The major cyclic peptide cyanopeptolin 1138, produced by Planktothrix strain NIVA CYA 116, was characterized and shown to be structurally very close to the earlier-characterized oscillapeptin E. A cyanopeptolin gene cluster likely to encode the corresponding peptide synthetase was sequenced from the same strain. The 30-kb oci gene cluster contains two novel domains previously not detected in nonribosomal peptide synthetase gene clusters (a putative glyceric acid-activating domain and a sulfotransferase domain), in addition to seven nonribosomal peptide synthetase modules. Unlike in two previously described cyanopeptolin gene clusters from Anabaena and Microcystis, a halogenase gene is not present. The three cyanopeptolin gene clusters show similar gene and domain arrangements, while the binding pocket signatures deduced from the adenylation domain sequences and the additional tailoring domains vary. This suggests loss and gain of tailoring domains within each genus, after the diversification of the three clades, as major events leading to the present diversity. The ABC transporter genes associated with the cyanopeptolin gene clusters form a monophyletic clade and accordingly are likely to have evolved as part of the functional unit. Phylogenetic analyses of adenylation and condensation domains, including domains from cyanopeptolins and microcystins, show a closer similarity between the Planktothrix and Microcystis cyanopeptolin domains than between these and the Anabaena domain. No clear evidence of recombination between cyanopeptolins and microcystins could be detected. There were no strong indications of horizontal gene transfer of cyanopeptolin gene sequences across the three genera, supporting independent evolution within each genus.     

Cloning and sequence analysis of the putative rifamycin polyketide synthase gene cluster from Amycolatopsis mediterranei

The 54-kbp Type I polyketide synthase gene cluster, most probably involved in rifamycin biosynthesis by Amycolatopsis mediterranei, was cloned in E. coli and completely sequenced. The DNA encodes five closely packed, very large open reading frames reading in one direction. As expected from the chemical structure of rifamycins, ten polyketide synthase modules and a CoA ligase domain were identified in the five open reading frames which contain one to three polyketide synthase modules each. The order of the functional domains on the DNA probably reflects the order in which they are used because each of the modules contains the predicted acetate or propionate transferase, dehydratase, and β-ketoacyl-ACP reductase functions, required for the respective step in rifamycin biosynthesis.     

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.     

Diversity and evolution of polyketide biosynthesis gene clusters in the Ceratocystidaceae

Polyketides are secondary metabolites with diverse biological activities. Polyketide synthases (PKS) are often encoded from genes clustered in the same genomic region. Functional analyses and genomic studies show that most fungi are capable of producing a repertoire of polyketides. We considered the potential of Ceratocystidaceae for producing polyketides using a comparative genomics approach. Our aims were to identify the putative polyketide biosynthesis gene clusters, to characterize them and predict the types of polyketide compounds they might produce. We used sequences from nineteen species in the genera, Ceratocystis, Endoconidiophora, Davidsoniella, Huntiella, Thielaviopsis and Bretziella, to identify and characterize PKS gene clusters, by employing a range of bioinformatics and phylogenetic tools. We showed that the genomes contained putative clusters containing a non-reducing type I PKS and a type III PKS. Phylogenetic analyses suggested that these genes were already present in the ancestor of the Ceratocystidaceae. By contrast, the various reducing type I PKS-containing clusters identified in these genomes appeared to have distinct evolutionary origins. Although one of the identified clusters potentially allows for the production of melanin, their functional characterization will undoubtedly reveal many novel and important compounds implicated in the biology of the Ceratocystidaceae.     

Evolution of Chemical Diversity in Echinocandin Lipopeptide Antifungal Metabolites

The echinocandins are a class of antifungal drugs that includes caspofungin, micafungin, and anidulafungin. Gene clusters encoding most of the structural complexity of the echinocandins provided a framework for hypotheses about the evolutionary history and chemical logic of echinocandin biosynthesis. Gene orthologs among echinocandin-producing fungi were identified. Pathway genes, including the nonribosomal peptide synthetases (NRPSs), were analyzed phylogenetically to address the hypothesis that these pathways represent descent from a common ancestor. The clusters share cooperative gene contents and linkages among the different strains. Individual pathway genes analyzed in the context of similar genes formed unique echinocandin-exclusive phylogenetic lineages. The echinocandin NRPSs, along with the NRPS from the inp gene cluster in Aspergillus nidulans and its orthologs, comprise a novel lineage among fungal NRPSs. NRPS adenylation domains from different species exhibited a one-to-one correspondence between modules and amino acid specificity that is consistent with models of tandem duplication and subfunctionalization. Pathway gene trees and Ascomycota phylogenies are congruent and consistent with the hypothesis that the echinocandin gene clusters have a common origin. The disjunct Eurotiomycete-Leotiomycete distribution appears to be consistent with a scenario of vertical descent accompanied by incomplete lineage sorting and loss of the clusters from most lineages of the Ascomycota. We present evidence for a single evolutionary origin of the echinocandin family of gene clusters and a progression of structural diversification in two fungal classes that diverged approximately 290 to 390 million years ago. Lineage-specific gene cluster evolution driven by selection of new chemotypes contributed to diversification of the molecular functionalities.     

Birth,death and horizontal transfer of the fumonisin biosynthetic gene cluster during the evolutionary diversification of Fusarium

Fumonisins are a family of carcinogenic secondary metabolites produced by members of the Fusarium fujikuroi species complex (FFSC) and rare strains of Fusarium oxysporum. In Fusarium, fumonisin biosynthetic genes (FUM) are clustered, and the cluster is uniform in gene organization. Here, sequence analyses indicated that the cluster exists in five different genomic contexts, defining five cluster types. In FUM gene genealogies, evolutionary relationships between fusaria with different cluster types were largely incongruent with species relationships inferred from primary‐metabolism (PM) gene genealogies, and FUM cluster types are not trans‐specific. In addition, synonymous site divergence analyses indicated that three FUM cluster types predate diversification of FFSC. The data are not consistent with balancing selection or interspecific hybridization, but they are consistent with two competing hypotheses: (i) multiple horizontal transfers of the cluster from unknown donors to FFSC recipients and (ii) cluster duplication and loss (birth and death). Furthermore, low levels of FUM gene divergence in F. bulbicola, an FFSC species, and F. oxysporum provide evidence for horizontal transfer of the cluster from the former, or a closely related species, to the latter. Thus, uniform gene organization within the FUM cluster belies a complex evolutionary history that has not always paralleled the evolution of Fusarium.