Identification of a Saxitoxin Biosynthesis Gene with a History of Frequent Horizontal Gene Transfers

The paralytic shellfish poisoning (PSP) toxins, saxitoxin, and its derivatives, are produced by a complex and unique biosynthetic pathway. It involves reactions that are rare in other metabolic pathways, however, distantly related organisms, such as dinoflagellates and cyanobacteria, produce these toxins by an identical pathway. Speculative explanations for the unusual phylogenetic distribution of this metabolic pathway have been proposed, including a polyphyletic origin, the involvement of symbiotic bacteria, and horizontal gene transfer. This study describes for the first time the identity of one gene, sxt1, that is involved in the biosynthesis of saxitoxin in cyanobacteria. It encoded an O-carbamoyltransferase (OCTASE) that was proposed to carbamoylate the hydroxymethyl side chain of saxitoxin precursor. Orthologues of sxt1 were exclusively present in PSP-toxic strains of cyanobacteria and had a high sequence similarity to each other. L. wollei had a naturally mutated sxt1 gene that encoded an inactive enzyme, and was incapable of producing carbamoylated PSP-toxin analogues, supporting the proposed function of Sxt1. Phylogenetic analysis revealed that OCATSE genes were present exclusively in prokaryotic organisms and were characterized by a high rate of horizontal gene transfer. OCTASE has most likely evolved from an ancestral O-sialoglycoprotein endopeptidase from proteobacteria, whereas the most likely phylogenetic origin of sxt1 was an ancestral alpha-proteobacterium. The phylogeny of sxt1 suggested that the entire set of genes required for saxitoxin biosynthesis may spread by horizontal gene transfer.     

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

Background  

Cyanobacteria produce a wealth of secondary metabolites, including the group of small cyclic heptapeptide hepatotoxins that constitutes the microcystin family. The enzyme complex that directs the biosynthesis of microcystin is encoded in a single large gene cluster (mcy). mcy genes have a widespread distribution among cyanobacteria and are likely to have an ancient origin. The notable diversity within some of the Mcy modules is generated through various recombination events including horizontal gene transfer.     

Evolutionary mechanisms underlying the functional divergence of duplicate genes involved in vertebrates' circadian rhythm pathway

Circadian rhythms, that are governed physiologically and behaviorally by endogenous clock, have been described in many species. Living organisms use this endogenous circadian clock to anticipate environmental transitions, perform activities at biologically advantageous times during the day, and undergo characteristic seasonal responses. Gene duplication is one of the most important mechanisms in the evolution of gene diversity. After duplication, one or both of duplicates can accumulate amino acid changes, thereby promoting functional divergence through the action of natural selection. The circadian system, like many other multigene families, has undergone this genetic revolution, and so circadian genes that are found in single copies in insects are duplicated in vertebrates. We analyzed six groups of genes involved in vertebrates' circadian rhythm pathway to find signatures of molecular evolutionary processes such as gene duplication, natural selection, recombination, and functional divergence. The obtained results, then, were used to determine what evolutionary forces have influenced the fates of duplicated genes of each group. We showed in this research that recombination has not been widespread during the evolution of circadian genes and that purifying selection has been the prominent natural pressure operating on circadian genes. We also showed that the evolution of circadian genes has been depended on gene duplication and functional divergence. Finally, we put forward models best describing the evolutionary fates of circadian duplicates.     

Potent toxins in Arctic environments – Presence of saxitoxins and an unusual microcystin variant in Arctic freshwater ecosystems

Cyanobacteria are the predominant phototrophs in freshwater ecosystems of the polar regions where they commonly form extensive benthic mats. Despite their major biological role in these ecosystems, little attention has been paid to their physiology and biochemistry. An important feature of cyanobacteria from the temperate and tropical regions is the production of a large variety of toxic secondary metabolites. In Antarctica, and more recently in the Arctic, the cyanobacterial toxins microcystin and nodularin (Antarctic only) have been detected in freshwater microbial mats. To date other cyanobacterial toxins have not been reported from these locations. Five Arctic cyanobacterial communities were screened for saxitoxin, another common cyanobacterial toxin, and microcystins using immunological, spectroscopic and molecular methods. Saxitoxin was detected for the first time in cyanobacteria from the Arctic. In addition, an unusual microcystin variant was identified using liquid chromatography–mass spectrometry. Gene expression analyses confirmed the analytical findings, whereby parts of the sxt and mcy operon involved in saxitoxin and microcystin synthesis, were detected and sequenced in one and five of the Arctic cyanobacterial samples, respectively. The detection of these compounds in the cryosphere improves the understanding of the biogeography and distribution of toxic cyanobacteria globally. The sequences of sxt and mcy genes provided from this habitat for the first time may help to clarify the evolutionary origin of toxin production in cyanobacteria.     

The role of recombination,niche‐specific gene pools and flexible genomes in the ecological speciation of bacteria

Bacteria diversify into genetic clusters analogous to those observed in sexual eukaryotes, but the definition of bacterial species is an ongoing problem. Recent work has focused on adaptation to distinct ecological niches as the main driver of clustering, but there remains debate about the role of recombination in that process. One view is that homologous recombination occurs too rarely for gene flow to constrain divergent selection. Another view is that homologous recombination is frequent enough in many bacterial populations that barriers to gene flow are needed to permit divergence. Niche‐specific gene pools have been proposed as a general mechanism to limit gene flow. We use theoretical models to evaluate additional hypotheses that evolving genetic architecture, specifically the effect sizes of genes and gene gain and loss, can limit gene flow between diverging populations. Our model predicts that (a) in the presence of gene flow and recombination, ecological divergence is concentrated in few loci of large effect and (b) high rates of gene flow plus recombination promote gene loss and favor the evolution of niche‐specific genes. The results show that changing genetic architecture and gene loss can facilitate ecological divergence, even without niche‐specific gene pools. We discuss these results in the context of recent studies of sympatric divergence in microbes.     

Extreme Recombination Frequencies Shape Genome Variation and Evolution in the Honeybee,Apis mellifera

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.     

Detection of cyanobacterial sxt genes and paralytic shellfish toxins in freshwater lakes and brackish waters on Åland Islands,Finland

Harmful cyanobacteria are a globally growing concern. They produce a large variety of toxic compounds, including saxitoxin and its many structural variants, a group of potent neurotoxins collectively called paralytic shellfish toxins or PST. Nucleic acid based detection methods, such as qPCR, have been proposed as potential screening and monitoring tools for toxic cyanobacteria, but it is not clear how well the presence and quantity of saxitoxin biosynthesis (sxt) genes can be used to predict the production of PST in the environment. In this study, the prevalence of three sxt genes and their co-occurrence with paralytic shellfish toxins in the environment was investigated. The sxtA, sxtG and sxtB genes were present on average in 31% of the samples collected from lakes and brackish coastal waters on Åland Islands, Finland, during the three-year monitoring period. PST detection frequency varied from 13% to 59% from year to year, and concentrations were generally low. On average higher sxtB copy numbers were associated with PST detection, and although a positive correlation between gene copy numbers and toxin concentrations was observed (Spearman rank correlation, ρ = 0.53, P = 0.012), sxt gene presence or quantity didn’t reliably predict PST production. Sequencing of sxtA fragments and identification of main cyanobacteria indicated that the likely candidate responsible for PST production in the samples belonged to the genus Anabaena.     

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.     

Iteratively improving natamycin production in Streptomyces gilvosporeus by a large operon-reporter based strategy

Many high-value secondary metabolites are assembled by very large multifunctional polyketide synthases or non-ribosomal peptide synthetases encoded by giant genes, for instance, natamycin production in an industrial strain of Streptomyces gilvosporeus. In this study, a large operon reporter-based selection system has been developed using the selectable marker gene neo to report the expression both of the large polyketide synthase genes and of the entire gene cluster, thereby facilitating the selection of natamycin-overproducing mutants by iterative random mutagenesis breeding. In three successive rounds of mutagenesis and selection, the natamycin titer was increased by 110%, 230%, and 340%, respectively, and the expression of the whole biosynthetic gene cluster was correspondingly increased. An additional copy of the natamycin gene cluster was found in one overproducer. These findings support the large operon reporter-based selection system as a useful tool for the improvement of industrial strains utilized in the production of polyketides and non-ribosomal peptides.     

Genetic tests of ancient asexuality in Root Knot Nematodes reveal recent hybrid origins

Background  

The existence of "ancient asexuals", taxa that have persisted for long periods of evolutionary history without sexual recombination, is both controversial and important for our understanding of the evolution and maintenance of sexual reproduction. A lack of sex has consequences not only for the ecology of the asexual organism but also for its genome. Several genetic signatures are predicted from long-term asexual (apomictic) reproduction including (i) large "allelic" sequence divergence (ii) lack of phylogenetic clustering of "alleles" within morphological species and (iii) decay and loss of genes specific to meiosis and sexual reproduction. These genetic signatures can be hard to assess since it is difficult to demonstrate the allelic nature of very divergent sequences, divergence levels may be complicated by processes such as inter-specific hybridization, and genes may have secondary roles unrelated to sexual reproduction. Apomictic species of Meloidogyne root knot nematodes have been suggested previously to be ancient asexuals. Their relatives reproduce by meiotic parthenogenesis or facultative sexuality, which in combination with the abundance of nematode genomic sequence data, makes them a powerful system in which to study the consequences of reproductive mode on genomic divergence.     

Evidence for Recombination in the Microcystin Synthetase (<Emphasis Type="Italic">mcy</Emphasis>) Genes ofToxic Cyanobacteria <Emphasis Type="Italic">Microcystis</Emphasis><Emphasis Type="Bold">spp</Emphasis>

Recombination has been suggested to be an important factor for the genetic variation of bacterial genes, but few studies have dealt with intragenic recombination between the same or closely related species of cyanobacteria. Here we provide strong evidence for recombination in the microcystin synthetase (mcy) gene cluster of the toxic cyanobacteria Microcystis spp. This gene cluster contains 10 genes (mcyA to J) that encode a mixed polyketide synthase (PKS)/nonribosomal peptide synthetase (NRPS) complex. mcy gene sequences were determined for four selected regions (within mcyA, D, G, and J) within the mcy gene cluster from 1 Canadian and 10 Asian toxic Microcystis and compared with previously published mcy sequences. Split decomposition analysis indicated a reticulate phylogeny of mcyA, and several potential recombination tracts of mcyA were identified by the RDP analysis and a runs test implemented in GENECONV. In contrast, no recombination was detected in the mcyD, G, and J sequences. However, discrepancies among the four mcy gene genealogies were evident from the results of independent split decomposition analyses, which were further supported by incongruence length difference (ILD) tests. Taken together, these findings suggest that both intragenic and intergenic recombination within the mcy gene cluster contributes to the genetic diversity of the mcy genes of Microcystis spp.This article contains online supplementary material.     

Single cell genome amplification accelerates identification of the apratoxin biosynthetic pathway from a complex microbial assemblage

Filamentous marine cyanobacteria are extraordinarily rich sources of structurally novel, biomedically relevant natural products. To understand their biosynthetic origins as well as produce increased supplies and analog molecules, access to the clustered biosynthetic genes that encode for the assembly enzymes is necessary. Complicating these efforts is the universal presence of heterotrophic bacteria in the cell wall and sheath material of cyanobacteria obtained from the environment and those grown in uni-cyanobacterial culture. Moreover, the high similarity in genetic elements across disparate secondary metabolite biosynthetic pathways renders imprecise current gene cluster targeting strategies and contributes sequence complexity resulting in partial genome coverage. Thus, it was necessary to use a dual-method approach of single-cell genomic sequencing based on multiple displacement amplification (MDA) and metagenomic library screening. Here, we report the identification of the putative apratoxin. A biosynthetic gene cluster, a potent cancer cell cytotoxin with promise for medicinal applications. The roughly 58 kb biosynthetic gene cluster is composed of 12 open reading frames and has a type I modular mixed polyketide synthase/nonribosomal peptide synthetase (PKS/NRPS) organization and features loading and off-loading domain architecture never previously described. Moreover, this work represents the first successful isolation of a complete biosynthetic gene cluster from Lyngbya bouillonii, a tropical marine cyanobacterium renowned for its production of diverse bioactive secondary metabolites.     

Duplication of partial spinosyn biosynthetic gene cluster in Saccharopolyspora spinosa enhances spinosyn production

Spinosyns, the secondary metabolites produced by Saccharopolyspora spinosa, are the active ingredients in a family of insect control agents. Most of the S. spinosa genes involved in spinosyn biosynthesis are found in a contiguous c. 74-kb cluster. To increase the spinosyn production through overexpression of their biosynthetic genes, part of its gene cluster (c. 18 kb) participating in the conversion of the cyclized polyketide to spinosyn was obtained by direct cloning via Red/ET recombination rather than by constructing and screening the genomic library. The resultant plasmid pUCAmT-spn was introduced into S. spinosa CCTCC M206084 from Escherichia coli S17-1 by conjugal transfer. The subsequent single-crossover homologous recombination caused a duplication of the partial gene cluster. Integration of this plasmid enhanced production of spinosyns with a total of 388 (± 25.0) mg L(-1) for spinosyns A and D in the exconjugant S. spinosa trans1 compared with 100 (± 7.7) mg L(-1) in the parental strain. Quantitative real time polymerase chain reaction analysis of three selected genes (spnH, spnI, and spnK) confirmed the positive effect of the overexpression of these genes on the spinosyn production. This study provides a simple avenue for enhancing spinosyn production. The strategies could also be used to improve the yield of other secondary metabolites.     

Mosaic origins of a complex chimeric mitochondrial gene in Silene vulgaris

Chimeric genes are significant sources of evolutionary innovation that are normally created when portions of two or more protein coding regions fuse to form a new open reading frame. In plant mitochondria astonishingly high numbers of different novel chimeric genes have been reported, where they are generated through processes of rearrangement and recombination. Nonetheless, because most studies do not find or report nucleotide variation within the same chimeric gene, evolution after the origination of these chimeric genes remains unstudied. Here we identify two alleles of a complex chimera in Silene vulgaris that are divergent in nucleotide sequence, genomic position relative to other mitochondrial genes, and expression patterns. Structural patterns suggest a history partially influenced by gene conversion between the chimeric gene and functional copies of subunit 1 of the mitochondrial ATP synthase gene (atp1). We identified small repeat structures within the chimeras that are likely recombination sites allowing generation of the chimera. These results establish the potential for chimeric gene divergence in different plant mitochondrial lineages within the same species. This result contrasts with the absence of diversity within mitochondrial chimeras found in crop species.     

Relationships of the Escherichia coli O157, O111, and O55 O-antigen gene clusters with those of Salmonella enterica and Citrobacter freundii, which express identical O antigens

Escherichia coli O157, Salmonella enterica O30, and Citrobacter freundii F90 have identical O-antigen structures, as do E. coli O55 and S. enterica O50. The O-antigen gene cluster sequences for E. coli O157 and E. coli O55 have been published, and the genes necessary for O-antigen biosynthesis have been identified, although transferase genes for glycosidic linkages are only generic and have not been allocated to specific linkages. We determined sequences for S. enterica O30 and C. freundii F90 O-antigen gene clusters and compared them to the sequence of the previously described E. coli O157 cluster. We also determined the sequence of the S. enterica O50 O-antigen gene cluster and compared it to the sequence of the previously described E. coli O55 cluster. For both the S. enterica O30-C. freundii F90-E. coli O157 group and the S. enterica O50-E. coli O55 group of O antigens, the gene clusters have identical or nearly identical organizations. The two sets of gene clusters had comparable overall levels of similarity in their genes, which were lower than the levels determined for housekeeping genes for these species, which were 55 to 65% for the genes encoding glycosyltransferases and O-antigen processing proteins and 75 to 93% for the nucleotide-sugar pathway genes. Nonetheless, the similarity of the levels of divergence in the five gene clusters required us to consider the possibility that the parent gene cluster for each structure was in the common ancestor of the species and that divergence is faster than expected for the common ancestor hypothesis. We propose that the identical O-antigen gene clusters originated from a common ancestor, and we discuss some possible explanations for the increased rate of divergence that is seen in these genes.     

Deletion of the Escherichia coli O14:K7 O antigen gene cluster

Escherichia coli O14:K7 is a rough strain, lacking a typical O antigen, in which the enterobacterial common antigen is attached to the lipopolysaccharide core. The rough phenotype was previously mapped to the O antigen gene cluster; however, the nature of the nonfunctional locus was not defined. In this study, we have shown that the O antigen gene cluster of an O14:K7 type strain (Su4411/41) was most likely deleted via homologous recombination between the GDP-mannose pathway genes (manB and manC) of the colanic acid and O antigen gene clusters. A similar recombination event has previously been inferred for the deletion of E. coli Sonnei chromosomal O antigen genes. Therefore, recombination between the GDP-mannose pathway genes provides a convenient mechanism for the deletion of O antigen genes, which may occur if the typical O antigen becomes redundant.