Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 03 February 2016
Sec. Food Microbiology

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

\r\nShaopeng Zhang&#x;Shaopeng Zhang1†Bingxiong Hu&#x;Bingxiong Hu1†Wei WeiWei Wei2Ying XiongYing Xiong3Wenjun Zhu*Wenjun Zhu1*Fang PengFang Peng1Yang YuYang Yu4Yonglian ZhengYonglian Zheng1Ping Chen*Ping Chen1*
  • 1College of Biology and Pharmaceutical Engineering, Wuhan Polytechnic University, Wuhan, China
  • 2Institute for Interdisciplinary Research, Jianghan University, Wuhan, China
  • 3Hefei Enzyme Information Technology Co., Ltd, Wuhan, China
  • 4College of Plant Protection, Southwest University, Chongqing, China

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.

Introduction

Wolfiporia cocos (Schwein.) Ryvarden et Gilb. (Basidiomycota, Polyporaceae) is a fungus that parasitizes the roots of diverse species of Pinus. The fungus has a wide distribution in East Asia, including China, Japan, and Korea, and other regions of the world (Dai et al., 2009; Esteban, 2009; Wang et al., 2013). The sclerotium of W. cocos, also known as Fuling in China, has been used as an edible mushroom and/or a traditional herbal medicine for centuries (Dai et al., 2009; Kobira et al., 2012). Phytochemical and pharmacological research demonstrates that the major pharmacologically active components from Fuling, such as polysaccharides and triterpenes, possess a wide spectrum of pharmacological activities, including anti-viral., anti-tumor, anti-oxidant, free-radical scavenging, anti-rejection, anti-hyperglycemic, antibacterial., anti-inflammatory, and anti-hypertonic stress effects (Huang et al., 2007; Esteban, 2009; Rios, 2011; Feng et al., 2013; Wang et al., 2013; Zhao et al., 2013a; Tang et al., 2014).

Wolfiporia cocos sclerotial formation is dependent on parasitism of the wood of Pinus species (Kubo et al., 2006). However, commercial production of W. cocos sclerotia is currently limited by severe habitat destruction, ineffective protection, and the lack of Pinus species materials (Wang et al., 2012). Thus, exploration of novel artificial cultivation techniques that are independent of pine resources and the sustainable production of W. cocos sclerotia have attracted much interest. Studying the parasitic mechanisms of W. cocos on Pinus species and the genetic basis of sclerotial development will improve our understanding of the overall biology of the fungus and may facilitate commercial production.

With the development of next-generation sequencing technology, RNA deep-sequencing, especially the de novo sequencing, has been widely used for transcriptomic and functional genes study in several mushrooms (Liu et al., 2012; Yang et al., 2012; Shu et al., 2013). Nevertheless, compared with other fungi, the mechanisms of W. cocos sclerotial formation and biosynthesis of polysaccharides and triterpenes remain largely unknown, especially the associated transcriptomic data at early stage of sclerotial growth. Moreover, in other fungi, many studies have demonstrated that a series of important genes involved in sclerotial development are upregulated at early stage of sclerotial development (Yu et al., 2012a; Xing et al., 2013; Zhu et al., 2013; Xiao et al., 2014). Thus, sequencing and studying the W. cocos transcriptome of early stage of sclerotial growth will authentically provide great help for identifying more significant genes regulating sclerotial formation and other important economic traits.

To infect and colonize the plants tissues successfully, both saprophytic and phytopathogenic fungi secrete a diverse range of carbohydrate-active enzymes (CAZymes), especially cell wall degrading enzymes (CWDEs), to degrade barrier of the plant cell wall, which are known to be mainly composed of cellulose, hemicellulose, and pectin (Cosgrove, 2005). Based on catalytic activities associated motifs and/or domains, the CAZymes are functionally classified into six classes: glycosyl hydrolases (GHs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), carbohydrate-binding modules (CBMs), glycosyltransferases (GTs), and auxiliary activities (AAs; Christian et al., 2014; Blackman et al., 2015).

The carbohydrates from degraded plant tissue, such as mono and oligosaccharides, can be utilized as nutrition for fungal growth, development, and reproduction. And a number of studies have also suggested that the nutrition conditions, such as optimal carbon source, play important roles on the sclerotial formation of different fungi (Rai and Agnihotri, 1971; Xing et al., 2011). Besides, when cultured on artificial medium, such as PDA, W. cocos could directly absorb and utilize small molecular nutrition source, such as glucose for mycelium growth, but it was unable to form sclerotia (Kubo et al., 2006). And the sclerotial formation and growth relied on colonization on the wood tissue of Pinus species to sustainably absorb and utilize the carbohydrates and other unknown nutrition source degraded from pine woods by the catalytic activities of CAZymes. In addition, several studies have also demonstrated that when fungi colonized on different host species or tissue components, hydrolytic enzymes showed activities preferences (King et al., 2011; Couturier et al., 2012). Thus, secretion of various CAZymes is a significant means by which W. cocos colonize and develop sclerotia successfully on the wood of Pinus species.

In this study, the transcript profiles from W. cocos mycelium and 2-months-old sclerotium, the early stage of sclerotial growth, were firstly analyzed using de novo sequencing technology, resulting in the identification of many differentially expressed genes (DEGs) during early stage of sclerotial formation, especially the CAZymes. And qRT-PCR was also used to confirm the expression profiles of these CAZymes genes. Our results will provide significant genes resource and facilitate future genes functional studies to gain a better understanding of sclerotial formation and growth mechanisms and other development processes of important economic traits. And data obtained in this study also provide a starting point to achieve a better understanding of the interaction system of W. cocosPinus species plants for efficient use of Pinus species wood and the steady development of Chinese W. cocos industry.

Materials and Methods

Strains and Culture Conditions

The W. cocos strain AS5.78 was obtained from the Agricultural Culture Collection of China, Institute of Soil and Fertilization, Chinese Academy of Agricultural Sciences, Beijing. W. cocos mycelia were grown on a cellophane membrane laid on potato dextrose agar (PDA) medium at 25°C for 7 days. Sclerotium of W. cocos strain AS5.78 at 2-months-old (Figures 1A,B) was obtained from the field in Luotian county, Hubei province, China. The mycelia and sclerotium of W. cocos were collected and immediately frozen in liquid nitrogen, then stored at –80°C for total RNA extraction.

FIGURE 1
www.frontiersin.org

FIGURE 1. The mycelium and sclerotia of W. cocos used for De Novo Sequencing. (A) Colony morphology W. cocos mycelium grown on PDA for 7 days at 25°C. (B) W. cocos sclerotium of 2-months-old, scale bar = 5 cm. (C) The total RNA isolated from mycelium and 2-months-old sclerotium.

RNA Extraction, cDNA Library Construction, and Sequencing

Total RNA of mycelia and 2-months-old sclerotium were isolated with TriZOL reagent (Invitrogen, Carlsbad, CA, USA) according to manufacturer’s instructions (Figure 1C). The RNA quantity and quality were exampled with UV spectrophotometer (DU800, Beckman Coulter, USA) and gel electrophoresis, respectively, and then treated with DNase I (RNase Free; Takara, Dalian, China) to remove residual DNA according to manufacturer’s protocols. The construction of cDNA library was performed as described (Shu et al., 2013). Finally, the cDNA library was sequenced with Illumina HiSeq 2000 (Illumina Inc, San Diego CA, USA).

Hybrid Assembly and Unigenes Function Annotation

We collected another RNA-seq data set from NCBI under accession number SRP018935 (Shu et al., 2013), which including two sample from mycelium and sclerotium. This dataset was included in all the following analysis. Before assembly, the raw reads contain adapters, and unknown (N) sequences were discarded. To obtain the high-quality clean reads, all reads were trimmed at the base with quality score lower than 20, and it will be discard if it is shorter than 50 bp. The individual reads will be kept if one of a pair reads is discarded. Transcript assembly was conducted in two steps. First, the reference genome of W. cocos was downloaded from JGI database1 and all clean reads were mapped on it using tophat2 (version 2.0.13) with default parameters, except that unmapped reads were output. Sequentially, cufflinks (v2.2.1) was used to assembly transcripts based on the map results. Second, all the unmapped reads from two datasets were combined and further assembly was conducted with trinity (version r20140413p1), and all the transcripts generated in this step were cluster with cd-hit-est (version 4.6), and both identity of 99% and coverage of 90% for shorter sequences were required and the longest sequence for each cluster were kept as the final unigenes together with transcripts generated by cufflinks.

Homologous sequences searching and functional annotation of all unigenes were performed by the using of BLASTx program against the NCBI NR database2, Swissprot3, and KEGG pathway4 (E-value <10-5).

Expression Difference Analysis

As previously described (Yu et al., 2012b; Shu et al., 2013), the RPKM method (Reads per kb per Million reads) was used to calculate the unigene expression. The statistical method false discovery rate (FDR) was used in multiple tests to correct the threshold of P-value for the identification of DEGs between different developmental stages. In our analysis, we chose those with ratio ≥ 2 and FDR ≤ 0.001. The DEGs was analyzed with DEGSeq (Wang et al., 2010).

CAZymes Annotation

The annotated unigenes from W. cocos transcriptome that are related to carbohydrate-active enzymes were further confirmed by the using of BLASTx program against CAZymes database5 at the threshold value of 1.0E-20. Unigenes possessing a sequence identity more than 60% with biochemically characterized CAZymes were considered as candidate.

qRT-PCR Confirmation of CAZymes Genes

Total RNA of mycelia and sclerotia were extracted and treated with DNase I as described above. And the RevertAidTM First Strand cDNA Synthesis Kit (MBI Fermentas, Lithuania) was used to generate the first strand cDNA according to the manufacturer’s protocols. Gene expression was analyzed by qRT-PCR using a Bio-Rad CFX96 Real Time System (Bio-Rad, USA) and QuantiTect SYBR Green PCR master mix (Bio-Rad, USA), according to the manufacturer’s instructions. The PCR conditions were as follows: denaturation at 94°C for 3 min; 40 rounds of 94°C for 20 s, 55°C for 20 s, and 72°C for 30 s; final step of 72°C for 10 min. The primers for qRT-PCR are listed in Supplementary Table S1. The expression of W. cocos alpha-tubulin gene (F:5′ACTCCAGCTTGGACTTCTTG3′ and R: 5′TCTTCGTCTTCCACTCCTTTG3′; Shu et al., 2013) was used to normalize the RNA sample for each qRT-PCR. For each gene, qRT-PCR assays were repeated at least three times, with each repetition having three technical replicates.

Results

Illumina Sequencing and De Novo Assembly

To obtain an overview of the W. cocos genes expression profiles during early growth stages, cDNA samples from mycelium and 2-months-old sclerotium were sequenced by using an Illumina HiSeq2000 platform, and the comparative analysis between mycelium and 2-months-old sclerotium was studied. After filtering the adapters, low quality and unknown (N) sequences, 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.

Then, by the using of Cufflinks and Trinity software, these high-quality reads were assembled together to totally generate 62,143 unigenes with length above 200 bp. For the all-unigenes set, the average length of the was 1114 bp, the max and min length were 16,339 bp and 201 bp, respectively. About 43.54% (27,059) of all-unigenes were longer than 1,000 bp (Figure 2). In summary, 57,686 unigenes from mycelium, 56,010 unigenes from 2-months-old sclerotium in our data set, and 42,357 unigenes from mycelium, 43,091 unigenes from matured sclerotium in data set of Shu et al. (2013) were obtained in (Figure 3), respectively.

FIGURE 2
www.frontiersin.org

FIGURE 2. Assembled unigenes length distribution of W. cocos transcriptome from mycelium and 2-months-old sclerotium.

FIGURE 3
www.frontiersin.org

FIGURE 3. The number of unigenes for mycelium and sclerotium from different two data sets. Mycelium_P0 and Sclerotia_P1 were from the data set sequence in this project and Mycelium_W0, Sclerotia_W1 were from Shu et al. (2013).

Functional Annotation of W. cocos Transcriptome

BLASTX alignments (E-value cut-off of 10-5) of all unigenes sequences against several protein databases, including NBCI non-redundant (NR), Swiss-Prot, GO, and KEGG, showed that, respectively 10,137, 4,983, 8,024, and 7,152 unigenes could match significantly to proteins with known biological functions (Table 1), whereas the remainder are currently lacking in the database. Most of them were short fragments, and some of them might be non-coding RNA sequences or new genes.

TABLE 1
www.frontiersin.org

TABLE 1. All-in-one list of annotations.

The species distribution of the top BLASTX match against the NR database for the W. cocos unigenes showed that about 39.32% of W. cocos unigenes match with genes of Fibroporia, followed by Fomitopsis pinicola (23.50%), Postia placenta (10.15%), Ceriporiopsis subvermispora (7.82%), Dichomitus squalens (2.94%), Trametes versicolor (2.79%), and Phanerochaete carnosa (1.34%). The rest of 12.15% W. cocos unigenes had first hits with other fungal species (Figure 4).

FIGURE 4
www.frontiersin.org

FIGURE 4. Distribution of the homology search of expressed sequence tags against the nr database.

Analysis of Differentially Expressed Genes During Sclerotial Formation

For identification of DEGs between W. cocos mycelium and early sclerotium growth stages, the expression data of unigenes were analyzed with DEGSeq. Compared with mycelium, 1,371 and 1,255 unigenes showed upregulated and downregulated, respectively, in 2-months-old sclerotium (Figure 5A), 3,052 and 3,208 unigenes showed upregulated and downregulated, respectively, in matured sclerotium were identified in the data set of Shu et al. (2013), (Figure 5B). One thousand and ninety-three unigenes were differentially expressed in both data sets (Figure 5C).

FIGURE 5
www.frontiersin.org

FIGURE 5. Gene transcription profile between mycelium and sclerotial stages. Scatter plot of total unigenes from the W. cocos transcriptome. The data was normalized as RPKM values and represented on a log10 scale. Red areas represent up-regulated unigenes in sclerotia, green areas represent down-regulated unigenes in the sclerotia and the blue areas represent no significant expression difference unigenes between mycelium and sclerotia. (A) Sclerotia_P1 and Mycelium_P0 were from the data set sequence in this project. (B) Sclerotia_W1 and Mycelium_W0 were from Shu et al. (2013). (C) Differentially expressed unigenes in both data sets.

CAZymes Genes In W. cocos

Carbohydrate-active enzymes, including Glycoside hydrolases (GHs), CEs, AAs, CBMs, GTs, and PLs, are involved in carbohydrate metabolism. In our study, the BLASTx results indicated that a total of 306 unigenes from W. cocos transcriptome was identified to CAZymes. Among of these unigenes, 181 homologs belonged to the GHs, 88 candidates to GTs, 28 candidates to CEs, four candidates to CBMs, three to AAs, and two to PLs (Supplementary Table S2). In the GHs, GH16 and GH5 family both contained 28 unigenes, followed by GH18 (16 unigenes), GH28 (13 unigenes), GH31 (11 unigenes), GH3 and GH13 (both nine unigenes). In the GTs class, GT8 and GT2 both possessed 12 unigenes and was the most dominant, followed by GT4 (seven unigenes), GT20 (six unigenes), and GT15 (five unigenes). In CEs, CE4, and CE16 both contained seven unigenes, followed by CE10 (six unigenes). In CBMs, three unigenes was CMB21 and one was CMB43. There were only two PL14 in PLs class and three AA6 in AAs class (Supplementary Table S2).

The differential expression of these CAZymes genes was evaluated. The results indicated that 69 CAZymes genes were significantly up-regulated (≥1.5-fold) in the 2-months-old sclerotium compared to that of in mycelium (Supplementary Table S2, yellow fluorescence label), containing 43 GHs, 17 GTs, five CEs, two PLs, and two AAs. Then, some of these up-regulated genes were selected to further examined the expression pattern by the using of qRT-PCR (Figure 6) and the results correspond with de novo transcriptome sequencing data. The results suggested that these up-regulated CAZymes unigenes is involved in W. cocos sclerotial formation possibly when colonize on pine wood at the early stage.

FIGURE 6
www.frontiersin.org

FIGURE 6. qRT-PCR validation of differentially expressed CAZymes genes. The relative expression of target genes in mycelium stage was set as level 1. Expression of W. cocos alpha-tubulin gene was used to normalize different samples. Bars represent means and standard deviations (three replications).

Discussion

Sclerotia of W. cocos possess important medicinal value and have been widely used in East Asia for centuries (Dai et al., 2009; Kobira et al., 2012). The sclerotial formation is dependent on parasitism of the wood of Pinus species (Kubo et al., 2006). At present, the lack of Pinus species materials currently limits the commercial production of W. cocos sclerotia. Thus, the efficient use of limited pine resources is significant for promoting steady development of Chinese W. cocos industry.

While the pharmacologically activities of polysaccharides and triterpenes have been extensively studied (Esteban, 2009; Wang et al., 2015a,b), the genomic and transcriptomic information remain largely unknown on this fungus, especially for the CAZymes genes which play important roles in carbohydrate degradation and plant tissue colonization (Kubicek et al., 2014). In this study, the transcriptome of W. cocos from mycelium and 2-months-old sclerotium were analyzed and numerous DEGs, especially the CAZymes genes relating to carbohydrates degradation, were identified. Comparative analysis with the data of Shu et al. (2013) revealed that W. cocos highly expressed CAZymes genes, especially in the early stage of sclerotial development, which supports with the facts that W. cocos needs to sustainably absorb and utilize the carbohydrates and other unknown nutrition source degraded from pine woods for sclerotial growth.

Among these CAZymes genes, the GHs family was the most dominant prevalent. These results are similar to those previously described by Yu et al. (2012b) in Ganoderma lucidum (Basidiomycota, Polyporaceae). GHs mainly hydrolyze the glycoside bond between two or more carbohydrates or between a carbohydrate and a non-carbohydrate (Cantarel et al., 2009). To date, based on amino acid sequence, GHs are classified into 127 families and 91 of them were found in fungi (Zhao et al., 2013b). According to our transcriptome data, families GH16, GH18, GH28, and GH5 were the dominant among these up-regulated CAZymes genes in the early stage of sclerotial development (Supplementary Table S2). GH16 family enzymes possess xyloglucanase activity and may degrade β-1,3-glucans or xyloglucans, which play significant roles in hemicellulose degrading (Baumann et al., 2007; Nakajima et al., 2012). GH18 may degrade chitin and cleave the β-1,4-linkage between N-acetylglucosamine residues at the base of the oligosaccharide chain while GH28 possess polygalacturonases activity and play a critical role in pectin degradation (Zhao et al., 2013b; Blackman et al., 2015). GH5 is one of the largest GH families and consists of a wide range of enzymes activity on different substrates, such as β-1,3-glucans, β-1,4-glucans in cellulose, and β-1,4-mannans in hemicellulose (Aspeborg et al., 2012; Blackman et al., 2014). In our study, six GH5, 10 GH16, five GH18, and six GH28 genes showed up-regulated expression in the early stage of sclerotial growth (2 months; Supplementary Table S2), suggesting that these enzymes play important roles in plant cell wall degradation to provide sufficient nutrition for the growth of W. cocos sclerotium.

In addition, as W. cocos is unable to form sclerotium in the absence of pine woods, thus we inferred that some special components from Pinus species plants could induce the formation and development of W. cocos sclerotium, and these unknown components will also definitely induce the differential expression of genes involved in W. cocos sclerotial development. Characterization of these potential candidates of sclerotial development associated genes will be helpful to elucidate the mechanism of sclerotial formation. And studying on the DEGs involved in chemical components synthesis and metabolism will also give a new clue to identify the components involved in inducing W. cocos sclerotial formation.

In summary, we have analyzed the W. cocos transcriptome in mycelium and early stages of sclerotial growth, resulting in identifying various DEGs, especially the CAZymes genes. Characterization of these genes will enhance our understanding on mechanism of W. cocos sclerotial development. However, our finding also arise more questions to be answered, such as, what kinds of components promote W. cocos sclerotial formation? Do these components only come from pine woods or from the interaction product of W. cocos-host? And how do these components induce W. cocos sclerotial formation?

Conclusion

This study contributes to a better understanding of the DEGs during the early stage of W. cocos sclerotial growth in comparison with that of mycelium stage, especially the CAZymes genes. And it also provides a valuable first step toward the understanding and analysis of the mechanism of W. cocos sclerotial development and interaction system of W. cocos-pine woods.

Author Contributions

Conceived and designed the experiments: WZ. Performed the experiments: SZ, BH, WW, YX, and WZ. Analyzed the transcription data: SZ, BH, WW, YX, WZ, FP, and YY. Contributed reagents/materials/analysis tools: WW, WZ, PC, and YZ. Wrote the paper: WW, WZ, PC, and YZ. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This research was supported by State of Traditional Chinese Medicine Industry Research Program (No. 201107009), Major Support Program granted from Wuhan Polytechnic University (No. 2013472) and National Natural Science Foundation of China (No. 31501587). We also thank the reviewers for their thoughtful suggestions.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.00083

TABLE S1 | Primers for qRT-PCR to verify W. cocos CAZymes genes.

TABLE S2 | Differentially expressed carbohydrate-active enzymes (CAZymes) genes between W. cocos mycelium and 2-months-old sclerotium. The up-regulated genes were labeled with yellow fluorescence label.

Footnotes

  1. ^ http://genome.jgi.doe.gov/Wolco1/Wolco1.home.html
  2. ^ http://www.ncbi.nlm.nih.gov/
  3. ^ http://www.expasy.ch/sprot/
  4. ^ http://www.genome.jp/kegg/
  5. ^ http://www.cazy.org/

References

Aspeborg, H., Coutinho, P. M., Wang, Y., Brumer, H., and Henrissat, B. (2012). Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5). BMC Evol. Biol. 212:186. doi: 10.1186/1471-2148-12-186

PubMed Abstract | CrossRef Full Text | Google Scholar

Baumann, M. J., Eklof, J. M., Michel, G., Kallas,Å. M., Teeri, T. T., Czjzek, M., et al. (2007). Structural evidence for the evolution of xyloglucanase activity from xyloglucan endo-transglycosylases: biological implications for cell wall metabolism. Plant Cell 19, 1947–1963. doi: 10.1105/tpc.107.051391

PubMed Abstract | CrossRef Full Text | Google Scholar

Blackman, L. M., Cullerne, D. P., and Hardham, A. R. (2014). Bioinformatic characterisation of genes encoding cell wall degrading enzymes in the Phytophthora parasitica genome. BMC Genomics 15:785. doi: 10.1186/1471-2164-15-785

PubMed Abstract | CrossRef Full Text | Google Scholar

Blackman, L. M., Cullerne, D. P., Torreñ, A. P., Taylor, J., and Hardham, A. R. (2015). RNA-Seq analysis of the expression of genes encoding cell wall degrading enzymes during infection of lupin (Lupinus angustifolius) by Phytophthora parasitica. PLoS ONE 10:e0136899. doi: 10.1371/journal.pone.0136899

PubMed Abstract | CrossRef Full Text | Google Scholar

Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The carbohydrate-active enzymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663

PubMed Abstract | CrossRef Full Text | Google Scholar

Christian, C. P., Starr, T. L., and Glass, N. L. (2014). Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Annu. Rev. Phytopathol. 52, 427–451. doi: 10.1146/annurev-phyto-102313-045831

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D. J. (2005). Growth of the plant cell wall. Nat. Rev. Mol. Cell. Biol. 6, 850–861. doi: 10.1038/nrm1746

PubMed Abstract | CrossRef Full Text | Google Scholar

Couturier, M., Navarro, D., Olive, C., Chevret, D., Haon, M., Favel, A., et al. (2012). Post-genomic analyses of fungal lignocellulosic biomass degradation reveal the unexpected potential of the plant pathogen Ustilago maydis. BMC Genomics 13:57. doi: 10.1186/1471-2164-13-57

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y. C., Yang, Z. L., Cui, B. K., Yu, C. J., and Zhou, L. W. (2009). Species diversity and utilization of medicinal mushrooms and fungi in China. Int. J. Med. Mushrooms 11, 287–302. doi: 10.1615/IntJMedMushr.v11.i3.80

CrossRef Full Text | Google Scholar

Esteban, C. I. (2009). Medicinal interest of Poria cocos (= Wolfiporia extensa). Rev. Iberoam. Micol. 26, 103–107. doi: 10.1016/S1130-1406(09)70019-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, Y. L., Lei, P., Tian, T., Yin, L., Chen, D. Q., Chen, H., et al. (2013). Diuretic activity of some fractions of the epidermis of Poria cocos. J. Ethnopharmacol. 150, 1114–1118. doi: 10.1016/j.jep.2013.10.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Q., Jin, Y., Zhang, L., Cheung, P. C. K., and Kennedy, J. F. (2007). Structure, molecular size and antitumor activities of polysaccharides from Poria cocos mycelia produced in fermenter. Carbohydr. Polym. 70, 324–333. doi: 10.1016/j.carbpol.2007.04.015

CrossRef Full Text | Google Scholar

King, B. C., Waxman, K. D., Nenni, N. V., Walker, L. P., Bergstrom, G. C., and Gibson, D. M. (2011). Arsenal of plant cell wall degrading enzymes reflects host preferenceamong plant pathogenic fungi. Biotechnol. Biofuels. 4:4. doi: 10.1186/1754-6834-4-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobira, S., Atsumi, T., Kakiuchi, N., and Mikage, M. (2012). Difference in cultivation characteristics and genetic polymorphism between Chinese and Japanese strains of Wolfiporia cocos Ryvarden et Gilbertson (Poria cocos Wolf). J. Nat. Med. 66, 493–499. doi: 10.1007/s11418-011-0612-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubicek, C. P., Starr, T. L., and Glass, N. L. (2014). Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. Annu. Rev. Phytopathol. 52, 427–451. doi: 10.1146/annurev-phyto-102313-045831

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubo, T., Terabayashi, S., Takeda, S., Sasaki, H., Aburada, M., and Miyamoto, K. (2006). Indoor cultivation and cultural characteristics of Wolfiporia cocos sclerotia using mushroom culture bottles. Biol. Pharm. Bull. 29, 1191–1196. doi: 10.1248/bpb.29.1191

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D., Gong, J., Dai, W., Kang, X., Huang, Z., Zhang, H. M., et al. (2012). The genome of Ganderma lucidum provide insights into triterpense biosynthesis and wood degradation. PLoS ONE 7:e36146. doi: 10.1371/journal.pone.0036146

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakajima, M., Yamashita, T., Takahashi, M., Nakano, Y., and Takeda, T. (2012). A novel glycosylphosphatidylinositolanchored glycoside hydrolase from Ustilago esculenta functions in β-1,3-glucan degradation. Appl. Environ. Microbiol. 78, 5682–5689. doi: 10.1128/AEM.00483-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Rai, R. A., and Agnihotri, J. P. (1971). Influence of nutrition and pH on growth and sclerotiaformation of Sclerotinia sclerotiorum (Lib.) De Bary from Gaillardia pulchella Foug. Mycopathol. Mycol. Appl. 43, 89–95. doi: 10.1007/BF02051505

CrossRef Full Text | Google Scholar

Rios, J. L. (2011). Chemical constituents and pharmacological properties of Poria cocos. Planta Med. 77, 681–691. doi: 10.1055/s-0030-1270823

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, S., Chen, B., Zhou, M., Zhao, X., Xia, H., and Wang, M. (2013). De novo sequencing and transcriptome analysis of Wolfiporia cocos to reveal genes related to biosynthesis of triterpenoids. PLoS ONE 8:e71350. doi: 10.1371/journal.pone.0071350

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, J., Nie, J., Li, D., Zhu, W., Zhang, S., Ma, F., et al. (2014). Characterization and antioxidant activities of degraded polysaccharides from Poria cocos sclerotium. Carbohydr. Polym. 105, 121–126. doi: 10.1016/j.carbpol.2014.01.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Mukerabigwi, J. F., Zhang, Y., Han, L., Jiayinaguli, T., Wang, Q., et al. (2015a). In vivo immunological activity of carboxymethylated-sulfated (1→3)-β-d-glucan from sclerotium of Poria cocos. Int. J. Biol. Macromol. 79, 511–517. doi: 10.1016/j.ijbiomac.2015.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Dong, H., Yan, R., Li, H., Li, P., Chen, P., et al. (2015b). Comparative study of lanostane-type triterpene acids in different parts of Poria cocos (Schw.) Wolf by UHPLC–Fourier transform MS and UHPLC-triple quadruple MS. J. Pharm. Biomed. 102, 203–214. doi: 10.1016/j.jpba.2014.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K. Q., Yin, X. R., Huang, H., Fu, J., Feng, H. G., Wang, Q., et al. (2012). Production status and industrialization development countermeasures of Poria in Hubei Province. Mod. Chin. Med. 14, 24–27.

Google Scholar

Wang, L., Feng, Z., Wang, X., Wang, X., and Zhang, X. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26, 136–138. doi: 10.1093/bioinformatics/btp612

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. Z., Zhang, J., Zhao, Y. L., Li, T., Shen, T., Li, J. Q., et al. (2013). Mycology, cultivation, traditional uses, phytochemistry and pharmacology of Wolfiporia cocos (Schwein.) Ryvarden et Gilb.: a review. J. Ethnopharmacol. 147, 265–276. doi: 10.1016/j.jep.2013.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, X., Xie, J., Cheng, J., Li, G., Yi, X., Jiang, D., et al. (2014). Novel secretory protein Ss-Caf1 of the plant-pathogenic fungus Sclerotinia sclerotiorum is required for host penetration and normal sclerotial development. Mol. Plant Microbe Interact. 27, 40–55. doi: 10.1094/MPMI-05-13-0145-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, Y. M., Chen, J., Lv, Y. L., Liang, H. Q., and Guo, S. X. (2011). Dertermination of optimal carbon source and pH value for sclerotial formation of Polyporus umbellatus under artificial conditions. Mycol. Prog. 10, 121–125. doi: 10.1007/s11557-010-0725-y

CrossRef Full Text | Google Scholar

Xing, Y. M., Chen, J., Song, C., Liu, Y. Y., Guo, S. X., and Wang, C. L. (2013). Nox gene expression and cytochemical localization of hydrogen peroxide in Polyporus umbellatus sclerotial formation. Int. J. Mol. Sci. 14, 22967–22981. doi: 10.3390/ijms141122967

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, F., Xu, B., Zhao, S., Li, J., Yang, Y., Tang, X., et al. (2012). De novo sequencing and analysis of the termite mushroom (Termitomyces albuminosus) transcriptome todiscover putative genes involved in bioactive component biosynthesis. J. Biosci. Bioeng. 114, 228–231. doi: 10.1016/j.jbiosc.2012.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y., Jiang, D., Xie, J., Cheng, J., Li, G., Yi, X., et al. (2012a). Ss-Sl2, a novel cell wall protein with pan modules, is essential for sclerotial development and cellular integrity of Sclerotinia sclerotiorum. PLoS ONE 7:e34962. doi: 10.1371/journal.pone.0034962

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G. J., Wang, M., Huang, J., Yin, Y. L., Chen, Y. J., Jiang, S., et al. (2012b). Deep insight into the Ganoderma lucidum by comprehensive analysis of its transcriptome. PLoS ONE 7:e44031. doi: 10.1371/journal.pone.0044031

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y. Y., Feng, Y. L., Bai, X., Tan, X. J., Lin, R. C., and Mei, Q. (2013a). Ultra performance liquid chromatography-based metabonomic study of therapeutic effect of the surface layer of Poria cocos on adenine-induced chronic kidney disease provides new insight into anti-fibrosis mechanism. PLoS ONE 8:e59617. doi: 10.1371/journal.pone.0059617

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Z., Liu, H., Wang, C., and Xu, J.-R. (2013b). Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics 14:274. doi: 10.1186/1471-2164-14-274

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, W., Wei, W., Fu, Y., Cheng, J., Xie, J., Li, G., et al. (2013). A secretory protein of necrotrophic fungus Sclerotinia sclerotiorum that suppresses host resistance. PLoS ONE 8:e53901. doi: 10.1371/journal.pone.0053901

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Wolfiporia cocos, sclerotial development, transcriptome, CAZymes, differentially expressed genes

Citation: Zhang S, Hu B, Wei W, Xiong Y, Zhu W, Peng F, Yu Y, Zheng Y and Chen P (2016) De Novo Analysis of Wolfiporia cocos Transcriptome to Reveal the Differentially Expressed Carbohydrate-Active Enzymes (CAZymes) Genes During the Early Stage of Sclerotial Growth. Front. Microbiol. 7:83. doi: 10.3389/fmicb.2016.00083

Received: 18 November 2015; Accepted: 18 January 2016;
Published: 03 February 2016.

Edited by:

Andrea Gomez-Zavaglia, Center for Research and Development in Food Cryotechnology, Argentina

Reviewed by:

Shu Shaohua, Huazhong Agricultural University, China
Fanyun Zeng, Chinese Academy of Tropical Agricultural Sciences, China

Copyright © 2016 Zhang, Hu, Wei, Xiong, Zhu, Peng, Yu, Zheng and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wenjun Zhu, ODI4NjIxMDhAcXEuY29t; Ping Chen, Y2hlbnBpbmd2aXAyNEAxNjMuY29t

These authors have contributed equally to this work.

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.