杜仲PAL基因cDNA全長序列特征分析

李鐵柱，杜紅巖，劉慧敏

(1.中國林業科學研究院 經濟林研究開發中心，河南 鄭州 450003；2.國家林業局杜仲工程技術研究中心，河南 鄭州 450003；3.中南林業科技大學 林學院，湖南 長沙 410004)

杜仲是我國傳統中藥材[1-3]，杜仲組織和器官中有11種苯丙素類化合物被發現，包括：綠原酸、咖啡酸、香草酸、松柏酸、間羥基苯丙酸、松柏苷、丁香苷、寇布拉苷、愈創木丙三醇、綠原酸甲酯、咖啡酸乙酯[4-9]。

苯丙素類化合物在植物體內生物合成，苯丙氨酸解氨酶(PAL)是這一過程中的第1個關鍵酶，該酶催化L-苯丙氨酸生成反式肉桂酸,再由反式肉桂酸生成芥子酸、阿魏酸、香豆酸等，在其它酶的作用下進一步轉化為香豆素、綠原酸，最終轉化為木質素類、黃酮類和酚類等次生代謝產物[10]。PAL基因一般都含有多個基因家族，PAL亞基又包括2～5個基因的小的多基因家族編碼，隨植物不同而不同[11-15]，如松樹PAL僅有1個PAL基因，菜豆有3個PAL基因：PAL1、PAL2、PAL3，歐芹和煙草的PAL基因有 4 個，番茄則有5個PAL基因，水稻有11個PAL基因，而馬鈴薯的PAL基因就更多，估計達到40～50個之多。PAL多基因家族中，各家族成員的表達活性并不相同，PAL基因是否表達及表達量的相對多少都受植物發育階段和環境信號的調節，表現為嚴格的組織和發育特異性[16-19]。文中對杜仲PAL基因cDNA和基因組全長序列特征進行系統的分析，旨在為杜仲苯丙素類高效合成的分子調控提供依據。

1 材料與方法

分別于2010年5月28日和10月28日，在國家林業局泡桐研究開發中心采集杜仲優良品種“華仲6號”的葉片和果實，提取RNA，送至深圳華大基因公司進行轉錄組測序和基因功能注釋[20-21]。轉錄組測序得到的Unigene，利用NCBI網站上的BLAST工具，和其它物種已知的CDS序列進行在線比較，用相似性得分比較高的已知基因確定并命名杜仲轉錄組中的待分析基因。

利用 ExPASyProtParam(http://cn.expasy.org)程序與ScanProsite程序對EuPAL基因氨基酸殘基數目、氨基酸組成比例、蛋白質相對分子質量、理論等電點以及蛋白質的親水疏水性等性能參數進行預測與分析；利用GOR4程序(http:// bioinf.cs.ucl.ac.uk/psipred)預測EuPAL編碼蛋白質的二級結構；利用Swiss-model(http://swissmodel.expasy.org/)程序，根據基因氨基酸序列，在PBD(Protein data bank，http://www.rcsb.org/pdb/home/home.do)數據庫搜索相似序列，采用alignment mode建模的方式進行同源建模，預測蛋白質的三級結構；用MEGA軟件進行遺傳系統進化分析，構建EuPAL系統發育樹。

2 結果與分析

2.1 EuPAL基因與其它植物種的相似性

EuPAL1(Unigene9151)氨基酸序列與長春 花Catharanthus roseus(BAA95629.1)、煙 草Nicotiana tabacum(BAA22948.1)、 香芹Petroselinum crispum(CAA57057.1)、胡蘿卜Daucus carota(BAC56977.1)、蓖 麻Ricinus communis(XP_002519521.1)等物種的苯丙氨酸脫氨酶氨基酸序列的相似性分別為88％、86％、87％、87％、87％，見圖 1。EuPAL1全長cDNA長度為2 133 bp，編碼 711個氨基酸。

圖1 EuPAL全長氨基酸序列與同源序列的多重比對Fig.1 Multiple alignment of full length amino acid sequence of EuPAL and homologous sequences

另外，EuPAL2(Unigene18458)的氨基酸序列與紫背天葵Gynura bicolor(BAJ17655.1)、葡 萄(XP_002285277.1)、萵 苣Lactuca sativa(AF299330_1)、金光菊Rudbeckia hirta(ABN79671.2)、木薯Manihot esculenta(AAK62030.1)的PAL氨基酸序列的相似性分別為83％、84％、83％、83％、82％，轉錄組數據中未得到全長cDNA序列，但與EuPAL1的氨基酸序列具有較高的相似性。

2.2 EuPAL編碼蛋白質的理化特性

ExPaSy ProtParam軟件預測EuPAL酶蛋白的氨基酸組成中以亮氨酸(11.0％)、丙氨酸(9.4％)、甘氨酸(8.7％)含量為最高，蛋白質相對分子質量約為77.14kD，理論等電點為6.24，蛋白不穩定系數為35.25，屬于穩定蛋白質。

2.3 EuPAL編碼蛋白的親水/疏水性分析

Expasy Protscale程序分析EuPAL編碼蛋白的親水性/疏水性，結果如圖2所示。EuPAL編碼蛋白氨基酸殘基中以A116親水性最強(-1.074)，以A696疏水性最強(1.502)，親水氨基酸殘基的數量的數量小于疏水氨基酸殘基數量，故推斷EuPAL編碼蛋白為疏水性蛋白質。

2.4 EuPAL編碼蛋白的二級結構與保守結構域預測

用GOR4軟件在線預測EuPAL基因編碼蛋白的二級結構，其中α-螺旋占44.30％，β-折疊占8.72％，無規卷曲結構占46.98％，屬于混合型結構的蛋白質(見圖3)。

蛋白質保守結構域分析結果見圖4，EuPAL為多結構域蛋白，屬Lyase-I-Like超家族。

2.5 EuPAL編碼蛋白的三級結構預測

采用Alignment Mode的方式對杜仲PAL基因編碼蛋白進行同源建模。經PBD(Protein data bank，http://www.rcsb.org/pdb/home/home.do) 數據庫搜索，發現與18 458號基因序列一樣，與香菜Petroselinum crispum苯丙氨酸氨脫氨酶晶體結構有較高的同源性，也可能涉及杜仲苯丙烷代謝途徑，在PBD編號為1W27(DOI：10.221 0/pdb1W27/pdb)，以此為模板進行同源性建模，獲得的蛋白模型QMEAN 得分為0.100(見圖5)。

圖2 EuPAL蛋白的親水/疏水性預測Fig.2 Prediction of hydrophobicity and rophilicity of the deduced EuPAL protein

圖3 EuPAL蛋白二級結構預測Fig.3 Prediction of secondary structure of the deduced EuPAL protein

圖4 EuPAL蛋白的保守結構域預測Fig.4 Prediction of conserved domains of the deduced EuPAL protein

2.6 PAL基因系統進化樹

圖5 EuPAL蛋白三級結構預測Fig.5 Prediction of tertiary structure of the deduced EuPAL protein

系統發育樹表明，EuPAL1和EuPAL2基因屬于不同類別(見圖6和圖7)。EuPAL1序列優先與芭蕉、山茶、咖啡等的PAL蛋白聚在一起。推測EuPAL基因很可能還參與到杜仲苯丙素類化合物代謝途徑，從而影響杜仲的生長發育。

圖6 PAL同源蛋白的系統進化樹分析Fig.6 Phylogenetic analysis of PAL homologous proteins

圖7 EuPAL1和EuPAL2系統發育樹Fig.7 EuPAL1 and EuPAL2 phylogenetic tree

3 結論與討論

杜仲苯丙氨酸脫氨酶基因(PAL)是槲皮素合成途徑中的關鍵酶，轉錄組數據中獲得了1條PAL基因全長cDNA，命名為EuPAL1，EuPAL1全長cDNA長度為2 133 bp，編碼 711個氨基酸，相對分子質量為77.14kD，理論等電點為6.24；EuPAL1的氨基酸組成中甘氨酸、丙氨酸、亮氨酸含量較高，屬于穩定的疏水性蛋白，蛋白質二級結構是以α-螺旋和無規卷曲結構為主的混合結構，屬于lyase-I-like蛋白質超家族。

綠原酸是重要的苯丙素類化合物之一，是杜仲重要的藥理功效生物活性成分，受到普遍關注，具有廣泛的生理活性，具有消炎抗菌以及興奮中樞神經等多種作用，對消化、生殖系統疾病均具顯著療效。葡萄糖在酶的催化下轉化成莽草酸，后者再轉化成苯丙氨酸，最后經合成酶的作用得到綠原酸，其中苯丙氨酸脫氨酶基因(PAL)為關鍵酶基因之一。不同植物中PAL基因存在豐富的變異[22-23]，PAL的相對分子質量為55～88kD，分布于細胞質和一些膜細胞器中，如葉綠體、白色體、線粒體、過氧化物體。不同來源的PAL相對分子質量不同，一般介于220～330kD。PAL多基因家族一般均含有2～3個家族成員。PAL基因是一個比較保守的基因，該基因在不同植物種編碼區變化范圍不大，PAL基因編碼區一般在2 100 bp左右，但存在一些差異，同一植物多基因家族中的各個基因編碼長度也不盡相同，PAL在植物次生代謝方面起重要作用，擬南芥、豌豆、松樹楊樹、水稻、煙草、茶樹的PAL基因已被克隆，結果表明，在多數植物中PAL基因以基因家族形式存在，但在不同植物中家族成員數量變化較大，如楊樹至少有2個、菜豆3個、歐芹4個、茶樹1個、馬鈴薯40多個[24]。從杜仲幼果和成熟葉片中鑒定PAL基因，發現至少存在12個家族成員，說明杜仲該基因的遺傳多樣性非常高，超過一般植物，這可能與杜仲除了木質素外，還富含苯丙素、黃酮類等活性成分密切相關。

[1] 杜紅巖,劉昌勇,李 欽,等.杜仲葉中3種主要活性成分含量的季節變化[J].中南林業科技大學學報,2011,31(8):6-9.

[2] 劉淑明,梁宗鎖,董娟娥.不同水分條件下皮葉兩用杜仲林的生長效應[J].中南林業科技大學學報,2007,27(5):49-53.

[3] 杜紅巖,李 欽,杜蘭英,等.杜仲雄花茶營養成分的測定分析[J].中南林業科技大學學報,2007,27(6):88-91.

[4] 薛 萍.杜仲研究現狀與發展前景[J].經濟林研究,1995,13(3): 56-58.

[5] 呂 強,彭密軍,彭 勝,等.不同栽培模式對杜仲葉及枝皮中多種活性成分含量的影響[J].經濟林研究,2012,30(1):73-76.

[6] Deyama T,Kawa T,Nishibe S.The constituent ofEucommia ulmoidesOliv Ⅱ-isolation and structure of three new lignan glycosides [J].Chem Pharm Bull,1985,33(9):3651-3657.

[7] Deyama T,Kawa T,Tagawa S,et al.The constituent ofEucommia ulmoidesOliv Ⅲ-isolation and structure of a new lignan glycosides [J].Chem Pharm Bull,1986,4(2):523-527.

[8] Deyama T,Kawa T,Tagawa S,et al.The constituent ofEucommia ulmoidesOliv Ⅳ-isolation and structure of a new sesquilignan and neoligna glycosides [J].Chem Pharm Bull,1987,35(5): 1803-1807.

[9] Nakamura T,Akazawa Y,Onizuka S,et al.Twelve phenolics from leaves ofEucommia ulmoides[J].Nat Med,1998,52(3):460-465.

[10] 程水源,陳昆松,劉衛紅,等.植物苯丙氨酸解氨酶(PAL)基因的表達調控與研究展望[J].果樹學報,2003,20(5):351-357.

[11] 江昌俊,余有本.苯丙氨酸解氨酶的研究進展(綜述)[J].安徽農業大學學報,2001,(4):425-430.

[12] 程水源,陳昆松,劉衛紅,等.植物苯丙氨酸解氨酶基因的表達調控與研究展望[J].果樹學報,2003,(5):351-357.

[13] 張淑紅,高寶嘉,溫秀軍.棗瘋病過氧化物酶及苯丙氨酸解氨酶的研究[J].植物保護,2004,(5):59-62.

[14] 蘇 艷.苯丙氨酸解氨酶菌株的選育及其全長cDNA序列的克隆[D].南寧:廣西大學,2008.

[15] 崔建東,李 艷,牟德華.苯丙氨酸解氨酶(PAL)的研究進展[J].食品工業科技,2008,(7):306-308.

[16] Hahl brock K,Scheel D.Physiology and molecular biology of phenyl propanoid metabolism [J].Annu Rev Plant Physiol Plant Mol Bio,1989,40(6): 347-369.

[17] Holton TA,Cornish E.Genetics and biochemistry of anthocyanin biosynthesis [J].Plant Cell,1995,7(9): 1071-1083.

[18] Ritter H,Schulz G E.Structural basis for the entrance into the phenyl propanoid metabolism catalyzed by phenylalanine ammonialyase [J].Plant Cell,2004,16(2): 3426-3436.

[19] Butland S L,Chow M L,Ellis B E.Adiverse family of phenylalanine ammonialyase genes expressed in pine trees and cell cultures [J].Plant Mol Bio,1998,37(1): 15-24.

[20] 李鐵柱,杜紅巖,劉慧敏,等.杜仲果實和葉片轉錄組數據組裝及基因功能注釋[J].中南林業科技大學學報,2012,32(11):122-130.

[21] 李鐵柱,杜紅巖,劉慧敏,等.杜仲幼果和成熟果實轉錄組數據組裝及基因功能注釋[J].中南林業科技大學學報,2012,32(10): 9-17.

[22] 徐曉梅,楊署光.苯丙氨酸解氨酶研究進展[J].安徽農業科學,2009,(31):15115-15119,15122.

[23] 宋 婕.丹參苯丙氨酸解氨酶基因(SmPAL1)的克隆及其功能初探[D].陜西師范大學,2007.

[24] 熊 彬.利用苯丙氨酸解氨酶轉化反式肉桂酸生產L-苯丙氨酸的研究[D].南寧:廣西大學,2005.