A Systematic Computational Analysis of Biosynthetic Gene Cluster Evolution: Lessons for Engineering Biosynthesis

Bacterial secondary metabolites are widely used as antibiotics, anticancer drugs, insecticides and food additives. Attempts to engineer their biosynthetic gene clusters (BGCs) to produce unnatural metabolites with improved properties are often frustrated by the unpredictability and complexity of the enzymes that synthesize these molecules, suggesting that genetic changes within BGCs are limited by specific constraints. Here, by performing a systematic computational analysis of BGC evolution, we derive evidence for three findings that shed light on the ways in which, despite these constraints, nature successfully invents new molecules: 1) BGCs for complex molecules often evolve through the successive merger of smaller sub-clusters, which function as independent evolutionary entities. 2) An important subset of polyketide synthases and nonribosomal peptide synthetases evolve by concerted evolution, which generates sets of sequence-homogenized domains that may hold promise for engineering efforts since they exhibit a high degree of functional interoperability, 3) Individual BGC families evolve in distinct ways, suggesting that design strategies should take into account family-specific functional constraints. These findings suggest novel strategies for using synthetic biology to rationally engineer biosynthetic pathways.     

Valinomycin Biosynthetic Gene Cluster in Streptomyces: Conservation,Ecology and Evolution

Many Streptomyces strains are known to produce valinomycin (VLM) antibiotic and the VLM biosynthetic gene cluster (vlm) has been characterized in two independent isolates. Here we report the phylogenetic relationships of these strains using both parsimony and likelihood methods, and discuss whether the vlm gene cluster shows evidence of horizontal transmission common in natural product biosynthetic genes. Eight Streptomyces strains from around the world were obtained and sequenced for three regions of the two large nonribosomal peptide synthetase genes (vlm1 and vlm2) involved in VLM biosynthesis. The DNA sequences representing the vlm gene cluster are highly conserved among all eight environmental strains. The geographic distribution pattern of these strains and the strict congruence between the trees of the two vlm genes and the housekeeping genes, 16S rDNA and trpB, suggest vertical transmission of the vlm gene cluster in Streptomyces with no evidence of horizontal gene transfer. We also explored the relationship of the sequence of vlm genes to that of the cereulide biosynthetic genes (ces) found in Bacillus cereus and found them highly divergent from each other at DNA level (genetic distance values≥95.6%). It is possible that the vlm gene cluster and the ces gene cluster may share a relatively distant common ancestor but these two gene clusters have since evolved independently.     

Organization of the Biosynthetic Gene Cluster for the Polyketide Antitumor Macrolide,Pladienolide, in Streptomyces platensis Mer-11107

Pladienolides are novel 12-membered macrolides produced by Streptomyces platensis Mer-11107. They show strong antitumor activity and are a potential lead in the search for novel antitumor agents. We sequenced the 65-kb region covering the biosynthetic gene cluster, and found four polyketide synthase genes (pldAI-pldAIV) composed of 11 modules, three genes involved in post-modifications (pldB-D), and a luxR-family regulatory gene (pldR). The thioesterase domain of pldAIV was more dissimilar to that of polyketide synthase systems synthesizing 12/14-membered macrolide polyketides than to that of systems synthesizing other cyclic polyketides. The pldB gene was identified as a 6-hydroxylase belonging to a cytochrome P450 of the CYP107 family. This was clarified by a disruption experiment on pldB, in which the disruptant produced 6-dehydroxy pladienolide B. Two genes located downstream of pldB, designated pldC and pldD, are thought to be a probable genes for 7-O-acetylase and 18, 19-epoxydase respectively.     

The Yersiniabactin Biosynthetic Gene Cluster of Yersinia enterocolitica: Organization and Siderophore-Dependent Regulation

The ability to synthesize and uptake the Yersinia siderophore yersiniabactin is a hallmark of the highly pathogenic, mouse-lethal species Yersinia pestis, Y. pseudotuberculosis, and Y. enterocolitica 1B. We have identified four genes, irp1, irp3, irp4, and irp5, on a 13-kb chromosomal DNA fragment of Y. enterocolitica O8, WA-314. These genes constitute the yersiniabactin biosynthetic gene cluster together with the previously defined irp2. The irp1 gene consists of 9,486 bp capable of encoding a 3,161-amino-acid high-molecular-weight protein 1 (HMWP1) polypeptide with a predicted mass of 384.6 kDa. The first 3,000 bp of irp1 show similarity to the corresponding regions of the polyketide synthase genes of Bacillus subtilis and Streptomyces antibioticus. The remaining part of irp1 is most similar to irp2, encoding HMWP2, which might be the reason for immunological cross-reactivity of the two polypeptides. Irp4 was found to have 41.7% similarity to thioesterase-like protein of the anguibactin biosynthetic genes of Vibrio anguillarum. Irp5 shows 41% similarity to EntE, the 2,3-dihydroxybenzoic acid-activating enzyme utilized in enterobactin synthesis of Escherichia coli. Irp4 and Irp5 are nearly identical to YbtT and YbtE, recently identified in Y. pestis. irp3 has no similarity to any known gene. Inactivation of either irp1 or irp2 abrogates yersiniabactin synthesis. Mutations in irp1 or fyuA (encoding yersiniabactin/pesticin receptor) result in downregulation of irp2 that can be upregulated by the addition of yersiniabactin. A FyuA-green fluorescent protein translational fusion was downregulated in an irp1 mutant. Upregulation was achieved by addition of yersiniabactin but not desferal, pesticin, or pyochelin, which indicates high specificity of the FyuA receptor and autoregulation of genes involved in synthesis and uptake of yersiniabactin.     

Isolation, Characterization, and Heterologous Expression of the Novel Lantibiotic Epicidin 280 and Analysis of Its Biosynthetic Gene Cluster

Epicidin 280 is a novel type A lantibiotic produced by Staphylococcus epidermidis BN 280. During C₁₈ reverse-phase high-performance liquid chromatography two epicidin 280 peaks were obtained; the two compounds had molecular masses of 3,133 ± 1.5 and 3,136 ± 1.5 Da, comparable antibiotic activities, and identical amino acid compositions. Amino acid sequence analysis revealed that epicidin 280 exhibits 75% similarity to Pep5. The strains that produce epicidin 280 and Pep5 exhibit cross-immunity, indicating that the immunity peptides cross-function in antagonization of both lantibiotics. The complete epicidin 280 gene cluster was cloned and was found to comprise at least five open reading frames (eciI, eciA, eciP, eciB, and eciC, in that order). The proteins encoded by these open reading frames exhibit significant sequence similarity to the biosynthetic proteins of the Pep5 operon of Staphylococcus epidermidis 5. A gene for an ABC transporter, which is present in the Pep5 gene cluster but is necessary only for high yields (G. Bierbaum, M. Reis, C. Szekat, and H.-G. Sahl, Appl. Environ. Microbiol. 60:4332–4338, 1994), was not detected. Instead, upstream of the immunity gene eciI we found an open reading frame, eciO, which could code for a novel lantibiotic modification enzyme involved in reduction of an N-terminally located oxopropionyl residue. Epicidin 280 produced by the heterologous host Staphylococcus carnosus TM 300 after introduction of eciIAPBC (i.e., no eciO was present) behaved homogeneously during reverse-phase chromatography.     

Cloning and Heterologous Expression of the Vibrioferrin Biosynthetic Gene Cluster from a Marine Metagenomic Library

A biosynthetic gene cluster of siderophore consisting of five open reading frames (ORFs) was cloned by functional screening of a metagenomic library constructed from tidal-flat sediment. Expression of the cloned biosynthetic genes in Escherichia coli led to the production of vibrioferrin, a siderophore originally reported for the marine bacterium Vibrio parahaemolyticus. To the best of our knowledge, this is the first example of heterologous production of a siderophore by biosynthetic genes cloned from a metagenomic library. The cloned cluster was one of the largest of the clusters obtained by functional screening. In this study, we demonstrated and extended the possibility of function-based metagenomic research.     

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.     

来源卤化二型聚酮类生物合成基因簇的两种关键功能基因鉴定与表达

为催化卤化产物,挖掘具有生物活性的新卤化物提供新的方法途径,拟构建卤化酶原核表达载体催化天然卤化物的卤化过程.对链霉菌Streptomyces sp.FJS31-2基因组中一个II型PKS类产物zunyimycin生物合成基因簇进行系列分析,针对其中的后修饰酶卤化酶基因和单加氧酶基因进行原核表达纯化.通过分子生物学手段...     

木霉菌T23胶毒素合成基因的生物信息学分析与克隆

胶毒素是生防木霉菌重要的次生代谢产物之一。本研究以生防木霉菌T23为供试材料,旨在通过生物信息学技术及表达分析,挖掘木霉菌T23中胶毒素合成候选基因,探索木霉菌胶毒素合成的分子调控机制,可为新型生物农药的开发及应用提供理论依据。研究表明,木霉菌T23中胶毒素合成候选基因簇全长28 kb,簇内包含了8个基因,分别与烟曲霉胶毒素合成基因簇内的gliP、gliC、gliN、gliK、gliI、gliG、gliF、gliM高度同源。提取培养2 d、3 d、4 d、5 d的木霉菌T23菌丝的RNA,通过半定量RT-PCR技术探索各候选基因在木霉菌T23不同生长时期的表达情况,显示各基因在不同生长时期均有表达,属于组成型表达基因。成功克隆得到木霉菌T23中的gliP-T23基因并完成基因结构分析,该基因全长6 339 bp,由4个外显子和3个内含子组成,为后续的基因功能验证提供基础。     

吸水链霉菌17997格尔德霉素部分生物合成基因簇的克隆和分析

吸水链霉菌17997(Streptomyceshy groscopicus17997)是我所从中国云南土壤中分离到的格尔德霉素(geldanamycin,GDM)产生菌,GDM具有良好的抗肿瘤和抗病毒活性,但其肝毒性和水溶性差的缺点限制了其在临床上的应用。为了实现对GDM结构的生物学改造,首先要获得GDM的生物合成基因。根据GDM后修饰基因——氨甲酰基转移酶基因(gdmN)的保守序列筛选S.hygroscopicus17997的柯斯质粒基因组文库,共获得6个阳性克隆,选择CT-4阳性柯斯质粒进行亚克隆和测序,又通过PCR延伸的方法获得了与CT4连锁的将近5kb的外源序列,共获得28.356kb的外源DNA序列,其中包含了13个可能阅读框架,通过同源比较证实该序列与S.hygroscopicusNRRL3602中的GDM生物合成基因有很高的同源性。为进一步研究GDM生物合成基因的功能,并通过组合生物学的方法改造GDM的结构奠定了基础。     

Gibberellin Biosynthetic Pathway inGibberella fujikuroi:Evidence for a Gene Cluster

Differential screening of aGibberella fujikuroicDNA library was used to successfully clone and identify genes involved in the pathway of gibberellin biosynthesis. Several cDNA clones that hybridized preferentially to a cDNA probe prepared from mycelium induced for gibberellin production were isolated and characterized. The deduced amino acid sequences of two (identical) clones contained the conserved heme-binding motif of cytochrome P450 monooxygenases (FXXGXXXCXG). One of these cDNA fragments was used as a homologous probe for the screening of a genomic library. A hybridizing 6.7-kb genomicSalI fragment was cloned into pUC19. The sequencing of this clone revealed that a second cytochrome P450 monooxygenase gene was closely linked to the first one. Since at least four cytochrome P450 monooxygenase-catalyzed steps are involved in the synthesis of gibberellins, chromosome walking was performed to find a further gene of this family or other genes involved in gibberellin pathway. Next to the two P450 monooxygenase genes, a putative geranylgeranyl diphosphate synthase gene, the copalyl diphosphate synthase gene, which is the first specific gene of the gibberellin pathway, and a third P450 monooxygenase gene were identified. These results suggest that at least some of the genes involved in the biosynthesis of gibberellins are closely linked in a gene cluster inG. fujikuroi,as has been recently found for other “dispensable” pathways in fungi.     

Identification of the Biosynthetic Gene Cluster and an Additional Gene for Resistance to the Antituberculosis Drug Capreomycin

Capreomycin (CMN) belongs to the tuberactinomycin family of nonribosomal peptide antibiotics that are essential components of the drug arsenal for the treatment of multidrug-resistant tuberculosis. Members of this antibiotic family target the ribosomes of sensitive bacteria and disrupt the function of both subunits of the ribosome. Resistance to these antibiotics in Mycobacterium species arises due to mutations in the genes coding for the 16S or 23S rRNA but can also arise due to mutations in a gene coding for an rRNA-modifying enzyme, TlyA. While Mycobacterium species develop resistance due to alterations in the drug target, it has been proposed that the CMN-producing bacterium, Saccharothrix mutabilis subsp. capreolus, uses CMN modification as a mechanism for resistance rather than ribosome modification. To better understand CMN biosynthesis and resistance in S. mutabilis subsp. capreolus, we focused on the identification of the CMN biosynthetic gene cluster in this bacterium. Here, we describe the cloning and sequence analysis of the CMN biosynthetic gene cluster from S. mutabilis subsp. capreolus ATCC 23892. We provide evidence for the heterologous production of CMN in the genetically tractable bacterium Streptomyces lividans 1326. Finally, we present data supporting the existence of an additional CMN resistance gene. Initial work suggests that this resistance gene codes for an rRNA-modifying enzyme that results in the formation of CMN-resistant ribosomes that are also resistant to the aminoglycoside antibiotic kanamycin. Thus, S. mutabilis subsp. capreolus may also use ribosome modification as a mechanism for CMN resistance.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Occurrence of UV-Absorbing, Mycosporine-Like Compounds among Cyanobacterial Isolates and an Estimate of Their Screening Capacity

A survey of 20 strains of cyanobacteria (belonging to 13 genera) isolated from habitats exposed to strong insolation revealed that 13 strains contained one or more water-soluble, UV-absorbing, mycosporine amino acid (MAA)-like compounds. Some of the compounds were identical in several strains. In all, 13 distinct compounds were found. The UV absorption spectra of MAAs complemented well that of the extracellular sunscreen pigment scytonemin, which many of the strains also produced. Even though the specific MAA contents were variable among strains, they were invariably higher when the cultures were grown with UV radiation than when it was absent. In five strains tested, the MAA complement accumulated as a solute in the cytoplasmic cell fraction. The sunscreen capacities of MAA and scytonemin and their combined capacity were estimated for each strain and condition on the basis of the specific contents, cell size, and cellular location of the compounds. The estimates suggested that significant, albeit not complete, protection from UV photodamage could be gained from the possession of either MAA or scytonemin but especially from simultaneous screening by both types of compounds.