Significance Microbial natural products are a major source of new drug leads, yet discovery efforts are constrained by the lack of information describing the diversity and distributions of the associated biosynthetic pathways among bacteria. Using the marine actinomycete genus Salinispora as a model, we analyzed genome sequence data from 75 closely related strains. The results provide evidence for high levels of pathway diversity, with most being acquired relatively recently in the evolution of the genus. The distributions and evolutionary histories of these pathways provide insight into the mechanisms that generate new chemical diversity and the strategies used by bacteria to maximize their population-level capacity to produce diverse secondary metabolites.

Recent genome sequencing efforts have led to the rapid accumulation of uncharacterized or “orphaned” secondary metabolic biosynthesis gene clusters (BGCs) in public databases. This increase in DNA-sequenced big data has given rise to significant challenges in the applied field of natural product genome mining, including (i) how to prioritize the characterization of orphan BGCs, and (ii) how to rapidly connect genes to biosynthesized small molecules. Here we show that by correlating putative antibiotic resistance genes that encode target-modified proteins with orphan BGCs, we predict the biological function of pathway specific small molecules before they have been revealed in a process we call target-directed genome mining. By querying the pan-genome of 86 Salinispora bacterial genomes for duplicated house-keeping genes co-localized with natural product BGCs, we prioritized an orphan polyketide synthase-nonribosomal peptide synthetase hybrid BGC (tlm) with a putative fatty acid synthase resistance gene. We employed a new synthetic double-stranded DNA-mediated cloning strategy based on transformation-associated recombination to efficiently capture tlm and the related ttm BGCs directly from genomic DNA and to heterologously express them in Streptomyces hosts. We show the production of a group of unusual thiotetronic acid natural products, including the well-known fatty acid synthase inhibitor thiolactomycin that was first described over 30 years ago, yet never at the genetic level in regards to biosynthesis and auto-resistance. This finding not only validates the target-directed genome mining strategy for the discovery of antibiotic producing gene clusters without a priori knowledge of the molecule synthesized, but also paves the way for the investigation of novel enzymology involved in thiotetronic acid natural product biosynthesis.

Summary Comparative genomics is providing new opportunities to address the diversity and distributions of genes encoding the biosynthesis of specialized metabolites. An analysis of 119 genome sequences representing three closely related species of the marine actinomycete genus Salinispora reveals extraordinary biosynthetic diversity in the form of 176 distinct biosynthetic gene clusters (BGCs) of which only 24 have been linked to their products. Remarkably, more than half of the BGCs were observed in only one or two strains, suggesting they were acquired relatively recently in the evolutionary history of the genus. These acquired gene clusters are concentrated in specific genomic islands, which represent hot spots for BGC acquisition. While most BGCs are stable in terms of their chromosomal position, others migrated to different locations or were exchanged with unrelated gene clusters suggesting a plug and play type model of evolution that provides a mechanism to test the relative fitness effects of specialized metabolites. Transcriptome analyses were used to address the relationships between BGC abundance, chromosomal position, and product discovery. The results indicate that recently acquired BGCs can be functional and that complex evolutionary processes shape the micro-diversity of specialized metabolism observed in closely related environmental bacteria.