Computer-aided, resistance gene-guided genome mining for proteasome and HMG-CoA reductase inhibitors

Secondary metabolites (SMs) are biologically active small molecules, many of which are medically valuable. Fungal genomes contain vast numbers of SM biosynthetic gene clusters (BGCs) with unknown products, suggesting that huge numbers of valuable SMs remain to be discovered. It is challenging, however, to identify SM BGCs, among the millions present in fungi, that produce useful compounds. One solution is resistance gene-guided genome mining, which takes advantage of the fact that some BGCs contain a gene encoding a resistant version of the protein targeted by the compound produced by the BGC. The bioinformatic signature of such BGCs is that they contain an allele of an essential gene with no SM biosynthetic function, and there is a second allele elsewhere in the genome. We have developed a computer-assisted approach to resistance gene-guided genome mining that allows users to query large databases for BGCs that putatively make compounds that have targets of therapeutic interest. Working with the MycoCosm genome database, we have applied this approach to look for SM BGCs that target the proteasome β6 subunit, the target of the proteasome inhibitor fellutamide B, or HMG-CoA reductase, the target of cholesterol reducing therapeutics such as lovastatin. Our approach proved effective, finding known fellutamide and lovastatin BGCs as well as fellutamide- and lovastatin-related BGCs with variations in the SM genes that suggest they may produce structural variants of fellutamides and lovastatin. Gratifyingly, we also found BGCs that are not closely related to lovastatin BGCs but putatively produce novel HMG-CoA reductase inhibitors.

ONE-SENTENCE SUMMARY

A new computer-assisted approach to resistance gene-directed genome mining is reported along with its use to identify fungal biosynthetic gene clusters that putatively produce proteasome and HMG-CoA reductase inhibitors.

Genetic co-option into plant-filamentous pathogen interactions

Plants are engaged in a coevolutionary arms race with their pathogens that drives rapid diversification and specialization of genes involved in resistance and virulence. However, some major innovations in plant-pathogen interactions, such as molecular decoys, trans-kingdom RNA interference, two-speed genomes, and receptor networks, evolved through the expansion of the functional landscape of genes. This is a typical outcome of genetic co-option, the evolutionary process by which available genes are recruited into new biological functions. Co-option into plant-pathogen interactions emerges generally from (i) cis-regulatory variation, (ii) horizontal gene transfer (HGT), (iii) mutations altering molecular promiscuity, and (iv) rewiring of gene networks and protein complexes. Understanding these molecular mechanisms is key for the functional and predictive biology of plant-pathogen interactions.

Genome-Based Analysis of Verticillium Polyketide Synthase Gene Clusters

Simple Summary

Fungi can produce many types of secondary metabolites, including mycotoxins. Poisonous mushrooms and mycotoxins that cause food spoilage have been known for a very long time. For example, Aspergillus flavus, which can grow on grains and nuts, produces highly toxic substances called Aflatoxins. Despite their menace to other living organisms, mycotoxins can be used for medicinal purposes, i.e., as antibiotics, growth-promoting compounds, and other kinds of drugs. These and other secondary metabolites produced by plant-pathogenic fungi may cause host plants to display disease symptoms and may play a substantial role in disease progression. Therefore, the identification and characterization of the genes involved in their biosynthesis are essential for understanding the molecular mechanism involved in their biosynthetic pathways and further promoting sustainable knowledge-based crop production.

Abstract

Polyketides are structurally diverse and physiologically active secondary metabolites produced by many organisms, including fungi. The biosynthesis of polyketides from acyl-CoA thioesters is catalyzed by polyketide synthases, PKSs. Polyketides play roles including in cell protection against oxidative stress, non-constitutive (toxic) roles in cell membranes, and promoting the survival of the host organisms. The genus Verticillium comprises many species that affect a wide range of organisms including plants, insects, and other fungi. Many are known as causal agents of Verticillium wilt diseases in plants. In this study, a comparative genomics approach involving several Verticillium species led us to evaluate the potential of Verticillium species for producing polyketides and to identify putative polyketide biosynthesis gene clusters. The next step was to characterize them and predict the types of polyketide compounds they might produce. We used publicly available sequences from ten species of Verticillium including V. dahliae, V. longisporum, V. nonalfalfae, V. alfalfae, V. nubilum, V. zaregamsianum, V. klebahnii, V. tricorpus, V. isaacii, and V. albo-atrum to identify and characterize PKS gene clusters by utilizing a range of bioinformatic and phylogenetic approaches. We found 32 putative PKS genes and possible clusters in the genomes of Verticillium species. All the clusters appear to be complete and functional. In addition, at least five clusters including putative DHN-melanin-, cytochalasin-, fusarielien-, fujikurin-, and lijiquinone-like compounds may belong to the active PKS repertoire of Verticillium. These results will pave the way for further functional studies to understand the role of these clusters.

Comparative Genomics of Fusarium circinatum Isolates Used to Screen Southern Pines for Pitch Canker Resistance

Pitch canker, caused by the fungal pathogen Fusarium circinatum, is a global disease affecting many Pinus spp. Often fatal, this disease causes significant mortality in both commercially grown and natural pine forests and is an issue of current and growing concern. F. circinatum isolates collected from three locations in the U.S. state of Florida were shown to be virulent on both slash and loblolly pine, with two of the isolates causing equivalent and significantly larger lesions than those caused by the third isolate during pathogenicity trials. In addition, significant genetic variation in lesion length in the pedigreed slash pine population was evident and rankings of parents for lesion length were similar across isolates. Experimental data demonstrate that both host and pathogen genetics contribute to disease severity. High-quality genomic assemblies of all three isolates were created and compared for structural differences and gene content. No major structural differences were observed among the isolates; however, missing or altered genes do contribute to genomic variation in the pathogen population. This work evaluates in planta virulence among three isolates of F. circinatum, provides genomic resources to facilitate study of this organism, and details comparative genomic methods that may be used to explore the pathogen's contribution to disease development.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Secondary Metabolite Gene Regulation in Mycotoxigenic Fusarium Species: A Focus on Chromatin

Fusarium is a species-rich group of mycotoxigenic plant pathogens that ranks as one of the most economically important fungal genera in the world. During growth and infection, they are able to produce a vast spectrum of low-molecular-weight compounds, so-called secondary metabolites (SMs). SMs often comprise toxic compounds (i.e., mycotoxins) that contaminate precious food and feed sources and cause adverse health effects in humans and livestock. In this context, understanding the regulation of their biosynthesis is crucial for the development of cropping strategies that aim at minimizing mycotoxin contamination in the field. Nevertheless, currently, only a fraction of SMs have been identified, and even fewer are considered for regular monitoring by regulatory authorities. Limitations to exploit their full chemical potential arise from the fact that the genes involved in their biosynthesis are often silent under standard laboratory conditions and only induced upon specific stimuli mimicking natural conditions in which biosynthesis of the respective SM becomes advantageous for the producer. This implies a complex regulatory network. Several components of these gene networks have been studied in the past, thereby greatly advancing the understanding of SM gene regulation and mycotoxin biosynthesis in general. This review aims at summarizing the latest advances in SM research in these notorious plant pathogens with a focus on chromatin structure.

Potentials of Endophytic Fungi in the Biosynthesis of Versatile Secondary Metabolites and Enzymes

World population growth and modernization have engendered multiple environmental problems: the propagation of humans and crop diseases and the development of multi-drug-resistant fungi, bacteria and viruses. Thus, a considerable shift towards eco-friendly products has been seen in medicine, pharmacy, agriculture and several other vital sectors. Nowadays, studies on endophytic fungi and their biotechnological potentials are in high demand due to their substantial, cost-effective and eco-friendly contributions in the discovery of an array of secondary metabolites. For this review, we provide a brief overview of plant–endophytic fungi interactions and we also state the history of the discovery of the untapped potentialities of fungal secondary metabolites. Then, we highlight the huge importance of the discovered metabolites and their versatile applications in several vital fields including medicine, pharmacy, agriculture, industry and bioremediation. We then focus on the challenges and on the possible methods and techniques that can be used to help in the discovery of novel secondary metabolites. The latter range from endophytic selection and culture media optimization to more in-depth strategies such as omics, ribosome engineering and epigenetic remodeling.

Fusarium chaquense, sp. nov, a novel type A trichothecene-producing species from native grasses in a wetland ecosystem in Argentina

The Chaco wetland is among the most biologically diverse regions in Argentina. In collections of fungi from asymptomatic native grasses (Poaceae) from the wetlands, we identified isolates of Fusarium that were morphologically similar to F. armeniacum, but distinct from it by their production of abundant microconidia. All the isolates had identical, or nearly identical, partial sequences of TEF1 and RPB2. But they were distinct from reference sequences from F. armeniacum and Fusarium species closely related to it. Phylogenetic analysis of 34 full-length housekeeping gene sequences retrieved from whole genome sequences of three Chaco wetland isolates, 29 genes resolved the isolates as an exclusive clade within the F. sambucinum species complex. Based on results of the morphological and phylogenetic analysis, we concluded that the Chaco wetland isolates are a distinct and novel species, herein described as Fusarium chaquense, sp. nov., which is closely related to F. armeniacum. F. chaquense in culture can produce the trichothecenes T-2 and HT-2 toxin, neosolaniol, diacetoxyscirpenol, and monoacetoxyscirpenol, as well as beauvericin and the pigment aurofusarin. Genome sequence analysis also revealed the presence of three previously described loci required for trichothecene biosynthesis. This research represents the first study of Fusarium in a natural ecosystem in Argentina.

Bioprospecting Trichoderma: A Systematic Roadmap to Screen Genomes and Natural Products for Biocontrol Applications

Natural products derived from microbes are crucial innovations that would help in reaching sustainability development goals worldwide while achieving bioeconomic growth. Trichoderma species are well-studied model fungal organisms used for their biocontrol properties with great potential to alleviate the use of agrochemicals in agriculture. However, identifying and characterizing effective natural products in novel species or strains as biological control products remains a meticulous process with many known challenges to be navigated. Integration of recent advancements in various "omics" technologies, next generation biodesign, machine learning, and artificial intelligence approaches could greatly advance bioprospecting goals. Herein, we propose a roadmap for assessing the potential impact of already known or newly discovered Trichoderma species for biocontrol applications. By screening publicly available Trichoderma genome sequences, we first highlight the prevalence of putative biosynthetic gene clusters and antimicrobial peptides among genomes as an initial step toward predicting which organisms could increase the diversity of natural products. Next, we discuss high-throughput methods for screening organisms to discover and characterize natural products and how these findings impact both fundamental and applied research fields.

Comparative analyses of the Hymenoscyphus fraxineus and Hymenoscyphus albidus genomes reveals potentially adaptive differences in secondary metabolite and transposable element repertoires

Background

The dieback epidemic decimating common ash (Fraxinus excelsior) in Europe is caused by the invasive fungus Hymenoscyphus fraxineus. In this study we analyzed the genomes of H. fraxineus and H. albidus, its native but, now essentially displaced, non-pathogenic sister species, and compared them with several other members of Helotiales. The focus of the analyses was to identify signals in the genome that may explain the rapid establishment of H. fraxineus and displacement of H. albidus.

Results

The genomes of H. fraxineus and H. albidus showed a high level of synteny and identity. The assembly of H. fraxineus is 13 Mb longer than that of H. albidus’, most of this difference can be attributed to higher dispersed repeat content (i.e. transposable elements [TEs]) in H. fraxineus. In general, TE families in H. fraxineus showed more signals of repeat-induced point mutations (RIP) than in H. albidus, especially in Long-terminal repeat (LTR)/Copia and LTR/Gypsy elements.

Comparing gene family expansions and 1:1 orthologs, relatively few genes show signs of positive selection between species. However, several of those did appeared to be associated with secondary metabolite genes families, including gene families containing two of the genes in the H. fraxineus-specific, hymenosetin biosynthetic gene cluster (BGC).

Conclusion

The genomes of H. fraxineus and H. albidus show a high degree of synteny, and are rich in both TEs and BGCs, but the genomic signatures also indicated that H. albidus may be less well equipped to adapt and maintain its ecological niche in a rapidly changing environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07837-2.

Mind the mushroom: natural product biosynthetic genes and enzymes of Basidiomycota

Covering: up to September 2020 Mushroom-forming fungi of the division Basidiomycota have traditionally been recognised as prolific producers of structurally diverse and often bioactive secondary metabolites, using the methods of chemistry for research. Over the past decade, -omics technologies were applied on these fungi, and sophisticated heterologous gene expression platforms emerged, which have boosted research into the genetic and biochemical basis of the biosyntheses. This review provides an overview on experimentally confirmed natural product biosyntheses of basidiomycete polyketides, amino acid-derived products, terpenoids, and volatiles. We also present challenges and solutions particular to natural product research with these fungi. 222 references are cited.

Origin and Evolution of Fusidane-Type Antibiotics Biosynthetic Pathway through Multiple Horizontal Gene Transfers

Fusidane-type antibiotics represented by fusidic acid, helvolic acid, and cephalosporin P1 have very similar core structures, but they are produced by fungi belonging to different taxonomic groups. The origin and evolution of fusidane-type antibiotics biosynthetic gene clusters (BGCs) in different antibiotics producing strains remained an enigma. In this study, we investigated the distribution and evolution of the fusidane BGCs in 1,284 fungal genomes. We identified 12 helvolic acid BGCs, 4 fusidic acid BGCs, and 1 cephalosporin P1 BGC in Pezizomycotina fungi. Phylogenetic analyses indicated six horizontal gene transfer (HGT) events in the evolutionary trajectory of the BGCs, including 1) three transfers across Eurotiomycetes and Sordariomycetes classes; 2) one transfer between genera under Sordariomycetes class; and 3) two transfers within Aspergillus genus under Eurotiomycetes classes. Finally, we proposed that the ancestor of fusidane BGCs would be originated from the Zoopagomycota by ancient HGT events according to the phylogenetic trees of key enzymes in fusidane BGCs (OSC and P450 genes). Our results extensively clarify the evolutionary trajectory of fusidane BGCs by HGT among distantly related fungi and provide new insights into the evolutionary mechanisms of metabolic pathways in fungi.

Biosynthetic gene clusters and the evolution of fungal chemodiversity

Covering: up to 2019Fungi produce a remarkable diversity of secondary metabolites: small, bioactive molecules not required for growth but which are essential to their ecological interactions with other organisms. Genes that participate in the same secondary metabolic pathway typically reside next to each other in fungal genomes and form biosynthetic gene clusters (BGCs). By synthesizing state-of-the-art knowledge on the evolution of BGCs in fungi, we propose that fungal chemodiversity stems from three molecular evolutionary processes involving BGCs: functional divergence, horizontal transfer, and de novo assembly. We provide examples of how these processes have contributed to the generation of fungal chemodiversity, discuss their relative importance, and outline major, outstanding questions in the field.

An endophyte of Macrochloa tenacissima (esparto or needle grass) from Tunisia is a novel species in the Fusarium redolens species complex

Here, we report on the morphological, molecular, and chemical characterization of a novel Fusarium species recovered from the roots and rhizosphere of Macrochloa tenacissima (halfa, esparto, or needle grass) in central Tunisia. Formally described here as F. spartum, this species is a member of the Fusarium redolens species complex but differs from the other two species within the complex, F. redolens and F. hostae, by its endophytic association with M. tenacissima and its genealogical exclusivity based on multilocus phylogenetic analyses. To assess their sexual reproductive mode, a polymerase chain reaction (PCR) assay was designed and used to screen the three strains of F. spartum, 51 of F. redolens, and 14 of F. hostae for mating type (MAT) idiomorph. Genetic architecture of the MAT locus in the former two species suggests that if they reproduce sexually, it is via obligate outcrossing. By comparison, results of the PCR assay indicated that 13/14 of the F. hostae strains possessed MAT1-1 and MAT1-2 idiomorphs and thus might be self-fertile or homothallic. However, when the F. hostae strains were selfed, 11 failed to produce perithecia and one only produced several small abortive perithecia. Cirrhi with ascospores, however, were only produced by 8/28 and 4/84 of the variable size perithecia, respectively, of F. hostae NRRL 29888 and 29890. The potential for the three F. redolens clade species to produce mycotoxins, pigments, and phytohormones was assessed by screening whole genome sequence data and by analyzing extracts on cracked maize kernel cultures via liquid chromatography-mass spectrometry.

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.

Rapid evolution in plant-microbe interactions - a molecular genomics perspective

Rapid (co-)evolution at multiple timescales is a hallmark of plant-microbe interactions. The mechanistic basis for the rapid evolution largely rests on the features of the genomes of the interacting partners involved. Here, we review recent insights into genomic characteristics and mechanisms that enable rapid evolution of both plants and phytopathogens. These comprise fresh insights in allelic series of matching pairs of resistance and avirulence genes, the generation of novel pathogen effectors, the recently recognised small RNA warfare, and genomic aspects of secondary metabolite biosynthesis. In addition, we discuss the putative contributions of permissive host environments, transcriptional plasticity and the role of ploidy on the interactions. We conclude that the means underlying the rapid evolution of plant-microbe interactions are multifaceted and depend on the particular nature of each interaction.

Gain and loss of a transcription factor that regulates late trichothecene biosynthetic pathway genes in Fusarium

Trichothecenes are among the mycotoxins of most concern to food and feed safety and are produced by species in two lineages of Fusarium: the F. incarnatum-equiseti (FIESC) and F. sambucinum (FSAMSC) species complexes. Previous functional analyses of the trichothecene biosynthetic gene (TRI) cluster in members of FSAMSC indicate that the transcription factor gene TRI6 activates expression of other TRI cluster genes. In addition, previous sequence analyses indicate that the FIESC TRI cluster includes TRI6 and another uncharacterized transcription factor gene (hereafter TRI21) that was not reported in FSAMSC. Here, gene deletion analysisindicated that in FIESC TRI6 functions in a manner similar to FSAMSC, whereas TRI21 activated expression of some genes that function late in the trichothecene biosynthetic pathway but not early-pathway genes. Consistent with this finding, TRI21 was required for formation of diacetoxyscripenol, a late-trichothecene-pathway product, but not for isotrichodermin, an early-pathway product. Although intact homologs of TRI21 were not detected in FSAMSC or other trichothecene-producing fungal genera, TRI21 fragments were detected in some FSAMSC species. This suggests that the gene was acquired by Fusarium after divergence from other trichothecene-producing fungi, was subsequently lost in FSAMSC, but was retained in FIESC. Together, our results indicate fundamental differences in regulation of trichothecene biosynthesis in FIESC and FSAMSC.

Complex Evolutionary Origins of Specialized Metabolite Gene Cluster Diversity among the Plant Pathogenic Fungi of the Fusarium graminearum Species Complex

Fungal genomes encode highly organized gene clusters that underlie the production of specialized (or secondary) metabolites. Gene clusters encode key functions to exploit plant hosts or environmental niches. Promiscuous exchange among species and frequent reconfigurations make gene clusters some of the most dynamic elements of fungal genomes. Despite evidence for high diversity in gene cluster content among closely related strains, the microevolutionary processes driving gene cluster gain, loss, and neofunctionalization are largely unknown. We analyzed the Fusarium graminearum species complex (FGSC) composed of plant pathogens producing potent mycotoxins and causing Fusarium head blight on cereals. We de novo assembled genomes of previously uncharacterized FGSC members (two strains of F. austroamericanum, F. cortaderiae, and F. meridionale). Our analyses of 8 species of the FGSC in addition to 15 other Fusarium species identified a pangenome of 54 gene clusters within FGSC. We found that multiple independent losses were a key factor generating extant cluster diversity within the FGSC and the Fusarium genus. We identified a modular gene cluster conserved among distantly related fungi, which was likely reconfigured to encode different functions. We also found strong evidence that a rare cluster in FGSC was gained through an ancient horizontal transfer between bacteria and fungi. Chromosomal rearrangements underlying cluster loss were often complex and were likely facilitated by an enrichment in specific transposable elements. Our findings identify important transitory stages in the birth and death process of specialized metabolism gene clusters among very closely related species.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

The birth, evolution and death of metabolic gene clusters in fungi

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

Fungal secondary metabolism: regulation, function and drug discovery

One of the exciting movements in microbial sciences has been a refocusing and revitalization of efforts to mine the fungal secondary metabolome. The magnitude of biosynthetic gene clusters (BGCs) in a single filamentous fungal genome combined with the historic number of sequenced genomes suggests that the secondary metabolite wealth of filamentous fungi is largely untapped. Mining algorithms and scalable expression platforms have greatly expanded access to the chemical repertoire of fungal-derived secondary metabolites. In this Review, I discuss new insights into the transcriptional and epigenetic regulation of BGCs and the ecological roles of fungal secondary metabolites in warfare, defence and development. I also explore avenues for the identification of new fungal metabolites and the challenges in harvesting fungal-derived secondary metabolites.

Chromosome rearrangements shape the diversification of secondary metabolism in the cyclosporin producing fungus Tolypocladium inflatum

BACKGROUND

Genes involved in production of secondary metabolites (SMs) in fungi are exceptionally diverse. Even strains of the same species may exhibit differences in metabolite production, a finding that has important implications for drug discovery. Unlike in other eukaryotes, genes producing SMs are often clustered and co-expressed in fungal genomes, but the genetic mechanisms involved in the creation and maintenance of these secondary metabolite biosynthetic gene clusters (SMBGCs) remains poorly understood.

RESULTS

In order to address the role of genome architecture and chromosome scale structural variation in generating diversity of SMBGCs, we generated chromosome scale assemblies of six geographically diverse isolates of the insect pathogenic fungus Tolypocladium inflatum, producer of the multi-billion dollar lifesaving immunosuppressant drug cyclosporin, and utilized a Hi-C chromosome conformation capture approach to address the role of genome architecture and structural variation in generating intraspecific diversity in SMBGCs. Our results demonstrate that the exchange of DNA between heterologous chromosomes plays an important role in generating novelty in SMBGCs in fungi. In particular, we demonstrate movement of a polyketide synthase (PKS) and several adjacent genes by translocation to a new chromosome and genomic context, potentially generating a novel PKS cluster. We also provide evidence for inter-chromosomal recombination between nonribosomal peptide synthetases located within subtelomeres and uncover a polymorphic cluster present in only two strains that is closely related to the cluster responsible for biosynthesis of the mycotoxin aflatoxin (AF), a highly carcinogenic compound that is a major public health concern worldwide. In contrast, the cyclosporin cluster, located internally on chromosomes, was conserved across strains, suggesting selective maintenance of this important virulence factor for infection of insects.

CONCLUSIONS

This research places the evolution of SMBGCs within the context of whole genome evolution and suggests a role for recombination between chromosomes in generating novel SMBGCs in the medicinal fungus Tolypocladium inflatum.

Functional analysis of polyketide synthase genes in the biocontrol fungus Clonostachys rosea

Clonostachys rosea is a mycoparasitic fungus used for biological control of plant diseases. Its genome contains 31 genes putatively encoding for polyketide synthases (PKSs), 75% of which are arranged in biosynthetic gene clusters. Gene expression analysis during C. rosea interactions with the fungal plant pathogens Botrytis cinerea and Fusarium graminearum showed common and species-specific induction of PKS genes. Our data showed a culture media dependent correlation between PKS gene expression and degree of antagonism in C. rosea. The pks22 and pks29 genes were highly induced during fungal-fungal interactions but not during pigmentation, and gene deletion studies revealed that PKS29 was required for full antagonism against B. cinerea, and for biocontrol of fusarium foot rot on barley. Metabolite analysis revealed that Δpks29 strains has a 50% reduced production (P = 0.001) of an unknown polyketide with molecular formula C₁₅H₂₈O₃, while Δpks22 strains lost the ability to produce four previously unknown polyketides named Clonorosein A-D. Clonorosein A and B were purified, their structures determined, and showed strong antifungal activity against B. cinerea and F. graminearum. These results show that PKS22 is required for production of antifungal polyketide Clonorosein A-D, and demonstrate the role of PKS29 in antagonism and biocontrol of fungal plant diseases.

Metabolic Gene Clusters in Eukaryotes

In bacteria, more than half of the genes in the genome are organized in operons. In contrast, in eukaryotes, functionally related genes are usually dispersed across the genome. There are, however, numerous examples of functional clusters of nonhomologous genes for metabolic pathways in fungi and plants. Despite superficial similarities with operons (physical clustering, coordinate regulation), these clusters have not usually originated by horizontal gene transfer from bacteria, and (unlike operons) the genes are typically transcribed separately rather than as a single polycistronic message. This clustering phenomenon raises intriguing questions about the origins of clustered metabolic pathways in eukaryotes and the significance of clustering for pathway function. Here we review metabolic gene clusters from fungi and plants, highlight commonalities and differences, and consider how these clusters form and are regulated. We also identify opportunities for future research in the areas of large-scale genomics, synthetic biology, and experimental evolution.

The MSDIN family in amanitin-producing mushrooms and evolution of the prolyl oligopeptidase genes

The biosynthetic pathway for amanitins and related cyclic peptides in deadly Amanita (Amanitaceae) mushrooms represents the first known ribosomal cyclic peptide pathway in the Fungi. Amanitins are found outside of the genus in distantly related agarics Galerina (Strophariaceae) and Lepiota (Agaricaceae). A long-standing question in the field persists: why is this pathway present in these phylogenetically disjunct agarics? Two deadly mushrooms, A. pallidorosea and A. subjunquillea, were deep sequenced, and sequences of biosynthetic genes encoding MSDINs (cyclic peptide precursor) and prolyl oligopeptidases (POPA and POPB) were obtained. The two Amanita species yielded 29 and 18 MSDINs, respectively. In addition, two MSDIN sequences were cloned from L. brunneoincarnata basidiomes. The toxin MSDIN genes encoding amatoxins or phallotoxins from the three genera were compared, and a phylogenetic tree constructed. Prolyl oligopeptidase B (POPB), a key enzyme in the biosynthetic pathway, was used in phylogenetic reconstruction to infer the evolutionary history of the genes. Phylogenies of POPB and POPA based on both coding and amino acid sequences showed very different results: while POPA genes clearly reflected the phylogeny of the host species, POPB did not; strikingly, it formed a well-supported monophyletic clade, despite that the species belong to different genera in disjunct families. POPA, a known house-keeping gene, was shown to be restricted in a branch containing only Amanita species and the phylogeny resembled that of those Amanita species. Phylogenetic analyses of MSDIN and POPB genes showed tight coordination and disjunct distribution. A POPB gene tree was compared with a corresponding species tree, and distances and substitution rates were compared. The result suggested POPB genes have significant smaller distances and rates than the house-keeping rpb2, discounting massive gene loss. Under this assumption, the incongruency between the gene tree and species tree was shown with strong support. Additionally, k-mer analyses consistently cluster Galerina and Amanita POPB genes, while Lepiota POPB is distinct. Our result suggests that horizontal gene transfer (HGT), at least between Amanita and Galerina, was involved in the acquisition of POPB genes, which may shed light on the evolution of the α-amanitin biosynthetic pathway.

Evolution and Diversity of Biosynthetic Gene Clusters in Fusarium

Plant pathogenic fungi in the Fusarium genus cause severe damage to crops, resulting in great financial losses and health hazards. Specialized metabolites synthesized by these fungi are known to play key roles in the infection process, and to provide survival advantages inside and outside the host. However, systematic studies of the evolution of specialized metabolite-coding potential across Fusarium have been scarce. Here, we apply a combination of bioinformatic approaches to identify biosynthetic gene clusters (BGCs) across publicly available genomes from Fusarium, to group them into annotated families and to study gain/loss events of BGC families throughout the history of the genus. Comparison with MIBiG reference BGCs allowed assignment of 29 gene cluster families (GCFs) to pathways responsible for the production of known compounds, while for 57 GCFs, the molecular products remain unknown. Comparative analysis of BGC repertoires using ancestral state reconstruction raised several new hypotheses on how BGCs contribute to Fusarium pathogenicity or host specificity, sometimes surprisingly so: for example, a gene cluster for the biosynthesis of hexadehydro-astechrome was identified in the genome of the biocontrol strain Fusarium oxysporum Fo47, while being absent in that of the tomato pathogen F. oxysporum f.sp. lycopersici. Several BGCs were also identified on supernumerary chromosomes; heterologous expression of genes for three terpene synthases encoded on the Fusarium poae supernumerary chromosome and subsequent GC/MS analysis showed that these genes are functional and encode enzymes that each are able to synthesize koraiol; this observed functional redundancy supports the hypothesis that localization of copies of BGCs on supernumerary chromosomes provides freedom for evolutionary innovations to occur, while the original function remains conserved. Altogether, this systematic overview of biosynthetic diversity in Fusarium paves the way for targeted natural product discovery based on automated identification of species-specific pathways as well as for connecting species ecology to the taxonomic distributions of BGCs.

Evolution of structural diversity of trichothecenes, a family of toxins produced by plant pathogenic and entomopathogenic fungi

Trichothecenes are a family of terpenoid toxins produced by multiple genera of fungi, including plant and insect pathogens. Some trichothecenes produced by the fungus Fusarium are among the mycotoxins of greatest concern to food and feed safety because of their toxicity and frequent occurrence in cereal crops, and trichothecene production contributes to pathogenesis of some Fusarium species on plants. Collectively, fungi produce over 150 trichothecene analogs: i.e., molecules that share the same core structure but differ in patterns of substituents attached to the core structure. Here, we carried out genomic, phylogenetic, gene-function, and analytical chemistry studies of strains from nine fungal genera to identify genetic variation responsible for trichothecene structural diversity and to gain insight into evolutionary processes that have contributed to the variation. The results indicate that structural diversity has resulted from gain, loss, and functional changes of trichothecene biosynthetic (TRI) genes. The results also indicate that the presence of some substituents has arisen independently in different fungi by gain of different genes with the same function. Variation in TRI gene duplication and number of TRI loci was also observed among the fungi examined, but there was no evidence that such genetic differences have contributed to trichothecene structural variation. We also inferred ancestral states of the TRI cluster and trichothecene biosynthetic pathway, and proposed scenarios for changes in trichothecene structures during divergence of TRI cluster homologs. Together, our findings provide insight into evolutionary processes responsible for structural diversification of toxins produced by pathogenic fungi.

Toxins produced by pathogens can contribute to infection and/or colonization of hosts. Some toxins consist of a family of metabolites with similar but distinct chemical structures. This structural variation can affect biological activity, which in turn likely contributes to adaptation to different environments, including to different hosts. Trichothecene toxins consist of over 150 structurally distinct molecules produced by certain fungi, including some plant and insect pathogens. In multiple systems that have been examined, trichothecenes contribute to pathogenesis on plants. To elucidate the evolutionary processes that have given rise to trichothecene structural variation, we conducted comparative analyses of nine fungal genera, most of which produce different trichothecene structures. Using genomic, molecular biology, phylogenetic, and analytical chemistry approaches, we obtained evidence that trichothecene structural variation has arisen primarily from gain, loss, and functional changes of trichothecene biosynthetic genes. Our results also indicate that some structural changes have arisen independently in different fungi. Our findings provide insight into genetic and biochemical changes that can occur in toxin biosynthetic pathways as fungi with the pathways adapt to different environmental conditions.

Horizontal gene cluster transfer increased hallucinogenic mushroom diversity

Secondary metabolites are a heterogeneous class of chemicals that often mediate interactions between species. The tryptophan‐derived secondary metabolite, psilocin, is a serotonin receptor agonist that induces altered states of consciousness. A phylogenetically disjunct group of mushroom‐forming fungi in the Agaricales produce the psilocin prodrug, psilocybin. Spotty phylogenetic distributions of fungal compounds are sometimes explained by horizontal transfer of metabolic gene clusters among unrelated fungi with overlapping niches. We report the discovery of a psilocybin gene cluster in three hallucinogenic mushroom genomes, and evidence for its horizontal transfer between fungal lineages. Patterns of gene distribution and transmission suggest that synthesis of psilocybin may have provided a fitness advantage in the dung and late wood‐decay fungal niches, which may serve as reservoirs of fungal indole‐based metabolites that alter behavior of mycophagous and wood‐eating invertebrates. These hallucinogenic mushroom genomes will serve as models in neurochemical ecology, advancing the (bio)prospecting and synthetic biology of novel neuropharmaceuticals.

Drivers of genetic diversity in secondary metabolic gene clusters within a fungal species

Filamentous fungi produce a diverse array of secondary metabolites (SMs) critical for defense, virulence, and communication. The metabolic pathways that produce SMs are found in contiguous gene clusters in fungal genomes, an atypical arrangement for metabolic pathways in other eukaryotes. Comparative studies of filamentous fungal species have shown that SM gene clusters are often either highly divergent or uniquely present in one or a handful of species, hampering efforts to determine the genetic basis and evolutionary drivers of SM gene cluster divergence. Here, we examined SM variation in 66 cosmopolitan strains of a single species, the opportunistic human pathogen Aspergillus fumigatus. Investigation of genome-wide within-species variation revealed 5 general types of variation in SM gene clusters: nonfunctional gene polymorphisms; gene gain and loss polymorphisms; whole cluster gain and loss polymorphisms; allelic polymorphisms, in which different alleles corresponded to distinct, nonhomologous clusters; and location polymorphisms, in which a cluster was found to differ in its genomic location across strains. These polymorphisms affect the function of representative A. fumigatus SM gene clusters, such as those involved in the production of gliotoxin, fumigaclavine, and helvolic acid as well as the function of clusters with undefined products. In addition to enabling the identification of polymorphisms, the detection of which requires extensive genome-wide synteny conservation (e.g., mobile gene clusters and nonhomologous cluster alleles), our approach also implicated multiple underlying genetic drivers, including point mutations, recombination, and genomic deletion and insertion events as well as horizontal gene transfer from distant fungi. Finally, most of the variants that we uncover within A. fumigatus have been previously hypothesized to contribute to SM gene cluster diversity across entire fungal classes and phyla. We suggest that the drivers of genetic diversity operating within a fungal species shown here are sufficient to explain SM cluster macroevolutionary patterns.

All organisms produce metabolites, which are small molecules important for growth, reproduction, and other essential functions. Some organisms, including fungi, plants, and bacteria, make specialized forms of metabolites known as “secondary” metabolites that are ecologically important and improve their producers’ chances of survival and reproduction. In fungi, the genes in pathways that synthesize secondary metabolites are typically located next to each other in the genome and organized in contiguous gene clusters. These gene clusters, along with the metabolites they produce, are highly distinct, even between otherwise similar fungi, and it is often difficult to reconstruct how these differences evolved. To understand how secondary metabolic pathways evolve in fungi, we compared secondary metabolic gene clusters in 66 strains of one species of filamentous fungus, the human pathogen Aspergillus fumigatus. We show that these gene clusters vary extensively within this species, and describe the genetic processes that cause these differences. We identify 5 types of variants: single nucleotide changes, gene and gene cluster gain and loss, different gene clusters at the same genomic position, and mobile gene clusters that “jump” around the genome. These results provide a road map to the types and frequencies of genomic changes underlying the extensive diversity of fungal secondary metabolites.

Fungal Gene Cluster Diversity and Evolution

Metabolic gene clusters (MGCs) have provided some of the earliest glimpses at the biochemical machinery of yeast and filamentous fungi. MGCs encode diverse genetic mechanisms for nutrient acquisition and the synthesis/degradation of essential and adaptive metabolites. Beyond encoding the enzymes performing these discrete anabolic or catabolic processes, MGCs may encode a range of mechanisms that enable their persistence as genetic consortia; these include enzymatic mechanisms to protect their host fungi from their inherent toxicities, and integrated regulatory machinery. This modular, self-contained nature of MGCs contributes to the metabolic and ecological adaptability of fungi. The phylogenetic and ecological patterns of MGC distribution reflect the broad diversity of fungal life cycles and nutritional modes. While the origins of most gene clusters are enigmatic, MGCs are thought to be born into a genome through gene duplication, relocation, or horizontal transfer, and analyzing the death and decay of gene clusters provides clues about the mechanisms selecting for their assembly. Gene clustering may provide inherent fitness advantages through metabolic efficiency and specialization, but experimental evidence for this is currently limited. The identification and characterization of gene clusters will continue to be powerful tools for elucidating fungal metabolism as well as understanding the physiology and ecology of fungi.