Deep sequencing of non-ribosomal peptide synthetases and polyketide synthases from the microbiomes of Australian marine sponges

The biosynthesis of non-ribosomal peptide and polyketide natural products is facilitated by multimodular enzymes that contain domains responsible for the sequential condensation of amino and carboxylic subunits. These conserved domains provide molecular targets for the discovery of natural products from microbial metagenomes. This study demonstrates the application of tag-encoded FLX amplicon pyrosequencing (TEFAP) targeting non-ribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) genes as a method for determining the identity and diversity of natural product biosynthesis genes. To validate this approach, we assessed the diversity of NRPS and PKS genes within the microbiomes of six Australian marine sponge species using both TEFAP and metagenomic whole-genome shotgun sequencing approaches. The TEFAP approach identified 100 novel ketosynthase (KS) domain sequences and 400 novel condensation domain sequences within the microbiomes of the six sponges. The diversity of KS domains within the microbiome of a single sponge species Scopalina sp. exceeded that of any previously surveyed marine sponge. Furthermore, this study represented the first to target the condensation domain from NRPS biosynthesis and resulted in the identification of a novel condensation domain lineage. This study highlights the untapped potential of Australian marine sponges for the isolation of novel bioactive natural products. Furthermore, this study demonstrates that TEFAP approaches can be applied to functional genes, involved in natural product biosynthesis, as a tool to aid natural product discovery. It is envisaged that this approach will be used across multiple environments, offering an insight into the biological processes that influence the production of secondary metabolites.     

Ketosynthase Domain Probes Identify Two Subclasses of Fungal Polyketide Synthase Genes

Analysis of fungal polyketide synthase gene sequences suggested that these might be divided into two subclasses, designated WA-type and MSAS-type. Two pairs of degenerate PCR primers (LC1 and LC2c, LC3 and LC5c) were designed for the amplification of ketosynthase domain fragments from fungal PKS genes in each of these subclasses. Both primer pairs were shown to amplify one or more PCR products from the genomes of a range of ascomycetous Deuteromycetes and Southern blot analysis confirmed that the products obtained with each pair of primers emanated from distinct genomic loci. PCR products obtained from Penicillium patulum and Aspergillus parasiticus with the LC1/2c primer pair and from Phoma sp. C2932 with both primer pairs were cloned and sequenced; the deduced protein sequences were highly homologous to the ketosynthase domains of other fungal PKS genes. Genes from which LC1/2c fragments were amplified (WA-type) were shown by a phylogenetic analysis to be closely related to fungal PKS genes involved in pigment and aflatoxin biosynthetic pathways, whereas the gene from which the LC3/5c fragment was amplified (MSAS-type) was shown to be closely related to genes encoding 6-methylsalicylic acid synthase (MSAS). The phylogenetic tree strongly supported the division of fungal PKS genes into two subclasses. The LC-series primers may be useful molecular tools to facilitate the cloning of novel fungal polyketide synthase genes.     

Diversity of type I polyketide synthase genes in the wood-decay fungus Xylaria sp. BCC 1067

Fungal type I polyketide (PK) compounds are highly valuable for medical treatment and extremely diverse in structure, partly because of the enzymatic activities of reducing domains in polyketide synthases (PKSs). We have cloned several PKS genes from the fungus Xylaria sp. BCC 1067, which produces two polyketides: depudecin (reduced PK) and 19,20-epoxycytochalasin Q (PK-nonribosomal peptide (NRP) hybrid). Two new degenerate primer sets, KA-series and XKS, were designed to amplify reducing PKS and PKS-NRP synthetase hybrid genes, respectively. Five putative PKS genes were amplified in Xylaria using KA-series primers and two more with the XKS primers. All seven are predicted to encode proteins homologous to highly reduced (HR)-type PKSs. Previously designed primers in LC-, KS-, and MT-series identified four additional PKS gene fragments. Selected PKS fragments were used as probes to identify PKS genes from the genomic library of this fungus. Full-length sequences for five PKS genes were obtained: pks12, pks3, pksKA1, pksMT, and pksX1. They are structurally diverse with 1-9 putative introns and products ranging from 2162 to 3654 amino acids in length. The finding of 11 distinct PKS genes solely by means of PCR cloning supports that PKS genes are highly diverse in fungi. It also indicates that our KA-series primers can serve as powerful tools to reveal the genetic potential of fungi in production of multiple types of HR PKs, which the conventional compound screening could underestimate.     

Towards Prediction of Metabolic Products of Polyketide Synthases: An In Silico Analysis

Sequence data arising from an increasing number of partial and complete genome projects is revealing the presence of the polyketide synthase (PKS) family of genes not only in microbes and fungi but also in plants and other eukaryotes. PKSs are huge multifunctional megasynthases that use a variety of biosynthetic paradigms to generate enormously diverse arrays of polyketide products that posses several pharmaceutically important properties. The remarkable conservation of these gene clusters across organisms offers abundant scope for obtaining novel insights into PKS biosynthetic code by computational analysis. We have carried out a comprehensive in silico analysis of modular and iterative gene clusters to test whether chemical structures of the secondary metabolites can be predicted from PKS protein sequences. Here, we report the success of our method and demonstrate the feasibility of deciphering the putative metabolic products of uncharacterized PKS clusters found in newly sequenced genomes. Profile Hidden Markov Model analysis has revealed distinct sequence features that can distinguish modular PKS proteins from their iterative counterparts. For iterative PKS proteins, structural models of iterative ketosynthase (KS) domains have revealed novel correlations between the size of the polyketide products and volume of the active site pocket. Furthermore, we have identified key residues in the substrate binding pocket that control the number of chain extensions in iterative PKSs. For modular PKS proteins, we describe for the first time an automated method based on crucial intermolecular contacts that can distinguish the correct biosynthetic order of substrate channeling from a large number of non-cognate combinatorial possibilities. Taken together, our in silico analysis provides valuable clues for formulating rules for predicting polyketide products of iterative as well as modular PKS clusters. These results have promising potential for discovery of novel natural products by genome mining and rational design of novel natural products.     

四种红树植物根际土壤微生物Ⅰ型和Ⅱ型PKS基因的检测与多样性分析

[目的]探究红树林土壤中聚酮合酶(Polyketide synthase,PKS)基因的多样性和新颖性.[方法]用Ⅰ型和Ⅱ型PKS基因酮基合成酶(Ketosynthase,KS)域的简并引物对海南清澜港红树林海莲、黄槿、银叶、老鼠簕4种红树根际土壤样品中DNA进行PCR扩增,之后利用PCR-限制性酶切片段多样性(PCR-RFLP)和测序分析法对Ⅰ型和Ⅱ型PKS基因的多样性进行探讨.[结果]对得到的72条Ⅰ型PKS基因的酮基合成酶(Ketosynthase,KS)域DNA序列进行PCR-RFLP分析,共得到51个可操作分类单元(Operational taxonomic unit,OTUs),其中37个OTUs为单克隆产生,没有明显的优势OTU.选取了26个代表不同OTU的克隆进行测序分析,这些序列与GenBank中已知序列的最大相似率均未超过85％. KS域氨基酸序列的系统发育分析显示,所得KS域来源广泛,包括蓝细菌门(Cyanobacteria)、变形杆菌门(Proteobacteria)、厚壁菌门(Firmicutes)、放线菌门(Actinobacteria)和一些未可培养细菌;对55条PKSⅡ基因KS域DNA序列的PCR-RFLP分析后共得到25个OTUs,有两个明显的优势OTUs,代表的克隆子数所占比例超过10％.[结论]PCR-RFLP分析表明红树林根际土壤中存在着丰富多样的Ⅰ型和Ⅱ型PKS基因,且前者多样性更高;低的序列相似度表明所获得的PKSⅠ基因KS域序列独特;系统发育分析表明得到的PKSⅠ基因来源广泛.     

Combinatorial polyketide biosynthesis by de novo design and rearrangement of modular polyketide synthase genes

Type I polyketide synthase (PKS) genes consist of modules approximately 3-6 kb long, which encode the structures of 2-carbon units in polyketide products. Alteration or replacement of individual PKS modules can lead to the biosynthesis of 'unnatural' natural products but existing techniques for this are time consuming. Here we describe a generic approach to the design of synthetic PKS genes where facile cassette assembly and interchange of modules and domains are facilitated by a repeated set of flanking restriction sites. To test the feasibility of this approach, we synthesized 14 modules from eight PKS clusters and associated them in 154 bimodular combinations spanning over 1.5-million bp of novel PKS gene sequences. Nearly half the combinations successfully mediated the biosynthesis of a polyketide in Escherichia coli, and all individual modules participated in productive bimodular combinations. This work provides a truly combinatorial approach for the production of polyketides.     

Phylogenetic analysis of polyketide synthase genes fromAspergillus ochraceus

A number of polyketide synthase gene sequences fromAspergillus ochraceus were isolated by both SSH-PCR and degenerate PCR. The deduced amino acid sequences of the corresponding clonedpks DNA fragments were then aligned with the amino acid sequences of other polyketide synthase enzymes. One of thesepks genes is essential for ochratoxin A biosynthesis (OTA-PKS). The OTA-PKS was most similar to methylsalicylic acid synthase (MSAS) type PKS proteins based on the alignment of the ketosynthase domains while if the acyl transferase domains were aligned it appeared to be more similar to PKS enzymes fromCochliobolus heterostrophus. The three PKS proteins identified by degenerate PCR were all from different PKS types, one was a MSAS type enzyme, the second was similar to the PKS proteins involved in lovastatin biosynthesis while the third was not similar to any of the other phylogenetic groupings. Data is presented which suggests that the use of phylogenetic analysis to predict the function of PKS proteins/genes is likely to be significantly enhanced by analyzing more than one domain of the protein. Presented at the EU-USA Bilateral Workshop on Toxigenic Fungi & Mycotoxins, New Orleans, USA, July 5–7, 2005 Financial support: Irish Government under the National Development Plan 2000–2006     

Polyketide synthase pathways identified from a metagenomic library are derived from soil Acidobacteria

Polyketides are structurally diverse secondary metabolites, many of which have antibiotic or anticancer activity. Type I modular polyketide synthase (PKS) genes are typically large and encode repeating enzymatic domains that elongate and modify the nascent polyketide chain. A fosmid metagenomic library constructed from an agricultural soil was arrayed and the macroarray was screened for the presence of conserved ketosynthase [β-ketoacyl synthase (KS)] domains, enzymatic domains present in PKSs. Thirty-four clones containing KS domains were identified by Southern hybridization. Many of the KS domains contained within metagenomic clones shared significant similarity to PKS or nonribosomal peptide synthesis genes from members of the Cyanobacteria or the Proteobacteria phyla. However, analysis of complete clone insert sequences indicated that the blast analysis for KS domains did not reflect the true phylogenetic origin of many of these metagenomic clones that had a %G+C content and significant sequence similarity to genes from members of the phylum Acidobacteria. This conclusion of an Acidobacteria origin for several clones was further supported by evidence that cultured soil Acidobacteria from different subdivisions have genetic loci closely related to PKS domains contained within metagenomic clones, suggesting that Acidobacteria may be a source of novel polyketides. This study also demonstrates the utility of combining data from culture-dependent and -independent investigations in expanding our collective knowledge of microbial genomic diversity.     

Biosynthesis of Dictyostelium discoideum differentiation-inducing factor by a hybrid type I fatty acid-type III polyketide synthase

Differentiation-inducing factors (DIFs) are well known to modulate formation of distinct communal cell types from identical Dictyostelium discoideum amoebas, but DIF biosynthesis remains obscure. We report complimentary in vivo and in vitro experiments identifying one of two approximately 3,000-residue D. discoideum proteins, termed 'steely', as responsible for biosynthesis of the DIF acylphloroglucinol scaffold. Steely proteins possess six catalytic domains homologous to metazoan type I fatty acid synthases (FASs) but feature an iterative type III polyketide synthase (PKS) in place of the expected FAS C-terminal thioesterase used to off load fatty acid products. This new domain arrangement likely facilitates covalent transfer of steely N-terminal acyl products directly to the C-terminal type III PKS active sites, which catalyze both iterative polyketide extension and cyclization. The crystal structure of a steely C-terminal domain confirms conservation of the homodimeric type III PKS fold. These findings suggest new bioengineering strategies for expanding the scope of fatty acid and polyketide biosynthesis.     

Identification of unique type II polyketide synthase genes in soil

Many bacteria, particularly actinomycetes, are known to produce secondary metabolites synthesized by polyketide synthases (PKS). Bacterial polyketides are a particularly rich source of bioactive molecules, many of which are of potential pharmaceutical relevance. To directly access PKS gene diversity from soil, we developed degenerate PCR primers for actinomycete type II KS(alpha) (ketosynthase) genes. Twenty-one soil samples were collected from diverse sources in New Jersey, and their bacterial communities were compared by terminal restriction fragment length polymorphism (TRFLP) analysis of PCR products generated using bacterial 16S rRNA gene primers (27F and 1525R) as well as an actinomycete-specific forward primer. The distribution of actinomycetes was highly variable but correlated with the overall bacterial species composition as determined by TRFLP. Two samples were identified to contain a particularly rich and unique actinomycete community based on their TRFLP patterns. The same samples also contained the greatest diversity of KS(alpha) genes as determined by TRFLP analysis of KS(alpha) PCR products. KS(alpha) PCR products from these and three additional samples with interesting TRFLP pattern were cloned, and seven novel clades of KS(alpha) genes were identified. Greatest sequence diversity was observed in a sample containing a moderate number of peaks in its KS(alpha) TRFLP. The nucleotide sequences were between 74 and 81% identical to known sequences in GenBank. One cluster of sequences was most similar to the KS(alpha) involved in ardacin (glycopeptide antibiotic) production by Kibdelosporangium aridum. The remaining sequences showed greatest similarity to the KS(alpha) genes in pathways producing the angucycline-derived antibiotics simocyclinone, pradimicin, and jasomycin.     

Evolution of ketosynthase domains of polyketide synthase genes in the Cladonia chlorophaea species complex (Cladoniaceae)

Lichen-forming fungi synthesize a diversity of polyketides, but only a few non-reducing polyketide synthase (PKS) genes from a lichen-forming fungus have been linked with a specific polyketide. While it is a challenge to link the large number of PKS paralogs in fungi with specific products, it might be expected that the PKS paralogs from closely related species would be similar because of recent evolutionary divergence. The objectives of this study were to reconstruct a PKS gene phylogeny of the Cladonia chlorophaea species complex based on the ketosynthase domain, a species phylogeny of the complex, and to explore the presence of PKS gene paralogs among members of the species complex. DNA was isolated from 51 individuals of C. chlorophaea and allies to screen for the presence of 13 PKS paralogs. A 128 sequence PKS gene phylogeny using deduced amino acid sequences estimated from the 13 PKS paralogs and sequences subjected to BLASTx comparisons showed losses of each of two PKS domains (reducing and methylation). This research provided insight into the evolution of PKS genes in the C. chlorophaea group, species evolution in the group, and it identified potential directions for further investigation of polyketide synthesis in the C. chlorophaea species complex.     

Diversity of Polyketide Synthase (PKS) genes in metagenomic community of freshwater sponge Lubomirskia baicalensis

The screening of metagenomic DNA of the microbial community associated with the Baikalian sponge Lubomirskia baicalensis was performed in order to investigate the presence of polyketide synthase (PKS) genes. PKS enzyme systems take part in the synthesis of a great number of biologically active substances. The cloning and sequencing of amplified products of the ketosynthase domain section of the PKS gene cluster revealed 15 fragments of PKS genes with amino acid sequences differing from each other by 35?C65%. A BLASTX analysis showed that all of these sequences belong to KS domains identified in various groups of microorganisms, i.e., Alpha-, Beta-, and Deltaproteobacteria; Verrucomicrobia; Cyanobacteria; and Chlorophyta. Some sequences were related to genes that participate in the biosynthesis of curacin A (CurI, CurJ), stigmatellin (StiC, StiG), nostophycin (NpnB), and cryptophycin (CrpB). The homology of the found sequences with those of the EMBL database lies in the range of 50?C82%, which indicates that the freshwater sponge community contains genes that encode new, not yet studied polyketide substances of potential biotechnological significance.     

PCR Detection of Type I Polyketide Synthase Genes in Myxobacteria

The diversity of type I modular polyketide synthase (PKS) was explored by PCR amplification of DNA encoding ketosynthase and acyltransferase domains in myxobacteria. The sequencing of the amplicons revealed that many PKS genes were distantly related to the published sequences. Thus, myxobacteria may be excellent resources for novel and diverse polyketides.     

Metagenomic analysis reveals diverse polyketide synthase gene clusters in microorganisms associated with the marine sponge Discodermia dissoluta

Sponge-associated bacteria are thought to produce many novel bioactive compounds, including polyketides. PCR amplification of ketosynthase domains of type I modular polyketide synthases (PKS) from the microbial community of the marine sponge Discodermia dissoluta revealed great diversity and a novel group of sponge-specific PKS ketosynthase domains. Metagenomic libraries totaling more than four gigabases of bacterial genomes associated with this sponge were screened for type I modular PKS gene clusters. More than 90% of the clones in total sponge DNA libraries represented bacterial DNA inserts, and 0.7% harbored PKS genes. The majority of the PKS hybridizing clones carried small PKS clusters of one to three modules, although some clones encoded large multimodular PKSs (more than five modules). The most abundant large modular PKS appeared to be encoded by a bacterial symbiont that made up < 1% of the sponge community. Sequencing of this PKS revealed 14 modules that, if expressed and active, is predicted to produce a multimethyl-branched fatty acid reminiscent of mycobacterial lipid components. Metagenomic libraries made from fractions enriched for unicellular or filamentous bacteria differed significantly, with the latter containing numerous nonribosomal peptide synthetase (NRPS) and mixed NRPS-PKS gene clusters. The filamentous bacterial community of D. dissoluta consists mainly of Entotheonella spp., an unculturable sponge-specific taxon previously implicated in the biosynthesis of bioactive peptides.     

A computational screen for type I polyketide synthases in metagenomics shotgun data

Background

Polyketides are a diverse group of biotechnologically important secondary metabolites that are produced by multi domain enzymes called polyketide synthases (PKS).

Methodology/Principal Findings

We have estimated frequencies of type I PKS (PKS I) – a PKS subgroup – in natural environments by using Hidden-Markov-Models of eight domains to screen predicted proteins from six metagenomic shotgun data sets. As the complex PKS I have similarities to other multi-domain enzymes (like those for the fatty acid biosynthesis) we increased the reliability and resolution of the dataset by maximum-likelihood trees. The combined information of these trees was then used to discriminate true PKS I domains from evolutionary related but functionally different ones. We were able to identify numerous novel PKS I proteins, the highest density of which was found in Minnesota farm soil with 136 proteins out of 183,536 predicted genes. We also applied the protocol to UniRef database to improve the annotation of proteins with so far unknown function and identified some new instances of horizontal gene transfer.

Conclusions/Significance

The screening approach proved powerful in identifying PKS I sequences in large sequence data sets and is applicable to many other protein families.     

Molecular characterization of a Dictyostelium discoideum gene encoding a multifunctional enzyme of the pyrimidine pathway

We have isolated and characterized a Dictyostelium discoideum gene (PYR1-3) encoding a multifunctional protein that carries the three first enzymatic activities of the de novo pyrimidine biosynthetic pathway. The PYR1-3 gene is adjacent to another gene of the pyrimidine biosynthetic pathway (PYR4); the two genes are separated by a 1.5-kb non-coding sequence and transcribed divergently. The PYR1-3 gene is transcribed to form a 7.5-kb polyadenylated mRNA. As with the other genes of the pyrimidine biosynthetic pathway, the PYR1-3 mRNA level is high during growth and decreases sharply during development. We have determined the nucleotide sequence of 63% of the coding region of the PYR1-3 gene. We have identified the activities of the protein encoded by the D. discoideum PYR1-3 gene by comparison of amino acid sequences with the products of genes of known function. The PYR1-3 gene contains four distinct regions that probably correspond to four domains in the protein. From the NH2 extremity to the COOH extremity, these domains are: glutamine amidotransferase, carbamoylphosphate synthetase, dihydroorotase and aspartate transcarbamylase. This organization is identical to the one found in the rudimentary gene of Drosophila. The evolutionary implications of this finding are discussed.     

Targeting polyketide synthase gene pool within actinomycetes: new degenerate primers

Natural products provide a unique element of molecular diversity and biological functionality and they are still indispensable for drug discovery. The polyketides, comprising a large and structurally diverse family of bioactive natural products, have been isolated from a group of mycelia-forming Gram-positive microorganisms, the actinomycetes. Relatively high amino acid sequence identity of the actinomycetes type I polyketide synthases (PKSs) was used to design three degenerate primer pairs for homology-based PCR detection of novel PKS genes, with particular interest into PKSs involved in biosynthesis of immunosuppressive-like metabolites. The stepdown PCR method, described here, enables fast insight into the PKS arsenal within actinomycetes. Designed primers and stepdown PCR were applied for the analysis of two natural isolates, Streptomyces sp. strains NP13 and MS405. Sequence analysis of chosen clones revealed the presence of two distinctive sequences in strain Streptomyces sp. NP13, but only one of these showed homology to PKS-related sequences. On analysing PCR amplicons derived from Streptomyces sp. strain MS405, three different PKS-related sequences were identified demonstrating a potential of designed primers to target PKS gene pool within single organism.     

Characterization of the biosynthetic gene cluster for the antifungal polyketide soraphen A from Sorangium cellulosum So ce26

A genomic DNA region of over 80 kb that contains the complete biosynthetic gene cluster for the synthesis of the antifungal polyketide metabolite soraphen A was cloned from Sorangium cellulosum So ce26. The nucleotide sequence of the soraphen A gene region, including 67,523 bp was determined. Examination of this sequence led to the identification of two adjacent type I polyketide synthase (PKS) genes that encode the soraphen synthase. One of the soraphen A PKS genes includes three biosynthetic modules and the second contains five additional modules for a total of eight. The predicted substrate specificities of the acyltransferase (AT) domains, as well as the reductive loop domains identified within each module, are consistent with expectations from the structure of soraphen A. Genes were identified in the regions flanking the two soraphen synthase genes that are proposed to have roles in the biosynthesis of soraphen A. Downstream of the soraphen PKS genes is an O-methyltransferase (OMT) gene. Upstream of the soraphen PKS genes there is a gene encoding a reductase and a group of genes that are postulated to have roles in the synthesis of methoxymalonyl-acyl carrier protein (ACP). This unusual extender unit is proposed to be incorporated in two positions of the soraphen polyketide chain. One of the genes in this group contains distinct domains for an AT, an ACP, and an OMT.