溫秀芳,戴亮芳,趙 俊,羅向東,鄧曉娟,張帆濤,謝建坤

(江西師范大學(xué) 生命科學(xué)學(xué)院,南昌 330022)

東鄉(xiāng)野生稻耐冷漸滲系DNA甲基化水平和模式變化研究

溫秀芳,戴亮芳,趙 俊,羅向東*,鄧曉娟,張帆濤,謝建坤

(江西師范大學(xué) 生命科學(xué)學(xué)院,南昌 330022)

根據(jù)水稻全基因組和特定位點的CpG島序列設(shè)計引物,采用McrBC酶酶切DNA,以東鄉(xiāng)野生稻耐冷漸滲系(IL5335和IL5423)及其雙親為試材,研究其基因組和特定位點的DNA甲基化水平和模式變化特征,探討甲基化變異對野生優(yōu)異基因漸滲的影響。結(jié)果顯示:(1)在覆蓋全基因組的83個CpG島中,IL5335和IL5423的基因組甲基化頻率分別為46.6%和53.8%,低于受體親本協(xié)青早B的62.6%;大部分(75.9%～80.7%)受體親本的甲基化模式在兩耐冷漸滲中能穩(wěn)定遺傳,另外一些位點發(fā)生了甲基化模式的改變,主要表現(xiàn)為脫甲基化(13.3%～18.1%)和過甲基化(4.4%～6.0%)。(2)在耐冷QTLs區(qū)間,兩耐冷漸滲系的甲基化水平為13.3%～26.7%,遠(yuǎn)低于受體親本的61.4%;它們在該區(qū)域的甲基化模式變異主要為脫甲基化(33.3%～40.0%),高于全基因組的平均變異率。(3)分析逆轉(zhuǎn)座子Houba和Osr14區(qū)域的51個CpG島發(fā)現(xiàn),耐冷漸滲系在該區(qū)域具有較高頻率的甲基化修飾和較低的甲基化模式變異。研究表明,種間雜交漸滲誘發(fā)了受體親本廣泛的甲基化水平和模式變異,為野生優(yōu)異基因的有效利用提供了新的機遇和挑戰(zhàn)。

東鄉(xiāng)野生稻;耐冷漸滲系;DNA甲基化;McrBC酶

DNA甲基化(DNA methylation)是真核生物DNA中一種可逆且能遺傳的表觀遺傳修飾,廣泛存在于細(xì)菌和動植物基因組中[1]。DNA甲基化通常被認(rèn)為是基因組DNA的一種防御型機制[2]。但事實上,在植物的生長發(fā)育過程中DNA甲基化修飾起著非常重要的作用,它能關(guān)閉或啟動某些基因的轉(zhuǎn)錄表達程序,參與基因的表達調(diào)控、組蛋白修飾、染色體印記、著絲粒和端粒結(jié)構(gòu)的維持等生命活動[3]。因此,DNA甲基化的遺傳與可塑性對生物體維持正常的生命活動至關(guān)重要。

研究表明,外界因素可誘發(fā)植物基因組甲基化變異,如遠(yuǎn)緣雜交、轉(zhuǎn)基因、離體培養(yǎng)、環(huán)境脅迫等[4-5]。前人對擬南芥、棉花、大米草、玉米等種間雜交后代的研究發(fā)現(xiàn),外源物種DNA的導(dǎo)入可誘發(fā)受體物種基因組廣泛的DNA甲基化水平和模式的變化,引發(fā)相關(guān)基因的表達或沉默,進而引起種間雜交后代的基因表達的表型與其雙親具有顯著差異[6-7]。因此,DNA甲基化水平和模式的改變一方面為種間雜交創(chuàng)造新的遺傳變異提供了新的育種途徑和思路,但同時對野生優(yōu)異基因的有效利用提出了新的挑戰(zhàn)。為此,闡明種間雜交后代中甲基化模式的變異的遺傳與分子機制具有重要的理論和實踐意義。

東鄉(xiāng)野生稻是一種全球分布最北的普通野生稻(OryzarufipogonGriff.),蘊藏大量栽培稻品種改良所沒有的優(yōu)異基因,如高產(chǎn)、強耐冷、抗旱等[8]。利用東鄉(xiāng)野生稻創(chuàng)制水稻(O.sativaL.)新種質(zhì)已成為改良栽培稻重要農(nóng)藝性狀和拓寬水稻遺傳基礎(chǔ)的重要課題。此前本課題組通過種間雜交和多次回交、自交獲得了一些攜帶野生稻優(yōu)異基因的抗旱、耐冷、耐低磷漸滲系[8-9]。隨后的研究發(fā)現(xiàn),這些高世代漸滲系中伴隨一些非孟德爾式的表觀遺傳變異,如高世代(BC1F9)強耐冷漸滲系IL5243和IL5335在株型、株高以及谷粒顏色等形態(tài)學(xué)性狀還存在一定變異或分離,其逆轉(zhuǎn)錄酶的序列結(jié)構(gòu)和熒光定量的相對表達量與供體親本有很大的差異[10],減數(shù)分裂也存在落后染色體及雙核仁等不穩(wěn)定現(xiàn)象[8]。為了更全面地了解野生稻優(yōu)異基因在向栽培稻漸滲過程中發(fā)生的遺傳與表觀遺傳變異特征,本研究基于水稻全基因組序列信息,對水稻全基因組的CpG島進行分析,利用CpG兩側(cè)的序列設(shè)計甲基化特異性PCR擴增引物,分析東鄉(xiāng)野生稻耐冷漸滲系的在種間雜交漸滲過程中所發(fā)生甲基化水平和模式變化特征。在此基礎(chǔ)上,基于我們此前鑒定的東鄉(xiāng)野生稻耐冷基因粗定位QTLs結(jié)果[11]和耐冷漸滲系中已發(fā)生逆轉(zhuǎn)錄酶活性變化的轉(zhuǎn)座子結(jié)果[10],對漸滲系基因組的特定區(qū)域的CpG島甲基化水平和模式變化特征作進一步分析,探討DNA甲基化水平和模式變化對外源優(yōu)異基因漸滲的影響,為東鄉(xiāng)野生稻耐冷基因在實際育種中的充分利用提供有益信息。

1 材料和方法

1.1 材 料

本課題組此前以東鄉(xiāng)野生稻(O.rufipogonGriff.)為供體親本,以栽培稻‘協(xié)青早B’(O.sativaL.cv.XieqingzaoB)為受體親本雜交獲得種間雜種F1。種間雜種F1隨后與受體親本進行回交得BC1F1,經(jīng)單粒傳連續(xù)自交得BC1F9代群體,該漸滲群體經(jīng)耐冷鑒定篩選,獲得強耐冷漸滲系IL5335和IL5243[8]。本文以東鄉(xiāng)野生稻耐冷漸滲系(IL5335和IL5423)及其雙親為試材,研究其基因組和特定位點DNA甲基化水平和模式變化特征。

1.2 CpG島甲基化特異擴增引物設(shè)計

利用Rice Genome Annotation Project(http://rice.plantbiology.msu.edu/)中水稻基因組信息,均勻地選擇分布在12條染色體上不同BAC中的DNA序列,結(jié)合我們對耐冷基因粗定位結(jié)果[11],利用程序Methyl Primer Express v1.0在這些序列中尋找CpG島,針對這些富含CG序列的DNA利用Primers 3 Input(version 4.0)設(shè)計了98對甲基化特異性PCR擴增引物,其中15對為耐冷基因QTL定位區(qū)間的引物,其余83條為覆蓋全基因組的甲基化特異性擴增引物。耐冷粗定位區(qū)間為第1號染色體上的RM246-RM1095,4號染色體上的S47M-S410M,12號染色體上的RM519-RM270[11]。另外,我們此前發(fā)現(xiàn)耐冷漸滲系中逆轉(zhuǎn)座子Houba和Osr15已發(fā)生逆轉(zhuǎn)錄酶活性變化[10],為此根據(jù)各染色體上Houba和Osr15轉(zhuǎn)座區(qū)域的DNA序列,設(shè)計特異引物51對。引物序列交于上海英俊公司合成。

1.3 甲基化引物擴增分析

采用CTAB法提取水稻基因組DNA。DNA檢測合格后用McrBC酶于37 ℃酶切(約16 h以上),隨后65 ℃溫育1 h使酶滅活終止酶切,McrBC酶酶切體系50 μL,包含10×NBE buffer 5 μL,100×BSA 0.5 μL,GTP(100 mmol/L)0.5 μL,McrBC酶(10 U/μL)1 μL,DNA 1 μL,滅菌ddH2O 42 μL。酶切產(chǎn)物用ddH2O稀釋10倍備用,并用陽性和陰性引物檢測酶切是否完全。陽性引物是針對Tos17基因設(shè)計的(正向引物:GCTACCCGTTCTTGGACTAT;反向引物:CTGAAATCGGAGCACTGACA);陰性引物是針對不易發(fā)生甲基化的激酶位點設(shè)計的(正向引物:CGACTAAACCACTCCAA-TCATC;反向引物:CCAATCAAAACTTCTCCTGTAA)。酶切時設(shè)空白對照,即酶切體系中McrBC酶用ddH2O代替,其他體系相同[12]。

McrBC酶切完全DNA樣品用于甲基化引物擴增,擴增反應(yīng)總體積為15 μL,包含10×buffer 1.5 μL,MgCl2(25 mmol/L)0.9 μL,dNTP(2.5 mmol/L)1.2 μL,甲基化引物(5 μmol/L)2.0 μL,Taqpolymerase(5 U/μL)0.2 μL,McrBC酶切后的模板DNA(25 ng/μL)2 μL,滅菌水7.2 μL。PCR擴增條件為96 ℃預(yù)變性5 min,96 ℃變性30 s,58 ℃退火30 s,72 ℃延伸45 s,35個循環(huán)后,72 ℃延伸10 min,最后4 ℃保存5 min。PCR產(chǎn)物用1.0%的瓊脂糖凝膠檢測,并利用凝膠成像分析系統(tǒng)進行拍照,每對引物對所有材料的擴增重復(fù)3次。對照組實驗中未擴增出來的材料同樣也進行3次重復(fù)擴增實驗。

1.4 數(shù)據(jù)統(tǒng)計分析

沒有擴增條帶的記為“0”,表示該DNA片段具有甲基化位點;擴增出目的條帶的記為“1”,表示該DNA片段無甲基化位點。CpG甲基化模式變化分為A、B、C三類,其中A類表示受體親本‘協(xié)青早B’的甲基化模式在漸滲后代中被穩(wěn)定遺傳,包括A1、A2兩個亞類;受體親本‘協(xié)青早B’中無甲基化的位點在漸滲后代變?yōu)榧谆稽c的情況定為B類,即發(fā)生了過甲基化變異;C類表示‘協(xié)青早B’為甲基化位點,而漸滲后代中變無甲基化位點,即發(fā)生了脫甲基化變異。

2 結(jié)果與分析

2.1 基因組DNA酶切效果檢測

為了確定McrBC酶對參試材料的基因組DNA酶切是否完全,本研究分別采用陽性和陰性引物對酶切產(chǎn)物進行擴增檢測。陽性引物可擴增基因組最易發(fā)生甲基化且被切開的Tos17基因,如果酶切不完全,材料可擴增出1條700 bp左右的條帶;陰性引物是針對不易發(fā)生甲基化的激酶位點設(shè)計的,可擴增出約200 bp左右的片段。本研究的McrBC酶切產(chǎn)物檢測結(jié)果表明(圖1,A),漸滲系及其雙親的酶切產(chǎn)物的陽性引物擴增均無700 bp的擴增產(chǎn)物,用陰性引物擴增時出現(xiàn)1條200 bp的條帶,表明漸滲系及其雙親的McrBC酶切產(chǎn)物可用于進一步甲基化檢測。

2.2 耐冷漸滲系基因組甲基化水平和模式變化特點

本研究利用覆蓋全基因組12條染色體上的83對特異性甲基化引物,對McrBC酶酶切后的雙親和東鄉(xiāng)野生稻耐冷漸滲系(IL5335和IL5423)進行特異擴增,分析其甲基化水平和模式變化特點。結(jié)果表明,受體親本‘協(xié)青早B’、供體親本東鄉(xiāng)野生稻和兩耐冷漸滲系都具有較高頻率的DNA甲基化(45.4%～62.6%)。其中,受體親本‘協(xié)青早B’的甲基化修飾程度要遠(yuǎn)高出供體親本東鄉(xiāng)野生稻,分別為62.6%和45.4%(表1和圖1,B);耐冷漸滲系(IL5335和IL5243)的甲基化水平介于雙親之間,且與受體親本存在明顯差異,分別為46.6%和53.8%。IL5335的甲基化水平與東鄉(xiāng)野生稻比較接近,兩者無顯著性差異。表明種間雜交漸滲誘發(fā)了受體親本的甲基化水平發(fā)生了顯著變化,而且在不同的漸滲系中甲基化水平變異也存在明顯差異。

進一步分析這83個CpG島甲基化模式的遺傳和變異特征發(fā)現(xiàn),漸滲系的CpG島甲基化模式變異主要分為A、B、C三類(表2),其中A類表示受體親本‘協(xié)青早B’的甲基化模式在漸滲中能穩(wěn)定遺傳,分為A1、A2兩個亞類。A1和A2分別表示未甲基化位點的遺傳和甲基化位點的遺傳模式。B類和C類分別表示耐冷漸滲系中發(fā)生了過甲基化變異和脫甲基化變異。統(tǒng)計分析表明,大部分受體親本‘協(xié)青早B’甲基化模式在兩耐冷漸滲系IL5335和IL5243中能穩(wěn)定遺傳,遺傳率分別為75.9%和80.7%。兩耐冷漸滲系與受體親本相比,分別有20和16個CpG島位點發(fā)生了甲基化模式變異,并主要表現(xiàn)為脫甲基化變異(C類)和過甲基化變異(B類)。其中IL5335和IL5243脫甲基化頻率分別為18.1%和13.3%;過甲基化變化(B類)分別為4.4%和6%。表明種間雜交漸滲誘發(fā)了受體親本的甲基化模式發(fā)生了顯著變異。

圖1 參試材料McrBC酶切效果檢測(A)與DNA甲基化分析部分電泳圖(B)M.Marker;P、N對應(yīng)的泳道分別表示陽性和陰性引物對參試材料的擴增;a～d對應(yīng)泳道的材料分別為栽培稻‘協(xié)青早B’、東鄉(xiāng)野生稻、IL5335和IL5243;0對應(yīng)的泳道表示無擴增條帶,表示該位點DNA含有甲基化,1對應(yīng)的泳道表示擴增出目的條帶,表示該位點無DNA甲基化

表1 耐冷漸滲系及其雙親基因組DNA甲基化水平特征

注:XB.協(xié)青早B;DWR.東鄉(xiāng)野生稻;Mean為平均值;同列不同字母表示0.05水平顯著性差異;下同。

Note:XB.OryzasativaL.cv.Xieqingzao B;DWR.Dongxiang wild rice.Mean:average values of the different material.The different normal letter within the same column represents significant difference at 0.05 level.The same as below.

表2 耐冷漸滲系基因組DNA甲基化模式的遺傳與變異特點

注:XB(-)和ILs(-)分別表示在XB和ILs中未擴增出目的條帶,表示在相應(yīng)的位點中發(fā)生了DNA甲基化;XB(+)和ILs(+)分別表示在XB和ILs中能擴增出目的條帶,表示在相應(yīng)的位點未發(fā)生DNA甲基化;每列括號前的數(shù)據(jù)表示位點數(shù),括號內(nèi)表示其頻率;下同。

Note:XB(-) and ILs(-) represent no band had been amplified in XB and ILs,respectively,which means that it’s DNA methylation site;XB(+) and ILs(+) represent bands could be amplified in XB and ILs,respectively,which means that it’s non-methylation site.The data before brackets was the number of locus,and in the brackets was its frequency;The same as below.

2.3 耐冷粗定位區(qū)間DNA甲基化變化特點

基于我們前期東鄉(xiāng)野生稻耐冷QTLs的粗定位結(jié)果,對漸滲系的耐冷粗定位區(qū)間的CpG進行甲基化水平和模式特征分析。結(jié)果發(fā)現(xiàn),受體親本‘協(xié)青早B’在耐冷粗定位區(qū)間的甲基化水平達53.3%,與全基因組CpG島甲基化水平(61.4%)無顯著性差異。東鄉(xiāng)野生稻和兩耐冷漸滲系(IL5335和IL5243)在該區(qū)間的甲基化水平分別為20.0%、13.3%和26.7%,低于其全基因組CpG島甲基化的平均水平(表1和表3)。受體親本‘協(xié)青早B’的甲基化模式(A1和A2)大部分在兩耐冷漸滲系中能得到較穩(wěn)定的遺傳,分別為70.0%和60.0%。與受體親本相比,IL5335和IL5243耐冷漸滲系中也發(fā)現(xiàn)一定比例的甲基化遺傳模式變異,并主要表現(xiàn)為脫甲基化(C類),分別為40%和33.3%,高于全基因甲基化模式的變異的平均水平(表1和表3)。

2.4 逆轉(zhuǎn)座子區(qū)域DNA甲基化變化特點

基于我們此前的研究結(jié)果,針對耐冷漸滲系(IL5335和IL5243)中已發(fā)生逆轉(zhuǎn)錄酶活性變化的逆轉(zhuǎn)座子Houba和Osr15區(qū)域設(shè)計了51對甲基化引物,進一步研究雙親及耐冷漸滲系在特定相關(guān)區(qū)域的甲基化水平和模式變化特征。研究發(fā)現(xiàn),在Houba和Osr15逆轉(zhuǎn)座子區(qū)域,‘協(xié)青早B’、東鄉(xiāng)野生稻、IL5335和IL5243的甲基化水平分別為62.7%、56.9%、66.7%和49.0%,均高于全基因組DNA甲基化的平均水平(表1和表4)。在Houba和Osr15區(qū)域,‘協(xié)青早B’、IL5335和IL5243之間的甲基化水平無顯著差異,但均與東鄉(xiāng)野生稻之間存在顯著差異。IL5335和IL5243甲基化模式大部分繼承(86.6%)大部分與受體親本’協(xié)青早B’一致,并主要表現(xiàn)為A1類型的遺傳模式。在轉(zhuǎn)座子區(qū)域,IL5335和IL5243的甲基化模式總變異率(包括過甲基化和脫甲基化)分別為15.7%和9.8%,均低于其他區(qū)域甲基化模式的變異率。表明耐冷漸滲系在Houba和Osr14轉(zhuǎn)座子區(qū)域具有較高頻率的甲基化修飾和較低的甲基化模式變異率,這利于轉(zhuǎn)座子和基因組結(jié)構(gòu)的穩(wěn)定,為漸滲系優(yōu)異性狀的轉(zhuǎn)移和性狀穩(wěn)定奠定基礎(chǔ)。

表3 耐冷漸滲系在耐冷區(qū)間DNA甲基化水平和模式的遺傳與變異特點

表4 耐冷漸滲系在逆轉(zhuǎn)座子區(qū)域DNA甲基化水平和模式的遺傳與變異特點

3 討 論

外源DNA漸滲會給受體基因組帶來一種強烈的“沖擊”(Genome shock),這種沖擊不僅表現(xiàn)在漸滲系后代中會產(chǎn)生各種基于基因組序列變化的遺傳學(xué)變異,同時還可能發(fā)生許多表觀遺傳學(xué)變異[13]。遠(yuǎn)緣雜交育種實踐告訴我們,表觀遺傳學(xué)變異在許多作物新品種改良創(chuàng)造了很多新的抗性和新表型,為作物遺傳育種提供了新途徑,因而對這些變異發(fā)生的分子機制深入研究將有助于作物品種改良[6,14]。作為表觀遺傳修飾的重要內(nèi)容,DNA甲基化可在阻斷遺傳信息傳遞而不改變DNA序列的情況下引起生物體形態(tài)性狀的變化。闡明遠(yuǎn)緣雜交漸滲過程中基因組的甲基化水平和模式變異特征及分子機制可為物種的多樣性維持機制和分子進化機制提供理論基礎(chǔ),同時也對新品種的培育提供有益信息,相關(guān)研究已成為表觀遺傳學(xué)研究的熱點[15-16]。

McrBC是一種甲基化內(nèi)切酶,可特異切割含有甲基化胞嘧啶的DNA單鏈或雙鏈,但不作用于非甲基化的DNA。結(jié)合PCR技術(shù),此方法不僅可以檢測絕大部分甲基化的胞嘧啶(如5-甲基胞嘧啶,5-羥甲基胞嘧啶或N4-甲基胞嘧啶),還可針對特異位點設(shè)計引物,研究特定染色體區(qū)域的甲基化特征,是檢測胞嘧啶甲基化的有利工具[17]。本研究利用該技術(shù)分析耐冷漸滲系全基因組83個CpG島和66個特定區(qū)域CpG島(15耐冷區(qū)間和51轉(zhuǎn)座子區(qū)間)發(fā)現(xiàn),兩耐冷漸滲后代穩(wěn)定遺傳了回交親本的大部分甲基化水平和模式,但也有一部分位點發(fā)生了DNA甲基化模式的變異,包括過甲基化和脫甲基化,但不同區(qū)域CpG位點甲基化水平和變異模式有明顯差異,這與水稻[14-15]、擬南芥[18]、大米草[19]以及玉米[20]等植物雜交后代的甲基化研究相似。與其他DNA甲基化檢測方法(如甲基化敏感多態(tài)性,MSAP)相比,McrBC酶切法檢測得到的DNA甲基化頻率略低一些[4,21],但該技術(shù)技術(shù)穩(wěn)定、重復(fù)性好,能精確研究特定位點、特定基因的甲基化水平和模式變異特征,是研究DNA甲基化水平和模式變化對外源優(yōu)異基因漸滲影響的有效手段。

通常認(rèn)為基因處于甲基化修飾狀態(tài)時表達抑制,而處于非甲基化狀態(tài)時表達活躍,但無論基因處于何種狀態(tài),都可能使?jié)u滲后代產(chǎn)生相關(guān)有利性狀的增強或者產(chǎn)生親本沒有的有利新表型,這取決于特異位點的甲基化增強或降低對這些性狀的影響[6-7]。本研究結(jié)果表明,兩耐冷漸滲系與受體親本‘協(xié)青早B’相比,漸滲后代甲基化變異包括脫甲基化和過甲基化兩種方式,并以脫甲基化為主,東鄉(xiāng)野生稻和兩耐冷漸滲后代IL5335和IL5243在耐冷QTLs粗定位區(qū)間的甲基化水平明顯低于普通區(qū)間的CpG島甲基化水平(分別為20.0%、13.3%和26.7%),并且發(fā)生甲基化模式變化的位點中有53.8%是在該區(qū)域,推測優(yōu)異漸滲系在其形成過程中發(fā)生了更為廣泛的DNA甲基化模式變化,但耐冷QTLs區(qū)間的優(yōu)異基因在雜交漸滲過程中的DNA甲基化變異特征和分子機制還需進一步深入的研究。另外,我們此前發(fā)現(xiàn)漸滲系中Houba和Osr14轉(zhuǎn)座子的活性較高,而本研究的結(jié)果表明,各參試材料在Houba和Osr14轉(zhuǎn)座子區(qū)域的甲基化水平較其他區(qū)域要略高,這是否與DNA甲基化沉默轉(zhuǎn)座子的活性、促進漸滲系性狀穩(wěn)定有關(guān)也還有待于進一步研究。

[1] DONG Z Y,WANG Y M,ZHANG Z J,etal.Extent and pattern of DNA methylation alteration in rice lines derived from introgressive hybridization of rice andZizanialatifoliaGriseb[J].TheoreticalandAppliedGenetics,2006,113:196-205.

[2] GOLL M G,BESTOR T H.Eukaryotic cytosine methyltransferases[J].AnnualReviewofBiochemistry,2005,74:481-51.

[3] NYBOM H,WEISING K,ROTTER B.DNA fingerprinting in botany:past,present,future[J].InvestigativeGenetics,2014,5:1-35.

[4] CHEN L Z,CHEN J F.Changes of cytosine methylation induced by wide hybridization and allopolyploid in cucumis[J].Genome,2008,51(10):789-799.

[5] LI L Q(李麗琴),FU CH H(付春華),ZHAO CH F(趙春芳),etal.Genetic and epigenetic variations ofTaxusmediaL.cv.hicksii during process of dedifferentiation[J].PlantPhysiologyCommunications(植物生理學(xué)通訊),2009,45(6):544-548(in Chinese).

[6] SANETOMO R,HOSAKA K.Reciprocal differences in DNA sequence and methylation status of the pollen DNA between F1hybrids ofSolanumtuberosum×S.demissum[J].Euphytica,2011,182:219-229.

[7] ABID G,MUHOVISKI Y,JACQUEMIN J M,etal.Changes in DNA-methylation during zygotic embryogenesis in interspecific hybrids of beans(Phaseolusssp.)[J].PlantCellTissOrganCulture,2011,105:383-393.

[8] LUO X D,JIAN S R,CAO J F,etal.Identification and molecular cytology analysis of cold tolerance introgression lines derived fromOryzasativaL.mating withO.rufipogonGriff[J].Euphytica,2012,187:461-469.

[9] JIAN SH R(簡水溶),WAN Y(萬 勇),LUO X D(羅向東),etal.Genetic analysis of cold tolerance at the seedling stage in Dongxiang wild rice(Oryzarufipogon)[J].ChineseBulletinofBotany(植物學(xué)報),2011,46(1):21-27(in Chinese).

[10] CHEN Y L(陳雅玲),LUO X D(羅向東),ZHANG F T(張帆濤),etal.Cloning and expression analysis of retrotransposon reverse transcriptase in introgression lines from Dongxiang wild rice[J].ChineseBulletinofBotany(植物學(xué)報),2013,48(2):138-144(in Chinese).

[11] 簡水溶.東鄉(xiāng)野生稻耐冷性狀的QTL定位[D].南昌:江西師范大學(xué),2011.

[12] SCHUMACHER A,KAPRANOV P,KAMINSKY Z,etal.Microarray-based DNA methylation profiling:technology and applications[J].NucleicAcidsResearch,2006,34:528-542.

[13] DENG X J,LUO X D,CHEN Y L,etal.Genetic diversity and genetic changes in the introgression lines derived fromOryzasativaL.mating withO.rufipogonGriff[J].JournalofIntegrativeAgriculture,2012,11(7):1 059-1 066.

[14] YASUDA K,TSUKIYAMA T,KARKI S,etal.Mobilization of the active transposon mPing in interspecific hybrid rice betweenOryzasativaandO.glaberrima[J].Euphytica,2013,192:17-24.

[15] KANTAMA L,JUNBUATHONG S,SAKULKOO J,etal.Epigenetic changes and transposon reactivation in Thai rice hybrids[J].MolecularBreeding,2013,31:815-827.

[16] WANG W S,PAN Y J,ZHAO X Q,etal.Drought-induced site-specific DNA methylation and its association with drought tolerance in rice(OryzasativaL.)[J].JournalofExperimentalBotany,2011,62:1 951-1 960.

[17] ZILBERMAN D,HENIKOFF S.Genome-wide analysis of DNA methylation patterns[J].Development,2007,134:3 959-3 965.

[18] MADLUNG A,MASUELLI R W,WATSON B,etal.Remodeling of DNA methylation and phenotypic and transcriptional changes in syntheticArabidopsisallotetraploids[J].PlantPhysiology,2002,129:733-746.

[19] SALMON A,AINOUCHE M L,WENDEL J F.Genetic and epigenetic consequences of recent hybridization and polyploidy inSpartina(Poaceae)[J].MolecularEcology,2005,14:1 163-1 175.

[20] SEKHON R S,CHOPRA S.Progressive loss of DNA methylation releases epigenetic gene silencing from a tandemly repeated maizeMybgene[J].Genetics,2009,181:81-91.

[21] DONG Z Y,WANG Y M,ZHANG Z J,etal.Extent and pattern of DNA methylation alteration in rice lines derived from introgressive hybridization of rice andZizanialatifoliaGriseb[J].TheoreticalandAppliedGenetics,2006,113:196-205.

(編輯:宋亞珍)

DNA Methylation Level and Pattern Changes in Cold ToleranceIntrogression Lines Derived fromOryzarufipogonGriff.

WEN Xiufang,DAI Liangfang,ZHAO Jun,LUO Xiangdong*,DENG Xiaojuan,ZHANG Fantao,XIE Jiankun

(College of Life Science,Jiangxi Normal University,Nanchang 330022,China)

In the present study,DNA methylation specific PCR primers were designed according to the sequences of CpG island region throughout 12 chromosomes and the specific genomic position.In order to investigate the changes of DNA methylation level and pattern in the cold tolerance introgression lines (ILs,IL5335 and IL5423) and their parents,we used these special primers to amplify the previously digested genomic DNA by McrBC enzyme.(1)Analysis of 83 pairs of primers distributed across the 12 chromosomes of rice showed that the frequency of DNA methylation in IL5335 and IL5423 was to 46.6% and 53.8%,respectively,which was lower than that of recurrent parent (OryzasativaL.cv.Xieqingzao B).Most of DNA methylation patterns of recurrent parent could be inherited in the ILs.At the same time,some of methylation patterns have changed including demethylation (13.3-18.1%) and hypermethylation (4.4%-6.0%).(2)In the QTLs region related to cold tolerance,the frequency of DNA methylation in IL5335 and IL5423 was 13.3%-26.7%,which was much lower than that of recurrent parent.The frequency of DNA methylation patterns changes of the two ILs in the QTLs region was mainly for demethylation (33.3%-40.0%),which was higher than that of the average variation frequency of whole-genome.(3)Further analysis on 51 pairs of primers based on two retrotransposons (HoubaandOsr14) demonstrated that there was higher level of frequency of DNA methylation and lower frequency of methylation pattern changes.These results suggested that the extensive DNA methylation level and pattern changes have been happened during the interspecific hybridization and alien gene introgression,which would provide useful information for efficiently exploiting and using the interesting wild genes for rice enhancement.

Dongxiang wild rice;cold tolerance introgression lines;DNA methylation;McrBC enzyme

1000-4025(2015)04-0694-07

10.7606/j.issn.1000-4025.2015.04.0694

2014-12-27;修改稿收到日期:2015-01-25

國家自然科學(xué)基金(31260255);江西省教育廳科技計劃項目(GJJ12184);江西省科技支撐計劃(20122BBF60064);江西省青年科學(xué)家培養(yǎng)計劃(20112BCB23007)

溫秀芳(1989-),女,在讀碩士研究生,主要從事植物遺傳與分子生物學(xué)研究。E-mail:1025108610@qq.com

*通信作者:羅向東,博士,教授,主要從事植物遺傳與分子生物學(xué)研究。E-mail:xdluolf@163.com

Q789

A