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

    植物位點特異性甲基化研究的引物設計及PCR體系優(yōu)化

    2016-10-20 08:56:16王鶴潼崔偉娜賈春云惠秀娟臺培東成智博
    農業(yè)環(huán)境科學學報 2016年9期
    關鍵詞:亞硫酸鹽甲基化特異性

    王鶴潼,宋 婕,崔偉娜,曹 霞,何 蕾,賈春云,惠秀娟,臺培東,成智博,劉 宛*

    (1.中國科學院沈陽應用生態(tài)研究所污染生態(tài)與環(huán)境工程重點實驗室,沈陽 110016;2.遼寧何氏醫(yī)學院,沈陽110163;3.遼寧大學,沈陽 110036;4.上海應用技術學院,上海 201418;5.沈陽農業(yè)大學,沈陽 110866;6.通遼市農業(yè)科學研究院蔬菜所,內蒙古通遼 028000)

    植物位點特異性甲基化研究的引物設計及PCR體系優(yōu)化

    王鶴潼1,2,宋婕3,崔偉娜4,曹霞5,6,何蕾3,賈春云1,惠秀娟1,臺培東1,成智博1,劉宛1*

    (1.中國科學院沈陽應用生態(tài)研究所污染生態(tài)與環(huán)境工程重點實驗室,沈陽 110016;2.遼寧何氏醫(yī)學院,沈陽110163;3.遼寧大學,沈陽 110036;4.上海應用技術學院,上海 201418;5.沈陽農業(yè)大學,沈陽 110866;6.通遼市農業(yè)科學研究院蔬菜所,內蒙古通遼 028000)

    通過設計大量引物,探索在植物甲基化引物設計中長、短引物的設計方法及其相應的PCR條件,同時比較不同Taq酶對甲基化率影響。結果表明:17~30 bp短引物適合使用Touch-down程序PCR,但特異性差、成功率低;31~50 bp長引物使用55℃退火60℃延伸的PCR條件成功率最高,其中反向引物的長度及Tm小于正向引物較佳;高保真酶因其3′→5′外切酶活性不宜用于甲基化研究,而EpiTaqTMHS在PCR結果中表現(xiàn)最佳。

    植物甲基化;引物設計;Taq酶;PCR條件

    王鶴潼,宋婕,崔偉娜,等.植物位點特異性甲基化研究的引物設計及PCR體系優(yōu)化[J].農業(yè)環(huán)境科學學報,2016,35(9):1686-1693.

    WANG He-tong,SONG Jie,CUI Wei-na,et al.Primer design and PCR condition optimization for plant site-specific methylation[J].Journal of Agro-Environment Science,2016,35(9):1686-1693.

    表觀遺傳修飾是一種不同于核苷酸序列的可遺傳信息,主要包括DNA甲基化和組蛋白修飾,研究發(fā)現(xiàn)其在植物生長發(fā)育、基因表達調控以及抗逆脅迫中起著至關重要的作用[1-2]。DNA甲基化廣泛存在于多種有機體中,且多數(shù)發(fā)生在胞嘧啶第5碳原子上形成5-甲基胞嘧啶。在植物中,DNA甲基化是一種普遍現(xiàn)象,20%~30%的基因組胞嘧啶處于甲基化狀態(tài),并分布于CG、CHG、CHH位點(H代表A、G或T)[3-5]。隨著人們對植物DNA甲基化研究愈來愈關注,開發(fā)出眾多檢測DNA甲基化的方法。根據(jù)研究內容和目的的不同,可分為檢測全基因組甲基化和位點特異性甲基化兩類方法,其中前者主要以高效液相色譜法[6]、基于甲基化敏感內切酶(如:Methylation sensitive amplification polymorphism,MSAP[7]和Methylation-sensitive arbitrarily primed PCR,MSAP-PCR[8-10]等)和基于South ern探針雜交(如:CGI array和cDNA microarray等)[2]的方法為主,由于芯片技術需要開發(fā)大量探針且成本較高,而高效液相色譜法條件優(yōu)化較復雜,因此在植物研究中多數(shù)使用基于甲基化敏感內切酶的方法。

    在位點特異性甲基化研究中,以使用重亞硫酸鹽限制酶組合分析[11](COBRA)、甲基化特異性PCR[12](MSP)和重亞硫酸鹽測序法(Bisulfite Sequencing)[13-15]為主。隨著DNA測序技術發(fā)展,由于植物基因組CG含量較高而無法使用焦磷酸測序,重亞硫酸鹽測序法的檢測效率和準確性均高于COBRA法和MSP法,目前位點特異性甲基化研究中重亞硫酸鹽測序法的使用最為廣泛[16-17]。然而,隨著該方法在植物上的不斷應用,發(fā)現(xiàn)其中仍然存在一些較棘手的技術性問題,即DNA目的片段重亞硫酸鹽修飾后擴增結果不理想和測序后甲基化率不穩(wěn)定,其中涉及引物的設計、PCR條件的優(yōu)化和Taq酶的選擇三個方面[18]。由于植物基因組CG含量與動物及人相比較高[19],基因組修飾后重復序列過多,導致應用該方法時在設計引物和選擇目標區(qū)域上會受到極大限制,同時還因所需目的基因區(qū)域的范圍限制,導致所設計引物的Tm、特異性和穩(wěn)定性不佳,使其在修飾后的DNA上不易擴增。另外,我們還發(fā)現(xiàn)不同的Taq酶對測序后的DNA甲基化率有明顯影響,而這些問題都讓植物甲基化研究受到阻礙。因此,本研究從完善方法著手,首次通過分析不同引物設計方式和不同引物擴增體系對甲基化結果的影響,優(yōu)化植物研究中重亞硫酸鹽測序法的實驗條件,最終建立一個兼顧實驗操作性和結果穩(wěn)定性的重亞硫酸鹽測序法實驗體系,為植物的位點特異性甲基化研究提供依據(jù)。

    1 材料與方法

    1.1植物材料的培養(yǎng)與處理

    選用哥倫比亞型擬南芥種子,經(jīng)10%次氯酸鈉消毒后,于滅菌超純水中4℃浸種24 h打破春化。浸種后在無菌臺上將種子轉移至滅菌后0.5×MS培養(yǎng)基(0.5%蔗糖)中,并于21℃、10 h光照培養(yǎng)15 d[9]。

    1.2DNA的提取

    將不同培養(yǎng)罐中的幼苗混合,使用北京康為世紀生物科技有限公司的新型植物基因組DNA提取試劑盒(CW0531)提取混合幼苗地上部DNA,然后使用Eppendorf Biophotometor Plus檢測DNA含量,根據(jù)測定結果上樣50 ng DNA進行瓊脂糖電泳,用TaKaRa DL2000 DNA Marker進行含量校準并分裝成每管1 μg,-70℃冰箱保存。

    1.3重亞硫酸鹽測序法

    按照說明書使用 EZ DNA Methylation-GoldTMKits(Zymo research)修飾 100 ng DNA。根據(jù)TAIR (www.tair.org)擬南芥DNA序列數(shù)據(jù)庫,使用Primer 5和Oligo 7自設計引物PCR擴增錯配修復基因MLH1、MSH2、MSH6和MSH7啟動子片段(引物序列及PCR條件詳見“結果與分析”);擴增產(chǎn)物使用1%瓊脂糖(TaKaRa)凝膠電泳檢測;目的片段使用TaKaRaMiniBEST Agarose Gel DNA Extraction Kit回收純化;純化后的PCR片段連接載體后克隆測序(上海英駿生物技術有限公司)。

    2 結果與分析

    2.1植物甲基化研究中短引物的設計與應用

    由于植物核基因組CG含量較高,修飾后造成基因組存在大量重復序列使引物特異性降低,因此本實驗室為研究脅迫誘導擬南芥幼苗錯配修復基因啟動子甲基化變異,設計了大量長度在17~30 bp的正、反向引物,發(fā)現(xiàn)其中僅有3對引物MLH1-1F/MLH1-1R、MLH1-410F/MLH1-410R和MSH6_P6F/MSH6_P6R3(表1)擴增出所需目的片段,且擴增穩(wěn)定(圖1)。通過比較引物相關信息,發(fā)現(xiàn)這3對成功擴增的引物,其正向引物的長度和Tm均大于反向引物(表1),說明在正、反向引物的設計上,需要反向引物的長度短于正向引物,前者可以優(yōu)先退火結合模板;而正引物較長,可增加引物的特異性,其Tm高于反向引物利于在反向引物擴增后,迅速結合反向引物所在的模板鏈,最終完成整個PCR過程。

    圖1表明,在PCR擴增條件優(yōu)化上,通過使用多種方式擴增后發(fā)現(xiàn),Touch-down對于這三種甲基化短引物的擴增均有較好的效果:引物MLH1-1F/MLH1-1R使用55℃退火30 s,且每個循環(huán)降0.5℃,運行17個循環(huán),然后47℃退火30 s運行21個循環(huán)為退火條件;引物MLH1-410F/MLH1-410R使用55℃退火30s,且每個循環(huán)降1℃,運行11個循環(huán),然后45℃退火30 s運行27個循環(huán)為退火條件;引物MSH6_P6F/ MSH6_P6R3使用56℃退火30 s,且每個循環(huán)降0.5℃,運行9個循環(huán),然后52℃退火30 s運行30個循環(huán)為退火條件。

    表1 實驗中自設計的甲基化短引物Table 1 Self-designed short primers for methylation in this study

    圖1 自設計短引物電泳驗證Figure 1 Electriphoresis validation of products amplified by 3 pairs of short,self-designed primers

    該結果說明,甲基化短引物設計的實現(xiàn)需對應優(yōu)化PCR條件,根據(jù)正、反向引物的Tm決定Touchdown程序的遞降范圍,原則上Touch-down程序退火最高溫度設定于正、反向引物Tm之間,不可超過反向引物Tm過多,最低退火溫度與最高退火溫度之差原則上小于10℃。而且,引物3′端的ΔG需要低于普通PCR引物,更利于甲基化引物退火后快速引發(fā)延伸反應。

    結果還發(fā)現(xiàn),雖然引物MSH6_P6F/MSH6_P6R3能夠擴增出較穩(wěn)定的PCR產(chǎn)物,且與DNA marker比較后其片段大小也與目的片段基本相似,但經(jīng)過克隆測序并與NCBI數(shù)據(jù)庫比對后發(fā)現(xiàn),所擴增的序列與葉綠體DNA高度同源(圖2),并非所需目的基因MSH6的PCR片段。這說明甲基化短引物仍存在弊端,無法避免基因組修飾后導致的引物特異性嚴重下降。

    2.2植物甲基化研究中長引物的設計與應用

    由于甲基化短引物無法避免基因組修飾后導致的引物特異性嚴重下降,致使PCR產(chǎn)物可能為非目的片段。為了克服該問題,本研究設計了大量長度為31~50 bp的長引物。表2中提供了部分該實驗中的自設計引物,其中覆蓋了擬南芥MSH2和MSH7基因的啟動子區(qū)0~-600,且每個基因含有2對有效引物,表2中帶“*”號引物為有效引物,其中dMSH2-11F/ dMSH2-5R和dMSH2b-7F/dMSH2b-4R為MSH2基因啟動子引物,dMSH7-3F/dMSH7-2R和dMSH7b-4F/dMSH7b-2R為MSH7基因啟動子引物。根據(jù)實驗結果該4對引物均為有效擴增(圖3),且經(jīng)測序比對驗證是目的片段。

    由于引物長度的加大可以使引物特異性得到很大的提升,正向引物已很難非特異性的結合非對應模板。根據(jù)短引物設計所得經(jīng)驗,在設計引物中依然盡量保持反向引物的長度小于正向引物,至少保證絕對不大于正向引物長度。結合表2和圖3結果分析,表明反向引物僅需不大于正向引物,兩者Tm僅需相差不大,且有較低的3′端ΔG,便可擴增出目的片段。在引物的PCR程序選擇上,通過對各種方式的嘗試,最終發(fā)現(xiàn),使用55℃退火30 s、60℃延伸2 min適合甲基化長引物擴增,而且其特異性和穩(wěn)定性較高。

    圖2 引物MSH6_P6F/MSH6_P6R3擴增產(chǎn)物Blast驗證Figure 2 Blast validation of PCR products amplified by primers MSH6_P6F/MSH6_P6R3

    表2顯示,隨著引物加長,不可避免的問題便是引物中覆蓋了許多可修飾位點,但這些位點是否被修飾未知,因此只能通過使用簡并堿基代替。過多的簡并會導致引物的特異性和穩(wěn)定性下降,因此本實驗在設計引物時,每個引物均設計成3~6對簡并形式(表2)。通過比較有效引物和無效引物,在我們成功擴增的4對引物中,發(fā)現(xiàn)并非所有可修飾位點全部簡并是最佳的,其3′端需要簡并,尤其是3′端的前1/3部分,然而其前3個堿基不能出現(xiàn)可修飾位點;其他位置的可修飾位點并無明顯的規(guī)律,一般情況下非5′端連續(xù)3個可修飾位點的堿基需要簡并,同時連續(xù)可修飾位點應避開3′端。

    2.3不同Taq酶對甲基化率的影響

    本實驗室對擬南芥錯配修復基因的甲基化研究發(fā)現(xiàn),不同的Taq酶對PCR擴增效果和甲基化率有明顯的差別。因此本實驗選用3種不同類型的Taq酶(表3)對相同DNA的MLH1基因啟動子進行擴增,結果如圖4所示,發(fā)現(xiàn)PrimeSTARHS所得結果的甲基化率遠大于其他兩種Taq酶,而專為重亞硫酸鹽修飾模板擴增開發(fā)的EpiTaqTMHS比常規(guī)的TaKaRaTaq的甲基化率略高。該結果表明,高保真酶PrimeSTARHS由于其結果過高并不適合于甲基化研究,而TaKaRaTaq由于其酶能力所限,無法檢測出目的片段內全部甲基化變異位點。

    3 討論

    3.1引物設計及優(yōu)化在植物位點特異性甲基化研究中的重要性

    圖3 自設計長引物電泳驗證Figure 3 Electriphoresis validation of products amplified by 4 pairs of short,self-designed primers

    隨著分子生物學和生物技術的高速發(fā)展,對于一般生物科學研究而言,實驗中使用的關鍵技術已不再是制約,而且各類儀器、技術服務以及試劑盒的開發(fā)極大地方便了科學研究并簡化了實驗過程。因此,實驗和研究人員所面臨的難題便是課題或實驗的設計,以及方法或技術的選擇和優(yōu)化[20]。對于植物位點特異性甲基化研究而言,其特異性位點的引物設計、選擇和優(yōu)化便是難點。

    表2 實驗中自設計的甲基化長引物Table 2 Self-designed long primers for methylation in this study

    表3 本研究所選用的Taq酶Table 3 Taq polymerase used in this study

    一方面,由于植物基因組CG含量較高[19],導致經(jīng)過重亞硫酸鹽修飾后基因組重復序列過多,引物結合模板的特異性下降,而且胞嘧啶過多也給設計引物造成較大的困難,因為在指定的區(qū)域或范圍內,很難找到較優(yōu)異的引物。另一方面,由于PCR擴增的模板是修飾過的,與原DNA差異較大,網(wǎng)上并沒有相關的數(shù)據(jù)庫,也無法對所設計的引物進行特異性分析(Primer-blast),只有通過PCR實踐檢測驗證。這些原因均極大地阻礙了實驗的進度和研究的開展,導致即使擁有高效的DNA提取、修飾及純化的試劑盒,實驗卻仍然停滯不前。因此,在植物位點特異性甲基化研究中,引物的設計和優(yōu)化尤為重要。

    3.2植物位點特異性甲基化的引物設計

    圖4 3種Taq酶擴增修飾后擬南芥MLH1基因啟動子的甲基化結果Figure 4 Methylation profile of Arabidopsis MLH1 promoter amplified by 3 Taq polymerase after sodium bisulfite modification

    本實驗室主要研究Cd脅迫下擬南芥的遺傳損傷[9-10,21-22],其中包括對錯配修復基因表觀遺傳損傷的研究,并為此設計了大量的引物。然而經(jīng)PCR擴增檢驗的結果并不理想,約20余對短引物(17~30 bp)中,僅有3對成功擴增(表1)。通過整理引物并總結經(jīng)驗,我們發(fā)現(xiàn)主要原因在于引物的特異性上,由于引物長度較短,修飾的DNA重復序列高,導致引物的非特異性極高,這便使目的片段很難擴增。然而根據(jù)這些所設計引物,我們對有效引物和失敗引物進行了比較分析,得到了一些有用的經(jīng)驗。我們發(fā)現(xiàn)正向引物的長度和Tm需要大于反向引物,且3′端ΔG也要保持較低(-6.0 kcal·mol-1左右),這樣引物的成功率會顯著提升(表1)。該現(xiàn)象是因為甲基化引物PCR與常規(guī)PCR有一重要區(qū)別,便是在實際擴增過程中只利用了一條鏈作為模板,即單鏈擴增,而與該模板鏈對應的引物便是反向引物,故在設計引物時,首先需考慮的便是讓反向引物更易退火且特異性較高。反向引物長度短于正向引物,便會先于正向引物退火結合到模板上;同時因3′端ΔG較低,會讓引物退火后快速有效的引發(fā)延伸反應,而正向引物盡管Tm高于反向引物,但因引物長度大導致其引物特異性略高,且沒有對應模板,不易退火。當反向引物擴增后,正向引物的高Tm會迅速讓其退火反向引物擴增鏈,并保持高于反向引物的退火能力及延伸能力,這樣可盡量保持兩條引物的同步擴增,又不至于單條鏈過于富集造成PCR效率下降。

    在成功擴增產(chǎn)物的引物中,經(jīng)測序比對驗證仍有非目的片段的擴增(圖2),因此短引物用作甲基化引物依然存在弊端。為此本實驗室又設計了大量長引物(31~50 bp)。根據(jù)表2發(fā)現(xiàn),除不同簡并形式外,長引物的成功率可在50%以上,其中設計引物原則基本與短引物保持一致,即反向引物更易于結合模板。另外,引物除5′端外要保持較低ΔG,即在引物內部穩(wěn)定性上(Internal stability)呈“W”型,并且不能存在ΔG過低的發(fā)夾或二聚體。另一方面,由于引物過長必然導致覆蓋較多的可修飾位點,引物當中的可變堿基需要簡并,然而不同的簡并形式會對PCR擴增造成直接影響。根據(jù)經(jīng)驗我們發(fā)現(xiàn)3′端若有可修飾位點則必須簡并,但該狀況越少越好,因為它直接影響引物延伸反應的引發(fā),但是3′末端的前3個堿基不允許存在可修飾位點。另外,需要避免2個堿基以上的連續(xù)可修飾位點,且不應出現(xiàn)在3′端,如含有則須簡并。終上所述,雖然長引物依然存在一些不可控因素,但與短引物相比,其引物設計成功率已得到極大提升,并且通過選擇適合的PCR程序和Taq酶也可在一定程度上提高其成功率。

    3.3PCR程序和Taq酶的選擇

    圖1結果表明,短引物擴增所使用的PCR程序是Touch-down模式。在該方式下,設定的起始退火溫度高于反向引物的Tm,使反向引物能有較高的特異性,而該溫度低于正向引物Tm,可在反向引物延伸后迅速以其為模板進行擴增。隨著前幾輪擴增,使有義模板富集,此時退火溫度已下降,可增加PCR擴增的整體效率以便得到較高含量的目的片段。但是,因為退火溫度只是影響PCR效率的次要因素,如果引物特異性差,PCR程序如何優(yōu)化也于事無補,所以長引物的選用成必然趨勢。

    表2和圖3結果表明,通過多種方式的嘗試和驗證,最終發(fā)現(xiàn)長引物的PCR程序如下:95℃變性30 s,55℃退火15~30 s,60℃延伸2~3 min。在該PCR程序中,選用長退火或延伸時間會造成PCR總時間過長而引起酶活力下降,因而需要根據(jù)不同Taq酶的擴增能力(主要指單位時間擴增片段長度)和活力的持久度,決定退火及延伸時間,此PCR條件與Henderson等[18]所用方法相似。該PCR條件的優(yōu)點是,對于Tm較高的長引物,55℃會讓引物迅速地在高特異性和穩(wěn)定性的位點上退火,然后在60℃進行邊退火邊延伸,這樣可使一些低穩(wěn)定性和可修飾的位點繼續(xù)退火,也可進行較穩(wěn)定的延伸。然而缺點便是PCR程序時間過長且對酶的性能要求較高,因為60℃并不是Taq最佳溫度,在該溫度下其延伸效率不足72℃的一半,致使延伸時間過長,而普通的Taq在2 h后性能會明顯下降,最終影響PCR的效果。

    圖4結果表明,在Taq酶的選擇上,PrimeSTAR HS并不適合甲基化PCR,因為所有的高保真酶都含有3′→5′核酸外切酶活性,而尿嘧啶與胸腺嘧啶相比少一個甲基,雖然在堿基配對時也可形成兩個氫鍵,但形成雙鏈螺旋結構時并不穩(wěn)定,可能會影響高保真酶的擴增,甚至會誤認其是錯配而將其外切掉,導致已修飾位點不易擴增而含甲基化位點的DNA可正常延伸,最終引起結果中高甲基化的極端現(xiàn)象。而且在部分品牌高保真酶的說明書上,雖然沒有強調卻有提示“不宜使用含dUTP或尿嘧啶的模板和引物”。因此,甲基化引物的PCR更注重對Taq酶的擴增能力。相對于普通的TaKaRaTaq,由于EpiTaqTMHS專為修飾DNA設計,其擴增能力得到加強,且在保真性能上不依賴外切酶活性,因此其結果優(yōu)于TaKaRaTaq。

    4 結論

    (1)17~30 bp短引物使用Touch-down程序PCR效果較好,但無法避免其特異性低的弊端;31~50 bp長引物可解決該問題,并以55℃退火60℃延伸PCR程序擴增,成功率最高,其中反向引物的長度及Tm低于正向引物較佳。

    (2)高保真酶因其外切酶活性所致的高甲基化現(xiàn)象,并不適宜用于甲基化研究,而EpiTaqTMHS由于專為修飾模板PCR設計,其有較強的擴增含尿嘧啶模板及不依賴3′→5′外切酶活性的保真能力,在實際應用中效果最佳。

    [1]Bird A.P.CpG-rich islands and the function of DNA methylation[J].Nature,1986,321(6067):209-213.

    [2]Ushijima T.Detection and interpretation of altered methylation patterns in cancer cells[J].Nat Rev Cancer,2005,5(3):223-231.

    [3]Richards E J.DNA methylation and plant development[J].Trends Genet,1997,13(8):319-323.

    [4]Zhang X,Yazaki J,Sundaresan A,et al.Genome-wide high-resolution mapping and functional analysis of DNA methylation in arabidopsis[J]. Cell,2006,126(6):1189-1201.

    [5]Wagner I,Capesius I.Determination of 5-methylcytosine from plant DNA byhigh-performanceliquidchromatography[J].BiochimicaEtBiophysica Acta,1981,654(1):52-56.

    [6]Johnston J W,Harding K,Bremner D H,et al.HPLC analysis of plant DNA methylation:A study of critical methodological factors[J].Plant Physiol Biochem,2005,43(9):844-853.

    [7]Reyna-Lopez G E,Simpson J,Ruiz-Herrera J.Differences in DNA methylation patterns are detectable during the dimorphic transition of fungi by amplification of restriction polymorphisms[J].Mol Gen Genet,1997,253(6):703-710.

    [8]Gonzalgo M L,Liang G,Spruck C H,et al.Identification and characterization of differentially methylated regions of genomic DNA by methylation-sensitive arbitrarily primed PCR[J].Cancer Res,1997,57(4):594-599.

    [9]王鶴潼,何蕾,宋杰,等.改進MSAP-PCR技術應用于Cd脅迫下擬南芥DNA甲基化分析[J].農業(yè)環(huán)境科學學報,2015,34(8):1618-1624.

    WANG He-tong,HE Lei,SONG Jie,et al.Assay of DNA methylation in Arabidopsis under Cd stress using improved MSAP-PCR technique[J]. Journal of Agro-Environment Science,2015,34(8):1618-1624.

    [10]Wang H,He L,Song J,et al.Cadmium-induced genomic instability in Arabidopsis:Molecular toxicological biomarkers for early diagnosis of cadmium stress[J].Chemosphere,2016,150:258-265.

    [11]Xiong Z,Laird P W.COBRA:A sensitive and quantitative DNA methylation assay[J].Nucleic Acids Res,1997,25(12):2532-2534.

    [12]Herman J G,Graff J R,Myohanen S,et al.Methylation-specific PCR:A novel PCR assay for methylation status of CpG islands[J].Proc Natl Acad Sci U S A,1996,93(18):9821-9826.

    [13]Sadri R,Hornsby P J.Rapid analysis of DNA methylation using new restriction enzyme sites created by bisulfite modification[J].Nucleic Acids Res,1996,24(24):5058-5059.

    [14]Cokus S J,F(xiàn)eng S,Zhang X,et al.Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning[J].Nature,2008,452(7184):215-219.

    [15]Lister R,O′Malley R C,Tonti-Filippini J,et al.Highly integrated single-base resolution maps of the epigenome in Arabidopsis[J].Cell,2008,133(3):523-536.

    [16]李利紅,儀慧蘭,王藝雯,等.二氧化硫脅迫誘導擬南芥NIT2基因DNA甲基化修飾[J].農業(yè)環(huán)境科學學報,2012,31(4):685-690.

    LI Li-hong,YI Hui-lan,WANG Yi-wen,et al.Sulfur dioxide induces DNA methylation alteration of a gene encoding nitrilase 2 protein in Arabidopsis plants[J].Journal of Agro-Environment Science,2012,31 (4):685-690.

    [17]王藝雯,李利紅,儀慧蘭.H2O2誘導擬南芥NIT2基因甲基化特征與轉錄水平改變[J].環(huán)境科學學報,2015,35(1):324-328.

    WANG Yi-wen,LI Li-hong,YI Hui-lan.DNA methylation alterations and gene transcription changes of NIT2 in Arabidopsis plants induced by hydrogen peroxide[J].Acta Scientiae Circumstantiae,2015,35(1):324-328.

    [18]Henderson I R,Chan S R,Cao X,et al.Accurate sodium bisulfite sequencing in plants[J].Epigenetics,2010,5(1):47-49.

    [19]侯雷平,李梅蘭.DNA甲基化與植物的生長發(fā)育[J].植物生理學報,2001,37(6):584-588.

    HOU Lei-ping,LI Mei-lan.DNA methylation and plant growth and development[J].Plant Physiology Journal,2001,37(6):584-588.

    [20]Fulnecek J,Kovarik A.How to interpret methylation sensitive amplified polymorphism(MSAP)profiles?[J].Bmc Genet,2014,15:2.

    [21]Liu W,Sun L,Zhong M,et al.Cadmium-induced DNA damage and mutations in Arabidopsis plantlet shoots identified by DNA fingerprinting[J].Chemosphere,2012,89(9):1048-1055.

    [22]Liu W,Yang Y S,Li P J,et al.Risk assessment of cadmium-contaminated soil on plant DNA damage using RAPD and physiological indices[J].J Hazard Mater,2009,161(2/3):878-883.

    Primer design and PCR condition optimization for plant site-specific methylation

    WANG He-tong1,2,SONG Jie3,CUI Wei-na4,CAO Xia5,6,HE Lei3,JIA Chun-yun1,HUI Xiu-juan1,TAI Pei-dong1,CHENG Zhi-bo1,LIU Wan1*
    (1.Key Laboratory of Pollution Ecology and Environmental Engineering,Institute of Applied Ecology,Chinese Academy of Sciences,Shenyang 110016,China;2.Liaoning He University,Shenyang 110163,China;3.Liaoning University,Shenyang 110036,China;4.Shanghai Institute of Technology,Shanghai 201418,China;5.Shenyang Agricultural University,Shenyang 110866,China;6.Institute of Vegetable Science,Tong Liao Academy of Agricultural Sciences,Tongliao 028000,China)

    The present study was to explore design methods for short and long primers in plant methylation research and their corresponding PCR conditions employing large number of designed primers and to analyze the effects on methylation rates using three different Taq polymerases.Results showed:a Touch-down PCR was suitable when using short primers of 17~30 bp,however,low specificity and success rates were observed;PCR performed the best when using long primers of 31~50 bp,provided that length and Tmof forward primers>those of reverse primers and 55℃and 60℃served as annealing and prolonging temperature,respectively;Super-Fidelity Taq polymerase like Prime STAR Risn′t suitable for methylation research because of its 3′→5′exonuclease activities.However,PCR performed the best when using EpiTaqTMHS.

    plant methylation;primer design;Taq polymerase;PCR condition

    Q943.2

    A

    1672-2043(2016)09-1686-08doi:10.11654/jaes.2016-0173

    2016-02-02

    國家自然科學基金項目(21347007,41272255,41472237);國家科技重大專項(2012ZX7505-001)

    王鶴潼(1985—),男,遼寧沈陽人,博士,從事植物分子毒理及表觀遺傳研究。E-mail:tony.w.china@hotmail.com

    劉宛E-mail:liuwan63@hotmail.com

    猜你喜歡
    亞硫酸鹽甲基化特異性
    亞硫酸鹽在食品中的研究進展
    核農學報(2021年7期)2021-12-07 12:00:17
    精確制導 特異性溶栓
    食品中亞硫酸鹽的使用及檢測
    BOPIM-dma作為BSA Site Ⅰ特異性探針的研究及其應用
    亞硫酸鹽在食品加工中的作用及其應用
    重復周圍磁刺激治療慢性非特異性下腰痛的臨床效果
    兒童非特異性ST-T改變
    鼻咽癌組織中SYK基因啟動子區(qū)的甲基化分析
    胃癌DNA甲基化研究進展
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    亚洲精品美女久久久久99蜜臀| 天天添夜夜摸| 人人妻人人澡人人看| 国内毛片毛片毛片毛片毛片| 电影成人av| 国产黄频视频在线观看| 欧美亚洲 丝袜 人妻 在线| 高清毛片免费观看视频网站 | 可以免费在线观看a视频的电影网站| 亚洲精品乱久久久久久| 老司机福利观看| 黑人猛操日本美女一级片| 久久久精品94久久精品| 青青草视频在线视频观看| 久久人人爽av亚洲精品天堂| 久久99热这里只频精品6学生| 美女福利国产在线| 久久国产亚洲av麻豆专区| 丁香六月欧美| 一区二区三区乱码不卡18| 咕卡用的链子| 热99国产精品久久久久久7| 午夜激情av网站| 亚洲欧美精品综合一区二区三区| 亚洲中文字幕日韩| 亚洲三区欧美一区| 久久久久久久国产电影| 色在线成人网| h视频一区二区三区| 国产av一区二区精品久久| 一二三四社区在线视频社区8| 国产三级黄色录像| 久久久水蜜桃国产精品网| 美女高潮到喷水免费观看| 亚洲九九香蕉| 男女边摸边吃奶| 18在线观看网站| 国产极品粉嫩免费观看在线| 精品一区二区三区四区五区乱码| 成人三级做爰电影| 成人18禁高潮啪啪吃奶动态图| 咕卡用的链子| 黄色视频在线播放观看不卡| 亚洲欧美激情在线| 国产精品久久久久久精品古装| 国产深夜福利视频在线观看| 亚洲av成人不卡在线观看播放网| 国产又爽黄色视频| 他把我摸到了高潮在线观看 | 久久ye,这里只有精品| 一区二区三区国产精品乱码| 9热在线视频观看99| 欧美日韩黄片免| 国产精品久久久av美女十八| kizo精华| 欧美黑人欧美精品刺激| 亚洲av电影在线进入| 国产精品 国内视频| 欧美日韩一级在线毛片| 精品熟女少妇八av免费久了| 操美女的视频在线观看| 精品欧美一区二区三区在线| 国产高清国产精品国产三级| 日日摸夜夜添夜夜添小说| 无限看片的www在线观看| 老司机在亚洲福利影院| 麻豆乱淫一区二区| 国产高清视频在线播放一区| 国产男靠女视频免费网站| 欧美午夜高清在线| 国产欧美亚洲国产| 99国产综合亚洲精品| 久久精品国产亚洲av高清一级| 日韩大码丰满熟妇| 视频区图区小说| 老司机福利观看| 精品一区二区三区av网在线观看 | 丝袜在线中文字幕| 丝袜喷水一区| 国产精品国产av在线观看| 五月开心婷婷网| 大型av网站在线播放| 91国产中文字幕| 国产成人免费观看mmmm| 国产日韩一区二区三区精品不卡| 国产av国产精品国产| 精品国产国语对白av| 自拍欧美九色日韩亚洲蝌蚪91| 性少妇av在线| 伦理电影免费视频| 99精品在免费线老司机午夜| 亚洲成人手机| 成年版毛片免费区| 动漫黄色视频在线观看| 午夜精品国产一区二区电影| 成年动漫av网址| 亚洲少妇的诱惑av| 精品一区二区三卡| 精品人妻在线不人妻| 美女高潮喷水抽搐中文字幕| 一本一本久久a久久精品综合妖精| 国产欧美亚洲国产| 国内毛片毛片毛片毛片毛片| 高清黄色对白视频在线免费看| 两性夫妻黄色片| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 淫妇啪啪啪对白视频| 高清欧美精品videossex| 在线观看免费视频日本深夜| 嫩草影视91久久| 欧美午夜高清在线| 露出奶头的视频| 久久精品91无色码中文字幕| 成人黄色视频免费在线看| 中文字幕色久视频| 欧美久久黑人一区二区| 午夜免费鲁丝| 欧美激情久久久久久爽电影 | 搡老岳熟女国产| 久久中文看片网| 黑丝袜美女国产一区| 汤姆久久久久久久影院中文字幕| 99精品久久久久人妻精品| bbb黄色大片| 免费久久久久久久精品成人欧美视频| 制服人妻中文乱码| 91精品国产国语对白视频| 欧美日韩av久久| 欧美 日韩 精品 国产| 啦啦啦免费观看视频1| 国产伦理片在线播放av一区| 在线天堂中文资源库| 性色av乱码一区二区三区2| 免费人妻精品一区二区三区视频| 两人在一起打扑克的视频| 狂野欧美激情性xxxx| 日本av手机在线免费观看| 香蕉久久夜色| 另类精品久久| 国产精品麻豆人妻色哟哟久久| 一本一本久久a久久精品综合妖精| 国产男女超爽视频在线观看| 麻豆av在线久日| 国产一区二区在线观看av| 国产成人啪精品午夜网站| 欧美乱码精品一区二区三区| 国产精品1区2区在线观看. | 国产老妇伦熟女老妇高清| 少妇粗大呻吟视频| 国产日韩欧美在线精品| 五月开心婷婷网| 男女边摸边吃奶| 夜夜夜夜夜久久久久| 最近最新中文字幕大全电影3 | 99riav亚洲国产免费| 两人在一起打扑克的视频| 国产一区二区 视频在线| 国产日韩欧美视频二区| 国产日韩欧美在线精品| 日韩熟女老妇一区二区性免费视频| 一级片免费观看大全| e午夜精品久久久久久久| 国产99久久九九免费精品| av不卡在线播放| 大型黄色视频在线免费观看| 亚洲精品美女久久久久99蜜臀| 高清黄色对白视频在线免费看| 日韩有码中文字幕| 久久国产精品人妻蜜桃| 亚洲自偷自拍图片 自拍| 悠悠久久av| 国产一卡二卡三卡精品| 国产福利在线免费观看视频| 香蕉国产在线看| 国产精品欧美亚洲77777| 一区二区三区乱码不卡18| 亚洲av成人不卡在线观看播放网| 大香蕉久久网| 成年人黄色毛片网站| 精品一区二区三区av网在线观看 | 久久天堂一区二区三区四区| 精品免费久久久久久久清纯 | 最新在线观看一区二区三区| 美女福利国产在线| 亚洲国产欧美日韩在线播放| 久久天堂一区二区三区四区| 99国产精品一区二区蜜桃av | 亚洲综合色网址| 9191精品国产免费久久| 男人舔女人的私密视频| 高清视频免费观看一区二区| 国产福利在线免费观看视频| 两人在一起打扑克的视频| 欧美另类亚洲清纯唯美| 日本av手机在线免费观看| 亚洲熟女毛片儿| 精品视频人人做人人爽| 精品国产一区二区三区久久久樱花| 99精品久久久久人妻精品| 久久ye,这里只有精品| 欧美精品亚洲一区二区| 国产真人三级小视频在线观看| 美女高潮喷水抽搐中文字幕| 纵有疾风起免费观看全集完整版| 国产成人免费观看mmmm| 国产免费视频播放在线视频| 国产成人欧美| 欧美黄色片欧美黄色片| 又大又爽又粗| av网站在线播放免费| 又黄又粗又硬又大视频| 99国产极品粉嫩在线观看| 亚洲男人天堂网一区| 香蕉国产在线看| 9色porny在线观看| 又紧又爽又黄一区二区| 日韩精品免费视频一区二区三区| 久久人妻福利社区极品人妻图片| 一区在线观看完整版| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品成人av观看孕妇| 两个人免费观看高清视频| 亚洲国产精品一区二区三区在线| 夜夜夜夜夜久久久久| 美女福利国产在线| 日韩中文字幕欧美一区二区| 精品视频人人做人人爽| 久久人妻熟女aⅴ| 国产单亲对白刺激| 高清欧美精品videossex| 涩涩av久久男人的天堂| 一个人免费在线观看的高清视频| 后天国语完整版免费观看| 久久热在线av| 麻豆成人av在线观看| h视频一区二区三区| 99国产精品一区二区蜜桃av | 午夜久久久在线观看| 欧美人与性动交α欧美软件| 精品人妻1区二区| 日韩视频一区二区在线观看| 精品福利永久在线观看| bbb黄色大片| 久久午夜综合久久蜜桃| 中文字幕高清在线视频| 黑人巨大精品欧美一区二区蜜桃| 精品国产乱码久久久久久男人| 成人国产av品久久久| 精品人妻在线不人妻| 国产精品 国内视频| 亚洲欧洲精品一区二区精品久久久| 中文亚洲av片在线观看爽 | 国产欧美日韩精品亚洲av| 成人手机av| 亚洲情色 制服丝袜| 欧美日韩黄片免| 亚洲精品国产精品久久久不卡| 99久久国产精品久久久| 精品少妇内射三级| 国产亚洲精品第一综合不卡| 国产精品一区二区在线观看99| 热99re8久久精品国产| av电影中文网址| 嫩草影视91久久| 99国产综合亚洲精品| 日韩大片免费观看网站| 老司机深夜福利视频在线观看| 国产又爽黄色视频| 久久精品91无色码中文字幕| 国产av国产精品国产| 99国产极品粉嫩在线观看| 亚洲国产精品一区二区三区在线| 午夜成年电影在线免费观看| 日本黄色视频三级网站网址 | 天天添夜夜摸| 在线av久久热| 美女国产高潮福利片在线看| 成年动漫av网址| 亚洲av国产av综合av卡| 9色porny在线观看| 中文字幕制服av| 叶爱在线成人免费视频播放| 午夜激情久久久久久久| 久久国产精品影院| 国产精品免费一区二区三区在线 | 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 一本大道久久a久久精品| 精品一品国产午夜福利视频| 久久性视频一级片| 在线观看免费高清a一片| 日韩成人在线观看一区二区三区| 亚洲国产精品一区二区三区在线| 欧美中文综合在线视频| 黄网站色视频无遮挡免费观看| 99精品欧美一区二区三区四区| 精品熟女少妇八av免费久了| 色播在线永久视频| 久久国产精品男人的天堂亚洲| 大型av网站在线播放| 国产人伦9x9x在线观看| 999精品在线视频| 成在线人永久免费视频| 一本—道久久a久久精品蜜桃钙片| 曰老女人黄片| 他把我摸到了高潮在线观看 | 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品久久久久5区| 成人18禁在线播放| 国产高清视频在线播放一区| 久久中文字幕人妻熟女| 国产亚洲欧美精品永久| av免费在线观看网站| 国产av一区二区精品久久| 91国产中文字幕| 亚洲av成人一区二区三| 巨乳人妻的诱惑在线观看| 又大又爽又粗| 啦啦啦中文免费视频观看日本| 波多野结衣av一区二区av| 中国美女看黄片| 亚洲三区欧美一区| 亚洲av美国av| 国产麻豆69| 亚洲欧美精品综合一区二区三区| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| 黑人操中国人逼视频| 久久狼人影院| 久久午夜亚洲精品久久| 午夜福利视频精品| tube8黄色片| 18禁美女被吸乳视频| 国产亚洲欧美精品永久| 久热爱精品视频在线9| 国产在线视频一区二区| 首页视频小说图片口味搜索| 国产精品免费视频内射| 国产三级黄色录像| 一级片'在线观看视频| 国产无遮挡羞羞视频在线观看| 丝袜美腿诱惑在线| 色尼玛亚洲综合影院| 999久久久精品免费观看国产| 国产欧美亚洲国产| 777米奇影视久久| 中文亚洲av片在线观看爽 | 欧美激情高清一区二区三区| 2018国产大陆天天弄谢| 老汉色∧v一级毛片| 国产亚洲精品久久久久5区| 国产高清国产精品国产三级| 黑人巨大精品欧美一区二区蜜桃| 母亲3免费完整高清在线观看| 精品国内亚洲2022精品成人 | 国产亚洲av高清不卡| 一个人免费在线观看的高清视频| tocl精华| 夫妻午夜视频| 亚洲精品乱久久久久久| 日本撒尿小便嘘嘘汇集6| 中文字幕人妻熟女乱码| 久久久久久久久久久久大奶| 亚洲av欧美aⅴ国产| 老司机影院毛片| 亚洲精品国产精品久久久不卡| 久久这里只有精品19| 国产精品熟女久久久久浪| 国产伦理片在线播放av一区| 亚洲成a人片在线一区二区| 麻豆成人av在线观看| 亚洲性夜色夜夜综合| 久久热在线av| 天天躁狠狠躁夜夜躁狠狠躁| 超色免费av| 老司机深夜福利视频在线观看| 99精品欧美一区二区三区四区| 性少妇av在线| 一边摸一边做爽爽视频免费| 少妇裸体淫交视频免费看高清 | 9色porny在线观看| 国产一区有黄有色的免费视频| 成人手机av| 亚洲av成人一区二区三| 国产一卡二卡三卡精品| 桃花免费在线播放| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 91老司机精品| 一个人免费看片子| 久久九九热精品免费| 老司机亚洲免费影院| 亚洲成av片中文字幕在线观看| 中文字幕色久视频| 蜜桃国产av成人99| 欧美精品一区二区大全| 91九色精品人成在线观看| 男女床上黄色一级片免费看| 成人国语在线视频| 老汉色∧v一级毛片| 男女免费视频国产| 亚洲 欧美一区二区三区| 一区福利在线观看| 国产麻豆69| 麻豆乱淫一区二区| 国产伦人伦偷精品视频| 国产精品电影一区二区三区 | 美女午夜性视频免费| 18禁黄网站禁片午夜丰满| 日韩中文字幕欧美一区二区| 亚洲精品中文字幕一二三四区 | 啦啦啦 在线观看视频| 亚洲综合色网址| 18禁观看日本| 国产黄频视频在线观看| 欧美国产精品va在线观看不卡| 国产一卡二卡三卡精品| 少妇的丰满在线观看| 久久久久久人人人人人| 亚洲精华国产精华精| 亚洲av成人一区二区三| 人妻一区二区av| 大片电影免费在线观看免费| 国产精品影院久久| 香蕉国产在线看| 又紧又爽又黄一区二区| 国产伦人伦偷精品视频| 丁香六月欧美| 黄片小视频在线播放| 成人三级做爰电影| 精品人妻1区二区| 男人操女人黄网站| 欧美日韩av久久| 精品国产一区二区三区久久久樱花| 久久青草综合色| 久久中文看片网| √禁漫天堂资源中文www| 亚洲成人免费av在线播放| 日韩有码中文字幕| 一级a爱视频在线免费观看| 少妇裸体淫交视频免费看高清 | 99九九在线精品视频| 十八禁网站网址无遮挡| 色尼玛亚洲综合影院| 少妇精品久久久久久久| 蜜桃国产av成人99| 五月开心婷婷网| 成人影院久久| 啦啦啦免费观看视频1| 99久久精品国产亚洲精品| 亚洲精品中文字幕在线视频| 丝袜人妻中文字幕| 色老头精品视频在线观看| 国产一区二区三区在线臀色熟女 | 人人妻人人澡人人看| 热re99久久精品国产66热6| 69av精品久久久久久 | 国产成人精品在线电影| 欧美大码av| 99re在线观看精品视频| 国产一区二区在线观看av| 成人av一区二区三区在线看| 久久精品91无色码中文字幕| 国产成+人综合+亚洲专区| 欧美日韩中文字幕国产精品一区二区三区 | 1024香蕉在线观看| 色综合欧美亚洲国产小说| 色在线成人网| 99热国产这里只有精品6| 水蜜桃什么品种好| 狠狠狠狠99中文字幕| 久久精品国产综合久久久| 亚洲精品粉嫩美女一区| 色尼玛亚洲综合影院| 日本欧美视频一区| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 亚洲av欧美aⅴ国产| 一本久久精品| 国产免费现黄频在线看| 91老司机精品| 亚洲视频免费观看视频| 欧美一级毛片孕妇| 97在线人人人人妻| 别揉我奶头~嗯~啊~动态视频| 欧美成人午夜精品| 男女下面插进去视频免费观看| 欧美乱码精品一区二区三区| 国产精品.久久久| 精品国内亚洲2022精品成人 | 国产精品亚洲一级av第二区| 国产成人精品无人区| 老司机深夜福利视频在线观看| 国产精品一区二区在线观看99| 高潮久久久久久久久久久不卡| 亚洲成人国产一区在线观看| 国产精品欧美亚洲77777| av在线播放免费不卡| 久久久久久亚洲精品国产蜜桃av| 国产深夜福利视频在线观看| 精品熟女少妇八av免费久了| 国产成人欧美| 日本av免费视频播放| 亚洲成人免费电影在线观看| 男女边摸边吃奶| 亚洲午夜精品一区,二区,三区| 亚洲av日韩在线播放| 国产精品久久久久成人av| 久久国产精品大桥未久av| 丝袜美足系列| 岛国毛片在线播放| 脱女人内裤的视频| 91精品三级在线观看| 日本av手机在线免费观看| 中文字幕人妻丝袜一区二区| 丰满迷人的少妇在线观看| 国产不卡av网站在线观看| 丰满人妻熟妇乱又伦精品不卡| 另类亚洲欧美激情| 免费av中文字幕在线| 国产av一区二区精品久久| 国产亚洲欧美精品永久| 国产xxxxx性猛交| 久久久久国产一级毛片高清牌| 国产男靠女视频免费网站| 亚洲国产av新网站| av线在线观看网站| 精品一区二区三区四区五区乱码| 亚洲av成人不卡在线观看播放网| 欧美日韩一级在线毛片| 12—13女人毛片做爰片一| 亚洲色图 男人天堂 中文字幕| 建设人人有责人人尽责人人享有的| 黄片大片在线免费观看| 久热爱精品视频在线9| 久久久久国内视频| 一区二区三区激情视频| 嫁个100分男人电影在线观看| 中文亚洲av片在线观看爽 | 欧美在线一区亚洲| 老熟妇乱子伦视频在线观看| av网站免费在线观看视频| 97人妻天天添夜夜摸| av网站免费在线观看视频| 欧美日韩精品网址| 麻豆av在线久日| 建设人人有责人人尽责人人享有的| 汤姆久久久久久久影院中文字幕| 国产视频一区二区在线看| 久久免费观看电影| 97人妻天天添夜夜摸| 电影成人av| 国产亚洲av高清不卡| 亚洲av欧美aⅴ国产| 岛国在线观看网站| 亚洲av欧美aⅴ国产| 黑人巨大精品欧美一区二区mp4| 18在线观看网站| 国产精品久久久久久精品古装| 99香蕉大伊视频| 黑人巨大精品欧美一区二区蜜桃| 夜夜夜夜夜久久久久| 国产成人系列免费观看| 亚洲欧洲日产国产| 黄片小视频在线播放| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| netflix在线观看网站| 女人久久www免费人成看片| 亚洲中文日韩欧美视频| 久久人妻福利社区极品人妻图片| 在线观看舔阴道视频| 国产亚洲精品久久久久5区| 国产精品久久久久久精品电影小说| 69av精品久久久久久 | 女人久久www免费人成看片| 日日爽夜夜爽网站| 99国产综合亚洲精品| 女性被躁到高潮视频| 日韩人妻精品一区2区三区| 日韩一区二区三区影片| 中文字幕人妻熟女乱码| 乱人伦中国视频| 成人国产av品久久久| 欧美日本中文国产一区发布| 亚洲成人手机| 精品免费久久久久久久清纯 | 精品少妇一区二区三区视频日本电影| 国产黄频视频在线观看| 国产精品久久久久久精品古装| 国产精品美女特级片免费视频播放器 | 中文欧美无线码| 亚洲专区中文字幕在线| 99久久精品国产亚洲精品| 日韩视频在线欧美| 搡老乐熟女国产| 欧美日本中文国产一区发布| 午夜日韩欧美国产| 18禁美女被吸乳视频| 又紧又爽又黄一区二区| 免费观看a级毛片全部| 国产在线一区二区三区精| av又黄又爽大尺度在线免费看| av视频免费观看在线观看| 高潮久久久久久久久久久不卡| 国产男女超爽视频在线观看| 欧美国产精品va在线观看不卡| 麻豆av在线久日| 国产高清视频在线播放一区| 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区视频在线观看| 丰满迷人的少妇在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产一区二区在线观看av|