Composition and expression of genes encoding carbohydrate-active enzymes in the straw-degrading mushroom Volvariella volvacea

Screening candidate genes related to volatile synthesis in shiitake mushrooms and construction of regulatory networks to effectively improve mushroom aroma

BACKGROUND

Metabolite formation is a physiological stress response during the growth and development of shiitake mushrooms (Lentinula edodes). The characteristic flavor metabolites are important quality components in shiitake mushrooms. To investigate the formation mechanisms of characteristic flavor metabolites, transcriptome analyses were performed on shiitake mushrooms harvested at different growth stages.

RESULTS

In total, 30 genes related to the synthesis of characteristic volatiles of mushrooms were identified via screening. Through KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of the selected genes and correlation analyses of gene expressions, the main volatile synthesis pathways were determined as histidine metabolism, glutathione metabolism and biosynthesis of unsaturated fatty acids. Gene cluster and correlation analyses were performed to clarify the combined effects of different genes in the enzymatic reactions. Further, a correlation network of candidate genes was built based on the gene expression levels.

CONCLUSION

The activities of flavor synthases and the content of characteristic flavor metabolites were analyzed; the enzyme activity changes and metabolic product distribution sites were clarified. A synthesis and regulation network was constructed for the candidate genes and characteristic volatiles, and information was obtained for 16 hub genes. Moreover, it was essential to identify and characterize the key genes and synthases involved in the synthesis of the characteristic volatiles of mushrooms. This information provides us with a better understanding of the biosynthesis and regulation of the volatiles, which will lay the foundation for improving the quality of shiitake mushrooms. © 2021 Society of Chemical Industry.

Analysis of the Genome Sequence of Strain GiC-126 of Gloeostereum incarnatum with Genetic Linkage Map

Gloeostereum incarnatum has edible and medicinal value and was first cultivated and domesticated in China. We sequenced the G. incarnatum monokaryotic strain GiC-126 on an Illumina HiSeq X Ten system and obtained a 34.52-Mb genome assembly sequence that encoded 16,895 predicted genes. We combined the GiC-126 genome with the published genome of G. incarnatum strain CCMJ2665 to construct a genetic linkage map (GiC-126 genome) that had 10 linkage groups (LGs), and the 15 assembly sequences of CCMJ2665 were integrated into 8 LGs. We identified 1912 simple sequence repeat (SSR) loci and detected 700 genes containing 768 SSRs in the genome; 65 and 100 of them were annotated with gene ontology (GO) terms and KEGG pathways, respectively. Carbohydrate-active enzymes (CAZymes) were identified in 20 fungal genomes and annotated; among them, 144 CAZymes were annotated in the GiC-126 genome. The A mating-type locus (MAT-A) of G. incarnatum was located on scaffold885 at 38.9 cM of LG1 and was flanked by two homeodomain (HD1) genes, mip and beta-fg. Fourteen segregation distortion markers were detected in the genetic linkage map, all of which were skewed toward the parent GiC-126. They formed three segregation distortion regions (SDR1-SDR3), and 22 predictive genes were found in scaffold1920 where three segregation distortion markers were located in SDR1. In this study, we corrected and updated the genomic information of G. incarnatum. Our results will provide a theoretical basis for fine gene mapping, functional gene cloning, and genetic breeding the follow-up of G. incarnatum.

A Comparative Transcriptome Analysis of Volvariella volvacea Identified the Candidate Genes Involved in Fast Growth at the Mycelial Growth Stage

The edible straw mushroom, Volvariella volvacea, is one of the most important cultivated mushrooms in tropical and sub-tropical regions. Strain improvement for V. volvacea is difficult because of the unknown mechanisms involved in its growth regulation and substrate utilization. A comparative physiological and transcriptomic study was conducted between two commercially available straw mushroom strains (v9 and v26) to explore their fast-growth regulation mechanism(s). The physiological study showed that V. volvacea v9 had a shorter growth cycle and higher biological efficiency (4% higher) than that in v26. At least 14,556 unigenes were obtained from the four cDNA libraries (two replicates per strain). Among them, the expression of 1597 unigenes was up-regulated while 1352 were down-regulated. Four heat-shock proteins were highly expressed in v9, showing that v9 has the better ability to handle stresses and/or environmental changes. Moreover, up to 14 putative transporter genes were expressed at a higher level in v9 than those in v26, implying that v9 has a better ability to transport nutrients or export xenobiotics efficiently. Our report allows to identify the candidate genes involved in the fast growth requirement of V. volvacea, which represents a valuable resource for strain improvement in this commercially important edible mushroom.

Genome Survey and Transcriptome Analysis on Mycelia and Primordia of Agaricus blazei

Agaricus blazei, a type of edible straw-rotting mushroom with somewhat sweet taste and fragrance of almonds, has attracted considerable scientific and practical attention. High-throughput Illumina PE150 and PacBio RSII platform were employed to generate a genomic sequence. De novo assembly generated 36 contigs with 38,686,133 bp in size, containing 10,119 putative predicted genes. Additionally, we also studied transcriptional regulation of the mycelia and the primordia for exploration of genes involved in fruiting body formation. Expression profiling analysis revealed that 2,164 genes were upregulated in mycelia and 1,557 in primordia. Functional enrichment showed that differentially expressed genes associated with response to stress, ribosome biogenesis, arginine biosynthesis, and steroid biosynthesis pathway were more active in fruiting body. The genome and transcriptome analysis of A. blazei provide valuable sequence resources and contribute to our understanding of genes related to the biosynthesis pathway of polysaccharide and benzaldehyde, as well as the fruiting body formation.

Gene Mutations in Ganoderma lucidum During Long-Term Preservation by Repeated Subculturing

Subculturing is frequently used for the preservation of basidiomycetes. In this study, to assess and verify the risks of repeated subculturing on the long-term preservation of strains of culture collections, we performed single nucleotide polymorphism (SNP) analysis in genes encoding enzymes of the mevalonate pathway, 1,3-β-glucan synthesis, lignin degradation, and the tricarboxylic acid (TCA) cycle of mycelia before and after preserving for a 4-year period by the subculturing 30 times every 45 days of Ganoderma lucidum NBRC 8346. As a result of analyzing 60 genes of the strain, SNPs were found in 18 genes, and 14 of them were found in the exon region. Nonsynonymous coding SNPs were found in two genes (atoB_2, hmgr) encoding enzymes of mevalonate pathway and five genes (lcc1_9, lcc1_11, lcc1_13, dslcc6, aa5_1_9) encoding enzymes of lignin degradation after subculturing of G. lucidum NBRC 8346.

Comparative genomics of 40 edible and medicinal mushrooms provide an insight into the evolution of lignocellulose decomposition mechanisms

Diversity comparison and phylogenetic analyses of carbohydrate-active enzymes (CAZymes), auxiliary activities (AAs) and cytochromes P450 among 40 fungi, which are based on different nutritional pathways, help clarify and explain their divergence and improvement of various life-styles. Molecular clock analyses allow us to understand the evolutionary and developmental rules in decomposition gene families. Our results suggested that fungi in different ecological types acquired an obvious preference on specific decomposing gene families during evolutionary selection. White rot and litter saprotrophic fungi possessed more complete types of varied degradation gene families and were superior in quantities. With evolution and development of lignocellulose decomposition mechanism, certain families (like CBM1, GH6, GH7, GH10, and CYP53) disappeared in brown rot fungi and symbiotic fungi. In addition, the earlier time of phylogenetic divergence determined the more integrated and larger decomposition families. And various gains and losses in gene quantity of varied decomposition families led in particularly phylogenetic clades or nodes, then accelerated in forming varied ecotypes of species.

Three complete mitochondrial genomes of straw-rotting edible fungus Volvariella volvacea using next generation sequencing

The straw-rotting edible fungus Volvariella volvacea is a widely cultivated edible fungus across China and Southeast Asian countries. Three complete mitochondrial genomes of V. volvacea from China, Thailand, and India were determined using the next-generation sequencing technology. The genome sizes of the three strains (China, Thailand, and India) were 62,541 bp, 64,531 bp, and 65,668 bp with GC contents of 38.46%, 38.56%, and 38.52%, respectively. All the genomes encoded 14 conserved protein-coding genes, the small ribosomal RNA subunits (rns), large ribosomal RNA subunits (rnl), and 23 tRNAs were located on the same strand. In the putative protein-coding genes, four introns were distributed in cox1 in the genomes of V23-1 and V8. 5 introns (four introns invaded into cox1and one intron invaded into cob) were detected in Tai8. The phylogenetic analysis confirmed that V. volvacea was a number of Agaricales. This mitochondrial genome may open new avenues for understanding the phylogeny and evolution of Pluteaceae and Agaricales.

Selection and Evaluation of Appropriate Reference Genes for RT-qPCR Normalization of Volvariella volvacea Gene Expression under Different Conditions

Volvariella volvacea (V. volvacea), commonly referred to as Chinese (paddy straw) mushroom, is a basidiomycete with a protein-rich volva and pileus. Selecting appropriate reference genes is a crucial step in the normalization of quantitative real-time PCR data. Therefore, 12 candidate reference genes were selected from the V. volvacea transcriptome based on previous studies and then BestKeeper, geNorm, and NormFinder were used to identify reference genes stably expressed during different developmental stages and conditions. Of the 12 candidate reference genes, SPRY domain protein (SPRYp), alpha-tubulin (TUBα), cyclophilin (CYP), L-asparaginase (L-asp), and MSF1-domain-containing protein (MSF1) were the most stably expressed under different experimental conditions, while 18S ribosomal RNA (18S), 28S ribosomal RNA (28S), and beta-actin (ACTB) were the least stably expressed. This investigation not only revealed potential factors influencing the suitability of reference genes, but also identified optimal reference genes from a pool of candidate genes under a wide range of conditions.

A genetic linkage map of Pleurotus tuoliensis integrated with physical mapping of the de novo sequenced genome and the mating type loci

Background

Pleurotus tuoliensis (Bailinggu) is a commercially cultivated mushroom species with an increasing popularity in China and other Asian countries. Commercial profits are now low, mainly due to a low yield, long cultivation period and sensitivity to diseases. Breeding efforts are thus required to improve agronomical important traits. Developing saturated genetic linkage and physical maps is a start for applying genetic and molecular approaches to accelerate the precise breeding programs.

Results

Here we present a genetic linkage map for P. tuoliensis constructed by using 115 haploid monokaryons derived from a hybrid strain H6. One thousand one hundred and eighty-two SNP markers developed by 2b–RAD (type IIB restriction-site associated DNA) approach were mapped to 12 linkage groups. The map covers 1073 cM with an average marker spacing of 1.0 cM. The genome of P. tuoliensis was de novo sequenced as 40.8 Mb and consisted of 500 scaffolds (>500 bp), which showed a high level of colinearity to the genome of P. eryngii var. eryngii. A total of 97.4% SNP markers (1151) were physically localized on 78 scaffolds, and the physical length of these anchored scaffolds were 33.9 Mb representing 83.1% of the whole genome. Mating type loci A and B were mapped on separate linkage groups and identified physically on the assembled genomes. Five putative pheromone receptors and two putative pheromone precursors were identified for the mating type B locus.

Conclusions

This study reported a first genetic linkage map integrated with physical mapping of the de novo sequenced genome and the mating type loci of an important cultivated mushroom in China, P. tuoliensis. The de novo sequenced and annotated genome, assembled using a 2b–RAD generated linkage map, provides a basis for marker-assisted breeding of this economic important mushroom species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4421-z) contains supplementary material, which is available to authorized users.

Small GTPases and Stress Responses of vvran1 in the Straw Mushroom Volvariella volvacea

Small GTPases play important roles in the growth, development and environmental responses of eukaryotes. Based on the genomic sequence of the straw mushroom Volvariella volvacea, 44 small GTPases were identified. A clustering analysis using human small GTPases as the references revealed that V. volvacea small GTPases can be grouped into five families: nine are in the Ras family, 10 are in the Rho family, 15 are in the Rab family, one is in the Ran family and nine are in the Arf family. The transcription of vvran1 was up-regulated upon hydrogen peroxide (H₂O₂) stress, and could be repressed by diphenyleneiodonium chloride (DPI), a NADPH oxidase-specific inhibitor. The number of vvran1 transcripts also increased upon cold stress. Diphenyleneiodonium chloride, but not the superoxide dismutase (SOD) inhibitor diethy dithiocarbamate (DDC), could suppress the up-regulation of vvran1 gene expression to cold stress. These results combined with the high correlations between gene expression and superoxide anion (O₂⁻) generation indicated that vvran1 could be one of the candidate genes in the downstream of O₂⁻ mediated pathways that are generated by NADPH oxidase under low temperature and oxidative stresses.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.

Fruiting Body Formation in Volvariella volvacea Can Occur Independently of Its MAT-A-Controlled Bipolar Mating System, Enabling Homothallic and Heterothallic Life Cycles

Volvariella volvacea is an important crop in Southeast Asia, but erratic fruiting presents a serious challenge for its production and breeding. Efforts to explain inconsistent fruiting have been complicated by the multinucleate nature, typical lack of clamp connections, and an incompletely identified sexual reproductive system. In this study, we addressed the life cycle of V. volvacea using whole genome sequencing, cloning of MAT loci, karyotyping of spores, and fruiting assays. Microscopy analysis of spores had previously indicated the possible coexistence of heterothallic and homothallic life cycles. Our analysis of the MAT loci showed that only MAT-A, and not MAT-B, controlled heterokaryotization. Thus, the heterothallic life cycle was bipolar. Karyotyping of single spore isolates (SSIs) using molecular markers supported the existence of heterokaryotic spores. However, most SSIs were clearly not heterokaryotic, yet contained structural variation (SV) markers relating to both alleles of both parents. Heterokaryons from crossed, self-sterile homokaryons could produce fruiting bodies, agreeing with bipolar heterothallism. Meanwhile, some SSIs with two different MAT-A loci also produced fruiting bodies, which supported secondary homothallism. Next, SSIs that clearly contained only one MAT-A locus (homothallism) were also able to fruit, demonstrating that self-fertile SSIs were not, per definition, secondary homothallic, and that a third life cycle or genetic mechanism must exist. Finally, recombination between SV markers was normal, yet 10 out of 24 SV markers showed 1:2 or 1:3 distributions in the spores, and large numbers of SSIs contained doubled SV markers. This indicated selfish genes, and possibly partial aneuploidy.

The Accordant Trend of Both Parameters (rgs Expression and cAMP Content) Follows the Pattern of Development of Fruiting Body in Volvariella volvacea

The formation of fruiting body in Volvariella volvacea is affected by endogenous genes and environmental factors. However, its regulation at a molecular level is still poorly understood. To study the genes involved in the formation of fruiting body, we cloned a new regulator of the G protein signaling (RGS) encoding gene (rgs) from V. volvacea. Phylogenetic analysis showed that RGS in V. volvacea and other basidiomycete RGS proteins from Schizophyllum commune and Coprinus cinereus belong to the same clade. In addition, we assayed intracellular cAMP content in the three developmental stages (mycelium, fruiting body primordia, and button). We also found that the expression of rgs was highly positively correlated to the content of intracellular cAMP during fruiting body formation. The conserved protein sequences and expression of rgs, together with high concent of cAMP at primordia tissue, suggested that rgs gene and cAMP may play a crucial role in fruiting body formation in V. volvacea.

Structural Variation (SV) Markers in the Basidiomycete Volvariella volvacea and Their Application in the Construction of a Genetic Map

Molecular markers and genetic maps are useful tools in genetic studies. Novel molecular markers and their applications have been developed in recent years. With the recent advancements in sequencing technology, the genomic sequences of an increasingly great number of fungi have become available. A novel type of molecular marker was developed to construct the first reported linkage map of the edible and economically important basidiomycete Volvariella volvacea by using 104 structural variation (SV) markers that are based on the genomic sequences. Because of the special and simple life cycle in basidiomycete, SV markers can be effectively developed by genomic comparison and tested in single spore isolates (SSIs). This stable, convenient and rapidly developed marker may assist in the construction of genetic maps and facilitate genomic research for other species of fungi.

Plant-polysaccharide-degrading enzymes from Basidiomycetes

SUMMARY

Basidiomycete fungi subsist on various types of plant material in diverse environments, from living and dead trees and forest litter to crops and grasses and to decaying plant matter in soils. Due to the variation in their natural carbon sources, basidiomycetes have highly varied plant-polysaccharide-degrading capabilities. This topic is not as well studied for basidiomycetes as for ascomycete fungi, which are the main sources of knowledge on fungal plant polysaccharide degradation. Research on plant-biomass-decaying fungi has focused on isolating enzymes for current and future applications, such as for the production of fuels, the food industry, and waste treatment. More recently, genomic studies of basidiomycete fungi have provided a profound view of the plant-biomass-degrading potential of wood-rotting, litter-decomposing, plant-pathogenic, and ectomycorrhizal (ECM) basidiomycetes. This review summarizes the current knowledge on plant polysaccharide depolymerization by basidiomycete species from diverse habitats. In addition, these data are compared to those for the most broadly studied ascomycete genus, Aspergillus, to provide insight into specific features of basidiomycetes with respect to plant polysaccharide degradation.

Crosses between monokaryons of Pleurotus sapidus or Pleurotus florida show an improved biotransformation of (+)-valencene to (+)-nootkatone

Several hundred monokaryotic and new dikaryotic strains derived thereof were established from (+)-valencene tolerant Pleurotus species. When grouped according to their growth rate on agar plates and compared to the parental of Pleurotus sapidus 69, the slowly growing monokaryons converted (+)-valencene more efficiently to the grapefruit flavour compound (+)-nootkatone. The fast growing monokaryons and the slow×slow and the fast×fast dikaryotic crosses showed similar or inferior yields. Some slow×fast dikaryons, however, exceeded the biotransformation capability of the parental dikaryon significantly. The activity of the responsible enzyme, lipoxygenase, showed a weak correlation with the yields of (+)-nootkatone indicating that the determination of enzyme activity using the primary substrate linoleic acid may be misleading in predicting the biotransformation efficiency. This exploratory study indicated that a classical genetics approach resulted in altered and partly improved terpene transformation capability (plus 60%) and lipoxygenase activity of the strains.

Whole genome and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose degradation

Flammulina velutipes is a fungus with health and medicinal benefits that has been used for consumption and cultivation in East Asia. F. velutipes is also known to degrade lignocellulose and produce ethanol. The overlapping interests of mushroom production and wood bioconversion make F. velutipes an attractive new model for fungal wood related studies. Here, we present the complete sequence of the F. velutipes genome. This is the first sequenced genome for a commercially produced edible mushroom that also degrades wood. The 35.6-Mb genome contained 12,218 predicted protein-encoding genes and 287 tRNA genes assembled into 11 scaffolds corresponding with the 11 chromosomes of strain KACC42780. The 88.4-kb mitochondrial genome contained 35 genes. Well-developed wood degrading machinery with strong potential for lignin degradation (69 auxiliary activities, formerly FOLymes) and carbohydrate degradation (392 CAZymes), along with 58 alcohol dehydrogenase genes were highly expressed in the mycelium, demonstrating the potential application of this organism to bioethanol production. Thus, the newly uncovered wood degrading capacity and sequential nature of this process in F. velutipes, offer interesting possibilities for more detailed studies on either lignin or (hemi-) cellulose degradation in complex wood substrates. The mutual interest in wood degradation by the mushroom industry and (ligno-)cellulose biomass related industries further increase the significance of F. velutipes as a new model.

First structural insights into α-L-arabinofuranosidases from the two GH62 glycoside hydrolase subfamilies

α-L-arabinofuranosidases are glycoside hydrolases that specifically hydrolyze non-reducing residues from arabinose-containing polysaccharides. In the case of arabinoxylans, which are the main components of hemicellulose, they are part of microbial xylanolytic systems and are necessary for complete breakdown of arabinoxylans. Glycoside hydrolase family 62 (GH62) is currently a small family of α-L-arabinofuranosidases that contains only bacterial and fungal members. Little is known about the GH62 mechanism of action, because only a few members have been biochemically characterized and no three-dimensional structure is available. Here, we present the first crystal structures of two fungal GH62 α-L-arabinofuranosidases from the basidiomycete Ustilago maydis (UmAbf62A) and ascomycete Podospora anserina (PaAbf62A). Both enzymes are able to efficiently remove the α-L-arabinosyl substituents from arabinoxylan. The overall three-dimensional structure of UmAbf62A and PaAbf62A reveals a five-bladed β-propeller fold that confirms their predicted classification into clan GH-F together with GH43 α-L-arabinofuranosidases. Crystallographic structures of the complexes with arabinose and cellotriose reveal the important role of subsites +1 and +2 for sugar binding. Intriguingly, we observed that PaAbf62A was inhibited by cello-oligosaccharides and displayed binding affinity to cellulose although no activity was observed on a range of cellulosic substrates. Bioinformatic analyses showed that UmAbf62A and PaAbf62A belong to two distinct subfamilies within the GH62 family. The results presented here provide a framework to better investigate the structure-function relationships within the GH62 family.

Using a grass substrate to compare decay among two clades of brown rot fungi

Interest in the mechanisms of wood-degrading fungi has grown in tandem with lignocellulose bioconversion efforts, yet many potential biomass feedstocks are non-woody. Using corn stover (Zea mays) as a substrate, we tracked degradative capacities among brown rot fungi from the Antrodia clade, including Postia placenta, the first brown rot fungus to have its genome sequenced. Decay dynamics were compared against Gloeophyllum trabeum from the Gloeophyllum clade. Weight loss induced by P. placenta (6.2 %) and five other Antrodia clade isolates (average 7.4 %) on corn stalk after 12 weeks demonstrated inefficiency among these fungi, relative to decay induced by G. trabeum (44.4 %). Using aspen (Populus sp.) as a woody substrate resulted in, on average, a fourfold increase in weight loss induced by Antrodia clade fungi, while G. trabeum results matched those on stover. The sequence and trajectories of chemical constituent losses differed as a function of substrate but not fungal clade. Instead, chemical data suggest that characters unique to stover limit decay by the Antrodia clade, rather than disparities in growth rate or extractives toxicity. High p-coumaryl lignin content, lacking the methoxy groups characteristically cleaved during brown rot, is among potential rate-distinguishing characters in grasses. This ineptitude among Antrodia clade fungi on grasses was supported by meta-analysis of other unrelated studies using grass substrates. Concerning application, results expose a problem if adopting the strategy of the model decay fungus P. placenta to treat corn stover, a widely available plant feedstock. Overall, the results insinuate phylogenetically distinct modes of brown rot and demonstrate the benefit of using non-woody substrates to probe wood degradation mechanisms.