DNA甲基化參與調(diào)控馬鈴薯干旱脅迫響應(yīng)

李鵬程 畢真真 梁文君 孫 超 張俊蓮 白江平,*

DNA甲基化參與調(diào)控馬鈴薯干旱脅迫響應(yīng)

李鵬程1,2,**畢真真1,2,**梁文君1,2孫 超1,2張俊蓮1白江平1,2,*

1甘肅省干旱生境作物學(xué)重點實驗室 / 甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點實驗室, 甘肅蘭州 730070;2甘肅農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 甘肅蘭州 730070

非生物脅迫下表觀遺傳對調(diào)控植物基因表達起重要作用, 但是有關(guān)馬鈴薯干旱脅迫下的表觀遺傳研究甚少。本研究以馬鈴薯品種大西洋、費烏瑞它、C119、C16和青薯9號為試驗材料, 以MS培養(yǎng)基為對照以及分別添加200 mmol L-1甘露醇、60 μmol L-1甲基化抑制劑(5-azadC)和60 μmol L-1甲基化抑制劑+200 mmol L-1甘露醇, 處理24 d后對試管苗表型性狀和生理指標進行綜合分析。結(jié)果發(fā)現(xiàn), 不同品種馬鈴薯對甘露醇和甲基化抑制劑響應(yīng)程度趨勢類似。在干旱和DNA甲基化抑制劑分別處理下, 馬鈴薯植株干鮮重、株高、葉片數(shù)和葉綠素含量均顯著減少(<0.05), SOD、POD、CAT活性和Pro、MDA含量均顯著增加(<0.05), 而分枝數(shù)、根長、平均根粗均無明顯變化, 表明馬鈴薯不同性狀指標在響應(yīng)干旱脅迫和DNA去甲基化時, 受到的調(diào)控通路可能不同。進一步比較干旱脅迫和DNA去甲基化共同處理與分別處理下表型性狀和生理指標的差異發(fā)現(xiàn), 共同處理使馬鈴薯植株表型性狀受到進一步抑制, 同時活化了SOD、POD、CAT, 并且使Pro、MDA含量增加, 表明馬鈴薯在響應(yīng)干旱脅迫過程中, 部分表型的形成與DNA甲基化調(diào)控相關(guān)。這將為深入研究馬鈴薯干旱脅迫響應(yīng)與表觀遺傳學(xué)之間的調(diào)控網(wǎng)絡(luò)通路提供初步的理論基礎(chǔ)。

馬鈴薯; 干旱脅迫; DNA甲基化; 形態(tài)生理; 綜合評價

馬鈴薯(L.)是茄科一年生草本塊莖植物, 其栽培范圍遍布全球。中國的馬鈴薯種植區(qū)域主要分布在西北干旱和半干旱地區(qū), 總產(chǎn)量位居全國第一, 現(xiàn)已發(fā)展成為這些地區(qū)的重要支柱產(chǎn)業(yè)[1], 但是其單產(chǎn)卻低于全國平均水平[2], 主要是因為馬鈴薯在塊莖形成及膨大期間的水分缺失嚴重影響了馬鈴薯植株的正常生長發(fā)育[3]。因此, 可通過提高馬鈴薯對水分缺失的耐受性, 進而提高馬鈴薯的產(chǎn)量。

對表觀遺傳機理的解析可以作為作物抗逆遺傳改良的重要途徑之一[4-5]。植物中的DNA甲基化分布具有物種、組織、器官、位點特異性和時空特異性[6-7], 在植物生長發(fā)育過程中具有十分重要的調(diào)控作用[8]。非生物脅迫會引起植物體內(nèi)DNA甲基化水平的改變, 激活脅迫相關(guān)基因的表達、轉(zhuǎn)座子的活性以及調(diào)控啟動子區(qū)域的開關(guān), 以增強植物適應(yīng)及抵御脅迫的能力[9]。潘雅姣等[10]發(fā)現(xiàn), 水稻在干旱脅迫處理下, 降低其葉和根的甲基化程度有利于激活基因的表達, 從而提高水稻耐旱性。研究DNA甲基化常用的方法是通過高分辨率的參考基因組裝配和亞硫酸氫鹽測序繪制基因組中現(xiàn)有的甲基化模式。這種方法富有成效, 但它是描述性的, 而不是試驗性的。另一種方法是通過試驗來改變其甲基化水平[11], 植物中常用5-Aza-2’-脫氧胞苷(5-azadC)來抑制其甲基化程度。5-azadC是胞嘧啶的類似物,不可逆地與DNA甲基轉(zhuǎn)移酶結(jié)合, 使基因組難以維持甲基化狀態(tài), 從而引起DNA的去甲基化[12]。毛銳濤[13]用DNA甲基化抑制劑處理馬鈴薯品種大西洋試管苗, 發(fā)現(xiàn)DNA去甲基化處理嚴重影響其生長發(fā)育。在小麥中, 5-氮雜胞苷通過降低小麥基因組DNA甲基化水平來提高耐鹽性[14]。以上研究表明甲基化變異在抗旱研究中存在潛在的應(yīng)用價值, 可通過發(fā)掘與作物抗旱相關(guān)的表觀遺傳調(diào)控相關(guān)基因, 對作物的抗旱性進行遺傳改良。

目前, 相對于模式植物擬南芥以及水稻、小麥、玉米等主要農(nóng)作物, 國內(nèi)外學(xué)者有關(guān)馬鈴薯DNA甲基化的研究主要集中在離體培養(yǎng)及脅迫環(huán)境下DNA甲基化的變化方面[15-18]。但有關(guān)DNA甲基化在調(diào)控馬鈴薯響應(yīng)干旱脅迫下生長發(fā)育和生理生化的研究卻少見報道。本研究主要通過甲基化抑制劑處理不同干旱敏感的馬鈴薯品種, 并比較在干旱脅迫下去甲基化對馬鈴薯試管苗表型性狀及生理生化的響應(yīng), 為進一步解析馬鈴薯干旱脅迫響應(yīng)的表觀遺傳學(xué)機理提供初步的表型基礎(chǔ)及實驗思路, 為開辟研究馬鈴薯抗旱分子育種的新方法提供一定的理論依據(jù)。

1 材料與方法

1.1 試驗材料

大西洋(Atl)和費烏瑞它(Fav)為干旱敏感型品種; C119 (CIP編號為398098.119)為抗旱中間型材料; C16 (CIP編號為397077.16)和青薯9號(QS9)為抗旱品種。其中C16和C119引自國際馬鈴薯中心(CIP), 其余由甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點實驗室提供。5-azadC (5-aza-2’-deoxycytidine)購于Sigma公司。

1.2 試驗方法

1.2.1 試驗設(shè)計 分別將5個品種的馬鈴薯試管苗莖段接種于4種不同的培養(yǎng)基上。即正常的MS培養(yǎng)基(對照, T0)、含200 mmol L-1甘露醇的MS培養(yǎng)基(干旱處理, T1)、含60 μmol L-15-azadC的MS培養(yǎng)基(去甲基化處理, T2)、含60 μmol L-15-azadC與200 mmol L-1甘露醇的MS培養(yǎng)基(去甲基化與干旱共同處理, T3)。培養(yǎng)條件為溫度(22±2)℃, 光強25.0~37.5 μmol m-2s-1, 光照16 h d-1。處理24 d后, 從每個處理隨機挑選10株較為一致的試管苗分析形態(tài)指標, 同時取3份樣品供生理生化指標的測定。

1.2.2 指標測定 (1)形態(tài)指標測定。株高、葉片數(shù)和分枝數(shù)直接觀察測量, 使用根系掃描儀(EPSON Scan 2 12000XL 2.2)和根系形態(tài)與結(jié)構(gòu)分析儀(WinRHIZO 2017a)測定馬鈴薯試管苗總根長和平均根粗等主要根系參數(shù)。使用佳能單反相機(Canon EOS 5D Mark II)進行植株形態(tài)拍照(圖1)。(2)植株干鮮重測定。將試管苗用無菌水清洗, 吸去表面的水分, 稱量鮮重; 然后將植株放入錫箔紙中于105℃殺青40 min, 85℃烘直至恒重, 稱量干重。(3)生理生化指標測定。葉綠素含量采用80%丙酮浸泡法[19], 超氧化物歧化酶(SOD)活性測定采用氮藍四唑光還原法[20], 過氧化物歧化酶(POD)活性測定采用愈創(chuàng)木酚法[20], 過氧化氫酶(CAT)活性測定采用紫外吸收法[20], 丙二醛(MDA)含量測定采用硫代巴比妥酸法[20], 游離脯氨酸(Pro)含量測定采用酸性茚三酮比色法[20]。所有生理生化指標均3次重復(fù), 取平均值。

1.3 數(shù)據(jù)分析

使用Genstat 19.0進行方差分析、主成分分析和聚類分析; GraphPad Prism 7.0數(shù)據(jù)作圖, 利用隸屬函數(shù)對5個馬鈴薯品種分別進行干旱及去甲基化響應(yīng)性綜合評價。

不同處理下的所有性狀數(shù)據(jù)均采用相對值。

= 1, 2, 3, …,(3)

A~E分別為大西洋、費烏瑞它、C119、C16和青薯9號處理24 d后結(jié)果。T0: CK; T1: 5-azadC; T2: 甘露醇; T3: 5-azadC+甘露醇。每個處理選取3株用于照相。標尺為2 cm。

A–E are the phenotype of Atl, FAV, C119, C16, and QS9 after 24 days of treatment. T0: CK; T1: 5-azadC; T2: mannitol; T3: 5-azadC+mannitol. Three plantlets of each treatment were selected for photography. Bar = 2 cm.

式(4)中,(X)指第個綜合指標的隸屬函數(shù)值,X指第個綜合指標值,max指第個綜合指標的最大值,min指第個綜合指標的最小值。

各綜合指標權(quán)重的計算:

式(5)中,W表示第個綜合指標在所有綜合指標中的重要程度及權(quán)重;P為各品種第個綜合指標的貢獻率。

綜合評價值的計算:

式(6)中,表示不同馬鈴薯品種抗旱性和對5-azadC的響應(yīng)性的綜合評價值。通過計算加權(quán)隸屬函數(shù)值所得值。值越大表明抗旱性越強以及對5-azadC的響應(yīng)性越大。

2 結(jié)果與分析

2.1 干旱脅迫對馬鈴薯試管苗生長的影響

2.1.1 表型性狀和生理生化指標分析 馬鈴薯試管苗在甘露醇模擬干旱脅迫下各表型性狀(圖1)和生理生化特性(圖2)均存在差異, 一般利用相對值來反映不同品種的抗旱性。各個性狀抗旱系數(shù)由高到低依次為丙二醛、過氧化氫酶、游離脯氨酸、過氧化物酶、超氧化物歧化酶、平均根粗、葉綠素、總根長、葉片數(shù)、干重、分枝數(shù)、株高、鮮重, 變化范圍在0.4175~2.3830 (表1)。表明參試馬鈴薯品種不同性狀的抗旱性存在明顯差異。各品種對干旱的綜合響應(yīng)系數(shù)CRC值介于0.9515~1.3897, 馬鈴薯品種的抗旱性由強到弱依次為青薯9號>C16>費烏瑞它>C119>大西洋。

2.1.2 主成分分析 由表2和表3可知, 與抗旱性相關(guān)的前3個主成分的累計貢獻率為98.26%, 其中PC1的特征值是8.595, 貢獻率為66.12%, 其較大的特征向量是分枝數(shù)和葉片數(shù); PC2的特征值是2.46, 貢獻率為18.92%, 載荷較高的是株高和游離脯氨酸; PC3的特征值是1.719, 貢獻率為13.22%, 載荷較高的是丙二醛。這樣可以分別把抗旱系數(shù)的13個單項指標轉(zhuǎn)化為3個相互獨立的綜合指標, 同時也可以代表原始數(shù)據(jù)的大部分信息。

2.1.3 綜合評價 根據(jù)對抗旱性的主成分分析分別得到的3個相互獨立的綜合指標, 利用公式(4)計算不同品種各綜合指標的隸屬函數(shù)值(X)。從表4可以看出, 干旱脅迫處理下, 同一綜合指標CI(1)中, 青薯9號(1)最大, 為1.0000, 表明青薯9號在CI(1)綜合指標上的抗旱性最強。根據(jù)各綜合指標貢獻率, 利用公式(5)計算其權(quán)重W, 計算得到與抗旱性相關(guān)的3個綜合指標權(quán)重分別為0.6729、0.1926、0.1345。進一步利用公式(6)計算抗旱性的綜合評價值, 根據(jù)值, 抗旱性為青薯9號>C16> C119>大西洋>費烏瑞它。

圖2 不同處理下馬鈴薯試管苗表型和生理生化指標的變化

顯著性分析使用Duncan’s法進行多重比較(< 0.05)。Multiple comparisons were performed with Duncan’s method (< 0.05).

表1 不同品種馬鈴薯各指標對甘露醇和5-azadC分別處理的響應(yīng)程度

(續(xù)表1)

FW: fresh weight; DW: dry weight; SH: shoot height; NB: number of branches; NL: number of leaves; RL: root length; ARD: average root diameter; SOD: superoxide dismutase; POD: peroxidase; CAT: catalase; Pro: proline; MDA: malonaldehyde; Chl: chlorophyll; CRC: comprehensive response coefficient.

表2 各綜合指標特征值及貢獻率

表3 各表型因子載荷矩陣

縮寫同表1。Abbreviations are the same as those given in Table 1.

表4 供試品種的綜合性狀指標、權(quán)重、μ(X)及綜合評價值(D)

2.2 DNA甲基化抑制劑對馬鈴薯試管苗生長的影響

2.2.1 表型性狀和生理生化指標分析 馬鈴薯試管苗在5-azadC處理下各表型性狀(圖1)和生理生化特性(圖2)也存在明顯差異, 由表1可知, 各個考察性狀對5-azadC響應(yīng)系數(shù)的變化范圍為0.9902~1.0797, 由高到低依次為丙二醛、超氧化物歧化酶、分枝數(shù)、干重、過氧化氫酶、過氧化物酶、株高、平均根粗、鮮重、葉綠素、葉片數(shù)、游離脯氨酸、總根長。利用公式(3)計算各品種對5-azadC綜合響應(yīng)系數(shù)的CRC值介于0.8753~1.1114。根據(jù)5-azadC綜合響應(yīng)系數(shù)對供試馬鈴薯進行排序, 響應(yīng)程度為青薯9號>C16>費烏瑞它>大西洋>C119。

2.2.2 主成分分析 由表2和表3可知, 與5-azadC響應(yīng)性相關(guān)的前3個主成分的累計貢獻率為92.01%, 其中PC1的特征值是6.750, 貢獻率為51.92%, 對應(yīng)較大的特征向量是SOD、POD、CAT、Pro; PC2的特征值是2.757, 貢獻率為21.20%, 載荷較高的是平均根粗; PC3的特征值是2.456, 貢獻率為18.89%, 載荷較高的是總根長、葉片數(shù)、Chl和株高。這樣可以把對5-azadC響應(yīng)系數(shù)的13個單項指標轉(zhuǎn)化為3個相互獨立的綜合指標, 同時也可以代表原始數(shù)據(jù)的大部分信息。

2.2.3 綜合評價 根據(jù)對5-azadC響應(yīng)性的主成分分析分別得到的3個相互獨立的綜合指標, 利用公式(4)計算不同品種各綜合指標的隸屬函數(shù)值(X)。從表4可以看出, 同一綜合指標CI(1)中, 青薯9號(1)最大, 為1.0000。表明青薯9號在CI(1)綜合指標上對5-azadC的響應(yīng)最弱。根據(jù)各綜合指標貢獻率, 利用公式(5)計算其權(quán)重W, 計算得到與對5-azadC的響應(yīng)性相關(guān)的3個綜合指標權(quán)重分別為0.5643、0.2304、0.2053。進一步利用公式(6)計算對5-azadC的響應(yīng)性的綜合評價值, 并根據(jù)值對5-azadC的響應(yīng)性為青薯9號>C16>C119>費烏瑞它>大西洋。

2.3 干旱脅迫及DNA甲基化抑制劑對馬鈴薯試管苗生長影響的差異比較

從圖2可以看出, 干旱和DNA去甲基化處理下, 不同基因型馬鈴薯試管苗的干鮮重、株高和總根長顯著減少(<0.05), 而青薯9號和費烏瑞它的分枝數(shù)無顯著變化。C16分枝數(shù)在去甲基化處理下無顯著變化, 而干旱處理下顯著減少(<0.05)。青薯9號葉片數(shù)在干旱處理下無顯著變化, 而去甲基化處理下顯著減少(<0.05)。平均根粗除C16在干旱處理下顯著減少(<0.05)外, 其余均無顯著變化。SOD、POD在干旱脅迫下均顯著增加(<0.05), 除C119 SOD無明顯變化外, 而去甲基化處理下青薯9號顯著增加(<0.05), 其余均無差異。CAT在干旱處理下均顯著增加(<0.05), 而去甲基化處理下大西洋、C119和C16無顯著變化, 費烏瑞它和青薯9號顯著增加(<0.05)。Pro在干旱處理下均顯著增加(<0.05), 而去甲基化處理下除大西洋和費烏瑞它無顯著變化外, 其余均顯著增加(<0.05)。MDA在干旱處理下均顯著增加(<0.05), 而去甲基化處理下均無顯著變化。大西洋、費烏瑞它、C16和青薯9號的葉綠素含量在去甲基化處理下顯著增加(<0.05), 而干旱處理下大西洋和費烏瑞它顯著減少(<0.05), C16和青薯9號無差異。C119在2種處理下的葉綠素含量均顯著減少(<0.05)。

從表5可以看出, 干旱處理與5-azadC處理相比, 使馬鈴薯植株干鮮重、株高、葉片數(shù)和葉綠素含量均顯著減少(<0.05), 其中葉片數(shù)和葉綠素含量減少幅度最大, 分別為33.35%和37.00%; SOD、POD、CAT活性和Pro、MDA含量均顯著增加(<0.05), MDA含量增加幅度最大, 為121.01%; 而分枝數(shù)、根長、平均根粗均無明顯變化。

2.4 干旱脅迫及DNA甲基化抑制劑共同處理對馬鈴薯試管苗生長的影響

以上結(jié)果表明, 甲基化抑制劑對馬鈴薯試管苗生長的影響與干旱脅迫的影響非常相似, 為了進一步分析DNA去甲基化與馬鈴薯響應(yīng)干旱脅迫響應(yīng)的關(guān)系, 比較了5-azadC和甘露醇共同處理與分別單獨處理24 d后各馬鈴薯植株形態(tài)和生理生化指標。從圖2可以看出, 與5-azadC和甘露醇單獨處理相比, 共同處理下, 除QS9的干鮮重和株高外, 其余馬鈴薯試管苗的表型性狀均呈下降趨勢(圖1); 而5種基因型馬鈴薯品種的SOD、POD、CAT活性、Pro和MDA含量呈增加趨勢, 其中響應(yīng)5-azadC和甘露醇越強的品種(如大西洋、費烏瑞它), 其SOD、POD、CAT活性以及Pro含量增加幅度越大, 而MDA含量則反之。

通過二因素方差分析對干旱處理/對照處理(T1/T0)、5-azadC處理/對照處理(T2/T0)、共同處理/對照處理(T3/T0)之間的差異性進行分析, 由表5可知, 與5-azadC處理相比, 5-azadC和甘露醇共同處理下, 馬鈴薯試管苗干鮮重、株高、葉片數(shù)、分枝數(shù)、總根長和葉綠素含量均顯著減少(<0.05), 其中分枝數(shù)和葉片數(shù)減少幅度最大, 分別為51.37%和52.41%; SOD、POD、CAT活性和Pro、MDA含量均顯著增加(<0.05), MDA含量增加幅度最大, 為139.82%, 而平均根粗則無明顯變化。與甘露醇處理相比, 兩者共同處理下馬鈴薯試管苗的干鮮重、株高、葉片數(shù)和總根長均顯著減少(<0.05), 其葉片數(shù)減少幅度最大, 為28.39%, SOD、POD活性和葉綠素含量均顯著增加(<0.05), 其中葉綠素含量增加幅度最大, 為46.54%, 而分枝數(shù)、平均根粗、CAT活性和Pro含量無明顯變化。

表5 甘露醇和5-azadC單獨與共同處理下馬鈴薯試管苗各指標比較分析

顯著性分析使用Duncan’s法進行多重比較(< 0.05)。縮寫同表1。

The significance analysis was performed by multiple comparisons using the Duncan’s method (< 0.05). Abbreviations are the same as those given in Table 1.

3 討論

干旱是影響植物正常生長發(fā)育的一個最重要的非生物因素。在長期進化中, 植物形成了系統(tǒng)的抗旱機制[21]。一些研究發(fā)現(xiàn), 植物通過基因組DNA的甲基化、去甲基化過程快速地應(yīng)對外界干旱環(huán)境, 表明DNA甲基化在植物響應(yīng)干旱脅迫中起重要作用[22]。而目前有關(guān)馬鈴薯抗旱性與DNA甲基化之間的研究較少。本研究利用DNA甲基化抑制劑、甘露醇模擬干旱脅迫分別處理不同品種馬鈴薯試管苗, 分析對其形態(tài)和生理生化的影響, 來初步預(yù)測DNA甲基化與馬鈴薯響應(yīng)干旱脅迫之間的關(guān)系。

植物暴露于外源5-azadC, 可通過限制DNA甲基化來誘導(dǎo)表型和發(fā)育性狀變異。許多植物已經(jīng)顯示出這種現(xiàn)象, 例如造成水稻株高降低以及馬鈴薯花和葉片形態(tài)的異常發(fā)育[23-24]。研究發(fā)現(xiàn), 去甲基化處理顯著抑制了馬鈴薯的生長和發(fā)育。植物對不同處理的響應(yīng)程度往往受多種因素影響。為了減少不同品種間各自的遺傳背景差異, 使評價更具有代表性, 往往用多種相關(guān)指標綜合評價及分析[25-27]。本試驗分別選取了馬鈴薯在5-azadC、甘露醇處理下的13個指標用于綜合評價。通過主成分分析可將13個單項指標轉(zhuǎn)換成3個新的相對獨立的綜合指標, 進一步對3個新的綜合指標利用隸屬函數(shù)分析法得到響應(yīng)5-azadC和干旱的綜合評價值。根據(jù)綜合評價值, 5種馬鈴薯品種對干旱和DNA甲基化抑制劑響應(yīng)程度趨勢類似。同時, 對干旱脅迫及DNA甲基化抑制劑差異比較發(fā)現(xiàn), 馬鈴薯植株干鮮重、株高、葉片數(shù)和葉綠素含量均顯著減少(<0.05); SOD、POD、CAT活性和Pro、MDA含量均顯著增加(<0.05), 而分枝數(shù)、根長、平均根粗均無明顯變化。這說明馬鈴薯不同指標對響應(yīng)干旱脅迫和DNA甲基化抑制的調(diào)控通路可能有所差異。

為了進一步研究DNA甲基化與馬鈴薯響應(yīng)干旱脅迫之間的關(guān)系。本研究比較了5-azadC和甘露醇共同處理與分別處理下, 不同馬鈴薯植株的表型性狀和生理特性差異, 發(fā)現(xiàn)雖然不同馬鈴薯品種響應(yīng)5-azadC和甘露醇的強弱不同, 但是共同處理使其表型性狀均受到抑制, 同時活化了SOD、POD、CAT活性, 增加了Pro、MDA的含量。通過二因素方差分析消除品種間差異, 比較不同處理下的表型和生理生化指標差異。與5-azadC、甘露醇單獨處理相比, 二者共同處理下馬鈴薯植株干鮮重、株高、葉片數(shù)、總根長均顯著降低, 而SOD、POD活性均顯著升高, 說明干旱脅迫與去甲基化可能相互作用, 通過不同調(diào)控途徑最終來影響植株的這些表型性狀和生理指標; 而平均根粗無明顯變化, 說明干旱脅迫與去甲基化可能通過相同的途徑調(diào)控這些性狀。在研究5-azaC對苜蓿耐鹽性中發(fā)現(xiàn), 5-azaC可能對植物生長有一些副作用, 導(dǎo)致生長明顯受到抑制, 表現(xiàn)為幼苗株高和干重減少。這表明DNA甲基化對苜蓿的耐鹽性是必要的, 并且與MSAP (methylation sensitive amplification polymorphism, MSAP)的結(jié)果相一致[28]。本研究也發(fā)現(xiàn)DNA去甲基化顯著抑制馬鈴薯的生長發(fā)育。雖然去甲基化影響植物的表型生長, 但植物卻可以同時通過DNA去甲基化來調(diào)控體內(nèi)相關(guān)基因表達、外源基因的防御及基因表達模式的改變等過程, 參與植物對干旱脅迫的響應(yīng), 從而使植物更加快速地應(yīng)對外界干旱環(huán)境[29]。鐘蘭[30]在研究發(fā)現(xiàn), 5-氮雜胞可通過參與活化小麥體內(nèi)響應(yīng)逆境脅迫的抗氧化酶系統(tǒng), 來顯著緩解鹽脅迫對小麥生長所造成的抑制作用。本研究結(jié)果也說明, 5-azadC可能通過DNA去甲基化參與馬鈴薯對干旱脅迫的響應(yīng), 提高馬鈴薯SOD、POD的活性, 最終顯著增強馬鈴薯體內(nèi)響應(yīng)干旱脅迫的抗氧化酶系統(tǒng)活性, 來緩解干旱脅迫對馬鈴薯植株造成的影響。

[1] 陳利洪, 舒幫榮, 李鑫. 中國干旱半干旱地區(qū)農(nóng)牧業(yè)廢棄物資源量評價. 干旱區(qū)資源與環(huán)境, 2018, 32(11): 9–14. Chen L H, Shu B R, Li X. Evaluation of agricultural and animal husbandry waste resources in arid and semi-arid areas of China., 2018, 32(11): 9–14 (in Chinese with English abstract).

[2] 張緒成, 馬一凡, 于顯楓, 侯慧芝, 王紅麗, 方彥杰. 西北半干旱區(qū)深旋松耕作對馬鈴薯水分利用和產(chǎn)量的影響. 應(yīng)用生態(tài)學(xué)報, 2018, 29: 3293–3301. Zhang X C, Ma Y F, Yu X F, Hou H Z, Wang H L, Fang Y J. Effects of vertically rotary sub-soiling tillage on water utilization and yield of potato in semi-arid area of northwest China., 2018, 29: 3293–3301 (in Chinese with English abstract).

[3] Ramírez D A, Yactayo W, Rens L R, Rolando J L, Palacios S, De Mendiburu F, Mares V, Barreda C, Loayza H, Monneveux P, Zotarelli L. Defining biological thresholds associated to plant water status for monitoring water restriction effects: stomatal conductance and photosynthesis recovery as key indicators in potato., 2016, 177: 369–378.

[4] 杜康兮, 沈文輝, 董愛武. 表觀遺傳調(diào)控植物響應(yīng)非生物脅迫的研究進展. 植物學(xué)報, 2018, 53: 581–593. Du K X, Shen W H, Dong A W. Advances in epigenetic regulation of abiotic stress response in plants., 2018, 53: 581–593 (in Chinese with English abstract).

[5] Ma C, Liang B, Chang B, Liu L, Yan J, Yang Y, Zhao Z. Transcriptome profiling reveals transcriptional regulation by DNA methyltransferase inhibitor 5-Aza-2’-deoxycytidine enhancing red pigmentation in bagged “Granny Smith” apples ()., 2018, 19: 3133.

[6] 黃韞宇, 張海軍, 邢燕霞, 齊艷, 孫倩倩, 周春蕾, 趙冰, 郭仰東. NaCl脅迫對黃瓜種子萌發(fā)的影響及DNA甲基化的MSAP分析. 中國農(nóng)業(yè)科學(xué), 2013, 46: 1646–1656. Huang W Y, Zhang H J, Xing Y X, Qi Y, Sun Q Q, Zhou C L, Zhao B, Guo Y D. Effects of NaCl stress on seed germination and DNA methylation status detected by MSAP analysis in cucumber., 2013, 46: 1646–1656 (in Chinese with English abstract).

[7] Kawakatsu T, Huang S S, Jupe F, Sasaki E, Schmitz R J, Urich M A, Castanon R, Nery J R, Barragan C, He Y, Chen H. Epigenomic diversity in a global collection ofaccessions., 2016, 166: 492–505.

[8] Rathore M S, Jha B. DNA methylation and methylation polymorphism in genetically stable in vitro regenerates ofL. using methylation-sensitive AFLP markers., 2016, 178: 1002–1014.

[9] Khan A, Yadav N S, Morgenstern Y, Zemach A, Grafi G. Activation of Tag1 transposable elements indedifferentiating cells and their regulation by CHROMOMETHYLASE 3-mediated CHG methylation., 2016, 1859: 1289–1298.

[10] 潘雅姣, 傅彬英, 王迪, 朱苓華, 黎志康. 水稻干旱脅迫誘導(dǎo)DNA甲基化時空變化特征分析. 中國農(nóng)業(yè)科學(xué), 2009, 42: 3009–3018. Pan Y J, Fu B Y, Wang D, Zhu L H, Li Z K. Spatial and temporal profiling of DNA methylation induced by drought stress in rice., 2009, 42: 3009–3018 (in Chinese with English abstract).

[11] Kottler E J, Vanwallendael A, Franks S J. Experimental treatment with a hypomethylating agent alters life history traits and fitness in., 2018, 2018: 1–10.

[12] 許嶸, 張文州, 陳琳琳, 丁小明, 謝志新, 林陽君, 辛軍. DNA去甲基化藥物作用機制的研究進展. 現(xiàn)代藥物與臨床, 2017, 32(1): 152–157. Xu R, Zhang W Z, Chen L L, Ding X M, Xie Z X, Lin Y J, Xin J. Research progress on action mechanism of DNA demethylating agents., 2017, 32(1): 152–157 (in Chinese with English abstract).

[13] 毛銳濤. 馬鈴薯甲基轉(zhuǎn)移酶和去甲基化酶預(yù)測及表達研究. 內(nèi)蒙古大學(xué)碩士學(xué)位論文, 內(nèi)蒙古呼和浩特, 2014. Mao R T. Phylogenetics and Expression of Methyltransferases and Demethylases in Potato. MS Thesis of Inner Mongolia University, Hohhot, Inner Mongolia, China, 2014 (in Chinese with English abstract).

[14] Zhong L, Xu Y H, Wang J B. The effect of 5-azacytidine on wheat seedlings responses to NaCl stress., 2010, 54: 753–756.

[15] Law R D, Suttle J C. Transient decreases in methylation at 5'-CCGG-3' sequences in potato (L.) meristem DNA during progression of tubers through dormancy precede the resumption of sprout growth., 2003, 51: 437–447.

[16] 李媛媛, 程鵬, 熊興耀, 洪亞輝. 干旱脅迫下馬鈴薯幼苗DNA甲基化變化研究. 中國馬鈴薯, 2012, 26: 11–16. Li Y Y, Cheng P, Xiong X Y, Hong Y H. DNA methylation in potato under drought stress., 2012, 26: 11–16 (in Chinese with English abstract).

[17] 李媛媛. 干旱脅迫下馬鈴薯苗基因組甲基化變化情況及生理生化變化研究. 湖南農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 湖南長沙, 2011. Li Y Y. Studies on Methylation and Physiological Alterations under Drought Stress in Potato. MS Thesis of Hunan Agricultural University, Changsha, China, 2011 (in Chinese with English abstract).

[18] 王芳, 石茹, 王艦. 馬鈴薯莖尖玻璃化法超低溫保存后DNA甲基化的遺傳變異. 分子植物育種, 2013, 11: 351–357. Wang F, Shi R, Wang J. Genetic variation of DNA methylation in potato shoot tips after the cryopreservation by vitrification approach., 2013, 11: 351–357 (in Chinese with English abstract).

[19] Arnon D I. Copper enzymes in isolated chloroplasts. Polyphenoloxidase in Beta vulgaris., 1949, 24: 1–15.

[20] 王學(xué)奎, 黃見良. 植物生理生化實驗原理與技術(shù)(第3版). 北京: 高等教育出版社, 2015. pp 131–287. Wang X K, Huang J L. Principles and Techniques of Plant Physiological Biochemical Experiment, 3rd edn. Beijing: Higher Education Press, 2015. pp 131–287 (in Chinese).

[21] Zhang H, Lang Z, Zhu J K. Dynamics and function of DNA methylation in plants., 2018, 19: 489–506.

[22] 李艷, 錢偉強. 植物中DNA甲基化及去甲基化研究進展. 生命科學(xué), 2017, (3): 302–309. Li Y, Qian W Q. Mechanisms of DNA methylation and demethylation in plants., 2017, (3): 302–309 (in Chinese with English abstract).

[23] Akimoto K, Katakami H, Kim H J, Ogawa E, Sano C M, Wada Y, Sano H. Epigenetic inheritance in rice plants., 2007, 100: 205–217.

[24] Marfil C F, Asurmendi S, Masuelli R W. Changes in micro RNA expression in a wild tuber-bearing Solanum species induced by 5-azacytidine treatment., 2012, 31: 1449–1461.

[25] 郝小琴, 姚鵬鶴, 高崢榮, 吳子愷. 低溫脅迫對微胚乳超甜超高油玉米耐寒性生理生化特性的影響. 作物學(xué)報, 2014, 40: 1470–1484. Hao X Q, Yao P H, Gao Z R, Wu Z K. Effects of low temperature stress on the physiological and biochemical characteristics of cold tolerance in micro-endosperm super sweet and super high oil maize., 2014, 40: 1470–1484 (in Chinese with English abstract).

[26] 徐寧, 陳冰嬬, 王明海, 包淑英, 王桂芳, 郭中校. 綠豆品種資源萌發(fā)期耐堿性鑒定. 作物學(xué)報, 2017, 43: 112–121. Xu N, Chen B R, Wang M H, Bao S Y, Wang G F, Guo Z X. Identification of alkali tolerance of mungbean germplasm resources during germination., 2017, 43: 112–121 (in Chinese with English abstract).

[27] 楊宏羽, 平海濤, 王蒂, 王麗, 劉玉匯, 白江平, 余斌, 張俊蓮. 不同倍性馬鈴薯品種的抗旱性. 中國沙漠, 2016, 36: 1041–1049. Yang H Y, Ping H T, Wang D, Wang L, Liu Y H, Bai J P, Yu B, Zhang J L. Drought resistance of different ploidy potato varieties., 2016, 36: 1041–1049 (in Chinese with English abstract).

[28] Al-Lawati A, Al-Bahry S, Victor R, Al-Lawati A H, Yaish M W. Salt stress alters DNA methylation levels in alfalfa (spp.).2016, 15: 15018299.

[29] 黃玉蘭, 殷奎德, 向君亮, 徐業(yè)梅, 劉兆兵, 豐經(jīng)書, 任海娟. 干旱脅迫下植物生理生化及DNA甲基化的研究進展. 玉米科學(xué), 2016, 24(2): 96–102. Huang Y L, Yin K D, Xiang Y L, Xu Y M, Liu Z B, Feng J S, Ren H J. Research progress on physiological reaction and DNA methylation under drought stress in plant., 2016, 24(2): 96–102 (in Chinese with English abstract).

[30] 鐘蘭. 鹽脅迫下小麥幼苗的生理生化特性及表觀遺傳學(xué)研究. 武漢大學(xué)博士學(xué)位論文, 湖北武漢, 2009. Zhong L. Studies on the Physio-biochemical Properties and Epigenetics of Seedlings of Wheat under Salt Stress. PhD Dissertation of Wuhan University. Wuhan, Hubei, China, 2009 (in Chinese with English abstract).

DNA methylation involved in regulating drought stress response of potato

LI Peng-Cheng1,2,**, BI Zhen-Zhen1,2,**, LIANG Wen-Jun1,2, SUN Chao1,2, ZHANG Jun-Lian1, and BAI Jiang-Ping1,2,*

1Gansu Provincial Key Lab of Aridland Crop Science / Gansu Key Lab of Crop Improvement & Germplasm Enhancement, Lanzhou 730070, Gansu, China;2College of Agronomy, Gansu Agricultural University, Lanzhou 730070, Gansu, China

Although epigenetics is essential for regulating gene expression in plant under abiotic stress, there are few reports regarding the epigenetics of potato under drought stress. In this study, five potato cultivars including Atlantic, Favorita, C119, C16, and Qingshu 9, were used to verify that DNA methylation is involved in potato drought stress response. The plantlets were cultured for 24 days on MS medium, MS medium supplemented with 60 μmol L-1methylation inhibitor (5-azadC), 200 mmol L-1mannitol, and 60 μmol L-1methylation inhibitor combined with 200 mmol L-1mannitol, respectively. Phenotypes, and physiological and biochemical traits were recorded and comprehensively analyzed. The response of five potato cultivars to mannitol and 5-azadC was similar.With treatments of mannitol or 5-azadC, the dry/fresh weight, shoot height, leaf number and chlorophyll content decreased significantly (< 0.05), while the activities of SOD, POD, CAT, and contents of proline and MDA increased significantly (< 0.05). No significant changes were observed in branch number, root length and average root diameter, indicating that potato traits might have different regulatory pathways in response to drought stress and DNA demethylation. Compared with the individual treatment, combined-treatment further inhibited the growth of plantlets, and increased activities of SOD, POD, CAT, and contents of MDA and Pro under drought stress and 5-azadC treatment, indicating that the formation of some phenotypes (not all) in response to drought was regulated through the DNA methylation. The results provide preliminary data for further study on epigenetics regulatory pathway of potato under drought stress.

potato; drought stress; DNA methylation; morphological physiology; comprehensive analysis

本研究由國家自然科學(xué)基金項目(31660432, 31460369), 國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系(馬鈴薯)建設(shè)專項(CARS-09-P14), 甘肅省馬鈴薯產(chǎn)業(yè)體系(GARS-03-P1), 中國科學(xué)院“西部之光”人才培養(yǎng)計劃(2014-01), 蘭州市科技發(fā)展計劃項目(2015-3-62), 甘肅省科技廳項目(18JR3RA170)和甘肅省教育廳高校科研項目(2018A-040)資助。

This study was supported by the National Natural Science Foundation of China (31660432, 31460369), the China Agriculture Research System (CARS-09-P14), the Gansu Potato Industry System (GARS-03-P1), the “Light of the West” Talent Training Program of the Chinese Academy of Sciences (2014-01), the Lanzhou Science and Technology Development Plan (2015-3-62), the Gansu Science and Technology Fund (18JR3RA170), and the Gansu Provincial Department of Education University Research Project (2018A-040).

白江平, E-mail:baijp@gsau.edu.cn

**同等貢獻(Contributed equally to this work)

李鵬程, E-mail:526040572@qq.com; 畢真真, E-mail: bizhen925@sina.com

2019-02-16;

2019-05-12;

2019-06-06.

10.3724/SP.J.1006.2019.94024

URL: http://kns.cnki.net/kcms/detail/11.1809.S.20190605.1800.006.html