李麗秀　賴鐘雄

摘 要 以早鐘6號(hào)枇杷試管苗葉片為試驗(yàn)材料，采用同源克隆方法從葉片mRNA中分離得到CAT基因的2個(gè)成員Ej-CAT1和Ej-CAT2，Ej-CAT1 cDNA序列全長(zhǎng)1 738 bp，包括1 479 bp ORF和257 bp 3′-UTR；Ej-CAT2序列全長(zhǎng)1 742 bp，包含1 479 bp ORF及261 bp 3′-UTR。2個(gè)CAT基因成員的核苷酸序列相似性為98.04%，其開(kāi)放閱讀框推導(dǎo)的氨基酸序列具有97.32%的相似性。生物信息學(xué)分析結(jié)果表明：2個(gè)蛋白均編碼492個(gè)氨基酸，是親水蛋白，不含信號(hào)肽，不存在螺旋卷曲結(jié)構(gòu)，是非跨膜蛋白，亞細(xì)胞定位于過(guò)氧化物酶體；Ej-CAT1及Ej-CAT2蛋白分別有18、17個(gè)功能位點(diǎn)，22、19個(gè)磷酸化位點(diǎn)。CAT基因在枇杷試管苗離體保存過(guò)程中不同時(shí)期具有不同的表達(dá)量，整體趨勢(shì)呈“W”形狀，在離體保存1個(gè)月時(shí)表達(dá)量最高，保存7個(gè)月時(shí)表達(dá)量最低。CAT表達(dá)量變化可能跟枇杷試管苗的生長(zhǎng)發(fā)育模式有關(guān)。

關(guān)鍵詞 枇杷；試管苗；CAT；基因克隆；表達(dá)分析

中圖分類號(hào) S667.3 文獻(xiàn)標(biāo)識(shí)碼 A

過(guò)氧化物酶（Hydrogen Peroxidase）又稱觸酶（CAT，Catalase），廣泛存在于動(dòng)物、植物及微生物中。過(guò)氧化物酶是最早發(fā)現(xiàn)與植物種子活力有關(guān)的氧化酶之一[1]，能參與活性氧代謝過(guò)程，催化H2O2分解為H2O和O2，防止細(xì)胞過(guò)氧化；能催化分解甲酸鹽、亞硝酸及過(guò)氧化乙醇[2-6]。CAT能與SOD、POD組成活性氧防御系統(tǒng)，在植物防御、延緩衰老等方面起著重要作用，能夠控制細(xì)胞內(nèi)的氧化還原平衡，與植物的抗旱、脅迫應(yīng)答（耐鹽脅迫、耐溫度脅迫）、抗病及減輕重金屬等有關(guān)[7-8]。此外，CAT在食品工業(yè)、紡織、造紙、環(huán)保、醫(yī)藥等行業(yè)中得到廣泛應(yīng)用[7-12]。

許多研究表明，CAT基因是多基因家族，不同成員在植物中的分布及功能有差異，具有組織差異性表達(dá)特點(diǎn)。1997年，Willekens等[12]從煙草中克隆得到3種CAT同工酶，并報(bào)道了cat1主要是清除光呼吸產(chǎn)生的H2O2，cat2可能特異的清除氧化脅迫產(chǎn)生的H2O2，而cat3主要是清除乙醛酸循環(huán)體中脂肪酸產(chǎn)生的H2O2。Shikanai[13]及Miyagawa[14]等報(bào)道了煙草葉綠體中的CAT表達(dá)增強(qiáng)了對(duì)光、干旱等氧化脅迫的抗性。John G等[15]報(bào)道了將CATs轉(zhuǎn)化其他物種，能提高植物對(duì)某些逆境的適應(yīng)性。程華等[16]從銀杏中分離得到GbCAT1，其在葉、根、莖、果中均有表達(dá)，且表達(dá)量受滲透壓、ABA、紫外、溫度脅迫等影響。梁國(guó)慶等[17]研究發(fā)現(xiàn)高鈣處理12 h和低鈣處理24 h，蘋果CAT1基因表達(dá)量顯著增加，且不完全受鈣調(diào)控。王鳳德等[18]報(bào)道了豌豆CAT在煙草中的過(guò)量表達(dá)，提高了煙草的抗旱性及抗氧化性。齊興柱等[19]從香蕉中分離得到2個(gè)CAT基因，并發(fā)現(xiàn)其在氧化脅迫環(huán)境中起主要作用。林星谷等[20]從冬棗果實(shí)中分離得到ZjCAT基因，該基因在冬棗果實(shí)成熟衰老中具調(diào)控作用。楊芳[21]研究發(fā)現(xiàn)豌豆CAT轉(zhuǎn)化玉米提高了玉米的抗旱性。余迪求[22]發(fā)現(xiàn)過(guò)氧化氫酶在轉(zhuǎn)基因馬鈴薯中的過(guò)量表達(dá)，能提高馬鈴薯的抗病性。在南瓜中得到3個(gè)CAT基因成員，擬南芥也有3個(gè)CAT基因成員，小麥至少有3個(gè)CAT成員等。

目前CAT基因的研究領(lǐng)域涉及植物、動(dòng)物及微生物，花椰菜、銀杏、蘋果、香蕉、冬棗、金絲桃、玉米、巴西橡膠、豌豆、小麥、番茄、甘蔗、蘿卜、南瓜、辣椒、水稻等植物[11，16-35]CATs基因均已被克隆出來(lái)，并做了相應(yīng)的分析研究，但在枇杷上還未見(jiàn)報(bào)道。為了進(jìn)一步研究CAT在枇杷試管苗離體保存以及其它生長(zhǎng)發(fā)育過(guò)程中的作用，本試驗(yàn)以枇杷試管苗為材料，采用同源克隆方法，從中克隆得到CAT基因的2個(gè)成員Ej-CAT1、Ej-CAT2的全長(zhǎng)序列，并對(duì)其進(jìn)行序列分析和結(jié)構(gòu)預(yù)測(cè)；同時(shí)，對(duì)該基因在離體保存過(guò)程中6個(gè)不同時(shí)期的表達(dá)進(jìn)行定量分析，為進(jìn)一步研究其功能和調(diào)控途徑奠定基礎(chǔ)。

1 材料與方法

1.1 材料

早鐘6號(hào)枇杷試管苗，由福建農(nóng)林大學(xué)園藝植物生物工程研究所保存。選取生長(zhǎng)旺盛的試管苗的葉片作為試驗(yàn)材料。試管苗離體保存材料接種在1/3 MS+PP333的培養(yǎng)基上，取材時(shí)間分別為保存1、4、5、7、9及11個(gè)月。

1.2 方法

1.2.1 引物的設(shè)計(jì)及合成 根據(jù)Gnebank上已登錄的其他物種CAT cDNA序列，設(shè)計(jì)保守區(qū)引物T1F、T1R；應(yīng)用得到的保守區(qū)序列，設(shè)計(jì)3′端、5′端擴(kuò)增引物P2F、P3F和P4R、P5R；根據(jù)拼接的cDNA全長(zhǎng)序列，在其起始密碼子和終止密碼子處分別設(shè)計(jì)引物ORF-F、ORF-R，驗(yàn)證開(kāi)放閱讀框。定量表達(dá)分析的引物在CAT ORF內(nèi)設(shè)計(jì)，設(shè)計(jì)的引物委托北京六合華大基因科技股份有限公司合成。各引物序列及其PCR擴(kuò)增時(shí)所用參數(shù)見(jiàn)表1所示。

1.2.2 目的片段的擴(kuò)增 以合成的cDNA第一鏈為模板，進(jìn)行保守區(qū)擴(kuò)增及3′ RACE；以cDNA第二鏈為模板，進(jìn)行5′端擴(kuò)增。PCR反應(yīng)體系：cDNA模板1μL，2×Taq PCR MasterMix 12.5 μL，上游引物1 μL，下游引物1 μL，加ddH2O 9.5 μL至總體積25 μL。反應(yīng)程序如下：94 ℃預(yù)變性3 min，94 ℃變性1 min，按表1所示溫度退火1 min，72 ℃延伸t min[表1所示時(shí)間（t）]，從變性到延伸共35個(gè)循環(huán)，最后再72 ℃延伸10 min結(jié)束反應(yīng)，

1.2.3 目的片段的回收、連接轉(zhuǎn)化及序列測(cè)定

用DNA回收試劑盒回收目的片段，將目的片段與PMD-18T載體連接并轉(zhuǎn)化DH5α感受態(tài)細(xì)胞，轉(zhuǎn)化產(chǎn)物涂在含有1 μL Amp ∶ 1 mL LB的平板中培養(yǎng)，12～16 h后挑起單菌落進(jìn)行重組子菌液鑒定，陽(yáng)性克隆子送到北京六合華大基因科技股份有限公司測(cè)序。

1.2.4 生物信息學(xué)分析 用GenBank上的Blast檢索已有序列與分析其保守結(jié)構(gòu)域，采用DNAMAN 6.0和primer 5.0軟件進(jìn)行引物設(shè)計(jì)、序列拼接及翻譯氨基酸。用ExPASy（http：//www.expasy.ch/toos.html）數(shù)據(jù)庫(kù)中的ProtParam工具和ProtScale工具預(yù)測(cè)分析推測(cè)蛋白質(zhì)的理化性質(zhì)，蛋白質(zhì)跨膜結(jié)構(gòu)預(yù)測(cè)與分析采用TMPRED（http：//www.ch.embnet.org/software/TMPRED_form.html），磷酸化位點(diǎn)預(yù)測(cè)用NetPhos2.0（http：//www.cbs.dtu.dk/services/NetPhos），利用MEGA 5.0構(gòu)建系統(tǒng)進(jìn)化樹。

1.2.5 熒光定量分析 將枇杷試管苗接種在1/3 MS+4.0 mg/L PP333+30 g/L白糖的培養(yǎng)基上，分別取保存時(shí)間為1、4、5、7、9、11個(gè)月的材料提取總RNA，按TaKaRa公司的PrimeScriptR RT reagent Kit說(shuō)明書逆轉(zhuǎn)錄成cDNA。熒光定量以Actin基因?yàn)閮?nèi)參對(duì)照，在枇杷CAT基因最保守的地方設(shè)計(jì)引物，引物序列見(jiàn)表2。采用兩步法qRT-PCR，反應(yīng)體系為SYBRR Primer Ex TaqTMⅡ（2×）10.0 μL，PCR Forward Primer 0.8 μL，PCR Reverse Primer 0.8 μL，cDNA模板1 μL，加ddH2O 7.4 μL。反應(yīng)程序?yàn)椋?5 ℃ 5 min；95 ℃ 10 s，58 ℃ 10 s，72 ℃ 10 s，共40個(gè)循環(huán)；95 ℃ 5 s，60 ℃ 1 min；40 ℃ 30 s。每個(gè)反應(yīng)重復(fù)3次，以6個(gè)時(shí)期的cDNA混樣稀釋5倍梯度制作標(biāo)準(zhǔn)曲線，并得到各個(gè)Ct值，最終計(jì)算CAT基因的相對(duì)表達(dá)量。

2 結(jié)果與分析

2.1 枇杷CAT基因保守區(qū)的克隆

以枇杷保守區(qū)cDNA為模板，用引物T1F、T1R進(jìn)行保守區(qū)PCR擴(kuò)增，得到2條731 bp的相似序列（圖1-A），與預(yù)期結(jié)果相符合，將其命名為CAT1、CAT2。經(jīng)NCBI比對(duì)，發(fā)現(xiàn)CAT1、CAT2核酸序列與碧桃、甜櫻桃、葡萄、棗、楊樹、蓖麻、荔枝、龍眼等CAT的相似性均在80%以上，與碧桃的相似性在90%以上；CAT1、CAT2推導(dǎo)氨基酸序列與桃、棗、楊樹、李、葡萄、金絲桃等相似性在85%以上，故推斷CAT1、CAT2均是枇杷CAT基因的保守區(qū)序列。

經(jīng)DNAMAN 6.0比對(duì)CAT1、CAT2核酸及氨基酸序列，結(jié)果見(jiàn)圖2。從圖2可知，兩核酸序列相似性為96.58%，在中間部位有25個(gè)不同堿基，氨基酸相似性為95.47%，含有11個(gè)不同氨基酸。

2.2 枇杷CAT基因的3′、5′末端的獲得與ORF分析

2.2.1 枇杷CAT基因的3′、5′末端的獲得 以保守區(qū)cDNA第一鏈為模板，將引物P2F、P3F與通用引物3P、3NP配對(duì)進(jìn)行巢式PCR擴(kuò)增。第二輪PCR擴(kuò)增產(chǎn)物經(jīng)膠回收、連接轉(zhuǎn)化及測(cè)序得到2條序列，其擴(kuò)增產(chǎn)物電泳圖見(jiàn)圖1-B，序列長(zhǎng)度分別為578、582 bp。經(jīng)DNAMAN6.0 比對(duì)，發(fā)現(xiàn)得到的2條序列均與保守區(qū)序列有重疊部分，能夠拼接起來(lái)，故推斷2條序列均為枇杷CAT的3′序列。以5′cDNA為模板，用引物P4F、P4R多次克隆枇杷CAT的5′RACE，均未能得到序列，重新比對(duì)GenBank中其它物種的CAT基因序列后發(fā)現(xiàn)，CAT基因在起始密碼子處高度同源，故在CAT基因的起始密碼子處設(shè)計(jì)1條順式上游引物T5F，以5′cDNA為模板，T4F、T5F為引物進(jìn)行PCR擴(kuò)增反應(yīng)，擴(kuò)增得到1條750 bp左右的目的條帶（圖1-C）。測(cè)序得到1條722 bp大小的序列，與保守區(qū)序列拼接有一段高度重疊部分，說(shuō)明得到的序列是CAT基因的5′端序列。

2.2.2 枇杷CAT基因ORF分析 將所得的5′端序列、保守區(qū)序列與3′端序列拼接起來(lái)，得到2條全長(zhǎng)序列，分別命名為Ej-CAT1、Ej-CAT2，GenBank登錄號(hào)分別為：JX307086、JX307085。Ej-CAT1 cDNA序列全長(zhǎng)1 738 bp，包括1 479 bp ORF及257 bp 3′-UTR，Poly（A）尾長(zhǎng)19 bp；Ej-CAT2 cDNA序列全長(zhǎng)1 742 bp，包括1 479 bp ORF及261 bp 3′-UTR，Poly（A）尾長(zhǎng)23 bp。在Ej-CAT1、Ej-CAT2序列的起始密碼子和終止密碼子處各設(shè)1條引物，進(jìn)行PCR擴(kuò)增，擴(kuò)增得到2條1 500 bp左右的條帶（圖1-D），經(jīng)測(cè)序得到2條1 479 bp的相似序列，分別與Ej-CAT1、Ej-CAT2的開(kāi)放閱讀框序列一致，進(jìn)一步驗(yàn)證了Ej-CAT1、Ej-CAT2序列的準(zhǔn)確性。

2個(gè)CAT基因成員的核苷酸序列相似性為98.04%，其開(kāi)放閱讀框推導(dǎo)的氨基酸序列具有97.32%的相似性。NCBI比對(duì)發(fā)現(xiàn)，2個(gè)CAT ORF序列與碧桃、甜櫻桃、棗、葡萄、毛果楊、橡膠樹、龍眼、棉花等其他物種的覆蓋度多在99%以上，相似性在80%以上。將CAT1、CAT2翻譯的氨基酸在線Blast，發(fā)現(xiàn)與碧桃CAT氨基酸的同源性為89%，與葡萄、棗、毛果楊、擬南芥的相似性分別為89%、88%、87%、83%，與其他物種的CAT氨基酸的相似性也在80%以上，屬于Catalase-like超級(jí)家族，由此判斷2個(gè)序列均為枇杷CAT的cDNA序列。

2.3 枇杷CAT的生物信息學(xué)分析

2.3.1 蛋白質(zhì)的理化性質(zhì)預(yù)測(cè)與分析 由預(yù)測(cè)結(jié)果可知Ej-CAT1蛋白翻譯492個(gè)氨基酸，分子總量為56 991.2，等電點(diǎn)pI為6.66，帶正電荷（Arg+Lys）與帶負(fù)電荷（Asp+Glu）氨基酸殘基總數(shù)分別為57、61。該蛋白含有20種氨基酸，其中Phe、Pro、Val 3種氨基酸的含量最多，均占7.1%；Cys、Met的含量最少，僅占全部氨基酸的1.2%。其分子式為C2 586H3 869N715O727S1，原子總數(shù)為7 909，不穩(wěn)定指標(biāo)為40.46，說(shuō)明該蛋白是一個(gè)不穩(wěn)定蛋白。其脂肪系數(shù)為71.71，總親水性值（GRAVY）為-0.520，預(yù)測(cè)該蛋白是親水蛋白。

Ej-CAT2蛋白分子總量為57 087.3，等電點(diǎn)為6.65，帶負(fù)電荷與帶正電荷氨基酸殘基總數(shù)分別為61、57。該蛋白的分子式為C2 594H3 869N715O727S12，原子總數(shù)為7 917，不穩(wěn)定系數(shù)為41.57，是不穩(wěn)定蛋白。其脂肪系數(shù)為70.33，總親水性值（GRAVY）為-0.539，說(shuō)明該蛋白是一個(gè)親水性蛋白，而且親水性比Ej-CAT1蛋白的親水性強(qiáng)。該蛋白翻譯的492個(gè)氨基酸中共含有20種氨基酸，其中Phe、Pro兩種氨基酸的含量最多，占總氨基酸數(shù)的7.3%；Cys的含量最少，僅占1.0%，而且其含量低于Ej-CAT1蛋白中的Cys含量。

2.3.2 蛋白質(zhì)親水性/疏水性預(yù)測(cè) ProtScale分析Ej-CAT1、Ej-CAT2蛋白的親疏水性，預(yù)測(cè)結(jié)果顯示Ej-CAT1蛋白親疏水性值低于-2.000的氨基酸有36個(gè)，最低值為-2.722，位于第375氨基酸處；分值大于+2.000的氨基酸有4個(gè)，最大值為第326處氨基酸的+2.667，親水性氨基酸明顯大于疏水性氨基酸，預(yù)測(cè)該蛋白為親水性蛋白，與ProtParam工具的預(yù)測(cè)結(jié)果一致。

Ej-CAT2蛋白的預(yù)測(cè)結(jié)果顯示，親疏水性值低于-2.000的氨基酸有33個(gè)，最低值是第427氨基酸處，分值為-2.678；分值大于+2.000的氨基酸只有2個(gè)，分別位于氨基酸103與186處，其中186處氨基酸的分值最大，分值為+2.233，其親水性氨基酸多于疏水性氨基酸，預(yù)測(cè)結(jié)果為親水性蛋白。

2.3.3 蛋白質(zhì)的磷酸化位點(diǎn)預(yù)測(cè) 經(jīng)NetPhos2.0 對(duì)枇杷Ej-CAT1蛋白、Ej-CAT2蛋白的磷酸化位點(diǎn)進(jìn)行預(yù)測(cè)，結(jié)果見(jiàn)圖3。從圖3-A可知，Ej-CAT1蛋白的Ser、Thr和Tyr都存在磷酸化現(xiàn)象，共具有22個(gè)磷酸化位點(diǎn)，其中10個(gè)Ser位點(diǎn)，分別位于氨基酸10、14、70、85、104、112、164、347、437、491處；6個(gè)Thr位點(diǎn)，分別位于氨基酸97、115、273、292、351、414處；6個(gè)Tyr位點(diǎn)，分別位于氨基酸4、221、256、348、392、453處。從圖3-B中可知，Ej-CAT2蛋白共有19個(gè)磷酸化位點(diǎn)，其中9個(gè)Ser位點(diǎn)，分別位于第10、14、70、85、104、112、164、437、491個(gè)氨基酸；5個(gè)Thr位點(diǎn)，分別位于第97、115、273、292、351個(gè)氨基酸；5個(gè)Tyr位點(diǎn)，分別位于第4、221、256、392、453個(gè)氨基酸。推測(cè)2個(gè)CAT蛋白均是以Ser為主，Thr、Tyr為輔的方式。

與Ej-CAT1蛋白的相比，Ej-CAT2蛋白少了3個(gè)磷酸化位點(diǎn)，分別為第347氨基酸處的Ser位點(diǎn)，第414氨基酸處的Thr位點(diǎn)與第348氨基酸處的Tyr位點(diǎn)。

2.3.4 蛋白質(zhì)功能位點(diǎn)預(yù)測(cè) 用ScaanProsite在線分析枇杷Ej-CAT1蛋白、Ej-CAT2蛋白的功能位點(diǎn)，結(jié)果顯示：Ej-CAT1蛋白含有18個(gè)功能位點(diǎn)，其中3個(gè)N端糖基化位點(diǎn)，蛋白激酶C磷酸化位點(diǎn)4個(gè)，酪蛋白激酶Ⅱ磷酸化位點(diǎn)有5個(gè)，酪氨酸磷酸化位點(diǎn)1個(gè)（386-392aa），N端豆蔻酰位點(diǎn)4個(gè)，酰胺化位點(diǎn)1個(gè)（414-417aa）。1個(gè)Catalase proximal heme-ligand singnature，1個(gè)Catalase proximal active site singnature。Ej-CAT2蛋白含有17個(gè)功能位點(diǎn)，其中3個(gè)N端糖基化位點(diǎn)，4個(gè)蛋白激酶C磷酸化位點(diǎn)，4個(gè)酪蛋白激酶Ⅱ磷酸化位點(diǎn)，1個(gè)酪氨酸磷酸化位點(diǎn)，4個(gè)N端豆蔻酰位點(diǎn)，1個(gè)酰胺化位點(diǎn)。1個(gè)Catalase proximal heme-ligand singnature，1個(gè)Catalase proximal active site singnature。

與枇杷Ej-CAT1蛋白功能位點(diǎn)相比，Ej-CAT2蛋白少了1個(gè)酪蛋白激酶Ⅱ磷酸化位點(diǎn)（347-350 aa），Ej-CAT2蛋白386-392 aa處酪氨酸磷酸化位點(diǎn)為Rde.Evd.Y，332-337 aa處N-豆蔻酰化位點(diǎn)為GvyySN，與Ej-CAT1蛋白有所不同。

2.3.5 蛋白質(zhì)亞細(xì)胞定位預(yù)測(cè) 用PSORT對(duì)Ej-CAT1、Ej-CAT2蛋白亞細(xì)胞定位進(jìn)行預(yù)測(cè)，發(fā)現(xiàn)2個(gè)CAT蛋白定位于過(guò)氧化物酶體的可能性較大，分值為0.748；其次是線粒體基質(zhì)，分值為0.100。而Ej-CAT1、Ej-CAT2 蛋白信號(hào)肽及跨膜結(jié)構(gòu)預(yù)測(cè)表明，2個(gè)蛋白均是非分泌蛋白，不含跨膜結(jié)構(gòu)，是非跨膜蛋白。

2.3.6 蛋白質(zhì)結(jié)構(gòu)預(yù)測(cè) PSIPRED軟件預(yù)測(cè)表明Ej-CAT1、Ej-CAT2蛋白的二級(jí)結(jié)構(gòu)是由α-螺旋、β-折疊和無(wú)規(guī)則卷曲3種結(jié)構(gòu)組成，其中無(wú)規(guī)則卷曲所占的比例較大，α-螺旋、β-折疊所占比例相當(dāng)。枇杷CAT兩個(gè)蛋白的二級(jí)結(jié)構(gòu)并沒(méi)有顯著差異，只有兩個(gè)細(xì)微差異。Ej-CAT2蛋白在氨基酸349-355處、475處形成α-螺旋，而Ej-CAT1蛋白則沒(méi)有形成α-螺旋。GOR4軟件預(yù)測(cè)結(jié)果與PSIPRED預(yù)測(cè)結(jié)果類似，Ej-CAT1二級(jí)結(jié)構(gòu)中無(wú)規(guī)則卷曲占59.55%，α-螺旋占18.50%，β-折疊占21.95%；枇杷Ej-CAT2二級(jí)結(jié)構(gòu)中無(wú)規(guī)則卷曲占59.76%，α-螺旋占16.46%，β-折疊占23.78%。

利用SWISS-MODEL工具對(duì)Ej-CAT1、Ej-CAT2蛋白的三級(jí)結(jié)構(gòu)進(jìn)行在線預(yù)測(cè)，預(yù)測(cè)結(jié)果如圖4所示。兩個(gè)CAT蛋白的三級(jí)結(jié)構(gòu)基本相似，主要是由無(wú)規(guī)則卷曲，α-螺旋，β-折疊組成。

2.3.7 蛋白質(zhì)系統(tǒng)進(jìn)化樹的構(gòu)建 經(jīng)NCBI BLASTn比對(duì)后，發(fā)現(xiàn)2個(gè)枇杷CAT氨基酸與碧桃（Prunus persica）、棗（Ziziphus jujuba）、毛果楊（Populus trichocarpa）、葡萄（Vitis vinifera）、甜櫻桃（Prunus avium）、貫葉金絲桃（Hypericum performatum）、陸地棉（Gossypium hirsutum）、胡麻（Sesamum indicum）、煙草（Nicotiana benthamiana）、荷花（Nelumbo nucifera）、木欖（Bruguiera gymnorhiza）、冰葉日中花（Mesem-bryanthemum crystallinum）、甜椒（Capsicum annuum）、木薯（Manihot esculenta）、擬南芥（Arabidopsis thaliana）、馬蹄蓮（Zantedeschia aethiopica）、蘿卜（Raphanus sativus）、蕓薹（Brassica juncea）、水稻（Oryza sativa）、玉米（Zea mays）的相似性分別均在80%以上。將枇杷兩個(gè)CAT成員與這22個(gè)氨基酸進(jìn)行分子進(jìn)化樹的構(gòu)建，結(jié)果見(jiàn)圖5。從圖5可見(jiàn)，整個(gè)進(jìn)化樹分成2個(gè)組群，水稻屬于小組群，Ej-CAT1、Ej-CAT2與其余的植物均屬于大組群。Ej-CAT1、Ej-CAT2與葡萄、碧桃、甜櫻桃、木薯、小麥同處一小分支，有著很近的進(jìn)化關(guān)系；與甜椒、荷花、馬蹄蓮、玉米有著相對(duì)較遠(yuǎn)的進(jìn)化關(guān)系，與水稻cat-c的進(jìn)化距離最遠(yuǎn)；與其他植物有著相對(duì)較近的進(jìn)化關(guān)系。

2.4 枇杷試管苗離體保存過(guò)程中CAT基因的相對(duì)定量表達(dá)分析

由Real Time PCR反應(yīng)得到的CAT及內(nèi)參基因Actin基因的溶解曲線，均只出現(xiàn)1個(gè)單峰。從擴(kuò)增曲線圖上看，兩個(gè)基因的曲線拐點(diǎn)清楚，指數(shù)期明顯，整體平行性較好，表明CAT及Actin基因的引物特異性較好。制作的標(biāo)準(zhǔn)曲線結(jié)果顯示，枇杷CAT和內(nèi)參基因Actin的標(biāo)準(zhǔn)曲線線性范圍較廣，起始模板稀釋倍數(shù)相關(guān)系數(shù)與Cq值良好。CAT及Actin標(biāo)準(zhǔn)曲線的擴(kuò)增效率分別為1.933、1.958，標(biāo)準(zhǔn)曲線的斜率均在-3.4～-3.6 之間。綜上，這兩個(gè)基因的引物特異性良好，符合試驗(yàn)要求，可進(jìn)行下一步試驗(yàn)。

6個(gè)不同時(shí)期材料進(jìn)行qPCR，得到所有Ct值，計(jì)算各時(shí)期的相對(duì)表達(dá)量，結(jié)果見(jiàn)圖6。CAT基因在6個(gè)不同保存時(shí)期都有表達(dá)，表達(dá)量的整體趨勢(shì)呈字母“W”形狀。保存1個(gè)月時(shí)，表達(dá)量最高；到保存4個(gè)月時(shí)，表達(dá)量急劇下降，其表達(dá)量?jī)H是保存1個(gè)月時(shí)表達(dá)量的1/3；保存5個(gè)月時(shí)表達(dá)量稍微升高，到保存7個(gè)月時(shí)表達(dá)量又開(kāi)始下降；保存9個(gè)月時(shí)，該基因表達(dá)量迅速上升，其表達(dá)量是保存7個(gè)月表達(dá)量的2倍；而保存9個(gè)月、11個(gè)月的表達(dá)量沒(méi)有明顯變化，表達(dá)相對(duì)穩(wěn)定。

3 討論與結(jié)論

3.1 枇杷CAT蛋白功能多樣性

枇杷葉片兩個(gè)CAT成員的氨基酸序列與多種植物的CAT蛋白氨基酸序列的相似性均在80%以上。在CAT氨基酸進(jìn)化樹分析中，枇杷不同CAT蛋白與葡萄、碧桃、甜櫻桃、木薯、聚類在一起，與水稻cat-c的進(jìn)化距離最遠(yuǎn)。根據(jù)生物學(xué)信息分析結(jié)果，枇杷葉片不同CAT成員均不含跨膜結(jié)構(gòu)，不含信號(hào)肽，沒(méi)有卷曲螺旋結(jié)構(gòu)，是不穩(wěn)定的親水蛋白。亞細(xì)胞定位于過(guò)氧化物酶體，能與SOD、POD組成活性氧防御系統(tǒng)，參與活性氧代謝過(guò)程，催化H2O2分解為H2O和O2，能控制細(xì)胞內(nèi)的氧化還原平衡，在枇杷試管苗離體保存中的植物防御、延緩衰老等方面起著重要作用。潛在磷酸化預(yù)測(cè)結(jié)果顯示，CAT1蛋白具有22個(gè)磷酸化位點(diǎn)，CAT2蛋白有19個(gè)磷酸化位點(diǎn)，2個(gè)CAT成員均是以絲氨酸為主，酪氨酸、蘇氨酸為輔的方式，調(diào)控枇杷試管苗離體保存中的信號(hào)轉(zhuǎn)導(dǎo)及基因表達(dá)等活動(dòng)。蛋白質(zhì)功能的預(yù)測(cè)結(jié)果顯示：CAT1和CAT2分別有18、17個(gè)功能位點(diǎn)，其中分別具有3個(gè)N-糖基化位點(diǎn)，4個(gè)蛋白激酶C磷酸化位點(diǎn)，4個(gè)N-豆蔻酰位點(diǎn)，1個(gè)酰胺化位點(diǎn)，5或4個(gè)酪蛋白激酶Ⅱ磷酸化位點(diǎn)，CAT1蛋白比CAT2蛋白多了1個(gè)酪蛋白激酶Ⅱ磷酸化位點(diǎn)。這些均說(shuō)明了枇杷CAT含有多種蛋白質(zhì)翻譯后修飾方式[36]，其在枇杷體內(nèi)擁有不同的功能作用。

枇杷2個(gè)CAT成員的氨基酸序列具有97.36%的相似性，擁有部分相似的功能和結(jié)構(gòu)，但2個(gè)CAT蛋白還存在不同的結(jié)構(gòu)與功能，因此，枇杷CAT兩個(gè)成員可能在枇杷試管苗離體保存中相互協(xié)調(diào)、各自分工的方式來(lái)發(fā)揮作用。

3.2 枇杷CAT基因在試管苗離體保存中的作用

CAT對(duì)H2O2具有親和力，能夠在各種脅迫下清除多余的H2O2，防止形成羥自由基。CAT是多基因家族，其表達(dá)受到發(fā)育和環(huán)境等因素影響[33]，且不同基因成員的表達(dá)具有時(shí)空差異。程華等[16]認(rèn)為銀杏CAT基因的表達(dá)受到ABA、滲透壓及溫度等環(huán)境因素的影響。

枇杷試管苗離體保存過(guò)程中，各個(gè)時(shí)期CAT基因的表達(dá)量不同，在培養(yǎng)基中加入植物生長(zhǎng)抑制劑PP333，保存1個(gè)月時(shí)CAT基因表達(dá)量最高，之后保存4個(gè)月、5個(gè)月及7個(gè)月的表達(dá)量均有所回落，在保存9個(gè)月和11個(gè)月時(shí)表達(dá)量又上升到一定水平。高表達(dá)量的CAT能夠有效減少試管苗體內(nèi)的活性氧積累，維持植物的生理代謝，能夠延緩衰老。試管苗保存1個(gè)月時(shí)，枇杷試管苗葉片慢慢發(fā)育成熟，此時(shí)會(huì)有大量的活性氧產(chǎn)生，高表達(dá)CAT能夠減少活性氧積累，降低氧化傷害。在之后的幾個(gè)月，CAT基因的表達(dá)量減少，可能跟試管苗側(cè)芽抽出新的嫩葉有關(guān)。到保存9個(gè)月、11個(gè)月時(shí)，CAT基因表達(dá)量維持在一定程度，可能是試管苗大部分枝葉進(jìn)入成熟和衰老階段，培養(yǎng)基中PP333也能夠有效的控制植物的衰老，一大部分減少了植物體內(nèi)一些活性氧的積累。枇杷CAT基因在整個(gè)離體保存中表達(dá)量不斷變化，推測(cè)其在試管苗離體保存中起著重要的作用。

參考文獻(xiàn)

[1] 張坤生， 田薈琳. 過(guò)氧化氫酶的功能研究[J]. 食品科技， 2007（1）： 8-11.

[2] 宋新華， 趙鳳云. 植物體內(nèi)過(guò)氧化氫酶的研究進(jìn)展[J]. 安微農(nóng)業(yè)科學(xué)， 2007， 35（31）： 9 824-9 827.

[3] 劉 冰， 梁嬋娟. 生物過(guò)氧化氫酶研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào)， 2005， 21（5）： 223-225.

[4] 南芝潤(rùn)， 范月仙. 植物過(guò)氧化氫酶的研究進(jìn)展[J]. 安徽農(nóng)學(xué)通報(bào)， 2008， 14（5）： 27-29.

[5] 楊冬業(yè)， 張麗珍， 徐淑慶. 西瓜不同生長(zhǎng)階段的過(guò)氧化氫酶活性研究[J]. 安徽農(nóng)業(yè)科學(xué)， 2012， 40（18）： 9 604-9 606.

[6] 劉靈芝， 鐘廣蓉， 熊 蓮， 等. 過(guò)氧化氫酶的研究與應(yīng)用進(jìn)展[J]. 化學(xué)與生物工程， 2009， 26（3）： 15-18.

[7] 陳金鋒， 王宮南， 程素滿. 過(guò)氧化氫酶在植物脅迫響應(yīng)中的功能研究進(jìn)展[J]. 西北植物學(xué)報(bào)， 2008， 28（1）： 188-193.

[8] 宋新華， 趙鳳云. 植物體內(nèi)過(guò)氧化氫酶的研究進(jìn)展[J]. 安微農(nóng)業(yè)科學(xué)， 2007， 35（31）： 9 824-9 827.

[9] 徐 剛. 過(guò)氧化氫酶在染紗生產(chǎn)中的應(yīng)用[J]. 紡織科技進(jìn)展， 2005 （2）： 29-30.

[10] 張永科， 趙景嬋， 郭治安， 等. 過(guò)氧化氫酶催化乙烯表面改性的研究[J]. 西北大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2007， 37（4）： 595-597.

[11] 李 蓉. 花椰菜Fe-SOD、 CAT、 Tub基因和青花菜Tub基因cDNA的克隆及分析[D]. 福州： 福建農(nóng)林大學(xué)， 2011.

[12] Willskens H， Chamnongpol S， Davey M， et al. Catalase is a sink for H2O2 and is indispensable for stress defence in C3 plants[J]. EMBOJ， 1997， 16（16）： 4 806-4 816.

[13] Shikanai T， Takeda T， Yamauchi H， et al. Inhibition of ascorbate peroxidase under oxidative syress in tobacco having bacterial catalase in chlorop[J]. Lasts FEBS Lett， 1998， 428： 47-51.

[14] Miyagawa Y， Tamoi M， Shigeoka. Evaluation of the defence system in chlorop lasts to photooxidative stress caused by paraquat using transgenic tobacco plants expressing catalase from Escherichiacoil[J]. Plant Cell Physiol， 2000， 41： 311-320.

[15] John G， Scandalios， Alberto Acev， et al. Catalase gene expression in responseto chronichigh temperature snbss in maize[J]. Plant Science， 2000， 156： 103-110.

[16] 程 華， 李琳玲， 許 鋒， 王 燕， 等. 銀杏過(guò)氧化氫酶基因（CAT1）的克隆及表達(dá)分析[J]. 林業(yè)科學(xué)研究， 2010， 23（4）： 493-499.

[17] 梁國(guó)慶， 孫靜文， 周 衛(wèi)， 等. 鈣對(duì)蘋果果實(shí)超氧化物歧化酶、 過(guò)氧化氫酶活性及其基因表達(dá)的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào)， 2011， 17（2）： 438-444.

[18] 王鳳德， 衣艷君， 王海慶， 等. 豌豆過(guò)氧化氫酶在煙草葉綠體中的過(guò)量表達(dá)提高了植物的抗逆性[J]. 生態(tài)學(xué)報(bào)， 2011， 31（4）： 1 058-1 063.

[19] 齊興柱， 楊臘英， 黃俊生. FOC4的2個(gè)過(guò)氧化氫酶基因的克隆及表達(dá)分析及其引起的香蕉苗活性氧迸發(fā)研究[J]. 中國(guó)農(nóng)學(xué)通報(bào)， 2012， 28（15）： 163-169.

[20] 林星谷， 孔德倉(cāng)， 龐曉明， 等. 冬棗過(guò)氧化氫酶基因的克隆及表達(dá)分析[J]. 西北植物學(xué)報(bào)， 2012， 32（6）： 1 086-1 092.

[21] 楊 芳. 豌豆過(guò)氧化氫酶基因在玉米中的轉(zhuǎn)化[D]. 濟(jì)南： 山東師范大學(xué)， 2006.

[22] 余迪求， 岑 川， 楊明蘭， 等. 玉米不同組織過(guò)氧化氫酶水楊酸敏感性的差異和外源水楊酸處理提高玉米抗病性的研究[J]. 植物學(xué)報(bào)， 1999， 41（12）： 1 293-1 298.

[23] 李夜雨， 楊 平. 貫葉金絲桃過(guò)氧化氫酶基因的生物信息學(xué)分析[J]. 安徽農(nóng)學(xué)通報(bào)， 2009， 15（13）： 34-35.

[24] Liu Hanmei， Zhang Huaiyu， Tan Zhenbo， et al. Molecular Cloning and Chilling Expression of Maize Catalase-3 Gene[J]. Journal of Sichuan Agriculture University， 2006， 24（3）： 272-275.

[25] 余乃通， 孫玉娟， 張雨良， 等. 海南島巴西橡膠和香蕉cat相關(guān)基因的克隆及分析[J]. 植物研究， 2012， 32（1）： 47-53.

[26] Garratt L， Cjanagoudar B， Lowe K， et al. Salinity tolerance and antioxidant status in cotton cultures[J]. Free Radic Biol Med， 2002， 33（4）： 502-511.

[27] 佘茂云. 小麥組織培養(yǎng)再生相關(guān)基因克隆及分子特性分析[D]. 北京： 中國(guó)農(nóng)業(yè)科學(xué)院， 2011.

[28] 佘茂云. 小麥TaCAT新基因克隆及分子生物學(xué)和生化特性分析[J]. 科技導(dǎo)報(bào)， 2011， 29（12）： 23-31.

[29] Polidoros A， Mylona P， Scandal J. Transgenic tobacco plants expressing the maize Cat2 gene have altered catalase levels that affect plant pathogen interactions and resistance to oxidative stress[J]. Transgenic Res， 2001， 10（6）： 555-569.

[30] 陳姍姍， 郭晉隆， 李國(guó)印， 等. 甘蔗過(guò)氧化氫酶基因的電子克隆及生物信息學(xué)分析[J]. 生物信息學(xué)， 2012， 10（1）： 65-70.

[31] Yang W L， Liu J M， Chen F， et al. Identification of Festuca arundinacea Schreb Cat1 gene and analysis of its expression under abiotic stresses[J]. Journal Integrative Plant Biol， 2006， 48（3）： 334-343.

[32] Kwon S I， Lee H， An C S. Differential expression of three catalase genes in the small radish（Rhaphanus sativus L. var.sativus）[J]. Mol Cells， 2007， 24（1）： 37-44.

[33] Esaka M， Yamada N， Kitabayashi M， et al. cDNA cloing and differential gene expression of three catalase in pumpkin[J]. Plant molecular biology， 1997， 33（1）： 141-155.

[34] Lee S H， An c s. Differential expression of three catalase genes in hot pepper（Capsicum annuum L）[J]. Molecules and Cells， 2005， 20（2）： 247-255.

[35] Chen Z， Iyer S， Caplan A， et al. Differential accumulation of salicylic acid and salicylic acid-sensitive catalase in different rice tissues[J]. Plant Physiology， 1997， 114（1）： 193-201.

[36] Taouatas N， Altelaar A F， Drugan M M， et al. Strong cation exchange-based fractionation of Lys-N-generated peptidws facilitates the targeted analysis of post-translational modifications[J]. Molecular & Cellular Proteomics， 2009， 8（1）： 190-200.