• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2015-02-28 01:58:13溫秀芳戴亮芳羅向東鄧曉娟張帆濤謝建坤
    西北植物學(xué)報(bào) 2015年4期
    關(guān)鍵詞:東鄉(xiāng)親本甲基化

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

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

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

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

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

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

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

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

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

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

    1 材料和方法

    1.1 材 料

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

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

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

    1.3 甲基化引物擴(kuò)增分析

    采用CTAB法提取水稻基因組DNA。DNA檢測(cè)合格后用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倍備用,并用陽(yáng)性和陰性引物檢測(cè)酶切是否完全。陽(yáng)性引物是針對(duì)Tos17基因設(shè)計(jì)的(正向引物:GCTACCCGTTCTTGGACTAT;反向引物:CTGAAATCGGAGCACTGACA);陰性引物是針對(duì)不易發(fā)生甲基化的激酶位點(diǎn)設(shè)計(jì)的(正向引物:CGACTAAACCACTCCAA-TCATC;反向引物:CCAATCAAAACTTCTCCTGTAA)。酶切時(shí)設(shè)空白對(duì)照,即酶切體系中McrBC酶用ddH2O代替,其他體系相同[12]。

    McrBC酶切完全DNA樣品用于甲基化引物擴(kuò)增,擴(kuò)增反應(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擴(kuò)增條件為96 ℃預(yù)變性5 min,96 ℃變性30 s,58 ℃退火30 s,72 ℃延伸45 s,35個(gè)循環(huán)后,72 ℃延伸10 min,最后4 ℃保存5 min。PCR產(chǎn)物用1.0%的瓊脂糖凝膠檢測(cè),并利用凝膠成像分析系統(tǒng)進(jìn)行拍照,每對(duì)引物對(duì)所有材料的擴(kuò)增重復(fù)3次。對(duì)照組實(shí)驗(yàn)中未擴(kuò)增出來(lái)的材料同樣也進(jìn)行3次重復(fù)擴(kuò)增實(shí)驗(yàn)。

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

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

    2 結(jié)果與分析

    2.1 基因組DNA酶切效果檢測(cè)

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

    2.2 耐冷漸滲系基因組甲基化水平和模式變化特點(diǎn)

    本研究利用覆蓋全基因組12條染色體上的83對(duì)特異性甲基化引物,對(duì)McrBC酶酶切后的雙親和東鄉(xiāng)野生稻耐冷漸滲系(IL5335和IL5423)進(jìn)行特異擴(kuò)增,分析其甲基化水平和模式變化特點(diǎn)。結(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)野生稻比較接近,兩者無(wú)顯著性差異。表明種間雜交漸滲誘發(fā)了受體親本的甲基化水平發(fā)生了顯著變化,而且在不同的漸滲系中甲基化水平變異也存在明顯差異。

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

    圖1 參試材料McrBC酶切效果檢測(cè)(A)與DNA甲基化分析部分電泳圖(B)M.Marker;P、N對(duì)應(yīng)的泳道分別表示陽(yáng)性和陰性引物對(duì)參試材料的擴(kuò)增;a~d對(duì)應(yīng)泳道的材料分別為栽培稻‘協(xié)青早B’、東鄉(xiāng)野生稻、IL5335和IL5243;0對(duì)應(yīng)的泳道表示無(wú)擴(kuò)增條帶,表示該位點(diǎn)DNA含有甲基化,1對(duì)應(yīng)的泳道表示擴(kuò)增出目的條帶,表示該位點(diǎn)無(wú)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甲基化模式的遺傳與變異特點(diǎn)

    注:XB(-)和ILs(-)分別表示在XB和ILs中未擴(kuò)增出目的條帶,表示在相應(yīng)的位點(diǎn)中發(fā)生了DNA甲基化;XB(+)和ILs(+)分別表示在XB和ILs中能擴(kuò)增出目的條帶,表示在相應(yīng)的位點(diǎn)未發(fā)生DNA甲基化;每列括號(hào)前的數(shù)據(jù)表示位點(diǎn)數(shù),括號(hào)內(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甲基化變化特點(diǎn)

    基于我們前期東鄉(xiāng)野生稻耐冷QTLs的粗定位結(jié)果,對(duì)漸滲系的耐冷粗定位區(qū)間的CpG進(jìn)行甲基化水平和模式特征分析。結(jié)果發(fā)現(xiàn),受體親本‘協(xié)青早B’在耐冷粗定位區(qū)間的甲基化水平達(dá)53.3%,與全基因組CpG島甲基化水平(61.4%)無(wú)顯著性差異。東鄉(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甲基化變化特點(diǎn)

    基于我們此前的研究結(jié)果,針對(duì)耐冷漸滲系(IL5335和IL5243)中已發(fā)生逆轉(zhuǎn)錄酶活性變化的逆轉(zhuǎn)座子Houba和Osr15區(qū)域設(shè)計(jì)了51對(duì)甲基化引物,進(jìn)一步研究雙親及耐冷漸滲系在特定相關(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之間的甲基化水平無(wú)顯著差異,但均與東鄉(xiāng)野生稻之間存在顯著差異。IL5335和IL5243甲基化模式大部分繼承(86.6%)大部分與受體親本’協(xié)青早B’一致,并主要表現(xiàn)為A1類型的遺傳模式。在轉(zhuǎn)座子區(qū)域,IL5335和IL5243的甲基化模式總變異率(包括過(guò)甲基化和脫甲基化)分別為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甲基化水平和模式的遺傳與變異特點(diǎn)

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

    3 討 論

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

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

    通常認(rèn)為基因處于甲基化修飾狀態(tài)時(shí)表達(dá)抑制,而處于非甲基化狀態(tài)時(shí)表達(dá)活躍,但無(wú)論基因處于何種狀態(tài),都可能使?jié)u滲后代產(chǎn)生相關(guān)有利性狀的增強(qiáng)或者產(chǎn)生親本沒(méi)有的有利新表型,這取決于特異位點(diǎn)的甲基化增強(qiáng)或降低對(duì)這些性狀的影響[6-7]。本研究結(jié)果表明,兩耐冷漸滲系與受體親本‘協(xié)青早B’相比,漸滲后代甲基化變異包括脫甲基化和過(guò)甲基化兩種方式,并以脫甲基化為主,東鄉(xiāng)野生稻和兩耐冷漸滲后代IL5335和IL5243在耐冷QTLs粗定位區(qū)間的甲基化水平明顯低于普通區(qū)間的CpG島甲基化水平(分別為20.0%、13.3%和26.7%),并且發(fā)生甲基化模式變化的位點(diǎn)中有53.8%是在該區(qū)域,推測(cè)優(yōu)異漸滲系在其形成過(guò)程中發(fā)生了更為廣泛的DNA甲基化模式變化,但耐冷QTLs區(qū)間的優(yōu)異基因在雜交漸滲過(guò)程中的DNA甲基化變異特征和分子機(jī)制還需進(jìn)一步深入的研究。另外,我們此前發(fā)現(xiàn)漸滲系中Houba和Osr14轉(zhuǎn)座子的活性較高,而本研究的結(jié)果表明,各參試材料在Houba和Osr14轉(zhuǎn)座子區(qū)域的甲基化水平較其他區(qū)域要略高,這是否與DNA甲基化沉默轉(zhuǎn)座子的活性、促進(jìn)漸滲系性狀穩(wěn)定有關(guān)也還有待于進(jì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(簡(jiǎn)水溶),WAN Y(萬(wàn) 勇),LUO X D(羅向東),etal.Genetic analysis of cold tolerance at the seedling stage in Dongxiang wild rice(Oryzarufipogon)[J].ChineseBulletinofBotany(植物學(xué)報(bào)),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é)報(bào)),2013,48(2):138-144(in Chinese).

    [11] 簡(jiǎn)水溶.東鄉(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

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

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

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

    Q789

    A

    猜你喜歡
    東鄉(xiāng)親本甲基化
    甘蔗親本農(nóng)藝性狀評(píng)價(jià)與分析
    東鄉(xiāng)五月
    東鄉(xiāng)記憶
    金山(2020年10期)2020-12-07 05:59:37
    一個(gè)我從未抵達(dá)過(guò)的地方——東鄉(xiāng)
    快樂(lè)東鄉(xiāng)娃
    幾種蘋(píng)果砧木實(shí)生后代與親本性狀的相關(guān)性
    云瑞10系列生產(chǎn)性創(chuàng)新親本2種方法評(píng)價(jià)
    鼻咽癌組織中SYK基因啟動(dòng)子區(qū)的甲基化分析
    胃癌DNA甲基化研究進(jìn)展
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    日韩av在线免费看完整版不卡| 久久久久久久久大av| 久久久久久久久久成人| 99久久精品一区二区三区| 美女被艹到高潮喷水动态| 亚洲久久久久久中文字幕| 赤兔流量卡办理| av福利片在线观看| eeuss影院久久| 综合色av麻豆| 老司机影院成人| 又爽又黄a免费视频| 成人高潮视频无遮挡免费网站| 91午夜精品亚洲一区二区三区| 成人国产av品久久久| 下体分泌物呈黄色| 久热久热在线精品观看| 国产探花极品一区二区| 免费看a级黄色片| 高清av免费在线| 嫩草影院入口| 国产乱人偷精品视频| 九草在线视频观看| 国内揄拍国产精品人妻在线| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 中国国产av一级| 国产精品久久久久久av不卡| 内地一区二区视频在线| 国产精品一区二区三区四区免费观看| 国产黄频视频在线观看| 噜噜噜噜噜久久久久久91| 人妻少妇偷人精品九色| 亚洲精华国产精华液的使用体验| 只有这里有精品99| 国产精品女同一区二区软件| 三级国产精品片| 欧美97在线视频| 在线天堂最新版资源| 亚洲怡红院男人天堂| 免费观看性生交大片5| 国产精品伦人一区二区| 国产成人freesex在线| 欧美日韩在线观看h| 久久影院123| 欧美xxxx性猛交bbbb| av线在线观看网站| 一级黄片播放器| 一级片'在线观看视频| av在线蜜桃| 99视频精品全部免费 在线| 精品99又大又爽又粗少妇毛片| av在线播放精品| 美女xxoo啪啪120秒动态图| 日本av手机在线免费观看| 精品久久久久久久末码| 亚洲国产精品成人久久小说| 在线亚洲精品国产二区图片欧美 | 国产精品国产三级国产av玫瑰| 美女内射精品一级片tv| 免费看日本二区| 久久久久久久久久成人| 精品一区二区免费观看| 丝袜脚勾引网站| 欧美潮喷喷水| 午夜福利高清视频| 日韩成人av中文字幕在线观看| 亚洲国产最新在线播放| 久久精品熟女亚洲av麻豆精品| 成人无遮挡网站| 在现免费观看毛片| 欧美丝袜亚洲另类| 最近手机中文字幕大全| av国产久精品久网站免费入址| 亚洲国产精品成人久久小说| 特大巨黑吊av在线直播| 亚洲不卡免费看| 赤兔流量卡办理| 亚洲欧美成人精品一区二区| 国产精品蜜桃在线观看| 免费看光身美女| 亚洲精品日本国产第一区| 五月开心婷婷网| 亚洲精品乱码久久久v下载方式| 在线播放无遮挡| 国产亚洲av片在线观看秒播厂| 97精品久久久久久久久久精品| 一本久久精品| 男女那种视频在线观看| 精品99又大又爽又粗少妇毛片| a级毛色黄片| 日日啪夜夜爽| 内射极品少妇av片p| 纵有疾风起免费观看全集完整版| 亚洲av.av天堂| 国产成人福利小说| 黄色视频在线播放观看不卡| 免费看光身美女| 26uuu在线亚洲综合色| 黄色怎么调成土黄色| 狠狠精品人妻久久久久久综合| 五月开心婷婷网| 亚洲av男天堂| 熟女人妻精品中文字幕| 成人高潮视频无遮挡免费网站| 少妇人妻 视频| videossex国产| 免费少妇av软件| 91久久精品电影网| 舔av片在线| 亚洲成人一二三区av| 日韩,欧美,国产一区二区三区| 少妇丰满av| 一区二区三区乱码不卡18| 最近最新中文字幕大全电影3| 全区人妻精品视频| 夫妻午夜视频| 欧美日韩亚洲高清精品| 亚洲人成网站高清观看| 国产亚洲av片在线观看秒播厂| 日本猛色少妇xxxxx猛交久久| 欧美激情在线99| 韩国高清视频一区二区三区| 精品人妻一区二区三区麻豆| 午夜免费观看性视频| 人人妻人人澡人人爽人人夜夜| 2022亚洲国产成人精品| 熟女电影av网| 亚洲精品乱码久久久久久按摩| 精品一区二区免费观看| 视频中文字幕在线观看| 亚洲欧美一区二区三区国产| 久久久久久久大尺度免费视频| 亚洲色图av天堂| 亚洲精品aⅴ在线观看| 国产午夜精品久久久久久一区二区三区| 欧美日韩视频高清一区二区三区二| 如何舔出高潮| 国产精品一区二区性色av| 一级毛片电影观看| 亚洲精品国产av蜜桃| 99精国产麻豆久久婷婷| 精品一区二区免费观看| 少妇猛男粗大的猛烈进出视频 | 色网站视频免费| www.色视频.com| 亚洲av免费高清在线观看| 国产伦精品一区二区三区四那| 亚洲精品视频女| 欧美zozozo另类| 乱码一卡2卡4卡精品| 国产女主播在线喷水免费视频网站| 99九九线精品视频在线观看视频| 国产 一区 欧美 日韩| 两个人的视频大全免费| 女人十人毛片免费观看3o分钟| av黄色大香蕉| 国产亚洲5aaaaa淫片| 国精品久久久久久国模美| 中文字幕人妻熟人妻熟丝袜美| 两个人的视频大全免费| 欧美日韩视频高清一区二区三区二| 别揉我奶头 嗯啊视频| 交换朋友夫妻互换小说| 大片免费播放器 马上看| 免费少妇av软件| 久久99精品国语久久久| 国产亚洲av嫩草精品影院| 欧美变态另类bdsm刘玥| 亚洲国产av新网站| 韩国高清视频一区二区三区| 在线看a的网站| 青春草亚洲视频在线观看| 中文字幕免费在线视频6| 天堂中文最新版在线下载 | 大又大粗又爽又黄少妇毛片口| 欧美亚洲 丝袜 人妻 在线| 亚洲aⅴ乱码一区二区在线播放| 精品久久久噜噜| 久久99精品国语久久久| 男人和女人高潮做爰伦理| 欧美老熟妇乱子伦牲交| 亚洲av免费在线观看| 青春草国产在线视频| 亚洲av二区三区四区| 一边亲一边摸免费视频| 亚洲国产精品国产精品| 久久国产乱子免费精品| 久久久久精品性色| 免费av不卡在线播放| 日本猛色少妇xxxxx猛交久久| 女的被弄到高潮叫床怎么办| 两个人的视频大全免费| 国模一区二区三区四区视频| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看| 国产美女午夜福利| 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 国产高清三级在线| 搡女人真爽免费视频火全软件| 成年人午夜在线观看视频| 日韩一区二区视频免费看| 精品酒店卫生间| 韩国高清视频一区二区三区| 亚洲va在线va天堂va国产| 免费观看的影片在线观看| 亚洲精品成人久久久久久| 国产精品一区二区三区四区免费观看| 亚洲在久久综合| 国产 一区精品| 色播亚洲综合网| 国产精品国产三级国产专区5o| 免费黄频网站在线观看国产| 白带黄色成豆腐渣| 久久久久精品久久久久真实原创| 丝袜美腿在线中文| 欧美高清性xxxxhd video| 亚洲国产日韩一区二区| 免费播放大片免费观看视频在线观看| 91久久精品国产一区二区成人| 真实男女啪啪啪动态图| 视频区图区小说| 亚洲国产色片| 国产精品久久久久久久电影| 亚洲色图综合在线观看| 国产伦理片在线播放av一区| 国产男女内射视频| 国产日韩欧美亚洲二区| 五月伊人婷婷丁香| 水蜜桃什么品种好| 精品熟女少妇av免费看| 亚洲av一区综合| 老师上课跳d突然被开到最大视频| 亚洲精品久久久久久婷婷小说| 能在线免费看毛片的网站| 中文资源天堂在线| 中文精品一卡2卡3卡4更新| 18禁动态无遮挡网站| 在线观看av片永久免费下载| 性色avwww在线观看| 欧美激情国产日韩精品一区| 亚洲精品国产色婷婷电影| 欧美最新免费一区二区三区| 亚洲天堂av无毛| 久久久久网色| 最近中文字幕2019免费版| 超碰av人人做人人爽久久| 丝袜喷水一区| 26uuu在线亚洲综合色| 免费黄色在线免费观看| 国内精品美女久久久久久| 国精品久久久久久国模美| 波多野结衣巨乳人妻| 菩萨蛮人人尽说江南好唐韦庄| 韩国高清视频一区二区三区| 午夜日本视频在线| 五月天丁香电影| 三级经典国产精品| av在线播放精品| 国产在线男女| 啦啦啦在线观看免费高清www| 欧美日韩亚洲高清精品| 久久99蜜桃精品久久| 亚洲国产精品成人久久小说| 亚洲欧美精品自产自拍| 国产成人免费无遮挡视频| 成人一区二区视频在线观看| 免费av毛片视频| 性色av一级| 国内精品美女久久久久久| 性插视频无遮挡在线免费观看| 国产精品国产av在线观看| 一级毛片 在线播放| 国产精品伦人一区二区| 街头女战士在线观看网站| 2021少妇久久久久久久久久久| 色综合色国产| 久久97久久精品| 亚洲人成网站在线观看播放| 国产一区亚洲一区在线观看| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 久久99热6这里只有精品| 国产高潮美女av| 国产精品嫩草影院av在线观看| 99久久九九国产精品国产免费| 亚洲内射少妇av| 又爽又黄a免费视频| 五月伊人婷婷丁香| 极品教师在线视频| 欧美日韩在线观看h| 狂野欧美激情性bbbbbb| 日本免费在线观看一区| 18禁在线无遮挡免费观看视频| 80岁老熟妇乱子伦牲交| 亚洲国产欧美人成| 99久久精品一区二区三区| 亚洲av免费在线观看| 久久精品国产亚洲av天美| 欧美日韩在线观看h| 久久精品夜色国产| 22中文网久久字幕| 亚洲成人久久爱视频| 一级毛片黄色毛片免费观看视频| 丰满少妇做爰视频| 国产男女超爽视频在线观看| 精品久久久久久久末码| 人人妻人人澡人人爽人人夜夜| 黄色怎么调成土黄色| 日韩在线高清观看一区二区三区| 久久久久久国产a免费观看| 亚洲在线观看片| 精品一区二区免费观看| 97在线人人人人妻| 日日撸夜夜添| 老师上课跳d突然被开到最大视频| 久久久久精品久久久久真实原创| 国产v大片淫在线免费观看| 自拍偷自拍亚洲精品老妇| 狂野欧美白嫩少妇大欣赏| 国产高潮美女av| 亚洲av免费高清在线观看| 精品国产一区二区三区久久久樱花 | 在线观看人妻少妇| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 日本色播在线视频| 国产精品一二三区在线看| 男女边摸边吃奶| 久久久色成人| 高清午夜精品一区二区三区| 国产成人免费无遮挡视频| 免费看光身美女| 超碰av人人做人人爽久久| 亚洲精华国产精华液的使用体验| 久久精品综合一区二区三区| 看黄色毛片网站| 亚洲内射少妇av| 成人综合一区亚洲| 一级a做视频免费观看| 亚洲三级黄色毛片| 黑人高潮一二区| 极品少妇高潮喷水抽搐| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 亚州av有码| 别揉我奶头 嗯啊视频| 又大又黄又爽视频免费| 亚洲图色成人| 国产精品久久久久久精品电影| 久久精品国产亚洲av天美| 精品国产三级普通话版| 亚洲美女搞黄在线观看| 一区二区av电影网| 乱码一卡2卡4卡精品| 久久久久久久久大av| 极品少妇高潮喷水抽搐| 亚洲av免费在线观看| 观看免费一级毛片| 亚洲精品乱久久久久久| 午夜精品一区二区三区免费看| 一级毛片 在线播放| 尤物成人国产欧美一区二区三区| 日韩大片免费观看网站| 中文精品一卡2卡3卡4更新| 久久久久久久午夜电影| 色视频www国产| 三级经典国产精品| 男女国产视频网站| 精品亚洲乱码少妇综合久久| 午夜日本视频在线| 插阴视频在线观看视频| 国产又色又爽无遮挡免| 精品久久久久久电影网| 亚洲人与动物交配视频| 女人久久www免费人成看片| 国产成人精品婷婷| 国产精品久久久久久久久免| 国内精品美女久久久久久| 最近中文字幕高清免费大全6| 亚洲在线观看片| 亚洲精品日韩av片在线观看| 在线 av 中文字幕| 国产白丝娇喘喷水9色精品| 我要看日韩黄色一级片| 久久久a久久爽久久v久久| 亚洲欧美精品专区久久| 欧美3d第一页| 亚洲天堂av无毛| 在线亚洲精品国产二区图片欧美 | 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 香蕉精品网在线| 久久久久精品久久久久真实原创| videos熟女内射| 色吧在线观看| 可以在线观看毛片的网站| 国产男女内射视频| 午夜免费男女啪啪视频观看| 日韩成人av中文字幕在线观看| 免费看光身美女| 国产成人免费无遮挡视频| 亚洲经典国产精华液单| 久久国内精品自在自线图片| 亚洲欧美精品自产自拍| 少妇 在线观看| 2022亚洲国产成人精品| 欧美日韩亚洲高清精品| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 九色成人免费人妻av| 少妇的逼水好多| 舔av片在线| 国产精品99久久久久久久久| 一二三四中文在线观看免费高清| 人人妻人人澡人人爽人人夜夜| 插阴视频在线观看视频| 成人毛片a级毛片在线播放| 五月天丁香电影| 91午夜精品亚洲一区二区三区| 在线免费观看不下载黄p国产| 国产精品福利在线免费观看| 亚洲av不卡在线观看| 精品久久久久久电影网| 亚洲三级黄色毛片| 综合色av麻豆| 免费大片黄手机在线观看| 久久午夜福利片| 高清午夜精品一区二区三区| 中文乱码字字幕精品一区二区三区| 人妻夜夜爽99麻豆av| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 高清毛片免费看| 各种免费的搞黄视频| 中文字幕久久专区| 久久人人爽人人爽人人片va| 久久这里有精品视频免费| av在线天堂中文字幕| 亚洲在久久综合| av卡一久久| 色网站视频免费| 在线 av 中文字幕| 亚洲欧美精品专区久久| 国产在线一区二区三区精| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 一级片'在线观看视频| 亚洲精品一二三| 久久99热这里只频精品6学生| 亚洲av二区三区四区| 日本欧美国产在线视频| av在线观看视频网站免费| 又爽又黄无遮挡网站| 免费播放大片免费观看视频在线观看| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 免费av观看视频| 亚洲,一卡二卡三卡| 午夜激情久久久久久久| 国产免费又黄又爽又色| 18禁在线播放成人免费| 国产精品一区二区三区四区免费观看| 毛片女人毛片| 亚洲av中文av极速乱| 欧美日韩亚洲高清精品| 国产成年人精品一区二区| 男人添女人高潮全过程视频| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 午夜老司机福利剧场| 夫妻性生交免费视频一级片| 国产人妻一区二区三区在| 国产黄频视频在线观看| 日日啪夜夜爽| 国产精品三级大全| 色5月婷婷丁香| 免费观看性生交大片5| 国产有黄有色有爽视频| 国产黄色视频一区二区在线观看| 亚洲熟女精品中文字幕| tube8黄色片| 有码 亚洲区| 国产探花在线观看一区二区| 日韩强制内射视频| 精品酒店卫生间| 国产精品熟女久久久久浪| 日韩av不卡免费在线播放| 久久99热6这里只有精品| 久久精品久久久久久噜噜老黄| 国产成人精品一,二区| 精品久久久久久久久av| 国产欧美亚洲国产| 美女cb高潮喷水在线观看| freevideosex欧美| 国产乱人偷精品视频| 久久久久久久午夜电影| 三级经典国产精品| 国产美女午夜福利| 日韩不卡一区二区三区视频在线| 欧美潮喷喷水| 亚洲自偷自拍三级| 99久久人妻综合| 你懂的网址亚洲精品在线观看| 男的添女的下面高潮视频| 日韩亚洲欧美综合| 2021天堂中文幕一二区在线观| 久久久久网色| 亚洲欧美精品自产自拍| 丝袜美腿在线中文| 十八禁网站网址无遮挡 | 久久久久久久久久久丰满| 亚洲精品,欧美精品| 免费看av在线观看网站| 精品久久久久久久久亚洲| 女人被狂操c到高潮| 欧美日韩精品成人综合77777| 国产精品蜜桃在线观看| 亚洲欧洲国产日韩| 国产精品伦人一区二区| 欧美日韩精品成人综合77777| 亚洲久久久久久中文字幕| 亚洲伊人久久精品综合| 免费黄频网站在线观看国产| 在线观看三级黄色| 日韩欧美一区视频在线观看 | 我的老师免费观看完整版| 2021少妇久久久久久久久久久| 国产精品一区二区性色av| 最后的刺客免费高清国语| 精品国产露脸久久av麻豆| 日韩免费高清中文字幕av| 成年女人看的毛片在线观看| 性插视频无遮挡在线免费观看| 国产成人精品福利久久| 欧美激情在线99| h日本视频在线播放| 亚洲四区av| 国产精品一区二区性色av| 免费大片18禁| 亚洲国产日韩一区二区| 少妇裸体淫交视频免费看高清| 高清视频免费观看一区二区| 超碰av人人做人人爽久久| 久久鲁丝午夜福利片| 一个人观看的视频www高清免费观看| www.av在线官网国产| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 一二三四中文在线观看免费高清| 亚洲经典国产精华液单| 亚洲欧美清纯卡通| 国产免费福利视频在线观看| 五月天丁香电影| 男女边摸边吃奶| 99热这里只有是精品在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品乱码久久久v下载方式| 中国三级夫妇交换| 国产成人a区在线观看| 一边亲一边摸免费视频| 成年女人在线观看亚洲视频 | 伊人久久精品亚洲午夜| 欧美xxxx黑人xx丫x性爽| av网站免费在线观看视频| 一级av片app| 一本色道久久久久久精品综合| 亚洲av中文字字幕乱码综合| 国产精品精品国产色婷婷| 伦精品一区二区三区| 啦啦啦中文免费视频观看日本| 人人妻人人澡人人爽人人夜夜| 日韩免费高清中文字幕av| 两个人的视频大全免费| 国产视频内射| 日韩精品有码人妻一区| 亚洲国产成人一精品久久久| 国产欧美另类精品又又久久亚洲欧美| kizo精华| 免费观看性生交大片5| 在线a可以看的网站| 亚洲国产精品成人久久小说| 国产 一区精品| 真实男女啪啪啪动态图| 国产男人的电影天堂91| 久久亚洲国产成人精品v| 免费观看性生交大片5| 欧美日韩一区二区视频在线观看视频在线 | 国产综合精华液| 五月开心婷婷网| 国产片特级美女逼逼视频| 亚洲天堂av无毛| 欧美一级a爱片免费观看看| 观看美女的网站| 国产精品.久久久| 麻豆成人午夜福利视频| 亚洲精品久久久久久婷婷小说| 99热这里只有精品一区| 啦啦啦在线观看免费高清www| 在线精品无人区一区二区三 | 亚洲精品中文字幕在线视频 | 特大巨黑吊av在线直播| 久久久成人免费电影| 男女无遮挡免费网站观看| 免费大片18禁| 亚洲电影在线观看av| 久久久久性生活片| 少妇人妻 视频| 成人综合一区亚洲| 国产精品爽爽va在线观看网站| 人妻制服诱惑在线中文字幕| 免费观看的影片在线观看| 中文欧美无线码| 日本一本二区三区精品|