巴西橡膠樹HbCYP450基因克隆與表達分析

李曉娜 肖厚貞 萬三連 張冬 張宇

摘 要 細胞色素P450氧化酶催化植物初生代謝和次生代謝等多種類型的氧化還原反應。本研究在橡膠樹葉片中克隆得到一個CYP450基因，編碼516個氨基酸，在3～511位氨基酸序列間存在一個P450超級保守結構域，且存在特征性的血紅素結合位點保守序列，與蓖麻、胡楊、毛楊果和大豆CYP450的序列一致性分別為74%、68%、68%和65%。HbCYP450基因在橡膠樹的葉片、花和膠乳中均表達，但在樹皮中不表達。該基因的表達量受干旱、機械傷害和白粉病的誘導，在植物激素吲哚乙酸（IAA）、茉莉酸（JA）、乙烯利（ET）、脫落酸（ABA）、水楊酸（SA）和過氧化氫（H2O2）作用下，HbCYP450基因表達量呈現先顯著上升后顯著下降的規律。說明HbCYP450是植物P450家族成員，其表達量受白粉菌和逆境誘導。

關鍵詞 巴西橡膠樹；基因克??；表達分析；CYP450

中圖分類號 S432.1 文獻標識碼 A

Abstract Cytochrome P450 monooxygenases catalyze a wide varieties of monooxygenation reactions in primary and secondary metabolism in plants. In this study， a CYP450 gene was cloned from rubber tree leaves. It encodes 516 amino acids which has one conserved domain at position 3-511 and a special conserved sequence of heme binding sites. HbCYP450 has identity of 74%， 68%， 68%， and 65% with Ricinus communis L.， Populus euphratica， Populus trichocarpa and Glycine max， respectively. HbCYP450 expressed in leaf， flower and latex， but no expressed in bark， which expression of gene induced by drought， wounding and powdery mildew disease.Under the treatment of IAA， ABA， SA， JA， ET and H2O2， the expression patterns of HbCYP450 were significantly upregulated， then downregulated dramatically. These results suggest that HbCYP450 is one of the members of P450 family， and the expression are induced by powdery mildew and stress induction.

Key words Hevea brasiliensis； gene cloning； expression analysis； CYP450

doi 10.3969/j.issn.1000-2561.2017.11.017

細胞色素P450作為植物中最大的酶蛋白家族，參與生物合成途徑和生物解毒途徑[1-2]。具有生物合成功能的CYP450在木質素中間物、植物激素（脫落酸、赤霉素、生長激素等）、淄醇、萜類、黃酮類等生物物質的合成中起重要作用[3-7]，具有生物解毒功能的CYP450能夠使污染物、農藥等有毒物質變成非毒性物質[8]。CYP450廣泛存在于動植物、細菌以及真菌等細胞體內，目前從植物中獲得的CYP450普遍存在于擬南芥[9]、番茄、小麥[10]、水稻、蓖麻[11]、玉米、大豆等植物中。目前，擬南芥中有40個CYP450基因功能已經確定，水稻中有39個CYP450基因功能已經確定[12]。根據CYP450蛋白序列的相似性和親緣關系將CYP450歸類為多個家族和亞家族，其中，同一家族成員的相似性大于40%，同一亞族成員的相似性大于55%[13]，同一亞家族的基因具有類似的功能[14]。

橡膠樹是我國重要的產膠植物，其生長發育經常受到各種生物脅迫和非生物脅迫，其中白粉病對橡膠樹的嫩葉、嫩芽和花絮均能造成危害，使植株葉片掉落，嚴重時會導致植株死亡。干旱和機械損傷等逆境因素會導致巴西橡膠樹病蟲害發生率和死皮率增加，引起巨大的經濟損失[15-16]。研究發現，木質素含量會影響植株抵抗脅迫和水分運輸能力[17-19]，并且高含量木質素有利于植物抗病蟲害。楊家書等[20]研究發現小麥的抗白粉病能力與木質素含量成正相關。橡膠樹膠乳的合成途徑中涉及諸多氧化還原反應[21]，每一個氧化酶反應都需要一個專一的CYP450酶[22]，CYP450對膠乳的生產可能具有一定的影響，但目前對橡膠樹中CYP450的研究甚少。筆者通過克隆橡膠樹中CYP450基因，分析其在不同處理下的表達情況，為探究其在橡膠樹生長發育和抗逆脅迫中的作用打下良好基礎。

1 材料與方法

1.1 材料

以中國熱帶農業科學院種植的橡膠樹品種CATAS 7-33-97正常割膠的成齡樹為材料，進行基因的克隆與不同組織的表達分析。選取穩定期生長均勻一致的橡膠樹CATAS 7-33-97芽接苗進行機械損傷、白粉菌生理小種 HO-73接種、脫落酸（ABA）、乙烯利（ETH）、茉莉酸（JA）和過氧化氫（H2O2）等處理[23]。

1.2 方法

1.2.1 樣品的采集和cDNA的合成 機械傷害處理：在芽接苗的嫩葉上用解剖刀割出橫切面，再將葉片置于水中，模擬機械損傷[24]；白粉菌處理采用Wang等[25]的方法；植物激素及過氧化氫處理方法參考Qin等[26]的方法：用乙醇溶解各種激素，噴于橡膠樹芽接苗上。在處理后0、0.5、2、6、10、24、48和72 h采集樣品。將不同處理下不同時間段采集到的葉片和不同組織放置在-80 ℃冰箱中冷凍保存。

提取橡膠樹樹皮、膠乳、葉片和花的總RNA，具體操作方法參考何鑫等[27]的方法。采用超微量核酸蛋白分析儀進行總RNA濃度和純度的檢測，并利用1%甲醛變性膠電泳進行RNA完整性的檢測[28]。cDNA的合成采用RNA反轉錄試劑盒（Takara，中國大連）完成。

1.2.2 橡膠樹HbCYP450基因克隆 根據Genbank上公布人類和一些模式植物的過氧化物酶蛋白序列，在EST數據庫中進行tblastn搜索，并通過軟件DNAman拼接同源片段，得到橡膠樹HbCYP450基因的cDNA序列[29]。利用軟件primer 6.0設計了基因特異引物HbCYP450-F（5′-TCGGTTGCTCATTG

TCTT-3′）和HbCYP450-R （5′-CAGTGGGCTACTAATACTATAC-3′），將反轉錄的cDNA第Ⅰ鏈作為模板，對橡膠樹細胞色素CYP450基因的cDNA序列進行擴增。

1.2.3 橡膠樹HbCYP450基因結構分析 使用NCBI ORF Finder預測基因的編碼區序列以及分析巴西橡膠樹HbCYP450的理化性質，通過數據庫NCBI Conserved Domains、在線分析軟件SMART和TMHMM Server v. 2.0分析其結構域， 用SignalP 4.1 Server預測氨基酸序列，再通過ExPASy ProtParam分析信號肽結構，亞細胞定位使用PSORT Prediction進行分析。最后在NCBI蛋白數據庫中利用blastp搜索其他植物的CYP450同源蛋白并獲取其序列，利用軟件DNAMAN 6.0進行序列比對，利用MAGE6.0分析細胞色素P450家族系統的進化關系。

1.2.4 橡膠樹HbCYP450基因表達分析 采用軟件primer 6.0設計HbCYP450基因熒光定量表達分析的引物HbCYP450F（5′-ATTCTATGGCTCCGCTA

TAC-3′）和HbCYP450R（5′-TGATTGCTGCTGCTGATT-3′）。以橡膠樹HbACTIN（F： GATGTGGATATCAGGAAGGA和R： CATACTGCTTGGAGCAAGA）為內參，進行熒光定量PCR反應[30]。

1.3 統計分析和作圖

采用單因素方差分析和Duncan檢驗分析，數據統計分析采用SPSS 21.0，使用Origin 2015科技繪圖軟件制圖。

2 結果與分析

2.1 橡膠樹HbCYP450基因克隆與生物信息學分析

利用模式植物擬南芥、水稻和小麥的CYP450蛋白序列在橡膠樹EST數據庫中做tblastn搜索后，克隆得到全長為1 610 bp的cDNA序列，編碼區長度為1 548 bp，編碼516個氨基酸，分子量為59.08 ku，等電點PI為8.59，分子式為C2 687H4 233N717O735S23，總原子數為8395，其編碼形成的是一種親水且不穩定的蛋白。在3～511位氨基酸處存在一個P450超級保守結構域（圖1劃線部分），在5～27位氨基酸處存在一個跨膜螺旋，通過細胞定位在線分析軟件可知，HbCYP450在線粒體、質膜、細胞核、內質網膜和內質網腔中定位的可能性分別是0.328 0、0.460、0.126、0.1和0.1。將同源序列進行比對分析發現，橡膠樹HbCYP450與蓖麻RcCYP450（XP_002524040.1）、胡楊PeCYP450（XP_011005628.1）、毛楊果PtCYP450（XP_006380437.1）、大豆GmCYP450（XP_003556213.1）的序列一致性分別為74%、68%、68%、65%。在HbCYP450蛋白上存在一些CYP450特有的保守結構，如C-螺旋序列WXXXR，在HbCYP450中位于132～136位氨基酸，氧結合結構域AGXXT，在HbCYP450中位于314～318位氨基酸，K-螺旋序列EXXR和PXRF，在HbCYP450分別位于372～375位和426～429位氨基酸，血紅素結合位點保守結構序列：FXXGXRXCXG，在HbCYP450中位于450～459位氨基酸，序列為：FGAGRRMCPG。將橡膠樹HbCYP450與擬南芥CYP450進行系統進化分析發現，擬南芥CYP89A2與HbCYP450親緣關系最近（圖2）。

2.2 橡膠樹HbCYP450基因的表達分析

HbCYP450基因在橡膠樹的花、葉和膠乳中均有表達，且在膠乳中的表達量是葉片中的65倍，但其在樹皮中不表達。白粉菌侵染橡膠樹葉片時，HbCYP450基因表達量在侵染初期便顯著上調，隨著病害程度的加重，表達量顯著增加，說明白粉菌能夠誘導HbCYP450基因正向表達。在機械傷害橡膠樹初期，HbCYP450基因的表達量顯著上調，處理后24 h表達量達到最高值（圖3）。說明HbCYP450基因參與橡膠樹生物和非生物脅迫響應過程。

在植物激素誘導下，HbCYP450基因的表達情況如圖4所示。在吲哚乙酸（IAA）的誘導下，HbCYP450基因的表達量在處理后6 h顯著上調，到24 h表達量又顯著下調。在茉莉酸（JA）、乙烯利（ET）、脫落酸（ABA）、水楊酸（SA）和過氧化氫（H2O2）分別處理下，橡膠樹葉片中HbCYP450基因表達量的變化趨勢與IAA誘導下的變化相同，其中IAA與ET、ABA、JA三者都是在處理后6 h表達量達到最高值，分別是處理前的6、25、11和20倍。H2O2和SA分別是在24 h和10 h達到表達最高值，分別是處理前的18和8倍。說明植物激素能夠有規律的誘導HbCYP450基因表達，同時也說明該基因參與橡膠樹植物激素處理響應。

3 討論

P450s是多組分氧化還原酶系統[31]，在真核生物的生長發育和新陳代謝中有至關重要的作用[32]。一般來說，每一個氧化酶反應都需要一個專一的P450酶，所以在真核生物中存在大量的P450s。生物體中有很多CYPs已經鑒定出來，并且功能得到了驗證。本研究在橡膠樹中克隆出一個P450基因，其編碼的蛋白在3～511位具有P450保守的超級家族結構域，在450～459位氨基酸存在一段血紅素結合位點，序列為FGAGRRMCPG，與P450家族保守的血紅素結合域FxxGxRxCxG序列相符。氨基酸序列比對結果顯示，橡膠樹HbCYP450的氨基酸序列與蓖麻、胡楊、毛楊果、大豆的一致性分別達到74%、68%、68%、65%，依據P450基因家族分類原則[13]，初步判斷橡膠樹HbCYP450屬于P450家族中CYP89A2亞家族，系統進化關系分析也發現，橡膠樹HbCYP450與擬南芥CYP89A2的親緣關系最近，兩個結果相一致。

橡膠樹HbCYP450基因在花、葉和膠乳中均有表達，但在樹皮中不表達，在膠乳中表達量最高，推測HbCYP450基因在橡膠樹合成膠乳的過程中起著至關重要的作用，而在木質素形成中不起作用，或作用很小。通過NCBI發現，擬南芥中CYP45089A2具有運輸電子的活性、單加氧酶活性、鐵離子結合、氧結合、亞鐵血紅素結合的功能，參與氧化還原反應。在同一亞族的P450基因具有相似的功能，因此HbCYP450在橡膠樹中可能也具有以上的作用。P450催化代謝途徑可以產生重要的次生代謝產物，該次生代謝物質可增強植物抵抗病蟲害和非生物脅迫的能力[33]。通過表達分析發現白粉菌侵染、機械傷害和干旱這些生物和非生物脅迫均能使HbCYP450基因表達量顯著上調。肉桂酸-4-羥化酶（CHC）作為CYP450的一種（CYP73），其活性和基因的表達受創傷、金屬離子、病原菌侵染等因素的調控[34-35]。HbCYP450與CYP73在進化上屬于同一分支，因此可以推測HbCYP450參與病害和逆境脅迫響應。P450參與許多萜類化合物的生物合成，包括赤霉素、脫落酸和一些植物防御物質[12]，還參與生長素的合成[36]。赤霉素、脫落酸和生長素都是常見且重要的植物激素，在本研究中HbCYP450基因在5類植物激素和H2O2作用下，表達量均顯著性變化。說明HbCYP450基因參與植物激素的響應。本研究為探究橡膠樹CYP450家族成員和各成員在橡膠樹抗病、抗逆機制奠定了一定的基礎。

參考文獻

[1] Ohkawa H， Imaishi H， Shiota N， et al. Molecular mechanisms of herbicide resistance with special emphasis on cytochrome P450 monooxygenases[J]. Plant Biotech， 1998， 15（4）： 173.

[2] 冷欣夫， 邱星輝. 細胞色素P450酶系的結構、 功能與應用前景[M]. 北京： 科學出版社， 2001.

[3] Ehlting J， Hamberger B， Million-Rousseau R， et al. Cyto-chromes P450 in phenolic metabolism[J]. Phytochem Rev， 2006， 5（2-3）： 239.

[4] Tanaka Y. Flower colour and cytochromes P450[J]. Phyto-chem Rev， 2006， 5（2-3）： 283.

[5] Ralston L， Yu O. Metabolons involving plant cytochromeP450s[J]. Phytochem Rev， 2006， 5（2-3）： 459.

[6] Bourgaud F， Hehn A， Larbat R， et al. Biosynthesis of cou-marins in plants： a major pathway still to be unravelled forcytochrome P450 enzymes[J]. Phytochem Rev， 2006， 5（2-3）： 293.

[7] Ayabe S I， Akashi T. Cytochrome P450s in flavonoid metab-olism[J]. Phytochem Rev， 2006， 5（2-3）： 271.

[8] Nomura T， Bishop G J. Cytochrome P450s in plant steroid hormone synthesis and metabolism[J]. Phytochem Rev， 2006， 5（2-3）： 421.

[9] Paquette S M， Bak S， Feyereisen R. Intron-exon organization and phylogeny in a large superfamily， the paralogous cytochrome P450 genes of Arabidopsis thaliana[J]. DNA Cell Biol， 2000， 19（5）： 307-317.

[10] Murphy P J， West C A. The role of mixed function oxidases in kaurene metabolism in Echinocystis macrocarpa Greene endosperm[J]. Arch Biochem Biophys， 1969， 133（2）： 395-407.

[11] Lew F L， West C A. （-）-kaur-16-en-7β-ol-19-oic acid， an intermediate in gibberellin biosynthesis[J]. Phytochemistry， 1971， 10（9）： 2 065-2 076.

[12] 賀麗虹， 趙淑娟， 胡之璧. 植物細胞色素P450基因與功能研究進展[J]. 藥物生物技術， 2008， 15（2）： 142-147.

[13] Nelson D R， Koymans L， Kamataki T， et al. P450 superfamily： update on new sequences， gene mapping， accession numbers and nomenclature[J]. Pharmacogenetics， 1996， 6（1）： 1-42.

[14] Nelson D R， Werck-Reichhart D. A P450-centric view of plant evolution[J]. Plant J， 2011， 66（1）： 194-211.

[15] 李新梅. 干旱期橡膠樹的管理[J]. 致富天地， 2010（5）： 31.

[16] 溫福光. 物候因子干旱災害長期預報[J]. 災害學， 1991， 6 （1）： 89-94.

[17] Zhong R， Morrison W H， Negrel J， et al. Dual methylation pathways in lignin biosynthesis[J]. Plant Cell， 1998， 10： 2 033-2 046.

[18] Jones L， Ennos A R， Turner S R. Cloning and characterization of irregular xylem4 （irx4）： a severely lignin-deficient mutant of A. thaliana[J]. Plant J， 2001， 26： 205-216.

[19] Goujon T， Ferret V， Mila I， et al. Down-regulation of AtCCR1 gene in A. thaliana： effects on phenotype， lignins and cel wa degradability[J]. Planta， 2003， 217： 218-228.

[20] 楊家書， 吳 畏， 吳友三， 等. 植物苯丙酸類代謝與小麥對白粉病抗性的關系[J]. 植物病理學報， 1986， （03）： 169-175.

[21] 段翠芳， 曾日中， 黎 瑜. 激素對巴西橡膠樹橡膠生物合成的調控[J]. 熱帶農業科學， 2004， （05）： 61-68.

[22] Qu X， Pu X， Chen F， et al. Molecular Cloning， Heterologous Expression， and Functional Characterization of an NADPH-Cytochrome P450 Reductase Gene from Camptotheca acuminata， a Camptothecin-Producing Plant[J]. PLoS One， 2015， 10（8）： e0135397.

[23] Wang L F. Physiological and molecular responses to drought stress in rubber tree （Hevea brasiliensis Muell. Arg.）[J]. Plant Physiol Biochem， 2014， 83： 243-249.

[24] Piffanelli P， Zhou F Z， Casais C， et al. The barley MLO modulator of defense and cell death is responsive to biotic and abiotic stress stimuli[J]. Plant Physiol， 2002， 129（3）： 1 076-1 085.

[25] Wang L F， Wang M， Zhang Y. Effects of powdery mildew infection on chloroplast and mitochondrial functions in rubber tree[J]. Trop Plant Pathol， 2014， 39（3）： 242-250.

[26] Qin B， Zheng F， Zhang Y， 2015. Molecular cloning and characterization of a Mlo gene in rubber tree （Hevea brasiliensis）[J]. Plant Physiol 175， 78-85.

[27] 何 鑫. 巴西橡膠樹JAZ和MYC家族幾個成員基因表達和產量相關性的研究[D]. ?？冢?海南大學， 2013.

[28] 莊海燕. 巴西橡膠樹水通道蛋白基因cDNA的克隆及其在乙烯利刺激下表達的初步分析[D]. 陜西： 西北農林科技大學， 2010.

[29] Chao J， Chen Y， Wu S， et al.， Comparative transcriptome analysis of latex from rubber tree clone CATAS 8-79 and PR107 reveals new cues for the regulation of latex regeneration and duration of latex flow[J]. BMC Plant Biol， 2015， 15： 104.

[30] 康桂娟， 黎 瑜， 曾日中. 巴西橡膠樹HbNAC24基因克隆和表達分析[J]. 西北植物學報， 2014， 34（12）： 2 374-2 381.

[31] Wang M， Roberts D L， Paschke R， et al. Three-dimensional structure of NADPH cytochrome P450 reductase： prototype for FMN- and FAD-containing enzymes[J]. Proc Natl Acad Sci USA. 1997， 94： 8 411-8 416.

[32] Rana S， Latto S K， Dhar N， et al. NADPH-cytochrome P450 reductase： molecular cloning and functional characterization of two paralogs from Withania somnifera （L.） Dunal[J]. PLoS ONE， 2013， 8： e57 068.

[33] 趙 劍， 楊文杰， 朱蔚華. 細胞色素與植物的次生代謝[J]. 生命科學， 1999， 11（3）： 127-131.

[34] Bell-Lelong D A， Cusumano JC， Meyer K， et al. Cinnamate-4-hydroxylase expression in Arabidopsis. Regulation in response to development and the environment[J]. Plant Physiol， 1997， 113（3）： 729-738.

[35] 王亞玲， 李 群， 何祖華. 一個水稻 P450 CYP72A基因簇受病原菌誘導和發育及組織特異性表達[J]. 科學通報， 2003， 48 （21）： 2 266-2 270.

[36] Kenneth A Feldmann. Cytochrome P450s as genes for crop improvement[J]. Curr Opin Plant Biol， 2001， 4（2）： 162-167.