Transcriptome Analysis of Clavariadelphus pistillaris Fruiting Bodies at Different Development Stages

doi:10.13560/j.cnki.biotech.bull.1985.2019-0233

Abstract

Abstract: This work aims to study the differential expressed genes（DEGs）in Clavariadelphus pistillaris fruiting bodies between the period of juvenile mushroom（CP-1）and mature mushroom（CP-2）and to explore the key genes and metabolic pathways during the growth and development of C. pistillaris. RNA-Seq sequencing technology was used to sequence the C. pistillaris fruiting bodies,and their bioinformatics was analyzed. The sequencing results demonstrated that a total of 8.78 G（CP-1）and 8.28 G（CP-2）clean reads were obtained. Functional annotations of gene KOG and GO enrichment analysis showed that the metabolic processes,cellular and cellular parts,and oxidoreductase activities were the key gene categories in the development of C. pistillaris fruiting bodies. The analysis of gene expression level revealed that the main high-expression genes in CP-1 and CP-2 were polyphenol oxidase,ribosome,and glycoside hydrolase family,etc,which played an important role in protein and cell wall synthesis. The results of DEGs showed that oxidative phosphorylation and citrate cycle（TCA cycle）were the key metabolic pathways for the growth and development of CP-1,and the environmental oxygen content had an important effect on the growth and development of C. pistillari;meanwhile,the molecular mechanism of regulating the concentration of C1 compound and the response to various environmental stresses at CP-2 changed. KEGG pathway enrichment showed that MAPK was a key signaling pathway involved in regulating the growth and development of C. pistillaris fruiting bodies,and it played an important role in regulating the polar growth of mycelia,cell cycle and the morphological transformation of fungal cells.

Key words: Clavariadelphus pistillaris, RNA-Seq sequence, differentially-expressed gene, metabolic pathway

TANG Xian, DING Xiang, DONG Ming-ming, ZHU Miao, SONG Zhi-qiang, HOU Yi-ling. Transcriptome Analysis of Clavariadelphus pistillaris Fruiting Bodies at Different Development Stages[J]. Biotechnology Bulletin, 2019, 35(10): 119-129.

References

[1] 李玉. 中国食用菌产业发展现状、机遇和挑战——走中国特色菇业发展之路, 实现食用菌产业强国之梦[J]. 菌物研究, 2018, 16(3):5-11.
[2] 王建瑞, 刘宇, 图力古尔. 山东省大型真菌物种濒危程度与优先保育评价[J]. 生态学报, 2015, 35(3):837-848.
[3] 贾昌路, 张瑶, 朱玲, 等. 转录组测序技术在生物测序中的应用研究进展[J]. 分子植物育种, 2015, 13(10):2388-2394.
[4] 荣成博, 赵爽, 宋爽, 等. 白灵侧耳子实体发育相关基因的挖掘[J]. 生物技术通报, 2018, 34(4):115-120.
[5] 吴小梅, 张昕, 李南羿. 双孢蘑菇子实体不同发育时期的转录组分析[J]. 菌物学报, 2017, 36(2):193-203.
[6] 施肖堃, 蔡志欣, 郭仲杰, 等. 双孢蘑菇As2796子实体发育转录组测序分析[J]. 福建农业学报, 2018, 33(3):282-287.
[7] 李帆, 陈利丁, 艾柳英, 等. 刺芹侧耳不同发育时期的转录组差异分析[J]. 菌物学报, 2018, 37(12):1586-1597.
[8] 卯晓岚. 中国大型真菌[M]. 郑州:河南科学技术出版社, 2000.
[9] Kranabetter JM, Friesen J, Gamiet S, et al.Epigeous fruiting bodies of ectomycorrhizal fungi as indicators of soil fertility and associated nitrogen status of boreal forests[J]. Mycorrhiza, 2009, 19(8):535-548.
[10] Pereira E, Barros L, Martins A, et al.Towards chemical and nutritional inventory of Portuguese wild edible mushrooms in different habitats[J]. Food Chemistry, 2012, 130(2):394-403.
[11] Lingham RB, Silverman KC, Jayasuriya H, et al.Clavaric acid and steroidal analogues as Ras- and FPP-directed inhibitors of human farnesyl-protein transferase[J]. Journal of Medicinal Chemistry, 1998, 41(23):4492-4501.
[12] Yamac M, Bilgili F.Antimicrobial activities of fruit bodies and/or mycelial cultures of some mushroom isolates[J]. Pharmaceutical Biology, 2006, 44(9):660-667.
[13] 赵洁, 杨淑达, 唐丽萍. 棒瑚菌系统分类研究进展[J]. 昆明医科大学学报, 2016, 37(3):231-233
[14] 袁明生, 孙佩琼. 中国大型真菌彩色图谱[M]. 成都:四川科学技术出版社, 2013.
[15] Cock PJA, Fields CJ, Goto N, et al.The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants[J]. Nucleic Acids Research, 2010, 38(6):1767-1771.
[16] Grabherr MG, Haas BJ, Yassour M, et al.Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology, 2011, 29(7):644-652.
[17] Chen G, Shi T, Shi L.Characterizing and annotating the genome using RNA-seq data[J]. Science China Life Sciences, 2017, 60(2):116-125.
[18] Mccarthy DJ, Chen Y, Smyth GK.Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation[J]. Nucleic Acids Research, 2012, 40(10):4288-4297.
[19] Young MD, Wakefield MJ, Smyth GK, et al.Gene ontology analysis for RNA-seq:accounting for selection bias[J]. Genome Biology, 2010, 11(2):R14.
[20] Kanehisa M, Sato Y, Kawashima M, et al.KEGG as a reference resource for gene and protein annotation[J]. Nucleic Acids Research, 2016, 44(D1):D457-D462.
[21] Wang J, Liu B, Xiao Q, et al.Cloning and expression analysis of litchi(Litchi Chinensis Sonn. )polyphenol oxidase gene and relationship with postharvest pericarp browning[J]. PLoS One, 2014, 9(4):e93982.
[22] LI NY, Cai WM, Jin QL, et al.Molecular cloning and expression of polyphenoloxidase genes from the mushroom, Agaricus bisporus[J]. Agricultural Sciences in China, 2011, 10(2):185-194.
[23] Kim HI, Lee CS, Park YJ.Further characterization of hydrophobin genes in genome of Flammulina velutipes[J]. Mycoscience, 2016, 57(5):320-325.
[24] 黄姗, 王中康, 陈环, 等. 莱氏野村菌疏水蛋白基因Nrhyd的克隆及其表达特征分析[J]. 菌物学报, 2012, 31(3):350-358.
[25] Karkehabadi S, Hansson H, Mikkelsen NE, et al.Structural studies of a glycoside hydrolase family 3 β-glucosidase from the model fungus Neurospora crassa[J]. Acta Crystallographica Section F:Structural Biology Communications, 2018, 74(12):787-796.
[26] Hayouka Z, Thomas NC, Mortenson DE, et al.Quasiracemate crystal structures of magainin 2 derivatives support the functional significance of the phenylalanine zipper motif[J]. Journal of the American Chemical Society, 2015, 137(37):11884-11887.
[27] Boczek NJ, Gomez-Hurtado N, Ye D, et al.Spectrum and prevalence of CALM1-, CALM2-, and CALM3-Encoded calmodulin variants in long QT syndrome and functional characterization of a novel long QT syndrome-associated calmodulin missense variant, E141G[J]. Circulation Cardiovascular Genetics, 2016, 9(2):136-146.
[28] 李天明, 朱灵桓, 刘天佳, 等. Pseudom onas sp. 101甲酸脱氢酶基因的合成、表达及定点饱和突变[J]. 生物技术通报, 2013(5):105-110.
[29] 朱钦士. 正转和反转的三羧酸循环[J]. 生物学通报, 2015, 50(1):16-19.
[30] Luke NS, Devito MJ, Portier CJ, et al.Employing a mechanistic model for the MAPK pathway to examine the impact of cellular all or none behavior on overall tissue response[J]. Dose Response, 2010, 8(3):347-367.
[31] Atay O, Skotheim JM.Spatial and temporal signal processing and decision making by MAPK pathways[J]. The Journal of Cell Biology, 2017, 216(2):317-330.
[32] Brand AC, Morrison E, Milne S, et al.Cdc42 GTPase dynamics control directional growth responses[J]. Proceedings of the National Academy of Sciences, 2014, 111(2):811-816.
[33] Woods B, Lew DJ.Polarity establishment by Cdc42:Key roles for positive feedback and differential mobility[J]. Small GTPases, 2016, 10(2):130-137.