MIBiG 3.0: a community-driven effort to annotate experimentally validated biosynthetic gene clusters.

With an ever-increasing amount of (meta)genomic data being deposited in sequence databases, (meta)genome mining for natural product biosynthetic pathways occupies a critical role in the discovery of novel pharmaceutical drugs, crop protection agents and biomaterials. The genes that encode these pathways are often organised into biosynthetic gene clusters (BGCs). In 2015, we defined the Minimum Information about a Biosynthetic Gene cluster (MIBiG): a standardised data format that describes the minimally required information to uniquely characterise a BGC. We simultaneously constructed an accompanying online database of BGCs, which has since been widely used by the community as a reference dataset for BGCs and was expanded to 2021 entries in 2019 (MIBiG 2.0). Here, we describe MIBiG 3.0, a database update comprising large-scale validation and re-annotation of existing entries and 661 new entries. Particular attention was paid to the annotation of compound structures and biological activities, as well as protein domain selectivities. Together, these new features keep the database up-to-date, and will provide new opportunities for the scientific community to use its freely available data, e.g. for the training of new machine learning models to predict sequence-structure-function relationships for diverse natural products. MIBiG 3.0 is accessible online at https://mibig.secondarymetabolites.org/.

Microbiome composition modulates secondary metabolism in a multispecies bacterial community.

Bacterial secondary metabolites are a major source of antibiotics and other bioactive compounds. In microbial communities, these molecules can mediate interspecies interactions and responses to environmental change. Despite the importance of secondary metabolites in human health and microbial ecology, little is known about their roles and regulation in the context of multispecies communities. In a simplified model of the rhizosphere composed of Bacillus cereus, Flavobacterium johnsoniae, and Pseudomonas koreensis, we show that the dynamics of secondary metabolism depend on community species composition and interspecies interactions. Comparative metatranscriptomics and metametabolomics reveal that the abundance of transcripts of biosynthetic gene clusters (BGCs) and metabolomic molecular features differ between monocultures or dual cultures and a tripartite community. In both two- and three-member cocultures, P. koreensis modified expression of BGCs for zwittermicin, petrobactin, and other secondary metabolites in B. cereus and F. johnsoniae, whereas the BGC transcriptional response to the community in P. koreensis itself was minimal. Pairwise and tripartite cocultures with P. koreensis displayed unique molecular features that appear to be derivatives of lokisin, suggesting metabolic handoffs between species. Deleting the BGC for koreenceine, another P. koreensis metabolite, altered transcript and metabolite profiles across the community, including substantial up-regulation of the petrobactin and bacillibactin BGCs in B. cereus, suggesting that koreenceine represses siderophore production. Results from this model community show that bacterial BGC expression and chemical output depend on the identity and biosynthetic capacity of coculture partners, suggesting community composition and microbiome interactions may shape the regulation of secondary metabolism in nature.

Genomic diversification of the specialized parasite of the fungus-growing ant symbiosis.

Fungi shape the diversity of life. Characterizing the evolution of fungi is critical to understanding symbiotic associations across kingdoms. In this study, we investigate the genomic and metabolomic diversity of the genus Escovopsis, a specialized parasite of fungus-growing ant gardens. Based on 25 high-quality draft genomes, we show that Escovopsis forms a monophyletic group arising from a mycoparasitic fungal ancestor 61.82 million years ago (Mya). Across the evolutionary history of fungus-growing ants, the dates of origin of most clades of Escovopsis correspond to the dates of origin of the fungus-growing ants whose gardens they parasitize. We reveal that genome reduction, determined by both genomic sequencing and flow cytometry, is a consistent feature across the genus Escovopsis, largely occurring in coding regions, specifically in the form of gene loss and reductions in copy numbers of genes. All functional gene categories have reduced copy numbers, but resistance and virulence genes maintain functional diversity. Biosynthetic gene clusters (BGCs) contribute to phylogenetic differences among Escovopsis spp., and sister taxa in the Hypocreaceae. The phylogenetic patterns of co-diversification among BGCs are similarly exhibited across mass spectrometry analyses of the metabolomes of Escovopsis and their sister taxa. Taken together, our results indicate that Escovopsis spp. evolved unique genomic repertoires to specialize on the fungus-growing ant-microbe symbiosis.

Metabolomics and Genomics Enable the Discovery of a New Class of Nonribosomal Peptidic Metallophores from a Marine Micromonospora.

Although microbial genomes harbor an abundance of biosynthetic gene clusters, there remain substantial technological gaps that impair the direct correlation of newly discovered gene clusters and their corresponding secondary metabolite products. As an example of one approach designed to minimize or bridge such gaps, we employed hierarchical clustering analysis and principal component analysis (hcapca, whose sole input is MS data) to prioritize 109 marine Micromonospora strains and ultimately identify novel strain WMMB482 as a candidate for in-depth "metabologenomics" analysis following its prioritization. Highlighting the power of current MS-based technologies, not only did hcapca enable the discovery of one new, nonribosomal peptide bearing an incredible diversity of unique functional groups, but metabolomics for WMMB482 unveiled 16 additional congeners via the application of Global Natural Product Social molecular networking (GNPS), herein named ecteinamines A-Q (1-17). The ecteinamines possess an unprecedented skeleton housing a host of uncommon functionalities including a menaquinone pathway-derived 2-naphthoate moiety, 4-methyloxazoline, the first example of a naturally occurring Ψ[CH2NH] "reduced amide", a methylsulfinyl moiety, and a d-cysteinyl residue that appears to derive from a unique noncanonical epimerase domain. Extensive in silico analysis of the ecteinamine (ect) biosynthetic gene cluster and stable isotope-feeding experiments helped illuminate the novel enzymology driving ecteinamine assembly as well the role of cluster collaborations or "duets" in producing such structurally complex agents. Finally, ecteinamines were found to bind nickel, cobalt, zinc, and copper, suggesting a possible biological role as broad-spectrum metallophores.

Nutrient Availability Shifts the Biosynthetic Potential of Soil-Derived Microbial Communities.

Secondary metabolites produced by microorganisms are the main source of antimicrobials and other pharmaceutical drugs. Soil microbes have been the primary discovery source for these secondary metabolites, often producing complex organic compounds with specific biological activities. Research suggests that secondary metabolism broadly shapes microbial ecological interactions, but little is known about the factors that shape the abundance, distribution, and diversity of biosynthetic gene clusters in the context of microbial communities. In this study, we investigate the role of nutrient availability on the abundance of biosynthetic gene clusters in soil-derived microbial consortia. Soil microbial consortia enriched in high sugar medium (150 mg/L of glucose and 200 mg/L of trehalose) had more biosynthetic gene clusters and higher inhibitory activity than those enriched in low sugar medium (15 mg/L of glucose + 20 mg/L of trehalose). Our results demonstrate that experimental microbial communities are a promising tool to study the ecology of specialized metabolites.

Bacillimidazoles A-F, Imidazolium-Containing Compounds Isolated from a Marine Bacillus.

Chemical investigations of a marine sponge-associated Bacillus revealed six new imidazolium-containing compounds, bacillimidazoles A-F (1-6). Previous reports of related imidazolium-containing natural products are rare. Initially unveiled by timsTOF (trapped ion mobility spectrometry) MS data, extensive HRMS and 1D and 2D NMR analyses enabled the structural elucidation of 1-6. In addition, a plausible biosynthetic pathway to bacillimidazoles is proposed based on isotopic labeling experiments and invokes the highly reactive glycolytic adduct 2,3-butanedione. Combined, the results of structure elucidation efforts, isotopic labeling studies and bioinformatics suggest that 1-6 result from a fascinating intersection of primary and secondary metabolic pathways in Bacillus sp. WMMC1349. Antimicrobial assays revealed that, of 1-6, only compound six displayed discernible antibacterial activity, despite the close structural similarities shared by all six natural products.

Needles in haystacks: reevaluating old paradigms for the discovery of bacterial secondary metabolites.

Covering: up to 2021Natural products research is in the midst of a renaissance ushered in by a modern understanding of microbiology and the technological explosions of genomics and metabolomics. As the exploration of uncharted chemical space expands into high-throughput discovery campaigns, it has become increasingly clear how design elements influence success: (bio)geography, habitat, community dynamics, culturing/induction methods, screening methods, dereplication, and more. We explore critical considerations and assumptions in natural products discovery. We revisit previous estimates of chemical rediscovery and discuss their relatedness to study design and producer taxonomy. Through frequency analyses of biosynthetic gene clusters in publicly available genomic data, we highlight phylogenetic biases that influence rediscovery rates. Through selected examples of how study design at each level determines discovery outcomes, we discuss the challenges and opportunities for the future of high-throughput natural product discovery.

The confluence of big data and evolutionary genome mining for the discovery of natural products.

This review covers literature between 2003-2021The development and application of genome mining tools has given rise to ever-growing genetic and chemical databases and propelled natural products research into the modern age of Big Data. Likewise, an explosion of evolutionary studies has unveiled genetic patterns of natural products biosynthesis and function that support Darwin's theory of natural selection and other theories of adaptation and diversification. In this review, we aim to highlight how Big Data and evolutionary thinking converge in the study of natural products, and how this has led to an emerging sub-discipline of evolutionary genome mining of natural products. First, we outline general principles to best utilize Big Data in natural products research, addressing key considerations needed to provide evolutionary context. We then highlight successful examples where Big Data and evolutionary analyses have been combined to provide bioinformatic resources and tools for the discovery of novel natural products and their biosynthetic enzymes. Rather than an exhaustive list of evolution-driven discoveries, we highlight examples where Big Data and evolutionary thinking have been embraced for the evolutionary genome mining of natural products. After reviewing the nascent history of this sub-discipline, we discuss the challenges and opportunities of genomic and metabolomic tools with evolutionary foundations and/or implications and provide a future outlook for this emerging and exciting field of natural product research.

Evolution and Ecology of Actinobacteria and Their Bioenergy Applications.

The ancient phylum Actinobacteria is composed of phylogenetically and physiologically diverse bacteria that help Earth's ecosystems function. As free-living organisms and symbionts of herbivorous animals, Actinobacteria contribute to the global carbon cycle through the breakdown of plant biomass. In addition, they mediate community dynamics as producers of small molecules with diverse biological activities. Together, the evolution of high cellulolytic ability and diverse chemistry, shaped by their ecological roles in nature, make Actinobacteria a promising group for the bioenergy industry. Specifically, their enzymes can contribute to industrial-scale breakdown of cellulosic plant biomass into simple sugars that can then be converted into biofuels. Furthermore, harnessing their ability to biosynthesize a range of small molecules has potential for the production of specialty biofuels.

Antileishmanial macrolides from ant-associated Streptomyces sp. ISID311.

Three antifungal macrolides cyphomycin (1), caniferolide C (2) and GT-35 (3) were isolated from Streptomyces sp. ISID311, a bacterial symbiont associated with Cyphomyrmex fungus-growing ants. The planar structures of these compounds were established by 1 and 2D NMR data and MS analysis. The relative configurations of 1-3 were established using Kishi's universal NMR database method, NOE/ROE analysis and coupling constants analysis assisted by comparisons with NMR data of related compounds. Detailed bioinformatic analysis of cyphomycin biosynthetic gene cluster confirmed the stereochemical assignments. Compounds 1-3 displayed high antagonism against different strains of Escovopsis sp., pathogen fungi specialized to the fungus-growing ant system. Compounds 1-3 also exhibited potent antiprotozoal activity against intracellular amastigotes of the human parasite Leishmania donovani with IC50 values of 2.32, 0.091 and 0.073 µM, respectively, with high selectivity indexes.

Bacillibactins E and F from a Marine Sponge-Associated Bacillus sp.

Chemical investigation of a marine sponge-associated Bacillus sp. led to the discovery of bacillibactins E and F (1 and 2). Despite containing the well-established cyclic triester core of iron-binding natural products such as enterobactin, bacillibactins E and F (1 and 2) are the first bacterial siderophores that contain nicotinic and benzoic acid moieties. The structures of the new compounds, including their absolute configurations, were determined by extensive spectroscopic analyses and Marfey's method. A plausible biosynthetic pathway to 1 and 2 is proposed; this route bears great similarity to other previously established bacillibactin-like pathways but appears to differentiate itself by a promiscuous DhbE, which likely installs the nicotinic moiety of 1 and the benzoic acid group of 2.

Evolutionary dynamics of natural product biosynthesis in bacteria.

Covering: 2008 up to 2019The forces of biochemical adaptive evolution operate at the level of genes, manifesting in complex phenotypes and the global biodiversity of proteins and metabolites. While evolutionary histories have been deciphered for some other complex traits, the origins of natural product biosynthesis largely remain a mystery. This fundamental knowledge gap is surprising given the many decades of research probing the genetic, chemical, and biophysical mechanisms of bacterial natural product biosynthesis. Recently, evolutionary thinking has begun to permeate this otherwise mechanistically dominated field. Natural products are now sometimes referred to as 'specialized' rather than 'secondary' metabolites, reinforcing the importance of their biological and ecological functions. Here, we review known evolutionary mechanisms underlying the overwhelming chemical diversity of bacterial secondary metabolism, focusing on enzyme promiscuity and the evolution of enzymatic domains that enable metabolic traits. We discuss the mechanisms that drive the assembly of natural product biosynthetic gene clusters and propose formal definitions for 'specialized' and 'secondary' metabolism. We further explore how biosynthetic gene clusters evolve to synthesize related molecular species, and in turn how the biological and ecological roles that emerge from metabolic diversity are acted on by selection. Finally, we reconcile chemical, functional, and genetic data into an evolutionary model, the dynamic chemical matrix evolutionary hypothesis, in which the relationships between chemical distance, biomolecular activity, and relative fitness shape adaptive landscapes.

Local Adaptation of Bacterial Symbionts within a Geographic Mosaic of Antibiotic Coevolution.

The geographic mosaic theory of coevolution (GMC) posits that coevolutionary dynamics go beyond local coevolution and are comprised of the following three components: geographic selection mosaics, coevolutionary hot spots, and trait remixing. It is unclear whether the GMC applies to bacteria, as horizontal gene transfer and cosmopolitan dispersal may violate theoretical assumptions. Here, we test key GMC predictions in an antibiotic-producing bacterial symbiont (genus Pseudonocardia) that protects the crops of neotropical fungus-farming ants (Apterostigma dentigerum) from a specialized pathogen (genus Escovopsis). We found that Pseudonocardia antibiotic inhibition of common Escovopsis pathogens was elevated in A. dentigerum colonies from Panama compared to those from Costa Rica. Furthermore, a Panama Canal Zone population of Pseudonocardia on Barro Colorado Island (BCI) was locally adapted, whereas two neighboring populations were not, consistent with a GMC-predicted selection mosaic and a hot spot of adaptation surrounded by areas of maladaptation. Maladaptation was shaped by incongruent Pseudonocardia-Escovopsis population genetic structure, whereas local adaptation was facilitated by geographic isolation on BCI after the flooding of the Panama Canal. Genomic assessments of antibiotic potential of 29 Pseudonocardia strains identified diverse and unique biosynthetic gene clusters in BCI strains despite low genetic diversity in the core genome. The strength of antibiotic inhibition was not correlated with the presence/absence of individual biosynthetic gene clusters or with parasite location. Rather, biosynthetic gene clusters have undergone selective sweeps, suggesting that the trait remixing dynamics conferring the long-term maintenance of antibiotic potency rely on evolutionary genetic changes within already-present biosynthetic gene clusters and not simply on the horizontal acquisition of novel genetic elements or pathways.IMPORTANCE Recently, coevolutionary theory in macroorganisms has been advanced by the geographic mosaic theory of coevolution (GMC), which considers how geography and local adaptation shape coevolutionary dynamics. Here, we test GMC in an ancient symbiosis in which the ant Apterostigma dentigerum cultivates fungi in an agricultural system analogous to human farming. The cultivars are parasitized by the fungus Escovopsis The ants maintain symbiotic actinobacteria with antibiotic properties that help combat Escovopsis infection. This antibiotic symbiosis has persisted for tens of millions of years, raising the question of how antibiotic potency is maintained over these time scales. Our study tests the GMC in a bacterial defensive symbiosis and in a multipartite symbiosis framework. Our results show that this multipartite symbiotic system conforms to the GMC and demonstrate that this theory is applicable in both microbes and indirect symbiont-symbiont interactions.

Competition among Nasal Bacteria Suggests a Role for Siderophore-Mediated Interactions in Shaping the Human Nasal Microbiota.

Resources available in the human nasal cavity are limited. Therefore, to successfully colonize the nasal cavity, bacteria must compete for scarce nutrients. Competition may occur directly through interference (e.g., antibiotics) or indirectly by nutrient sequestration. To investigate the nature of nasal bacterial competition, we performed coculture inhibition assays between nasal Actinobacteria and Staphylococcus spp. We found that isolates of coagulase-negative staphylococci (CoNS) were sensitive to growth inhibition by Actinobacteria but that Staphylococcus aureus isolates were resistant to inhibition. Among Actinobacteria, we observed that Corynebacterium spp. were variable in their ability to inhibit CoNS. We sequenced the genomes of 10 Corynebacterium species isolates, including 3 Corynebacterium propinquum isolates that strongly inhibited CoNS and 7 other Corynebacterium species isolates that only weakly inhibited CoNS. Using a comparative genomics approach, we found that the C. propinquum genomes were enriched in genes for iron acquisition and harbored a biosynthetic gene cluster (BGC) for siderophore production, absent in the noninhibitory Corynebacterium species genomes. Using a chrome azurol S assay, we confirmed that C. propinquum produced siderophores. We demonstrated that iron supplementation rescued CoNS from inhibition by C. propinquum, suggesting that inhibition was due to iron restriction through siderophore production. Through comparative metabolomics and molecular networking, we identified the siderophore produced by C. propinquum as dehydroxynocardamine. Finally, we confirmed that the dehydroxynocardamine BGC is expressed in vivo by analyzing human nasal metatranscriptomes from the NIH Human Microbiome Project. Together, our results suggest that bacteria produce siderophores to compete for limited available iron in the nasal cavity and improve their fitness.IMPORTANCE Within the nasal cavity, interference competition through antimicrobial production is prevalent. For instance, nasal Staphylococcus species strains can inhibit the growth of other bacteria through the production of nonribosomal peptides and ribosomally synthesized and posttranslationally modified peptides. In contrast, bacteria engaging in exploitation competition modify the external environment to prevent competitors from growing, usually by hindering access to or depleting essential nutrients. As the nasal cavity is a nutrient-limited environment, we hypothesized that exploitation competition occurs in this system. We determined that Corynebacterium propinquum produces an iron-chelating siderophore, and this iron-sequestering molecule correlates with the ability to inhibit the growth of coagulase-negative staphylococci. Furthermore, we found that the genes required for siderophore production are expressed in vivo Thus, although siderophore production by bacteria is often considered a virulence trait, our work indicates that bacteria may produce siderophores to compete for limited iron in the human nasal cavity.

Complexity of Complement Resistance Factors Expressed by Acinetobacter baumannii Needed for Survival in Human Serum.

Acinetobacter baumannii is a bacterial pathogen with increasing impact in healthcare settings, due in part to this organism's resistance to many antimicrobial agents, with pneumonia and bacteremia as the most common manifestations of disease. A significant proportion of clinically relevant A. baumannii strains are resistant to killing by normal human serum (NHS), an observation supported in this study by showing that 12 out of 15 genetically diverse strains of A. baumannii are resistant to NHS killing. To expand our understanding of the genetic basis of A. baumannii serum resistance, a transposon (Tn) sequencing (Tn-seq) approach was used to identify genes contributing to this trait. An ordered Tn library in strain AB5075 with insertions in every nonessential gene was subjected to selection in NHS. We identified 50 genes essential for the survival of A. baumannii in NHS, including already known serum resistance factors, and many novel genes not previously associated with serum resistance. This latter group included the maintenance of lipid asymmetry genetic pathway as a key determinant in protecting A. baumannii from the bactericidal activity of NHS via the alternative complement pathway. Follow-up studies validated the role of eight additional genes identified by Tn-seq in A. baumannii resistance to killing by NHS but not by normal mouse serum, highlighting the human species specificity of A. baumannii serum resistance. The identification of a large number of genes essential for serum resistance in A. baumannii indicates the degree of complexity needed for this phenotype, which might reflect a general pattern that pathogens rely on to cause serious infections.

antiSMASH 4.0-improvements in chemistry prediction and gene cluster boundary identification.

Many antibiotics, chemotherapeutics, crop protection agents and food preservatives originate from molecules produced by bacteria, fungi or plants. In recent years, genome mining methodologies have been widely adopted to identify and characterize the biosynthetic gene clusters encoding the production of such compounds. Since 2011, the 'antibiotics and secondary metabolite analysis shell-antiSMASH' has assisted researchers in efficiently performing this, both as a web server and a standalone tool. Here, we present the thoroughly updated antiSMASH version 4, which adds several novel features, including prediction of gene cluster boundaries using the ClusterFinder method or the newly integrated CASSIS algorithm, improved substrate specificity prediction for non-ribosomal peptide synthetase adenylation domains based on the new SANDPUMA algorithm, improved predictions for terpene and ribosomally synthesized and post-translationally modified peptides cluster products, reporting of sequence similarity to proteins encoded in experimentally characterized gene clusters on a per-protein basis and a domain-level alignment tool for comparative analysis of trans-AT polyketide synthase assembly line architectures. Additionally, several usability features have been updated and improved. Together, these improvements make antiSMASH up-to-date with the latest developments in natural product research and will further facilitate computational genome mining for the discovery of novel bioactive molecules.

Emerging evolutionary paradigms in antibiotic discovery.

Antibiotics revolutionized medicine and remain its cornerstone. Despite their global importance and the continuous threat of resistant pathogens, few antibiotics have been discovered in recent years. Natural products, especially the secondary metabolites of Actinobacteria, have been the traditional discovery source of antibiotics. In nature, the chemistry of antibiotic natural products is shaped by the unique evolution and ecology of their producing organisms, yet these influences remain largely unknown. Here, we highlight the ecology of antibiotics employed by microbes in defensive symbioses and review the evolutionary processes underlying the chemical diversity and activity of microbe-derived antibiotics, including the dynamics of vertical and lateral transmission of biosynthetic pathways and the evolution of efficacy, targeting specificity, and toxicity. We argue that a deeper understanding of the ecology and evolution of microbial interactions and the metabolites that mediate them will allow for an alternative, rational approach to discover new antibiotics.

SANDPUMA: ensemble predictions of nonribosomal peptide chemistry reveal biosynthetic diversity across Actinobacteria.

SUMMARY: Nonribosomally synthesized peptides (NRPs) are natural products with widespread applications in medicine and biotechnology. Many algorithms have been developed to predict the substrate specificities of nonribosomal peptide synthetase adenylation (A) domains from DNA sequences, which enables prioritization and dereplication, and integration with other data types in discovery efforts. However, insufficient training data and a lack of clarity regarding prediction quality have impeded optimal use. Here, we introduce prediCAT, a new phylogenetics-inspired algorithm, which quantitatively estimates the degree of predictability of each A-domain. We then systematically benchmarked all algorithms on a newly gathered, independent test set of 434 A-domain sequences, showing that active-site-motif-based algorithms outperform whole-domain-based methods. Subsequently, we developed SANDPUMA, a powerful ensemble algorithm, based on newly trained versions of all high-performing algorithms, which significantly outperforms individual methods. Finally, we deployed SANDPUMA in a systematic investigation of 7635 Actinobacteria genomes, suggesting that NRP chemical diversity is much higher than previously estimated. SANDPUMA has been integrated into the widely used antiSMASH biosynthetic gene cluster analysis pipeline and is also available as an open-source, standalone tool. AVAILABILITY AND IMPLEMENTATION: SANDPUMA is freely available at https://bitbucket.org/chevrm/sandpuma and as a docker image at https://hub.docker.com/r/chevrm/sandpuma/ under the GNU Public License 3 (GPL3). CONTACT: chevrette@wisc.edu or marnix.medema@wur.nl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.