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

Targeting lignocellulolytic gene clusters in novel Trichoderma atroviride and Trichoderma harzianum strains through bacterial artificial chromosome-guided analysis

Lignocellulosic biomass is a complex carbon source with recalcitrant properties whose degradation via industrial enzymatic hydrolysis is challenging, directly affecting the cost of reliable energy production. In nature, filamentous fungi, including Trichoderma species, degrade lignocellulose via an arsenal of hydrolytic and oxidative enzymes that act synergistically to process it into soluble sugar monomers. This work explored the genomic content of Trichoderma atroviride and Trichoderma harzianum strains with hydrolytic abilities by identifying regions possessing degradative enzyme-encoding genes, namely, hydrolytic clusters. We employed bacterial artificial chromosome (BAC) methodology to target specific genomic regions and explore their genetic organization, proximal gene context, and gene expression under degradative conditions. With this tool, it was possible to inspect the linear structure and expression profile of target hydrolytic-rich genomic regions. The present work offers a perspective on the organization of genome regions related to carbohydrate metabolism. This study revealed novel genes and genome regions that are positively regulated during cellulose degradation, contributing to elucidating differences in gene organization that potentially impact hydrolysis among Trichoderma species.

EvoWeaver: large-scale prediction of gene functional associations from coevolutionary signals

The known universe of uncharacterized proteins is expanding far faster than our ability to annotate their functions through laboratory study. Computational annotation approaches rely on similarity to previously studied proteins, thereby ignoring unstudied proteins. Coevolutionary approaches hold promise for injecting new information into our knowledge of the protein universe by linking proteins through 'guilt-by-association'. However, existing coevolutionary algorithms have insufficient accuracy and scalability to connect the entire universe of proteins. We present EvoWeaver, a method that weaves together 12 signals of coevolution to quantify the degree of shared evolution between genes. EvoWeaver accurately identifies proteins involved in protein complexes or separate steps of a biochemical pathway. We show the merits of EvoWeaver by partly reconstructing known biochemical pathways without any prior knowledge other than that available from genomic sequences. Applying EvoWeaver to 1545 gene groups from 8564 genomes reveals missing connections in popular databases and potentially undiscovered links between proteins.

Adaptive evolution of sesquiterpene deoxyphomenone in mycoparasitism by Hansfordia pulvinata associated with horizontal gene transfer from Aspergillus species

Leaf mold caused by the ascomycete fungus Cladosporium fulvum is a devastating disease of tomato plants. The mycoparasitic fungus Hansfordia pulvinata is an effective biocontrol agent that parasitizes C. fulvum hyphae on leaves and secretes 13-deoxyphomenone, an eremophilane-type sesquiterpene, which was also identified as a sporulation-inducing factor in Aspergillus oryzae. Here, we identified deoxyphomenone biosynthesis (DPH) gene clusters conserved in both H. pulvinata and Aspergillus section Flavi, including A. oryzae and A. flavus. Functional disruption of DPH1 orthologous genes encoding sesquiterpene cyclase in H. pulvinata, A. oryzae, and its close relative A. flavus revealed that deoxyphomenone in H. pulvinata had exogenic antifungal activity against C. fulvum and controlled endogenic sporulation in Aspergillus species. Complete DPH clusters, highly similar to those in H. pulvinata, were exclusive to Aspergillus section Flavi, while species in other Aspergillus sections contained fragmented DPH clusters. A comparative genomics analysis revealed that these DPH gene clusters share a common origin and are horizontally transferred from an ancestor of Aspergillus to H. pulvinata. Our results suggest that after horizontal transfer, H. pulvinata maintained the DPH cluster as the inhibitory effect of deoxyphomenone on spore germination and mycelial growth contributed to its mycoparasitism on the host fungus C. fulvum.

IMPORTANCE

Tomato leaf mold disease caused by C. fulvum poses a significant economic threat to tomato production globally. Breeders have developed tomato cultivars with Cf resistance genes. C. fulvum frequently evolves new races that overcome these genetic defenses, complicating control efforts. Additionally, the pathogen has developed resistance to chemical fungicides, prompting the need for sustainable alternatives like biocontrol agents. The mycoparasitic fungus H. pulvinata is crucial as an effective agent against C. fulvum. Clarifying the mechanism of mycoparasitism is significant, as it enhances its application as a biocontrol agent against plant pathogens. This study revealed how H. pulvinata produces deoxyphomenone, an antifungal compound, through horizontal gene transfer from Aspergillus species. It is hypothesized that mycoparasitism could be one of the mechanisms that facilitated horizontal gene transfer between fungi. These insights facilitate the development of eco-friendly, sustainable agricultural practices by reducing dependence on chemical fungicides and promoting natural pathogen control methods.

Fungi: Pioneers of chemical creativity - Techniques and strategies to uncover fungal chemistry

Natural product discovery from fungi for drug development and description of novel chemistry has been a tremendous success. This success is expected to accelerate even further, owing to the advent of sophisticated technical advances of technical advances that recently led to the discovery of an unparalleled biodiversity in the fungal kingdom. This review aims to give an overview on i) important secondary metabolite-derived drugs or drug leads, ii) discuss the analytical and strategic framework of how natural product discovery and drug lead identification transformed from earlier days to the present, iii) how knowledge of fungal biology and biodiversity facilitates the discovery of new compounds, and iv) point out endeavors in understanding fungal secondary metabolite chemistry in order to systematically explore fungal genomes by utilizing synthetic biology. An outlook is given, underlining the necessity for a collaborative and cooperative scenario to harness the full potential of the fungal secondary metabolome.

Overexpression of the RNA-binding protein NrdA affects global gene expression and secondary metabolism in Aspergillus species

RNA-binding protein Nrd1 plays a role in RNA polymerase II transcription termination. In this study, we showed that the orthologous NrdA is important in global mRNA expression and secondary metabolism in Aspergillus species. We constructed an nrdA conditional expression strain using the Tet-On system in Aspergillus luchuenesis mut. kawachii. Downregulation of nrdA caused a severe growth defect, indicating that NrdA is essential for the proliferation of A. kawachii. Parallel RNA-sequencing and RNA immunoprecipitation-sequencing analysis identified potential NrdA-interacting transcripts, corresponding to 32% of the predicted protein-coding genes of A. kawachii. Subsequent gene ontology analysis suggested that overexpression of NrdA affects the production of secondary metabolites. To clarify this, we constructed Aspergillus nidulans, Aspergillus fumigatus, and Aspergillus oryzae strains overexpressing NrdA in the early developmental stage. Overexpression of NrdA reduced the production of sterigmatocystin and penicillin in A. nidulans, as well as that of helvolic acid and pyripyropene A in A. fumigatus. Moreover, it increased the production of kojic acid and reduced the production of penicillin in A. oryzae. These effects were accompanied by almost consistent changes in the mRNA levels of relevant genes. Collectively, these results suggest that NrdA is the essential RNA-binding protein, which plays a significant role in global gene expression and secondary metabolism in Aspergillus species.IMPORTANCENrd1, a component of the Nrd1-Nab3-Sen1 complex, is an essential RNA-binding protein involved in transcriptional termination in yeast. However, its role in filamentous fungi has not been studied. In this study, we characterized an orthologous NrdA in the Aspergillus species, identified potential NrdA-interacting mRNA, and investigated the effect of overexpression of NrdA on mRNA expression in Aspergillus luchuensis mut. kawachii. The results indicated that NrdA controls global gene expression involved in versatile metabolic pathways, including the secondary metabolic process, at least in the early developmental stage. We demonstrated that NrdA overexpression significantly affected the production of secondary metabolites in Aspergillus nidulans, Aspergillus oryzae, and Aspergillus fumigatus. Our findings are of importance to the fungal research community because the secondary metabolism is an industrially and clinically important aspect for the Aspergillus species.

Uncovering the transcriptional landscape of Fomes fomentarius during fungal-based material production through gene co-expression network analysis

BACKGROUND

Fungal-based composites have emerged as renewable, high-performance biomaterials that are produced on lignocellulosic residual streams from forestry and agriculture. Production at an industrial scale promises to revolutionize the world humans inhabit by generating sustainable, low emission, non-toxic and biodegradable construction, packaging, textile, and other materials. The polypore Fomes fomentarius is one of the basidiomycete species used for biomaterial production, yet nothing is known about the transcriptional basis of substrate decomposition, nutrient uptake, or fungal growth during composite formation. Co-expression network analysis based on RNA-Seq profiling has enabled remarkable insights into a range of fungi, and we thus aimed to develop such resources for F. fomentarius.

RESULTS

We analysed gene expression from a wide range of laboratory cultures (n = 9) or biomaterial formation (n = 18) to determine the transcriptional landscape of F. fomentarius during substrate decomposition and to identify genes important for (i) the enzymatic degradation of lignocellulose and other plant-based substrates, (ii) the uptake of their carbon monomers, and (iii) genes guiding mycelium formation through hyphal growth and cell wall biosynthesis. Simple scripts for co-expression network construction were generated and tested, and harnessed to identify a fungal-specific transcription factor named CacA strongly co-expressed with multiple chitin and glucan biosynthetic genes or Rho GTPase encoding genes, suggesting this protein is a high-priority target for engineering adhesion and branching during composite growth. We then updated carbohydrate activated enzymes (CAZymes) encoding gene annotation, used phylogenetics to assign putative uptake systems, and applied network analysis to predict repressing/activating transcription factors for lignocellulose degradation. Finally, we identified entirely new types of co-expressed contiguous clusters not previously described in fungi, including genes predicted to encode CAZymes, hydrophobins, kinases, lipases, F-box domains, chitin synthases, amongst others.

CONCLUSION

The systems biology data generated in this study will enable us to understand the genetic basis of F. fomentarius biomaterial formation in unprecedented detail. We provided proof-of-principle for accurate network-derived predictions of gene function in F. fomentarius and generated the necessary data and scripts for analysis by any end user. Entirely new classes of contiguous co-expressed gene clusters were discovered, and multiple transcription factor encoding genes which are high-priority targets for genetic engineering were identified.

Multivariate Divergence in Wild Microbes: No Evidence for Evolution along a Genetic Line of Least Resistance

AbstractTrait evolution depends both on the direct fitness effects of specific traits and on indirect selection arising from genetically correlated traits. Although well established in plants and animals, the role of trait correlations in microbial evolution remains a major open question. Here, we tested whether genetic correlations in a suite of metabolic traits are conserved between two sister lineages of fungal endophytes and whether phenotypic divergence between lineages occurred in the direction of the multivariate trait combination containing the most genetic variance within lineages (i.e., the genetic lines of least resistance). We found that while one lineage grew faster across nearly all substrates, lineages differed in their mean response to specific substrates and in their overall multivariate metabolic trait means. The structure of the genetic variance-covariance (G) matrix was conserved between lineages, yet to our surprise divergence in metabolic phenotypes between lineages was nearly orthogonal to the major axis of genetic variation within lineages, indicating that divergence did not occur along the genetic lines of least resistance. Our findings suggest that the evolutionary genetics of trait correlations in microorganisms warrant further research and highlight the extensive functional variation that exists at very fine taxonomic scales in host-associated microbial communities.

Comparative genome analysis of endophytic Bacillus amyloliquefaciens MR4: a potential biocontrol agent isolated from wild medicinal plant root tissue

In this investigation, 396 endophytic bacterial strains from six indigenous medicinal plant species within the Xinjiang Tumor Peak National Nature Reserve were subjected to screening. The strain MR4 emerged as a noteworthy contender, demonstrating pronounced biocontrol capabilities coupled with exceptional cold tolerance. Through morphological scrutiny and comprehensive genomic sequencing, MR4 was identified as Bacillus amyloliquefaciens. Antagonistic assays revealed MR4's efficacy in suppressing the causative agents of cotton wilt and verticillium wilt, achieving inhibition rates surpassing 50%. Analyses, underpinned by PCR methodologies, indicated MR4's capacity to biosynthesize a minimum of eight distinct antimicrobial agents. The whole-genome sequencing data indicated that B. amyloliquefaciens MR4 had the genome size and GC content of 4,017,872 bp and 47.14%, respectively, and 4191 coding genes were identified. The genome consists of a single chromosome and one plasmid. Moreover, it was augmented by annotations from various databases, including GO, KEGG, and COG. The pathogenicity of MR4 undergoes evaluation, while predictions concerning MR4's secondary metabolites have disclosed gene clusters for 13 varieties of these compounds, with particular emphasis on surfactins and fengycin. Comparative analyses with four paradigmatic strains shed light on MR4's genomic composition and its phylogenetic lineage within the Bacillus genus. The genomic data pertaining to MR4 have been duly submitted to the NCBI GenBank, bearing the accession numbers CP146236 (Chr1) and CP146237 (plas1). This study endeavors to furnish potent microbial resources for the biocontrol and enhancement of plant growth, thereby providing a theoretical groundwork for MR4's agronomic utilization.

A unique metabolic gene cluster regulates lactose and galactose metabolism in the yeast Candida intermedia

Lactose assimilation is a relatively rare trait in yeasts, and Kluyveromyces yeast species have long served as model organisms for studying lactose metabolism. Meanwhile, the metabolic strategies of most other lactose-assimilating yeasts remain unknown. In this work, we have elucidated the genetic determinants of the superior lactose-growing yeast Candida intermedia. Through genomic and transcriptomic analyses, we identified three interdependent gene clusters responsible for the metabolism of lactose and its hydrolysis product galactose: the conserved LAC cluster (LAC12, LAC4) for lactose uptake and hydrolysis, the conserved GAL cluster (GAL1, GAL7, and GAL10) for galactose catabolism through the Leloir pathway, and a "GALLAC" cluster containing the transcriptional activator gene LAC9, second copies of GAL1 and GAL10, and a XYL1 gene encoding an aldose reductase involved in carbon overflow metabolism. Bioinformatic analysis suggests that the GALLAC cluster is unique to C. intermedia and has evolved through gene duplication and divergence, and deletion mutant phenotyping proved that the cluster is indispensable for C. intermedia's growth on lactose and galactose. We also show that the regulatory network in C. intermedia, governed by Lac9 and Gal1 from the GALLAC cluster, differs significantly from the galactose and lactose regulons in Saccharomyces cerevisiae, Kluyveromyces lactis, and Candida albicans. Moreover, although lactose and galactose metabolism are closely linked in C. intermedia, our results also point to important regulatory differences.IMPORTANCEThis study paves the way to a better understanding of lactose and galactose metabolism in the non-conventional yeast C. intermedia. Notably, the unique GALLAC cluster represents a new, interesting example of metabolic network rewiring and likely helps to explain how C. intermedia has evolved into an efficient lactose-assimilating yeast. With the Leloir pathway of budding yeasts acting like a model system for understanding the function, evolution, and regulation of eukaryotic metabolism, this work provides new evolutionary insights into yeast metabolic pathways and regulatory networks. In extension, the results will facilitate future development and use of C. intermedia as a cell-factory for conversion of lactose-rich whey into value-added products.

CLOCI: unveiling cryptic fungal gene clusters with generalized detection

Gene clusters are genomic loci that contain multiple genes that are functionally and genetically linked. Gene clusters collectively encode diverse functions, including small molecule biosynthesis, nutrient assimilation, metabolite degradation, and production of proteins essential for growth and development. Identifying gene clusters is a powerful tool for small molecule discovery and provides insight into the ecology and evolution of organisms. Current detection algorithms focus on canonical 'core' biosynthetic functions many gene clusters encode, while overlooking uncommon or unknown cluster classes. These overlooked clusters are a potential source of novel natural products and comprise an untold portion of overall gene cluster repertoires. Unbiased, function-agnostic detection algorithms therefore provide an opportunity to reveal novel classes of gene clusters and more precisely define genome organization. We present CLOCI (Co-occurrence Locus and Orthologous Cluster Identifier), an algorithm that identifies gene clusters using multiple proxies of selection for coordinated gene evolution. Our approach generalizes gene cluster detection and gene cluster family circumscription, improves detection of multiple known functional classes, and unveils non-canonical gene clusters. CLOCI is suitable for genome-enabled small molecule mining, and presents an easily tunable approach for delineating gene cluster families and homologous loci.

The gradual establishment of complex coumarin biosynthetic pathway in Apiaceae

Complex coumarins (CCs) represent characteristic metabolites found in Apiaceae plants, possessing significant medical value. Their essential functional role is likely as protectants against pathogens and regulators responding to environmental stimuli. Utilizing genomes and transcriptomes from 34 Apiaceae plants, including our recently sequenced Peucedanum praeruptorum, we conduct comprehensive phylogenetic analyses to reconstruct the detailed evolutionary process of the CC biosynthetic pathway in Apiaceae. Our results show that three key enzymes - p-coumaroyl CoA 2'-hydroxylase (C2'H), C-prenyltransferase (C-PT), and cyclase - originated successively at different evolutionary nodes within Apiaceae through various means of gene duplications: ectopic and tandem duplications. Neofunctionalization endows these enzymes with novel functions necessary for CC biosynthesis, thus completing the pathway. Candidate genes are cloned for heterologous expression and subjected to in vitro enzymatic assays to test our hypothesis regarding the origins of the key enzymes, and the results precisely validate our evolutionary inferences. Among the three enzymes, C-PTs are likely the primary determinant of the structural diversity of CCs (linear/angular), due to divergent activities evolved to target different positions (C-6 or C-8) of umbelliferone. A key amino acid variation (Ala161/Thr161) is identified and proven to play a crucial role in the alteration of enzymatic activity, possibly resulting in distinct binding forms between enzymes and substrates, thereby leading to different products. In conclusion, this study provides a detailed trajectory for the establishment and evolution of the CC biosynthetic pathway in Apiaceae. It explains why only a portion, not all, of Apiaceae plants can produce CCs and reveals the mechanisms of CC structural diversity among different Apiaceae plants.

Bioactivity-driven fungal metabologenomics identifies antiproliferative stemphone analogs and their biosynthetic gene cluster

INTRODUCTION

Fungi biosynthesize chemically diverse secondary metabolites with a wide range of biological activities. Natural product scientists have increasingly turned towards bioinformatics approaches, combining metabolomics and genomics to target secondary metabolites and their biosynthetic machinery. We recently applied an integrated metabologenomics workflow to 110 fungi and identified more than 230 high-confidence linkages between metabolites and their biosynthetic pathways.

OBJECTIVES

To prioritize the discovery of bioactive natural products and their biosynthetic pathways from these hundreds of high-confidence linkages, we developed a bioactivity-driven metabologenomics workflow combining quantitative chemical information, antiproliferative bioactivity data, and genome sequences.

METHODS

The 110 fungi from our metabologenomics study were tested against multiple cancer cell lines to identify which strains produced antiproliferative natural products. Three strains were selected for further study, fractionated using flash chromatography, and subjected to an additional round of bioactivity testing and mass spectral analysis. Data were overlaid using biochemometrics analysis to predict active constituents early in the fractionation process following which their biosynthetic pathways were identified using metabologenomics.

RESULTS

We isolated three new-to-nature stemphone analogs, 19-acetylstemphones G (1), B (2) and E (3), that demonstrated antiproliferative activity ranging from 3 to 5 µM against human melanoma (MDA-MB-435) and ovarian cancer (OVACR3) cells. We proposed a rational biosynthetic pathway for these compounds, highlighting the potential of using bioactivity as a filter for the analysis of integrated-Omics datasets.

CONCLUSIONS

This work demonstrates how the incorporation of biochemometrics as a third dimension into the metabologenomics workflow can identify bioactive metabolites and link them to their biosynthetic machinery.

Tools to make Stachybotrys chartarum genetically amendable: Key to unlocking cryptic biosynthetic gene clusters

The soil and indoor fungus Stachybotrys chartarum can induce respiratory disorders, collectively referred to as stachybotryotoxicosis, owing to its prolific production of diverse bioactive secondary metabolites (SMs) or mycotoxins. Although many of these toxins responsible for the harmful effects on animals and humans have been identified in the genus Stachybotrys, however a number of SMs remain elusive. Through in silico analyses, we have identified 37 polyketide synthase (PKS) genes, highlighting that the chemical profile potential of Stachybotrys is far from being fully explored. Additionally, by leveraging phylogenetic analysis of known SMs produced by non-reducing polyketide synthases (NR-PKS) in other filamentous fungi, we showed that Stachybotrys possesses a rich reservoir of untapped SMs. To unravel natural product biosynthesis in S. chartarum, genetic engineering methods are crucial. For this purpose, we have developed a reliable protocol for the genetic transformation of S. chartarum and applied it to the ScPKS14 biosynthetic gene cluster. This cluster is homologous to the already known Claviceps purpurea CpPKS8 BGC, responsible for the production of ergochromes. While no novel SMs were detected, we successfully applied genetic tools, such as the generation of deletionand overexpression strains of single cluster genes. This toolbox can now be readily employed to unravel not only this particular BGC but also other candidate BGCs present in S. chartarum, making this fungus accessible for genetic engineering.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize diverse secondary metabolites, small organic bioactive molecules with key roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce a cadre of secondary metabolites, some of which are useful (e.g. the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and others harmful (e.g. the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. Fungal genomes often also encode resistance genes that confer protection against toxic secondary metabolites. Some Penicillium species, such as Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity. To investigate the evolutionary conservation of homologs of the gliotoxin BGC and of genes involved in gliotoxin resistance in Penicillium, we analyzed 35 Penicillium genomes from 23 species. Homologous, lesser fragmented gliotoxin BGCs were found in 12 genomes, mostly fragmented remnants of the gliotoxin BGC were found in 21 genomes, whereas the remaining 2 Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes was observed. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes, suggestive of constraint and potential functional significance or more recent decay. Gene tree-species tree reconciliation analyses suggested that the history of homologs in the gliotoxin BGC across the genus Penicillium likely involved multiple duplications, losses, and horizontal gene transfers. Our analyses suggest that genes encoded in BGCs can have complex evolutionary histories and be retained in genomes long after the loss of secondary metabolite biosynthesis.

Genomic and Transcriptomic Comparisons of the Twig Blight Pathogen, Passalora sequoiae, with Mycosphaerellaceae Foliar and Conifer Pathogens

Passalora sequoiae is a foliar pathogen to conifer tree species. In this study, we conducted whole-genome and transcriptome analyses on isolates of P. sequoiae collected from symptomatic Leyland cypress leaves from a Christmas tree farm in Mississippi. The objectives for this research were to elucidate the pathogenicity mechanisms of P. sequoiae by characterizing the genome and transcriptome and possibly identify unique and shared predicted genes in comparison with non-conifer/canker and foliar pathogens in the family Mycosphaerellaceae. P. sequoiae was found to be similar to other foliar Mycosphaerellaceae pathogens and likely represents a hemibiotrophic lifestyle based on comparisons across pathogens. The genome and in planta transcriptome highlighted some unique features of P. sequoiae: the significant presence of chitin synthases and fructose-degrading carbohydrate-degrading enzymes, trans-AT PKS genes, and antibiotic gene clusters that were unique to P. sequoiae compared with the other Mycosphaerellaceae species genomes. Several transcripts that were highly expressed in planta were identified as effectors, yet the functions were not characterized. These targets provide ample resources to continue to characterize pathogen-conifer host interactions in conifer foliar pathogens. Furthermore, this research helps build genomic resources for an important plant pathogen on Leyland cypress that will further our ability to develop novel management practices that could begin with breeding for resistance.

Predicting fungal secondary metabolite activity from biosynthetic gene cluster data using machine learning

Fungal secondary metabolites (SMs) contribute to the diversity of fungal ecological communities, niches, and lifestyles. Many fungal SMs have one or more medically and industrially important activities (e.g., antifungal, antibacterial, and antitumor). The genes necessary for fungal SM biosynthesis are typically located right next to each other in the genome and are known as biosynthetic gene clusters (BGCs). However, whether fungal SM bioactivity can be predicted from specific attributes of genes in BGCs remains an open question. We adapted machine learning models that predicted SM bioactivity from bacterial BGC data with accuracies as high as 80% to fungal BGC data. We trained our models to predict the antibacterial, antifungal, and cytotoxic/antitumor bioactivity of fungal SMs on two data sets: (i) fungal BGCs (data set comprised of 314 BGCs) and (ii) fungal (314 BGCs) and bacterial BGCs (1,003 BGCs). We found that models trained on fungal BGCs had balanced accuracies between 51% and 68%, whereas training on bacterial and fungal BGCs had balanced accuracies between 56% and 68%. The low prediction accuracy of fungal SM bioactivities likely stems from the small size of the data set; this lack of data, coupled with our finding that including bacterial BGC data in the training data did not substantially change accuracies currently limits the application of machine learning approaches to fungal SM studies. With >15,000 characterized fungal SMs, millions of putative BGCs in fungal genomes, and increased demand for novel drugs, efforts that systematically link fungal SM bioactivity to BGCs are urgently needed.IMPORTANCEFungi are key sources of natural products and iconic drugs, including penicillin and statins. DNA sequencing has revealed that there are likely millions of biosynthetic pathways in fungal genomes, but the chemical structures and bioactivities of >99% of natural products produced by these pathways remain unknown. We used artificial intelligence to predict the bioactivities of diverse fungal biosynthetic pathways. We found that the accuracies of our predictions were generally low, between 51% and 68%, likely because the natural products and bioactivities of only very few fungal pathways are known. With >15,000 characterized fungal natural products, millions of putative biosynthetic pathways present in fungal genomes, and increased demand for novel drugs, our study suggests that there is an urgent need for efforts that systematically identify fungal biosynthetic pathways, their natural products, and their bioactivities.

Endosymbiotic Fungal Diversity and Dynamics of the Brown Planthopper across Developmental Stages, Tissues, and Sexes Revealed Using Circular Consensus Sequencing

Endosymbiotic fungi play an important role in the growth and development of insects. Understanding the endosymbiont communities hosted by the brown planthopper (BPH; Nilaparvata lugens Stål), the most destructive pest in rice, is a prerequisite for controlling BPH rice infestations. However, the endosymbiont diversity and dynamics of the BPH remain poorly studied. Here, we used circular consensus sequencing (CCS) to obtain 87,131 OTUs (operational taxonomic units), which annotated 730 species of endosymbiotic fungi in the various developmental stages and tissues. We found that three yeast-like symbionts (YLSs), Polycephalomyces prolificus, Ophiocordyceps heteropoda, and Hirsutella proturicola, were dominant in almost all samples, which was especially pronounced in instar nymphs 4-5, female adults, and the fat bodies of female and male adult BPH. Interestingly, honeydew as the only in vitro sample had a unique community structure. Various diversity indices might indicate the different activity of endosymbionts in these stages and tissues. The biomarkers analyzed using LEfSe suggested some special functions of samples at different developmental stages of growth and the active functions of specific tissues in different sexes. Finally, we found that the incidence of occurrence of three species of Malassezia and Fusarium sp. was higher in males than in females in all comparison groups. In summary, our study provides a comprehensive survey of symbiotic fungi in the BPH, which complements the previous research on YLSs. These results offer new theoretical insights and practical implications for novel pest management strategies to understand the BPH-microbe symbiosis and devise effective pest control strategies.

Endophyte genomes support greater metabolic gene cluster diversity compared with non-endophytes in Trichoderma

Trichoderma is a cosmopolitan genus with diverse lifestyles and nutritional modes, including mycotrophy, saprophytism, and endophytism. Previous research has reported greater metabolic gene repertoires in endophytic fungal species compared to closely-related non-endophytes. However, the extent of this ecological trend and its underlying mechanisms are unclear. Some endophytic fungi may also be mycotrophs and have one or more mycoparasitism mechanisms. Mycotrophic endophytes are prominent in certain genera like Trichoderma, therefore, the mechanisms that enable these fungi to colonize both living plants and fungi may be the result of expanded metabolic gene repertoires. Our objective was to determine what, if any, genomic features are overrepresented in endophytic fungi genomes in order to undercover the genomic underpinning of the fungal endophytic lifestyle. Here we compared metabolic gene cluster and mycoparasitism gene diversity across a dataset of thirty-eight Trichoderma genomes representing the full breadth of environmental Trichoderma's diverse lifestyles and nutritional modes. We generated four new Trichoderma endophyticum genomes to improve the sampling of endophytic isolates from this genus. As predicted, endophytic Trichoderma genomes contained, on average, more total biosynthetic and degradative gene clusters than non-endophytic isolates, suggesting that the ability to create/modify a diversity of metabolites potential is beneficial or necessary to the endophytic fungi. Still, once the phylogenetic signal was taken in consideration, no particular class of metabolic gene cluster was independently associated with the Trichoderma endophytic lifestyle. Several mycoparasitism genes, but no chitinase genes, were associated with endophytic Trichoderma genomes. Most genomic differences between Trichoderma lifestyles and nutritional modes are difficult to disentangle from phylogenetic divergences among species, suggesting that Trichoderma genomes maybe particularly well-equipped for lifestyle plasticity. We also consider the role of endophytism in diversifying secondary metabolism after identifying the horizontal transfer of the ergot alkaloid gene cluster to Trichoderma.

Wheat stripe rust resistance locus YR63 is a hot spot for evolution of defence genes - a pangenome discovery

BACKGROUND

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), poses a threat to global wheat production. Deployment of widely effective resistance genes underpins management of this ongoing threat. This study focused on the mapping of stripe rust resistance gene YR63 from a Portuguese hexaploid wheat landrace AUS27955 of the Watkins Collection.

RESULTS

YR63 exhibits resistance to a broad spectrum of Pst races from Australia, Africa, Asia, Europe, Middle East and South America. It was mapped to the short arm of chromosome 7B, between two single nucleotide polymorphic (SNP) markers sunCS_YR63 and sunCS_67, positioned at 0.8 and 3.7 Mb, respectively, in the Chinese Spring genome assembly v2.1. We characterised YR63 locus using an integrated approach engaging targeted genotyping-by-sequencing (tGBS), mutagenesis, resistance gene enrichment and sequencing (MutRenSeq), RNA sequencing (RNASeq) and comparative genomic analysis with tetraploid (Zavitan and Svevo) and hexaploid (Chinese Spring) wheat genome references and 10+ hexaploid wheat genomes. YR63 is positioned at a hot spot enriched with multiple nucleotide-binding and leucine rich repeat (NLR) and kinase domain encoding genes, known widely for defence against pests and diseases in plants and animals. Detection of YR63 within these gene clusters is not possible through short-read sequencing due to high homology between members. However, using the sequence of a NLR member we were successful in detecting a closely linked SNP marker for YR63 and validated on a panel of Australian bread wheat, durum and triticale cultivars.

CONCLUSIONS

This study highlights YR63 as a valuable source for resistance against Pst in Australia and elsewhere. The closely linked SNP marker will facilitate rapid introgression of YR63 into elite cultivars through marker-assisted selection. The bottleneck of this study reinforces the necessity for a long-read sequencing such as PacBio or Oxford Nanopore based techniques for accurate detection of the underlying resistance gene when it is part of a large gene cluster.

Functional Crosstalk between Discrete Indole Terpenoid Gene Clusters in Tolypocladium album

Indole terpenoids make up a large group of secondary metabolites that display an enticing array of bioactivities. While indole diterpene (IDT) and rarely indole sesquiterpene (IST) pathways have been found individually in filamentous fungi, here we show that both cluster types are encoded within the genome of Tolypocladium album. Through heterologous reconstruction, we demonstrate the SES cluster encodes for IST biosynthesis and can tailor IDT substrates produced by the TER cluster.

Structural diversification of natural substrates modified by the O-methyltransferase AurJ from Fusarium Graminearum

Aromatic polyketide and phenylpropanoid derivatives are a large class of natural products produced by bacteria, fungi, and plants. The O-methylation is a unique decoration that can increase structural diversity of aromatic compounds and improve their pharmacological properties, but the substrate specificity of O-methyltransferase hinders the discovery of more natural products with O-methylation through biosynthesis. Here, we reported that the O-methyltransferase AurJ from plant pathogenic fungus Fusarium graminearum could methylate a broad range of natural substrates of monocyclic, bicyclic, and tricyclic aromatic precursors, exhibiting excellent substrate tolerance. This finding will partly change our stereotype about the specificity of traditional methyltransferases, and urge us to mine more O-methyltransferases with good substrate tolerance and discover more methylated natural products for drug discovery and development through directed evolution and combinatorial biosynthesis.

Unraveling the evolutionary dynamics of the TPS gene family in land plants

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases (TPSs) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family.

Unveiling the hidden role of aquatic viruses in hydrocarbon pollution bioremediation

Hydrocarbon pollution poses substantial environmental risks to water and soil. Bioremediation, which utilizes microorganisms to manage pollutants, offers a cost-effective solution. However, the role of viruses, particularly bacteriophages (phages), in bioremediation remains unexplored. This study examines the diversity and activity of hydrocarbon-degradation genes encoded by environmental viruses, focusing on phages, within public databases. We identified 57 high-quality phage-encoded auxiliary metabolic genes (AMGs) related to hydrocarbon degradation, which we refer to as virus-encoded hydrocarbon degradation genes (vHYDEGs). These genes are encoded by taxonomically diverse aquatic phages and highlight the under-characterized global virosphere. Six protein families involved in the initial alkane hydroxylation steps were identified. Phylogenetic analyses revealed the diverse evolutionary trajectories of vHYDEGs across habitats, revealing previously unknown biodegraders linked evolutionarily with vHYDEGs. Our findings suggest phage AMGs may contribute to alkane and aromatic hydrocarbon degradation, participating in the initial, rate-limiting hydroxylation steps, thereby aiding hydrocarbon pollution bioremediation and promoting their propagation. To support future research, we developed vHyDeg, a database containing identified vHYDEGs with comprehensive annotations, facilitating the screening of hydrocarbon degradation AMGs and encouraging their bioremediation applications.

Biosynthetic gene cluster synteny: Orthologous polyketide synthases in Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus), which together generate a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here, we provide a comparative view of the biosynthetic gene clusters of three lichen mycobionts derived from Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata. In addition, we present a high-quality PacBio metagenome of Parmelia sulcata, from which we extracted the mycobiont bin containing 214 biosynthetic gene clusters. Most biosynthetic gene clusters in these genomes were associated with T1PKSs, followed by NRPSs and terpenes. This study focused on biosynthetic gene clusters related to polyketide synthesis. Based on ketosynthase homology, we identified nine highly syntenic clusters present in all three species. Among the four clusters belonging to nonreducing PKSs, two are putatively linked to lichen substances derived from orsellinic acid (orcinol depsides and depsidones, e.g., lecanoric acid, physodic acid, lobaric acid), one to compounds derived from methylated forms of orsellinic acid (beta orcinol depsides, e.g., atranorin), and one to melanins. Five clusters with orthologs in all three species are linked to reducing PKSs. Our study contributes to sorting and dereplicating the vast PKS diversity found in lichenized fungi. High-quality sequences of biosynthetic gene clusters of these three common species provide a foundation for further exploration into biotechnological applications and the molecular evolution of lichen substances.

Predicting fungal secondary metabolite activity from biosynthetic gene cluster data using machine learning

Fungal secondary metabolites (SMs) play a significant role in the diversity of ecological communities, niches, and lifestyles in the fungal kingdom. Many fungal SMs have medically and industrially important properties including antifungal, antibacterial, and antitumor activity, and a single metabolite can display multiple types of bioactivities. The genes necessary for fungal SM biosynthesis are typically found in a single genomic region forming biosynthetic gene clusters (BGCs). However, whether fungal SM bioactivity can be predicted from specific attributes of genes in BGCs remains an open question. We adapted previously used machine learning models for predicting SM bioactivity from bacterial BGC data to fungal BGC data. We trained our models to predict antibacterial, antifungal, and cytotoxic/antitumor bioactivity on two datasets: 1) fungal BGCs (dataset comprised of 314 BGCs), and 2) fungal (314 BGCs) and bacterial BGCs (1,003 BGCs); the second dataset was our control since a previous study using just the bacterial BGC data yielded prediction accuracies as high as 80%. We found that the models trained only on fungal BGCs had balanced accuracies between 51-68%, whereas training on bacterial and fungal BGCs yielded balanced accuracies between 61-74%. The lower accuracy of the predictions from fungal data likely stems from the small number of BGCs and SMs with known bioactivity; this lack of data currently limits the application of machine learning approaches in studying fungal secondary metabolism. However, our data also suggest that machine learning approaches trained on bacterial and fungal data can predict SM bioactivity with good accuracy. With more than 15,000 characterized fungal SMs, millions of putative BGCs present in fungal genomes, and increased demand for novel drugs, efforts that systematically link fungal SM bioactivity to BGCs are urgently needed.

Emergence of the fungal immune system

Investigation of fungal biology has been frequently motivated by the fact that many fungal species are important plant and animal pathogens. Such efforts have contributed significantly toward our understanding of fungal pathogenic lifestyles (virulence factors and strategies) and the interplay with host immune systems. In parallel, work on fungal allorecognition systems leading to the characterization of fungal regulated cell death determinants and pathways, has been instrumental for the emergent concept of fungal immunity. The uncovered evolutionary trans-kingdom parallels between fungal regulated cell death pathways and innate immune systems incite us to reflect further on the concept of a fungal immune system. Here, I briefly review key findings that have shaped the fungal immunity paradigm, providing a perspective on what I consider its most glaring knowledge gaps. Undertaking to fill such gaps would establish firmly the fungal immune system inside the broader field of comparative immunology.

Biosynthetic Potential of Hypogymnia Holobionts: Insights into Secondary Metabolite Pathways

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus). They are known to produce a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here we provide a comprehensive view of the biosynthetic gene clusters of all organisms comprising a lichen thallus: fungi, green algae, and bacteria. We present two high-quality PacBio metagenomes, in which we identified a total of 460 biosynthetic gene clusters. Lichen mycobionts yielded 73-114 clusters, other lichen associated ascomycetes 8-40, green algae of the genus Trebouxia 14-19, and lichen-associated bacteria 101-105 clusters. The mycobionts contained mainly T1PKSs, followed by NRPSs, and terpenes; Trebouxia reads harbored mainly clusters linked to terpenes, followed by NRPSs and T3PKSs. Other lichen-associated ascomycetes and bacteria contained a mix of diverse biosynthetic gene clusters. In this study, we identified for the first time the biosynthetic gene clusters of entire lichen holobionts. The yet untapped biosynthetic potential of two species of the genus Hypogymnia is made accessible for further research.

A genetic tool to express long fungal biosynthetic genes

BACKGROUND

Secondary metabolites (SMs) from mushroom-forming fungi (Basidiomycota) and early diverging fungi (EDF) such as Mucoromycota are scarcely investigated. In many cases, production of SMs is induced by unknown stress factors or is accompanied by seasonable developmental changes on fungal morphology. Moreover, many of these fungi are considered as non-culturable under laboratory conditions which impedes investigation into SM. In the post-genomic era, numerous novel SM genes have been identified especially from EDF. As most of them encode multi-module enzymes, these genes are usually long which limits cloning and heterologous expression in traditional hosts.

RESULTS

An expression system in Aspergillus niger is presented that is suitable for the production of SMs from both Basidiomycota and EDF. The akuB gene was deleted in the expression host A. niger ATNT∆pyrG, resulting in a deficient nonhomologous end-joining repair mechanism which in turn facilitates the targeted gene deletion via homologous recombination. The ∆akuB mutant tLK01 served as a platform to integrate overlapping DNA fragments of long SM genes into the fwnA locus required for the black pigmentation of conidia. This enables an easy discrimination of correct transformants by screening the transformation plates for fawn-colored colonies. Expression of the gene of interest (GOI) is induced dose-dependently by addition of doxycycline and is enhanced by the dual TetON/terrein synthase promoter system (ATNT) from Aspergillus terreus. We show that the 8 kb polyketide synthase gene lpaA from the basidiomycete Laetiporus sulphureus is correctly assembled from five overlapping DNA fragments and laetiporic acids are produced. In a second approach, we expressed the yet uncharacterized > 20 kb nonribosomal peptide synthetase gene calA from the EDF Mortierella alpina. Gene expression and subsequent LC-MS/MS analysis of mycelial extracts revealed the production of the antimycobacterial compound calpinactam. This is the first report on the heterologous production of a full-length SM multidomain enzyme from EDF.

CONCLUSIONS

The system allows the assembly, targeted integration and expression of genes of > 20 kb size in A. niger in one single step. The system is suitable for evolutionary distantly related SM genes from both Basidiomycota and EDF. This uncovers new SM resources including genetically intractable or non-culturable fungi.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

The evolution of the gliotoxin biosynthetic gene cluster in Penicillium fungi

Fungi biosynthesize a diversity of secondary metabolites, small organic bioactive molecules that play diverse roles in fungal ecology. Fungal secondary metabolites are often encoded by physically clustered sets of genes known as biosynthetic gene clusters (BGCs). Fungi in the genus Penicillium produce diverse secondary metabolites that have been both useful (e.g., the antibiotic penicillin and the cholesterol-lowering drug mevastatin) and harmful (e.g., the mycotoxin patulin and the immunosuppressant gliotoxin) to human affairs. BGCs often also encode resistance genes that confer self-protection to the secondary metabolite-producing fungus. Some Penicillium species, such as Penicillium lilacinoechinulatum and Penicillium decumbens, are known to produce gliotoxin, a secondary metabolite with known immunosuppressant activity; however, an evolutionary characterization of the BGC responsible for gliotoxin biosynthesis among Penicillium species is lacking. Here, we examine the conservation of genes involved in gliotoxin biosynthesis and resistance in 35 Penicillium genomes from 23 species. We found homologous, less fragmented gliotoxin BGCs in 12 genomes, mostly fragmented remnants of the gliotoxin BGC in 21 genomes, whereas the remaining two Penicillium genomes lacked the gliotoxin BGC altogether. In contrast, we observed broad conservation of homologs of resistance genes that reside outside the BGC across Penicillium genomes. Evolutionary rate analysis revealed that BGCs with higher numbers of genes evolve slower than BGCs with few genes. Even though the gliotoxin BGC is fragmented to varying degrees in nearly all genomes examined, ancestral state reconstruction suggests that the ancestor of Penicillium species possessed the gliotoxin BGC. Our analyses suggest that genes that are part of BGCs can be retained in genomes long after the loss of secondary metabolite biosynthesis.

Metagenomics Shines Light on the Evolution of "Sunscreen" Pigment Metabolism in the Teloschistales (Lichen-Forming Ascomycota)

Fungi produce a vast number of secondary metabolites that shape their interactions with other organisms and the environment. Characterizing the genes underpinning metabolite synthesis is therefore key to understanding fungal evolution and adaptation. Lichenized fungi represent almost one-third of Ascomycota diversity and boast impressive secondary metabolites repertoires. However, most lichen biosynthetic genes have not been linked to their metabolite products. Here we used metagenomic sequencing to survey gene families associated with production of anthraquinones, UV-protectant secondary metabolites present in various fungi, but especially abundant in a diverse order of lichens, the Teloschistales (class Lecanoromycetes, phylum Ascomycota). We successfully assembled 24 new, high-quality lichenized-fungal genomes de novo and combined them with publicly available Lecanoromycetes genomes from taxa with diverse secondary chemistry to produce a whole-genome tree. Secondary metabolite biosynthetic gene cluster (BGC) analysis showed that whilst lichen BGCs are numerous and highly dissimilar, core enzyme genes are generally conserved across taxa. This suggests metabolite diversification occurs via re-shuffling existing enzyme genes with novel accessory genes rather than BGC gains/losses or de novo gene evolution. We identified putative anthraquinone BGCs in our lichen dataset that appear homologous to anthraquinone clusters from non-lichenized fungi, suggesting these genes were present in the common ancestor of the subphylum Pezizomycotina. Finally, we identified unique transporter genes in Teloschistales anthraquinone BGCs that may explain why these metabolites are so abundant and ubiquitous in these lichens. Our results support the importance of metagenomics for understanding the secondary metabolism of non-model fungi such as lichens.

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.

A putative terpene cyclase gene (CcPtc1) is required for fungal development and virulence in Cytospora chrysosperma

Cytospora chrysosperma is a destructive plant pathogenic fungus, which causes canker disease on numerous woody plants. However, knowledge concerning the interaction between C. chrysosperma and its host remains limited. Secondary metabolites produced by phytopathogens often play important roles in their virulence. Terpene cyclases (TC), polyketide synthases (PKS) and non-ribosomal peptide synthetases (NRPS) are the key components for the synthesis of secondary metabolites. Here, we characterized the functions of a putative terpene type secondary metabolite biosynthetic core gene CcPtc1 in C. chrysosperma, which was significantly up-regulated in the early stages of infection. Importantly, deletion of CcPtc1 greatly reduced fungal virulence to the poplar twigs and they also showed significantly reduced fungal growth and conidiation compared with the wild-type (WT) strain. Furthermore, toxicity test of the crude extraction from each strain showed that the toxicity of crude extraction secreted by ΔCcPtc1 were strongly compromised in comparison with the WT strain. Subsequently, the untargeted metabolomics analyses between ΔCcPtc1 mutant and WT strain were conducted, which revealed 193 significantly different abundant metabolites (DAMs) inΔCcPtc1 mutant compared to the WT strain, including 90 significantly downregulated metabolites and 103 significantly up-regulated metabolites, respectively. Among them, four key metabolic pathways that reported to be important for fungal virulence were enriched, including pantothenate and coenzyme A (CoA) biosynthesis. Moreover, we also detected significant alterations in a series of terpenoids, among which (+)-ar-turmerone, pulegone, ethyl chrysanthemumate, and genipin were significantly down-regulated, while cuminaldehyde and (±)-abscisic acid were significantly up-regulated. In conclusion, our results demonstrated that CcPtc1 acts as a virulence-related secondary metabolism factor and provides new insights into the pathogenesis of C. chrysosperma.

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates,1,2,3 bilaterian animals,4,5 and filamentous fungi6). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ∼400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization.7 We found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ∼100 million years, the microsynteny of many genes-especially genes in metabolic clusters (e.g., in the GAL gene cluster8)-was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

Transperons: RNA operons as effectors of coordinated gene expression in eukaryotes

Coordinated gene expression allows spatiotemporal control of cellular processes and is achieved by the cotranscription/translation of functionally related genes/proteins. Prokaryotes evolved polycistronic messages (operons) to confer expression from a single promoter to efficiently cotranslate proteins functioning on the same pathway. Yet, despite having far greater diversity (e.g., gene number, distribution, modes of expression), eukaryotic cells employ individual promoters and monocistronic messages. Although gene expression is modular, it does not account for how eukaryotes achieve coordinated localized translation. The RNA operon theory states that mRNAs derived from different chromosomes assemble into ribonucleoprotein particles (RNPs) that act as functional operons to generate protein cohorts upon cotranslation. Work in yeast has now validated this theory and shown that intergenic associations and noncanonical histone functions create pathway-specific RNA operons (transperons) that regulate cell physiology. Herein the involvement of chromatin organization in transperon formation and programmed gene coexpression is discussed.

A widespread inversion polymorphism conserved among Saccharomyces species is caused by recurrent homogenization of a sporulation gene family

Saccharomyces genomes are highly collinear and show relatively little structural variation, both within and between species of this yeast genus. We investigated the only common inversion polymorphism known in S. cerevisiae, which affects a 24-kb 'flip/flop' region containing 15 genes near the centromere of chromosome XIV. The region exists in two orientations, called reference (REF) and inverted (INV). Meiotic recombination in this region is suppressed in crosses between REF and INV orientation strains such as the BY x RM cross. We find that the inversion polymorphism is at least 17 million years old because it is conserved across the genus Saccharomyces. However, the REF and INV isomers are not ancient alleles but are continually being re-created by re-inversion of the region within each species. Inversion occurs due to continual homogenization of two almost identical 4-kb sequences that form an inverted repeat (IR) at the ends of the flip/flop region. The IR consists of two pairs of genes that are specifically and strongly expressed during the late stages of sporulation. We show that one of these gene pairs, YNL018C/YNL034W, codes for a protein that is essential for spore formation. YNL018C and YNL034W are the founder members of a gene family, Centroid, whose members in other Saccharomycetaceae species evolve fast, duplicate frequently, and are preferentially located close to centromeres. We tested the hypothesis that Centroid genes are a meiotic drive system, but found no support for this idea.

Functional Guilds, Community Assembly, and Co-occurrence Patterns of Fungi in Metalliferous Mine Tailings Ponds in Mainland China

Metalliferous mine tailings ponds are generally characterized by low levels of nutrient elements, sustained acidic conditions, and high contents of toxic metals. They represent one kind of extreme environments that are believed to resemble the Earth's early environmental conditions. There is increasing evidence that the diversity of fungi inhabiting mine tailings ponds is much higher than previously thought. However, little is known about functional guilds, community assembly, and co-occurrence patterns of fungi in such habitats. As a first attempt to address this critical knowledge gap, we employed high-throughput sequencing to characterize fungal communities in 33 mine tailings ponds distributed across 18 provinces of mainland China. A total of 5842 fungal phylotypes were identified, with saprotrophic fungi being the major functional guild. The predictors of fungal diversity in whole community and sub-communities differed considerably. Community assembly of the whole fungal community and individual functional guilds were primarily governed by stochastic processes. Total soil nitrogen and total phosphorus mediated the balance between stochastic and deterministic processes of the fungal community assembly. Co-occurrence network analysis uncovered a high modularity of the whole fungal community. The observed main modules largely consisted of saprotrophic fungi as well as various phylotypes that could not be assigned to known functional guilds. The richness of core fungal phylotypes, occupying vital positions in co-occurrence network, was positively correlated with edaphic properties such as soil enzyme activity. This indicates the important roles of core fungal phylotypes in soil organic matter decomposition and nutrient cycling. These findings improve our understanding of fungal ecology of extreme environments.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

RimO (SrrB) is required for carbon starvation signaling and production of secondary metabolites in Aspergillus nidulans

Depending on the prevailing environmental, developmental and nutritional conditions, fungi activate biosynthetic gene clusters (BGCs) to produce condition-specific secondary metabolites (SMs). For activation, global chromatin-based de-repression must be integrated with pathway-specific induction signals. Here we describe a new global regulator needed to activate starvation-induced SMs. In our transcriptome dataset, we found locus AN7572 strongly transcribed solely under conditions of starvation-induced SM production. The predicted AN7572 protein is most similar to the stress and nutritional regulator Rim15 of Saccharomyces cerevisiae, and to STK-12 of Neurospora crassa. Based on this similarity and on stress and nutritional response phenotypes of A. nidulans knock-out and overexpression strains, AN7572 is designated rimO. In relation to SM production, we found that RimO is required for the activation of starvation-induced BGCs, including the sterigmatocystin (ST) gene cluster. Here, RimO regulates the pathway-specific transcription factor AflR both at the transcriptional and post-translational level. At the transcriptional level, RimO mediates aflR induction following carbon starvation and at the post-translational level, RimO is required for nuclear accumulation of the AflR protein. Genome-wide transcriptional profiling showed that cells lacking rimO fail to adapt to carbon starvation that, in the wild type, leads to down-regulation of genes involved in basic metabolism, membrane biogenesis and growth. Consistently, strains overexpressing rimO are more resistant to oxidative and osmotic stress, largely insensitive to glucose repression and strongly overproduce several SMs. Our data indicate that RimO is a positive regulator within the SM and stress response network, but this requires nutrient depletion that triggers both, rimO gene transcription and activation of the RimO protein.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Examination of Genome-Wide Ortholog Variation in Clinical and Environmental Isolates of the Fungal Pathogen Aspergillus fumigatus

Aspergillus fumigatus is both an environmental saprobe and an opportunistic human fungal pathogen. Knowledge of genomic variation across A. fumigatus isolates is essential for understanding the evolution of pathogenicity, virulence, and resistance to antifungal drugs. Here, we investigated 206 A. fumigatus isolates (133 clinical and 73 environmental isolates), aiming to identify genes with variable presence across isolates and test whether this variation was related to the clinical or environmental origin of isolates. The PanOrtho genome of A. fumigatus consists of 13,085 ortholog groups, of which 7,773 (59.4%) are shared by all isolates (core groups) and 5,312 (40.6%) vary in their gene presence across isolates (accessory groups plus singletons). Despite differences in the distribution of orthologs across all isolates, no significant differences were observed among clinical versus environmental isolates when phylogeny was accounted for. Orthologs that differ in their distribution across isolates tend to occur at low frequency and/or be restricted to specific isolates; thus, the degree of genomic conservation between orthologs of A. fumigatus is high. These results suggest that differences in the distribution of orthologs within A. fumigatus cannot be associated with the clinical or environmental origin of isolates. IMPORTANCE Aspergillus fumigatus is a cosmopolitan species of fungus responsible for thousands of cases of invasive disease annually. Clinical and environmental isolates of A. fumigatus exhibit extensive phenotypic differences, including differences related to virulence and antifungal drug resistance. A comprehensive survey of the genomic diversity present in A. fumigatus and its relationship to the clinical or environmental origin of isolates can contribute to the prediction of the mechanisms of evolution and infection of the species. Our results suggest that there is no significant variation in ortholog distribution between clinical and environmental isolates when accounting for evolutionary history. The work supports the hypothesis that environmental and clinical isolates of A. fumigatus do not differ in their gene contents.

Evolution of pathogenicity in obligate fungal pathogens and allied genera

Obligate fungal pathogens (ascomycetes and basidiomycetes) and oomycetes are known to cause diseases in cereal crop plants. They feed on living cells and most of them have learned to bypass the host immune machinery. This paper discusses some of the factors that are associated with pathogenicity drawing examples from ascomycetes, basidiomycetes and oomycetes, with respect to their manifestation in crop plants. The comparisons have revealed a striking similarity in the three groups suggesting convergent pathways that have arisen from three lineages independently leading to an obligate lifestyle. This review has been written with the intent, that new information on adaptation strategies of biotrophs, modifications in pathogenicity strategies and population dynamics will improve current strategies for breeding with stable resistance.

Functional Characterization of the GlcNAc Catabolic Pathway in Cryptococcus deneoformans

The amino sugar N-acetyl-d-glucosamine (GlcNAc) is the key constituent of cell wall components and plays an important role in pathogenesis in a wide range of fungi. However, catabolism of GlcNAc has not been studied in basidiomycete fungi. In this study, we identified and characterized a gene cluster essential for GlcNAc utilization in Cryptococcus deneoformans, an environmental human fungal pathogen. The C. deneoformans genome contains a GlcNAc transporter (Ngt1), a GlcNAc kinase (Hxk3), a GlcNAc-6-phosphate deacetylase (Dac1), and a glucosamine-6-phosphate deaminase (Nag1). Their expression levels were highly induced in cultures containing GlcNAc as the sole carbon source, and the corresponding mutants showed severe growth defects in the presence of GlcNAc. Functional and biochemical analyses revealed that HXK3 encodes a novel GlcNAc kinase. Site-directed mutations of conserved residues of Hxk3 indicated that ATP binding and GlcNAc binding are essential for GlcNAc kinase activities. Taken together, the results from this study provide crucial insights into basidiomycete GlcNAc catabolism. IMPORTANCEN-Acetylglucosamine (GlcNAc) is recognized as not only the building block of chitin but also an important signaling molecule in fungi. The catabolic pathway of GlcNAc also plays an important role in vital biological processes in fungi. However, the utilization pathway of GlcNAc in the phylum Basidiomycota, which contains more than 41,000 species, remains unknown. Cryptococcus deneoformans is a representative basidiomycetous pathogen that causes life-threatening meningitis. In this study, we characterized a gene cluster essential for GlcNAc utilization in C. deneoformans and identified a novel GlcNAc kinase. The results of this study provide important insights into basidiomycete GlcNAc catabolism and offer a starting point for revealing its role in pathogenesis.

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.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker's yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

Fungal-derived compounds and mycogenic nanoparticles with antimycobacterial activity: a review

Tuberculosis (TB) is a persistent lung infection caused by Mycobacterium tuberculosis. The disease is characterized by high mortality rates of over 1 million per year. Unfortunately, the potency and effectiveness of currently used anti-TB drugs is gradually decreasing due to the constant development of persistence and resistance by M. tuberculosis. The adverse side effects associated with current anti-TB drugs, along with anti-TB drug resistance, present an opportunity to bio-prospect novel potent anti-TB drugs from unique sources. Fundamentally, fungi are a rich source of bioactive secondary metabolites with valuable therapeutic potential. Enhancing the potency and effectiveness of fungal-based anti-TB drug leads by chemical synthesis and/or modification with nanomaterials, may result in the discovery of novel anti-TB drugs. In this review, the antimycobacterial activity of fungal-derived compounds and mycogenic nanoparticles are summarized. Numerous fungal-derived compounds as well as some mycogenic nanoparticles that exhibit strong antimycobacterial activity that is comparable to that of approved drugs, were found. If fully explored, fungi holds the promise to become key drivers in the generation of lead compounds in TB-drug discovery initiatives.

Contrasting Strategies for Sucrose Utilization in a Floral Yeast Clade

Yeast species in the Wickerhamiella and Starmerella genera (W/S clade) thrive in the sugar-rich floral niche. We have previously shown that species belonging to this clade harbor an unparalleled number of genes of bacterial origin, among which is the SUC2 gene, encoding a sucrose-hydrolyzing enzyme. In this study, we used complementary in silico and experimental approaches to examine sucrose utilization in a broader cohort of species representing extant diversity in the W/S clade. Distinct strategies and modes of sucrose assimilation were unveiled, involving either extracellular sucrose hydrolysis through secreted bacterial Suc2 or intracellular assimilation using broad-substrate-range α-glucoside/H⁺ symporters and α-glucosidases. The intracellular pathway is encoded in two types of gene clusters reminiscent of the MAL clusters in Saccharomyces cerevisiae, where they are involved in maltose utilization. The genes composing each of the two types of MAL clusters found in the W/S clade have disparate evolutionary histories, suggesting that they formed de novo. Both transporters and glucosidases were shown to be functional and additionally involved in the metabolization of other disaccharides, such as maltose and melezitose. In one Wickerhamiella species lacking the α-glucoside transporter, maltose assimilation is accomplished extracellularly, an attribute which has been rarely observed in fungi. Sucrose assimilation in Wickerhamiella generally escaped both glucose repression and the need for an activator and is thus essentially constitutive, which is consistent with the abundance of both glucose and sucrose in the floral niche. The notable plasticity associated with disaccharide utilization in the W/S clade is discussed in the context of ecological implications and energy metabolism. IMPORTANCE Microbes usually have flexible metabolic capabilities and are able to use different compounds to meet their needs. The yeasts belonging to the Wickerhamiella and Starmerella genera (forming the so-called W/S clade) are usually found in flowers or insects that visit flowers and are known for having acquired many genes from bacteria by a process called horizontal gene transfer. One such gene, dubbed SUC2, is used to assimilate sucrose, which is one of the most abundant sugars in floral nectar. Here, we show that different lineages within the W/S clade used different solutions for sucrose utilization that dispensed SUC2 and differed in their energy requirements, in their capacity to scavenge small amounts of sucrose from the environment, and in the potential for sharing this resource with other microbial species. We posit that this plasticity is possibly dictated by adaptation to the specific requirements of each species.

The Evolution of the Satratoxin and Atranone Gene Clusters of Stachybotrys chartarum

Stachybotrys chartarum is frequently isolated from damp building materials or improperly stored animal forage. Human and animal exposure to the secondary metabolites of this mold is linked to severe health effects. The mutually exclusive production of either satratoxins or atranones defines the chemotypes A and S. Based upon the genes (satratoxin cluster, SC1-3, sat or atranone cluster, AC1, atr) that are supposed to be essential for satratoxin and atranone production, S. chartarum can furthermore be divided into three genotypes: the S-type possessing all sat- but no atr-genes, the A-type lacking the sat- but harboring all atr-genes, and the H-type having only certain sat- and all atr-genes. We analyzed the above-mentioned gene clusters and their flanking regions to shed light on the evolutionary relationship. Furthermore, we performed a deep re-sequencing and LC-MS/MS (Liquid chromatography-mass spectrometry) analysis. We propose a first model for the evolution of the S. chartarum genotypes. We assume that genotype H represents the most ancient form. A loss of the AC1 and the concomitant acquisition of the SC2 led to the emergence of the genotype S. According to our model, the genotype H also developed towards genotype A, a process that was accompanied by a loss of SC1 and SC3.

Positively selected genes in the hoary bat (Lasiurus cinereus) lineage: prominence of thymus expression, immune and metabolic function, and regions of ancient synteny

Background

Bats of the genus Lasiurus occur throughout the Americas and have diversified into at least 20 species among three subgenera. The hoary bat (Lasiurus cinereus) is highly migratory and ranges farther across North America than any other wild mammal. Despite the ecological importance of this species as a major insect predator, and the particular susceptibility of lasiurine bats to wind turbine strikes, our understanding of hoary bat ecology, physiology, and behavior remains poor.

Methods

To better understand adaptive evolution in this lineage, we used whole-genome sequencing to identify protein-coding sequence and explore signatures of positive selection. Gene models were predicted with Maker and compared to seven well-annotated and phylogenetically representative species. Evolutionary rate analysis was performed with PAML.

Results

Of 9,447 single-copy orthologous groups that met evaluation criteria, 150 genes had a significant excess of nonsynonymous substitutions along the L. cinereus branch (P < 0.001 after manual review of alignments). Selected genes as a group had biased expression, most strongly in thymus tissue. We identified 23 selected genes with reported immune functions as well as a divergent paralog of Steep1 within suborder Yangochiroptera. Seventeen genes had roles in lipid and glucose metabolic pathways, partially overlapping with 15 mitochondrion-associated genes; these adaptations may reflect the metabolic challenges of hibernation, long-distance migration, and seasonal variation in prey abundance. The genomic distribution of positively selected genes differed significantly from background expectation by discrete Kolmogorov-Smirnov test (P < 0.001). Remarkably, the top three physical clusters all coincided with islands of conserved synteny predating Mammalia, the largest of which shares synteny with the human cat-eye critical region (CECR) on 22q11. This observation coupled with the expansion of a novel Tbx1-like gene family may indicate evolutionary innovation during pharyngeal arch development: both the CECR and Tbx1 cause dosage-dependent congenital abnormalities in thymus, heart, and head, and craniodysmorphy is associated with human orthologs of other positively selected genes as well.