Genome expansion and lineage-specific genetic innovations in the forest pathogenic fungi Armillaria

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Composition and diversity of soil microbial communities change by introducing Phallus impudicus into a Gastrodia elata Bl.-based soil

BACKGROUND

The Gastrodia elata Bl. is an orchid, and its growth demands the presence of Armillaria species. The strong competitiveness of Armillaria species has always been a concern of major threat to other soil organisms, thus disrupting the equilibrium of soil biodiversity. Introducing other species to where G. elata was cultivated, could possibly alleviate the problems associated with the disequilibrium of soil microenvironment; however, their impacts on the soil microbial communities and the underlying mechanisms remain unclear. To reveal the changes of microbial groups associated with soil chemical properties responding to different cultivation species, the chemical property measurements coupled with the next-generation pyrosequencing analyses were applied with soil samples collected from fallow land, cultivation of G. elata and Phallus impudicus, respectively.

RESULTS

The cultivation of G. elata induced significant increases (p < 0.05) in soil pH and NO3-N content compared with fallow land, whereas subsequent cultivation of P. impudicus reversed these G. elata-induced increases and was also found to significantly increase (p < 0.05) the content of soil NH4+-N and AP. The alpha diversities of soil microbial communities were significantly increased (p < 0.01) by cultivation of G. elata and P. impudicus as indicated with Chao1 estimator and Shannon index. The structure and composition of soil microbial communities differed responding to different cultivation species. In particular, the relative abundances of Bacillus, norank_o_Gaiellales, Mortierella and unclassified_k_Fungi were significantly increased (p < 0.05), while the abundances of potentially beneficial genera such as Acidibacter, Acidothermus, Cryptococcus, and Penicillium etc., were significantly decreased (p < 0.05) by cultivation of G. elata. It's interesting to find that cultivation of P. impudicus increased the abundances of these genera that G. elata decreased before, which contributed to the difference of composition and structure. The results of CCA and heatmap indicated that the changes of soil microbial communities had strong correlations with soil nutrients. Specifically, among 28 genera presented, 50% and 42.9% demonstrated significant correlations with soil pH and NO3-N in response to cultivation of G. elata and P. impudicus.

CONCLUSIONS

Our findings suggested that the cultivation of P. impudicus might have potential benefits as result of affecting soil microorganisms coupled with changes in soil nutrient profile.

Snowball: a novel gene family required for developmental patterning of fruiting bodies of mushroom-forming fungi (Agaricomycetes)

The morphogenesis of sexual fruiting bodies of fungi is a complex process determined by a genetically encoded program. Fruiting bodies reached the highest complexity levels in the Agaricomycetes; yet, the underlying genetics is currently poorly known. In this work, we functionally characterized a highly conserved gene termed snb1, whose expression level increases rapidly during fruiting body initiation. According to phylogenetic analyses, orthologs of snb1 are present in almost all agaricomycetes and may represent a novel conserved gene family that plays a substantial role in fruiting body development. We disrupted snb1 using CRISPR/Cas9 in the agaricomycete model organism Coprinopsis cinerea. snb1 deletion mutants formed unique, snowball-shaped, rudimentary fruiting bodies that could not differentiate caps, stipes, and lamellae. We took advantage of this phenotype to study fruiting body differentiation using RNA-Seq analyses. This revealed differentially regulated genes and gene families that, based on wild-type RNA-Seq data, were upregulated early during development and showed tissue-specific expression, suggesting a potential role in differentiation. Taken together, the novel gene family of snb1 and the differentially expressed genes in the snb1 mutants provide valuable insights into the complex mechanisms underlying developmental patterning in the Agaricomycetes.

IMPORTANCE

Fruiting bodies of mushroom-forming fungi (Agaricomycetes) are complex multicellular structures, with a spatially and temporally integrated developmental program that is, however, currently poorly known. In this study, we present a novel, conserved gene family, Snowball (snb), termed after the unique, differentiation-less fruiting body morphology of snb1 knockout strains in the model mushroom Coprinopsis cinerea. snb is a gene of unknown function that is highly conserved among agaricomycetes and encodes a protein of unknown function. A comparative transcriptomic analysis of the early developmental stages of differentiated wild-type and non-differentiated mutant fruiting bodies revealed conserved differentially expressed genes which may be related to tissue differentiation and developmental patterning fruiting body development.

Genome Sequencing and Analysis of Nigrospora oryzae, a Rice Leaf Disease Fungus

Nigrospora oryzae is one of several fungal pathogens known to cause brown streaks, leaf spots, and latent infections in rice. In this study, the entire 42.09-Mb genome of N. oryzae was sequenced at a depth of 169× using the Oxford Nanopore Technologies platform. The draft genome sequence was comprised of 26 scaffolds, possessed an average GC content of 58.83%, and contained a total of 10,688 protein-coding genes. Analysis of the complete genome sequence revealed that CAZyme-encoding genes account for 6.11% of all identified genes and that numerous transcription factors (TFs) associated with diverse biological processes belong predominantly to Zn-clus (22.20%) and C2H2 (10.59%) fungal TF classes. In addition, genes encoding 126 transport proteins and 3307 pathogen-host interaction proteins were identified. Comparative analysis of the previously reported N. oryzae reference strain GZL1 genome and the genome of a representative strain ZQ1 obtained here revealed 9722 colinear genes. Collectively, these findings provide valuable insights into N. oryzae genetic mechanisms and phenotypic characteristics.

The role of symbiotic fungi in the life cycle of Gastrodia elata Blume (Orchidaceae): a comprehensive review

Gastrodia elata Blume, a fully mycoheterotrophic perennial plant of the family Orchidaceae, is a traditional Chinese herb with medicinal and edible value. Interestingly, G. elata requires symbiotic relationships with Mycena and Armillaria strains for seed germination and plant growth, respectively. However, there is no comprehensive summary of the symbiotic mechanism between fungi and G. elata. Here, the colonization and digestion of hyphae, the bidirectional exchange of nutrients, the adaptation of fungi and G. elata to symbiosis, and the role of microorganisms and secondary metabolites in the symbiotic relationship between fungi and G. elata are summarized. We comprehensively and deeply analyzed the mechanism of symbiosis between G. elata and fungi from three perspectives: morphology, nutrition, and molecules. The aim of this review was to enrich the understanding of the mutualistic symbiosis mechanisms between plants and fungi and lay a theoretical foundation for the ecological cultivation of G. elata.

Enzymatic machinery of wood-inhabiting fungi that degrade temperate tree species

Deadwood provides habitat for fungi and serves diverse ecological functions in forests. We already have profound knowledge of fungal assembly processes, physiological and enzymatic activities, and resulting physico-chemical changes during deadwood decay. However, in situ detection and identification methods, fungal origins, and a mechanistic understanding of the main lignocellulolytic enzymes are lacking. This study used metaproteomics to detect the main extracellular lignocellulolytic enzymes in 12 tree species in a temperate forest that have decomposed for 8 ½ years. Mainly white-rot (and few brown-rot) Basidiomycota were identified as the main wood decomposers, with Armillaria as the dominant genus; additionally, several soft-rot xylariaceous Ascomycota were identified. The key enzymes involved in lignocellulolysis included manganese peroxidase, peroxide-producing alcohol oxidases, laccase, diverse glycoside hydrolases (cellulase, glucosidase, xylanase), esterases, and lytic polysaccharide monooxygenases. The fungal community and enzyme composition differed among the 12 tree species. Ascomycota species were more prevalent in angiosperm logs than in gymnosperm logs. Regarding lignocellulolysis as a function, the extracellular enzyme toolbox acted simultaneously and was interrelated (e.g. peroxidases and peroxide-producing enzymes were strongly correlated), highly functionally redundant, and present in all logs. In summary, our in situ study provides comprehensive and detailed insight into the enzymatic machinery of wood-inhabiting fungi in temperate tree species. These findings will allow us to relate changes in environmental factors to lignocellulolysis as an ecosystem function in the future.

Fungal Alcohol Dehydrogenases: Physiological Function, Molecular Properties, Regulation of Their Production, and Biotechnological Potential

Fungal alcohol dehydrogenases (ADHs) participate in growth under aerobic or anaerobic conditions, morphogenetic processes, and pathogenesis of diverse fungal genera. These processes are associated with metabolic operation routes related to alcohol, aldehyde, and acid production. The number of ADH enzymes, their metabolic roles, and their functions vary within fungal species. The most studied ADHs are associated with ethanol metabolism, either as fermentative enzymes involved in the production of this alcohol or as oxidative enzymes necessary for the use of ethanol as a carbon source; other enzymes participate in survival under microaerobic conditions. The fast generation of data using genome sequencing provides an excellent opportunity to determine a correlation between the number of ADHs and fungal lifestyle. Therefore, this review aims to summarize the latest knowledge about the importance of ADH enzymes in the physiology and metabolism of fungal cells, as well as their structure, regulation, evolutionary relationships, and biotechnological potential.

Armillaria altimontana in North America: Biology and Ecology

Armillaria altimontana is a fungus (Basidiomycota, Agaricomycetes, Agaricales, and Physalacriaceae) that is generally considered as a weak/opportunistic pathogen or saprophyte on many tree hosts. It widely occurs across the northwestern USA to southern British Columbia, Canada, but relatively little is known about its ecological role in the diverse forest ecosystems where it occurs. This review summarizes the biology and ecology of A. altimontana, including its identification, life cycle, distribution, host associations, and bioclimatic models under climate change.

Vertical and horizontal gene transfer shaped plant colonization and biomass degradation in the fungal genus Armillaria

The fungal genus Armillaria contains necrotrophic pathogens and some of the largest terrestrial organisms that cause tremendous losses in diverse ecosystems, yet how they evolved pathogenicity in a clade of dominantly non-pathogenic wood degraders remains elusive. Here we show that Armillaria species, in addition to gene duplications and de novo gene origins, acquired at least 1,025 genes via 124 horizontal gene transfer events, primarily from Ascomycota. Horizontal gene transfer might have affected plant biomass degrading and virulence abilities of Armillaria, and provides an explanation for their unusual, soft rot-like wood decay strategy. Combined multi-species expression data revealed extensive regulation of horizontally acquired and wood-decay related genes, putative virulence factors and two novel conserved pathogenicity-induced small secreted proteins, which induced necrosis in planta. Overall, this study details how evolution knitted together horizontally and vertically inherited genes in complex adaptive traits of plant biomass degradation and pathogenicity in important fungal pathogens.

A practical approach to genome assembly and annotation of Basidiomycota using the example of Armillaria

Technological advancements in genome sequencing, assembly and annotation platforms and algorithms that resulted in several genomic studies have created an opportunity to further our understanding of the biology of phytopathogens, including Armillaria species. Most Armillaria species are facultative necrotrophs that cause root- and stem-rot, usually on woody plants, significantly impacting agriculture and forestry worldwide. Genome sequencing, assembly and annotation in terms of samples used and methods applied in Armillaria genome projects are evaluated in this review. Infographic guidelines and a database of resources to facilitate future Armillaria genome projects were developed. Knowledge gained from genomic studies of Armillaria species is summarized and prospects for further research are provided. This guide can be applied to other diploid and dikaryotic fungal genomics.

The mycoremediation potential of the armillarioids: a comparative genomics analysis

Genes involved in mycoremediation were identified by comparative genomics analysis in 10 armillarioid species and selected groups of white-rot Basidiomycota (14) and soft-rot Ascomycota (12) species to confine the distinctive bioremediation capabilities of the armillarioids. The genomes were explored using phylogenetic principal component analysis (pPCA), searching for genes already documented in a biocatalysis/biodegradation database. The results underlined a distinct, increased potential of aromatics-degrading genes/enzymes in armillarioids, with particular emphasis on a high copy number and diverse spectrum of benzoate 4-monooxygenase [EC:1.14.14.92] homologs. In addition, other enzymes involved in the degradation of various monocyclic aromatics were more abundant in the armillarioids than in the other white-rot basidiomycetes, and enzymes involved in the degradation of polycyclic aromatic hydrocarbons (PAHs) were more prevailing in armillarioids and other white-rot species than in soft-rot Ascomycetes. Transcriptome profiling of A. ostoyae and A. borealis isolates confirmed that several genes involved in the degradation of benzoates and other monocyclic aromatics were distinctively expressed in the wood-invading fungal mycelia. Data were consistent with armillarioid species offering a more powerful potential in degrading aromatics. Our results provide a reliable, practical solution for screening the likely fungal candidates for their full biodegradation potential, applicability, and possible specialization based on their genomics data.

Insights into the Ecological Diversification of the Hymenochaetales based on Comparative Genomics and Phylogenomics With an Emphasis on Coltricia

To elucidate the genomic traits of ecological diversification in the Hymenochaetales, we sequenced 15 new genomes, with attention to ectomycorrhizal (EcM) Coltricia species. Together with published data, 32 genomes, including 31 Hymenochaetales and one outgroup, were comparatively analyzed in total. Compared with those of parasitic and saprophytic members, EcM species have significantly reduced number of plant cell wall degrading enzyme genes, and expanded transposable elements, genome sizes, small secreted proteins, and secreted proteases. EcM species still retain some of secreted carbohydrate-active enzymes (CAZymes) and have lost the key secreted CAZymes to degrade lignin and cellulose, while possess a strong capacity to degrade a microbial cell wall containing chitin and peptidoglycan. There were no significant differences in secreted CAZymes between fungi growing on gymnosperms and angiosperms, suggesting that the secreted CAZymes in the Hymenochaetales evolved before differentiation of host trees into gymnosperms and angiosperms. Nevertheless, parasitic and saprophytic species of the Hymenochaetales are very similar in many genome features, which reflect their close phylogenetic relationships both being white rot fungi. Phylogenomic and molecular clock analyses showed that the EcM genus Coltricia formed a clade located at the base of the Hymenochaetaceae and divergence time later than saprophytic species. And Coltricia remains one to two genes of AA2 family. These indicate that the ancestors of Coltricia appear to have originated from saprophytic ancestor with the ability to cause a white rot. This study provides new genomic data for EcM species and insights into the ecological diversification within the Hymenochaetales based on comparative genomics and phylogenomics analyses.

Hallmarks of Comparative Transcriptome between Rhizomorphs and Hyphae of Armillaria sp. 541 Participating in Fungal Symbiosis with Emphasis on LysM Domains

Armillaria sp. 541, a genus of root-infecting fungi, forms a symbiosis with traditional Chinese medicine Gastrodia elata (Orchid) and Polyporus umbellatus via extensive networks of durable rhizomorphs. It is not clear the hallmarks of comparative transcriptome between the rhizomorphs and hyphae of Armillaria sp. 541. In the present study, transcriptomic analysis of Armillaria sp. 541 identified 475 differentially expressed genes (DEGs) between Armillaria rhizomorphs (AR) and hyphae (AH). Of them, 285 genes were upregulated and 190 were downregulated. Bioinformatics analyses and tests demonstrated DEGs involved in oxidoreductase activity and peptidoglycan binding were significantly enriched in this process when rhizomorph formed from hyphae. We accordingly obtained 14 gene-encoding proteins containing the LysM domain, and further consensus pattern and phylogenetic analysis indicated that their amino acid sequences were conserved and their biological functions may be peptidoglycan binding for recognition between the fungus and host. Among these genes, one, named Armillaria LysM domain recognition gene (aLDRG), was expressed significantly when rhizomorphs were differentiated from hyphae. It was located in the cortical cells of the rhizomorph by in situ hybridization. Furthermore, biolayer interferometry (BLI) assay demonstrated that aLDRG can bind specifically to chitin oligosaccharide of the fungal cell wall, including N,N',N″-Triacetylchitotriose (CO3) and N,N',N″,N'″,N″″-Pentaacetylchitopentaose (CO5). Therefore, we deduced that Armillaria sp. 541 expressed higher levels of LysM protein aLDRG for better binding of oligosaccharide after rhizomorphs were generated. This study provides functional genes for further studies on the interaction between Armillaria sp. 541 and its host.

BASIDIN as a New Protein Effector of the Phytopathogen Causing Witche's Broom Disease in Cocoa

The fungus Moniliophthora perniciosa secretes protein effectors that manipulate the physiology of the host plant, but few effectors of this fungus have had their functions confirmed. We performed functional characterization of a promising candidate effector of M. perniciosa. The inoculation of rBASIDIN at 4 µmol L^-1 in the mesophyll of leaflets of Solanum lycopersicum caused symptoms of shriveling within 6 h without the presence of necrosis. However, when sprayed on the plant at a concentration of 11 µmol L^-1, it caused wilting symptoms only 2 h after application, followed by necrosis and cell death at 48 h. rBASIDIN applied to Theobroma cacao leaves at the same concentration caused milder symptoms. rBASIDIN caused hydrogen peroxide production in leaf tissue, damaging the leaf membrane and negatively affecting the photosynthetic rate of Solanum lycopersicum plants. Phylogenetic analysis indicated that BASIDIN has orthologs in other phytopathogenic basidiomycetes. Analysis of the transcripts revealed that BASIDIN and its orthologs are expressed in different fungal species, suggesting that this protein is differentially regulated in these basidiomycetes. Therefore, the results of applying BASIDIN allow the inference that it is an effector of the fungus M. perniciosa, with a strong potential to interfere in the defense system of the host plant.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Dual RNA-Seq Profiling Unveils Mycoparasitic Activities of Trichoderma atroviride against Haploid Armillaria ostoyae in Antagonistic Interaction Assays

Armillaria ostoyae, a species among the destructive forest pathogens from the genus Armillaria, causes root rot disease on woody plants worldwide. Efficient control measures to limit the growth and impact of this severe underground pathogen are under investigation. In a previous study, a new soilborne fungal isolate, Trichoderma atroviride SZMC 24276 (TA), exhibited high antagonistic efficacy, which suggested that it could be utilized as a biocontrol agent. The dual culture assay results indicated that the haploid A. ostoyae-derivative SZMC 23085 (AO) (C18/9) is highly susceptible to the mycelial invasion of TA. In the present study, we analyzed the transcriptome of AO and that of TA in in vitro dual culture assays to test the molecular arsenal of Trichoderma antagonism and the defense mechanisms of Armillaria. We conducted time-course analysis and functional annotation and analyzed enriched pathways and differentially expressed genes including biocontrol-related candidate genes from TA and defense-related candidate genes from AO. The results indicated that TA deployed several biocontrol mechanisms when confronted with AO. In response, AO initiated multiple defense mechanisms to protect against the fungal attack. To our knowledge, the present study offers the first transcriptome analysis of a biocontrol fungus attacking AO. Overall, this study provides insights that aid the further exploration of plant pathogen-biocontrol agent interaction mechanisms. IMPORTANCE Armillaria species can survive for decades in the soil on dead woody debris, develop rapidly under favorable conditions, and harmfully infect newly planted forests. Our previous study found Trichoderma atroviride to be highly effective in controlling Armillaria growth; therefore, our current work explored the molecular mechanisms that might play a key role in Trichoderma-Armillaria interactions. Direct confrontation assays combined with time course-based dual transcriptome analysis provided a reliable system for uncovering the interactive molecular dynamics between the fungal plant pathogen and its mycoparasitic partner. Furthermore, using a haploid Armillaria isolate allowed us to survey the deadly prey-invading activities of the mycoparasite and the ultimate defensive strategies of its prey. Our current study provides detailed insights into the essential genes and mechanisms involved in Armillaria defense against Trichoderma and the genes potentially involved in the efficiency of Trichoderma to control Armillaria. In addition, using a sensitive haploid Armillaria strain (C18/9), with its complete genome data already available, also offers the opportunity to test possible variable molecular responses of Armillaria ostoyae toward diverse Trichoderma isolates with various biocontrol abilities. Initial molecular tests of the dual interactions may soon help to develop a targeted biocontrol intervention with mycoparasites against plant pathogens.

Genomic and Secretomic Analyses of the Newly Isolated Fungus Perenniporia fraxinea SS3 Identified CAZymes Potentially Related to a Serious Pathogenesis of Hardwood Trees

Perenniporia fraxinea can colonize living trees and cause severe damage to standing hardwoods by secreting a number of carbohydrate-activate enzymes (CAZymes), unlike other well-studied Polyporales. However, significant knowledge gaps exist in understanding the detailed mechanisms for this hardwood-pathogenic fungus. To address this issue, five monokaryotic P. fraxinea strains, SS1 to SS5, were isolated from the tree species Robinia pseudoacacia, and high polysaccharide-degrading activities and the fastest growth were found for P. fraxinea SS3 among the isolates. The whole genome of P. fraxinea SS3 was sequenced, and its unique CAZyme potential for tree pathogenicity was determined in comparison to the genomes of other nonpathogenic Polyporales. These CAZyme features are well conserved in a distantly related tree pathogen, Heterobasidion annosum. Furthermore, the carbon source-dependent CAZyme secretions of P. fraxinea SS3 and a nonpathogenic and strong white-rot Polyporales member, Phanerochaete chrysosporium RP78, were compared by activity measurements and proteomic analyses. As seen in the genome comparisons, P. fraxinea SS3 exhibited higher pectin-degrading activities and higher laccase activities than P. chrysosporium RP78, which were attributed to the secretion of abundant glycoside hydrolase family 28 (GH28) pectinases and auxiliary activity family 1_1 (AA1_1) laccases, respectively. These enzymes are possibly related to fungal invasion into the tree lumens and the detoxification of tree defense substances. Additionally, P. fraxinea SS3 showed secondary cell wall degradation capabilities at the same level as that of P. chrysosporium RP78. Overall, this study suggested mechanisms for how this fungus can attack the cell walls of living trees as a serious pathogen and differs from other nonpathogenic white-rot fungi. IMPORTANCE Many studies have been done to understand the mechanisms underlying the degradation of plant cell walls of dead trees by wood decay fungi. However, little is known about how some of these fungi weaken living trees as pathogens. P. fraxinea belongs to the Polyporales, a group of strong wood decayers, and is known to aggressively attack and fell standing hardwood trees all over the world. Here, we report CAZymes potentially related to plant cell wall degradation and pathogenesis factors in a newly isolated fungus, P. fraxinea SS3, by genome sequencing in conjunction with comparative genomic and secretomic analyses. The present study provides insights into the mechanisms of the degradation of standing hardwood trees by the tree pathogen, which will contribute to the prevention of this serious tree disease.

Whole genome sequencing and analysis of Armillaria gallica Jzi34 symbiotic with Gastrodia elata

BACKGROUND

Armillaria species are plant pathogens, but a few Armillaria species can establish a symbiotic relationship with Gastrodia elata, a rootless and leafless orchid, that is used as a Chinese herbal medicine. Armillaria is a nutrient source for the growth of G. elata. However, there are few reports on the molecular mechanism of symbiosis between Armillaria species and G. elata. The genome sequencing and analysis of Armillaria symbiotic with G. elata would provide genomic information for further studying the molecular mechanism of symbiosis.

RESULTS

The de novo genome assembly was performed with the PacBio Sequel platform and Illumina NovaSeq PE150 for the A. gallica Jzi34 strain, which was symbiotic with G. elata. Its genome assembly contained ~ 79.9 Mbp and consisted of 60 contigs with an N50 of 2,535,910 bp. There were only 4.1% repetitive sequences in the genome assembly. Functional annotation analysis revealed a total of 16,280 protein coding genes. Compared with the other five genomes of Armillaria, the carbohydrate enzyme gene family of the genome was significantly contracted, while it had the largest set of glycosyl transferase (GT) genes. It also had an expansion of auxiliary activity enzymes AA3-2 gene subfamily and cytochrome P450 genes. The synteny analysis result of P450 genes reveals that the evolutionary relationship of P450 proteins between A. gallica Jzi34 and other four Armillaria was complex.

CONCLUSIONS

These characteristics may be beneficial for establishing a symbiotic relationship with G. elata. These results explore the characteristics of A. gallica Jzi34 from a genomic perspective and provide an important genomic resource for further detailed study of Armillaria. This will help to further study the symbiotic mechanism between A. gallica and G. elata.

Genomic and Transcriptomic Approaches Provide a Predictive Framework for Sesquiterpenes Biosynthesis in Desarmillaria tabescens CPCC 401429

Terpenoids constitute a structurally diverse class of secondary metabolites with wide applications in the pharmaceutical, fragrance and flavor industries. Desarmillaria tabescens CPCC 401429 is a basidiomycetous mushroom that could produce anti-tumor melleolides. To date, no studies have been conducted to thoroughly investigate the sesquiterpenes biosynthetic potential in Desarmillaria or related genus. This study aims to unravel the phylogeny, terpenome, and functional characterization of unique sesquiterpene biosynthetic genes of the strain CPCC 401429. Herein, we report the genome of the fungus containing 15,145 protein-encoding genes. MLST-based phylogeny and comparative genomic analyses shed light on the precise reclassification of D. tabescens suggesting that it belongs to the genus Desarmillaria. Gene ontology enrichment and pathway analyses uncover the hidden capacity for producing polyketides and terpenoids. Genome mining directed predictive framework reveals a diverse network of sesquiterpene synthases (STSs). Among twelve putative STSs encoded in the genome, six ones are belonging to the novel minor group: diverse Clade IV. In addition, RNA-sequencing based transcriptomic profiling revealed differentially expressed genes (DEGs) of the fungus CPCC 401429 in three different fermentation conditions, that of which enable us to identify noteworthy genes exemplified as STSs coding genes. Among the ten sesquiterpene biosynthetic DEGs, two genes including DtSTS9 and DtSTS10 were selected for functional characterization. Yeast cells expressing DtSTS9 and DtSTS10 could produce diverse sesquiterpene compounds, reinforced that STSs in the group Clade IV might be highly promiscuous producers. This highlights the potential of Desarmillaria in generating novel terpenoids. To summarize, our analyses will facilitate our understanding of phylogeny, STSs diversity and functional significance of Desarmillaria species. These results will encourage the scientific community for further research on uncharacterized STSs of Basidiomycota phylum, biological functions, and potential application of this vast source of secondary metabolites.

Haplotype-phased and chromosome-level genome assembly of Puccinia polysora, a giga-scale fungal pathogen causing southern corn rust

Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P. polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P. polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P. polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P. polysora.

Nonribosomal peptide synthetase gene clusters and characteristics of predicted NRPS-dependent siderophore synthetases in Armillaria and other species in the Physalacriaceae

Fungal secondary metabolites are often pathogenicity or virulence factors synthesized by genes contained in secondary metabolite gene clusters (SMGCs). Nonribosomal polypeptide synthetase (NRPS) clusters are SMGCs which produce peptides such as siderophores, the high affinity ferric iron chelating compounds required for iron uptake under aerobic conditions. Armillaria spp. are mostly facultative necrotrophs of woody plants. NRPS-dependent siderophore synthetase (NDSS) clusters of Armillaria spp. and selected Physalacriaceae were investigated using a comparative genomics approach. Siderophore biosynthesis by strains of selected Armillaria spp. was evaluated using CAS and split-CAS assays. At least one NRPS cluster and other clusters were detected in the genomes studied. No correlation was observed between the number and types of SMGCs and reported pathogenicity of the species studied. The genomes contained one NDSS cluster each. All NDSSs were multi-modular with the domain architecture (ATC)₃(TC)₂. NDSS clusters of the Armillaria spp. showed a high degree of microsynteny. In the genomes of Desarmillaria spp. and Guyanagaster necrorhizus, NDSS clusters were more syntenic with NDSS clusters of Armillaria spp. than to those of the other Physalacriaceae species studied. Three A-domain orthologous groups were identified in the NDSSs, and atypical Stachelhaus codes were predicted for the A3 orthologous group. In vitro biosynthesis of mainly hydroxamate and some catecholate siderophores was observed. Hence, Armillaria spp. generally contain one highly conserved, NDSS cluster although some interspecific variations in the products of these clusters is expected. Results from this study lays the groundwork for future studies to elucidate the molecular biology of fungal phyto-pathogenicity.

Genomic Comparisons of Two Armillaria Species with Different Ecological Behaviors and Their Associated Soil Microbial Communities

Armillaria species show considerable variation in ecological roles and virulence, from mycorrhizae and saprophytes to important root pathogens of trees and horticultural crops. We studied two Armillaria species that can be found in coniferous forests of northwestern USA and southwestern Canada. Armillaria altimontana not only is considered as a weak, opportunistic pathogen of coniferous trees, but it also appears to exhibit in situ biological control against A. solidipes, formerly North American A. ostoyae, which is considered a virulent pathogen of coniferous trees. Here, we describe their genome assemblies and present a functional annotation of the predicted genes and proteins for the two Armillaria species that exhibit contrasting ecological roles. In addition, the soil microbial communities were examined in association with the two Armillaria species within a 45-year-old plantation of western white pine (Pinus monticola) in northern Idaho, USA, where A. altimontana was associated with improved tree growth and survival, while A. solidipes was associated with reduced growth and survival. The results from this study reveal a high similarity between the genomes of the beneficial/non-pathogenic A. altimontana and pathogenic A. solidipes; however, many relatively small differences in gene content were identified that could contribute to differences in ecological lifestyles and interactions with woody hosts and soil microbial communities.

Unlocking the magic in mycelium: Using synthetic biology to optimize filamentous fungi for biomanufacturing and sustainability

Filamentous fungi drive carbon and nutrient cycling across our global ecosystems, through its interactions with growing and decaying flora and their constituent microbiomes. The remarkable metabolic diversity, secretion ability, and fiber-like mycelial structure that have evolved in filamentous fungi have been increasingly exploited in commercial operations. The industrial potential of mycelial fermentation ranges from the discovery and bioproduction of enzymes and bioactive compounds, the decarbonization of food and material production, to environmental remediation and enhanced agricultural production. Despite its fundamental impact in ecology and biotechnology, molds and mushrooms have not, to-date, significantly intersected with synthetic biology in ways comparable to other industrial cell factories (e.g. Escherichia coli,Saccharomyces cerevisiae, and Komagataella phaffii). In this review, we summarize a suite of synthetic biology and computational tools for the mining, engineering and optimization of filamentous fungi as a bioproduction chassis. A combination of methods across genetic engineering, mutagenesis, experimental evolution, and computational modeling can be used to address strain development bottlenecks in established and emerging industries. These include slow mycelium growth rate, low production yields, non-optimal growth in alternative feedstocks, and difficulties in downstream purification. In the scope of biomanufacturing, we then detail previous efforts in improving key bottlenecks by targeting protein processing and secretion pathways, hyphae morphogenesis, and transcriptional control. Bringing synthetic biology practices into the hidden world of molds and mushrooms will serve to expand the limited panel of host organisms that allow for commercially-feasible and environmentally-sustainable bioproduction of enzymes, chemicals, therapeutics, foods, and materials of the future.

Comparative transcriptome analysis of Armillaria gallica 012m in response to ethephon treatment

Background

Gastrodia elata, known as a rootless, leafless, achlorophyllous and fully mycoheterotrophic orchid, needs to establish symbionts with particular Armillaria species to acquire nutrition and energy. Previous research findings had approved that ethylene (ET) played an important role in plant-fungi interaction and some receptors of ET had been discovered in microorganisms. However, the molecular mechanisms underlying the role of ET in the interaction between G. elata and Armillaria species remain unknown.

Methods

Exiguous ethephon (ETH) was added to agar and liquid media to observe the morphological features of mycelium and count the biomass respectively. Mycelium cultured in liquid media with exiguous ETH (0.1 ppm, 2.0 ppm, 5.0 ppm) were chosen to perform whole-transcriptome profiling through the RNA-seq technology (Illumina NGS sequencing). The DEGs of growth-related genes and candidate ET receptor domains were predicted on SMART.

Results

ETH-0.1 ppm and ETH-2 ppm could significantly improve the mycelium growth of A. gallica 012m, while ETH-5 ppm inhibited the mycelium growth in both solid and liquid media. The number of up-regulated or down-regulated genes increased along with the concentrations of ETH. The growth of mycelia might benefit from the up-regulated expression of Pyr_redox (Pyridine nucleotide-disulphide oxidoreductase), GAL4 (C6 zinc finger) and HMG (High Mobility Group) genes in the ETH-0.1 ppm and ETH-2 ppm. Therefore, the growth of mycelia might be impaired by the down-regulated expression of ZnF_C2H2 and ribosomal protein S4 proteins in the ETH-5 ppm. Seven ET receptor domains were predicted in A. gallica 012m. Based on cluster analysis and comparative studies of proteins, the putative ETH receptor domains of A. gallica 012m have a higher homologous correlation with fungi.

Conclusions

The responses of A. gallica 012m to ETH had a concentration effect similar to the plants’ responses to ET. Therefore, the number of up-regulated or down-regulated genes are increased along with the concentrations of ETH. Seven ET receptor protein domains were predicted in the genome and transcriptome of A. gallica 012m. We speculate that ETH receptors exist in A. gallica 012m and ethylene might play an important role in the plant-fungi interaction.

A facultative ectomycorrhizal association is triggered by organic nitrogen

Increasing nitrogen (N) deposition often tends to negatively impact the functions of belowground ectomycorrhizal networks, although the exact molecular mechanisms underlying this trait are still unclear. Here, we assess how the root-associated fungus Clitopilus hobsonii establishes an ectomycorrhiza-like association with its host tree Populus tomentosa and how this interaction is favored by organic N over mineral N. The establishment of a functional symbiosis in the presence of organic N promotes plant growth and the transfer of ¹⁵N from the fungus to above ground plant tissues. Genomic traits and in planta transcriptional signatures suggest that C. hobsonii may have a dual lifestyle with saprotrophic and mutualistic traits. For example, several genes involved in the digestion of cellulose and hemicellulose are highly expressed during the interaction, whereas the expression of multiple copies of pectin-digesting genes is tightly controlled. Conversely, the nutritional mutualism is dampened in the presence of ammonium (NH₄⁺) or nitrate (NO₃^-). Increasing levels of NH₄⁺ led to a higher expression of pectin-digesting genes and a continuous increase in hydrogen peroxide production in roots, whereas the presence of NO₃^- resulted in toxin production. In summary, our results suggest that C. hobsonii is a facultative ectomycorrhizal fungus. Access to various forms of N acts as an on/off switch for mutualism caused by large-scale fungal physiological remodeling. Furthermore, the abundance of pectin-degrading enzymes with distinct expression patterns during functional divergence after exposure to NH₄⁺ or organic N is likely to be central to the transition from parasitism to mutualism.

Comparative genomic analysis reveals contraction of gene families with putative roles in pathogenesis in the fungal boxwood pathogens Calonectria henricotiae and C. pseudonaviculata

Background

Boxwood blight disease caused by Calonectria henricotiae and C. pseudonaviculata is of ecological and economic significance in cultivated and native ecosystems worldwide. Prior research has focused on understanding the population genetic and genomic diversity of C. henricotiae and C. pseudonaviculata, but gene family evolution in the context of host adaptation, plant pathogenesis, and trophic lifestyle is poorly understood. This study applied bioinformatic and phylogenetic methods to examine gene family evolution in C. henricotiae, C. pseudonaviculata and 22 related fungi in the Nectriaceae that vary in pathogenic and saprobic (apathogenic) lifestyles.

Results

A total of 19,750 gene families were identified in the 24 genomes, of which 422 were rapidly evolving. Among the six Calonectria species, C. henricotiae and C. pseudonaviculata were the only species to experience high levels of rapid contraction of pathogenesis-related gene families (89% and 78%, respectively). In contrast, saprobic species Calonectria multiphialidica and C. naviculata, two of the closest known relatives of C. henricotiae and C. pseudonaviculata, showed rapid expansion of pathogenesis-related gene families.

Conclusions

Our results provide novel insight into gene family evolution within C. henricotiae and C. pseudonaviculata and suggest gene family contraction may have contributed to limited host-range expansion of these pathogens within the plant family Buxaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02035-4.

Research Status and Application Prospects of the Medicinal Mushroom Armillaria mellea

Armillaria is one of the most common diseases underlying chronic root rot in woody plants. Although there is no particularly effective way to prevent it, soil disinfection is a common effective protective measure. However, Armillaria itself has important medicinal value and is a symbiotic fungus in the cultivation of Gastrodia elata and Polyporus umbellatus. Therefore, researching Armillaria is of great practical significance. In this review, the biological characteristics, cultivation methods, chemical components, food and medicinal value and efficacy of Armillaria were all reviewed, and its development and utilization direction were analyzed and discussed.

A Combination of Genome Mining with an OSMAC Approach Facilitates the Discovery of and Contributions to the Biosynthesis of Melleolides from the Basidiomycete Armillaria tabescens

Genome mining revealed that the genomes of basidiomycetes may include a considerable number of biosynthetic gene clusters (BGCs), yet numerous clusters remain unidentified. Herein, we report a combination of genome mining with an OSMAC (one strain, many compounds) approach to characterize the spectrum of melleolides produced by Armillaria tabescens CPCC 401429. Using F1 fermentation medium, the metabolic pathway of the gene cluster mel was successfully upregulated. From the extracts of the wild-type strain, two new melleolides (1 and 2), along with five new orsellinic acid-derived lactams (10-14), were isolated, and their structures were elucidated by LC-HR-ESIMS/MS and 2D-NMR. Several melleolides exhibited moderate anti-carcinoma (A549, NCI-H520, and H1299) effects with IC₅₀ values of 4.0-48.8 μM. RNA-sequencing based transcriptomic profiling broadened our knowledge of the genetic background, regulation, and mechanisms of melleolide biosynthesis. These results may promote downstream metabolic engineering studies of melleolides. Our study demonstrates the approach is effective for discovering new secondary metabolites from Armillaria sp. and will facilitate the mining of the unexploited biosynthetic potential in other basidiomycetes.

SNP4OrphanSpecies: A bioinformatics pipeline to isolate molecular markers for studying genetic diversity of orphan species

Background

For several decades, an increase in disease or pest emergences due to anthropogenic introduction or environmental changes has been recorded. This increase leads to serious threats to the genetic and species diversity of numerous ecosystems. Many of these events involve species with poor or no genomic resources (called here "orphan species"). This lack of resources is a serious limitation to our understanding of the origin of emergent populations, their ability to adapt to new environments and to predict future consequences to biodiversity. Analyses of genetic diversity are an efficient method to obtain this information rapidly, but require available polymorphic genetic markers.

New information

We developed a generic bioinformatics pipeline to rapidly isolate such markers with the goal for the pipeline to be applied in studies of invasive taxa from different taxonomic groups, with a special focus on forest fungal pathogens and insect pests. This pipeline is based on: 1) an automated de novo genome assembly obtained from shotgun whole genome sequencing using paired-end Illumina technology; 2) the isolation of single-copy genes conserved in species related to the studied emergent organisms; 3) primer development for multiplexed short sequences obtained from these conserved genes. Previous studies have shown that intronic regions of these conserved genes generally contain several single nucleotide polymorphisms within species. The pipeline's functionality was evaluated with sequenced genomes of five invasive or expanding pathogen and pest species in Europe (Armillariaostoyae (Romagn.) Herink 1973, Bursaphelenchusxylophilus Steiner & Buhrer 1934, Sphaeropsissapinea (fr.) Dicko & B. Sutton 1980, Erysiphealphitoides (Griffon & Maubl.) U. Braun & S. Takam. 2000, Thaumetopoeapityocampa Denis & Schiffermüller, 1775). We successfully isolated several pools of one hundred short gene regions for each assembled genome, which can be amplified in multiplex. The bioinformatics pipeline is user-friendly and requires little computational resources. This easy-to-set-up and run method for genetic marker identification will be useful for numerous laboratories studying biological invasions, but with limited resources and expertise in bioinformatics.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga’s genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host–pathogen interactions.

Sesquiterpenes with diverse skeletons from histone deacetylase inhibitor modified cultures of the basidiomycete Cyathus stercoreus (Schwein.) De Toni HFG134

Epigenetic modifiers are proved to be effective specialized products-mining tools by rationally regulating the gene expression of fungal biosynthetic pathways. Chemical investigation on the histone deacetylase inhibitor (HDI) vorinostat (also known as SAHA)-modified cultures of the basidiomycete Cyathus stercoreus (Schwein.) De Toni (Nidulariaceae) led to the isolation of nine previously undescribed sesquiterpenes, and four previously described ones. The structures of the nine undescribed compounds were determined by extensive NMR spectroscopic analysis, HRESIMS analysis, as well as ECD and NMR calculations. Notably, the isolated sesquiterpenes are exclusive or overproduced from the epigenetic modified cultures compared to the negative control cultures. Additionally, the skeleton types of the isolated sesquiterpenes include protoilludalane, illudalane, 1,11-seco-protoilludalane, 10,11-seco-illudalane, and 14(11→10)abeo-illudalane. It is noteworthy that the 14(11→10)abeo-illudalane skeleton is reported for the first time. Cystercorodiol A, 4-O-acetylcybrodol, cystercorotone, and cybrodol showed weak inhibitory activity against the bacterium Escherichia coli ATCC25922 with the inhibitory rates 34.7%, 33.0%, 32.3%, and 29.6% at the concentration 200 μM, respectively. This study suggested that epigenetic modifiers are also an effective tool for specialized metabolite-mining in basidiomycetes.

White rot fungal impact on the evolution of simple phenols during decay of silver fir wood by UHPLC-HQOMS

INTRODUCTION

Silver fir (Abies alba Mill.) is one of the most valuable conifer wood species in Europe. Among the main opportunistic pathogens that cause root and butt rot on silver fir are Armillaria ostoyae and Heterobasidion abietinum. Due to the different enzymatic pools of these wood-decay fungi, different strategies in metabolizing the phenols were available.

OBJECTIVE

This work explores the changes in phenolic compounds during silver fir wood degradation.

METHODOLOGY

Phenols were analyzed before and after fungus inoculation in silver fir macerated wood after 2, 4 and 6 months. All samples were analyzed using high-performance liquid chromatography coupled to a hybrid quadrupole-orbitrap mass spectrometer.

RESULTS

Thirteen compounds, including simple phenols, alkylphenyl alcohols, hydroxybenzoketones, hydroxycinnamaldehydes, hydroxybenzaldehydes, hydroxyphenylacetic acids, hydroxycinnamic acids, hydroxybenzoic acids and hydroxycoumarins, were detected. Pyrocatechol, coniferyl alcohol, acetovanillone, vanillin, benzoic acid, 4-hydroxybenzoic acid and vanillic acid contents decreased during the degradation process. Methyl vanillate, ferulic acid and p-coumaric were initially produced and then degraded. Scopoletin was accumulated. Pyrocatechol, acetovanillone and methyl vanillate were found for the first time in both degrading and non-degrading wood of silver fir.

CONCLUSIONS

Despite differences in the enzymatic pool, both fungi caused a significant decrease in the amounts of phenolic compounds with the accumulation of the only scopoletin. Principal component analysis revealed an initial differentiation between the degradation activity of the two fungal species during degradation, but similar phenolic contents at the end of wood degradation.

New Findings on the Biology and Ecology of the Ecuadorian Amazon Fungus Polyporus leprieurii var. yasuniensis

Polyporus leprieurii var. yasuniensis is a prolific wood-decay fungus inhabiting the forest floor of one of the most biodiverse places on earth, the Yasuní National Park in Ecuador. Basidiocarps and aerial rhizomorphs are commonly found growing on woody debris distributed along the floor of this forest ecosystem. Because of the extraordinary abundance of this fungus in the tropical rainforest, we carried out investigations to better understand the biological and ecological aspects contributing to its prolific distribution. Data on growth inhibition in paired competition studies with sixteen fungal isolates exemplifies defense mechanisms used to defend its territory, including pseudosclerotial plates and the development of a melanized rhizomorphic mat. Results of biomass loss on eleven types of tropical wood in microcosm experiments demonstrated the broad decay capacity of the fungus. In and ex situ observations provided information on how long rhizomorphs can prevail in highly competitive ecosystems as well as stressful conditions in the laboratory. Finally, high concentrations of metal ions occur on rhizomorphs as compared to colonized wood. Sequestration of metal ions from the environment by the melanized rhizomorphs may offer protection against competitors. The development of melanized rhizomorphs is key to find and colonize new substrates and resist changing environmental conditions.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Molecular Evolution of Lysine Biosynthesis in Agaricomycetes

As an indispensable essential amino acid in the human body, lysine is extremely rich in edible mushrooms. The α-aminoadipic acid (AAA) pathway is regarded as the biosynthetic pathway of lysine in higher fungal species in Agaricomycetes. However, there is no deep understanding about the molecular evolutionary relationship between lysine biosynthesis and species in Agaricomycetes. Herein, we analyzed the molecular evolution of lysine biosynthesis in Agaricomycetes. The phylogenetic relationships of 93 species in 34 families and nine orders in Agaricomycetes were constructed with six sequences of LSU, SSU, ITS (5.8 S), RPB1, RPB2, and EF1-α datasets, and then the phylogeny of enzymes involved in the AAA pathway were analyzed, especially homocitrate synthase (HCS), α-aminoadipate reductase (AAR), and saccharopine dehydrogenase (SDH). We found that the evolution of the AAA pathway of lysine biosynthesis is consistent with the evolution of species at the order level in Agaricomycetes. The conservation of primary, secondary, predicted tertiary structures, and substrate-binding sites of the enzymes of HCS, AAR, and SDH further exhibited the evolutionary conservation of lysine biosynthesis in Agaricomycetes. Our results provide a better understanding of the evolutionary conservation of the AAA pathway of lysine biosynthesis in Agaricomycetes.

Phylogenetic Relationships, Speciation, and Origin of Armillaria in the Northern Hemisphere: A Lesson Based on rRNA and Elongation Factor 1-Alpha

Armillaria species have a global distribution and play various roles in the natural ecosystems, e.g., pathogens, decomposers, and mycorrhizal associates. However, their taxonomic boundaries, speciation processes, and origin are poorly understood. Here, we used a phylogenetic approach with 358 samplings from Europe, East Asia, and North America to delimit the species boundaries and to discern the evolutionary forces underpinning divergence and evolution. Three species delimitation methods indicated multiple unrecognized phylogenetic species, and biological species recognition did not reflect the natural evolutionary relationships within Armillaria; for instance, biological species of A. mellea and D. tabescens are divergent and cryptic species/lineages exist associated with their geographic distributions in Europe, North America, and East Asia. While the species-rich and divergent Gallica superclade might represent three phylogenetic species (PS I, PS II, and A. nabsnona) that undergo speciation. The PS II contained four lineages with cryptic diversity associated with the geographic distribution. The genus Armillaria likely originated from East Asia around 21.8 Mya in early Miocene when Boreotropical flora (56–33.9 Mya) and the Bering land bridge might have facilitated transcontinental dispersal of Armillaria species. The Gallica superclade arose at 9.1 Mya and the concurrent vicariance events of Bering Strait opening and the uplift of the northern Tibetan plateau might be important factors in driving the lineage divergence.

A Transcriptomic Atlas of the Ectomycorrhizal Fungus Laccaria bicolor

Trees are able to colonize, establish and survive in a wide range of soils through associations with ectomycorrhizal (EcM) fungi. Proper functioning of EcM fungi implies the differentiation of structures within the fungal colony. A symbiotic structure is dedicated to nutrient exchange and the extramatricular mycelium explores soil for nutrients. Eventually, basidiocarps develop to assure last stages of sexual reproduction. The aim of this study is to understand how an EcM fungus uses its gene set to support functional differentiation and development of specialized morphological structures. We examined the transcriptomes of Laccaria bicolor under a series of experimental setups, including the growth with Populus tremula x alba at different developmental stages, basidiocarps and free-living mycelium, under various conditions of N, P and C supply. In particular, N supply induced global transcriptional changes, whereas responses to P supply seemed to be independent from it. Symbiosis development with poplar is characterized by transcriptional waves. Basidiocarp development shares transcriptional signatures with other basidiomycetes. Overlaps in transcriptional responses of L. bicolor hyphae to a host plant and N/C supply next to co-regulation of genes in basidiocarps and mature mycorrhiza were detected. Few genes are induced in a single condition only, but functional and morphological differentiation rather involves fine tuning of larger gene sets. Overall, this transcriptomic atlas builds a reference to study the function and stability of EcM symbiosis in distinct conditions using L. bicolor as a model and indicates both similarities and differences with other ectomycorrhizal fungi, allowing researchers to distinguish conserved processes such as basidiocarp development from nutrient homeostasis.

Evolutionary Morphogenesis of Sexual Fruiting Bodies in Basidiomycota: Toward a New Evo-Devo Synthesis

The development of sexual fruiting bodies is one of the most complex morphogenetic processes in fungi. Mycologists have long been fascinated by the morphological and developmental diversity of fruiting bodies; however, evolutionary developmental biology of fungi still lags significantly behind that of animals or plants. Here, we summarize the current state of knowledge on fruiting bodies of mushroom-forming Basidiomycota, focusing on phylogenetic and developmental biology. Phylogenetic approaches have revealed a complex history of morphological transformations and convergence in fruiting body morphologies. Frequent transformations and convergence is characteristic of fruiting bodies in contrast to animals or plants, where main body plans are highly conserved. At the same time, insights into the genetic bases of fruiting body development have been achieved using forward and reverse genetic approaches in selected model systems. Phylogenetic and developmental studies of fruiting bodies have each yielded major advances, but they have produced largely disjunct bodies of knowledge. An integrative approach, combining phylogenetic, developmental, and functional biology, is needed to achieve a true fungal evolutionary developmental biology (evo-devo) synthesis for fungal fruiting bodies.

Global Distribution and Richness of Armillaria and Related Species Inferred From Public Databases and Amplicon Sequencing Datasets

Armillaria is a globally distributed fungal genus most notably composed of economically important plant pathogens that are found predominantly in forest and agronomic systems. The genus sensu lato has more recently received attention for its role in woody plant decomposition and in mycorrhizal symbiosis with specific plants. Previous phylogenetic analyses suggest that around 50 species are recognized globally. Despite this previous work, no studies have analyzed the global species richness and distribution of the genus using data derived from fungal community sequencing datasets or barcoding initiatives. To assess the global diversity and species richness of Armillaria, we mined publicly available sequencing datasets derived from numerous primer regions for the ribosomal operon, as well as ITS sequences deposited on Genbank, and clustered them akin to metabarcoding studies. Our estimates reveal that species richness ranges from 50 to 60 species, depending on whether the ITS1 or ITS2 marker is used. Eastern Asia represents the biogeographic region with the highest species richness. We also assess the overlap of species across geographic regions and propose some hypotheses regarding the drivers of variability in species diversity and richness between different biogeographic regions.

Evolution of the Mode of Nutrition in Symbiotic and Saprotrophic Fungi in Forest Ecosystems

In this review, we highlight the main insights that have been gathered from recent developments using large-scale genomics of fungal saprotrophs and symbiotrophs (including ectomycorrhizal and orchid and ericoid mycorrhizal fungi) inhabiting forest ecosystems. After assessing the goals and motivations underlying our approach, we explore our current understanding of the limits and future potential of using genomics to understand the ecological roles of these forest fungi. Comparative genomics unraveled the molecular machineries involved in lignocellulose decomposition in wood decayers, soil and litter saprotrophs, and mycorrhizal symbionts. They also showed that transitions from saprotrophy to mutualism entailed widespread losses of lignocellulose-degrading enzymes; diversification of novel, lineage-specific symbiosis-induced genes; and convergent evolution of genetic innovations that facilitate the accommodationof mutualistic symbionts within their plant hosts. We also identify the major questions that remain unanswered and propose new avenues of genome-based research to understand the role of soil fungi in sustainable forest ecosystems.

The melanized layer of Armillaria ostoyae rhizomorphs: Its protective role and functions

Armillaria ostoyae (Romagn.) Herink is a highly pathogenic fungus that uses exploratory, cordlike structures called rhizomorphs to seek out new sources of nutrition, posing a parasitic threat to natural stands of trees, orchards, and vineyards. Rhizomorphs are notoriously difficult to destroy, and this resilience is due in large part to a melanized layer that protects the rhizomorph. While this structure has been previously observed, its structural and chemical defenses are yet to be discerned. Research was conducted on both lab-cultured and wild-harvested rhizomorph samples. While both environments produce rhizomorphs, only the wild-harvested rhizomorphs produced the melanized layer, allowing for direct investigation of its structure and properties. Imaging, chemical analysis, mechanical testing, and finite element modeling were used to understand the defense mechanisms provided by the melanized layer. Imaging showed a porous outer layer in both types of rhizomorphs, though the pores were smaller in the harvested melanized layer. This melanized layer contained calcium, which provides chemical defense against both human and natural control methods, but was absent from cultured samples. Nanoindentation resulted in a larger variance of hardness values for cultured rhizomorphs than for wild-harvested. Finite element analysis proved that the smaller pore structure of the melanized porous layer had the best balance between maximum deformation and resulting permanent deformation. These results allow for a better understanding of the defenses of this pathogenic fungus, which may lead to better control methods.

Transcriptomics Reveals the Putative Mycoparasitic Strategy of the Mushroom Entoloma abortivum on Species of the Mushroom Genus Armillaria

During mycoparasitism, a fungus—the host—is parasitized by another fungus—the mycoparasite. The genetic underpinnings of these relationships have been best characterized in ascomycete fungi. However, within basidiomycete fungi, there are rare instances of mushroom-forming species parasitizing the reproductive structures, or sporocarps, of other mushroom-forming species, which have been rarely investigated on a genetic level. One of the most enigmatic of these occurs between Entoloma abortivum and species of Armillaria, where hyphae of E. abortivum are hypothesized to disrupt the development of Armillaria sporocarps, resulting in the formation of carpophoroids. However, it remains unknown whether carpophoroids are the direct result of a mycoparasitic relationship. To address the nature of this unique interaction, we analyzed gene expression of field-collected Armillaria and E. abortivum sporocarps and carpophoroids. Transcripts in the carpophoroids are primarily from E. abortivum, supporting the hypothesis that this species is parasitizing Armillaria. Most notably, we identified differentially upregulated E. abortivum β-trefoil-type lectins in the carpophoroid, which we hypothesize bind to Armillaria cell wall galactomannoproteins, thereby mediating recognition between the mycoparasite and the host. The most differentially upregulated E. abortivum transcripts in the carpophoroid code for oxalate decarboxylases—enzymes that degrade oxalic acid. Oxalic acid is a virulence factor in many plant pathogens, including Armillaria species; however, E. abortivum has evolved a sophisticated strategy to overcome this defense mechanism. The number of gene models and genes that code for carbohydrate-active enzymes in the E. abortivum transcriptome was reduced compared to other closely related species, perhaps as a result of the specialized nature of this interaction.

IMPORTANCE By studying fungi that parasitize other fungi, we can understand the basic biology of these unique interactions. Studies focused on the genetic mechanisms regulating mycoparasitism between host and parasite have thus far concentrated on a single fungal lineage within the Ascomycota. The work presented here expands our understanding of mycoparasitic relationships to the Basidiomycota and represents the first transcriptomic study to our knowledge that examines fungal-fungal relationships in their natural setting. The results presented here suggest that even distantly related mycoparasites utilize similar mechanisms to parasitize their host. Given that species of the mushroom-forming pathogen Armillaria cause plant root-rot diseases in many agroecosystems, an enhanced understanding of this interaction may contribute to better control of these diseases through biocontrol applications.

Genomic and Experimental Investigations of Auriscalpium and Strobilurus Fungi Reveal New Insights into Pinecone Decomposition

Saprophytic fungi (SPF) play vital roles in ecosystem dynamics and decomposition. However, because of the complexity of living systems, our understanding of how SPF interact with each other to decompose organic matter is very limited. Here we studied their roles and interactions in the decomposition of highly specialized substrates between the two genera Auriscalpium and Strobilurus fungi-colonized fallen pinecones of the same plant sequentially. We obtained the genome sequences from seven fungal species with three pairs: A. orientale-S. luchuensis, A. vulgare-S. stephanocystis and A. microsporum-S. pachcystidiatus/S. orientalis on cones of Pinus yunnanensis, P. sylvestris and P. armandii, respectively, and the organic profiles of substrate during decomposition. Our analyses revealed evidence for both competition and cooperation between the two groups of fungi during decomposition, enabling efficient utilization of substrates with complementary profiles of carbohydrate active enzymes (CAZymes). The Auriscalpium fungi are highly effective at utilizing the primary organic carbon, such as lignin, and hemicellulose in freshly fallen cones, facilitated the invasion and colonization by Strobilurus fungi. The Strobilurus fungi have genes coding for abundant CAZymes to utilize the remaining organic compounds and for producing an arsenal of secondary metabolites such as strobilurins that can inhibit other fungi from colonizing the pinecones.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates—namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose–methanol–choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases—we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Symbiotic nitrogen fixation in the reproductive structures of a basidiomycete fungus

Nitrogen (N) fixation is a driving force for the formation of symbiotic associations between N₂-fixing bacteria and eukaryotes.¹ Limited examples of these associations are known in fungi, and none with sexual structures of non-lichenized species.^2-6 The basidiomycete Guyanagaster necrorhizus is a sequestrate fungus endemic to the Guiana Shield.⁷ Like the root rot-causing species in its sister genera Armillaria and Desarmillaria, G. necrorhizus sporocarps fruit from roots of decaying trees (Figures 1A-1C),⁸ and genome sequencing is consistent with observations that G. necrorhizus is a white-rotting decomposer. This species also represents the first documentation of an arthropod-dispersed sequestrate fungus. Numerous species of distantly related wood-feeding termites, which scavenge for N-rich food, feed on the mature spore-bearing tissue, or gleba, of G. necrorhizus. During feeding, mature spores adhere to termites for subsequent dispersal.⁹ Using chemical assays, isotope analysis, and high-throughput sequencing, we show that the sporocarps harbor actively N₂-fixing Enterobacteriaceae species and that the N content within fungal tissue increases with maturation. Untargeted proteomic profiling suggests that ATP generation in the gleba is accomplished via fermentation. The use of fermentation-an anaerobic process-indicates that the sporocarp environment is anoxic, likely an adaptation to protect the oxygen-sensitive nitrogenase enzyme. Sporocarps also have a thick outer covering, possibly to limit oxygen diffusion. The enriched N content within mature sporocarps may offer a dietary inducement for termites in exchange for spore dispersal. These results show that the flexible metabolic capacity of fungi may facilitate N₂-fixing associations, as well as higher-level organismal associations.

Draft Genome Sequence of the Ectomycorrhizal Fungus Astraeus odoratus from Northern Thailand

We report the draft genome sequence of Astraeus odoratus, an edible ectomycorrhizal fungus from northern Thailand. The assembled genome has a size of 45.1 Mb and 13,403 annotated protein-coding genes. This reference genome will provide a better understanding of the biology of mushroom-forming ectomycorrhizal fungi in the family Diplocystidiaceae.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.