Deep insight into the Ganoderma lucidum by comprehensive analysis of its transcriptome

Combined NMR and MS-based metabonomics and real-time PCR analyses reveal dynamic metabolic changes of Ganoderma lucidum during fruiting body growing

Ganoderma lucidum (G. lucidum) is a rare medicinal fungus with various beneficial properties. One of its main components, ganoderic acids (GAs), are important triterpenoids known for their sedative and analgesic, hepatoprotective, and anti-tumor activities. Understanding the growth and development of the G. lucidum fruiting body is crucial for determining the optimal time to harvest them. In this study, we used nuclear magnetic resonance (NMR) spectroscopy to systematically characterize the metabolites of G. lucidum at seven distinct developmental stages. We also measured the contents of seven kinds of GAs using LC-MS/MS. A total of 49 metabolites were detected in G. lucidum, including amino acids, sugars, organic acids and GAs. During the transition from the bud development period (I) to the budding period (II), we observed a rapid accumulation of glucose, tyrosine, nicotinamide ribotide, inosine and GAs. After the budding period, the contents of most metabolites decreased until the mature period (VII). In addition, the contents of GAs showed an initial raising, followed by a decline during the elongation period, except for GAF, which exhibited a rapid raise during the mature stage. We also detected the expression of several genes involved in GA synthesis, finding that most genes including 16 cytochrome P450 monooxygenase were all down-regulated during periods IV and VII compared to period I. These findings provide valuable insights into the dynamic metabolic profiles of G. lucidum throughout its growth stage, and it is recommended to harvest G. lucidum at period IV.

Comparative transcriptome analysis on candidate genes associated with fruiting body growth and development in Lyophyllum decastes

Lyophyllum decastes is a mushroom that is highly regarded for its culinary and medicinal properties. Its delectable taste and texture make it a popular choice for consumption. To gain a deeper understanding of the molecular mechanisms involved in the development of the fruiting body of L. decastes, we used RNA sequencing to conduct a comparative transcriptome analysis. The analysis encompassed various developmental stages, including the vegetative mycelium, primordial initiation, young fruiting body, medium-size fruiting body, and mature fruiting body stages. A range of 40.1 to 60.6 million clean reads were obtained, and de novo assembly generated 15,451 unigenes with an average length of 1,462.68 bp. Functional annotation of transcriptomes matched 76.84% of the unigenes to known proteins available in at least one database. The gene expression analysis revealed a significant number of differentially expressed genes (DEGs) between each stage. These genes were annotated and subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Highly differentially expressed unigenes were also identified, including those that encode extracellular enzymes, transcription factors, and signaling pathways. The accuracy of the RNA-Seq and DEG analyses was validated using quantitative PCR. Enzyme activity analysis experiments demonstrated that the extracellular enzymes exhibited significant differences across different developmental stages. This study provides valuable insights into the molecular mechanisms that underlie the development of the fruiting body in L. decastes.

Comparative analysis of proteomes and transcriptomes revealed the molecular mechanism of development and nutrition of Pleurotus giganteus at different fruiting body development stages

Pleurotus giganteus is a commercially cultivated high-temperature mushroom. Investigating the molecular mechanism of fruiting body development will help us to better understand the regulation of substrates and energy in this process. However, little information has been reported on the development and nutrients of the P. giganteus fruiting body. In the present study, P. giganteus is cultivated in a climate chamber, and comparative transcriptome, proteome, and nutritional analysis of P. giganteus fruiting bodies were performed. Our results revealed that Cytochrome P450 monooxygenases and hydrophobin proteins play important roles during the differentiation in the elongation stage. Later, carbon metabolism dominate the fruiting body metabolism and genes related to the carbohydrate metabolic process, glycolytic process, and gluconeogenesis were up-regulated in the mature fruiting bodies. The up-regulation of carbohydrate substrates utilization CAZymes genes and inconsistent protein expression in pileus indicated a reverse transportation of mRNA from the fruiting body to vegetative mycelia. In addition, protein concentration in the pileus is higher than that in the stem, while the stem is the major nitrogen metabolic and amino acid synthetic location. The integrated transcriptomic, proteomic, and nutritional analysis indicated a two-way transportation of substrates and mRNAs in P. giganteus. Stem synthesizes amino acids and transported them to pileus with reducing sugars, while pileus induces the expression of substrate degradation mRNA according to the needs of growth and development and transports them in the other direction.

Metabolic Processes and Biological Macromolecules Defined the Positive Effects of Protein-Rich Biostimulants on Sugar Beet Plant Development

Protein-based biostimulants (PBBs) have a positive effect on plant development, although the biological background for this effect is not well understood. Here, hydrolyzed wheat gluten (HWG) and potato protein film (PF) in two levels (1 and 2 g/kg soil) and in two different soils (low and high nutrient; LNC and HNC) were used as PBBs. The effect of these PBBs on agronomic traits, sugars, protein, and peptides, as well as metabolic processes, were evaluated on sugar beet in comparison with no treatment (control) and treatment with nutrient solution (NS). The results showed a significant growth enhancement of the plants using HWG and PF across the two soils. Sucrose and total sugar content in the roots were high in NS-treated plants and correlated to root growth in HNC soil. Traits related to protein composition, including nitrogen, peptide, and RuBisCO contents, were enhanced in PBB-treated plants (mostly for HWG and PF at 2 g/kg soil) by 100% and >250% in HNC and LNC, respectively, compared to control. The transcriptomic analysis revealed that genes associated with ribosomes and photosynthesis were upregulated in the leaf samples of plants treated with either HWG or PP compared to the control. Furthermore, genes associated with the biosynthesis of secondary metabolites were largely down-regulated in root samples of HWG or PF-treated plants. Thus, the PBBs enhanced protein-related traits in the plants through a higher transcription rate of genes related to protein- and photosynthesis, which resulted in increased plant growth, especially when added in certain amounts (2 g/kg soil). However, sucrose accumulation in the roots of sugar beet seemed to be related to the easy availability of nitrogen.

Transcriptome Analysis of Antrodia cinnamomea Mycelia from Different Wood Substrates

Antrodia cinnamomea, an edible and medicinal fungus with significant economic value and application prospects, is rich in terpenoids, benzenoids, lignans, polysaccharides, and benzoquinone, succinic and maleic derivatives. In this study, the transcriptome of A. cinnamomea cultured on the wood substrates of Cinnamomum glanduliferum (YZM), C. camphora (XZM), and C. kanehirae (NZM) was sequenced using the high-throughput sequencing technology Illumina HiSeq 2000, and the data were assembled by de novo strategy to obtain 78,729 Unigenes with an N50 of 4,463 bp. Compared with public databases, about 11,435, 6,947, and 5,994 Unigenes were annotated to the Non-Redundant (NR), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genome (KEGG), respectively. The comprehensive analysis of the mycelium terpene biosynthesis-related genes in A. cinnamomea revealed that the expression of acetyl-CoA acetyltransferase (AACT), acyl-CoA dehydrogenase (MCAD), 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA), mevalonate pyrophosphate decarboxylase (MVD), and isopentenyl diphosphate isomerase (IDI) was significantly higher on NZM compared to the other two wood substrates. Similarly, the expression of geranylgeranyltransferase (GGT) was significantly higher on YZM compared to NZM and XZM, and the expression of farnesyl transferase (FTase) was significantly higher on XZM. Furthermore, the expressions of 2,3-oxidized squalene cyclase (OCS), squalene synthase (SQS), and squalene epoxidase (SE) were significantly higher on NZM. Overall, this study provides a potential approach to explore the molecular regulation mechanism of terpenoid biosynthesis in A. cinnamomea.

Increasing the production of the bioactive compounds in medicinal mushrooms: an omics perspective

Macroscopic fungi, mainly higher basidiomycetes and some ascomycetes, are considered medicinal mushrooms and have long been used in different areas due to their pharmaceutically/nutritionally valuable bioactive compounds. However, the low production of these bioactive metabolites considerably limits the utilization of medicinal mushrooms both in commerce and clinical trials. As a result, many attempts, ranging from conventional methods to novel approaches, have been made to improve their production. The novel strategies include conducting omics investigations, constructing genome-scale metabolic models, and metabolic engineering. So far, genomics and the combined use of different omics studies are the most utilized omics analyses in medicinal mushroom research (both with 31% contribution), while metabolomics (with 4% contribution) is the least. This article is the first attempt for reviewing omics investigations in medicinal mushrooms with the ultimate aim of bioactive compound overproduction. In this regard, the role of these studies and systems biology in elucidating biosynthetic pathways of bioactive compounds and their contribution to metabolic engineering will be highlighted. Also, limitations of omics investigations and strategies for overcoming them will be provided in order to facilitate the overproduction of valuable bioactive metabolites in these valuable organisms.

Physiological changes and gene responses during Ganoderma lucidum growth with selenium supplementation

Ganoderma lucidum basidiomycota is highly appreciated for its health and nutrition value. In the present study, Ganoderma lucidum was cultivated as selenium transformation carrier, and the physiological changes and gene responses by selenium supplementation were revealed through high-throughput RNA-Seq technology. As a result, selenium supplementation increased the stipe length and the cap size, but decreased the cap thickness of G. lucidum. Mineral salt supplementation could greatly promote the formation of triterpene acids and selenium in G. lucidum. The highest yield was gained in the treatment with selenium content of 200 µg/g. Subsequently, the tissues of G. lucidum at budding and mature stages in this treatment group were sampled for transcriptome analysis and compared to those of a control group without selenium supplementation. A total of 16,113 expressed genes were obtained from the transcriptome of G. lucidum, and GO-annotated unigenes were mainly involved in molecular functions and KEGG-annotated ones were highly expressed in ribosomal pathway. Furthermore, genes involved in carbon metabolism pathway were most promoted by selenium at budding stage of G. lucidum, while gene expression was the highest in the pathway of amino acid biosynthesis at mature stage of G. lucidum. Specially, selenium-related genes in G. lucidum, such as GL23172-G, GL29881-G and GL28298-G, played a regulatory role in oxidoreductase, antioxidant activity and tryptophan synthesis. The results provide a theoretical basis for further study of selenium-enriched mushrooms and aid to development of Se-enriched foodstuff and health products made from fungi.

Comparative transcriptome and WGCNA reveal key genes involved in lignocellulose degradation in Sarcomyxa edulis

The developmental transcriptomes of Sarcomyxa edulis were assessed to explore the molecular mechanisms underlying lignocellulose degradation. Six stages were analyzed, spanning the entire developmental process: growth of mycelium until occupying half the bag (B1), mycelium under low-temperature stimulation after occupying the entire bag (B2), appearance of mycelium in primordia (B3), primordia (B4), mycelium at the harvest stage (B5), and mature fruiting body (B6). Samples from all six developmental stages were used for transcriptome sequencing, with three biological replicates for all experiments. A co-expression network of weighted genes associated with extracellular enzyme physiological traits was constructed using weighted gene co-expression network analysis (WGCNA). We obtained 19 gene co-expression modules significantly associated with lignocellulose degradation. In addition, 12 key genes and 8 kinds of TF families involved in lignocellulose degradation pathways were discovered from the four modules that exhibited the highest correlation with the target traits. These results provide new insights that advance our understanding of the molecular genetic mechanisms of lignocellulose degradation in S. edulis to facilitate its utilization by the edible mushroom industry.

De novo transcriptome assembly and comprehensive assessment provide insight into fruiting body formation of Sparassis latifolia

The genes associated with fruiting body formation of Sparasis latifolia are valuable for improving mushroom breeding. To investigate this process, 4.8 × 108 RNA-Seq reads were acquired from three stages: hyphal knot (SM), primordium (SP), and primordium differentiation (SPD). The de novo assembly generated a total of 48,549 unigenes, of which 71.53% (34,728) unigenes could be annotated by at least one of the KEGG (Kyoto Encyclopedia of Genes and Genomes), GO (Gene Ontology), and KOG (Eukaryotic Orthologous Group) databases. KEGG and KOG analyses respectively mapped 32,765 unigenes to 202 pathways and 19,408 unigenes to 25 categories. KEGG pathway enrichment analysis of DEGs (differentially expressed genes) indicated primordium initiation was significantly related to 66 pathways, such as "Ribosome", "metabolism of xenobiotics by cytochrome P450", and "glutathione metabolism" (among others). The MAPK and mTOR signal transduction pathways underwent significant adjustments during the SM to SP transition. Further, our research revealed the PI3K-Akt signaling pathway related to cell proliferation could play crucial functions during the development of SP and SPD. These findings provide crucial candidate genes and pathways related to primordium differentiation and development in S. latifolia, and advances our knowledge about mushroom morphogenesis.

Transcriptome and Differentially Expressed Gene Profiles in Mycelium, Primordium and Fruiting Body Development in Stropharia rugosoannulata

Stropharia rugosoannulata uses straw as a growth substrate during artificial cultivation and has been widely promoted in China. However, its fruiting body formation and development processes have not been elucidated. In this study, the developmental transcriptomes were analyzed at three stages: the mycelium (G-S), primordium (P-S) and fruiting body (M-F) stages. A total of 9690 differentially expressed genes (DEGs) were identified in the different developmental stages. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that these DEGs were involved mainly in hydrolase activity, structural molecule activity and oxidoreductase activity as well as xenobiotic biodegradation and metabolism and energy metabolism pathways. We further found that the higher expression of most carbohydrate enzyme (i.e., GH, CE, CBM, AA and PL) genes in the hyphal (i.e., G-S) stage was related mainly to substrate degradation, while the upregulation of glycosyltransferase (GT) gene expression in the P-S and M-F stages may be related to cell wall synthesis. In addition, we found that CO₂-sensing-related genes (i.e., CA-2, CA-3, PKA-1 and PKA-2) were upregulated in the P-S and M-F stages, heat shock protein genes (HSP60 and HSP90) were significantly downregulated in the P-S stage and upregulated in the M-F stage and the transcription factors (i.e., steA, MYB, nosA, HAP1, and GATA-4/5/6) involved in growth and development were significantly upregulated in the P-S stage. These results suggest that environmental factors (i.e., CO₂ and temperature) and transcription factors may play a key role in primordium formation. In short, this study provides new insights into the study of stimulating primordia formation affecting the development of fruiting bodies of S. rugosoannulata.

Efficient expression of heterologous genes by the introduction of the endogenous glyceraldehyde-3-phosphate dehydrogenase gene intron 1 in Ganoderma lucidum

Background

Ganoderma lucidum, a well-known medicinal mushroom, has received wide attention as a promising cell factory for producing bioactive compounds. However, efficient expression of heterologous genes remains a major challenge in Ganoderma, hindering metabolic regulation research and molecular breeding of this species.

Results

We show that the presence of glyceraldehyde-3-phosphate dehydrogenase gene (gpd) intron 1 at the 5′ end of, the 3′ end of, or within the heterologous phosphinothricin-resistant gene (bar) is efficient for its expression in G. lucidum. The enhanced expression of bar is exhibited by the higher accumulation of mRNA and increased amounts of protein. Moreover, the insertion of the gpd intron 1 in the β-glucuronidase gene (gus) elevates its mRNA accumulation and enzyme activity, which facilitates the use of this reporter gene in Ganoderma.

Conclusions

This study has demonstrated the importance of the introduction of gpd intron 1 for the efficient expression of bar and gus in G. lucidum. The presence of the gpd intron 1 in heterologous genes increases levels of mRNA accumulation and protein expression in basidiomycete Ganoderma. The developed method may be utilized in upregulating the expression of other heterologous genes in Ganoderma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01654-8.

Transcriptome of different fruiting stages in the cultivated mushroom Cyclocybe aegerita suggests a complex regulation of fruiting and reveals enzymes putatively involved in fungal oxylipin biosynthesis

BACKGROUND

Cyclocybe aegerita (syn. Agrocybe aegerita) is a commercially cultivated mushroom. Its archetypal agaric morphology and its ability to undergo its whole life cycle under laboratory conditions makes this fungus a well-suited model for studying fruiting body (basidiome, basidiocarp) development. To elucidate the so far barely understood biosynthesis of fungal volatiles, alterations in the transcriptome during different developmental stages of C. aegerita were analyzed and combined with changes in the volatile profile during its different fruiting stages.

RESULTS

A transcriptomic study at seven points in time during fruiting body development of C. aegerita with seven mycelial and five fruiting body stages was conducted. Differential gene expression was observed for genes involved in fungal fruiting body formation showing interesting transcriptional patterns and correlations of these fruiting-related genes with the developmental stages. Combining transcriptome and volatilome data, enzymes putatively involved in the biosynthesis of C8 oxylipins in C. aegerita including lipoxygenases (LOXs), dioxygenases (DOXs), hydroperoxide lyases (HPLs), alcohol dehydrogenases (ADHs) and ene-reductases could be identified. Furthermore, we were able to localize the mycelium as the main source for sesquiterpenes predominant during sporulation in the headspace of C. aegerita cultures. In contrast, changes in the C8 profile detected in late stages of development are probably due to the activity of enzymes located in the fruiting bodies.

CONCLUSIONS

In this study, the combination of volatilome and transcriptome data of C. aegerita revealed interesting candidates both for functional genetics-based analysis of fruiting-related genes and for prospective enzyme characterization studies to further elucidate the so far barely understood biosynthesis of fungal C8 oxylipins.

Transcriptional Dynamics of Genes Purportedly Involved in the Control of Meiosis, Carbohydrate, and Secondary Metabolism during Sporulation in Ganoderma lucidum

Ganoderma lucidum spores (GLS), the mature germ cells ejected from the abaxial side of the pileus, have diverse pharmacological effects. However, the genetic regulation of sporulation in this fungus remains unknown. Here, samples corresponding to the abaxial side of the pileus were collected from strain YW-1 at three sequential developmental stages and were then subjected to a transcriptome assay. We identified 1598 differentially expressed genes (DEGs) and found that the genes related to carbohydrate metabolism were strongly expressed during spore morphogenesis. In particular, genes involved in trehalose and malate synthesis were upregulated, implying the accumulation of specific carbohydrates in mature G. lucidum spores. Furthermore, the expression of genes involved in triterpenoid and ergosterol biosynthesis was high in the young fruiting body but gradually decreased with sporulation. Finally, spore development-related regulatory pathways were explored by analyzing the DNA binding motifs of 24 transcription factors that are considered to participate in the control of sporulation. Our results provide a dataset of dynamic gene expression during sporulation in G. lucidum. They also shed light on genes potentially involved in transcriptional regulation of the meiotic process, metabolism pathways in energy provision, and ganoderic acids and ergosterol biosynthesis.

A putative transcription factor LFC1 negatively regulates development and yield of winter mushroom

Basidioma is the fruiting body of mushroom species. The deep understanding on the mechanism of basidioma development is valuable for mushroom breeding and cultivation. From winter mushroom (Flammulina velutipes), one of the top five industrially cultivated mushrooms, a novel putative Zn(II)2Cys6 transcription factor LFC1 with negative regulatory function in basidioma development was identified. The transcript level of lfc1 was dramatically decreased during basidioma development. Neither overexpression nor knockdown of lfc1 affected hyphal vegetative growth. However, knockdown of lfc1 could promote basidioma development and shorten cultivation time by 2 days, while overexpression of lfc1 delayed the optimal harvest time by 3 days. In the lfc1 knockdown strain, in which the lfc1 expression was reduced by 72%, mushroom yield and biological efficiency could be increased at least by 24%. Knockdown of lfc1 did not affect the shape of caps but significantly increased basidioma length and number, while its overexpression did not affect basidioma length but dramatically reduced basidioma number. In addition, rather than producing basidiomata with round caps as in wild type, the caps of basidiomata in the lfc1 overexpression mutants were significantly larger and the cap edge was wrinkled. RNA-seq analysis revealed that 455 genes had opposite transcriptional responses to lfc1 overexpression and knockdown. Some of them were previously reported as genes involved in basidioma development, including 3 hydrophobin encoding genes, 2 lectin encoding genes, FVFD16, an Eln2 ortholog encoding gene, and 3 genes encoding membrane components. As LFC1 homologs are widely present in mushroom species, lfc1 can be useful in mushroom breeding.Key Points• A novel transcription factor LFC1 negatively regulates fruiting in winter mushroom• LFC1 regulated transcription of more than 400 genes.• Reduction of LFC1 expression could shorten cultivation time and increase yield.• lfc1 could be a potentially useful reference gene for mushroom breeding.

A WD40 Protein Encoding Gene Fvcpc2 Positively Regulates Mushroom Development and Yield in Flammulina velutipes

Ascomycota and Basidiomycota are two closely related phyla and fungi in two phyla share some common morphological developmental process during fruiting body formation. In Neurospora crassa, the Gβ-like protein CPC-2 with a seven-WD40 repeat domain was previously reported. By transforming CPC-2 ortholog encoding genes, from 7 different fungal species across Ascomycota and Basidiomycota, into the cpc-2 deletion mutant of N. crassa, we demonstrate that all tested CPC-2 ortholog genes were able to complement the defects of the cpc-2 deletion mutant in sexual development, indicating that CPC-2 proteins from Ascomycota and Basidiomycota have the similar cellular function. Using Flammulina velutipes as a model system for mushroom species, the CPC-2 ortholog FvCPC2 was characterized. Fvcpc2 increased transcription during fruiting body development. Knockdown of Fvcpc2 by RNAi completely impaired fruiting body formation. In three Fvcpc2 knockdown mutants, transcriptional levels of genes encoding adenylate cyclase and protein kinase A catalytic subunit were significantly lower and colony growth became slower than wild type. The addition of cAMP or the PKA-activator 8-Bromo-cAMP into the medium restored the Fvcpc2 knockdown mutants to the wild-type colony growth phenotype, suggesting that the involvement of cAMP production in the regulatory mechanisms of FvCPC2. Knockdown of Fvcpc2 also weakened transcriptional responses to sexual development induction by some genes related to fruiting body development, including 4 jacalin-related lectin encoding genes, 4 hydrophobin encoding genes, and 3 functionally-unknown genes, suggesting the participation of these genes in the mechanisms by which FvCPC2 regulates fruiting body development. All three Fvcpc2 overexpression strains displayed increased mushroom yield and shortened cultivation time compared to wild type, suggesting that Fvcpc2 can be a promising reference gene for Winter Mushroom breeding. Since the orthologs of FvCPC2 were highly conserved and specifically expressed during fruiting body development in different edible mushrooms, genes encoding FvCPC2 orthologs in other mushroom species also have potential application in breeding.

Enabling community-based metrology for wood-degrading fungi

BACKGROUND

Lignocellulosic biomass could support a greatly-expanded bioeconomy. Current strategies for using biomass typically rely on single-cell organisms and extensive ancillary equipment to produce precursors for downstream manufacturing processes. Alternative forms of bioproduction based on solid-state fermentation and wood-degrading fungi could enable more direct means of manufacture. However, basic methods for cultivating wood-degrading fungi are often ad hoc and not readily reproducible. Here, we developed standard reference strains, substrates, measurements, and methods sufficient to begin to enable reliable reuse of mycological materials and products in simple laboratory settings.

RESULTS

We show that a widely-available and globally-regularized consumer product (Pringles™) can support the growth of wood-degrading fungi, and that growth on Pringles™-broth can be correlated with growth on media made from a fully-traceable and compositionally characterized substrate (National Institute of Standards and Technology Reference Material 8492 Eastern Cottonwood Whole Biomass Feedstock). We also establish a Relative Extension Unit (REU) framework that is designed to reduce variation in quantification of radial growth measurements. So enabled, we demonstrate that five laboratories were able to compare measurements of wood-fungus performance via a simple radial extension growth rate assay, and that our REU-based approach reduced variation in reported measurements by up to ~ 75%.

CONCLUSIONS

Reliable reuse of materials, measures, and methods is necessary to enable distributed bioproduction processes that can be adopted at all scales, from local to industrial. Our community-based measurement methods incentivize practitioners to coordinate the reuse of standard materials, methods, strains, and to share information supporting work with wood-degrading fungi.

Insight into plant cell wall degradation and pathogenesis of Ganoderma boninense via comparative genome analysis

Background

G. boninense is a hemibiotrophic fungus that infects oil palms (Elaeis guineensis Jacq.) causing basal stem rot (BSR) disease and consequent massive economic losses to the oil palm industry. The pathogenicity of this white-rot fungus has been associated with cell wall degrading enzymes (CWDEs) released during saprophytic and necrotrophic stage of infection of the oil palm host. However, there is a lack of information available on the essentiality of CWDEs in wood-decaying process and pathogenesis of this oil palm pathogen especially at molecular and genome levels.

Methods

In this study, comparative genome analysis was carried out using the G. boninense NJ3 genome to identify and characterize carbohydrate-active enzyme (CAZymes) including CWDE in the fungal genome. Augustus pipeline was employed for gene identification in G. boninense NJ3 and the produced protein sequences were analyzed via dbCAN pipeline and PhiBase 4.5 database annotation for CAZymes and plant-host interaction (PHI) gene analysis, respectively. Comparison of CAZymes from G. boninense NJ3 was made against G. lucidum, a well-studied model Ganoderma sp. and five selected pathogenic fungi for CAZymes characterization. Functional annotation of PHI genes was carried out using Web Gene Ontology Annotation Plot (WEGO) and was used for selecting candidate PHI genes related to cell wall degradation of G. boninense NJ3.

Results

G. boninense was enriched with CAZymes and CWDEs in a similar fashion to G. lucidum that corroborate with the lignocellulolytic abilities of both closely-related fungal strains. The role of polysaccharide and cell wall degrading enzymes in the hemibiotrophic mode of infection of G. boninense was investigated by analyzing the fungal CAZymes with necrotrophic Armillaria solidipes, A. mellea, biotrophic Ustilago maydis, Melampsora larici-populina and hemibiotrophic Moniliophthora perniciosa. Profiles of the selected pathogenic fungi demonstrated that necrotizing pathogens including G. boninense NJ3 exhibited an extensive set of CAZymes as compared to the more CAZymes-limited biotrophic pathogens. Following PHI analysis, several candidate genes including polygalacturonase, endo β-1,3-xylanase, β-glucanase and laccase were identified as potential CWDEs that contribute to the plant host interaction and pathogenesis.

Discussion

This study employed bioinformatics tools for providing a greater understanding of the biological mechanisms underlying the production of CAZymes in G. boninense NJ3. Identification and profiling of the fungal polysaccharide- and lignocellulosic-degrading enzymes would further facilitate in elucidating the infection mechanisms through the production of CWDEs by G. boninense. Identification of CAZymes and CWDE-related PHI genes in G. boninense would serve as the basis for functional studies of genes associated with the fungal virulence and pathogenicity using systems biology and genetic engineering approaches.

De novo transcriptome and proteome analysis of Dictyophora indusiata fruiting bodies provides insights into the changes during morphological development

De novo transcriptome assembly and shotgun proteome analysis of Dictyophora indusiata fruiting bodies were performed. A total of 19,704 unigenes were sequenced, and 4380 proteins were identified. Annotation and functional analysis of the identified proteins were significantly enriched in small molecule synthetic and metabolic processes, protein modification regulation (phosphorylation and ubiquitination), and vesicle transport. Furthermore, quantitative developmental transcriptome analysis was performed between the peach-shaped and mature fruiting bodies, and the results revealed that the metabolism and transport activities were upregulated in the mature stage, while protein translation was downregulated; this regulation is likely the main reason for the significant changes in the nutrients of fruiting bodies. Furthermore, the cell wall stress-dependent MAPK sub-pathway was activated in the mature stage, and fungal cell wall degradation-related genes were upregulated, which could promote reconstruction of the cell wall and might play a key role in the morphological development of D. indusiata fruiting bodies.

A Comparative Transcriptome Analysis Reveals Physiological Maturation Properties of Mycelia in Pleurotus tuoliensis

Pleurotus tuoliensis is a precious edible fungus with extremely high nutritive and medicinal value. The cultivation period of P. tuoliensis is longer than those of other Pleurotus species, which is mainly due to a longer mycelium physiological maturation period (30–60 days). Currently, the molecular processes underlying physiological maturation of the mycelium remain unclear. We performed a comparative transcriptomic analysis of immature and mature mycelia using RNA-seq. De novo transcriptome assembly resulted in identification of 17,030 unigenes. 451 differentially expressed genes—including those encoding nucleoside diphosphate kinase (NDPK), glycoside hydrolase family proteins, exopolygalacturonase, and versatile peroxidases—were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that nucleotide synthesis and energy metabolism are highly active during the physiological maturation of mycelia, and genes related to these pathways were significantly upregulated in mature mycelia. NDPK is predicted to be essential for mycelia maturation. Our findings contribute to a comprehensive understanding of mycelia maturation in a commercially important fungal species. Future efforts will focus on the function of NDPK and the mechanism by which it regulates mycelia maturation.

Transcriptional profiling of Auricularia cornea in selenium accumulation

Auricularia cornea is a widely cultivated edible fungus with substantial nutritive value. This study aimed to enrich the multifunctional bionutrient element selenium in A. cornea to improve its quality and explore the accumulation of selenium in the fungus using high-throughput RNA-Seq technology. In general, the treatment group with a 100 µg/g supply of selenium outperformed the other treatment groups in terms of high yield, rich crude polysaccharides and a high total selenium concentration. Additional evidences demonstrated the budding and mature phases were two typical growth stages of A. cornea and were important for the accumulation of selenium. Therefore, the budding and mature phase tissues of A. cornea in the treatment group with a 100 µg/g supply of selenium were used for transcriptome analysis and compared to those of a control group that lacked additional selenium. A total of 2.56 × 105 unigenes from A. cornea transcriptome were assembled and annotated to five frequently used databases including NR, GO, KEGG, eggNOG and SwissProt. GO and KEGG pathway analysis revealed that genes involved in metabolic process and translation were up-expressed at the budding stage in response to selenium supplementation, including amino acid metabolism, lipid metabolism, ribosome. In addition, the differential gene expression patterns of A. cornea suggested that the up-expressed genes were more likely to be detected at the budding stage than at the mature stage. These results provide insights into the transcriptional response of A. cornea to selenium accumulation.

RNA Sequencing Reveals Differential Gene Expression of Cerrena Unicolor in Response to Variable Lighting Conditions

To elucidate the light-dependent gene expression in Cerrena unicolor FCL139, the transcriptomes of the fungus growing in white, blue, green, and red lighting conditions and darkness were analysed. Among 10,413 all-unigenes detected in C. unicolor, 7762 were found to be expressed in all tested conditions. Transcripts encoding putative fungal photoreceptors in the C. unicolor transcriptome were identified. The number of transcripts uniquely produced by fungus ranged from 20 during its growth in darkness to 112 in the green lighting conditions. We identified numerous genes whose expression differed substantially between the darkness (control) and each of the light variants tested, with the greatest number of differentially expressed genes (DEGs) (454 up- and 457 down-regulated) observed for the white lighting conditions. The DEGs comprised those involved in primary carbohydrate metabolism, amino acid metabolism, autophagy, nucleotide repair systems, signalling pathways, and carotenoid metabolism as defined using Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The analysis of the expression profile of genes coding for lignocellulose-degrading enzymes suggests that the wood-degradation properties of C. unicolor may be independent of the lighting conditions and may result from the overall stimulation of fungal metabolism by daylight.

Comparative transcriptome analysis reveals relationship of three major domesticated varieties of Auricularia auricula-judae

Auricularia auricula-judae is an edible mushroom and a traditional medicine in China as well as the fourth largest cultivated mushroom species in the world. Here for the first time, we present comparative transcriptome analyses of the fruiting bodies of three morphologically distinguishable A. auricula-judae cultivated varieties (Wujin, smooth; Banjin, partially wrinkled; and Quanjin, fully wrinkled) collected from Jilin Province, China. Biological triplicates were performed to determine the expression levels of 13,937 unigenes. Among them, only 13 unigenes were annotated to A. auricula-judae, highlighting the lack of publicly available reference sequences for this economically important species. Principal component analysis (PCA) determined that the gene expression profile of Quanjin was unique when compared to those of Banjin and Wujin. Such relationships were further supported by analyses of annotated and unannotated unigenes, differentially expressed unigenes, gene ontology functions, and the family of peroxidase genes. Using the KEGG database, significant alternations in biological pathways were detected among the three cultivars. This work contributes a large set of A. auricula-judae sequences to public database, establishes the relationships among major cultivars, and provides molecular guidance for breeding and cultivation.

Comparative transcriptome analysis of dikaryotic mycelia and mature fruiting bodies in the edible mushroom Lentinula edodes

Lentinula edodes is a popular cultivated edible mushroom with high nutritional and medicinal value. To understand the regulation of gene expression in the dikaryotic mycelium and mature fruiting body in the commercially important Korean L. edodes strain, we first performed comparative transcriptomic analysis, using Illumina HiSeq platform. De novo assembly of these sequences revealed 11,675 representative transcripts in two different stages of L. edodes. A total of 9,092 unigenes were annotated and subjected to Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Gene expression analysis revealed that 2,080 genes were differentially expressed, with 1,503 and 577 upregulated in the mycelium and a mature fruiting body, respectively. Analysis of 18 KEGG categories indicated that fruiting body-specific transcripts were significantly enriched in ‘replication and repair’ and ‘transcription’ pathways, which are important for premeiotic replication, karyogamy, and meiosis during maturation. We also searched for fruiting body-specific proteins such as aspartic protease, gamma-glutamyl transpeptidase, and cyclohexanone monooxygenase, which are involved in fruiting body maturation and isolation of functional substances. These transcriptomes will be useful in elucidating the molecular mechanisms of mature fruiting body development and beneficial properties, and contribute to the characterization of novel genes in L. edodes.

Investigation of lignocellulolytic enzymes during different growth phases of Ganoderma lucidum strain G0119 using genomic, transcriptomic and secretomic analyses

Ganoderma lucidum is a medicinal mushroom that is well known for its ability to enhance human health, and products made from this fungus have been highly profitable. The substrate-degrading ability of G. lucidum could be related to its growth. CAZy proteins were more abundant in its genome than in the other white rot fungi models. Among these CAZy proteins, changes in lignocellulolytic enzymes during growth have not been well studied. Using genomic, transcriptomic and secretomic analyses, this study focuses on the lignocellulolytic enzymes of G. lucidum strain G0119 to determine which of these degradative enzymes contribute to its growth. From the genome sequencing data, genes belonging to CAZy protein families, especially genes involved in lignocellulose degradation, were investigated. The gene expression, protein abundance and enzymatic activity of lignocellulolytic enzymes in mycelia over a growth cycle were analysed. The overall expression cellulase was higher than that of hemicellulase and lignin-modifying enzymes, particularly during the development of fruiting bodies. The cellulase and hemicellulase abundances and activities increased after the fruiting bodies matured, when basidiospores were produced in massive quantities till the end of the growth cycle. Additionally, the protein abundances of the lignin-modifying enzymes and the expression of their corresponding genes, including laccases and lignin-degrading heme peroxidases, were highest when the mycelia fully spread in the compost bag. Type I cellobiohydrolase was observed to be the most abundant extracellular lignocellulolytic enzyme produced by the G. lucidum strain G0119. The AA2 family haem peroxidases were the dominant lignin-modifying enzyme expressed during the mycelial growth phase, and several laccases might play roles during the formation of the primordium. This study provides insight into the changes in the lignocellulose degradation ability of G. lucidum during its growth and will facilitate the discovery of new approaches to accelerate the growth of G. lucidum in culture.

Comparative transcriptomic analysis of Cerrena unicolor revealed differential expression of genes engaged in degradation of various kinds of wood

To explore the number of enzymes engaged by Cerrena unicolor FCL139 for wood degradation, the transcriptomes of the fungus growing on birch, ash, maple sawdust and the control liquid medium were analyzed. Among 12,966 gene models predicted for the C. unicolor genome, 10,396 all-unigenes were detected, of which 9567 were found to be expressed in each of the tested growth media. The highest number (107) of unique transcripts was detected during fungus growth in the control liquid medium, while the lowest number (11) - in the fungal culture comprising maple saw dust. Analysis of C. unicolor transcriptomes identified numerous genes whose expression differed substantially between the mycelia growing in control medium and each of the sawdust media used, with the highest number (828) of upregulated transcripts observed during the fungus growth on the ash medium. Among the 294 genes that were potentially engaged in wood degradation, the expression of 59 was significantly (p < .01) changed in the tested conditions. The transcripts of 37 of those genes were at least four times more abundant in the cells grown in all sawdust media when compared to the control medium. Upregulated genes coding for cellulases and, to a lower extent, hemicellulases predominated during fungus growth on sawdust. Transcripts encoding cellulolytic enzymes were the most abundant in mycelia grown on birch and maple while lower number of such transcripts was detected in fungus growing on ash. The expression pattern of lignolytic activities-coding genes was strongly dependent on the type of sawdust applied for fungus growth medium.

De novo transcriptome assembly and characterization of the 10-hydroxycamptothecin-producing Xylaria sp. M71 following salicylic acid treatment

In the present study, we identified genes that are putatively involved in the production of fungal 10-hydroxycamptothecin via transcriptome sequencing and characterization of the Xylaria sp. M71 treated with salicylic acid (SA). A total of 60,664,200 raw reads were assembled into 26,044 unigenes. BLAST assigned 8,767 (33.7%) and 10,840 (41.6%) unigenes to 40 Gene Ontology (GO) annotations and 108 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 3,713 unigenes comprising 1,504 upregulated and 2,209 downregulated unigenes were found to be differentially expressed between SA-induced and control fungi. Based on the camptothecin biosynthesis pathway in plants, 13 functional genes of Xylaria sp. M71 were mapped to the mevalonate (MVA) pathway, suggesting that the fungal 10-hydroxycamptothecin is produced via the MVA pathway. In summary, analysis of the Xylaria sp. M71 transcriptome allowed the identification of unigenes that are putatively involved in 10-hydroxycamptothecin biosynthesis in fungi.

A Lentinus edodes polysaccharide induces mitochondrial-mediated apoptosis in human cervical carcinoma HeLa cells

In this study, a homogeneous polysaccharide (LEP1) with an average molecular weight of 53kDa was successfully purified from the fruiting bodies of Lentinus edodes and its anticancer efficacy on human cervical carcinoma HeLa cells in vitro and associated possible molecular mechanism were also evaluated. MTT assay showed that LEP1 exhibited a dose-dependent inhibitory effect on the proliferation of HeLa cells and caused apoptotic death. Our present findings provided the first evidence that LEP1 induced the apoptosis of HeLa cells via a mitochondria dependent pathway, as indicated by an increase in Bax/Bcl-2 ratio, a loss of mitochondrial membrane potential (Δym), the release of cytochrome c from the mitochondria to the cytosol, activation of caspase-9 and caspase-3, and cleavage of poly (ADP-ribose) polymerase (PARP) in HeLa cells. These combined results unequivocally indicated that the involvement of mitochondria-mediated signaling pathway in LEP1-induced apoptosis and strongly provided experimental evidence for the use of LEP1 as a potential therapeutic agent in the prevention and treatment of human cervical carcinoma.

The transcriptional landscape of basidiosporogenesis in mature Pisolithus microcarpus basidiocarp

Background

Pisolithus microcarpus (Cooke & Massee) G. Cunn is a gasteromycete that produces closed basidiocarps in symbiosis with eucalypts and acacias. The fungus produces a complex basidiocarp composed of peridioles at different developmental stages and an upper layer of basidiospores free of the hyphae and ready for wind dispersal upon the rupture of the basidiocarp pellis. During basidiosporogenesis, a process that takes place inside the basidiocarp peridioles, a conspicuous reserve of fatty acids is present throughout development. While several previous studies have described basidiosporogenesis inside peridioles, very little is known about gene expression changes that may occur during this part of the fungal life cycle. The objective of this work was to analyze gene transcription during peridiole and basidiospore development, while focusing specifically on cell cycle progression and lipid metabolism.

Results

Throughout different developmental stages of the peridioles we analyzed, 737 genes were regulated between adjacent compartments (>5 fold, FDR-corrected p-value < 0.05) corresponding to 3.49% of the genes present in the P. microcarpus genome. We identified three clusters among the regulated genes which showed differential expression between the peridiole developmental stages and the basidiospores. During peridiole development, transcripts for proteins involved in cellular processes, signaling, and information storage were detected, notably those for coding transcription factors, DNA polymerase subunits, DNA repair proteins, and genes involved in chromatin structure. For both internal embedded basidiospores (hereto referred to as “Internal spores”, IS) and external free basidiospores (hereto referred to as “Free spores”, FS), upregulated transcripts were found to involve primary metabolism, particularly fatty acid metabolism (FA). High expression of transcripts related to β-oxidation and the glyoxylate shunt indicated that fatty acids served as a major carbon source for basidiosporogenesis.

Conclusion

Our results show that basidiocarp formation in P. microcarpus involves a complex array of genes that are regulated throughout peridiole development. We identified waves of transcripts with coordinated regulation and identified transcription factors which may play a role in this regulation. This is the first work to describe gene expression patterns during basidiocarp formation in an ectomycorrhizal gasteromycete fungus and sheds light on genes that may play important roles in the developmental process.

Electronic supplementary material

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

Genetic engineering of Ganoderma lucidum for the efficient production of ganoderic acids

Ganoderma lucidum is a well-known traditional medicinal mushroom that produces ganoderic acids with numerous interesting bioactivities. Genetic engineering is an efficient approach to improve ganoderic acid biosynthesis. However, reliable genetic transformation methods and appropriate genetic manipulation strategies remain underdeveloped and thus should be enhanced. We previously established a homologous genetic transformation method for G. lucidum; we also applied the established method to perform the deregulated overexpression of a homologous 3-hydroxy-3-methyl-glutaryl coenzyme A reductase gene in G. lucidum. Engineered strains accumulated more ganoderic acids than wild-type strains. In this report, the genetic transformation systems of G. lucidum are described; current trends are also presented to improve ganoderic acid production through the genetic manipulation of G. lucidum.

Tricking Arthrinium malaysianum into Producing Industrially Important Enzymes Under 2-Deoxy D-Glucose Treatment

This study catalogs production of industrially important enzymes and changes in transcript expression caused by 2-deoxy D-glucose (2-DG) treatment in Arthrinium malaysianum cultures. Carbon Catabolite Repression (CCR) induced by 2-DG in this species is cAMP independent unlike many other organisms. Higher levels of secreted endoglucanase (EG), β-glucosidase (BGL), β-xylosidase (BXL), and filter paper activity assay (FPase) enzymes under 2-DG treatment can be exploited for commercial purposes. An integrated RNA sequencing and quantitative proteomic analysis was performed to investigate the cellular response to 2-DG in A. malaysianum. Analysis of RNASeq data under 2-DG treated and control condition reveals that 56% of the unigenes do not have any known similarity to proteins in non-redundant database. Gene Ontology IDs were assigned to 36% of the transcripts (13260) and about 5207 (14%) were mapped to Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). About 1711 genes encoding 2691 transcripts were differentially expressed in treated vs. control samples. Out of the 2691 differentially expressed transcripts, only 582 have any known function. The most up regulated genes belonged to Pentose Phosphate Pathways and carbohydrate degradation class as expected. In addition, genes involved in protein folding, binding, catalytic activity, DNA repair, and secondary metabolites were up-regulated under 2-DG treatment. Whereas genes encoding glycosylation pathways, growth, nutrient reservoir activity was repressed. Gene ontology analysis of the differentially expressed genes indicates metabolic process (35%) is the pre-dominant class followed by carbohydrate degradation (11%), protein folding, and trafficking (6.2%) and transport (5.3%) classes. Unlike other organisms, conventional unfolded protein response (UPR) was not activated in either control or treated conditions. Major enzymes secreted by A. malaysianum are those degrading plant polysaccharides, the most dominant ones being β-glucosidase, as demonstrated by the 2D gel analysis. A set of 7 differentially expressed mRNAs were validated by qPCR. Transmission electron microscopy analyses demonstrated that the 2-DG treated cell walls of hyphae showed significant differences in the cell-wall thickness. Overall 2-DG treatment in A. malaysianum induced secretion of large amount of commercially viable enzymes compared to other known species.

De Novo Analysis of Wolfiporia cocos Transcriptome to Reveal the Differentially Expressed Carbohydrate-Active Enzymes (CAZymes) Genes During the Early Stage of Sclerotial Growth

The sclerotium of Wolfiporia cocos has been used as an edible mushroom and/or a traditional herbal medicine for centuries. W. cocos sclerotial formation is dependent on parasitism of the wood of Pinus species. Currently, the sclerotial development mechanisms of W. cocos remain largely unknown and the lack of pine resources limit the commercial production. The CAZymes (carbohydrate-active enzymes) play important roles in degradation of the plant cell wall to provide carbohydrates for fungal growth, development, and reproduction. In this study, the transcript profiles from W. cocos mycelium and 2-months-old sclerotium, the early stage of sclerotial growth, were specially analyzed using de novo sequencing technology. A total of 142,428,180 high-quality reads of mycelium and 70,594,319 high-quality reads of 2-months-old sclerotium were obtained. Additionally, differentially expressed genes from the W. cocos mycelium and 2-months-old sclerotium stages were analyzed, resulting in identification of 69 CAZymes genes which were significantly up-regulated during the early stage of sclerotial growth compared to that of in mycelium stage, and more than half of them belonged to glycosyl hydrolases (GHs) family, indicating the importance of W. cocos GHs family for degrading the pine woods. And qRT-PCR was further used to confirm the expression pattern of these up-regulated CAZymes genes. Our results will provide comprehensive CAZymes genes expression information during W. cocos sclerotial growth at the transcriptional level and will lay a foundation for functional genes studies in this fungus. In addition, our study will also facilitate the efficient use of limited pine resources, which is significant for promoting steady development of Chinese W. cocos industry.

Transcript profiling of aquaporins during basidiocarp development in Laccaria bicolor ectomycorrhizal with Picea glauca

Sporocarp formation is part of the reproductive stage in the life cycle of many mycorrhizal macrofungi. Sporocarp formation is accompanied by a transcriptomic switch and profound changes in regulation of the gene families that play crucial roles in the sporocarp initiation and maturation. Since sporocarp growth requires efficient water delivery, in the present study, we investigated changes in transcript abundance of six fungal aquaporin genes that could be cloned from the ectomycorrhizal fungus Laccaria bicolor strain UAMH8232, during the initiation and development of its basidiocarp. Aquaporins are intrinsic membrane proteins facilitating the transmembrane transport of water and other small neutral molecules. In controlled-environment experiments, we induced basidiocarp formation in L. bicolor, which formed ectomycorrhizal associations with white spruce (Picea glauca) seedlings. We profiled transcript abundance corresponding to six fungal aquaporin genes at six different developmental stages of basidiocarp growth and development. We also compared physiological parameters of non-inoculated to mycorrhizal seedlings with and without the presence of basidiocarps. Two L. bicolor aquaporins--JQ585592, a functional channel for CO2, NO and H2O2, and JQ585595, a functional water channel--showed the greatest degree of upregulation during development of the basidiocarp. Our findings point to the importance of aquaporin-mediated transmembrane water and CO2 transport during distinct stages of basidiocarp development.

Recent Developments in Using Advanced Sequencing Technologies for the Genomic Studies of Lignin and Cellulose Degrading Microorganisms

Lignin is a complex polyphenyl aromatic compound which exists in tight associations with cellulose and hemicellulose to form plant primary and secondary cell wall. Lignocellulose is an abundant renewable biomaterial present on the earth. It has gained much attention in the scientific community in recent years because of its potential applications in bio-based industries. Microbial degradation of lignocellulose polymers was well studied in wood decaying fungi. Based on the plant materials they degrade these fungi were classified as white rot, brown rot and soft rot. However, some groups of bacteria belonging to the actinomycetes, α-proteobacteria and β-proteobacteria were also found to be efficient in degrading lignocellulosic biomass but not well understood unlike the fungi. In this review we focus on recent advancements deployed for finding and understanding the lignocellulose degradation by microorganisms. Conventional molecular methods like sequencing 16s rRNA and Inter Transcribed Spacer (ITS) regions were used for identification and classification of microbes. Recent progression in genomics mainly next generation sequencing technologies made the whole genome sequencing of microbes possible in a great ease. The whole genome sequence studies reveals high quality information about genes and canonical pathways involved in the lignin and other cell wall components degradation.

De novo transcriptome analysis of Inonotus baumii by RNA-seq

Inonotus baumii, a basidiomycete white rot fungus, has been widely used as traditional herbal medicine in China, Korea, Japan and other Asian countries for many years. Its extract is of great medicinal importance and plays a valuable role in the immune response and disease resistance. However, limited genetic resources for I. baumii have hindered exploration of this species. In order to gain a molecular understanding of this fungus, Illumina high-throughput technology was used to sequence and analyze the transcriptome of I. baumii, and 280,691 contigs, 43,890 scaffolds and 30,051 unigenes were obtained. Additionally, based on similarity search with known proteins, unigenes were annotated with gene descriptions, gene ontology (GO), clusters of orthologous group (COG), and database of protein families (Pfam) terms. According to the annotation of unigenes, a total of 12 candidate genes involved in the triterpenoid biosynthesis pathway and 21 putative FOLymes (fungal oxidative lignin enzymes) and 176 CAZymes (carbohydrate-active enzymes) were obtained using homology-based BlastX. Moreover, for better understanding of the transcripts function, the BlastX algorithm was used to search for homologous sequences against the Yeast genome. This is the first study on transcriptome analyses of I. baumii, which provided a dataset for functional gene mining and laid a basis for further functional genomics studies of I. baumii.

Genome analysis of medicinal Ganoderma spp. with plant-pathogenic and saprotrophic life-styles

Ganoderma is a fungal genus belonging to the Ganodermataceae family and Polyporales order. Plant-pathogenic species in this genus can cause severe diseases (stem, butt, and root rot) in economically important trees and perennial crops, especially in tropical countries. Ganoderma species are white rot fungi and have ecological importance in the breakdown of woody plants for nutrient mobilization. They possess effective machineries of lignocellulose-decomposing enzymes useful for bioenergy production and bioremediation. In addition, the genus contains many important species that produce pharmacologically active compounds used in health food and medicine. With the rapid adoption of next-generation DNA sequencing technologies, whole genome sequencing and systematic transcriptome analyses become affordable approaches to identify an organism's genes. In the last few years, numerous projects have been initiated to identify the genetic contents of several Ganoderma species, particularly in different strains of Ganoderma lucidum. In November 2013, eleven whole genome sequencing projects for Ganoderma species were registered in international databases, three of which were already completed with genomes being assembled to high quality. In addition to the nuclear genome, two mitochondrial genomes for Ganoderma species have also been reported. Complementing genome analysis, four transcriptome studies on various developmental stages of Ganoderma species have been performed. Information obtained from these studies has laid the foundation for the identification of genes involved in biological pathways that are critical for understanding the biology of Ganoderma, such as the mechanism of pathogenesis, the biosynthesis of active components, life cycle and cellular development, etc. With abundant genetic information becoming available, a few centralized resources have been established to disseminate the knowledge and integrate relevant data to support comparative genomic analyses of Ganoderma species. The current review carries out a detailed comparison of the nuclear genomes, mitochondrial genomes and transcriptomes from several Ganoderma species. Genes involved in biosynthetic pathways such as CYP450 genes and in cellular development such as matA and matB genes are characterized and compared in detail, as examples to demonstrate the usefulness of comparative genomic analyses for the identification of critical genes. Resources needed for future data integration and exploitation are also discussed.

MicroRNA-like small RNAs prediction in the development of Antrodia cinnamomea

Antrodia cinnamomea, a precious, host-specific brown-rot fungus that has been used as a folk medicine in Taiwan for centuries is known to have diverse bioactive compounds with potent pharmaceutical activity. In this study, different fermentation states of A. cinnamomea (wild-type fruiting bodies and liquid cultured mycelium) were sequenced using the next-generation sequencing (NGS) technique. A 45.58 Mb genome encoding 6,522 predicted genes was obtained. High quality reads were assembled into a total of 13,109 unigenes. Using a previously constructed pipeline to search for microRNAs (miRNAs), we then identified 4 predicted conserved miRNA and 63 novel predicted miRNA-like small RNA (milRNA) candidates. Target prediction revealed several interesting proteins involved in tri-terpenoid synthesis, mating type recognition, chemical or physical sensory protein and transporters predicted to be regulated by the miRNAs and milRNAs.

Transcriptome analysis and its application in identifying genes associated with fruiting body development in basidiomycete Hypsizygus marmoreus

To elucidate the mechanisms of fruit body development in H. marmoreus, a total of 43609521 high-quality RNA-seq reads were obtained from four developmental stages, including the mycelial knot (H-M), mycelial pigmentation (H-V), primordium (H-P) and fruiting body (H-F) stages. These reads were assembled to obtain 40568 unigenes with an average length of 1074 bp. A total of 26800 (66.06%) unigenes were annotated and analyzed with the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Eukaryotic Orthologous Group (KOG) databases. Differentially expressed genes (DEGs) from the four transcriptomes were analyzed. The KEGG enrichment analysis revealed that the mycelium pigmentation stage was associated with the MAPK, cAMP, and blue light signal transduction pathways. In addition, expression of the two-component system members changed with the transition from H-M to H-V, suggesting that light affected the expression of genes related to fruit body initiation in H. marmoreus. During the transition from H-V to H-P, stress signals associated with MAPK, cAMP and ROS signals might be the most important inducers. Our data suggested that nitrogen starvation might be one of the most important factors in promoting fruit body maturation, and nitrogen metabolism and mTOR signaling pathway were associated with this process. In addition, 30 genes of interest were analyzed by quantitative real-time PCR to verify their expression profiles at the four developmental stages. This study advances our understanding of the molecular mechanism of fruiting body development in H. marmoreus by identifying a wealth of new genes that may play important roles in mushroom morphogenesis.

Proteome exploration to provide a resource for the investigation of Ganoderma lucidum

Ganoderma lucidum is a basidiomycete white rot fungus that has been used for medicinal purposes worldwide. Although information concerning its genome and transcriptome has recently been reported, relatively little information is available for G. lucidum at the proteomic level. In this study, protein fractions from G. lucidum at three developmental stages (16-day mycelia, and fruiting bodies at 60 and 90 days) were prepared and subjected to LC-MS/MS analysis. A search against the G. lucidum genome database identified 803 proteins. Among these proteins, 61 lignocellulose degrading proteins were detected, most of which (49 proteins) were found in the 90-day fruiting bodies. Fourteen TCA-cycle related proteins, 17 peptidases, two argonaute-like proteins, and two immunomodulatory proteins were also detected. A majority (470) of the 803 proteins had GO annotations and were classified into 36 GO terms, with "binding", "catalytic activity", and "hydrolase activity" having high percentages. Additionally, 357 out of the 803 proteins were assigned to at least one COG functional category and grouped into 22 COG classifications. Based on the results from the proteomic and sequence alignment analyses, a potentially new immunomodulatory protein (GL18769) was expressed and shown to have high immunomodulatory activity. In this study, proteomic and biochemical analyses of G. lucidum were performed for the first time, revealing that proteins from this fungus can play significant bioactive roles and providing a new foundation for the further functional investigations that this fungus merits.

Traversing the fungal terpenome

Fungi (Ascomycota and Basidiomycota) are prolific producers of structurally diverse terpenoid compounds. Classes of terpenoids identified in fungi include the sesqui-, di- and triterpenoids. Biosynthetic pathways and enzymes to terpenoids from each of these classes have been described. These typically involve the scaffold generating terpene synthases and cyclases, and scaffold tailoring enzymes such as e.g. cytochrome P450 monoxygenases, NAD(P)+ and flavin dependent oxidoreductases, and various group transferases that generate the final bioactive structures. The biosynthesis of several sesquiterpenoid mycotoxins and bioactive diterpenoids has been well-studied in Ascomycota (e.g. filamentous fungi). Little is known about the terpenoid biosynthetic pathways in Basidiomycota (e.g. mushroom forming fungi), although they produce a huge diversity of terpenoid natural products. Specifically, many trans-humulyl cation derived sesquiterpenoid natural products with potent bioactivities have been isolated. Biosynthetic gene clusters responsible for the production of trans-humulyl cation derived protoilludanes, and other sesquiterpenoids, can be rapidly identified by genome sequencing and bioinformatic methods. Genome mining combined with heterologous biosynthetic pathway refactoring has the potential to facilitate discovery and production of pharmaceutically relevant fungal terpenoids.

Comparative transcriptomics of the model mushroom Coprinopsis cinerea reveals tissue-specific armories and a conserved circuitry for sexual development

Background

It is well known that mushrooms produce defense proteins and secondary metabolites against predators and competitors; however, less is known about the correlation between the tissue-specific expression and the target organism (antagonist) specificity of these molecules. In addition, conserved transcriptional circuitries involved in developing sexual organs in fungi are not characterized, despite the growing number of gene expression datasets available from reproductive and vegetative tissue. The aims of this study were: first, to evaluate the tissue specificity of defense gene expression in the model mushroom Coprinopsis cinerea and, second, to assess the degree of conservation in transcriptional regulation during sexual development in basidiomycetes.

Results

In order to characterize the regulation in the expression of defense loci and the transcriptional circuitries controlling sexual reproduction in basidiomycetes, we sequenced the poly (A)-positive transcriptome of stage 1 primordia and vegetative mycelium of C. cinerea A43mutB43mut. Our data show that many genes encoding predicted and already characterized defense proteins are differentially expressed in these tissues. The predicted specificity of these proteins with regard to target organisms suggests that their expression pattern correlates with the type of antagonists these tissues are confronted with. Accordingly, we show that the stage 1 primordium-specific protein CC1G_11805 is toxic to insects and nematodes. Comparison of our data to analogous data from Laccaria bicolor and Schizophyllum commune revealed that the transcriptional regulation of nearly 70 loci is conserved and probably subjected to stabilizing selection. A Velvet domain-containing protein was found to be up-regulated in all three fungi, providing preliminary evidence of a possible role of the Velvet protein family in sexual development of basidiomycetes. The PBS-soluble proteome of C. cinerea primordia and mycelium was analyzed by shotgun LC-MS. This proteome data confirmed the presence of intracellular defense proteins in primordia.

Conclusions

This study shows that the exposure of different tissues in fungi to different types of antagonists shapes the expression pattern of defense loci in a tissue-specific manner. Furthermore, we identify a transcriptional circuitry conserved among basidiomycetes during fruiting body formation that involves, amongst other transcription factors, the up-regulation of a Velvet domain-containing protein.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-492) contains supplementary material, which is available to authorized users.

De novo assembly of Auricularia polytricha transcriptome using Illumina sequencing for gene discovery and SSR marker identification

Auricularia polytricha (Mont.) Sacc., a type of edible black-brown mushroom with a gelatinous and modality-specific fruiting body, is in high demand in Asia due to its nutritional and medicinal properties. Illumina Solexa sequenceing technology was used to generate very large transcript sequences from the mycelium and the mature fruiting body of A. polytricha for gene discovery and molecular marker development. De novo assembly generated 36,483 ESTs with an N50 length of 636 bp. A total of 28,108 ESTs demonstrated significant hits with known proteins in the nr database, and 94.03% of the annotated ESTs showed the greatest similarity to A. delicata, a related species of A. polytricha. Functional categorization of the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways revealed the conservation of genes involved in various biological processes in A. polytricha. Gene expression profile analysis indicated that a total of 2,057 ESTs were differentially expressed, including 1,020 ESTs that were up-regulated in the mycelium and 1,037 up-regulated in the fruiting body. Functional enrichment showed that the ESTs associated with biosynthesis, metabolism and assembly of proteins were more active in fruiting body development. The expression patterns of homologous transcription factors indicated that the molecular mechanisms of fruiting body formation and development were not exactly the same as for other agarics. Interestingly, an EST encoding tyrosinase was significantly up-regulated in the fruiting body, indicating that melanins accumulated during the processes of the formation of the black-brown color of the fruiting body in A. polytricha development. In addition, a total of 1,715 potential SSRs were detected in this transcriptome. The transcriptome analysis of A. polytricha provides valuable sequence resources and numerous molecular markers to facilitate further functional genomics studies and genetic researches on this fungus.

Abundant and selective RNA-editing events in the medicinal mushroom Ganoderma lucidum

RNA editing is a widespread, post-transcriptional molecular phenomenon that diversifies hereditary information across various organisms. However, little is known about genome-scale RNA editing in fungi. In this study, we screened for fungal RNA editing sites at the genomic level in Ganoderma lucidum, a valuable medicinal fungus. On the basis of our pipeline that predicted the editing sites from genomic and transcriptomic data, a total of 8906 possible RNA-editing sites were identified within the G. lucidum genome, including the exon and intron sequences and the 5'-/3'-untranslated regions of 2991 genes and the intergenic regions. The major editing types included C-to-U, A-to-G, G-to-A, and U-to-C conversions. Four putative RNA-editing enzymes were identified, including three adenosine deaminases acting on transfer RNA and a deoxycytidylate deaminase. The genes containing RNA-editing sites were functionally classified by the Kyoto Encyclopedia of Genes and Genomes enrichment and gene ontology analysis. The key functional groupings enriched for RNA-editing sites included laccase genes involved in lignin degradation, key enzymes involved in triterpenoid biosynthesis, and transcription factors. A total of 97 putative editing sites were randomly selected and validated by using PCR and Sanger sequencing. We presented an accurate and large-scale identification of RNA-editing events in G. lucidum, providing global and quantitative cataloging of RNA editing in the fungal genome. This study will shed light on the role of transcriptional plasticity in the growth and development of G. lucidum, as well as its adaptation to the environment and the regulation of valuable secondary metabolite pathways.

Biosynthesis of terpenoid natural products in fungi

Tens of thousands of terpenoid natural products have been isolated from plants and microbial sources. Higher fungi (Ascomycota and Basidiomycota) are known to produce an array of well-known terpenoid natural products, including mycotoxins, antibiotics, antitumor compounds, and phytohormones. Except for a few well-studied fungal biosynthetic pathways, the majority of genes and biosynthetic pathways responsible for the biosynthesis of a small number of these secondary metabolites have only been discovered and characterized in the past 5-10 years. This chapter provides a comprehensive overview of the current knowledge on fungal terpenoid biosynthesis from biochemical, genetic, and genomic viewpoints. Enzymes involved in synthesizing, transferring, and cyclizing the prenyl chains that form the hydrocarbon scaffolds of fungal terpenoid natural products are systematically discussed. Genomic information and functional evidence suggest differences between the terpenome of the two major fungal phyla--the Ascomycota and Basidiomycota--which will be illustrated for each group of terpenoid natural products.

A first insight into Pycnoporus sanguineus BAFC 2126 transcriptome

Fungi of the genus Pycnoporus are white-rot basidiomycetes widely studied because of their ability to synthesize high added-value compounds and enzymes of industrial interest. Here we report the sequencing, assembly and analysis of the transcriptome of Pycnoporus sanguineus BAFC 2126 grown at stationary phase, in media supplemented with copper sulfate. Using the 454 pyrosequencing platform we obtained a total of 226,336 reads (88,779,843 bases) that were filtered and de novo assembled to generate a reference transcriptome of 7,303 transcripts. Putative functions were assigned for 4,732 transcripts by searching similarities of six-frame translated sequences against a customized protein database and by the presence of conserved protein domains. Through the analysis of translated sequences we identified transcripts encoding 178 putative carbohydrate active enzymes, including representatives of 15 families with roles in lignocellulose degradation. Furthermore, we found many transcripts encoding enzymes related to lignin hydrolysis and modification, including laccases and peroxidases, as well as GMC oxidoreductases, copper radical oxidases and other enzymes involved in the generation of extracellular hydrogen peroxide and iron homeostasis. Finally, we identified the transcripts encoding all of the enzymes involved in terpenoid backbone biosynthesis pathway, various terpene synthases related to the biosynthesis of sesquiterpenoids and triterpenoids precursors, and also cytochrome P450 monooxygenases, glutathione S-transferases and epoxide hydrolases with potential functions in the biodegradation of xenobiotics and the enantioselective biosynthesis of biologically active drugs. To our knowledge this is the first report of a transcriptome of genus Pycnoporus and a resource for future molecular studies in P. sanguineus.

Illumina-based de novo transcriptome sequencing and analysis of Amanita exitialis basidiocarps

Amanita exitialis is a lethal mushroom that was first discovered in Guangdong Province, China. The high content of amanitin in its basidiocarps makes it lethal to humans. To comprehensively characterize the A. exitialis transcriptome and analyze the Amanita toxins as well as their related gene family, transcriptome sequencing of A. exitialis was performed using Illumina HiSeq 2000 technology. A total of 25,563,688 clean reads were collected and assembled into 62,137 cDNA contigs with an average length of 481 bp and N50 length of 788 bp. A total of 27,826 proteins and 39,661 unigenes were identified among the assembled contigs. All of the unigenes were classified into 166 functional categories for understanding the gene functions and regulation pathways. The genes contributing to toxic peptide biosynthesis were analyzed. From this set, eleven gene sequences encoding the toxins or related cyclic peptides were discovered in the transcriptome. Three of these sequences matched the peptide toxins α-amanitin, β-amanitin, and phallacidin, while others matched amanexitide and seven matched unknown peptides. All of the genes encoding peptide toxins were confirmed by polymerase chain reaction (PCR) in A. exitialis, and the phylogenetic relationships among these proprotein sequences were discussed. The gene polymorphism and degeneracy of the toxin encoding sequences were found and analyzed. This study provides the first primary transcriptome of A. exitialis, which provided comprehensive gene expression information on the lethal amanitas at the transcriptional level, and could lay a strong foundation for functional genomics studies in those fungi.

Antidiabetic, antihyperlipidemic and antioxidant activities of a novel proteoglycan from ganoderma lucidum fruiting bodies on db/db mice and the possible mechanism

Previously, we screened a proteoglycan for anti-hyperglycemic, named FYGL, from Ganoderma Lucidum. For further research of the antidiabetic mechanisms of FYGL in vivo, the glucose homeostasis, activities of insulin-sensitive enzymes, glucose transporter expression and pancreatic function were analyzed using db/db mice as diabetic models in the present work. FYGL not only lead to a reduction in glycated hemoglobin level, but also an increase in insulin and C-peptide level, whereas a decrease in glucagons level and showed a potential for the remediation of pancreatic islets. FYGL also increased the glucokinase activities, and simultaneously lowered the phosphoenol pyruvate carboxykinase activities, accompanied by a reduction in the expression of hepatic glucose transporter protein 2, while the expression of adipose and skeletal glucose transporter protein 4 was increased. Moreover, the antioxidant enzyme activities were also increased by FYGL treatment. Thus, FYGL was an effective antidiabetic agent by enhancing insulin secretion and decreasing hepatic glucose output along with increase of adipose and skeletal muscle glucose disposal in the late stage of diabetes. Furthermore, FYGL is beneficial against oxidative stress, thereby being helpful in preventing the diabetic complications.

Profiling and quantifying differential gene transcription provide insights into ganoderic acid biosynthesis in Ganoderma lucidum in response to methyl jasmonate

Ganoderma lucidum is a mushroom with traditional medicinal properties that has been widely used in China and other countries in Eastern Asia. Ganoderic acids (GA) produced by G. lucidum exhibit important pharmacological activities. Previous studies have demonstrated that methyl jasmonate (MeJA) is a potent inducer of GA biosynthesis and the expression of genes involved in the GA biosynthesis pathway in G. lucidum. To further explore the mechanism of GA biosynthesis, cDNA-Amplified Fragment Length Polymorphism (cDNA-AFLP) was used to identify genes that are differentially expressed in response to MeJA. Using 64 primer combinations, over 3910 transcriptionally derived fragments (TDFs) were obtained. Reliable sequence data were obtained for 390 of 458 selected TDFs. Ninety of these TDFs were annotated with known functions through BLASTX searching the GenBank database, and 12 annotated TDFs were assigned into secondary metabolic pathways by searching the KEGGPATHWAY database. Twenty-five TDFs were selected for qRT-PCR analysis to confirm the expression patterns observed with cDNA-AFLP. The qRT-PCR results were consistent with the altered patterns of gene expression revealed by the cDNA-AFLP technique. Additionally, the transcript levels of 10 genes were measured at the mycelium, primordia, and fruiting body developmental stages of G. lucidum. The greatest expression levels were reached during primordia for all of the genes except cytochrome b2 reached its highest expression level in the mycelium stage. This study not only identifies new candidate genes involved in the regulation of GA biosynthesis but also provides further insight into MeJA-induced gene expression and secondary metabolic response in G. lucidum.