A Natural Product Chemist's Guide to Unlocking Silent Biosynthetic Gene Clusters.

Microbial natural products have provided an important source of therapeutic leads and motivated research and innovation in diverse scientific disciplines. In recent years, it has become evident that bacteria harbor a large, hidden reservoir of potential natural products in the form of silent or cryptic biosynthetic gene clusters (BGCs). These can be readily identified in microbial genome sequences but do not give rise to detectable levels of a natural product. Herein, we provide a useful organizational framework for the various methods that have been implemented for interrogating silent BGCs. We divide all available approaches into four categories. The first three are endogenous strategies that utilize the native host in conjunction with classical genetics, chemical genetics, or different culture modalities. The last category comprises expression of the entire BGC in a heterologous host. For each category, we describe the rationale, recent applications, and associated advantages and limitations.

A Pathogen-Responsive Gene Cluster for Highly Modified Fatty Acids in Tomato.

In response to biotic stress, plants produce suites of highly modified fatty acids that bear unusual chemical functionalities. Despite their chemical complexity and proposed roles in pathogen defense, little is known about the biosynthesis of decorated fatty acids in plants. Falcarindiol is a prototypical acetylenic lipid present in carrot, tomato, and celery that inhibits growth of fungi and human cancer cell lines. Using a combination of untargeted metabolomics and RNA sequencing, we discovered a biosynthetic gene cluster in tomato (Solanum lycopersicum) required for falcarindiol production. By reconstituting initial biosynthetic steps in a heterologous host and generating transgenic pathway mutants in tomato, we demonstrate a direct role of the cluster in falcarindiol biosynthesis and resistance to fungal and bacterial pathogens in tomato leaves. This work reveals a mechanism by which plants sculpt their lipid pool in response to pathogens and provides critical insight into the complex biochemistry of alkynyl lipid production.

Discovery of Reactive Microbiota-Derived Metabolites that Inhibit Host Proteases.

The gut microbiota modulate host biology in numerous ways, but little is known about the molecular mediators of these interactions. Previously, we found a widely distributed family of nonribosomal peptide synthetase gene clusters in gut bacteria. Here, by expressing a subset of these clusters in Escherichia coli or Bacillus subtilis, we show that they encode pyrazinones and dihydropyrazinones. At least one of the 47 clusters is present in 88% of the National Institutes of Health Human Microbiome Project (NIH HMP) stool samples, and they are transcribed under conditions of host colonization. We present evidence that the active form of these molecules is the initially released peptide aldehyde, which bears potent protease inhibitory activity and selectively targets a subset of cathepsins in human cell proteomes. Our findings show that an approach combining bioinformatics, synthetic biology, and heterologous gene cluster expression can rapidly expand our knowledge of the metabolic potential of the microbiota while avoiding the challenges of cultivating fastidious commensals.

Neofunctionalization of an OMT cluster dominates polymethoxyflavone biosynthesis associated with the domestication of citrus.

Polymethoxyflavones (PMFs) are a class of abundant specialized metabolites with remarkable anticancer properties in citrus. Multiple methoxy groups in PMFs are derived from methylation modification catalyzed by a series of hydroxylases and O-methyltransferases (OMTs). However, the specific OMTs that catalyze the systematic O-methylation of hydroxyflavones remain largely unknown. Here, we report that PMFs are highly accumulated in wild mandarins and mandarin-derived accessions, while undetectable in early-diverging citrus species and related species. Our results demonstrated that three homologous genes, CreOMT3, CreOMT4, and CreOMT5, are crucial for PMF biosynthesis in citrus, and their encoded methyltransferases exhibit multisite O-methylation activities for hydroxyflavones, producing seven PMFs in vitro and in vivo. Comparative genomic and syntenic analyses indicated that the tandem CreOMT3, CreOMT4, and CreOMT5 may be duplicated from CreOMT6 and contributes to the genetic basis of PMF biosynthesis in the mandarin group through neofunctionalization. We also demonstrated that N17 in CreOMT4 is an essential amino acid residue for C3-, C5-, C6-, and C3'-O-methylation activity and provided a rationale for the functional deficiency of OMT6 to produce PMFs in early-diverging citrus and some domesticated citrus species. A 1,041-bp deletion in the CreOMT4 promoter, which is found in most modern cultivated mandarins, has reduced the PMF content relative to that in wild and early-admixture mandarins. This study provides a framework for reconstructing PMF biosynthetic pathways, which may facilitate the breeding of citrus fruits with enhanced health benefits.

Evolution of the insect Hox gene cluster: Comparative analysis across 243 species.

The Hox gene cluster is an iconic example of evolutionary conservation between divergent animal lineages, providing evidence for ancient similarities in the genetic control of embryonic development. However, there are differences between taxa in gene order, gene number and genomic organisation implying conservation is not absolute. There are also examples of radical functional change of Hox genes; for example, the ftz, zen and bcd genes in insects play roles in segmentation, extraembryonic membrane formation and body polarity, rather than specification of anteroposterior position. There have been detailed descriptions of Hox genes and Hox gene clusters in several insect species, including important model systems, but a large-scale overview has been lacking. Here we extend these studies using the publicly-available complete genome sequences of 243 insect species from 13 orders. We show that the insect Hox cluster is characterised by large intergenic distances, consistently extreme in Odonata, Orthoptera, Hemiptera and Trichoptera, and always larger between the 'posterior' Hox genes. We find duplications of ftz and zen in many species and multiple independent cluster breaks, although certain modules of neighbouring genes are rarely broken apart suggesting some organisational constraints. As more high-quality genomes are obtained, a challenge will be to relate structural genomic changes to phenotypic change across insect phylogeny.

Multi-omic analysis tools for microbial metabolites prediction.

How to resolve the metabolic dark matter of microorganisms has long been a challenging problem in discovering active molecules. Diverse omics tools have been developed to guide the discovery and characterization of various microbial metabolites, which make it gradually possible to predict the overall metabolites for individual strains. The combinations of multi-omic analysis tools effectively compensates for the shortcomings of current studies that focus only on single omics or a broad class of metabolites. In this review, we systematically update, categorize and sort out different analysis tools for microbial metabolites prediction in the last five years to appeal for the multi-omic combination on the understanding of the metabolic nature of microbes. First, we provide the general survey on different updated prediction databases, webservers, or software that based on genomics, transcriptomics, proteomics, and metabolomics, respectively. Then, we discuss the essentiality on the integration of multi-omics data to predict metabolites of different microbial strains and communities, as well as stressing the combination of other techniques, such as systems biology methods and data-driven algorithms. Finally, we identify key challenges and trends in developing multi-omic analysis tools for more comprehensive prediction on diverse microbial metabolites that contribute to human health and disease treatment.

A bacterial-like Pictet-Spenglerase drives the evolution of fungi to produce ß-carboline glycosides together with separate genes.

Diverse ß-carboline (ßC) alkaloids are produced by microbes, plants, and animals with myriad bioactivities and drug potentials. However, the biosynthetic mechanism of ßCs remains largely elusive, especially regarding the hydroxyl and glucosyl modifications of ßCs. Here, we report the presence of the bacterial-like Pictet-Spenglerase gene Fcs1 in the entomopathogenic Beauveria fungi that can catalyze the biosynthesis of the ßC skeleton. The overexpression of Fcs1 in Beauveria bassiana led to the identification of six ßC methyl glycosides, termed bassicarbosides (BCSs) A-F. We verified that the cytochrome P450 (CYP) genes adjacent to Fcs1 cannot oxidize ßCs. Alternatively, the separated CYP684B2 family gene Fcs2 was identified to catalyze ßC hydroxylation together with its cofactor gene Fcs3. The functional homologue of Fcs2 is only present in the Fcs1-containing fungi and highly similar to the Fcs1-connected yet nonfunctional CYP. Both evolved quicker than those from fungi without Fcs1 homologues. Finally, the paired methyl/glucosyl transferase genes were verified to mediate the production of BCSs from hydroxy-ßCs. All these functionally verified genes are located on different chromosomes of Beauveria, which is in contrast to the typical content-clustered feature of fungal biosynthetic gene clusters (BGCs). We also found that the production of BCSs selectively contributed to fungal infection of different insect species. Our findings shed light on the biosynthetic mechanism of ßC glycosides, including the identification of a ßC hydroxylase. The results of this study also propose an evolving process of fungal BGC formation following the horizontal transfer of a bacterial gene to fungi.

Alternative splicing of lncRNA LAIR fine-tunes the regulation of neighboring yield-related gene LRK1 expression in rice.

The diversity in alternative splicing of long noncoding RNAs (lncRNAs) poses a challenge for functional annotation of lncRNAs. Moreover, little is known on the effects of alternatively spliced lncRNAs on crop yield. In this study, we cloned nine isoforms resulting from the alternative splicing of the lncRNA LAIR in rice. The LAIR isoforms are generated via alternative 5'/3' splice sites and different combinations of specific introns. All LAIR isoforms activate the expression of the neighboring LRK1 gene and enhance yield-related rice traits. In addition, there are slight differences in the binding ability of LAIR isoforms to the epigenetic modification-related proteins OsMOF and OsWDR5, which affect the enrichment of H4K16ac and H3K4me3 at the LRK1 locus, and consequently fine-tune the regulation of LRK1 expression and yield-related traits. These differences in binding may be caused by polymorphic changes to the RNA secondary structure resulting from alternative splicing. It was also observed that the composition of LAIR isoforms was sensitive to abiotic stress. These findings suggest that the alternative splicing of LAIR leads to the formation of a functional transcript population that precisely regulates yield-related gene expression, which may be relevant for phenotypic polymorphism-based crop breeding under changing environmental conditions.

The evolution of the GALactose utilization pathway in budding yeasts.

The Leloir galactose utilization or GAL pathway of budding yeasts, including that of the baker's yeast Saccharomyces cerevisiae and the opportunistic human pathogen Candida albicans, breaks down the sugar galactose for energy and biomass production. The GAL pathway has long served as a model system for understanding how eukaryotic metabolic pathways, including their modes of regulation, evolve. More recently, the physical linkage of the structural genes GAL1, GAL7, and GAL10 in diverse budding yeast genomes has been used as a model for understanding the evolution of gene clustering. In this review, we summarize exciting recent work on three different aspects of this iconic pathway's evolution: gene cluster organization, GAL gene regulation, and the population genetics of the GAL pathway.

Genes and evolutionary fates of the amanitin biosynthesis pathway in poisonous mushrooms.

The deadly toxin α-amanitin is a bicyclic octapeptide biosynthesized on ribosomes. A phylogenetically disjunct group of mushrooms in Agaricales (Amanita, Lepiota, and Galerina) synthesizes α-amanitin. This distribution of the toxin biosynthetic pathway is possibly related to the horizontal transfer of metabolic gene clusters among taxonomically unrelated mushrooms with overlapping habitats. Here, our work confirms that two biosynthetic genes, P450-29 and FMO1, are oxygenases important for amanitin biosynthesis. Phylogenetic and genetic analyses of these genes strongly support their origin through horizontal transfer, as is the case for the previously characterized biosynthetic genes MSDIN and POPB. Our analysis of multiple genomes showed that the evolution of the α-amanitin biosynthetic pathways in the poisonous agarics in the Amanita, Lepiota, and Galerina clades entailed distinct evolutionary pathways including gene family expansion, biosynthetic genes, and genomic rearrangements. Unrelated poisonous fungi produce the same deadly amanitin toxins using variations of the same pathway. Furthermore, the evolution of the amanitin biosynthetic pathway(s) in Amanita species generates a much wider range of toxic cyclic peptides. The results reported here expand our understanding of the genetics, diversity, and evolution of the toxin biosynthetic pathway in fungi.

Domoic acid biosynthesis in the red alga Chondria armata suggests a complex evolutionary history for toxin production.

Domoic acid (DA), the causative agent of amnesic shellfish poisoning, is produced by select organisms within two distantly related algal clades: planktonic diatoms and red macroalgae. The biosynthetic pathway to isodomoic acid A was recently solved in the harmful algal bloom-forming diatom Pseudonitzschia multiseries, establishing the genetic basis for the global production of this potent neurotoxin. Herein, we sequenced the 507-Mb genome of Chondria armata, the red macroalgal seaweed from which DA was first isolated in the 1950s, identifying several copies of the red algal DA (rad) biosynthetic gene cluster. The rad genes are organized similarly to the diatom DA biosynthesis cluster in terms of gene synteny, including a cytochrome P450 (CYP450) enzyme critical to DA production that is notably absent in red algae that produce the simpler kainoid neurochemical, kainic acid. The biochemical characterization of the N-prenyltransferase (RadA) and kainoid synthase (RadC) enzymes support a slightly altered DA biosynthetic model in C. armata via the congener isodomoic acid B, with RadC behaving more like the homologous diatom enzyme despite higher amino acid similarity to red algal kainic acid synthesis enzymes. A phylogenetic analysis of the rad genes suggests unique origins for the red macroalgal and diatom genes in their respective hosts, with native eukaryotic CYP450 neofunctionalization combining with the horizontal gene transfer of N-prenyltransferases and kainoid synthases to establish DA production within the algal lineages.

Identification and functional validation of super-enhancers in Arabidopsis thaliana.

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

A systematic comparison of natural product potential, with an emphasis on RiPPs, by mining of bacteria of three large ecosystems.

The implementation of several global microbiome studies has yielded extensive insights into the biosynthetic potential of natural microbial communities. However, studies on the distribution of several classes of ribosomally synthesized and post-translationally modified peptides (RiPPs), non-ribosomal peptides (NRPs) and polyketides (PKs) in different large microbial ecosystems have been very limited. Here, we collected a large set of metagenome-assembled bacterial genomes from marine, freshwater and terrestrial ecosystems to investigate the biosynthetic potential of these bacteria. We demonstrate the utility of public dataset collections for revealing the different secondary metabolite biosynthetic potentials among these different living environments. We show that there is a higher occurrence of RiPPs in terrestrial systems, while in marine systems, we found relatively more terpene-, NRP-, and PK encoding gene clusters. Among the many new biosynthetic gene clusters (BGCs) identified, we analyzed various Nif-11-like and nitrile hydratase leader peptide (NHLP) containing gene clusters that would merit further study, including promising products, such as mersacidin-, LAP- and proteusin analogs. This research highlights the significance of public datasets in elucidating the biosynthetic potential of microbes in different living environments and underscores the wide bioengineering opportunities within the RiPP family.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy.

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Plant terpene specialized metabolism: complex networks or simple linear pathways?

From the perspectives of pathway evolution, discovery and engineering of plant specialized metabolism, the nature of the biosynthetic routes represents a critical aspect. Classical models depict biosynthesis typically from an end-point angle and as linear, for example, connecting central and specialized metabolism. As the number of functionally elucidated routes increased, the enzymatic foundation of complex plant chemistries became increasingly well understood. The perception of linear pathway models has been severely challenged. With a focus on plant terpenoid specialized metabolism, we review here illustrative examples supporting that plants have evolved complex networks driving chemical diversification. The completion of several diterpene, sesquiterpene and monoterpene routes shows complex formation of scaffolds and their subsequent functionalization. These networks show that branch points, including multiple sub-routes, mean that metabolic grids are the rule rather than the exception. This concept presents significant implications for biotechnological production.

Unveiling biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2.

BACKGROUND: A growing number of studies have demonstrated that the polar regions have the potential to be a significant repository of microbial resources and a potential source of active ingredients. Genome mining strategy plays a key role in the discovery of bioactive secondary metabolites (SMs) from microorganisms. This work highlighted deciphering the biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2 by a combination of whole genome analysis and antiSMASH as well as feature-based molecular networking (MN) in the Global Natural Products Social Molecular Networking (GNPS). RESULTS: In this study, a high-quality whole genome sequence of an Arctic marine strain MNP-2, with a size of 34.9 Mb was successfully obtained. Its total number of genes predicted by BRAKER software was 13,218, and that of non-coding RNAs (rRNA, sRNA, snRNA, and tRNA) predicted by using INFERNAL software was 204. AntiSMASH results indicated that strain MNP-2 harbors 56 biosynthetic gene clusters (BGCs), including 18 NRPS/NRPS-like gene clusters, 10 PKS/PKS-like gene clusters, 8 terpene synthse gene clusters, 5 indole synthase gene clusters, 10 hybrid gene clusters, and 5 fungal-RiPP gene clusters. Metabolic analyses of strain MNP-2 grown on various media using GNPS networking revealed its great potential for the biosynthesis of bioactive SMs containing a variety of heterocyclic and bridge-ring structures. For example, compound G-8 exhibited a potent anti-HIV effect with an IC50 value of 7.2 nM and an EC50 value of 0.9 nM. Compound G-6 had excellent in vitro cytotoxicities against the K562, MCF-7, Hela, DU145, U1975, SGC-7901, A549, MOLT-4, and HL60 cell lines, with IC50 values ranging from 0.10 to 3.3 µM, and showed significant anti-viral (H1N1 and H3N2) activities with IC50 values of 15.9 and 30.0 µM, respectively. CONCLUSIONS: These findings definitely improve our knowledge about the molecular biology of genus A. sydowii and would effectively unveil the biosynthetic potential of strain MNP-2 using genomics and metabolomics techniques.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns.

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Restoration of the Functional nif Gene Cluster by Complex Recombination Events during Heterocyst Development in the Nitrogen-Fixing Cyanobacterium Calothrix sp. NIES-4101.

In the genome of the heterocystous cyanobacterium Calothrix sp. NIES-4101 (NIES-4101), the four genes essential for nitrogen fixation (nifB, nifH, nifD and nifK) are highly fragmented into 13 parts in a 350-kb chromosomal region, and four of these parts are encoded in the reverse strand. Such a complex fragmentation feature makes it difficult to restore the intact nifBHDK genes by the excision mechanism found in the nifD gene of the Anabaena sp. PCC 7120 heterocyst. To examine the nitrogen-fixing ability of NIES-4101, we confirmed that NIES-4101 grew well on a combined nitrogen-free medium and showed high nitrogenase activity, which strongly suggested that the complete nifBHDK genes are restored by a complex recombination process in heterocysts. Next, we resequenced the genome prepared from cells grown under nitrogen-fixing conditions. Two contigs covering the complete nifHDK and nifB genes were found by de novo assembly of the sequencing reads. In addition, the DNA fragments covering the nifBHDK operon were successfully amplified by PCR. We propose that the process of nifBHDK restoration occurs as follows. First, the nifD-nifK genes are restored by four excision events. Then, the complete nifH and nifB genes are restored by two excision events followed by two successive inversion events between the inverted repeat sequences and one excision event, forming the functional nif gene cluster, nifB-fdxN-nifS-nifU-nifH-nifD-nifK. All genes coding recombinases responsible for these nine recombination events are located close to the terminal repeat sequences. The restoration of the nifBHDK genes in NIES-4101 is the most complex genome reorganization reported in heterocystous cyanobacteria.

Creating large chromosomal segment deletions in Aspergillus flavus by a dual CRISPR/Cas9 system: Deletion of gene clusters for production of aflatoxin, cyclopiazonic acid, and ustiloxin B.

Aspergillus flavus produces hepatocarcinogenic aflatoxin that adversely impacts human and animal health and international trade. A promising means to manage preharvest aflatoxin contamination of crops is biological control, which employs non-aflatoxigenic A. flavus isolates possessing defective aflatoxin gene clusters to outcompete field toxigenic populations. However, these isolates often produce other toxic metabolites. The CRISPR/Cas9 technology has greatly advanced genome editing and gene functional studies. Its use in deleting large chromosomal segments of filamentous fungi is rarely reported. A system of dual CRISPR/Cas9 combined with a 60-nucleotide donor DNA that allowed removal of A. flavus gene clusters involved in production of harmful specialized metabolites was established. It efficiently deleted a 102-kb segment containing both aflatoxin and cyclopiazonic acid gene clusters from toxigenic A. flavus morphotypes, L-type and S-type. It further deleted the 27-kb ustiloxin B gene cluster of a resulting L-type mutant. Overall efficiencies of deletion ranged from 66.6 % to 85.6 % and efficiencies of deletions repaired by a single copy of donor DNA ranged from 50.5 % to 72.7 %. To determine the capacity of this technique, a pigment-screening setup based on absence of aspergillic acid gene cluster was devised. Chromosomal segments of 201 kb and 301 kb were deleted with efficiencies of 57.7 % to 69.2 %, respectively. This system used natural A. flavus isolates as recipients, eliminated a forced-recycling step to produce recipients for next round deletion, and generated maker-free deletants with sequences predefined by donor DNA. The research provides a method for creating genuine atoxigenic biocontrol strains friendly for field trial release.

The porcine skin microbiome exhibits broad fungal antagonism.

The skin and its microbiome function to protect the host from pathogen colonization and environmental stressors. In this study, using the Wisconsin Miniature Swine™ model, we characterize the porcine skin fungal and bacterial microbiomes, identify bacterial isolates displaying antifungal activity, and use whole-genome sequencing to identify biosynthetic gene clusters encoding for secondary metabolites that may be responsible for the antagonistic effects on fungi. Through this comprehensive approach of paired microbiome sequencing with culturomics, we report the discovery of novel species of Corynebacterium and Rothia. Further, this study represents the first comprehensive evaluation of the porcine skin mycobiome and the evaluation of bacterial-fungal interactions on this surface. Several diverse bacterial isolates exhibit potent antifungal properties against opportunistic fungal pathogens in vitro. Genomic analysis of inhibitory species revealed a diverse repertoire of uncharacterized biosynthetic gene clusters suggesting a reservoir of novel chemical and biological diversity. Collectively, the porcine skin microbiome represents a potential unique source of novel antifungals.