張濤，司福會，張玉喜，蓋樹鵬

?

外源GA3影響牡丹花芽DNA甲基化水平和相關(guān)酶基因的表達分析

張濤，司福會，張玉喜，蓋樹鵬

（青島農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院/山東省高校植物生物技術(shù)重點實驗室，山東青島 266109）

【目的】外施赤霉素（GA）是生產(chǎn)上輔助解除牡丹花芽休眠進行反季節(jié)盆花生產(chǎn)的常用手段，赤霉素如何發(fā)揮作用解除牡丹花芽休眠是未解的科學(xué)問題，本文旨在明確赤霉素是否通過表觀遺傳方式參與內(nèi)休眠調(diào)控?！痉椒ā繌腄NA甲基化入手，以4年生‘魯菏紅’為材料，500 mg·L-1的GA3處理花芽，并在處理后0.5、1、3和5 d后取健壯頂芽，CTAB法提取花芽DNA，HPLC技術(shù)測定GA3處理后牡丹花芽DNA甲基化水平的變化；利用轉(zhuǎn)錄組分析結(jié)合RACE技術(shù)得到牡丹3種DNA甲基轉(zhuǎn)移酶基因、、和DNA去甲基化酶基因的序列，生物信息學(xué)方法比對基因序列并分析可能的結(jié)構(gòu)域，MEGA 5.0構(gòu)建系統(tǒng)發(fā)育樹；以為內(nèi)參，利用qPCR方法分析其組織表達特性和對GA3的響應(yīng)?！窘Y(jié)果】3種DNA甲基轉(zhuǎn)移酶基因（）的開放閱讀框長短不一，、、開放閱讀框長度分別為1 118、1 056、2 175 bp，編碼的氨基酸數(shù)目分別為372、351、724。3種DNA甲基轉(zhuǎn)移酶均有甲基轉(zhuǎn)移酶結(jié)構(gòu)域。而DNA去甲基化酶基因的序列較長，開放閱讀框為6 636 bp，編碼2 211個氨基酸。上存在保守的DNA糖基化結(jié)構(gòu)域。系統(tǒng)發(fā)育分析表明，牡丹中的PsCMT和PsDRM蛋白與葡萄的VvCMT2蛋白親緣關(guān)系最近，PsMET和PsROS1與大部分木本植物聚為一支，與煙草等距離較遠。、、在牡丹初花期的各個組織（根、莖、葉、苞片、萼片、花瓣、雄蕊、心皮）中均有表達，其中在牡丹根中表達量較高，而在心皮中表達量最高。GA3處理5 d，牡丹花芽迅速萌動，60 d時萌動率為97.5%，成花率為92.5%；對照20 d時才有少量花芽萌動，至60 d時萌動率僅為23.1%。GA3顯著降低了牡丹花芽DNA甲基化水平（＜0.01），從處理前的38.9%降至處理5 d時的28.7%；GA3處理提高了和的表達水平，降低了的表達水平，而的表達水平差異不顯著?！窘Y(jié)論】GA3處理通過誘導(dǎo)表達、抑制表達導(dǎo)致了DNA低甲基化，進而促進了牡丹休眠解除，可能是赤霉素發(fā)揮作用的一種重要方式。

牡丹；DNA甲基轉(zhuǎn)移酶；；赤霉素；DNA甲基化；休眠解除

0 引言

【研究意義】牡丹（Andr.）是中國著名的傳統(tǒng)觀賞花卉之一，同其他多年生木本植物一樣，其在自然條件下必須經(jīng)過冬天一段時間低溫，才能解除休眠來年正常開花。足夠低溫累積是解除休眠的有效措施。低溫解除牡丹休眠進程中，甲基化水平是下降的[1]，暗示表觀遺傳修飾在其中發(fā)揮了重要作用。生產(chǎn)上，常用赤霉素輔助解除休眠，進行反季節(jié)花期調(diào)控生產(chǎn)。因此，明確赤霉素解除休眠是否與甲基化調(diào)控有關(guān)對于牡丹反季節(jié)生產(chǎn)具有重要意義?！厩叭搜芯窟M展】真核生物中，DNA胞嘧啶的甲基化與基因表達密切相關(guān)[2]，5-胞嘧啶的甲基化是最普遍的DNA修飾。胞嘧啶甲基化在多種生物發(fā)育過程中發(fā)揮關(guān)鍵作用，包括基因沉默、轉(zhuǎn)座子的失活和印記基因表達[3-4]。DNA甲基轉(zhuǎn)移酶（Dnmt）家族是植物進行甲基化修飾的工具，它能催化S-adenosylmethionine（SAM）上的甲基轉(zhuǎn)移到胞嘧啶嘧啶環(huán)的C-5位點，在表觀遺傳基因調(diào)控中起中心作用[5]。植物中的DNA甲基轉(zhuǎn)移酶主要包括甲基轉(zhuǎn)移酶、染色質(zhì)甲基轉(zhuǎn)移酶、結(jié)構(gòu)域重排甲基轉(zhuǎn)移酶這3類[6]。擬南芥中，主要維持重復(fù)序列和單拷貝DNA序列中CpG位點的DNA甲基化[7]。能維持重復(fù)和單拷貝DNA序列中CG類型的甲基化，并影響其形態(tài)特征[8]，是植物中特有的一類甲基化酶，可以在異染色質(zhì)發(fā)生DNA甲基化，其功能主要是維持CHG位點的DNA甲基化[2,9]。另外，植物DNA去甲基化酶（repressor of silencing 1）是一種DNA糖基化酶，它可以去除DNA5-甲基胞嘧啶上的甲基而發(fā)動DNA的去甲基化過程。目前，其功能和特性僅在擬南芥中進行了較為詳細的研究[10]?！颈狙芯壳腥朦c】關(guān)于赤霉素對甲基化影響的報道大都集中在模式植物上，在木本植物牡丹中的研究未見報道?！緮M解決的關(guān)鍵問題】探究赤霉素促進休眠解除與DNA甲基化變化的潛在聯(lián)系，為進一步研究赤霉素解除休眠的機制打下基礎(chǔ)。

1 材料與方法

1.1 試驗材料

供試材料為4年生‘魯菏紅’（‘Luhehong’）健壯植株，由菏澤牡丹研究所提供。2016年10月下旬定植于盆中，11月12日移入溫室，用500 mg·L-1的GA3處理花芽（棉花纏繞花芽，滴加GA3溶液1 mL），并在處理后0.5、1、3和5 d后取健壯頂芽，每組設(shè)3個重復(fù)，清水為對照；每處理15株，共90株，剝?nèi)[片后，液氮速凍后于-80℃保存?zhèn)溆?。另取上述處?0株，移入溫室后60 d調(diào)查植株生長情況。取初花期（大風(fēng)鈴期）牡丹花芽的各個組織（根、莖、葉、花瓣、萼片、苞片、雄蕊、心皮），每組取3個重復(fù)，用于組織表達分析。

1.2 總RNA提取與cDNA合成

采用MiniBEST Plant RNA Extraction Kit（TaKaRa，大連）提取總RNA，-80℃保存。cDNA的合成按照PrimeScriptTMRT reagent Kit with gDNA Eraser 試劑盒（TaKaRa，大連）說明書進行。RACE- Ready-cDNA第一鏈合成按照SMARTer? RACE 5′/3′ Kit說明書進行。

1.3 RACE擴增和實時定量PCR

根據(jù)已知序列設(shè)計特異性引物進行RACE擴增。

反應(yīng)條件：94℃ 30 s，72℃ 3 min，共5個循環(huán)；94℃ 30 s，70℃30 s，72℃3 min，共5個循環(huán)；94℃30 s，68℃30 s，72℃3 min，共25個循環(huán)。

將反轉(zhuǎn)錄得到的cDNA稀釋3倍作為實時定量PCR反應(yīng)的模板，設(shè)定3個技術(shù)重復(fù)。為內(nèi)參，引物見表1，數(shù)據(jù)分析采用2-ΔΔCT法。Real-time PCR反應(yīng)體系為20 μL，包括ddH2O 6.8 μL，SYBR Premix DimerEraser（2×）10.0 μL，PCR Forward Primer（10 μmol·L-1）0.6 μL，PCR Reverse Primer（10 μmol·L-1）0.6 μL，cDNA 2.0 μL。反應(yīng)在Light Cycler 480 real-time PCR System上進行，反應(yīng)條件參考司福會等[11]。

表1 引物相關(guān)信息

1.4 生物信息學(xué)分析

測序結(jié)果利用DNAMAN6.0進行分析和拼接，利用Blast（http://blast.ncbi.nlm.nih.gov/Blast.cgi）進行蛋白質(zhì)Blast對比并分析結(jié)構(gòu)域。利用IBS 1.0.2繪制蛋白質(zhì)結(jié)構(gòu)域分布圖。系統(tǒng)發(fā)育樹采用MEGA 5.0中的鄰接法（neighbor-joining，NJ）構(gòu)建，置信度檢驗設(shè)置為Bootstrap=1 000。

1.5 甲基化水平的測定

利用CTAB法提取植物總DNA[12]，100℃加熱2 min使約5 μg DNA變性，迅速置于冰上。甲基化測定用Hasbun等的方法[13]，變性的DNA在37℃經(jīng)核酸酶P1（Sigma）和堿性磷酸酶（Sigma）處理后，15 000×離心20 min，上清液經(jīng)0.025 μm濾紙（Millipore，Bedford，MA）過濾，上樣Eclipse XDB-C18柱（5 μm，150×4.6 nm，Aglient）用于HPLC分析（Aglient 1100 Series），流動相為KH2PO3/ methanol（98﹕2，v/v），約20 min。標樣胞嘧啶（C）和甲基胞嘧啶（mC）購自Sigma。使用下列公式計算甲基胞嘧啶百分比：mC（%）=（mC/（C+mC））×100。

2 結(jié)果

2.1 赤霉素處理解除牡丹花芽休眠的效應(yīng)

GA3處理后，牡丹花芽迅速萌動，處理5 d時，大部分花芽即已萌動；至60 d時調(diào)查，萌動率為97.5%，成花率為92.5%。對照萌動較晚，20 d時才有少量花芽萌動，至60 d時萌動率僅為23.1%，未見花朵開放，僅有5.1%的成花率（圖1）?？梢?，赤霉素處理可在一定程度上代替低溫，促進休眠解除，顯著提高萌動率和成花率。

2.2 赤霉素對甲基化水平的影響

利用HPLC分析了GA3處理后牡丹花芽的甲基化水平。由圖2可以看出，GA3處理顯著降低了花芽的甲基化水平，且隨著施加GA3天數(shù)的增加，花芽甲基化水平呈下降趨勢。從38.9%降至處理5 d時的28.7%，而對照甲基化水平維持在37.5%—41.5%，相對穩(wěn)定。

A：對照 Control；B：500 mg·L-1 GA3處理 Treated with 500 mg·L-1 GA3

*P＜0.05；** P＜0.01。圖6同 The same as Fig. 6

2.3 甲基化相關(guān)酶基因的克隆

根據(jù)DNA甲基化相關(guān)酶基因保守序列，從轉(zhuǎn)錄組數(shù)據(jù)庫中篩選出相關(guān)片段[12]，其中（Accession number：MH038045）和（Accession number：MH087469）已包含完整的ORF，而（ID：JI448809）和（ID：JI448935）編碼序列不完整，進行了RACE擴增。根據(jù)目的基因的序列設(shè)計了特異性引物（表1），經(jīng)過RACE PCR擴增，得到了997 bp的5′ RACE PCR產(chǎn)物和1 245 bp的5′ RACE PCR產(chǎn)物（圖3）。測序后利用DNAMAN與已知cDNA進行序列拼接，最終獲得1 506 bp的cDNA序列（Accession number：MH038046）、1 582 bp的cDNA序列（Accession number：MH038044）。開放閱讀框的長度為1 118 bp，編碼的氨基酸數(shù)目為372，蛋白質(zhì)分子量為41.469 kD。開放閱讀框的長度為1 056 bp，編碼351個氨基酸，蛋白質(zhì)分子量為39.853 kD。開放閱讀框的長度為2 175 bp，編碼724個氨基酸，蛋白質(zhì)分子量為81.477 kD。含有6 636 bp的開放閱讀框，編碼2 211個氨基酸，分子量為246.675 kD。

在線分析牡丹和的序列表明，在152—269 aa間含有一個BAH結(jié)構(gòu)域（The Bromo-adjacent homology domain），在311—375 aa間存在Cyt_C5_DNA_methylase結(jié)構(gòu)域；在22—351 aa間存在一個SAM-dependent MTase C5-type結(jié)構(gòu)域；在序列上找到兩個結(jié)構(gòu)域，一個是位于232—276 aa間的UBA（Ubiquitin-associated domain）結(jié)構(gòu)域，另一個是在393—723 aa的SAM-dependent MTase DRM-type結(jié)構(gòu)域；上存在DNA糖基化結(jié)構(gòu)域、HhH-GPD 結(jié)構(gòu)域、Perm-CXXC結(jié)構(gòu)域和RRM-DME結(jié)構(gòu)域（圖3）。

為了解牡丹PsDRM、PsMET、PsCMT、PsROS1蛋白與其他植物相應(yīng)基因的進化關(guān)系，對牡丹中4種DNA甲基化相關(guān)酶進行系統(tǒng)發(fā)育分析（圖4）。結(jié)果表明，PsCMT和PsDRM與葡萄的VvCMT2親緣關(guān)系最近，與模式植物擬南芥親緣關(guān)系較遠；PsMET和PsROS1與大部分木本植物聚為一支，比如胡楊，與煙草等距離較遠。

圖3 DNA甲基化酶（A）和去甲基化酶（B）的結(jié)構(gòu)域

圖4 牡丹甲基化相關(guān)酶的系統(tǒng)發(fā)育分析

2.4 牡丹甲基化相關(guān)酶基因的組織表達特性

在牡丹的初花期（大風(fēng)鈴期）取各個組織（根、莖、葉、苞片、萼片、花瓣、雄蕊、心皮），利用熒光定量PCR對DNA甲基化酶基因的表達模式進行了分析。結(jié)果表明，、、在各個組織中均能檢測到表達，但表達量不同。和在根、莖、雄蕊和苞片中表達量相對較高，在葉片和心皮中較低；而在根中的表達量最高，在心皮中的表達量相對較高，其次為莖和萼片，在雄蕊中最低；在心皮中表達量最高，其次是萼片，在葉片和花瓣中表達較少（圖5）。

2.5 赤霉素對甲基化相關(guān)酶基因表達的影響

植物中的甲基化水平是由甲基化酶和去甲基化酶共同維持的[14]。因此，對相關(guān)酶基因進行熒光定量PCR分析（圖6）。結(jié)果表明，GA3處理的牡丹花芽與對照相比表達量差異不顯著（＞0.05）。用GA3處理后，的表達量在0.5 d內(nèi)就有了大幅度的提高，這表明GA3顯著地提高了的表達水平，可快速響應(yīng)GA3。GA3處理0.5—3 d，的表達明顯減弱。與此相反，經(jīng)GA3處理后，去甲基化酶表達量顯著提高，且有逐漸增加的趨勢。

3 討論

DNA胞嘧啶甲基化、組蛋白修飾和RNA干擾（RNAi）是最常見的分子表觀遺傳學(xué)機制，在基因調(diào)控中發(fā)揮重要作用[14]。DNA胞嘧啶甲基化大約占總胞嘧啶的30%[15]，越來越受到人們的重視。DNA甲基化水平主要是由DNA甲基轉(zhuǎn)移酶和去甲基化酶創(chuàng)建和維持的。相關(guān)基因克隆和功能分析主要見于模式植物，但木本植物中相關(guān)基因克隆的報道較少。本研究從牡丹花芽中克隆了DNA甲基化酶基因、、和DNA去甲基化酶，序列分析表明具有相應(yīng)保守結(jié)構(gòu)域，并在不同組織中檢測到了它們的表達，暗示它們可能參與了牡丹細胞中DNA甲基化的調(diào)控。

1：根root；2：莖stem；3：葉Leaf；4：萼片sepal；5：苞片bract；6：花瓣petal；7：雄蕊stamen；8：心皮carpel

圖6 GA3處理后甲基化相關(guān)酶基因的表達水平

赤霉素是一種重要的植物激素，參與調(diào)控多種生長和發(fā)育過程[16]。目前已對赤霉素的生物合成和代謝途徑有了比較明確的了解[17]。然而，赤霉素在表觀遺傳調(diào)控中發(fā)揮怎樣的作用仍未可知。李梅蘭等[18]研究表明，在春化誘導(dǎo)不結(jié)球白菜開花期間，DNA甲基化水平降低與芽尖分生組織中內(nèi)源GA含量增加存在間接關(guān)系。本研究表明，赤霉素導(dǎo)致牡丹的甲基化水平下降。已有研究表明，休眠解除進程中伴隨著DNA低甲基化[1,19]。因此，赤霉素處理導(dǎo)致DNA低甲基化，可能是其發(fā)揮作用促進休眠解除的一種重要方式。

另外，赤霉素處理可能通過調(diào)節(jié)DNA甲基化相關(guān)酶基因的表達而影響甲基化水平。從小麥幼苗和發(fā)芽胚的細胞核中純化的DNMT的活性受到GA3抑制，但GA3添加到小麥胚的核提取物中似乎刺激了DNA甲基化酶活性[20]。然而，在水稻中，外施GA3不會引起表達的顯著變化[21]。近期在煙草中的試驗表明，和響應(yīng)GA3信號表達水平下降[14]。顯示了不同生物中的（、和）對GA響應(yīng)是有差異的[19-23]。ROS1介導(dǎo)的主動去甲基化可以導(dǎo)致內(nèi)外源基因的轉(zhuǎn)錄沉默[24]，引起5S rDNA染色質(zhì)的解聚及植物對生物脅迫和環(huán)境脅迫的響應(yīng)[25-27]。本研究中，GA3處理導(dǎo)致的表達量提高。CMT主要功能是維持CG位點甲基化，GA3處理促進了細胞分裂，需要CMT大量表達維持甲基化水平。DRM可導(dǎo)致CHG位點重頭甲基化，GA3處理后表達受到抑制，而表達量顯著提高，二者聯(lián)合作用，導(dǎo)致總體甲基化水平的下降。

4 結(jié)論

利用轉(zhuǎn)錄組分析和RACE技術(shù)，得到了牡丹DNA甲基轉(zhuǎn)移酶基因（、和）和DNA去甲基化酶基因，具有相應(yīng)保守結(jié)構(gòu)域。外施GA3顯著降低DNA甲基化水平，導(dǎo)致DNA低甲基化，進而促進了休眠解除，是赤霉素發(fā)揮作用的一種方式。DNA甲基化水平的降低可能與赤霉素誘導(dǎo)表達水平下降和表達量提高有關(guān)。

[1] 蓋樹鵬, 張風(fēng), 張玉喜, 鄭國生. 低溫解除牡丹休眠進程中基因組DNA甲基化敏感擴增多態(tài)性 (MSAP)分析. 農(nóng)業(yè)生物技術(shù)學(xué)報, 2012, 20(3): 261-267.

GAI S P, ZHANG F, ZHANG Y X, ZHENG G S. Analysis of genomic DNA methylation during chilling induced endo-dormancy release by methylation sensitive amplified polymorphism (MSAP) technology in tree peony ()., 2012, 20(3): 261-267. (in Chinese)

[2] PAPA C M, SPRINGER N M, MUSZYNSKI M G, MEELEY R, KAEPPLER S M. Maize chromomethylaseis required for CpNpG methylation., 2001, 13(8): 1919-1928.

[3] CHAN S W L, HENDERSON I R, JACOBSEN S E. Gardening the genome: DNA methylation in., 2005, 6(5): 351-360.

[4] GEHRING M, HENIKOFF S. DNA methylation dynamics in plant genomes., 2007, 1769(5/6): 276-286.

[5] LYKO F. The DNA methyltransferase family: a versatile toolkit for epigenetic regulation., 2018, 19(2): 81-92.

[6] FINNEGAN E J, KOVAC K A. Plant DNA methyltransferases., 2000, 43(2/3): 189-201.

[7] KISHIMOTO N, SAKAI H, JACKSON J, JACOBSEN S E, MEYEROWITZ E M, DENNIS E S, FINNEGAN E J. Site specificity of theMETI DNA methyltransferase demonstrated through hypermethylation of thelocus., 2001, 46(2): 171-183.

[8] GENGER R K, KOVAC K A, DENNIS E S, PEACOCK W J, FINNEGAN E J. Multiple DNA methyltransferase genes in., 1999, 41(2): 269-278.

[9] LINDROTH A M, Cao X, Jackson J P, Zilberman D, McCallum C M, Henikoff S, JACOBSEN S E. Requirement offor maintenance of CpXpG methylation., 2001, 292(5524): 2077-2080.

[10] GONG Z, MORALES-RUIZ T, ARIZA R R, ROLDAN-ARJONA T, DAVID L, ZHU J K., a repressor of transcriptional gene silencing in, encodes a DNA glycosylase/lyase., 2002, 111(6): 803-814.

[11] 司福會, 張玉喜, 劉春英, 蓋偉玲, 戰(zhàn)新梅, 蓋樹鵬. 牡丹DNA甲基結(jié)合域蛋白基因的克隆和表達特性分析. 華北農(nóng)學(xué)報, 2017, 32(2): 66-70.

SI F H, ZHANG Y X, LIU C Y, GAI W L, ZHAN X M, GAI S P. Cloning and expression pattern of DNA methyl-binding domain gene-, 2017, 32(2): 66-70. (in Chinese)

[12] GAI S P, ZHANG Y X, MU P, LIU C Y, LIU S, DONG L, ZHENG G S. Transcriptome analysis of tree peony during chilling requirement fulfillment: assembling, annotation and markers discovering., 2012, 497(2): 256-262.

[13] HASBUN R, VALLEDOR L, RODRIGUEZ J, SANTAMARIA E, RIOS D, SANCHEZ M, CANAL M J, RODRIGUEZ R. HPCE quantification of 5-methyl-2′-deoxycytidine in genomic DNA: methodological optimization for chestnut and other woody species., 2008, 46(8): 815-822.

[14] MANOHARLAL R, SAIPRASAD G V S, ULLAGADDI C, KOVA?íK A. Gibberellin A3as an epigenetic determinant of global DNA hypo-methylation in tobacco., 2018, 62(1): 11-23.

[15] FINNEGAN E J, DENNIS E S. Isolation and identification by sequence homology of a putative cytosine methyltransferase from., 1993, 21(10): 2383-2388.

[16] GOCAL G F W, SHELDON C C, GUBLER F, MORITZ T, BAGNALL D J, MACMILLAN C P, LI S F, PARISH R W, DENNIS E S, WEIGEL D, KING R W.genes, flowering, and gibberellin signaling in., 2001, 127(4): 1682-1693.

[17] SUN T P. Gibberellin-GID1-DELLA: a pivotal regulatory module for plant growth and development., 2010, 154(2): 567-570.

[18] 李梅蘭, 汪俏梅, 朱祝軍, 曾廣文. 春化對白菜DNA甲基化、GA含量及蛋白質(zhì)的影響. 園藝學(xué)報, 2002, 29(4): 353-357.

LI M L, WANG Q M, ZHU Z J, ZENG G W. Studies on the changes of DNA methylation level, GA content and protein in non-heading Chinese cabbage during vernalization., 2002, 29(4): 353-357. (in Chinese)

[19] AHMAD F, HUANG X, LAN H X, HUMA T, BAO Y M, HUANG J, ZHANG H S. Comprehensive gene expression analysis of the DNA (cytosine-5) methyltransferase family in rice (L.)., 2014, 13(3): 5159-5172.

[20] VLASOVA T I, DEMIDENKO Z N, KIRNOS M D, VANYUSHIN B F.DNA methylation by wheat nuclear cytosine DNA methyltransferase: effect of phytohormones., 1995, 157(1/2): 279-281.

[21] NAKANO Y, STEWARD N, SEKINE M, KUSANO T, SANO H. A tobaccocDNA encoding a DNA methyltransferase: molecular characterization and abnormal phenotypes of transgenic tobacco plants., 2000, 41(4): 448-457.

[22] SHARMA R, SINGH R K M, MALIK G, Deveshwar P, TYAGI A K, KAPOOR S, KAPOOR M. Rice cytosine DNA methyltransferases- gene expression profiling during reproductive development and abiotic stress., 2009, 276(21): 6301-6311.

[23] LIN Y T, WEI H M, LU H Y, LEI Y L, FU S F. Developmental- and tissue-specific expression ofencoding a chromomethylase in., 2015, 56(6): 1124-1143.

[24] TANAKA J, ISHIDA T, CHOI B I, YASUDA J, WATANABE T, LWAKURA Y. Latent HIV-1 reactivation in transgenic mice requires cell cycle-dependent demethylation of CREB/ATF sites in the LTR., 2003, 17(2): 167-175.

[25] DOUET J, BLANCHARD B, CUVILLIER C, TOURMENTE S. Interplay of RNA Pol IV and ROS1 during post-embryonic 5S rDNA chromatin remodeling., 2008, 49(12): 1783-1791.

[26] CHOI C S, SANO H. Abiotic-stress induces demethylation and transcriptional activation of a gene encoding a glycerophosphodiesterase- like protein in tobacco plants., 2007, 277(5): 589-600.

[27] PAVET V, QUINTERO C, CECCHINI N M, ROSA A L, ALVAREZ M E.displays centromeric DNA hypomethylation and cytological alterations of heterochromatin upon attack by-, 2006, 19(6): 577-587.

（責(zé)任編輯 岳梅）

Effect of Exogenous Gibberellin on DNA Methylation Level and Expression of Related Enzyme Genes in Tree Peony Floral Buds

ZHANG Tao, SI FuHui, ZHANG YuXi, GAI ShuPeng

(College of Life Sciences, Qingdao Agricultural University/Key Lab of Plant Biotechnology in Universities of Shandong Province, Qingdao 266109, Shandong)

【Objective】 Application of exogenous gibberellin (GA) is a commonly method to break floral bud dormancy in the anti-season production of tree peony. However, the mechanism of how gibberellin plays a role in breaking dormancy of tree peony is an unsolved problem. The objective of this study is to clarify whether GA participates in dormancy regulation through epigenetic mode. 【Method】 In this study, 4-year-old ‘Luhehong’ was used as material, and 500 mg·L-1exogenous GA3was applied. Terminal buds were harvested and DNA was extracted at 0.5, 1, 3 and 5 d after GA3applied. The DNA of peony floral buds was extracted using CTAB method and DNA methylation level of floral buds was analyzed by HPLC technology. Three types of DNA methyltransferase genes,,and, and one DNA demethylase genewere obtained using RACE amplification based on transcriptome analysis. Bioinformatics method was performed to analyze gene sequences and possible domains, and MEGA 5.0 was used to construct phylogenetic tree. Usingas reference gene, the tissue expression characteristics and response to GA3of these genes were analyzed by qPCR method. 【Result】The open reading frame (orf) of the three DNA methyltransferases (Dnmts) varies in length. The orf lengthof,andis 1 118, 1 056 and 2 175 bp, respectively. The number of amino acids encoded is 372, 351 and 724, respectively. All three DNA methyltransferases have methyltransferase domains. The DNA demethylase genehas a long sequence with an ORF of 6 636 bp and encodes 2 211 amino acids, and there is a conserved DNA glycosylation domain in. Phylogenetic analysis indicated that the PsCMT and PsDRM proteins in peony were closely related to the VvCMT2 protein of grape, PsMET and PsROS1 were grouped with most woody plants and located far away from tobacco.,andwere expressed in all tissues (roots, stems, leaves, bracts, sepals, petals, stamens, carpels) at the early flowering stage of peony, and the expression level of,andwas higher in the root of peony than that of other tissues, while the expression level ofwas the highest in carpels. Application of GA3dramatically accelerated buds sprouting and flowering. Bud sprouting could be observed 5 d after GA3application, and the germinating rate of bud was 97.5% with a 92.5% flowering rate after 60 d of GA3application. But in the control group, only a few buds germinated at 20 d, and the germinating rate was only 23.1% at 60 d. GA3significantly decreased the DNA methylation level of peony bud (＜0.01), from 38.9% before treatment to 28.7% at 5 d after treatment. GA3treatment increased the expression ofand, decreased the expression level of, but there was no significant difference inexpression. 【Conclusion】GA3treatment inducesexpression and inhibitsexpression, which leads to DNA hypomethylation and thus promotes dormancy release of peony, which may be an important way for gibberellin to play its role.

; DNA methyltransferase;; gibberellin; DNA methylation; dormancy release

10.3864/j.issn.0578-1752.2018.18.012

2018-04-16；

2018-06-19

國家自然科學(xué)基金面上項目（31372104，31672194）

張濤，E-mail：zhang1995tao@163.com。通信作者蓋樹鵬，E-mail：spgai@qau.edu.cn