陳恩發(fā),劉 輝,丁海兵,潘 牧,王啟富
(貴州省農業(yè)科學院 生物技術研究所,貴州 貴陽 550006)
【研究意義】研究魔芋軟腐病病原菌的致病機理,對進一步探明該病原菌的侵染過程和該病的防治提供理論依據?!厩叭搜芯窟M展】魔芋(Konjac)為天南星科魔芋屬多年生草本植物,在我國南部地區(qū)廣泛種植。魔芋含有豐富的葡甘聚糖,因其具有保水性、成膜性、膠凝性、乳化性、穩(wěn)定性等良好的特性,被廣泛地應用于食品、農業(yè)、環(huán)保、醫(yī)藥、化工、生物等領域[1-2]。近年來,隨著魔芋種植面積的不斷擴大,各種病害便隨之發(fā)生,且日漸嚴重,嚴重制約了魔芋產業(yè)的發(fā)展。魔芋軟腐病是危害魔芋生產的主要病害之一[3],發(fā)病率高,死亡率大且極易傳染,常大面積、連片爆發(fā),給魔芋產業(yè)帶來毀滅性打擊,極度增加了農民的種植風險[4]。劉佩瑛等[5]研究發(fā)現,魔芋軟腐病是由胡蘿卜軟腐歐氏桿菌胡蘿卜軟腐亞種(Eriwiniacarotovorapv.carotovora)侵染所致。修建華等[6]將軟腐病菌株分為3個組群,分別為菊歐文氏桿菌(Erwiniachrysanthemi,E.ch.)、胡蘿卜軟腐歐文氏桿菌胡蘿卜軟腐亞種(Erwiniacarotovorasubsp.carotovora,E.c.c.)和一組尚不能確定的菌株。魔芋感染軟腐病菌后,葉片或塊莖變軟、倒伏、發(fā)黑腐爛,并帶有刺激的酒臭味[7]。胡蘿卜軟腐歐文氏菌致病因子主要依賴其分泌的植物細胞壁降解酶,包括果膠酶、纖維素酶、以及多種蛋白水解酶等[8]。其中,果膠酶是致病過程中的主要致病因子,果膠酶是一種以多種形式存在的同功酶,由果膠酸鹽裂解酶(pel)、果膠甲基酯酶(pme)、果膠質裂解酶(pnl)、聚半乳糖醛酸酶等組成,有不同的基因編碼,并成簇排列。果膠酸鹽裂解酶(pel)是果膠酶中最重要的酶,由5個主要果膠裂解酶和至少3個次要果膠裂解酶組成[9]。在病原菌侵染的不同時期均有果膠酶基因的表達[10],有研究表明,果膠酶基因通過與一些基因的協同作用,可使病原菌穿透植物組織的能力提高,加速植物細胞壁的降解,使細胞壁成分被釋放出來為病原菌提供營養(yǎng)[11]?!颈狙芯壳腥朦c】國內對于魔芋軟腐病致病機理的研究報道較少,未見魔芋軟腐病果膠酸鹽裂解酶相關研究報道?!緮M解決的關鍵問題】通過克隆技術得到魔芋軟腐病中2個重要的果膠鹽裂解酶基因pelD和pelE,并對2個基因進行生物信息學、原核表達和酶活性分析,以期獲得正常表達的pelD和pelE基因。
1.1.1 菌株 研究采用的魔芋軟腐病病原菌菌株由課題組分離自貴州省臺江縣。
1.1.2 試劑EscherichiacoliDH5α、pET28a(+)和EscherichiacoliBL21(DE3)均購自寶生物工程(大連)有限公司;RNA提取試劑盒、DNA提取試劑盒、T4DNA連接酶、2×TaqPCR MasterMix、TaqDNA Polymerase、6×buffer、10×buffer、SYBR RT-PCT Kit、蛋白酶K、DLMarker 2000、DLMarker 15000、PMD-18T和限制性內切酶均購于TaKaRa生物公司;cDNA第一鏈合成試劑盒、SM0431 maker購于Fermentas公司;質粒小提試劑盒、瓊脂糖凝膠DNA回收試劑盒、PCR產物純化試劑盒均購自Tiangen公司;卡那霉素、氨芐青霉素、潮霉素B(Sangon Biotech)購于sigma公司;瓊脂糖購自Biowest公司;其他試劑均為進口或國產分析純。
1.1.3 培養(yǎng)基 NA 培養(yǎng)基:酵母浸膏 1 g、牛肉浸膏 3 g、葡萄糖 10 g、蛋白胨 5 g、瓊脂 18 g,加蒸餾水定容至1000 mL,調pH至7.0。NA 液體培養(yǎng)基:酵母浸膏1 g、葡萄糖 10 g、牛肉浸膏 3 g、蛋白胨5 g,加蒸餾水定容至 1000 mL,調節(jié)pH至7.0。LB 培養(yǎng)基:酵母粉5 g、胰蛋白胨10 g、NaCl 10 g,加適量蒸餾水緩慢加熱使其充分溶解,用蒸餾水定容至 1000 mL,加熱混勻。
1.2.1 基因組 DNA 的提取 提取參照DNA提取試劑盒說明書進行,-20 ℃保存?zhèn)溆谩?/p>
1.2.2pelD和pelE基因引物設計 從GenBank中下載已注冊的pelD和pelE基因同源序列,利用Primer 5 設計2個基因的擴增引物(表1)。
表1 pelD和pelE基因擴增引物Table 1 The amplified primers of pelD and pelE genes
1.2.3 PCR 擴增及產物克隆 以提取的DNA為模板,設計的引物進行PCR 擴增,參數:94 ℃預變性3 min;94 ℃變性30 s,64 ℃(pelD)、60 ℃(pelE)退火1 min,72 ℃延伸10 min,30個循環(huán);72 ℃再延伸10 min,12 ℃保存。電泳檢測(1 %瓊脂糖凝膠),膠回收目的片段。將目的片段與 PMD-19T載體的連接,反應體系10 μl:PMD-19TVector DNA 1 μl,PCR 目的片段產物 4 μl,Solution I 5 μl,16 ℃連接過夜。連接產物轉化至大腸桿菌E.coliDH5α感受態(tài)細胞,37 ℃培養(yǎng)過夜,挑選單菌落,搖床培養(yǎng)擴繁,參考Tiangen公司的質粒小提試劑盒說明書提取質粒,雙酶切驗證,送出測序。
1.2.4pelD和pelE基因序列生物信息學分析 用DNAMAN (Virsion 6.0)軟件和Protparam軟件預測分析目的基因片段序列??缒そY構預測在網頁http://www.cbs.dtu.dk/services/TMHMM/上進行;亞細胞定位分析在網頁http://www.cbs.dtu.dk/services/TargetP/上進行;二級結構分析在網頁http://npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html上進行,三級結構分析在網頁SWISS-MODLE (http://swissmodel.expasy.org/interactive)上進行;利用Doolittle 算法對目的蛋白進行親水/疏水性分析。
1.2.5 目的蛋白的原核表達分析 分別用相應的酶雙酶切連有目的基因片段(pelD和pelE基因)的T-easy質粒和pET28a載體,膠回收目的基因和載體片段,將目的基因片段和載體片段按1∶1混合,加入T4連接酶,室溫連接2~3 h,轉化至大腸桿菌E.coliDH5α感受態(tài)細胞中,37 ℃培養(yǎng)過夜,挑取單菌落進行搖菌培養(yǎng)、菌液PCR及酶切驗證,將重組質粒轉入大腸桿菌BL21(DE3)受體菌,37 ℃ 培養(yǎng)過夜,挑取單菌落,置于LB液體培養(yǎng)基中37 ℃過夜培養(yǎng)。取上清液接種于新鮮的LB 液體培養(yǎng)基中(菌液上清液:LB 液體培養(yǎng)基=1∶100)培養(yǎng),在A600下測吸光值在0.6~0.8時,加入IPTG調節(jié)終濃度為1 mmol/L,于37 ℃誘導反應5 h,分別取1 mL菌液,經超聲波破碎和離心后,將菌體沉淀用0.2 mmol/L PBS(Na2HPO4-Na2HPO4)重懸,室溫搖床培養(yǎng)3 h, 4 ℃、12 000 r/min離心30 min,12 %的SDS-PAGE凝膠電泳檢測上清液,用鎳瓊脂糖親和層析法洗脫。
1.2.6 酶活性分析 超聲波處理1.2.5中誘導5 h及未誘導的菌體(功率80 %,超聲波振蕩2 s,間歇4 s,總時間為10 min,制冷),15 000 r/min離心5 min,取上清液用于測定酶活,酶活測定參照李洋等[11]的方法進行。
從圖1看出,對pelD基因和pelE基因進行常規(guī) PCR 擴增,pelD基因在 1209 bp處有一特異條帶,pelE基因在1080 bp 處有一特異條帶,均與預測條帶大小一致。
將pelD和pelE基因的目的片段經過凝膠回收,克隆至 PMD19-T 載體,轉化E.coliDH5α感受態(tài)細胞。挑取陽性菌落進行 PCR 擴增(圖2)和提取質粒雙酶切(圖3),均得到與預期條帶大小相當的條帶。
測序結果表明,pelD基因序列大小為1209 bp,pelE基因序列大小為1080 bp,與目的基因片段大小一致。Blast比對分析發(fā)現,這2個基因的序列與已報道的果膠酶基因pelD和pelE具有較高的同源性,且有完整的閱讀框。
通過Protparam預測分析發(fā)現(圖4),pelD基因編碼一個完整的開放閱讀框,閱讀框架全長為1209 bp,編碼氨基酸殘基402個,具有起始密碼子ATG和終止密碼子TAA?;蚓幋a的蛋白分子式為 C1891H2938N536O588S6,分子量為42.78 kDa,等電點為5.79,半衰期為30 h。不穩(wěn)定系數為22.12,為穩(wěn)定性蛋白。含39個負電荷氨基酸(Asp + Glu),35個正電荷氨基酸(Arg + Lys),沒有賴氨酸拉鏈,果膠酸鹽裂解酶折疊和致病因子的蛋白保守區(qū)在第94~301個氨基酸區(qū)域,在1~16 個氨基酸處存在信號肽;pelE基因含有一個完整的開放閱讀框,閱讀框架全長為1080 bp,編碼氨基酸殘基359個,含有起始密碼子ATG 和終止密碼子TAA。基因編碼的蛋白分子式為C1891H2938N536O588S6,分子量為37.9 kDa,等電點為9.43,半衰期是30 h,不穩(wěn)定系數為31.50,為穩(wěn)定性蛋白。帶負電荷氨基酸 (Asp + Glu)24個,正電荷氨基酸(Arg + Lys)41個,不含有賴氨酸拉鏈,果膠酸鹽裂解酶折疊和致病因子的蛋白保守區(qū)為第94~301個氨基酸區(qū)域,在1~21個氨基酸處存在信號肽。
經跨膜結構分析發(fā)現,pelD蛋白無跨膜區(qū),pelE蛋白在50~72個氨基酸處存在一跨膜結構域。經Protscale分析 pelD蛋白的親水/疏水性發(fā)現,該蛋白含有165個疏水氨基酸,237個親水氨基酸,親水性氨基酸58.96 %,且在整個肽鏈中均勻分布,推測pelD蛋白為親水性蛋白。pelE蛋白含有144個疏水氨基酸,215個親水氨基酸,親水性氨基酸59.89 %,且在整個肽鏈中均勻分布,推測pelE蛋白也為親水性蛋白。對pelD蛋白質進行亞細胞定位分析得出,pelD蛋白定位于線粒體的可能性為4.4 %,作為信號肽的可能性為92.2 %。pelE蛋白作為信號肽的可能性為16.8.2 %,為其他成分的可能行為87.9 %。對pelD蛋白質二級結構預測發(fā)現,該蛋白含α螺旋24.88 %,β轉角9.20 %,無規(guī)則卷曲39.05 %,延伸鏈26.87 %。pelE蛋白質含α螺旋17.55 %,β轉角11.98 %,39.83 %無規(guī)則卷曲,延伸鏈30.64 %。并用 SWISS-MODLE預測得pelD 蛋白質的三級結構和pelE 蛋白質的三級結構(圖5),pelD蛋白三級結構折疊方式為單鏈右旋β螺旋,解析度為1.60A,扭曲值為-2.88,溶解值為-2.74,Cβ值為-1.51,Oligo-state為單體,GMQE值為0.64,QMEAN值為-3.68。pelE蛋白三級結構折疊方式為單鏈右旋β螺旋,解析度為2.10A,扭曲值為-2.91,溶解值為-1.56,Cβ值為-0.75。Oligo-state為單體,GMQE值為0.42,QMEAN值為-3.53。
對轉化重組表達載體 pET28a-pelD和pET28a-pelE菌株轉化至感受態(tài)大腸桿菌BL21(DE3),經IPTG誘導表達5 h后,SDS-PAGE 電泳分析看出(圖6),pelD 和 pelE蛋白均在在約45 ku處出現明顯的條帶。
從圖7看出,誘導后pelD和pelE的酶活性均高于誘導前,pelD的酶活是誘導前的1.79倍,pelE是誘導前的1.75倍。表明,pelD及pelE基因均能正常表達,且經誘導后活性均顯著上升,即研究獲得了具有正常功能的魔芋pelD和pelE基因。
果膠鹽裂解酶是是植物細胞壁降解酶中重要的酶類,是植物病原菌的主要致病因子,其有關克隆和分析報道較多,1995年Heikinheimo等[12]研究Ecc SCC3193中的果膠酶發(fā)現,其編碼的蛋白含有信號肽且與Ecc71菌株中的pel-3編碼的蛋白相似,但與其他細菌分泌的pel蛋白差異較大。筆者等研究獲得的pelD和pelE基因與李琦等[13]的研究結果相同,均含有信號肽。Hla等[14]研究表明,不同菌株中的同一種果膠酶同源性較低,可能是由細菌本身的特殊性導致。研究得到的pelD和pelE基因與其他果膠酶基因同源性較高,但這一結果是否證明這2個基因在進化上保守有待進一步研究。
果膠酶基因是軟腐病侵染植物組織過程中的重要致病基因,作為果膠酶基因家族中的重要成員,pelD和pelE基因主要在降解植物組織細胞壁果膠物質中發(fā)揮作用。有研究表明,作為獨立的轉錄單元,pelA、pelE和pelD基因對外界環(huán)境條件比較敏感,pelD和pelE基因的表達受pH 的影響較大,在堿性條件下,二者的活力均較低[15];在酸性條件下,pelD基因主要表達于病菌侵染過程的早期,而pelE基因則在正常情況下轉錄[16]。在cAMP途徑中,pelD和pelE基因的轉錄表達受CRP蛋白的抑制,但低濃度的Fe2+有利于pelD和pelE基因的表達。本研究對魔芋軟腐病菌的pelD和pelE基因進行克隆分析和酶活性分析獲得pelD和pelE基因的完整序列。經原核表達分析pelD和pelE的表達產物均在約45 ku 處有明顯的條帶,與Dickeya[17]的研究結果較吻合,IPTG誘導后pelD和pelE基因的酶活性分別是誘導前的1.79和1.75倍。表明,研究獲得了具有正常功能的魔芋軟腐病的pelD和pelE基因。為下一步研究魔芋軟腐病菌果膠酶基因的功能提供了參考,為魔芋軟腐病軟腐病的防治提供了重要依據。
通過PCR克隆技術獲得具有完整閱讀框的魔芋軟腐病pelD和pelE基因,2個基因均在約45 ku處出現明顯的條帶,pelD和pelE誘導后酶活分別是誘導前的1.79和1.75倍,2個基因具有正常功能。