Genome Analysis of Hypomyces perniciosus, the Causal Agent of Wet Bubble Disease of Button Mushroom (Agaricus bisporus)

Genome Sequencing of Three Pathogenic Fungi Provides Insights into the Evolution and Pathogenic Mechanisms of the Cobweb Disease on Cultivated Mushrooms

Fungal diseases not only reduce the yield of edible mushrooms but also pose potential threats to the preservation and quality of harvested mushrooms. Cobweb disease, caused primarily by fungal pathogens from the Hypocreaceae family, is one of the most significant diseases affecting edible mushrooms. Deciphering the genomes of these pathogens will help unravel the molecular basis of their evolution and identify genes responsible for pathogenicity. Here, we present high-quality genome sequences of three cobweb disease fungi: Hypomyces aurantius Cb-Fv, Cladobotryum mycophilum CB-Ab, and Cladobotryum protrusum CB-Mi, isolated from Flammulina velutipes, Agaricus bisporus, and Morchella importuna, respectively. The assembled genomes of H. aurantius, C. mycophilum, and C. protrusum are 33.19 Mb, 39.83 Mb, and 38.10 Mb, respectively. This is the first report of the genome of H. aurantius. Phylogenetic analysis revealed that cobweb disease pathogens are closely related and diverged approximately 17.51 million years ago. CAZymes (mainly chitinases, glucan endo-1,3-beta-glucosidases, and secondary metabolite synthases), proteases, KP3 killer proteins, lipases, and hydrophobins were found to be conserved and strongly associated with pathogenicity, virulence, and adaptation in the three cobweb pathogens. This study provides insights into the genome structure, genome organization, and pathogenicity of these three cobweb disease fungi, which will be a valuable resource for comparative genomics studies of cobweb pathogens and will help control this disease, thereby enhancing mushroom quality.

Genomic insights into the evolution and adaptation of secondary metabolite gene clusters in fungicolous species Cladobotryum mycophilum ATHUM6906

Mycophilic or fungicolous fungi can be found wherever fungi exist since they are able to colonize other fungi, which occupy a diverse range of habitats. Some fungicolous species cause important diseases on Basidiomycetes, and thus, they are the main reason for the destruction of mushroom cultivations. Nonetheless, despite their ecological significance, their genomic data remain limited. Cladobotryum mycophilum is one of the most aggressive species of the genus, destroying the economically important Agaricus bisporus cultivations. The 40.7 Mb whole genome of the Greek isolate ATHUM6906 is assembled in 16 fragments, including the mitochondrial genome and 2 small circular mitochondrial plasmids, in this study. This genome includes a comprehensive set of 12,282 protein coding, 56 rRNA, and 273 tRNA genes. Transposable elements, CAZymes, and pathogenicity related genes were also examined. The genome of C. mycophilum contained a diverse arsenal of genes involved in secondary metabolism, forming 106 biosynthetic gene clusters, which renders this genome as one of the most BGC abundant among fungicolous species. Comparative analyses were performed for genomes of species of the family Hypocreaceae. Some BGCs identified in C. mycophilum genome exhibited similarities to clusters found in the family Hypocreaceae, suggesting vertical heritage. In contrast, certain BGCs showed a scattered distribution among Hypocreaceae species or were solely found in Cladobotryum genomes. This work provides evidence of extensive BGC losses, horizontal gene transfer events, and formation of novel BGCs during evolution, potentially driven by neutral or even positive selection pressures. These events may increase Cladobotryum fitness under various environmental conditions and potentially during host-fungus interaction.

Characterization and Genome Analysis of Cladobotryum mycophilum, the Causal Agent of Cobweb Disease of Morchella sextelata in China

Cobweb disease is a fungal disease that can cause serious damage to edible mushrooms worldwide. To investigate cobweb disease in Morchella sextelata in Guizhou Province, China, we isolated and purified the pathogen responsible for the disease. Through morphological and molecular identification and pathogenicity testing on infected M. sextelata, we identified Cladobotryum mycophilum as the cause of cobweb disease in this region. This is the first known occurrence of this pathogen causing cobweb disease in M. sextelata anywhere in the world. We then obtained the genome of C. mycophilum BJWN07 using the HiFi sequencing platform, resulting in a high-quality genome assembly with a size of 38.56 Mb, 10 contigs, and a GC content of 47.84%. We annotated 8428 protein-coding genes in the genome, including many secreted proteins, host interaction-related genes, and carbohydrate-active enzymes (CAZymes) related to the pathogenesis of the disease. Our findings shed new light on the pathogenesis of C. mycophilum and provide a theoretical basis for developing potential prevention and control strategies for cobweb disease.

Edible mushroom industry in China: current state and perspectives

China is one of the largest producers and exporters of wild edible fungi in the world. Cultivation mushroom production value ranks within the top five after grain, vegetable, fruit, and edible oil plantation, greater than sugar, cotton, and tobacco business. More than 40 new varieties of high market value mushrooms from our group were highlighted in this article. Mushroom cultivations have a high impact on China's poverty alleviation program, with earnings at least ten times higher than rice and corn. The products were exported to 137 countries and regions, mainly to Japan, South Korea, ASEAN, the USA, and the European Union, among these, Japan is the biggest import market for cultivated mushrooms from China. Rapid development in the market and an increased demand for edible fungi generally enhance the economy of domestic edible fungi. We are the leading research group in logical farm design that is HACCP-certified to reduce the cost of investment for agriculture, thus broadening the consumption market of edible mushrooms and forming a demand-oriented leading industry for the promotion of human health. The enterprise needs to re-examine the operation plan and the strategic thinking to improve the fundamental drivers based on the available resources of the locality. Mushrooms growing intergrate with upgraded technologies and equipment to become smart agriculture which have smart production and intelligent factories. The purpose of suitable products will not change: delicious, nutritious, healthy, and modern. KEY POINTS: • Cultivation mushroom production value ranks within the top five after grain, vegetable, fruit, and edible oil plantation, greater than sugar, cotton, and tobacco business. • Mushroom cultivations have a high impact on China's poverty alleviation program, with earnings at least ten times higher than rice and corn. • The development of transportation and industrialization of mushrooms facilitate the modernization of mushroom industry in China.

Infection Process and Genome Assembly Provide Insights into the Pathogenic Mechanism of Destructive Mycoparasite Calcarisporium cordycipiticola with Host Specificity

Calcarisporium cordycipiticola is the pathogen in the white mildew disease of Cordyceps militaris, one of the popular mushrooms. This disease frequently occurs and there is no effective method for disease prevention and control. In the present study, C. militaris is found to be the only host of C. cordycipiticola, indicating strict host specificity. The infection process was monitored by fluorescent labeling and scanning and transmission electron microscopes. C. cordycipiticola can invade into the gaps among hyphae of the fruiting bodies of the host and fill them gradually. It can degrade the hyphae of the host by both direct contact and noncontact. The parasitism is initially biotrophic, and then necrotrophic as mycoparasitic interaction progresses. The approximate chromosome-level genome assembly of C. cordycipiticola yielded an N50 length of 5.45 Mbp and a total size of 34.51 Mbp, encoding 10,443 proteins. Phylogenomic analysis revealed that C. cordycipiticola is phylogenetically close to its specific host, C. militaris. A comparative genomic analysis showed that the number of CAZymes of C. cordycipiticola was much less than in other mycoparasites, which might be attributed to its host specificity. Secondary metabolite cluster analysis disclosed the great biosynthetic capabilities and potential mycotoxin production capability. This study provides insights into the potential pathogenesis and interaction between mycoparasite and its host.

Complete Genomic Characterization and Identification of Saccharomycopsisphalluae sp. nov., a Novel Pathogen Causes Yellow Rot Disease on Phallus rubrovolvatus

“Hongtuozhusun” (Phallus rubrovolvatus) is an important edible and medicinal mushroom endemic to Southwest China. However, yellow rot disease is a severe disease of P. rubrovolvatus that occurs extensively in Guizhou Province. It has caused major economic losses and hinders the development of the P. rubrovolvatus industry. In this study, 28 microorganism strains were isolated from diseased fruiting bodies of P. rubrovolvatus at various stages, two of which were confirmed to be pathogenic based on Koch’s postulates. These two strains are introduced herein as Saccharomycopsisphalluae sp. nov. based on morphological, physiological, and molecular analysis. We reported a high-quality de novo sequencing and assembly of the S. phalluae genome using single-molecule real-time sequencing technology. The whole genome was approximately 14.148 Mb with a G+C content of 43.55%. Genome assembly generated 8 contigs with an N50 length of 1,822,654 bp. The genome comprised 5966 annotated protein-coding genes. This is the first report of mushroom disease caused by Saccharomycopsis species. We expect that the information on genome properties, particularly in pathogenicity-related genes, assist in developing effective control measures in order to prevent severe losses and make amendments in management strategies.

De novo genome sequencing of mycoparasite Mycogone perniciosa strain MgR1 sheds new light on its biological complexity

Mycogone perniciosa is a mycoparasite causing Wet Bubble Diseases (WBD) of Agaricus bisporus. In the present study, the whole genome of M. perniciosa strain MgR1 was sequenced using Illumina NextSeq500 platform. This sequencing generated 8.03 Gb of high-quality data and a draft genome of 39 Mb was obtained through a de novo assembly of the high-quality reads. The draft genome resulted into prediction of 9276 genes from the 1597 scaffolds. NCBI-based homology analysis revealed the identification of 8660 genes. Notably, non-redundant protein database analysis of the M. perniciosa strain MgR1 revealed its close relation with the Trichoderma arundinaceum. Moreover, ITS-based phylogenetic analysis showed the highest similarity of M. perniciosa strain MgR1 with Hypomyces perniciosus strain CBS 322.22 and Mycogone perniciosa strain PPRI 5784. Annotation of the 3917 genes of M. perniciosa strain MgR1 grouped in three major categories viz. biological process (2583 genes), cellular component (2013 genes), and molecular function (2919 genes). UniGene analysis identified 2967 unique genes in M. perniciosa strain MgR1. In addition, prediction of the secretory and pathogenicity-related genes based on the fungal database indicates that 1512 genes (16% of predicted genes) encode for secretory proteins. Moreover, out of 9276 genes, 1296 genes were identified as pathogenesis-related proteins matching with 51 fungal and bacterial genera. Overall, the key pathogenic genes such as lysine M protein domain genes, G protein, hydrophobins, and cytochrome P450 were also observed. The draft genome of MgR1 provides an understanding of pathogenesis of WBD in A. bisporus and could be utilized to develop novel management strategies.

Intraspecific Mitochondrial DNA Comparison of Mycopathogen Mycogone perniciosa Provides Insight Into Mitochondrial Transfer RNA Introns

Mycogone perniciosa is the main causative agent of wet bubble disease, which causes severe damage to the production of the cultivated mushroom Agaricus bisporus around the world. Whole-genome sequencing of 12 isolates of M. perniciosa was performed using the Illumina sequencing platform, and the obtained paired-end reads were used to assemble complete mitochondrial genomes. Intraspecific comparisons of conserved protein-coding genes, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes, introns, and intergenic regions were conducted. Five different mitochondrial DNA (mtDNA) haplotypes were detected among the tested isolates, ranging from 89,080 to 93,199 bp in length. All of the mtDNAs contained the same set of 14 protein-coding genes and 2 rRNA and 27 tRNA genes, which shared high sequence similarity. In contrast, the number, insertion sites, and sequences of introns varied greatly among the mtDNAs. Eighteen of 43 intergenic regions differed among the isolates, reflecting 65 single nucleotide polymorphisms, 76 indels, and the gain/loss of nine long fragments. Intraspecific comparison revealed that two introns were located within tRNA genes, which is the first detailed description of mitochondrial tRNA introns. Intronic sequence comparison within the same insertion sites revealed the formation process of two introns, which also illustrated a fast evolutionary rate of introns among M. perniciosa isolates. Based on the intron distribution pattern, a pair of universal primers and four pairs of isolate-specific primers were designed and were used to identify the five mtDNA types. In summary, the rapid gain or loss of mitochondrial introns could be an ideal marker for population genetics analysis.

Control of Fungal Diseases in Mushroom Crops while Dealing with Fungicide Resistance: A Review

Mycoparasites cause heavy losses in commercial mushroom farms worldwide. The negative impact of fungal diseases such as dry bubble (Lecanicillium fungicola), cobweb (Cladobotryum spp.), wet bubble (Mycogone perniciosa), and green mold (Trichoderma spp.) constrains yield and harvest quality while reducing the cropping surface or damaging basidiomes. Currently, in order to fight fungal diseases, preventive measurements consist of applying intensive cleaning during cropping and by the end of the crop cycle, together with the application of selective active substances with proved fungicidal action. Notwithstanding the foregoing, the redundant application of the same fungicides has been conducted to the occurrence of resistant strains, hence, reviewing reported evidence of resistance occurrence and introducing unconventional treatments is worthy to pave the way towards the design of integrated disease management (IDM) programs. This work reviews aspects concerning chemical control, reduced sensitivity to fungicides, and additional control methods, including genomic resources for data mining, to cope with mycoparasites in the mushroom industry.

Genome-Wide Identification and Analysis of Chitinase GH18 Gene Family in Mycogone perniciosa

Mycogone perniciosa causes wet bubble disease in Agaricus bisporus and various Agaricomycetes species. In a previous work, we identified 41 GH18 chitinase genes and other pathogenicity-related genes in the genome of M. perniciosa Hp10. Chitinases are enzymes that degrade chitin, and they have diverse functions in nutrition, morphogenesis, and pathogenesis. However, these important genes in M. perniciosa have not been fully characterized, and their functions remain unclear. Here, we performed a genome-wide analysis of M. perniciosa GH18 genes and analyzed the transcriptome profiles and GH18 expression patterns in M. perniciosa during the time course of infection in A. bisporus. Phylogenetic analysis of the 41 GH18 genes with those of 15 other species showed that the genes were clustered into three groups and eight subgroups based on their conserved domains. The GH18 genes clustered in the same group shared different gene structures but had the same protein motifs. All GH18 genes were localized in different organelles, were unevenly distributed on 11 contigs, and had orthologs in the other 13 species. Twelve duplication events were identified, and these had undergone both positive and purifying selection. The transcriptome analyses revealed that numerous genes, including transporters, cell wall degrading enzymes (CWDEs), cytochrome P450, pathogenicity-related genes, secondary metabolites, and transcription factors, were significantly upregulated at different stages of M. perniciosa Hp10 infection of A. bisporus. Twenty-three out of the 41 GH18 genes were differentially expressed. The expression patterns of the 23 GH18 genes were different and were significantly expressed from 3 days post-inoculation of M. perniciosa Hp10 in A. bisporus. Five differentially expressed GH18 genes were selected for RT-PCR and gene cloning to verify RNA-seq data accuracy. The results showed that those genes were successively expressed in different infection stages, consistent with the previous sequencing results. Our study provides a comprehensive analysis of pathogenicity-related and GH18 chitinase genes’ influence on M. perniciosa mycoparasitism of A. bisporus. Our findings may serve as a basis for further studies of M. perniciosa mycoparasitism, and the results have potential value for improving resistance in A. bisporus and developing efficient disease-management strategies to mitigate wet bubble disease.

Genome Sequencing of Paecilomyces Penicillatus Provides Insights into Its Phylogenetic Placement and Mycoparasitism Mechanisms on Morel Mushrooms

Morels (Morchella spp.) are popular edible fungi with significant economic and scientific value. However, white mold disease, caused by Paecilomyces penicillatus, can reduce morel yield by up to 80% in the main cultivation area in China. Paecilomyces is a polyphyletic genus and the exact phylogenetic placement of P. penicillatus is currently still unclear. Here, we obtained the first high-quality genome sequence of P. penicillatus generated through the single-molecule real-time (SMRT) sequencing platform. The assembled draft genome of P. penicillatus was 40.2 Mb, had an N50 value of 2.6 Mb and encoded 9454 genes. Phylogenetic analysis of single-copy orthologous genes revealed that P. penicillatus is in Hypocreales and closely related to Hypocreaceae, which includes several genera exhibiting a mycoparasitic lifestyle. CAZymes analysis demonstrated that P. penicillatus encodes a large number of fungal cell wall degradation enzymes. We identified many gene clusters involved in the production of secondary metabolites known to exhibit antifungal, antibacterial, or insecticidal activities. We further demonstrated through dual culture assays that P. penicillatus secretes certain soluble compounds that are inhibitory to the mycelial growth of Morchella sextelata. This study provides insights into the correct phylogenetic placement of P. penicillatus and the molecular mechanisms that underlie P. penicillatus pathogenesis.

Sphaerostilbellins, New Antimicrobial Aminolipopeptide Peptaibiotics from Sphaerostilbella toxica

Sphaerostilbella toxica is a mycoparasitic fungus that can be found parasitizing wood-decay basidiomycetes in the southern USA. Organic solvent extracts of fermented strains of S. toxica exhibited potent antimicrobial activity, including potent growth inhibition of human pathogenic yeasts Candida albicans and Cryptococcus neoformans, the respiratory pathogenic fungus Aspergillus fumigatus, and the Gram-positive bacterium Staphylococcus aureus. Bioassay-guided separations led to the purification and structure elucidation of new peptaibiotics designated as sphaerostilbellins A and B. Their structures were established mainly by analysis of NMR and HRMS data, verification of amino acid composition by Marfey's method, and by comparison with published data of known compounds. They incorporate intriguing structural features, including an N-terminal 2-methyl-3-oxo-tetradecanoyl (MOTDA) residue and a C-terminal putrescine residue. The minimal inhibitory concentrations for sphaerostilbellins A and B were measured as 2 μM each for C. neoformans, 1 μM each for A. fumigatus, and 4 and 2 μM, respectively, for C. albicans. Murine macrophage cells were unaffected at these concentrations.

Genomic Features of Cladobotryum dendroides, Which Causes Cobweb Disease in Edible Mushrooms, and Identification of Genes Related to Pathogenicity and Mycoparasitism

Cladobotryum dendroides, which causes cobweb disease in edible mushrooms, is one of the major fungal pathogens. Our previous studies focused on the genetic and morphological characterization of this fungus, as well as its pathogenicity and the identification of appropriate fungicides. However, little is known about the genome characters, pathogenic genes, and molecular pathogenic mechanisms of C. dendroides. Herein, we reported a high-quality de novo genomic sequence of C. dendroides and compared it with closely-related fungi. The assembled C. dendroides genome was 36.69 Mb, consisting of eight contigs, with an N50 of 4.76 Mb. This genome was similar in size to that of C. protrusum, and shared highly conserved syntenic blocks and a few inversions with C. protrusum. Phylogenetic analysis revealed that, within the Hypocreaceae, Cladobotryum was closer to Mycogone than to Trichoderma, which is consistent with phenotypic evidence. A significant number of the predicted expanded gene families were strongly associated with pathogenicity, virulence, and adaptation. Our findings will be instrumental for the understanding of fungi-fungi interactions, and for exploring efficient management strategies to control cobweb disease.

Genetic and Pathogenic Variability of Mycogone perniciosa Isolates Causing Wet Bubble Disease on Agaricus bisporus in China

Wet bubble disease, caused by Mycogone perniciosa, is a major threat to Agaricus bisporus production in China. In order to understand the variability in genetic, pathogenicity, morphology, and symptom production of the fungus, 18 isolates of the pathogen were collected from diseased A. bisporus in different provinces in China. The isolates were characterized by a combination of morphological, cultural, molecular, and pathogenicity testing on different strains of A. bisporus and amplified fragment length polymorphism (AFLP) analysis. The 18 isolates were identified by Koch’s postulate and confirmed different pathogenic variability among them. The yellow to brown isolates were more virulent than the white isolates. AFLP markers clustered the isolates into two distinct groups based on their colony color, with a high level of polymorphism of Jaccard similarities range from 0.39% to 0.64%. However, there was no evidence of an association between the genetic diversity and the geographical origin of the isolates. Through knowledge of the genetic diversity, phenotypic virulence of M. perniciosa is a key factor for successful breeding of resistant strains of A. bisporus and developing of an integrated disease management strategy to manage wet bubble disease of A. bisporus.