The Natural Product Domain Seeker version 2 (NaPDoS2) webtool relates ketosynthase phylogeny to biosynthetic function.

The Natural Product Domain Seeker (NaPDoS) webtool detects and classifies ketosynthase (KS) and condensation domains from genomic, metagenomic, and amplicon sequence data. Unlike other tools, a phylogeny-based classification scheme is used to make broader predictions about the polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) genes in which these domains are found. NaPDoS is particularly useful for the analysis of incomplete biosynthetic genes or gene clusters, as are often observed in poorly assembled genomes and metagenomes, or when loci are not clustered, as in eukaryotic genomes. To help support the growing interest in sequence-based analyses of natural product biosynthetic diversity, here we introduce version 2 of the webtool, NaPDoS2, available at http://napdos.ucsd.edu/napdos2. This update includes the addition of 1417 KS sequences, representing a major expansion of the taxonomic and functional diversity represented in the webtool database. The phylogeny-based KS classification scheme now recognizes 41 class and subclass assignments, including new type II PKS subclasses. Workflow modifications accelerate run times, allowing larger datasets to be analyzed. In addition, default parameters were established using statistical validation tests to maximize KS detection and classification accuracy while minimizing false positives. We further demonstrate the applications of NaPDoS2 to assess PKS biosynthetic potential using genomic, metagenomic, and PCR amplicon datasets. These examples illustrate how NaPDoS2 can be used to predict biosynthetic potential and detect genes involved in the biosynthesis of specific structure classes or new biosynthetic mechanisms.

Structure and Candidate Biosynthetic Gene Cluster of a Manumycin-Type Metabolite from Salinispora pacifica.

A new manumycin-type natural product named pacificamide (1) and its candidate biosynthetic gene cluster (pac) were discovered from the marine actinobacterium Salinispora pacifica CNT-855. The structure of the compound was determined using NMR, electronic circular dichroism, and bioinformatic predictions. The pac gene cluster is unique to S. pacifica and found in only two of the 119 Salinispora genomes analyzed across nine species. Comparative analyses of biosynthetic gene clusters encoding the production of related manumycin-type compounds revealed genetic differences in accordance with the unique pacificamide structure. Further queries of manumycin-type gene clusters from public databases revealed their limited distribution across the phylum Actinobacteria and orphan diversity that suggests additional products remain to be discovered in this compound class. Production of the known metabolite triacsin D is also reported for the first time from the genus Salinispora. This study adds two classes of compounds to the natural product collective isolated from the genus Salinispora, which has proven to be a useful model for natural product research.

Extending the Salinilactone Family.

Five new members of the salinilactone family, salinilactones D-H, are reported. These bicyclic lactones are produced by Salinispora bacteria and display extended or shortened alkyl side chains relative to the recently reported salinilactones A-C. They were identified by GC/MS, gas chromatographic retention index, and comparison with synthetic samples. We further investigated the occurrence of salinilactones across six newly proposed Salinispora species to gain insight into how compound production varies among taxa. The growth-inhibiting effect of this compound family on multiple biological systems including non-Salinispora actinomycetes was analyzed. Additionally, we found strong evidence for significant cytotoxicity of the title compounds.

Experimental Evolution of Escherichia coli K-12 at High pH and with RpoS Induction.

Experimental evolution of Escherichia coli K-12 W3110 by serial dilutions for 2,200 generations at high pH extended the range of sustained growth from pH 9.0 to pH 9.3. pH 9.3-adapted isolates showed mutations in DNA-binding regulators and envelope proteins. One population showed an IS1 knockout of phoB (encoding the positive regulator of the phosphate regulon). A phoB::kanR knockout increased growth at high pH. phoB mutants are known to increase production of fermentation acids, which could enhance fitness at high pH. Mutations in pcnB [poly(A) polymerase] also increased growth at high pH. Three out of four populations showed deletions of torI, an inhibitor of TorR, which activates expression of torCAD (trimethylamine N-oxide respiration) at high pH. All populations showed point mutations affecting the stationary-phase sigma factor RpoS, either in the coding gene or in genes for regulators of RpoS expression. RpoS is required for survival at extremely high pH. In our microplate assay, rpoS deletion slightly decreased growth at pH 9.1. RpoS protein accumulated faster at pH 9 than at pH 7. The RpoS accumulation at high pH required the presence of one or more antiadaptors that block degradation (IraM, IraD, and IraP). Other genes with mutations after high-pH evolution encode regulators, such as those encoded by yobG (mgrB) (PhoPQ regulator), rpoN (nitrogen starvation sigma factor), malI, and purR, as well as envelope proteins, such as those encoded by ompT and yahO Overall, E. coli evolution at high pH selects for mutations in key transcriptional regulators, including phoB and the stationary-phase sigma factor RpoS.IMPORTANCEEscherichia coli in its native habitat encounters high-pH stress such as that of pancreatic secretions. Experimental evolution over 2,000 generations showed selection for mutations in regulatory factors, such as deletion of the phosphate regulator PhoB and mutations that alter the function of the global stress regulator RpoS. RpoS is induced at high pH via multiple mechanisms.

Acid Evolution of Escherichia coli K-12 Eliminates Amino Acid Decarboxylases and Reregulates Catabolism.

Acid-adapted strains of Escherichia coli K-12 W3110 were obtained by serial culture in medium buffered at pH 4.6 (M. M. Harden, A. He, K. Creamer, M. W. Clark, I. Hamdallah, K. A. Martinez, R. L. Kresslein, S. P. Bush, and J. L. Slonczewski, Appl Environ Microbiol 81:1932-1941, 2015, https://doi.org/10.1128/AEM.03494-14). Revised genomic analysis of these strains revealed insertion sequence (IS)-driven insertions and deletions that knocked out regulators CadC (acid induction of lysine decarboxylase), GadX (acid induction of glutamate decarboxylase), and FNR (anaerobic regulator). Each acid-evolved strain showed loss of one or more amino acid decarboxylase systems, which normally help neutralize external acid (pH 5 to 6) and increase survival in extreme acid (pH 2). Strains from populations B11, H9, and F11 had an IS5 insertion or IS-mediated deletion in cadC, while population B11 had a point mutation affecting the arginine activator adiY The cadC and adiY mutants failed to neutralize acid in the presence of exogenous lysine or arginine. In strain B11-1, reversion of an rpoC (RNA polymerase) mutation partly restored arginine-dependent neutralization. All eight strains showed deletion or downregulation of the Gad acid fitness island. Strains with the Gad deletion lost the ability to produce GABA (gamma-aminobutyric acid) and failed to survive extreme acid. Transcriptome sequencing (RNA-seq) of strain B11-1 showed upregulated genes for catabolism of diverse substrates but downregulated acid stress genes (the biofilm regulator ariR, yhiM, and Gad). Other strains showed downregulation of H2 consumption mediated by hydrogenases (hya and hyb) which release acid. Strains F9-2 and F9-3 had a deletion of fnr and showed downregulation of FNR-dependent genes (dmsABC, frdABCD, hybABO, nikABCDE, and nrfAC). Overall, strains that had evolved in buffered acid showed loss or downregulation of systems that neutralize unbuffered acid and showed altered regulation of catabolism.IMPORTANCE Experimental evolution of an enteric bacterium under a narrow buffered range of acid pH leads to loss of genes that enhance fitness above or below the buffered pH range, including loss of enzymes that may raise external pH in the absence of buffer. Prominent modes of evolutionary change involve IS-mediated insertions and deletions that knock out key regulators. Over generations of acid stress, catabolism undergoes reregulation in ways that differ for each evolving strain.

Benzoate- and Salicylate-Tolerant Strains of Escherichia coli K-12 Lose Antibiotic Resistance during Laboratory Evolution.

Escherichia coli K-12 W3110 grows in the presence of membrane-permeant organic acids that can depress cytoplasmic pH and accumulate in the cytoplasm. We conducted experimental evolution by daily diluting cultures in increasing concentrations of benzoic acid (up to 20 mM) buffered at external pH 6.5, a pH at which permeant acids concentrate in the cytoplasm. By 2,000 generations, clones isolated from evolving populations showed increasing tolerance to benzoate but were sensitive to chloramphenicol and tetracycline. Sixteen clones grew to stationary phase in 20 mM benzoate, whereas the ancestral strain W3110 peaked and declined. Similar growth occurred in 10 mM salicylate. Benzoate-evolved strains grew like W3110 in the absence of benzoate, in media buffered at pH 4.8, pH 7.0, or pH 9.0, or in 20 mM acetate or sorbate at pH 6.5. Genomes of 16 strains revealed over 100 mutations, including single-nucleotide polymorphisms (SNPs), large deletions, and insertion knockouts. Most strains acquired deletions in the benzoate-induced multiple antibiotic resistance (Mar) regulon or in associated regulators such as rob and cpxA, as well as the multidrug resistance (MDR) efflux pumps emrA, emrY, and mdtA Strains also lost or downregulated the Gad acid fitness regulon. In 5 mM benzoate or in 2 mM salicylate (2-hydroxybenzoate), most strains showed increased sensitivity to the antibiotics chloramphenicol and tetracycline; some strains were more sensitive than a marA knockout strain. Thus, our benzoate-evolved strains may reveal additional unknown drug resistance components. Benzoate or salicylate selection pressure may cause general loss of MDR genes and regulators. IMPORTANCE: Benzoate is a common food preservative, and salicylate is the primary active metabolite of aspirin. In the gut microbiome, genetic adaptation to salicylate may involve loss or downregulation of inducible multidrug resistance systems. This discovery implies that aspirin therapy may modulate the human gut microbiome to favor salicylate tolerance at the expense of drug resistance. Similar aspirin-associated loss of drug resistance might occur in bacterial pathogens found in arterial plaques.

Metagenomic Data Reveal Type I Polyketide Synthase Distributions Across Biomes.

Microbial polyketide synthase (PKS) genes encode the biosynthesis of many biomedically important natural products, yet only a small fraction of nature's polyketide biosynthetic potential has been realized. Much of this potential originates from type I PKSs (T1PKSs), which can be delineated into different classes and subclasses based on domain organization and structural features of the compounds encoded. Notably, phylogenetic relationships among PKS ketosynthase (KS) domains provide a method to classify the larger and more complex genes in which they occur. Increased access to large metagenomic datasets from diverse habitats provides opportunities to assess T1PKS biosynthetic diversity and distributions through the analysis of KS domain sequences. Here, we used the webtool NaPDoS2 to detect and classify over 35,000 type I KS domains from 137 metagenomic data sets reported from eight diverse biomes. We found biome-specific separation with soils enriched in modular cis -AT and hybrid cis -AT KSs relative to other biomes and marine sediments enriched in KSs associated with PUFA and enediyne biosynthesis. By extracting full-length KS domains, we linked the phylum Actinobacteria to soil-specific enediyne and cis -AT clades and identified enediyne and monomodular KSs in phyla from which the associated compound classes have not been reported. These sequences were phylogenetically distinct from those associated with experimentally characterized PKSs suggesting novel structures or enzyme functions remain to be discovered. Lastly, we employed our metagenome-extracted KS domains to evaluate commonly used type I KS PCR primers and identified modifications that could increase the KS sequence diversity recovered from amplicon libraries. Importance: Polyketides are a crucial source of medicines, agrichemicals, and other commercial products. Advances in our understanding of polyketide biosynthesis coupled with the accumulation of metagenomic sequence data provide new opportunities to assess polyketide biosynthetic potential across biomes. Here, we used the webtool NaPDoS2 to assess type I PKS diversity and distributions by detecting and classifying KS domains across 137 metagenomes. We show that biomes are differentially enriched in KS domain classes, providing a roadmap for future biodiscovery strategies. Further, KS phylogenies reveal both biome-specific clades that do not include biochemically characterized PKSs, highlighting the biosynthetic potential of poorly explored environments. The large metagenome-derived KS dataset allowed us to identify regions of commonly used type I KS PCR primers that could be modified to capture a larger extent of KS diversity. These results facilitate both the search for novel polyketides and our understanding of the biogeographical distribution of PKSs across earth's major biomes.

Metagenomic data reveals type I polyketide synthase distributions across biomes.

Microbial polyketide synthase (PKS) genes encode the biosynthesis of many biomedically or otherwise commercially important natural products. Despite extensive discovery efforts, metagenomic analyses suggest that only a small fraction of nature's polyketide biosynthetic potential has been realized. Much of this potential originates from type I PKSs (T1PKSs), which can be further delineated based on their domain organization and the structural features of the compounds they encode. Notably, phylogenetic relationships among ketosynthase (KS) domains provide an effective method to classify the larger and more complex T1PKS genes in which they occur. Increased access to large metagenomic data sets from diverse habitats provides opportunities to assess T1PKS biosynthetic diversity and distributions through their smaller and more tractable KS domain sequences. Here, we used the web tool NaPDoS2 to detect and classify over 35,000 type I KS domains from 137 metagenomic data sets reported from eight diverse, globally distributed biomes. We found biome-specific separation with soils enriched in KSs from modular cis-acetyltransferase (AT) and hybrid cis-AT KSs relative to other biomes and marine sediments enriched in KSs associated with polyunsaturated fatty acid and enediyne biosynthesis. We linked the phylum Actinobacteria to soil-derived enediyne and cis-AT KSs while marine-derived KSs associated with enediyne and monomodular PKSs were linked to phyla from which the compounds produced by these biosynthetic enzymes have not been reported. These KSs were phylogenetically distinct from those associated with experimentally characterized PKSs suggesting they may be associated with novel structures or enzyme functions. Finally, we employed our metagenome-extracted KS domains to evaluate the PCR primers commonly used to amplify type I KSs and identified modifications that could increase the KS sequence diversity recovered from amplicon libraries. IMPORTANCE Polyketides are a crucial source of medicines, agrichemicals, and other commercial products. Advances in our understanding of polyketide biosynthesis, coupled with the increased availability of metagenomic sequence data, provide new opportunities to assess polyketide biosynthetic potential across biomes. Here, we used the web tool NaPDoS2 to assess type I polyketide synthase (PKS) diversity and distributions by detecting and classifying ketosynthase (KS) domains across 137 metagenomes. We show that biomes are differentially enriched in type I KS domains, providing a roadmap for future biodiscovery strategies. Furthermore, KS phylogenies reveal biome-specific clades that do not include biochemically characterized PKSs, highlighting the biosynthetic potential of poorly explored environments. The large metagenome-derived KS data set allowed us to identify regions of commonly used type I KS PCR primers that could be modified to capture a larger extent of environmental KS diversity. These results facilitate both the search for novel polyketides and our understanding of the biogeographical distribution of PKSs across Earth's major biomes.

Expansion of Gamma-Butyrolactone Signaling Molecule Biosynthesis to Phosphotriester Natural Products.

Bacterial hormones, such as the iconic gamma-butyrolactone A-factor, are essential signaling molecules that regulate diverse physiological processes, including specialized metabolism. These low molecular weight compounds are common in Streptomyces species and display species-specific structural differences. Recently, unusual gamma-butyrolactone natural products called salinipostins were isolated from the marine actinomycete genus Salinispora based on their antimalarial properties. As the salinipostins possess a rare phosphotriester motif of unknown biosynthetic origin, we set out to explore its construction by the widely conserved 9-gene spt operon in Salinispora species. We show through a series of in vivo and in vitro studies that the spt gene cluster dually encodes the salinipostins and newly identified natural A-factor-like gamma-butyrolactones (Sal-GBLs). Remarkably, homologous biosynthetic gene clusters are widely distributed among many actinomycete genera, including Streptomyces, suggesting the significance of this operon in bacteria.