光照脅迫下廣葉繡球菌菌絲基因差異表達的初步分析*

肖冬來，張 迪，馬 璐，王宏雨，林衍銓

（福建省農業科學院食用菌研究所，福建 福州 350014）

光照脅迫下廣葉繡球菌菌絲基因差異表達的初步分析*

肖冬來，張 迪，馬 璐，王宏雨，林衍銓**

（福建省農業科學院食用菌研究所，福建 福州 350014）

利用高通量測序技術對黑暗和光照下廣葉繡球菌（Sparassis latifolia）菌絲進行RNA-Seq分析，比較光照脅迫下的差異表達基因，并對差異基因進行GO功能分析和KEGG pathway分析。結果表明，光照脅迫下顯著性差異表達的基因328個，其中，上調、下調表達的基因數分別為157個和171個。差異表達基因的GO注釋表明，生物學過程主要為細胞過程、代謝過程、單一生物過程，而分子功能主要為催化、結合和轉運子活性。KEGG pathway功能分析結果表明，差異基因主要集中在次級代謝產物合成、淀粉和蔗糖代謝和谷胱甘肽代謝通路中。

廣葉繡球菌；轉錄組；GO功能分析；差異表達基因

廣葉繡球菌（Sparassis latifolia Y.C.Dai&Zheng Wang）[1]，隸屬于多孔菌目（Polyporales） 繡球菌科（Sparassidacea） 繡球菌屬（Sparassis），是1種珍稀食（藥）用真菌，廣葉繡球菌中的β-葡聚糖能提高人體免疫力，具有抗癌、防癌特殊功效[2]。光照在食用菌菌絲生長、原基形成和產量、子實體形態等[3-4]過程都發揮著重要的調控作用。

廣葉繡球菌對光照的需求要高于其它的食用菌，適宜的光照是原基分化與子實體生長發育的必要條件，不同的光質和光量影響廣葉繡球菌菌絲生長速度和原基形成率[5]。Yang等[6]克隆了1個廣葉繡球菌藍光受體基因，光照可誘導該基因的表達，為光照誘導子實體發育的研究提供參考。本研究利用轉錄組測序技術分析了光照條件下廣葉繡球菌基因差異表達，以期進一步為廣葉繡球菌的光響應機制研究提供科學參考。

1 材料和方法

1.1 供試菌株和參考基因組

供試廣葉繡球菌菌株“閩繡1號”為福建省農業科學院食用菌研究所保藏。參考基因組為本單位測序完成的廣葉繡球菌菌株“閩繡1號”基因組，該數據尚未公開。

1.2 樣品制備

利用打孔器在活化的繡球菌PDA培養皿上，呈圓周形選取菌齡一致，直徑約7 mm的菌塊接種至表面鋪有玻璃紙的9 cm PDA培養皿。25℃黑暗培養20 d后，將培養皿置于裝有LED白色光源的培養箱中25℃培養48 h（調整透過培養皿蓋的光強度為300 lx），黑暗培養的菌絲為對照組。無菌玻璃棒刮取菌絲后液氮急凍保存備用。

1.3 差異表達基因分析

光照組和黑暗處理組樣品送至北京安諾優達基因科技有限公司提取總RNA、構建文庫，利用Illumina HiSeq平臺測序。基因表達量利用FPKM（fragments per kilobase per millon mapped fragments）值進行評估。采用DEGSeq v1.18.0[7]進行基因差異表達分析，比較光照組與黑暗組，并選取|log2Ratio|≥1和q＜0.05的基因作為差異表達基因，得到上調、下調基因個數。根據差異表達基因在每個樣品里的表達量，取以2為底的對數后，利用R軟件（v3.1.1）對差異表達基因進行分層聚類分析。

1.4 GO功能分類注釋與顯著性富集分析

根據GO數據庫（Gene Ontology，http://www.geneontology.org/），利用軟件Blast2GO得到差異表達基因對應的GO條目。針對GO數據庫中第2層的條目，統計差異表達基因在該條目里的個數，并計算百分比。根據計算每個條目的基因數目，然后應用超幾何檢驗，找出與整個基因組背景相比，差異表達基因顯著富集的GO條目。計算得到的p-value通過校正之后，以q＜0.05為閾值，滿足此條件的GO條目定義為在差異表達基因中顯著富集的GO條目。

1.5 代謝通路分析與顯著性富集分析

利用 KEGG（Kyoto Encyclopedia of Genes and Genomes）數據庫（http://www.kegg.jp/）分析差異基因參與的pathway，對KEGG中每個Pathway應用超幾何檢驗進行富集分析，找出差異表達基因顯著性富集的Pathway，計算得到的p-value通過校正之后，以q＜0.05為閾值，滿足此條件的pathway定義為在差異表達基因中顯著富集的pathway。

2 結果與分析

2.1 差異基因的篩選

光照組和黑暗組共有11 029個轉錄本可定位到參考基因組（廣葉繡球菌基因組測序結果預測為12 471個基因，數據尚未公開）。以黑暗組為對照，共篩選出差異表達基因328個，其中，上調、下調表達的基因數分別為157個和171個，差異表達顯著的部分基因見表1。

表1 部分顯著差異表達的基因Tab.1 Some significantly differentially expressed genes

2.2 差異基因的GO功能分析

差異表達基因的GO功能分析和GO注釋富集分析見圖1、表2。

差異表達基因的GO功能分析顯示，在生物學過程 （biological process） 分類中，細胞過程（cellular process）、代謝過程 （metabolic process） 和單一生物過程（single-organism process）所占比例較高，分別占總差異蛋白的19.8%、21.6%和16.5%。在細胞組分分類中差異表達基因在細胞組分（cell part）、細胞器（organelle）和膜組分（membrane part）類型所占比例較高，分別占總差異蛋白的16.5%、7.3%和5.5%。在分子功能組分中，催化（catalytic）、結合（binding）和轉運子（transporter）活性類型所占比例較高，分別占總差異蛋白的28.7%、16.2%和3.0%。

GO功能顯著性富集分析表明，顯著富集的功能有：生物學過程下的多糖分解過程（polysaccharide catabolic process）、碳水化合物代謝過程（carbohy－drate metabolic process） 和多糖代謝過程（polysaccharide metabolic process）；分子功能下的水解酶活性（hydrolase activity）、氧化還原酶活性（oxidoreductase activity）、催化活性 （catalytic activity）、β-甘露糖苷酶活性（beta-mannosidase activity）、單加氧酶活性（monooxygenase activity） 和半乳糖苷酶的活性（galactosidase activity）；細胞組分下的胞外區域 （extracellular region）。

圖1 差異表達基因的GO功能分析Fig.1 Functional annotation of differentially expressed genes according to GO analysis

表2 GO注釋富集分析Tab.2 GO enrichment analysis

2.3 差異基因的KEGG pathway分析

對差異基因進行pathway分析發現，差異基因可注釋到30條pathway上，差異基因數目較多的pathway分別為次級代謝產物合成（biosynthesis of secondary metabolite，涉及10個差異基因）、淀粉和蔗糖代謝（starch and sucrose metabolism，涉及4個差異基因）、谷胱甘肽代謝（glutathione metabolism，涉及3個差異基因）、半胱氨酸和蛋氨酸代謝（cysteine and methionine metabolism，涉及3個差異基因）。經KEGG pathway富集分析顯示，未檢測到顯著富集的pathway。

3 結論與討論

本研究通過轉錄組測序分析，得到廣葉繡球菌在光照和黑暗培養下差異基因328個，這些基因可能參與了廣葉繡球菌對光照的響應調控。廣葉繡球菌的發育過程是一個多因素調控的生命過程，除了自身的遺傳特性作為其調控發育的基礎外，外界栽培條件如栽培基質、溫度、光照等均可影響其生長發育。目前，關于光照調控廣葉繡球菌發育的分子機制研究還很少[6]，本研究將所有差異表達基因進行GO注釋和功能聚類分析，在顯著富集的生物學功能分類中，多糖分解、碳水化合物代謝和多糖代謝涉及的基因下調數量大于上調數量，推測光照可能影響了繡球菌菌絲碳水化合物降解能力。胞外酶活性在一定程度上可反映菌株生長的快慢和對基質的利用率[8]。應正河等[5]的研究表明繡球菌菌絲在黑暗條件下生長速度最快，顯著高于不同光質下菌絲的生長速度。

谷胱甘肽、蛋氨酸和甲硫氨酸參與了多種環境脅迫的響應，同時蛋氨酸和甲硫氨酸對細胞的分化和功能的產生有著重要的影響[9-11]。KEGG代謝通路分析結果顯示谷胱甘肽代謝、半胱氨酸和蛋氨酸代謝通路中均含有3個關鍵差異基因。其中2個谷胱甘肽S-轉移酶在光照誘導下表現上調，而半胱氨酸和蛋氨酸代謝通路中編碼3個關鍵基因腺苷高半胱氨酸酶（adenosylhomocysteinase）、s-腺苷甲硫氨酸合成酶（s-adenosylmethionine synthetase）、氨基環丙基甲酸甲酯氧化酶（aminocyclopropanecarboxylate oxidase） 均表現下調，相關差異基因在光響應中的分子機制還需進一步研究。

綜上所述，本研究應用轉錄組測序技術，初步分析了光照誘導下廣葉繡球菌菌絲基因表達變化情況、差異表達基因的GO分類和代謝過程，研究結果為進一步分析廣葉繡球菌光響應的分子機制，探討差異基因與菌絲發育的調控關系提供科學數據。

[1]Dai YC,Wang Z,Binder M,et al.Phylogeny and a new species of Sparassis(Polyporales,Basidiomycota):evidence from mitochondrial atp6,nuclear rDNA and rpb2 genes[J].Mycologia,2006,98(4):584-92.

[2]Kimura T.Natural products and biological activity of the pharmacologically active cauliflower mushroom Sparassis crispa[J].Bio Med Research International,2013,2013 (11):167-169.

[3]黃兵.不同光質LED對白靈菇商品性狀及產量的影響[D].長春：吉林農業大學，2015.

[4]白雪萍，趙麗，周軍，等.杏鮑菇子實體發育形態的光調控機制研究[J].湖北農業科學，2016（9）：2275-2278.

[5]應正河，林衍銓，馬璐，等.不同光質光量對繡球菌菌絲生長及原基形成的影響[J].福建農業學報，2013，28（6）：538-540.

[6]Yang C,Ma L,Ying ZH,et al.Sequence analysis and expression of a blue-light photoreceptor gene,slwc-1 from the cauliflower mushroom Sparassis latifolia[J].Current Microbiology,2017,74(4):469-475.

[7]Wang L,Feng Z,Wang X,et al.DEGseq:an r package for identifying differentially expressed genes from RNA-seq data[J].Bioinformatics,2010,26(1):136-138.

[8]安琪，吳雪君，吳冰，等.不同碳源和氮源對金針菇降解木質纖維素酶活性的影響[J].菌物學報，2015，34（4）：761-771.

[9]肖清鐵，王經源，鄭新宇，等.水稻根系響應鎘脅迫的蛋白質差異表達[J].生態學報，2015，35（24）：8276-8283.

[10]呂煥青，王志敏，湯青林，等.多胺生物合成途徑中兩個關鍵酶基因研究進展[J].生物技術通報，2015，31（2）：61-64.

[11]Jiang HW,Liu M,Chen IC,et al.A glutathione S-transferase regulated by light and hormones participates in the modulation of Arabidopsis seedling development[J].Plant P-hysiology,2010,154(4):1646-1658.

英國：橡樹下發現“黑色鉆石”之稱的珍稀黑松露

前不久，英國德文郡史蒂夫先生無意在一棵橡樹下挖到了10多顆有“廚師的黃金”之稱的珍稀黑松露。據了解，黑松露是1種生長于地下的野生食用真菌，含有豐富的蛋白質、維生素、18種氨基酸以及人體必需微量元素，具有極高的營養保健價值。黑松露每年在全世界產量非常稀少，因此價格相當昂貴，堪比黃金。在中國西南藏族、彝族等少數民族地區，黑松露也被當地人作為增強體力、治療疾病的“黑色鉆石”。

在國外，黑松露常見于在法國、意大利和西班牙的森林中，每千克價格高達2 000英鎊，除此品種之外，原產于于意大利的白松露也尤為珍稀，2007年一場拍賣會上，一塊重達3.3磅（約1.5 kg）的白松露以16.5萬英鎊價格成交。

中國食用菌商務網

2017.09.04

Preliminary Study on Differentially Expressed Genes of Sparassis latifolia under Light Inducing

XIAO Dong-lai,ZHANG Di,MA Lu,WANG Hong-yu,LIN Yan-quan
(Institute of Edible＆Medicinal Fungi,Fujian Academy of Agricultural Sciences,Fuzhou 350014,China)

The transcriptomes of Sparassis latifolia mycelium under light or dark inducing were compared using the whole genome sequence of Sparassis latifolia as a reference.Both light and dark inducing mycelium were subjected to RNA sequencing and exposed to GO function and KEGG pathway analysis.The results showed that a total of 328 genes were identified as believable.Among them,157 genes were up-regulated and 171 genes were down-regulated.The GO classification showed that,the differentially expressed genes were mostly enriched in cellular process,metabolic process,single-organism process biological pathway and catalytic,binding,transporter molecular functions.KEGG pathway analysis showed that the differentially expressed genes were involved in biosynthesis of secondary metabolites,starch and sucrose metabolism,glutathione metabolism.

Sparassis latifolia;transcriptome;GO function analysis;differentially expressed gene

S646.9

A

1003-8310（2017）05-0060-04

10.13629/j.cnki.53-1054.2017.05.014

福建省屬公益類科研院所基本科研專項（2015R1020-5）；福建省自然科學基金（2016J01133）。

肖冬來（1981-），男，博士，助理研究員，主要從事食用菌栽培生理研究。E-mail:xdljiangsu@163.com

**通信作者：林衍銓（1963-），男，大專，研究員，主要從事食用菌栽培研究。E-mail:lyq-406@163.com

2017-07-13