棒瑚菌子实体不同发育时期的转录组分析

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

摘要/Abstract

摘要： 旨为研究棒瑚菌子实体幼菇期（CP-1）和成熟期（CP-2）基因表达的差异,探索棒瑚菌在生长发育过程中的关键基因和关键代谢通路。采用RNA-Seq测序技术对棒瑚菌子实体进行测序并进行生物信息学分析。测序结果显示,共获得8.78 G（CP-1）和8.28 G（CP-2）的分析数据（Clean reads）。KOG功能注释和GO富集分析显示新陈代谢过程、细胞和细胞组分以及氧化还原酶活性是棒瑚菌子实体发育过程中的关键基因门类。基因表达水平分析结果显示,多酚氧化酶、核糖体、糖苷水解酶家族等基因为棒瑚菌子实体幼菇期（CP-1）和成熟期（CP-2）的主要高表达基因,在蛋白质和细胞壁合成中起重要作用。差异表达基因分析结果显示氧化磷酸化和三羧酸循环是棒瑚菌子实体幼菇期（CP-1）生长发育的关键代谢途径,且环境含氧量对棒瑚菌的生长发育有重要影响;同时棒瑚菌成熟期（CP-2）调节一碳化合物浓度和应对各种环境胁迫的分子机制发生了改变。KEGG 通路富集结果显示,MAPK信号通路是参与调节棒瑚菌子实体生长发育的关键信号通路,在调控菌丝的极性生长、细胞周期和真菌细胞形态转换中起重要作用。

关键词: 棒瑚菌, RNA-Seq测序, 差异表达基因, 代谢通路

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

唐贤, 丁祥, 董明明, 朱淼, 宋志强, 侯怡铃. 棒瑚菌子实体不同发育时期的转录组分析[J]. 生物技术通报, 2019, 35(10): 119-129.

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.

参考文献

[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.