Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

Structural diversity of decalin forming Diels-Alderase

The Diels-Alder (DA) reaction, specifically referring to the [4 + 2] cycloaddition reaction in pericyclic reactions, is a process that forms two carbon-carbon covalent bonds in a single step via an electron ring transition state. Among the secondary metabolites produced by microorganisms, numerous compounds are biosynthesized through DA reactions, most of which are enzymatic. Our research group has discovered an enzyme named Diels-Alderase (DAase) that catalyzes the DA reaction in filamentous fungi, and we have been investigating its catalytic mechanism. This review describes the reported microbial DAase enzymes, with a particular focus on those involved in the construction of the decalin ring.

Revealing Hidden Genes in Botrytis cinerea: New Insights into Genes Involved in the Biosynthesis of Secondary Metabolites

Utilizing bioinformatics tools, this study expands our understanding of secondary metabolism in Botrytis cinerea, identifying novel genes within polyketide synthase (PKS), non-ribosomal peptide synthetase (NRPS), sesquiterpene cyclase (STC), diterpene cyclase (DTC), and dimethylallyltryptophan synthase (DMATS) families. These findings enrich the genetic framework associated with B. cinerea's pathogenicity and ecological adaptation, offering insights into uncharted metabolic pathways. Significantly, the discovery of previously unannotated genes provides new molecular targets for developing targeted antifungal strategies, promising to enhance crop protection and advance our understanding of fungal biochemistry. This research not only broadens the scope of known secondary metabolites but also opens avenues for future exploration into B. cinerea's biosynthetic capabilities, potentially leading to novel antifungal compounds. Our work underscores the importance of integrating bioinformatics and genomics for fungal research, paving the way for sustainable agricultural practices by pinpointing precise molecular interventions against B. cinerea. This study sets a foundation for further investigations into the fungus's secondary metabolism, with implications for biotechnology and crop disease management.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. In this manuscript, we present a genome for a laboratory-tractable Entomophthora muscae isolate that infects fruit flies. Our E. muscae assembly is 1.03 Gb, consists of 7810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with recent datasets from entomophthoralean fungi, we show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity. In addition, we find that E. muscae and its closest allies possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Finally, we offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex that support the need for taxonomic revision within this group. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. This is in part due to their large genome sizes and difficulty in obtaining sequenceable material. In this manuscript, we leveraged a recently-isolated, laboratory-tractable Entomophthora muscae isolate and improved long-read sequencing to obtain a largely-complete entomophthoralean genome. Our E. muscae assembly is 1.03 Gb, consists of 7,810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with other available (transcriptomic and genomic) datasets from entomophthoralean fungi, we provide new insight into the biology of these understudied pathogens. We offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex. Our findings suggest that substantial taxonomic revision is needed to define species within this group and we provide recommendations for differentiating strains and species in the context of the existing body of E. muscae scientific literature. We show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity, despite the presence of machinery to defend against transposable elements(RNAi). In addition, we find that E. muscae and its closest allies are enriched for M16A peptidases and possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We find that E. muscae has an expanded group of acid-trehalases, consistent with trehalose being the primary sugar component of fly (and insect) hemolymph. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Natural trait variation across Saccharomycotina species

Among molecular biologists, the group of fungi called Saccharomycotina is famous for its yeasts. These yeasts in turn are famous for what they have in common-genetic, biochemical, and cell-biological characteristics that serve as models for plants and animals. But behind the apparent homogeneity of Saccharomycotina species lie a wealth of differences. In this review, we discuss traits that vary across the Saccharomycotina subphylum. We describe cases of bright pigmentation; a zoo of cell shapes; metabolic specialties; and species with unique rules of gene regulation. We discuss the genetics of this diversity and why it matters, including insights into basic evolutionary principles with relevance across Eukarya.

Genomic Insights into Paucibacter aquatile DH15, a Cyanobactericidal Bacterium, and Comparative Genomics of the Genus Paucibacter

Microcystis blooms threaten ecosystem function and cause substantial economic losses. Microorganism-based methods, mainly using cyanobactericidal bacteria, are considered one of the most ecologically sound methods to control Microcystis blooms. This study focused on gaining genomic insights into Paucibacter aquatile DH15 that exhibited excellent cyanobactericidal effects against Microcystis. Additionally, a pan-genome analysis of the genus Paucibacter was conducted to enhance our understanding of the ecophysiological significance of this genus. Based on phylogenomic analyses, strain DH15 was classified as a member of the species Paucibacter aquatile. The genome analysis supported that strain DH15 can effectively destroy Microcystis, possibly due to the specific genes involved in the flagellar synthesis, cell wall degradation, and the production of cyanobactericidal compounds. The pan-genome analysis revealed the diversity and adaptability of the genus Paucibacter, highlighting its potential to absorb external genetic elements. Paucibacter species were anticipated to play a vital role in the ecosystem by potentially providing essential nutrients, such as vitamins B7, B12, and heme, to auxotrophic microbial groups. Overall, our findings contribute to understanding the molecular mechanisms underlying the action of cyanobactericidal bacteria against Microcystis and shed light on the ecological significance of the genus Paucibacter.

Comparative genomic and phenotypic analyses of pathogenic fungi Neoscytalidium dimidiatum and Bipolaris papendorfii isolated from human skin scraping

Neoscytalidium dimidiatum and Bipolaris species are fungal plant pathogens that have been reported to cause human diseases. Recently, we have isolated numerous N. dimidiatum and Bipolaris species from the skin scrapings and nails of different patients. In this work, we have sequenced the genome of one strain of N. dimidiatum. The sequenced genome was compared to that of a previously reported Bipolaris papendorfii genome for a better understanding of their complex lifestyle and broad host-range pathogenicity. Both N. dimidiatum UM 880 (~ 43 Mb) and B. papendorfii UM 226 (~ 33 Mb) genomes include 11,015-12,320 putative coding DNA sequences, of which 0.51-2.49% are predicted transposable elements. Analysis of secondary metabolism gene clusters revealed several genes involved in melanin biosynthesis and iron uptake. The arsenal of CAZymes related to plants pathogenicity is comparable between the species, including genes involved in hemicellulose and pectin decomposition. Several important gene encoding keratinolytic peptidases were identified in N. dimidiatum and B. papendorfii, reflecting their potential pathogenic role in causing skin and nail infections. In this study, additional information on the metabolic features of these two species, such as nutritional profiling, pH tolerance, and osmotolerant, are revealed. The genomic characterization of N. dimidiatum and B. papendorfii provides the basis for the future functional studies to gain further insights as to what makes these fungi persist in plants and why they are pathogenic to humans.

Biomimetic engineering of nonribosomal peptide synthesis

Nonribosomal peptides (NRPs) have gained attention due to their diverse biological activities and potential applications in medicine and agriculture. The natural diversity of NRPs is a result of evolutionary processes that have occurred over millions of years. Recent studies have shed light on the mechanisms by which nonribosomal peptide synthetases (NRPSs) evolve, including gene duplication, recombination, and horizontal transfer. Mimicking natural evolution could be a useful strategy for engineering NRPSs to produce novel compounds with desired properties. Furthermore, the emergence of antibiotic-resistant bacteria has highlighted the urgent need for new drugs, and NRPs represent a promising avenue for drug discovery. This review discusses the engineering potential of NRPSs in light of their evolutionary history.

Nonribosomal Peptide Synthetases in Animals

Nonribosomal peptide synthetases (NRPSs) are a class of cytosolic enzymes that synthesize a range of bio-active secondary metabolites including antibiotics and siderophores. They are widespread among both prokaryotes and eukaryotes but are considered rare among animals. Recently, several novel NRPS genes have been described in nematodes, schistosomes, and arthropods, which led us to investigate how prevalent NRPS genes are in the animal kingdom. We screened 1059 sequenced animal genomes and showed that NRPSs were present in 7 out of the 19 phyla analyzed. A phylogenetic analysis showed that the identified NRPSs form clades distinct from other adenylate-forming enzymes that contain similar domains such as fatty acid synthases. NRPSs show a remarkably scattered distribution over the animal kingdom. They are especially abundant in rotifers and nematodes. In rotifers, we found a large variety of domain architectures and predicted substrates. In the nematode Plectus sambesii, we identified the beta-lactam biosynthesis genes L-δ-(α-aminoadipoyl)-L-cysteinyl-D-valine synthetase, isopenicillin N synthase, and deacetoxycephalosporin C synthase that catalyze the formation of beta-lactam antibiotics in fungi and bacteria. These genes are also present in several species of Collembola, but not in other hexapods analyzed so far. In conclusion, our survey showed that NRPS genes are more abundant and widespread in animals than previously known.

Mining for a new class of fungal natural products: the evolution, diversity, and distribution of isocyanide synthase biosynthetic gene clusters

The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We amalgamated a pipeline of tools to predict BGCs based on shared promoter motifs and located 3800 ICS BGCs in 3300 genomes, making ICS BGCs the fifth largest class of specialized metabolites compared to canonical classes found by antiSMASH. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1/2 gene cluster family (GCF), which was prior only studied in yeast, is present in ∼30% of all Ascomycetes. The dit variety ICS exhibits greater similarity to bacterial ICS than other fungal ICS, suggesting a potential convergence of the ICS backbone domain. The evolutionary origins of the dit GCF in Ascomycota are ancient and these genes are diversifying in some lineages. Our results create a roadmap for future research into ICS BGCs. We developed a website (https://isocyanides.fungi.wisc.edu/) that facilitates the exploration and downloading of all identified fungal ICS BGCs and GCFs.

Inventory of the Secondary Metabolite Biosynthetic Potential of Members within the Terminal Clade of the Fusarium solani Species Complex

The Fusarium solani species complex (FSSC) constitutes at least 77 phylogenetically distinct species including several agriculturally important and clinically relevant opportunistic pathogens. As with other Fusaria, they have been well documented to produce many secondary metabolites-compounds that are not required for the fungus to grow or develop but may be beneficial to the organism. An analysis of ten genomes from fungi within the terminal clade (clade 3) of the FSSC revealed each genome encoded 35 (F. cucurbitcola) to 48 (F. tenucristatum) secondary metabolite biosynthetic gene clusters (BGCs). A total of seventy-four different BGCs were identified from the ten FSSC genomes including seven polyketide synthases (PKS), thirteen nonribosomal peptide synthetases (NRPS), two terpene synthase BGCs, and a single dimethylallytryptophan synthase (DMATS) BGC conserved in all the genomes. Some of the clusters that were shared included those responsible for producing naphthoquinones such as fusarubins, a red pigmented compound, squalestatin, and the siderophores malonichrome, ferricrocin, and triacetylfusarinine. Eight novel NRPS and five novel PKS BGCs were identified, while BGCs predicted to produce radicicol, gibberellin, and fusaoctaxin were identified, which have not previously described in members of the FSSC. The diversity of the secondary metabolite repertoire of the FSSC may contribute to the expansive host range of these fungi and their ability to colonize broad habitats.

From Genes to Molecules, Secondary Metabolism in Botrytis cinerea: New Insights into Anamorphic and Teleomorphic Stages

The ascomycete Botrytis cinerea Pers. Fr., classified within the family Sclerotiniaceae, is the agent that causes grey mould disease which infects at least 1400 plant species, including crops of economic importance such as grapes and strawberries. The life cycle of B. cinerea consists of two phases: asexual (anamorph, Botrytis cinerea Pers. Fr.) and sexual (teleomorph, Botryotinia fuckeliana (de Bary) Wetzel). During the XVI International Symposium dedicated to the Botrytis fungus, which was held in Bari in June 2013, the scientific community unanimously decided to assign the most widely used name of the asexual form, Botrytis, to this genus of fungi. However, in the literature, we continue to find articles referring to both morphic stages. In this review, we take stock of the genes and metabolites reported for both morphic forms of B. cinerea between January 2015 and October 2022.

Classification and Multifaceted Potential of Secondary Metabolites Produced by Bacillus subtilis Group: A Comprehensive Review

Despite their remarkable biosynthetic potential, Bacillus subtilis have been widely overlooked. However, their capability to withstand harsh conditions (extreme temperature, Ultraviolet (UV) and γ-radiation, and dehydration) and the promiscuous metabolites they synthesize have created increased commercial interest in them as a therapeutic agent, a food preservative, and a plant-pathogen control agent. Nevertheless, the commercial-scale availability of these metabolites is constrained due to challenges in their accessibility via synthesis and low fermentation yields. In the context of this rising in interest, we comprehensively visualized the antimicrobial peptides produced by B. subtilis and highlighted their prospective applications in various industries. Moreover, we proposed and classified these metabolites produced by the B. subtilis group based on their biosynthetic pathways and chemical structures. The biosynthetic pathway, bioactivity, and chemical structure are discussed in detail for each class. We believe that this review will spark a renewed interest in the often disregarded B. subtilis and its remarkable biosynthetic capabilities.

The Botrytis cinerea Gene Expression Browser

For comprehensive gene expression analyses of the phytopathogenic fungus Botrytis cinerea, which infects a number of plant taxa and is a cause of substantial agricultural losses worldwide, we developed BEB, a web-based B. cinerea gene Expression Browser. This computationally inexpensive web-based application and its associated database contain manually curated RNA-Seq data for B. cinerea. BEB enables expression analyses of genes of interest under different culture conditions by providing publication-ready heatmaps depicting transcript levels, without requiring advanced computational skills. BEB also provides details of each experiment and user-defined gene expression clustering and visualization options. If needed, tables of gene expression values can be downloaded for further exploration, including, for instance, the determination of differentially expressed genes. The BEB implementation is based on open-source computational technologies that can be deployed for other organisms. In this case, the new implementation will be limited only by the number of transcriptomic experiments that are incorporated into the platform. To demonstrate the usability and value of BEB, we analyzed gene expression patterns across different conditions, with a focus on secondary metabolite gene clusters, chromosome-wide gene expression, previously described virulence factors, and reference genes, providing the first comprehensive expression overview of these groups of genes in this relevant fungal phytopathogen. We expect this tool to be broadly useful in B. cinerea research, providing a basis for comparative transcriptomics and candidate gene identification for functional assays.

Diversity and Biosynthetic Potential of Fungi Isolated from St. John's Island, Singapore

Adaptation to a wide variety of habitats allows fungi to develop unique abilities to produce diverse secondary metabolites with diverse bioactivities. In this study, 30 Ascomycetes fungi isolated from St. John's Island, Singapore were investigated for their general biosynthetic potential and their ability to produce antimicrobial secondary metabolites (SMs). All the 30 fungal isolates belong to the Phylum Ascomycota and are distributed into 6 orders and 18 genera with Order Hypocreales having the highest number of representative (37%). Screening for polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) genes using degenerate PCR led to the identification of 23 polyketide synthases (PKSs) and 5 nonribosomal peptide synthetases (NRPSs) grouped into nine distinct clades based on their reduction capabilities. Some of the identified PKSs genes share high similarities between species and known reference genes, suggesting the possibility of conserved biosynthesis of closely related compounds from different fungi. Fungal extracts were tested for their antimicrobial activity against S. aureus, Methicillin-resistant S. aureus (MRSA), and Candida albicans. Bioassay-guided fractionation of the active constituents from two promising isolates resulted in the isolation of seven compounds: Penilumamides A, D, and E from strain F4335 and xanthomegnin, viomellein, pretrichodermamide C and vioxanthin from strain F7180. Vioxanthin exhibited the best antibacterial activity with IC₅₀ values of 3.0 μM and 1.6 μM against S. aureus and MRSA respectively. Viomellein revealed weak antiproliferative activity against A549 cells with an IC₅₀ of 42 μM. The results from this study give valuable insights into the diversity and biosynthetic potential of fungi from this unique habitat and forms a background for an in-depth analysis of the biosynthetic capability of selected strains of interest with the aim of discovering novel fungal natural products.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Identification of a novel secreted metabolite cyclo(phenylalanyl-prolyl) from Batrachochytrium dendrobatidis and its effect on Galleria mellonella

BACKGROUND

The fungus, Batrachochytrium dendrobatidis, is the causative agent of chytridiomycosis and a leading cause of global decline in amphibian populations. The first stages of chytridiomycosis include: inflammation, hyperkeratosis, lethargy, loss of righting reflex, and disruption of internal electrolyte levels leading to eventual death of the host. Previous work indicates that B. dendrobatidis can produce immunomodulatory compounds and other secreted molecules that regulate the growth of the fungus. In this study, filtrates of the fungus grown in media and water were subjected to ultra-performance liquid chromatography-mass spectrometry and analyzed using Compound Discoverer 3.0.

RESULTS

Identification of cyclo(phenylalanyl-prolyl), chitobiose, and S-adenosylmethionine were verified by their retention times and fragmentation patterns from B. dendrobatidis supernatants. Previous studies have analyzed the effects of B. dendrobatidis on amphibian models, in vitro, or in cell culture. We studied the effects of live B. dendrobatidis cells, spent culture filtrates containing secreted metabolites, and cyclo(pheylalanyl-prolyl) on wax moth larvae (Galleria mellonella). Concentrated filtrates caused melanization within 24 h, while live B. dendrobatidis caused melanization within 48 h.

CONCLUSIONS

Here we show B. dendrobatidis produces secreted metabolites previously unreported. The impacts of these chemicals were tested on an alternate non-amphibian model system that has been used for other fungi to study pathogenicity traits in this fungus.

Evolutionary relationships of adenylation domains in fungi

Non-ribosomal peptide synthetases (NRPSs) and NRPS-like enzymes are abundant in microbes as they are involved in the production of primary and secondary metabolites. In contrast to the well-studied NRPSs, known to produce non-ribosomal peptides, NRPS-like enzymes exhibit more diverse activities and their evolutionary relationships are unclear. Here, we present the first in-depth phylogenetic analysis of fungal NRPS-like A domains from functionally characterized pathways, and their relationships to characterized A domains found in fungal NRPSs. This study clearly differentiated amino acid reductases, including NRPSs, from CoA/AMP ligases, which could be divided into 10 distinct phylogenetic clades that reflect their conserved domain organization, substrate specificity and enzymatic activity. In particular, evolutionary relationships of adenylate forming reductases could be refined and explained the substrate specificity difference. Consistent with their phylogeny, the deduced amino acid code of A domains differentiated amino acid reductases from other enzymes. However, a diagnostic code was found for α-keto acid reductases and clade 7 CoA/AMP ligases only. Comparative genomics of loci containing these enzymes revealed that they can be independently recruited as tailoring genes in diverse secondary metabolite pathways. Based on these results, we propose a refined and clear phylogeny-based classification of A domain-containing enzymes, which will provide a robust framework for future functional analyses and engineering of these enzymes to produce new bioactive molecules.

Subcellular localization of fungal specialized metabolites

Fungal specialized metabolites play an important role in the environment and have impacted human health and survival significantly. These specialized metabolites are often the end product of a series of sequential and collaborating biosynthetic enzymes that reside within different subcellular compartments. A wide variety of methods have been developed to understand fungal specialized metabolite biosynthesis in terms of the chemical conversions and the biosynthetic enzymes required, however there are far fewer studies elucidating the compartmentalization of the same enzymes. This review illustrates the biosynthesis of specialized metabolites where the localization of all, or some, of the biosynthetic enzymes have been determined and describes the methods used to identify the sub-cellular localization.

Characterization of NRPS and PKS genes involved in the biosynthesis of SMs in Alternaria dauci including the phytotoxic polyketide aldaulactone

Alternaria dauci is a Dothideomycete fungus, causal agent of carrot leaf blight. As a member of the Alternaria genus, known to produce a lot of secondary metabolite toxins, A. dauci is also supposed to synthetize host specific and non-host specific toxins playing a crucial role in pathogenicity. This study provides the first reviewing of secondary metabolism genetic basis in the Alternaria genus by prediction of 55 different putative core genes. Interestingly, aldaulactone, a phytotoxic benzenediol lactone from A. dauci, was demonstrated as important in pathogenicity and in carrot partial resistance to this fungus. As nothing is known about aldaulactone biosynthesis, bioinformatic analyses on a publicly available A. dauci genome data set that were reassembled, thanks to a transcriptome data set described here, allowed to identify 19 putative secondary metabolism clusters. We exploited phylogeny to pinpoint cluster 8 as a candidate in aldaulactone biosynthesis. This cluster contains AdPKS7 and AdPKS8, homologs with genes encoding a reducing and a non-reducing polyketide synthase. Clusters containing such a pair of PKS genes have been identified in the biosynthesis of resorcylic acid lactones or dihydroxyphenylacetic acid lactones. AdPKS7 and AdPKS8 gene expression patterns correlated with aldaulactone production in different experimental conditions. The present results highly suggest that both genes are responsible for aldaulactone biosynthesis.

Recent advances in the chemo-biological characterization of decalin natural products and unraveling of the workings of Diels-Alderases

The Diels-Alder (DA) reaction refers to a [4 + 2] cycloaddition reaction that falls under the category of pericyclic reactions. It is a reaction that allows regio- and stereo-selective construction of two carbon-carbon bonds simultaneously in a concerted manner to generate a six-membered ring structure through a six-electron cyclic transition state. The DA reaction is one of the most widely applied reactions in organic synthesis, yet its role in biological systems has been debated intensely over the last four decades. A survey of secondary metabolites produced by microorganisms suggests strongly that many of the compounds possess features that are likely formed through DA reactions, and most of them are considered to be catalyzed by enzymes that are commonly referred to as Diels-Alderases (DAases). In recent years, especially over the past 10 years or so, we have seen an accumulation of a substantial body of work that substantiates the argument that DAases indeed exist and play a critical role in the biosynthesis of complex metabolites. This review will cover the DAases involved in the biosynthesis of decalin moieties, which are found in many of the medicinally important natural products, especially those produced by fungi. In particular, we will focus on a subset of secondary metabolites referred to as pyrrolidine-2-one-bearing decalin compounds and discuss the decalin ring stereochemistry and the biological activities of those compounds. We will also look into the genes and enzymes that drive the biosynthetic construction of those complex natural products, and highlight the recent progress made on the structural and mechanistic understanding of DAases, especially regarding how those enzymes exert stereochemical control over the [4 + 2] cycloaddition reactions they catalyze.

A Secondary Metabolism Pathway Involved in the Production of a Putative Toxin Is Expressed at Early Stage of Monilinia laxa Infection

The necrotrophic pathogenic fungus Monilinia laxa causes brown rot disease on stone fruit generating significant yield losses. So far, a limited number of pathogenesis-related virulence factors, such as cell wall degrading enzymes and potential phytotoxins, have been described in Monilinia spp. Using RNA-sequencing data from highly virulent M. laxa ML8L strain at early stages of the infection process (6, 14, 24, and 48 h post-inoculation, hpi) on nectarine and the Pathogen-Host-Interactions (PHI) database, we selected a number of genes for further study and ranked them according to their transcription levels. We identified a class of genes highly expressed at 6 hpi and that their expression decreased to almost undetectable levels at 14 to 48 hpi. Among these genes we found Monilinia__061040 encoding a non-ribosomal peptide synthase (NRPS). Monilinia__061040 together with other five co-regulated genes, forms a secondary metabolism cluster potentially involved in the production of epipolythiodioxopiperazine (ETP) toxin. Quantitative-PCR data confirmed previous RNA sequencing results from the virulent ML8L strain. Interestingly, in a less virulent M. laxa ML5L strain the expression levels of this pathway were reduced compared to the ML8L strain during nectarine infection. In vitro experiments showed that liquid medium containing peach extract mimicked the results observed using nectarines. In fact, upregulation of the NRPS coding gene was also observed in minimal medium suggesting the existence of a fruit-independent mechanism of regulation for this putative toxin biosynthetic pathway that is also downregulated in the less virulent strain. These results emphasize the role of this secondary metabolism pathway during the early stage of brown rot disease development and show alternative models to study the induction of virulence genes in this fungus.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Analysis of Biosynthetic Gene Clusters, Secretory, and Antimicrobial Peptides Reveals Environmental Suitability of Exiguobacterium profundum PHM11

Halotolerant bacteria produce a wide range of bioactive compounds with important applications in agriculture for abiotic stress amelioration and plant growth promotion. In the present study, 17 biosynthetic gene clusters (BGCs) were identified in Exiguobacterium profundum PHM11 belonging to saccharides, desmotamide, pseudaminic acid, dipeptide aldehydes, and terpene biosynthetic pathways representing approximately one-sixth of genomes. The terpene biosynthetic pathway was conserved in Exiguobacterium spp. while the E. profundum PHM11 genome confirms the presence of the 1-deoxy-d-xylulose 5-phosphate (DXP) pathway for the isopentenyl diphosphate (IPP) synthesis. Further, 2,877 signal peptides (SPs) were identified using the PrediSi server, out of which 592 proteins were prophesied for the secretion having a transmembrane helix (TMH). In addition, antimicrobial peptides (AMPs) were also identified using BAGEL4. The transcriptome analysis of PHM11 under salt stress reveals the differential expression of putative secretion and transporter genes having SPs and TMH. Priming of the rice, wheat and maize seeds with PHM11 under salt stress led to improvement in the root length, root diameters, surface area, number of links and forks, and shoot length. The study shows that the presence of BGCs, SPs, and secretion proteins constituting TMH and AMPs provides superior competitiveness in the environment and make E. profundum PHM11 a suitable candidate for plant growth promotion under salt stress.

A Bioinformatics Workflow for Investigating Fungal Biosynthetic Gene Clusters

Predicting secondary metabolite biosynthetic gene clusters is a routine analysis performed for each newly sequenced fungal genome. Yet, the usefulness of such predictions remains restricted as they provide total numbers of biosynthetic pathways with only very limited biological significance. In this chapter, we describe a workflow to predict and analyze biosynthetic gene clusters in fungal genomes. It relies on similarity networking and phylogeny to perform genetic dereplication and to prioritize candidate gene clusters that potentially produce new compounds. This basic workflow includes the generation of high-quality figures for publication.

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.

Genome-Wide Analysis of Biosynthetic Gene Cluster Reveals Correlated Gene Loss with Absence of Usnic Acid in Lichen-Forming Fungi

Lichen-forming fungi are known to produce a large number of secondary metabolites. Some metabolites are deposited in the cortical layer of the lichen thallus where they exert important ecological functions, such as UV filtering. The fact that closely related lineages of lichen-forming fungi can differ in cortical chemistry suggests that natural product biosynthesis in lichens can evolve independent from phylogenetic constraints. Usnic acid is one of the major cortical pigments in lichens. Here we used a comparative genomic approach on 46 lichen-forming fungal species of the Lecanoromycetes to elucidate the biosynthetic gene content and evolution of the gene cluster putatively responsible for the biosynthesis of usnic acid. Whole-genome sequences were gathered from taxa belonging to different orders and families of Lecanoromycetes, where Parmeliaceae is the most well-represented taxon, and analyzed with a variety of genomic tools. The highest number of biosynthetic gene clusters was found in Evernia prunastri, Pannoparmelia angustata, and Parmotrema austrosinense, respectively, and lowest in Canoparmelia nairobiensis, Bulbothrix sensibilis, and Hypotrachyna scytodes. We found that all studied species producing usnic acid contain the putative usnic acid biosynthetic gene cluster, whereas the cluster was absent in all genomes of species lacking usnic acid. The absence of the gene cluster was supported by an additional unsuccessful search for ß-ketoacylsynthase, the most conserved domain of the gene cluster, in the genomes of species lacking usnic acid. The domain architecture of this PKS cluster—homologous to the already known usnic acid PKS cluster (MPAS) and CYT450 (MPAO)—varies within the studied species, whereas the gene arrangement is highly similar in closely related taxa. We hypothesize that the ancestor of these lichen-forming fungi contained the putative usnic acid producing PKS cluster and that the gene cluster was lost repeatedly during the evolution of these groups. Our study provides insight into the genomic adaptations to the evolutionary success of these lichen-forming fungal species and sets a baseline for further exploration of biosynthetic gene content and its evolutionary significance.

Deletion of the Bcnrps1 Gene Increases the Pathogenicity of Botrytis cinerea and Reduces Its Tolerance to the Exogenous Toxic Substances Spermidine and Pyrimethanil

During the infection of grapevine (Vitis vinifera) by the fungus Botrytis cinerea, the concentration of polyamines, which are toxic substances for the phytopathogen, increases in the grape. Nine NRPS genes have been identified in the genome of B. cinerea, yet the function of five of them remains unknown. For this reason, we have studied the expression of the 9 NRPS genes by RT-qPCR in a medium supplemented with sublethal concentrations of three polyamines (1,3-diaminopropane (1,3-DAP), spermidine (SPD), and spermine (SPM)). Our results show that the presence of polyamines in the culture medium triggered the overexpression of the Bcnrps1 gene in the pathogen. Deleting Bcnrps1 did not affect mycelial growth or adaptation to osmotic stress, and we show that its expression is not essential for the cycle of infection of the B. cinerea. However, mutating the Bcnrps1 gene resulted in overexpression of the Bcnrps6 gene, which encodes for the excretion of siderophores of the coprogen family. Moreover, gene deletion has reduced the tolerance of B. cinerea B05.10 to toxic substances such as the polyamine SPD and the fungicide pyrimethanil, and its virulence has increased. Our findings provide new insights into the function of the Bcnrps1 gene and its involvement in the tolerance of B. cinerea against exogenous toxic compounds.

Cyclic, Hydrophobic Hexapeptide Fusahexin Is the Product of a Nonribosomal Peptide Synthetase in Fusarium graminearum

The plant pathogenic fungus Fusarium graminearum is known to produce a wide array of secondary metabolites during plant infection. This includes several nonribosomal peptides. Recently, the fusaoctaxin (NRPS5/9) and gramilin (NRPS8) gene clusters were shown to be induced by host interactions. To widen our understanding of this important pathogen, we investigated the involvement of the NRPS4 gene cluster during infection and oxidative and osmotic stress. Overexpression of NRPS4 led to the discovery of a new cyclic hexapeptide, fusahexin (1), with the amino acid sequence cyclo-(d-Ala-l-Leu-d-allo-Thr-l-Pro-d-Leu-l-Leu). The structural analyses revealed an unusual ether bond between a proline C^δ to C^β of the preceding threonine resulting in an oxazine ring system. The comparative genomic analyses showed that the small gene cluster only encodes an ABC transporter in addition to the five-module nonribosomal peptide synthetase (NRPS). Based on the structure of fusahexin and the domain architecture of NRPS4, we propose a biosynthetic model in which the terminal module is used to incorporate two leucine units. So far, iterative use of NRPS modules has primarily been described for siderophore synthetases, which makes NRPS4 a rare example of a fungal nonsiderophore NRPS with distinct iterative module usage.

Rhodococcus comparative genomics reveals a phylogenomic-dependent non-ribosomal peptide synthetase distribution: insights into biosynthetic gene cluster connection to an orphan metabolite

Natural products (NPs) are synthesized by biosynthetic gene clusters (BGCs), whose genes are involved in producing one or a family of chemically related metabolites. Advances in comparative genomics have been favourable for exploiting huge amounts of data and discovering previously unknown BGCs. Nonetheless, studying distribution patterns of novel BGCs and elucidating the biosynthesis of orphan metabolites remains a challenge. To fill this knowledge gap, our study developed a pipeline for high-quality comparative genomics for the actinomycete genus Rhodococcus , which is metabolically versatile, yet understudied in terms of NPs, leading to a total of 110 genomes, 1891 BGCs and 717 non-ribosomal peptide synthetases (NRPSs). Phylogenomic inferences showed four major clades retrieved from strains of several ecological habitats. BiG-SCAPE sequence similarity BGC networking revealed 44 unidentified gene cluster families (GCFs) for NRPS, which presented a phylogenomic-dependent evolution pattern, supporting the hypothesis of vertical gene transfer. As a proof of concept, we analysed in-depth one of our marine strains, Rhodococcus sp. H-CA8f, which revealed a unique BGC distribution within its phylogenomic clade, involved in producing a chloramphenicol-related compound. While this BGC is part of the most abundant and widely distributed NRPS GCF, corason analysis unveiled major differences regarding its genetic context, co-occurrence patterns and modularity. This BGC is composed of three sections, two well-conserved right/left arms flanking a very variable middle section, composed of nrps genes. The presence of two non-canonical domains in H-CA8f’s BGC may contribute to adding chemical diversity to this family of NPs. Liquid chromatography-high resolution MS and dereplication efforts retrieved a set of related orphan metabolites, the corynecins, which to our knowledge are reported here for the first time in Rhodococcus . Overall, our data provide insights to connect BGC uniqueness with orphan metabolites, by revealing key comparative genomic features supported by models of BGC distribution along phylogeny.

Genomic Characteristics and Comparative Genomics Analysis of Two Chinese Corynespora cassiicola Strains Causing Corynespora Leaf Fall (CLF) Disease

Rubber tree Corynespora leaf fall (CLF) disease, caused by the fungus Corynespora cassiicola, is one of the most damaging diseases in rubber tree plantations in Asia and Africa, and this disease also threatens rubber nurseries and young rubber plantations in China. C. cassiicola isolates display high genetic diversity, and virulence profiles vary significantly depending on cultivar. Although one phytotoxin (cassicolin) has been identified, it cannot fully explain the diversity in pathogenicity between C. cassiicola species, and some virulent C. cassiicola strains do not contain the cassiicolin gene. In the present study, we report high-quality gapless genome sequences, obtained using short-read sequencing and single-molecule long-read sequencing, of two Chinese C. cassiicola virulent strains. Comparative genomics of gene families in these two stains and a virulent CPP strain from the Philippines showed that all three strains experienced different selective pressures, and metabolism-related gene families vary between the strains. Secreted protein analysis indicated that the quantities of secreted cell wall-degrading enzymes were correlated with pathogenesis, and the most aggressive CCP strain (cassiicolin toxin type 1) encoded 27.34% and 39.74% more secreted carbohydrate-active enzymes (CAZymes) than Chinese strains YN49 and CC01, respectively, both of which can only infect rubber tree saplings. The results of antiSMASH analysis showed that all three strains encode ~60 secondary metabolite biosynthesis gene clusters (SM BGCs). Phylogenomic and domain structure analyses of core synthesis genes, together with synteny analysis of polyketide synthase (PKS) and non-ribosomal peptide synthetase (NRPS) gene clusters, revealed diversity in the distribution of SM BGCs between strains, as well as SM polymorphisms, which may play an important role in pathogenic progress. The results expand our understanding of the C. cassiicola genome. Further comparative genomic analysis indicates that secreted CAZymes and SMs may influence pathogenicity in rubber tree plantations. The findings facilitate future exploration of the molecular pathogenic mechanism of C. cassiicola.

Rational engineering of specialized metabolites in bacteria and fungi

Bacteria and fungi have a high potential to produce compounds that display large structural change and diversity, thus displaying an extensive range of biological activities. Secondary metabolism or specialized metabolism is a term for pathways and small molecule products of metabolism that are not mandatory for the subsistence of the organism but improve and control their phenotype. Their interesting biological activities have occasioned their application in the fields of agriculture, food, and pharmaceuticals. Metabolic engineering is a powerful approach to improve access to these treasured molecules or to rationally engineer new ones. A thorough overview of engineering methods in secondary metabolism is presented, both in heterologous and epigenetic modification. Engineering methods to modify the structure of some secondary metabolite classes in their host are also intensively assessed.

Terpenoid Biosynthesis Dominates among Secondary Metabolite Clusters in Mucoromycotina Genomes

Early-diverging fungi harbour unprecedented diversity in terms of living forms, biological traits and genome architecture. Before the sequencing era, non-Dikarya fungi were considered unable to produce secondary metabolites (SM); however, this perspective is changing. The main classes of secondary metabolites in fungi include polyketides, nonribosomal peptides, terpenoids and siderophores that serve different biological roles, including iron chelation and plant growth promotion. The same classes of SM are reported for representatives of early-diverging fungal lineages. Encouraged by the advancement in the field, we carried out a systematic survey of SM in Mucoromycotina and corroborated the presence of various SM clusters (SMCs) within the phylum. Among the core findings, considerable representation of terpene and nonribosomal peptide synthetase (NRPS)-like candidate SMCs was found. Terpene clusters with diverse domain composition and potentially highly variable products dominated the landscape of candidate SMCs. A uniform low-copy distribution of siderophore clusters was observed among most assemblies. Mortierellomycotina are highlighted as the most potent SMC producers among the Mucoromycota and as a source of novel peptide products. SMC identification is dependent on gene model quality and can be successfully performed on a batch scale with genomes of different quality and completeness.

Bacterial-Like Nonribosomal Peptide Synthetases Produce Cyclopeptides in the Zygomycetous Fungus Mortierella alpina

Fungi are traditionally considered a reservoir of biologically active natural products. However, an active secondary metabolism has long not been attributed to early-diverging fungi such as Mortierella Here, we report on the biosynthesis of two series of cyclic pentapeptides, the malpicyclins and malpibaldins, as products of Mortierella alpina ATCC 32222. The molecular structures of malpicyclins were elucidated by high-resolution tandem mass spectrometry (HR-MS/MS), Marfey's method, and one-dimensional (1D) and 2D nuclear magnetic resonance (NMR) spectroscopy. In addition, malpibaldin biosynthesis was confirmed by HR-MS. Genome mining and comparative quantitative real-time PCR (qRT-PCR) expression analysis pointed at two pentamodular nonribosomal peptide synthetases (NRPSs), malpicyclin synthetase MpcA and malpibaldin synthetase MpbA, as candidate biosynthetic enzymes. Heterologous production of the respective adenylation domains and substrate specificity assays proved promiscuous substrate selection and confirmed their respective biosynthetic roles. In stark contrast to known fungal NRPSs, MpbA and MpcA contain bacterial-like dual epimerase/condensation domains allowing the racemization of enzyme-tethered l-amino acids and the subsequent incorporation of d-amino acids into the metabolites. Phylogenetic analyses of both NRPS genes indicated a bacterial origin and a horizontal gene transfer into the fungal genome. We report on the as-yet-unexplored nonribosomal peptide biosynthesis in basal fungi which highlights this paraphylum as a novel and underrated resource of natural products.IMPORTANCE Fungal natural compounds are industrially produced, with application in antibiotic treatment, cancer medications, and crop plant protection. Traditionally, higher fungi have been intensively investigated concerning their metabolic potential, but reidentification of already known compounds is frequently observed. Hence, alternative strategies to acquire novel bioactive molecules are required. We present the genus Mortierella as representative of the early-diverging fungi as an underestimated resource of natural products. Mortierella alpina produces two families of cyclopeptides, designated malpicyclins and malpibaldins, respectively, via two pentamodular nonribosomal peptide synthetases (NRPSs). These enzymes are much more closely related to bacterial than to other fungal NRPSs, suggesting a bacterial origin of these NRPS genes in Mortierella Both enzymes were biochemically characterized and are involved in as-yet-unknown biosynthetic pathways of natural products in basal fungi. Hence, this report establishes early-diverging fungi as prolific natural compound producers and sheds light on the origin of their biosynthetic capacity.

BiG-SLiCE: A highly scalable tool maps the diversity of 1.2 million biosynthetic gene clusters

BACKGROUND

Genome mining for biosynthetic gene clusters (BGCs) has become an integral part of natural product discovery. The >200,000 microbial genomes now publicly available hold information on abundant novel chemistry. One way to navigate this vast genomic diversity is through comparative analysis of homologous BGCs, which allows identification of cross-species patterns that can be matched to the presence of metabolites or biological activities. However, current tools are hindered by a bottleneck caused by the expensive network-based approach used to group these BGCs into gene cluster families (GCFs).

RESULTS

Here, we introduce BiG-SLiCE, a tool designed to cluster massive numbers of BGCs. By representing them in Euclidean space, BiG-SLiCE can group BGCs into GCFs in a non-pairwise, near-linear fashion. We used BiG-SLiCE to analyze 1,225,071 BGCs collected from 209,206 publicly available microbial genomes and metagenome-assembled genomes within 10 days on a typical 36-core CPU server. We demonstrate the utility of such analyses by reconstructing a global map of secondary metabolic diversity across taxonomy to identify uncharted biosynthetic potential. BiG-SLiCE also provides a "query mode" that can efficiently place newly sequenced BGCs into previously computed GCFs, plus a powerful output visualization engine that facilitates user-friendly data exploration.

CONCLUSIONS

BiG-SLiCE opens up new possibilities to accelerate natural product discovery and offers a first step towards constructing a global and searchable interconnected network of BGCs. As more genomes are sequenced from understudied taxa, more information can be mined to highlight their potentially novel chemistry. BiG-SLiCE is available via https://github.com/medema-group/bigslice.

Fusarium Cyclodepsipeptide Mycotoxins: Chemistry, Biosynthesis, and Occurrence

Most of the fungi from the Fusarium genus are pathogenic to cereals, vegetables, and fruits and the products of their secondary metabolism mycotoxins may accumulate in foods and feeds. Non-ribosomal cyclodepsipeptides are one of the main mycotoxin groups and include beauvericins (BEAs), enniatins (ENNs), and beauvenniatins (BEAEs). When ingested, even small amounts of these metabolites significantly affect human and animal health. On the other hand, in view of their antimicrobial activities and cytotoxicity, they may be used as components in drug discovery and processing and are considered as suitable candidates for anti-cancer drugs. Therefore, it is crucial to expand the existing knowledge about cyclodepsipeptides and to search for new analogues of these compounds. The present manuscript aimed to highlight the extensive variability of cyclodepsipeptides by describing chemistry, biosynthesis, and occurrence of BEAs, ENNs, and BEAEs in foods and feeds. Moreover, the co-occurrence of Fusarium species was compared to the amounts of toxins in crops, vegetables, and fruits from different regions of the world.

Divergence of Beauvericin Synthase Gene among Fusarium and Trichoderma Species

Beauvericin (BEA) is a cyclodepsipeptide mycotoxin, showing insecticidal, antibiotic and antimicrobial activities, as well as inducing apoptosis of cancer cell lines. BEA can be produced by multiple fungal species, including saprotrophs, plant, insect and human pathogens, particularly belonging to Fusarium, Beauveria and Isaria genera. The ability of Trichoderma species to produce BEA was until now uncertain. Biosynthesis of BEA is governed by a non-ribosomal peptide synthase (NRPS), known as beauvericin synthase (BEAS), which appears to present considerable divergence among different fungal species. In the present study we compared the production of beauvericin among Fusarium and Trichoderma strains using UPLC methods. BEAS fragments were sequenced and analyzed to examine the level of the gene's divergence between these two genera and confirm the presence of active BEAS copy in Trichoderma. Seventeen strains of twelve species were studied and phylogenetic analysis showed distinctive grouping of Fusarium and Trichoderma strains. The highest producers of beauvericin were F. proliferatum and F. nygamai. Trichoderma strains of three species (T. atroviride, T. viride, T. koningiopsis) were minor BEA producers. The study showed beauvericin production by Fusarium and Trichoderma species and high variance of the non-ribosomal peptide synthase gene among fungal species from the Hypocreales order.

Natural variation of root lesion nematode antagonism in the biocontrol fungus Clonostachys rosea and identification of biocontrol factors through genome-wide association mapping

Biological control is a promising approach to reduce plant diseases caused by nematodes to ensure high productivity in agricultural production. Large-scale analyses of genetic variation in fungal species used for biocontrol can generate knowledge regarding interaction mechanisms that can improve efficacy of biocontrol applications. In this study, we performed a genome-wide association study (GWAS) for in vitro antagonism against the root lesion nematode Pratylenchus penetrans in 53 previously genome re-sequenced strains of the biocontrol fungus Clonostachys rosea. Nematode mortality in C. rosea potato dextrose broth (PDB) culture filtrates was highly variable and showed continuous variation (p < .001) between strains, indicating a polygenic inheritance. Twenty-one strains produced culture filtrates with higher (p ≤ .05) nematode mortality compared with the PDB control treatment, while ten strains lowered (p ≤ .05) the mortality. The difference in in vitro antagonism against P. penetrans correlated with antagonism against the soybean cyst nematode Heterodera glycines, indicating lack of host specificity in C. rosea. An empirical Bayesian multiple hypothesis testing approach identified 279 single nucleotide polymorphism markers significantly (local false sign rate < 10-10) associated with the trait. Genes present in the genomic regions associated with nematicidal activity included several membrane transporters, a chitinase and genes encoding proteins predicted to biosynthesize secondary metabolites. Gene deletion strains of the predicted nonribosomal peptide synthetase genes nps4 and nps5 were generated and showed increased (p ≤ .001) fungal growth and conidiation rates compared to the wild type. Deletion strains also exhibited reduced (p < .001) nematicidal activity and reduced (p ≤ .05) biocontrol efficacy against nematode root disease and against fusarium foot rot on wheat. In summary, we show that the GWAS approach can be used to identify biocontrol factors in C. rosea, specifically the putative nonribosomal peptide synthetases NPS4 and NPS5.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus. Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Defects in the Ferroxidase That Participates in the Reductive Iron Assimilation System Results in Hypervirulence in Botrytis Cinerea

The plant-pathogenic fungus B. cinerea causes enormous economic losses, estimated at anywhere between $10 billion and $100 billion worldwide, under both pre- and postharvest conditions. Here, we present the characterization of a loss-of-function mutant in a component involved in iron acquisition that displays hypervirulence. While in different microbial systems iron uptake mechanisms appear to be critical to achieve full pathogenic potential, we found that the absence of the ferroxidase that is part of the reductive iron assimilation system leads to hypervirulence in this fungus. This is an unusual and rather underrepresented phenotype, which can be modulated by iron levels in the plant and provides an unexpected link between iron acquisition, reactive oxygen species (ROS) production, and pathogenesis in the Botrytis-plant interaction.

The plant pathogen Botrytis cinerea is responsible for gray-mold disease, which infects a wide variety of species. The outcome of this host-pathogen interaction, a result of the interplay between plant defense and fungal virulence pathways, can be modulated by various environmental factors. Among these, iron availability and acquisition play a crucial role in diverse biological functions. How B. cinerea obtains iron, an essential micronutrient, during infection is unknown. We set out to determine the role of the reductive iron assimilation (RIA) system during B. cinerea infection. This system comprises the BcFET1 ferroxidase, which belongs to the multicopper oxidase (MCO) family of proteins, and the BcFTR1 membrane-bound iron permease. Gene knockout and complementation studies revealed that, compared to the wild type, the bcfet1 mutant displays delayed conidiation, iron-dependent sclerotium production, and significantly reduced whole-cell iron content. Remarkably, this mutant exhibited a hypervirulence phenotype, whereas the bcftr1 mutant presents normal virulence and unaffected whole-cell iron levels and developmental programs. Interestingly, while in iron-starved plants wild-type B. cinerea produced slightly reduced necrotic lesions, the hypervirulence phenotype of the bcfet1 mutant is no longer observed in iron-deprived plants. This suggests that B. cinerea bcfet1 knockout mutants require plant-derived iron to achieve larger necrotic lesions, whereas in planta analyses of reactive oxygen species (ROS) revealed increased ROS levels only for infections caused by the bcfet1 mutant. These results suggest that increased ROS production, under an iron sufficiency environment, at least partly underlie the observed infection phenotype in this mutant.

Microorganisms populating the water-related indoor biome

Modernisation of our households created novel opportunities for microbial growth and thus changed the array of microorganisms we come in contact with. While many studies have investigated microorganisms in the air and dust, tap water, another major input of microbial propagules, has received far less attention. The quality of drinking water in developed world is strictly regulated to prevent immediate danger to human health. However, fungi, algae, protists and bacteria of less immediate concern are usually not screened for. These organisms can thus use water as a vector of transmission into the households, especially if they are resistant to various water treatment procedures. Good tolerance of unfavourable abiotic conditions is also important for survival once microbes enter the household. Limitation of water availability, high or low temperatures, application of antimicrobial chemicals and other measures are taken to prevent indoor microbial overgrowth. These conditions, together with a large number of novel chemicals in our homes, shape the diversity and abundance of indoor microbiota through constant selection of the most resilient species, resulting in a substantial overlap in diversity of indoor and natural extreme environments. At least in fungi, extremotolerance has been linked to human pathogenicity, explaining why many species found in novel indoor habitats (such as dishwasher) are notable opportunistic pathogens. As a result, microorganisms that often enter our households with water and are then enriched in novel indoor habitats might have a hitherto underestimated impact on the well-being of the increasingly indoor-bound human population. KEY POINTS: Domestic environment harbours a large diversity of microorganisms. Microbiota of water-related indoor habitats mainly originates from tap water. Bathrooms, kitchens and household appliances select for polyextremotolerant species. Many household-related microorganisms are human opportunistic pathogens.

Cyclodepsipeptide Biosynthesis in Hypocreales Fungi and Sequence Divergence of The

Fungi from the Hypocreales order synthesize a range of toxic non-ribosomal cyclic peptides with antimicrobial, insecticidal and cytotoxic activities. Entomopathogenic Beauveria, Isaria and Cordyceps as well as phytopathogenic Fusarium spp. are known producers of beauvericins (BEAs), beauvenniatins (BEAEs) or enniatins (ENNs). The compounds are synthesized by beauvericin/enniatin synthase (BEAS/ESYN1), which shows significant sequence divergence among Hypocreales members. We investigated ENN, BEA and BEAE production among entomopathogenic (Beauveria, Cordyceps, Isaria) and phytopathogenic (Fusarium) fungi; BEA and ENNs were quantified using an LC-MS/MS method. Phylogenetic analysis of partial sequences of putative BEAS/ESYN1 amplicons was also made. Nineteen fungal strains were identified based on sequence analysis of amplified ITS and tef-1α regions. BEA was produced by all investigated fungi, with F. proliferatum and F. concentricum being the most efficient producers. ENNs were synthesized mostly by F. acuminatum, F. avenaceum and C. confragosa. The phylogeny reconstruction suggests that ancestral BEA biosynthesis independently diverged into biosynthesis of other compounds. The divergent positioning of three Fusarium isolates raises the possibility of parallel acquisition of cyclic depsipeptide synthases in ancient complexes within Fusarium genus. Different fungi have independently evolved NRPS genes involved in depsipeptide biosynthesis, with functional adaptation towards biosynthesis of overlapping yet diversified metabolite profiles.

Changes in Peptaibol Production of Trichoderma Species during In Vitro Antagonistic Interactions with Fungal Plant Pathogens

Trichoderma species are widely used as biofungicides for the control of fungal plant pathogens. Several studies have been performed to identify the main genes and compounds involved in Trichoderma–plant–microbial pathogen cross-talks. However, there is not much information about the exact mechanism of this profitable interaction. Peptaibols secreted mainly by Trichoderma species are linear, 5–20 amino acid residue long, non-ribosomally synthesized peptides rich in α-amino isobutyric acid, which seem to be effective in Trichoderma–plant pathogenic fungus interactions. In the present study, reversed phase (RP) high-performance liquid chromatography (HPLC) coupled with electrospray ionization (ESI) mass spectrometry (MS) was used to detect peptaibol profiles of Trichoderma strains during interactions with fungal plant pathogens. MS investigations of the crude extracts deriving from in vitro confrontations of Trichoderma asperellum and T. longibrachiatum with different plant pathogenic fungi (Fusarium moniliforme, F. culmorum, F. graminearum, F. oxysporum species complex, Alternaria solani and Rhizoctonia solani) were performed to get a better insight into the role of these non-ribosomal antimicrobial peptides. The results revealed an increase in the total amount of peptaibols produced during the interactions, as well as some differences in the peptaibol profiles between the confrontational and control tests. Detection of the expression level of the peptaibol synthetase tex1 by qRT-PCR showed a significant increase in T. asperellum/R. solani interaction in comparison to the control. In conclusion, the interaction with plant pathogens highly influenced the peptaibol production of the examined Trichoderma strains.

Acrophiarin (antibiotic S31794/F-1) from Penicillium arenicola shares biosynthetic features with both Aspergillus- and Leotiomycete-type echinocandins

The antifungal echinocandin lipopeptide, acrophiarin, was circumscribed in a patent in 1979. We confirmed that the producing strain NRRL 8095 is Penicillium arenicola and other strains of P. arenicola produced acrophiarin and acrophiarin analogues. Genome sequencing of NRRL 8095 identified the acrophiarin gene cluster. Penicillium arenicola and echinocandin-producing Aspergillus species belong to the family Aspergillaceae of the Eurotiomycetes, but several features of acrophiarin and its gene cluster suggest a closer relationship with echinocandins from Leotiomycete fungi. These features include hydroxy-glutamine in the peptide core instead of a serine or threonine residue, the inclusion of a non-heme iron, α-ketoglutarate-dependent oxygenase for hydroxylation of the C3 of the glutamine, and a thioesterase. In addition, P. arenicola bears similarity to Leotiomycete echinocandin-producing species because it exhibits self-resistance to exogenous echinocandins. Phylogenetic analysis of the genes of the echinocandin biosynthetic family indicated that most of the predicted proteins of acrophiarin gene cluster exhibited higher similarity to the predicted proteins of the pneumocandin gene cluster of the Leotiomycete Glarea lozoyensis than to those of the echinocandin B gene cluster from A. pachycristatus. The fellutamide gene cluster and related gene clusters are recognized as relatives of the echinocandins. Inclusion of the acrophiarin gene cluster into a comprehensive phylogenetic analysis of echinocandin gene clusters indicated the divergent evolutionary lineages of echinocandin gene clusters are descendants from a common ancestral progenitor. The minimal 10-gene cluster may have undergone multiple gene acquisitions or losses and possibly horizontal gene transfer after the ancestral separation of the two lineages.

Acurin A, a novel hybrid compound, biosynthesized by individually translated PKS- and NRPS-encoding genes in Aspergillus aculeatus

This work presents the identification and proposed biosynthetic pathway for a compound of mixed polyketide-nonribosomal peptide origin that we named acurin A. The compound was isolated from an extract of the filamentous fungus Aspergillus aculeatus, and its core structure resemble that of the mycotoxin fusarin C produced by several Fusarium species. Based on bioinformatics in combination with RT-qPCR experiments and gene-deletion analysis, we identified a biosynthetic gene cluster (BGC) in A. aculeatus responsible for the biosynthesis of acurin A. Moreover, we were able to show that a polyketide synthase (PKS) and a nonribosomal peptide synthetase (NRPS) enzyme separately encoded by this BGC are responsible for the synthesis of the PK-NRP compound, acurin A, core structure. In comparison, the production of fusarin C is reported to be facilitated by a linked PKS-NRPS hybrid enzyme. Phylogenetic analyses suggest the PKS and NRPS in A. aculeatus resulted from a recent fission of an ancestral hybrid enzyme followed by gene duplication. In addition to the PKS- and NRPS-encoding genes of acurin A, we show that six other genes are influencing the biosynthesis including a regulatory transcription factor. Altogether, we have demonstrated the involvement of eight genes in the biosynthesis of acurin A, including an in-cluster transcription factor. This study highlights the biosynthetic capacity of A. aculeatus and serves as an example of how the CRISPR/Cas9 system can be exploited for the construction of fungal strains that can be readily engineered.

A detailed in silico analysis of secondary metabolite biosynthesis clusters in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

BACKGROUND

The broad host range pathogen Sclerotinia sclerotiorum infects over 400 plant species and causes substantial yield losses in crops worldwide. Secondary metabolites are known to play important roles in the virulence of plant pathogens, but little is known about the secondary metabolite repertoire of S. sclerotiorum. In this study, we predicted secondary metabolite biosynthetic gene clusters in the genome of S. sclerotiorum and analysed their expression during infection of Brassica napus using an existing transcriptome data set. We also investigated their sequence diversity among a panel of 25 previously published S. sclerotiorum isolate genomes.

RESULTS

We identified 80 putative secondary metabolite clusters. Over half of the clusters contained at least three transcriptionally coregulated genes. Comparative genomics revealed clusters homologous to clusters in the closely related plant pathogen Botrytis cinerea for production of carotenoids, hydroxamate siderophores, DHN melanin and botcinic acid. We also identified putative phytotoxin clusters that can potentially produce the polyketide sclerin and an epipolythiodioxopiperazine. Secondary metabolite clusters were enriched in subtelomeric genomic regions, and those containing paralogues showed a particularly strong association with repeats. The positional bias we identified was borne out by intraspecific comparisons that revealed putative secondary metabolite genes suffered more presence / absence polymorphisms and exhibited a significantly higher sequence diversity than other genes.

CONCLUSIONS

These data suggest that S. sclerotiorum produces numerous secondary metabolites during plant infection and that their gene clusters undergo enhanced rates of mutation, duplication and recombination in subtelomeric regions. The microevolutionary regimes leading to S. sclerotiorum secondary metabolite diversity have yet to be elucidated. Several potential phytotoxins documented in this study provide the basis for future functional analyses.