李輝亮　郭冬　汪晗　楊子平　彭世清

摘 要 通過(guò)同源克隆和RACE方法獲得了巴西橡膠樹(shù)DNA甲基轉(zhuǎn)移酶（DNA Methyltransferase， MET）基因（HbMET）。HbMET長(zhǎng)4 896 bp，含有4 635 bp的閱讀框架，編碼1 545個(gè)氨基酸。蛋白分子量174.36 ku，等電點(diǎn)為6.36。氨基酸序列與可可、毛果楊、黃瓜、擬南芥、碧桃、葡萄和煙草等物種的MET家族成員的同源性分別為73%、77%、74%、56%、72%和68%。進(jìn)化樹(shù)分析表明，HbMET氨基酸序列與毛果楊的MET親緣關(guān)系最近。HbMET在巴西橡膠樹(shù)的根、樹(shù)皮、葉、膠乳中均有表達(dá)，其中在膠乳中表達(dá)量最低；HbMET在橡膠樹(shù)自根幼態(tài)無(wú)性系的膠乳中表達(dá)比老態(tài)無(wú)性系膠乳中的表達(dá)低。

關(guān)鍵詞 巴西橡膠樹(shù)；自根幼態(tài)無(wú)性系；DNA甲基轉(zhuǎn)移酶；表達(dá)；DNA甲基化

中圖分類(lèi)號(hào) S794.1 文獻(xiàn)標(biāo)識(shí)碼 A

目前生產(chǎn)中種植的橡膠樹(shù)是通過(guò)對(duì)優(yōu)良品種嫁接獲得的，由于在嫁接過(guò)程中所用的砧木是未經(jīng)選擇的實(shí)生苗，同時(shí)嫁接后由于砧木和接穗的相互作用，造成種植的橡膠樹(shù)會(huì)出現(xiàn)不同個(gè)體之間生長(zhǎng)的差異和產(chǎn)膠的差異[1-3]。中國(guó)學(xué)者首次以橡膠樹(shù)花藥為外植體通過(guò)體細(xì)胞發(fā)生途徑獲得了再生植株[4-5]。由于再生植株具有自己的根系，擺脫了砧木的影響并恢復(fù)了幼態(tài)的特性，稱(chēng)為自根幼態(tài)無(wú)性系[2，6]。與嫁接苗相比，自根幼態(tài)無(wú)性系具有生長(zhǎng)速度快、干膠產(chǎn)量高、抗逆性強(qiáng)的特點(diǎn)[6-9]。自根幼態(tài)無(wú)性系與供體（老態(tài)無(wú)性系）差異產(chǎn)生的原因尚不清楚。以前對(duì)差異產(chǎn)生的原因研究主要集中在形態(tài)學(xué)觀察、排膠生理參數(shù)等方面，如樹(shù)的莖圍、乳管層數(shù)、堵塞指數(shù)和黃色體破裂指數(shù)、膠乳中糖、無(wú)機(jī)磷含量和干膠含量等，綜合相關(guān)研究結(jié)果，自根幼態(tài)無(wú)性系高產(chǎn)的構(gòu)成因子主要有：較快的生長(zhǎng)速度、較多的乳管列數(shù)目、較好的排膠和產(chǎn)膠生理特性[6-9]。本課題組研究發(fā)現(xiàn)在自根幼態(tài)無(wú)性系的膠乳和老態(tài)無(wú)性系的膠乳中存在一些基因和蛋白的差異表達(dá)，初步證實(shí)了在自根幼態(tài)無(wú)性系中膠乳上調(diào)表達(dá)的基因有95個(gè)，下調(diào)表達(dá)的基因有26個(gè)，同時(shí)發(fā)現(xiàn)在橡膠樹(shù)自根幼態(tài)無(wú)性系中13個(gè)蛋白上調(diào)表達(dá)，11個(gè)蛋白下調(diào)表達(dá)[10-15]。推測(cè)這些基因和蛋白的差異表達(dá)是自根幼態(tài)無(wú)性系和老態(tài)無(wú)性系的產(chǎn)量差異產(chǎn)生的根本原因。橡膠樹(shù)自根幼態(tài)無(wú)性系是由供體通過(guò)無(wú)性繁殖獲得的，自根幼態(tài)無(wú)性系和供體的遺傳背景是相同的，所含遺傳物質(zhì)DNA在堿基序列是相同的，正常條件下橡膠樹(shù)的無(wú)性繁殖不可能改變遺傳物質(zhì)DNA的堿基序列，因此自根幼態(tài)無(wú)性系與其供體差異的產(chǎn)生很可能與表觀遺傳相關(guān)。李輝亮[16]證實(shí)了橡膠樹(shù)自根幼態(tài)無(wú)性系與供體基因組DNA之間存在甲基化多態(tài)性的差異，進(jìn)一步表明自根幼態(tài)無(wú)性系與其供體差異的產(chǎn)生與表觀遺傳相關(guān)。

表觀遺傳是指不因DNA序列的變化而形成基因功能的改變，這種改變是細(xì)胞內(nèi)除了遺傳信息以外的其他可遺傳物質(zhì)發(fā)生的改變，且這種改變?cè)诎l(fā)育和細(xì)胞增殖過(guò)程中能穩(wěn)定傳遞[17]。表觀遺傳的分子機(jī)制包括DNA甲基化，組蛋白修飾，染色質(zhì)改型和小RNA干擾等。甲基化是DNA化學(xué)修飾的一種形式，能夠在不改變DNA序列的前提下，改變遺傳表現(xiàn)，是調(diào)節(jié)基因功能的重要手段[18]。它通過(guò)一系列DNA甲基轉(zhuǎn)移酶來(lái)建立和維持DNA甲基化模式[19]。在植物中目前發(fā)現(xiàn)至少有MET甲基轉(zhuǎn)移酶（Methyltransferase， MET）、染色質(zhì)甲基化酶（Chromomethylase， CMT）和域重排甲基轉(zhuǎn)移酶（Domains Rearranged Methyltransferase， DRM）3類(lèi)在結(jié)構(gòu)和功能上不同的胞嘧啶甲基轉(zhuǎn)移酶[19-27]。為了研究橡膠樹(shù)自根幼態(tài)無(wú)性系高產(chǎn)與DNA甲基化的關(guān)系，本課題組對(duì)橡膠樹(shù)的3類(lèi)胞嘧啶甲基轉(zhuǎn)移酶基因進(jìn)行了克隆和表達(dá)分析等相關(guān)研究。本研究將報(bào)道巴西橡膠樹(shù)DNA甲基轉(zhuǎn)移酶（Methyltransferase， MET）基因（HbMET）的克隆和表達(dá)，該結(jié)果為進(jìn)一步研究橡膠樹(shù)自根幼態(tài)無(wú)性系高產(chǎn)與DNA甲基化的關(guān)系奠定基礎(chǔ)。

1 材料與方法

1.1 材料

1.1.1 植物材料 巴西橡膠樹(shù)（Hevea brasiliensis）海墾2號(hào)自根幼態(tài)無(wú)性系及其供體老態(tài)無(wú)性系種植在中國(guó)熱帶農(nóng)業(yè)科學(xué)院試驗(yàn)場(chǎng)。所用膠乳的采集按照Tang等[28]的方法進(jìn)行。根、葉和樹(shù)皮采集后用液氮研磨，保存于-70 ℃?zhèn)溆谩?/p>

1.1.2 質(zhì)粒、菌種及試劑 3′-full RACE試劑盒和RNA PCR（AMV）Ver.3.0試劑盒購(gòu)自TaKaRa公司。Taq DNA Polymerase 和T4-DNA連接酶dNTPs、IPTG、X-gal購(gòu)于Promega公司。E.coli DH5α菌株由本實(shí)驗(yàn)室保存。

1.2 方法

1.2.1 RNA的提取 參照Tang等[28]的方法進(jìn)行。

1.2.2 HbMET的擴(kuò)增 根據(jù)已報(bào)道的植物MET基因家族保守序列設(shè)計(jì)簡(jiǎn)并引物，通過(guò)PCR獲得HbMET片段序列，再根據(jù)HbMET片段序列設(shè)計(jì)3′RACE引物。cDNA的合成和3′RACE擴(kuò)增按照3′-full RACE試劑盒使用說(shuō)明書(shū)進(jìn)行。實(shí)驗(yàn)所用引物見(jiàn)表1。

1.2.3 HbMET同源性分析和進(jìn)化樹(shù)構(gòu)建 利用DNAMAN軟件分析HbMET蛋白與其它植物MET家族的同源性和進(jìn)化樹(shù)構(gòu)建。

1.2.4 HbMET定量PCR分析 以HbACT基因?yàn)閮?nèi)參對(duì)照，采用SYBR Premix Ex TaqⅡ進(jìn)行定量PCR分析，基因的相對(duì)表達(dá)量按BIO-RAD CFX96熒光定量?jī)x使用說(shuō)明進(jìn)行分析。反應(yīng)條件：95 ℃預(yù)變性30 s；95 ℃ 5 s， 60 ℃ 20 s，72 ℃ 20 s，35個(gè)循環(huán)。PCR引物見(jiàn)表2。

2 結(jié)果與分析

2.1 HbMET的克隆

根據(jù)已有植物的MET基因家族保守序列設(shè)計(jì)簡(jiǎn)并引物，通過(guò)PCR獲得了HbMET的3個(gè)片段，分別為1 725 bp（HbMET1）、1 584 bp（HbMET2）和1 525 bp（HbMET3）（圖1）。利用DNAMAN軟件對(duì)這3個(gè)片段進(jìn)行拼接。根據(jù)拼接獲得的片段設(shè)計(jì)3′RACE引物，通過(guò)3′-RACE擴(kuò)增技術(shù)得到一個(gè)551 bp的 3′片段（圖1）。將保守區(qū)片段與3′端片段進(jìn)行拼接，得到一個(gè)長(zhǎng)4 896 bp的序列（命名為HbMET）。該序列含有4 635 bp的開(kāi)放式閱讀框架，3′端有248 bp的非編碼區(qū)，其中含有一條18個(gè)堿基的polyA尾，預(yù)測(cè)該基因編碼1 545個(gè)氨基酸（圖2），分子量174.36 ku，等電點(diǎn)為6.36。

2.2 HbMET生物信息學(xué)分析

將HbMET氨基酸序列與其它植物的MET比較，發(fā)現(xiàn)HbMET與可可（Theobroma cacao）、毛果楊（Populus trichocarpa）、黃瓜（Cucumis sativus）、擬南芥（Arabidopsis thaliana）、碧桃（Prunus persica）、葡萄（Vitis vinifera）和煙草（Nicotiana tabacum）等物種的MET家族成員的同源性分別達(dá)到為73%、77%、74%、56%、72%和68%（圖3）。構(gòu)建的進(jìn)化樹(shù)表明，HbMET與毛果楊的MET親緣關(guān)系最近（圖4）。

2.3 HbMET在自根幼態(tài)無(wú)性系組織特異性表達(dá)

對(duì)HbMET在自根幼態(tài)無(wú)性系的根、皮、花、葉和膠乳中的表達(dá)進(jìn)行了定量PCR分析，結(jié)果表明，HbMET在根中的表達(dá)最高，其次是在葉和樹(shù)皮中表達(dá)，膠乳中表達(dá)量最低（圖5）。

2.4 HbMET在自根幼態(tài)無(wú)性系與其供體膠乳的差異表達(dá)

天然橡膠是在橡膠樹(shù)的乳管中合成的，膠乳實(shí)質(zhì)上是乳管細(xì)胞的細(xì)胞質(zhì)[29]。乳管細(xì)胞中基因的表達(dá)與天然橡膠生物合成直接相關(guān)。為了解HbMET在自根幼態(tài)無(wú)性系與其供體膠乳的表達(dá)是否存在差異，對(duì)HbMET在自根幼態(tài)無(wú)性系和其供體膠乳的進(jìn)行定量PCR分析， 結(jié)果表明，HbMET在自根幼態(tài)無(wú)性系的表達(dá)比其供體膠乳的表達(dá)低（圖6）。

3 討論與結(jié)論

DNA甲基化是最早發(fā)現(xiàn)的表觀遺傳修飾途徑之一，大量研究結(jié)果表明，DNA甲基化能引起染色質(zhì)結(jié)構(gòu)、DNA構(gòu)象、DNA穩(wěn)定性及DNA與蛋白質(zhì)相互作用方式的改變，從而控制基因表達(dá)[17]。DNA甲基化主要是通過(guò)DNA甲基轉(zhuǎn)移酶家族來(lái)催化完成的[19-21]。目前一些植物的MET基因已經(jīng)克隆并進(jìn)行了相關(guān)研究[21-23]，但是到目前為止還未見(jiàn)橡膠樹(shù)MET基因研究的報(bào)道。本研究首次從橡膠樹(shù)中克隆了MET基因HbMET，同源性和分子進(jìn)化分析表明，推導(dǎo)的HbMET是一種與植物DNA甲基化相關(guān)的甲基化酶。過(guò)去關(guān)于巴西橡膠樹(shù)自根幼態(tài)無(wú)性系和供體差異產(chǎn)生的原因的研究主要集中在對(duì)自根幼態(tài)無(wú)性系高產(chǎn)構(gòu)成的形態(tài)和生理指標(biāo)[6-9]，以及膠乳中差異表達(dá)基因和蛋白的鑒定[10-15]。而從表觀遺傳的角度來(lái)分析自根幼態(tài)無(wú)性系和供體差異有可能闡明自根幼態(tài)無(wú)性系高產(chǎn)產(chǎn)生的根本原因。天然橡膠是橡膠樹(shù)乳管細(xì)胞中合成的，膠乳是乳管細(xì)胞的細(xì)胞質(zhì)[29]。本研究結(jié)果表明，HbMET在自根幼態(tài)無(wú)性系膠乳中的表達(dá)比其供體膠乳的表達(dá)低，可能在自根幼態(tài)無(wú)性系乳管細(xì)胞中DNA甲基化程度較供體乳管細(xì)胞的低，DNA甲基化可調(diào)控特定的基因表達(dá)，而植物啟動(dòng)子區(qū)的DNA甲基化通常抑制轉(zhuǎn)錄[21-22]。推測(cè)由于自根幼態(tài)無(wú)性系和供體乳管細(xì)胞中DNA甲基化的差異，使得自根幼態(tài)無(wú)性系和供體乳管細(xì)胞中基因的表達(dá)出現(xiàn)差異，從而導(dǎo)致巴西橡膠樹(shù)自根幼態(tài)無(wú)性系和供體在產(chǎn)膠能力等方面的差異，但這一推測(cè)還需進(jìn)一步深入研究。

參考文獻(xiàn)

[1] Clément-Demange A， Priyadarshan P M， Hoa T T T， et al. Hevea rubber breeding and genetics[J]. Plant Breeding Reviews， 2007， 29： 177-281.

[2] Chandrashekar T R， Mydin K K， Alice J， et al. Intraclonal variability for yield in rubber（Hevea brasiliensis）[J]. Ind J Nat Rubb Res， 1997， 10： 43-47.

[3] Hua Y W， Huang T D， Huang H S. Micropropagation of self-rooting juvenile clones by secondary somatic embryogenesis in Hevea brasiliensis[J]. Plant Breeding， 2010， 129： 202-207.

[4] 中國(guó)科學(xué)院遺傳研究所， 廣東省農(nóng)墾總局保亭熱帶作物研究所. 三葉橡膠（Hevea brasiliensis）花藥育株成功（簡(jiǎn)訊）[J]. 遺傳學(xué)報(bào)， 1977， 4： 186.

[5] 王澤云， 曾憲松， 陳傳琴. 從離體花藥誘導(dǎo)巴西橡膠植株[J]. 熱帶作物科技通訊， 1978， 4： 1-7.

[6] 陳雄庭， 王澤云， 吳蝴蝶， 等. 橡膠樹(shù)新種植材料-幼態(tài)自根無(wú)性系[J]. 熱帶作物學(xué)報(bào)， 2002， 23： 19-23.

[7] 幼態(tài)課題組. 幼態(tài)無(wú)性系在生產(chǎn)上的應(yīng)用[J]. 熱帶農(nóng)業(yè)科學(xué)， 1994， 14（2）： 1-3.

[8] 郝秉忠， 吳繼林. 巴西橡膠樹(shù)幼態(tài)無(wú)性系的高產(chǎn)特性[J]. 熱帶農(nóng)業(yè)科學(xué)， 1996， 16（2）： 1-8.

[9] 王澤云， 陳雄庭， 吳胡蝶. 橡膠樹(shù)新型種植材料--體胚植株[J]. 熱帶農(nóng)業(yè)科學(xué)， 2001， 21（6）： 11-15.

[10] 魯惠中， 李輝亮， 郭 冬， 等. 巴西橡膠樹(shù)HbHDT1基因的克隆及表達(dá)分析[J]. 熱帶作物學(xué)報(bào)， 2010， 31（1）： 59-64.

[11] 李輝亮， 郭冬， 彭世清. 巴西橡膠樹(shù)14-3-3蛋白基因的克隆及表達(dá)分析[J]. 熱帶作物學(xué)報(bào)， 2010， 31（8）： 1 452-1 457.

[12] 梁小蓮， 李輝亮， 彭世清. 巴西橡膠樹(shù)TCTP基因的克隆及表達(dá)[J]. 分子植物育種， 2009， 7： 188-193.

[13] Li H L， Guo D， Peng S Q. Differential gene expression profiles in latex from Hevea brasiliensis between self-rooting juvenile and donor clones[J]. Plant Growth Regulation， 2014， 74： 65-71.

[14] Li H L， Guo D， Lan F Y， et al. Protein differential expression in the latex from Hevea brasiliensis between self-rooting juvenile clones and donor clones[J]. Acta Physiol Plant， 2011， 33： 1 853-1 859.

[15] Li H L， Lu H Z， Guo D， et al. Molecular characterization of a thioredoxin h gene（HbTRX1）from Hevea brasiliensis showing differential expression in latex between self-rooting juvenile clones and donor clones[J]. Mol Biol Rep， 2011， 38： 1 989-1 994.

[16] 李輝亮. 橡膠樹(shù)自根幼態(tài)無(wú)性系與其供體（老態(tài)無(wú)性系）差異的分子基礎(chǔ)研究[D]. ?？冢?海南大學(xué)， 2010.

[17] Russo V E A， Martienssen R A， Riggs A D. Epigenetic mechanisms of gene regulation[M]. New York： Cold Spring Harbor Laboratory Press， 1996： 123-130.

[18] Zemach A， McDaniel I E， Silve P， et al. Genome-wide evolutionary analysis of eukaryotic DNA methylation[J]. Science， 2010， 328： 916-919.

[19] Goll M G， Bestor T H. Eukaryotic cytosine methyltransferases[J]. Annual Review Biochem， 2005， 74： 481-514.

[20] Finnegan E J， Kovac K A. Plant DNA methyltransferases[J]. Plant Mole Biol， 2000， 43： 189-201.

[21] Pavlopoulou A， Kossida S. Plant cytosine-5 DNA methyltransferases： structure， function and molecular evolution[J]. Genomics， 2007， 90： 530-541.

[22] Kishimoto N， Sakai H， Jackson J， et al. Site specificity of the Arabidopsis METI DNA methyltransferase demonstrated through hypermethylation of the SUPERMAN locus[J]. Plant Mol Biol， 2001， 46： 171-183.

[23] Chan S W， Henderson I R， Jacobsen S E. Gardening the genome： DNA methylation in Arabidopsis thaliana[J]. Nature Reviews Genetics， 2005， 6： 351-360.

[24] Henikoff S， Comai L. A DNA methyltransferase homolog with a chromodomain exists in multiple polymorphic forms in Arabidopsis[J]. Genetics， 1998， 149： 3 007-3 018.

[25] Cao X， Springer N M， Muszynski M， et al. Conserved plant genes with similarity to mammalian de novo DNA methyltranferases[J]. Proc Nat Acad Sci USA ， 2000， 97： 4 979-4 984.

[26] Cao X， Jacobsen S E. Role of the Arabidopsis DRM methyltransfeases in de novo methylation and gene silencing[J]. Current Biol， 2002， 12： 1 138-1 144.

[27] Chan S W， Zilberman D， Xie Z， et al. RNA silencing genes control de novo DNA methylation[J]. Science， 2004， 303： 1 336.

[28] Tang C， Qi J， Li H， et al. A convenient and efficient protocol for isolating high-quality RNA from latex of Hevea brasiliensis（para rubber tree）[J]. Biochem Biophys Methods， 2007， 70： 749-754.

[29] Kush A. Isoprenoid biosynthesis： the Hevea factory[J]. Plant Physiol Biochem， 1994， 32： 761-767.