銀耳二型態(tài)編碼MAPK關(guān)鍵差異基因的克隆及生物信息學(xué)分析

徐濤 孔旭強(qiáng) 占觀平 陳雅芳 曹繼璇 李亞星 孫淑靜

銀耳二型態(tài)編碼MAPK關(guān)鍵差異基因的克隆及生物信息學(xué)分析

徐濤 孔旭強(qiáng) 占觀平 陳雅芳 曹繼璇 李亞星 孫淑靜*

（福建農(nóng)林大學(xué)生命科學(xué)學(xué)院，福建 福州 350002）

銀耳是我國特色的食用菌品種，其二型態(tài)轉(zhuǎn)換與實際生產(chǎn)密切相關(guān)并制約著育種工作的進(jìn)展，而絲裂原活化蛋白激酶（Mitogen activated protein kinase，MAPK）途徑在真菌二型態(tài)轉(zhuǎn)換中具有重要調(diào)控作用。根據(jù)轉(zhuǎn)錄組測序數(shù)據(jù)克隆了銀耳二型態(tài)編碼MAPK的關(guān)鍵差異基因，包括基因、基因和基因分析表明，、、基因的基因組DNA序列全長分別為1 721 bp、1 865 bp、1 903 bp，ORF序列全長分別為1 191 bp、1 281 bp、1 149 bp，分別含有10、12、16個內(nèi)含子，編碼396、426、382個氨基酸，相對分子質(zhì)量在43.61 kDa～48.15 kDa間，理論等電點在5.82～6.89。序列同源性分析表明，、、基因所編碼的氨基酸序列分別與中同源蛋白的相似性最高。結(jié)構(gòu)域分析表明，TfFus3、TfSlt2、Kss1蛋白都具有典型的蛋白激酶保守結(jié)構(gòu)域且同屬于PKc_like super family。本研究結(jié)果為進(jìn)一步探究MAPK信號通路在銀耳二型態(tài)轉(zhuǎn)換中的作用奠定了基礎(chǔ)。

銀耳；絲裂原活化蛋白激酶；二型態(tài)；生物信息學(xué)

銀耳（Berk）隸屬于真菌門、擔(dān)子菌綱、銀耳目、銀耳科、銀耳屬，為我國久負(fù)盛名的食藥用菌[1]，素有“菌中之冠”的美稱。銀耳是典型的二型態(tài)真菌，具有酵母型和菌絲型2種生長階段。其酵母型包括由擔(dān)孢子芽殖生成的單核酵母狀孢子和由雙核菌絲或子實體在游離水過多、高溫、機(jī)械刺激等不良條件下形成的雙核酵母狀孢子[2]。而其菌絲型則包括單核菌絲、雙核菌絲和假菌絲3種類型[3]。具有不同極性的銀耳擔(dān)孢子可在和雙因子控制的四極性交配型系統(tǒng)的控制下形成具鎖狀聯(lián)合的雙核菌絲[4]，雙核酵母狀孢子也可在適宜條件下萌發(fā)形成雙核菌絲和假菌絲。銀耳酵母狀孢子直接與香灰菌配對時萌發(fā)率低、很難形成白毛團(tuán)，更無法獲得子實體，而用銀耳菌絲體與香灰菌配對時則可大大提高銀耳的成品率，因此探究銀耳酵母狀孢子向菌絲轉(zhuǎn)變的調(diào)控機(jī)制對今后銀耳的育種工作及生產(chǎn)至關(guān)重要。

真菌二型態(tài)轉(zhuǎn)變受外界環(huán)境因子和內(nèi)部信號轉(zhuǎn)導(dǎo)途徑的影響，其中絲裂原活化蛋白激酶（Mitogen activated protein kinase，MAPK）信號通路在真核細(xì)胞中較為保守，與真菌二型態(tài)的轉(zhuǎn)變密切相關(guān)[5]。1994年Liu等[6]最先發(fā)現(xiàn)MAPK途徑在釀酒酵母（）形態(tài)轉(zhuǎn)化中有重要作用，后又發(fā)現(xiàn)該途徑與白念珠菌（）的二型態(tài)也有關(guān)。在模式真菌釀酒酵母中，MAPK信號通路主要包括由Fus3介導(dǎo)的信息素反應(yīng)（pheromone response），Kss1介導(dǎo)的絲化和侵襲性生長（filamentation and invasive growth），Slt2介導(dǎo)的細(xì)胞壁完整性（cell wall integrity），Hog1介導(dǎo)的高滲透壓應(yīng)激反應(yīng)（high-osmolarity stress response）[7]。、作為2個同源基因主要調(diào)控釀酒酵母有性生殖和菌絲生長過程[8]。新生隱球菌（）的二型態(tài)轉(zhuǎn)變一直被認(rèn)為是一個依賴信息素的過程，其交配反應(yīng)通過信息素受體級聯(lián)反應(yīng)介導(dǎo)，其中與釀酒酵母、分別具有52%和51%的相似性[9]。白念珠菌中與菌絲生長有關(guān)的MAPK通路為Cek1和Cek2信號轉(zhuǎn)導(dǎo)通路，Cek1與Kss1同屬一個家族，Cek2與Fus3有55%的相似性，其可恢復(fù)釀酒酵母Fus3/Kss1缺失菌株的交配缺陷[10]。此外，有研究發(fā)現(xiàn)白念珠菌中（同源基因）雙敲除菌株無法在Spider培養(yǎng)基上形成侵襲性菌絲，而在釀酒酵母中過表達(dá)基因則可促進(jìn)在氮饑餓條件下假菌絲的形成[11]。/同源基因在許多植物病原體中也同樣起重要作用，禾谷鐮刀菌（）中同源基因的缺失會導(dǎo)致氣生菌絲生長減少，減弱了對植株的侵染能力，并且它對于有性生殖過程中的雌性生殖能力是必不可少的[12]。玉米小斑病菌（）中的和釀酒酵母（）的氨基酸序列有59%的相似性。缺失突變體的氣生菌絲發(fā)育不良，無附著孢，侵染能力嚴(yán)重降低，不產(chǎn)分生孢子，菌絲內(nèi)部有自溶現(xiàn)象[13]。

銀耳作為典型的二型態(tài)真菌可在酵母態(tài)與菌絲態(tài)之間相互轉(zhuǎn)變，且一旦形成酵母狀芽孢便很難形成菌絲。目前有關(guān)銀耳二型態(tài)的研究多集中于外界環(huán)境條件上，對二者轉(zhuǎn)變過程中的關(guān)鍵基因及其調(diào)節(jié)網(wǎng)絡(luò)知之甚少。本課題組前期研究發(fā)現(xiàn)，銀耳在二型態(tài)轉(zhuǎn)變過程中MAPK信號通路中有多個基因出現(xiàn)了差異性表達(dá)，本研究根據(jù)銀耳酵母狀芽孢與菌絲體的轉(zhuǎn)錄組數(shù)據(jù)設(shè)計引物，對銀耳二型態(tài)編碼MAPK的關(guān)鍵差異基因、、進(jìn)行克隆及生物信息學(xué)分析，為研究該通路在銀耳二型態(tài)轉(zhuǎn)變中的作用奠定基礎(chǔ)。

1 材料與方法

1.1 材料

（1）供試菌株。銀耳“繡銀1號”，編號TWW01-AX，由福建農(nóng)林大學(xué)生命科學(xué)學(xué)院特色食用菌品種創(chuàng)新與代謝工程團(tuán)隊選育并保存。

（2）供試培養(yǎng)基。PDA加富液體培養(yǎng)基：馬鈴薯200 g，葡萄糖20 g，蛋白胨3 g，酵母粉3 g，MgSO4·7H2O 1.5 g，KH2PO41.5 g，水1 L。

（3）試劑。Ezup柱式真菌基因組DNA抽提試劑盒、DEPC水、氨芐青霉素均購自生工生物工程（上海）股份有限公司；?DNA Polymerase、2K DNA Marker、Trans1-T1感受態(tài)細(xì)胞均購自北京全式金生物技術(shù)有限公司；TRIzol（Invitrogen）；Prime Script?RT regent Kit With gDNA Eraser反轉(zhuǎn)錄試劑盒（TaKaRa）； Hieff Clone? Zero TOPO-Blunt Simple Cloning Kit（上海翊圣生物科技有限公司）；D2500 Gel Extraction Kit（OMEGA），其余試劑均為國產(chǎn)分析純。

1.2 方法

（1）銀耳基因組DNA的提取。將銀耳酵母狀芽孢接種于PDA加富液體培養(yǎng)基中，25 ℃、150 r/min培養(yǎng)7 d，12 000 r/min離心10 min收集菌體，根據(jù)Ezup 柱式真菌基因組 DNA 抽提試劑盒說明書提取基因組DNA。

（2）銀耳總RNA的提取及cDNA合成。用Trizol法提取銀耳酵母狀芽孢總RNA并根據(jù)PrimeScript?RT regent Kit With gDNA Eraser反轉(zhuǎn)錄試劑盒說明書合成ORF模板。

（3）銀耳二型態(tài)差異基因的克隆及測序。根據(jù)銀耳酵母狀芽孢與菌絲體轉(zhuǎn)錄組數(shù)據(jù)庫中各基因在NR及KEGG數(shù)據(jù)庫的注釋結(jié)果，運用NCBI的ORFfinder功能搜尋MAPK信號通路中富集到的含有完整開放閱讀框的各差異基因序列，用Primer Premier 6軟件設(shè)計引物（表1）。分別以銀耳基因組DNA及總ORF為模板，根據(jù)DNA Polymerase說明書設(shè)置各基因的PCR反應(yīng)體系及程序（表2），并擴(kuò)增各基因序列。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳后，用D2500 Gel Extraction Kit試劑盒對目的片段進(jìn)行回收純化，連接至pESI-Blunt simple vector中并轉(zhuǎn)入Trans1-T1感受態(tài)細(xì)胞，挑選陽性菌落，送生工生物工程（上海）股份有限公司進(jìn)行測序。

表1 引物序列

表2 Hieff Canace? High-Fidelity DNA Polymerase PCR反應(yīng)體系

注：基因組DNA 50 ng～200 ng，質(zhì)粒DNA 10 pg～20 ng，cDNA 1～5 μL；PCR反應(yīng)程序為98 ℃預(yù)變性3 min，98 ℃變性10 s，60 ℃退火20 s，72 ℃延伸（速度為30 s/kb），35個循環(huán)，72 ℃繼續(xù)延伸5 min，4 ℃保存。

（4）差異基因的生物信息學(xué)分析。將克隆得到的各基因序列，通過在線網(wǎng)站（http://gsds.cbi. pku.edu.cn/index.php）繪制基因的結(jié)構(gòu)圖；采用ProtParam（http://www.expasy.org/）分析蛋白質(zhì)的分子量大小、等電點等理化性質(zhì)；利用Expasy在線服務(wù)器上的SOPMA（https://npsa-prabi.ibcp.fr/cgi-bin/ npsa_automat.pl?page=/NPSA/ npsa_sopma.html）分析蛋白質(zhì)的二級結(jié)構(gòu)；利用在線網(wǎng)站（https:// swissmodel.expasy.org/interactive）預(yù)測蛋白質(zhì)的三級結(jié)構(gòu)；蛋白質(zhì)信號肽預(yù)測與分析采用SignalP（http://www.cbs.dtu.dk/services/SignalP/）；蛋白質(zhì)跨膜結(jié)構(gòu)分析采用TMHMM程序（http://www.cbs.dtu.dk/services/TMHMM/）；根據(jù)各基因的Blastp結(jié)果下載不同物種的同源蛋白序列，同時在NCBI搜索下載該基因酵母的蛋白序列。用ClustalX軟件將上述蛋白質(zhì)序列進(jìn)行同源比對，參數(shù)為默認(rèn)值，以同源比對數(shù)據(jù)為基礎(chǔ)，采用MEGA7軟件中的近鄰法，Bootstrap值設(shè)置為1 000，其余參數(shù)為默認(rèn)值，構(gòu)建系統(tǒng)進(jìn)化樹。

2 結(jié)果與分析

2.1 TfFus3、TfSlt2、Kss1基因組DNA及ORF的擴(kuò)增

分別以銀耳酵母狀芽孢基因組DNA及cDNA為模板，以引物-F、-R克隆銀耳基因、以引物-F、-R克隆銀耳基因、以引物-F、-R克隆銀耳基因。cDNA擴(kuò)增產(chǎn)物電泳結(jié)果（圖1），擴(kuò)增得到的基因組DNA序列全長1 721 bp，開放閱讀框（ORF）序列全長1 191 bp；基因組DNA序列全長1 865 bp，ORF序列全長1 281 bp；基因組DNA序列全長1 903 bp，ORF序列全長1 149 bp。

2.2 TfFus3、TfSlt2、Kss1基因結(jié)構(gòu)及蛋白質(zhì)理化性質(zhì)

基因測序結(jié)果（圖2）顯示，3個基因的基因組序列全長分別為1 721 bp、1 865 bp、1 903 bp，ORF全長分別為1 191 bp、1 281 bp、1 149 bp，分別含有10、12、16個內(nèi)含子。其中，基因共編碼396個氨基酸，分子式為C2030H3156N542O591S16，相對分子質(zhì)量45.12 kDa，其中帶負(fù)電荷的殘基（Asp + Glu）49個，帶正電荷的殘基（Arg + Lys）48個，理論等電點為（pI）6.89，不穩(wěn)定性系數(shù)45.09（＞40）為不穩(wěn)定蛋白，脂肪指數(shù)78.46，總平均親水系數(shù)（GRAVY）? 0.352＜0，推測該蛋白為親水性蛋白，無信號肽及跨膜結(jié)構(gòu)域。Blast結(jié)果顯示，該基因編碼的氨基酸序列與隱球菌屬MAPK蛋白激酶相似性最高，為86.41%?；蚬簿幋a426個氨基酸，分子式為C2155H3331N595O623S19，相對分子質(zhì)量48.15 kDa，其中帶負(fù)電荷的殘基（Asp + Glu）52個，帶正電荷的殘基（Arg + Lys）44個，理論等電點（pI）為6.06，不穩(wěn)定性系數(shù)37.31（＜40）為穩(wěn)定蛋白，脂肪指數(shù)82.89，總平均親水系數(shù)（GRAVY）? 0.344＜0，推測該蛋白為親水性蛋白，無跨膜結(jié)構(gòu)域及信號肽?；蚬簿幋a382個氨基酸，分子式為C1953H3035N529O571S17，相對分子質(zhì)量43.61 kDa，其中帶負(fù)電荷的氨基酸殘基（Asp + Glu）52個，帶正電荷的氨基酸殘基（Arg + Lys）43個，理論等電點為（pI）5.82，不穩(wěn)定性系數(shù)38.03（＜40），為穩(wěn)定蛋白，脂肪指數(shù)為84.79，總平均親水系數(shù)（GRAVY）? 0.358＜0，推測該蛋白為親水性蛋白。無信號肽，無跨膜螺旋區(qū)。Blast結(jié)果顯示，該基因編碼的氨基酸序列與黃金銀耳絲裂原活化蛋白激酶相似性最高，一致性為82.4%。

圖1 TfFus3、TfSlt2、Kss1基因ORF序列擴(kuò)增結(jié)果

圖2 TfFus3、TfSlt2、Kss1基因結(jié)構(gòu)

2.3 TfFus3、TfSlt2、Kss1蛋白質(zhì)結(jié)構(gòu)預(yù)測

TfFus3蛋白質(zhì)二級結(jié)構(gòu)預(yù)測結(jié)果如圖3a所示，該蛋白質(zhì)主要由41.41%α-螺旋、13.13%延伸鏈、4.55%β-轉(zhuǎn)角和40.91%無規(guī)則卷曲構(gòu)成，分別由164個、52個、18個和162個氨基酸參與形成。TfFus3蛋白質(zhì)三級結(jié)構(gòu)預(yù)測結(jié)果如圖3 b所示，該蛋白與模板2b9h.1.A一致性達(dá)54.89%，GMQE為0.73（可信度較好），QMEAN為? 2.04（與模板的匹配度較好），該模板蛋白為絲裂源活化蛋白激酶Fus3。

TfSlt2蛋白質(zhì)二級結(jié)構(gòu)預(yù)測結(jié)果如圖4 a所示，其主要由43.9%α-螺旋、11.74%延伸鏈、5.87%β-轉(zhuǎn)角和38.5%無規(guī)則卷曲構(gòu)成，分別由187個、50個、25個和164個氨基酸參與形成。TfSlt2蛋白質(zhì)三級結(jié)構(gòu)預(yù)測結(jié)果如圖4 b所示，該蛋白與模板5z33.1.A一致性達(dá)65.36%，GMQE為0.75（可信度較好），QMEAN為? 1.41（與模板的匹配度較好），該模板蛋白為絲裂源活化蛋白激酶Slt2。

Kss1蛋白質(zhì)二級結(jié)構(gòu)預(yù)測結(jié)果如圖5 a所示，該蛋白質(zhì)主要由46.34%α-螺旋、14.14%延伸鏈、4.97%β-轉(zhuǎn)角和34.55%無規(guī)則卷曲構(gòu)成，分別由177個、54個、19個、132個氨基酸參與形成。Kss1蛋白質(zhì)三級結(jié)構(gòu)預(yù)測結(jié)果如圖5 b所示，該蛋白與模板2b9h.1.A一致性達(dá)54.49%，GMQE為0.71（可信度較好），QMEAN為? 2.41（與模板的匹配度較好），該模板為絲裂原活化蛋白激酶Kss1。

2.4 TfFus3、TfSlt2、Kss1蛋白保守結(jié)構(gòu)域預(yù)測

TfFus3TfSlt2Kss1蛋白保守結(jié)構(gòu)域預(yù)測結(jié)果如圖6、7、8所示。TfFus3蛋白在第44～390位氨基酸間含有一個胞外信號調(diào)節(jié)激酶1和2樣絲氨酸/蘇氨酸激酶的催化結(jié)構(gòu)域（STKc_ERK1_2_ like結(jié)構(gòu)域），屬于PKc蛋白超家族。該結(jié)構(gòu)域主要是將ATP上的γ-磷?；D(zhuǎn)移到蛋白質(zhì)底物上的絲氨酸/蘇氨酸殘基上，該亞家族由絲裂原活化蛋白激酶（MAPK）ERK1、ERK2、釀酒酵母Fus3和類似蛋白質(zhì)組成。TfSlt2蛋白在第21～358位氨基酸間含有一個真菌絲裂原活化蛋白激酶MPK1結(jié)構(gòu)域（STKc_MPK1），屬于PKc_like超家族。同樣催化將γ-磷?；鶑腁TP轉(zhuǎn)移到蛋白質(zhì)底物上的絲氨酸/蘇氨酸殘基上，該亞家族由釀酒酵母的MPK1、粟酒裂殖酵母（）的Pmk1和相似的蛋白質(zhì)組成。Kss1蛋白在第32～375位氨基酸間含有一個細(xì)胞外信號調(diào)節(jié)激酶1和2類絲氨酸/蘇氨酸激酶的催化結(jié)構(gòu)域（STKc_ERK1_2_like），屬于PKc_like超家族。STK催化γ-磷?；鶑腁TP轉(zhuǎn)移到蛋白質(zhì)底物上的絲氨酸/蘇氨酸殘基，這個亞家族由絲裂原活化蛋白激酶（MAPKs）ERK1、ERK2、釀酒酵母Fus3和類似的蛋白質(zhì)組成。

圖6 TfFus3蛋白保守結(jié)構(gòu)域預(yù)測結(jié)果

圖7 TfSlt2蛋白保守結(jié)構(gòu)域預(yù)測結(jié)果

圖8 Kss1蛋白保守結(jié)構(gòu)域預(yù)測結(jié)果

2.5 TfFus3、TfSlt2、Kss1蛋白系統(tǒng)進(jìn)化分析

（1）TfFus3蛋白系統(tǒng)進(jìn)化分析。用ClustalX2及MEGA7軟件將TfFus3蛋白與釀酒酵母中已知的相關(guān)蛋白進(jìn)行同源性比對，結(jié)果顯示該蛋白與酵母中的Fus3蛋白聚為一支（圖9）。

圖9 TfFus3蛋白與釀酒酵母（Saccharomyces cerevisiae）中相關(guān)蛋白序列同源比對結(jié)果

TfFus3蛋白系統(tǒng)進(jìn)化分析如圖10所示，銀耳TfFus3蛋白與、、隱球酵母菌屬的新型隱球酵母（）、蓋替隱球酵母（）、親緣關(guān)系較近。而與角狀膠角耳（）、草菇（、毛韌革菌（）、斑玉蕈（）、灰樹花（）等擔(dān)子菌親緣關(guān)系次之，與白地霉（）、黑曲霉（）釀酒酵母（）等子囊菌親緣性較遠(yuǎn)。

（2）TfSlt2蛋白系統(tǒng)進(jìn)化分析。用ClustalX2及MEGA7軟件對TfSlt2蛋白與釀酒酵母中已知的相關(guān)蛋白進(jìn)行同源性比對，結(jié)果顯示該蛋白與酵母中的Slt2蛋白聚為一支（圖11）。

TfSlt2蛋白系統(tǒng)進(jìn)化分析如圖12所示，銀耳TfSlt2蛋白序列與親緣關(guān)系最近，并與、新型隱球酵母（）、蓋替隱球酵母（）聚為一類。而與立枯絲核菌（）、膠角耳（）、、陶蘭柱擔(dān)菌（）、絨毛栓菌（）、玉米黑粉菌（）等擔(dān)子菌親緣關(guān)系次之，與黑酵母菌（）釀酒酵母()、核盤菌（）等子囊菌親緣性較遠(yuǎn)。

圖10 TfFus3蛋白系統(tǒng)進(jìn)化樹

圖11 TfSlt2蛋白與釀酒酵母（Saccharomyces cerevisiae）中相關(guān)蛋白序列同源比對結(jié)果

圖12 TfSlt2蛋白系統(tǒng)進(jìn)化樹

（3）Kss1蛋白系統(tǒng)進(jìn)化分析。用ClustalX2及MEGA7軟件將銀耳Kss1蛋白與釀酒酵母中已知的相關(guān)蛋白質(zhì)進(jìn)行同源性比對，結(jié)果顯示該蛋白與釀酒酵母中的Kss1蛋白聚為一個分支（圖13）。

Kss1蛋白系統(tǒng)進(jìn)化如圖14所示，銀耳Kss1蛋白與親緣關(guān)系最近，與、、、、、、的親緣關(guān)系次之，與釀酒酵母（）、光滑念珠菌（）親緣關(guān)系較遠(yuǎn)。

3 討論

在真核生物中，異三聚體G蛋白可介導(dǎo)Ca2+、cAMP 和 MAPK 3條信號轉(zhuǎn)導(dǎo)途徑將胞外環(huán)境信號傳遞到細(xì)胞內(nèi)[14]。在MAPK信號通路中主要通過交配信息素途徑介導(dǎo)兩單倍體的交配反應(yīng)[15]。當(dāng)交配信息素與相應(yīng)的受體結(jié)合后激活Gα亞基使GDP轉(zhuǎn)變?yōu)镚TP，此時Gαβγ三聚體解離為Gα和Gβγ復(fù)合體，他們分別通過cAMP、MAPK信號轉(zhuǎn)導(dǎo)途徑將信號傳遞至下游的效應(yīng)分子調(diào)節(jié)相應(yīng)的細(xì)胞活動[16]。G蛋白繼續(xù)通過支架蛋白將信號傳遞至MAPK信號通路中，再經(jīng)MAPKKK-MAPKK-MAPK三級磷酸化過程最終將信號傳遞至TfFus3從而使釀酒酵母的單倍體細(xì)胞發(fā)生交配反應(yīng)，除此之外該通路還與菌絲生長、孢子生成、致病性有關(guān)[17]。

圖13 Kss1蛋白與釀酒酵母（Saccharomyces cerevisiae）中相關(guān)蛋白序列同源比對結(jié)果

圖14 Kss1蛋白系統(tǒng)進(jìn)化樹

銀耳酵母狀芽孢在非適宜條件下，往往通過芽殖、裂殖的方式進(jìn)行無性增殖，當(dāng)其向菌絲態(tài)轉(zhuǎn)變時與白念珠菌侵染宿主的過程類似，即首先會在一端形成芽管，芽管不斷伸長后成為菌絲，在這一過程中細(xì)胞以單一或特定的方向進(jìn)行生長延伸，稱為細(xì)胞的極性生長[18]。酵母中極性生長主要發(fā)生在三個時期：酵母細(xì)胞極性位點的選擇、酵母細(xì)胞向菌絲的轉(zhuǎn)變過程、兩單倍體細(xì)胞的接合生長[19]。菌絲頂端這種高度極性生長的方式依賴于內(nèi)質(zhì)網(wǎng)的功能及細(xì)胞壁的有序合成及重塑[20]。研究發(fā)現(xiàn)，白念珠菌中同源基因敲除后與野生型相比其在42 ℃條件下培養(yǎng)時生長速率降低，且對熱敏感。在高Ca2+濃度、高滲透壓條件下突變體菌落變得光滑，菌絲生長受到抑制[21]。同時Slt2對球孢白僵菌、的毒力、侵襲性生長至關(guān)重要[22-23]。

在釀酒酵母中，有5條分別調(diào)控高滲環(huán)境生長的逆境反應(yīng)、交配反應(yīng)、菌絲生長、細(xì)胞壁完整性和孢子形成的MAPK途徑已被鑒定研究[24-26]。MAPKK是MAPK三級級聯(lián)信號途徑的中間蛋白，其高度保守的SxxxS/T位點受上一級雙特異性MAPKKK的雙重磷酸化誘導(dǎo)激活，其中Fus3和Kss1受同一個MAPKK蛋白的激活，共同調(diào)控交配反應(yīng)和菌絲生長[23]。所以克隆基因有助于研究銀耳二型態(tài)轉(zhuǎn)變的關(guān)鍵基因。分析結(jié)果顯示，銀耳基因與釀酒酵母中的基因具有高度同源性，銀耳Kss1蛋白與親緣關(guān)系最近，本研究結(jié)合銀耳不同形態(tài)下的基因表達(dá)差異分析，發(fā)現(xiàn)基因在菌絲形態(tài)中相較酵母形態(tài)上調(diào)表達(dá)。推測基因可能在銀耳菌絲形成過程中起正調(diào)控作用。

4 結(jié)論

將獲得的銀耳酵母狀芽孢及菌絲體的轉(zhuǎn)錄組數(shù)據(jù)，通過NR注釋、KEGG富集后查找含有完整開放閱讀框的差異基因序列?？寺〔⒎治隽司幋a絲裂原活化蛋白激酶的基因、和。

序列分析發(fā)現(xiàn)，基因組序列全長1 721 bp，ORF全長1 191 bp，含有10個內(nèi)含子，共編碼396個氨基酸，含有絲氨酸/蘇氨酸激酶催化結(jié)構(gòu)域（STKc_ERK1_2_like結(jié)構(gòu)域），與、、新型隱球酵母、蓋替隱球酵母的絲裂原活化蛋白激酶編碼基因親緣關(guān)系最近?；蚪M序列全長1 865 bp，ORF全長1 281 bp，含有12個內(nèi)含子，共編碼426個氨基酸，含有絲裂原活化蛋白激酶MPK1結(jié)構(gòu)域（STKc_MPK1），與親緣關(guān)系最近，并與、新型隱球酵母、蓋替隱球酵母聚為一類?；蚧蚪M序列全長1 903 bp，ORF全長1 149 bp，包含16個內(nèi)含子，共編碼382個氨基酸，含有絲氨酸/蘇氨酸激酶的催化結(jié)構(gòu)域（STKc_ERK1_2_like），與親緣關(guān)系最近。

[1] 彭衛(wèi)紅, 王勇, 黃忠乾, 等. 我國銀耳研究現(xiàn)狀與存在問題[J]. 食用菌學(xué)報, 2005, 12(1) : 51-56.

[2] 孫淑靜. 畫說銀耳優(yōu)質(zhì)高效生產(chǎn)技術(shù)[M]. 北京: 中國農(nóng)業(yè)科學(xué)技術(shù)出版社, 2020: 10-11.

[3] 朱季紅. 群體效應(yīng)因子影響銀耳二型態(tài)轉(zhuǎn)變的研究[D]. 武漢: 華中農(nóng)業(yè)大學(xué), 2015: 7.

[4] 陳明, 陳立國, 汪國蓮. 銀耳極性測定的初步研究[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報, 2000, 19(2): 138-141.

[5] 胡鵬, 張雅芬, 崔海峰, 等. 菰黑粉菌MAPKK同源基因UeMkk1的全長克隆及表達(dá)[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報, 2016, 24(3): 406-415.

[6] LIU H, KOHLER J, FINK G R. Suppression of hyphal formation inby mutation of ahomolog[J]. Science, 1994, 266(5191): 1723-1726.

[7] CHEN R E, THORNER J. Function and regulation in MAPK signaling pathways: lessons learned from the yeast[J]. Biochim Biophys Acta, 2007, 1773(8) : 1311-1340.

[8] 陳江野, 陳曦. 酵母交配過程的極性生長及其調(diào)控機(jī)制[J]. 生命科學(xué), 2002(3): 159-162.

[9] 張科貴, 何光軍. 新生隱球菌細(xì)胞通訊研究進(jìn)展[J]. 微生物學(xué)通報, 2020, 47(1): 272-281.

[10] 陳靜, 王勤, 陳江野., 一個白色念珠菌新的能互補fus/kss突變株的交配缺陷[J]. 生物化學(xué)和生物物理學(xué)報, 2000, 32(3) : 299-304.

[11] NAVARRO-GARCIA F, Alonso-Monge R, Rico H, et al. A role for the MAP kinase genein cell wall construction and morphological transitions in[J]. Microbiology, 1998, 144 (2): 411-424.

[12] HOU Z, XUE C, PENG Y, et al. A mitogen-activated protein kinase gene () inis required for female fertility, heterokaryon formation, and plant infection [J]. Molecular Plant-Microbe Interactions?, 2002, 15(11): 1119-1127.

[13] LEV S, SHARON A, HADAR R,. A mitogen-activated protein kinase of the corn leaf pathogenis involved in conidiation, appressorium formation, and pathogenicity: diverse roles for mitogen-activated protein kinase homologs in foliar pathogens[J]. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96(23): 13542-13547.

[14] NEVES S R, RAM P T, IYENGAR R. G protein pathways [J]. Science, 2002, 296(5573): 1636-1639.

[15] 張海峰. 稻瘟病菌G蛋白及MAPK信號途徑相關(guān)基因的功能分析[D]. 南京: 南京農(nóng)業(yè)大學(xué), 2011: 1-3.

[16] 劉娟, 吳堯, 馬愛民, 等. 真菌二型態(tài)環(huán)境因子及信號轉(zhuǎn)導(dǎo)途徑研究進(jìn)展[J]. 微生物學(xué)通報, 2008, 35(7): 1102-1106.

[17] 彭靜靜. 植物病原真菌中MAPK級聯(lián)通路研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué), 2014, 42(9) : 11-15.

[18] GOW N A R, VAN DE VEERDONK F L, BROWN A J, et al.morphogenesis and host defence: discriminating invasion from colonization[J]. Nature Reviews Microbiology, 2013, 10(2): 112-122.

[19] PRUYNE D, BRETSCHER A. Polarization of cell growth in yeast[J]. Journal of Cell Science, 2000, 113 (Pt 4): 571-585.

[20] 張冰. 白念珠菌Sac1蛋白的鑒定及其功能的研究[D]. 天津: 南開大學(xué), 2015: 79-82.

[21] NAVARRO-GARCIA F, Sánchez M, PLA J. Functional characterization of thegene of, which encodes a mitogen-activated protein kinase homolog related to cell integrity[J]. Molecular and Cellular Biology, 1995, 15(4): 2197-2206.

[22] LUOP X, KEYHANI N O, YU X, et al. The MAP kinase Bbslt2 controls growth, conidiation, cell wall integrity, and virulence in the insect pathogenic fungus[J]. Fungal Genetics and Biology, 2012, 49(7): 544-555.

[23] ONYILO F, TUSIIME G, TRIPATHI J N, et al. Silencing of the mitogen-activated protein kinases (MAPK)andinreduces growth and virulence on host plants[J]. Frontiers in Plant Science, 2018, 9: 291.

[24] ZHAO X, MEHRABI R, XU J R. Mitogen-activated protein kinase pathways and fungal pathogenesis [J]. Eukaryot Cell, 2007, 6(10): 1701-1714.

[25] GUSTIN M C, ALBERTYN J, ALEXANDER M, et al. MAP kinase pathways in the yeast[J]. Microbiology and Molecular Biology Reviews, 1998, 62(4): 1264-1300.

[26] HERSKOWITZ I. MAP kinase pathways in yeast: for mating and more[J]. Cell, 1995, 80(2): 187-197.

Cloning and bioinformatics analysis of key differential-expression genes encoding MAPK indiformis

XU Tao KONG Xuqiang ZHAN Guanping CHEN Yafang CAO Jixuan LI Yaxing SUN Shujing*

（College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China）

is a characteristic edible fungus in China. The diformis transformation is closely related to the actual production and restricts the progress of breeding. Mitogen activated protein kinase (MAPK) pathway plays an important role in the regulation of diformis transformation. Based on transcriptome sequencing data, key differential-expression genes encoding MAPK indiformis were cloned, includinggene,gene andgene. Analysis showed that the total length of genomic DNA sequence of,andgenes was 1 721 bp, 1 865 bp and 1 903 bp, and the total length of ORF sequence was 1 191 bp, 1 281 bp and 1 149 bp, containing 10, 12 and 16 introns, respectively. Coding 396, 426, 382 amino acids, the molecular weight ranging from 43.61 kDa to 48.15 kDa, the theoretical isoelectric point between 5.82 and 6.89. Sequence homology analysis showed that the amino acid sequences encoded by,andwere the most similar to homologous proteins in,and, respectively. Domain analysis showed that TfFus3, TfSlt2 and Kss1 all had typical conserved kinase domains and belonged to PKc_like super family. The results of this study laid a foundation for further exploring the role of MAPK signaling pathway in tremella dimorphism transformation.

; Mitogen activated protein kinase; dimorphic; bioinformatics

S646

A

2095-0934（2024）01-031-09

國家自然科學(xué)基金青年科學(xué)基金項目（31801556）；福建省種業(yè)工程特色食用菌育種與產(chǎn)業(yè)化開發(fā)（zycxny2021011）；福建省食用菌現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)體系

徐濤（2003—），男，在讀本科生，研究方向為食藥用真菌學(xué)。E-mail：2071463989@qq.com。

孫淑靜（1978—），女，博士，教授，主要從事食藥用真菌學(xué)研究。E-mail：shjsun2004@126.com。