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

    PCR產(chǎn)物的高通量測序方法及優(yōu)化

    2016-03-07 01:57:25李桂瀾胥傳來
    關(guān)鍵詞:建庫上機(jī)數(shù)據(jù)量

    李桂瀾, 胥傳來

    (1.北京大學(xué) 生命科學(xué)學(xué)院,北京100871;2.食品科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,江南大學(xué),江蘇 無錫214122)

    PCR產(chǎn)物的高通量測序方法及優(yōu)化

    李桂瀾1, 胥傳來2

    (1.北京大學(xué) 生命科學(xué)學(xué)院,北京100871;2.食品科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,江南大學(xué),江蘇 無錫214122)

    PCR產(chǎn)物的高通量測序被廣泛應(yīng)用于功能基因篩選、腫瘤相關(guān)基因的突變和甲基化檢測等。在高通量測序技術(shù)中,建庫和上機(jī)測序?qū)嶒?yàn)一直是決定最終DNA數(shù)據(jù)質(zhì)量的關(guān)鍵。作者優(yōu)化了PCR產(chǎn)物樣品的DNA文庫制備條件和體系,設(shè)計(jì)了適合PCR產(chǎn)物特點(diǎn)上機(jī)測序方法,將已建庫的PCR樣品中混入一定量的基因組標(biāo)準(zhǔn)品后再上機(jī)測序,由此保證了測序數(shù)據(jù)的高質(zhì)量和低冗余度。為低多樣性DNA樣品高通量測序技術(shù)方法運(yùn)用提供方法上的參考。

    PCR產(chǎn)物;高通量測序技術(shù);樣品文庫制備;上機(jī)測序

    DNA測序是常用的分子生物學(xué)研究技術(shù),通過測序分析能提供最真實(shí)可靠的基因序列信息。從1977年第一代傳統(tǒng) DNA雙脫氧鏈末端終止的sanger測序法問世以來,DNA測序技術(shù)經(jīng)歷了快速發(fā)展[1-2]?,F(xiàn)在基于邊合成邊測序(sequencing by synthesis)為主的第二代高通量測序技術(shù)逐漸被廣泛運(yùn)用于基礎(chǔ)研究、臨床醫(yī)學(xué)診斷和腸道微生物群與營養(yǎng)及代謝研究領(lǐng)域[2-5]。二代測序以單個讀長較短、通量大為特點(diǎn),實(shí)現(xiàn)了對單個物種的轉(zhuǎn)錄組和基因組細(xì)致全貌分析,也被稱為深度測序或高通量測序[6]。二代測序研究平臺從最初Roche公司的454焦磷酸測序和Illumina公司早先的Solexa測序技術(shù)發(fā)展到目前廣泛使用的Hiseq、Miseq系列,和Life Technoloiges的Ion PGM和Ion Proton的測序,測序技術(shù)無論從速度和通量上得到了前所未有的發(fā)展[7]。目前以Illumina為技術(shù)平臺的Hiseq X 5系統(tǒng)、Hiseq3000/4000最新的測序儀器也相繼發(fā)布,使二代測序樣品建庫更加高效,讀長(reads)更長,通量更高。

    PCR產(chǎn)物測序通過設(shè)計(jì)恰當(dāng)引物探針,運(yùn)用各種PCR技術(shù)將待測序的目的基因片段擴(kuò)增產(chǎn)物進(jìn)行測序[8-10]。PCR產(chǎn)物的高通量測序技術(shù)被輔助運(yùn)用于細(xì)胞文庫中功能基因的篩選,腫瘤相關(guān)基因突變及甲基化檢測和微生物菌群的營養(yǎng)及代謝相關(guān)探索研究[11-12]。以Illumina為平臺的高通量測序,對片段大小在150~350 bp的PCR混合物樣品,與基因組 DNA和微量的 ChIP-seq(Chromatin Immuno Precipitation,ChIP)樣品的文庫制備方法不同[13-14]。雖然目前建庫技術(shù)整體一直在被優(yōu)化,測序通量和測序長度也在不斷的進(jìn)行技術(shù)革新使之達(dá)到更高的測序要求,但I(xiàn)llumina的邊合成邊測序技術(shù)核心仍無法改變,也決定了Illumina二代測序的技術(shù)局限性。所以根據(jù)PCR產(chǎn)物樣品自身的特點(diǎn),文庫制備和上機(jī)測序?qū)嶒?yàn)都需要進(jìn)行特別的設(shè)計(jì)和處理。作者主要以Illumina系列的二代測序平臺,在普通的ChIP建庫方法,基于片段已知的DNA樣品通過兩步末端修飾反應(yīng)后加DNA通用接頭的實(shí)驗(yàn)原理,設(shè)計(jì)建立了一種適合PCR混合物樣品性質(zhì)特點(diǎn)的文庫制備及上機(jī)操作實(shí)驗(yàn),以確保高效低成本的建庫技術(shù)及上機(jī)測序中得到最佳的測序數(shù)據(jù)質(zhì)量。作者以Illumina為代表的高通量測序技術(shù)特點(diǎn)及具體實(shí)驗(yàn)方法,為不同的研究領(lǐng)域中選擇不同測序方法提供參考。

    1 材料與方法

    1.1 試劑與材料

    1.1.1 試劑 ChIP-Seq Sample Prep Master Mix試劑盒:美國 New England Biolabs公司;NEBNext Multiplex Oligos;Agencourt AMPure XP的 DNA純化磁珠:美國BECKMAN COULTER公司;Brilliant SYBR Green QPCR試劑盒:美國Agilent公司;上機(jī)雙端測序試劑 Miseq Reagent Kit V3 (2×150 cycles): 美國 Illumina公司;100 bp plus DNA Marker相對分子質(zhì)量標(biāo)準(zhǔn)品:北京全式金生物技術(shù)公司;100%乙醇、Tris、Tween20和 NaOH:美國Sigma公司。1.5 mL低DNA吸附(LoBind)微量離心管:美國Eppendorf。

    1.1.2 DNA分子 New England Biolabs公司提供:接頭分子(Adaptor):5’-pGATCGGAAGAGCACACG TCTGAACTCCAGTC/ideoxyU/ACACTCTTTCCCTAC ACGACGCTCTTCCGATCT-3’;上游引物:5’-AATG ATACGGCGACCACCGAGATCTACACTCTTTCCCTA CACGACGCTCTTCCGATCT-3’;下游引物:5’-CAA GCAGAAGACGGCATACGAGATAAGCTAGTGACTG GAGTTCAGACGTGTGCTCTTCCGATCT-3’;下游引物中的下劃線加粗的6個堿基為第10號標(biāo)簽(index)DNA分子。Takara合成熒光定量PCR引物:上游引物P1:5’-AATGATACGGCGACCACCGA-3’;下游引物P2:5’-CAAGCAGAAGACGGCATACG A-3’。

    1.2 儀器

    Miseq第二代高通量測序儀:美國Illumina公司;美國Bio-Rad核酸電泳儀、凝膠成像儀和PCR核酸擴(kuò)增儀;電子天平和pH計(jì):瑞士METTLER TOLEDO公司;NanoDrop2000超微量分光光度計(jì):美國Thermo Scientific公司;Qubit 2.0熒光定量儀:美國Life Technologies公司;Mx300P熒光定量PCR儀:美國安捷倫公司;Fragment Analyzer全自動毛細(xì)管電泳儀及其配套試劑:美國Advanced Analytical。ThermoMixer恒溫混勻儀:美國Eppendorf。

    1.3 實(shí)驗(yàn)方法

    1.3.1 測序樣品文庫制備 切膠純化回收后的PCR產(chǎn)物樣品T1、T2和T3,瓊脂糖凝膠電泳圖中條帶是分布在200~300 bp寬峰,見圖1。用NanoDrop2000測定質(zhì)量濃度分別為24.2、25.2、26.2 ng/μL;A260/A280的OD值為1.68、1.59和1.76。

    建庫的初始樣品用 Qubit 2.0熒光定量儀對PCR產(chǎn)物樣品準(zhǔn)確定量后各取2 μL(DNA總量在10~50 ng)于1.5 mL的Lobind離心管,準(zhǔn)備第一步的平末端修復(fù)反應(yīng):1 μL末端修復(fù)反應(yīng)聚合酶和5μL的10×末端修復(fù)反應(yīng)緩沖液,用RNase free H2O補(bǔ)足至反應(yīng)終體積50 μL,恒溫混勻儀中30℃反應(yīng)20 min。 結(jié)束后加入80 uL的AMPure XP磁珠純化反應(yīng)后樣品,用44 uL RNase free H2O洗脫。第二步直接向樣品溶液加入1 μL Klenow片段酶和5 μL的10×反應(yīng)緩沖液(含0.2 mmol/L dATP),在37℃反應(yīng)30 min,完成末端加“A”修飾反應(yīng)。同樣用80 μL磁珠純化,19 μL的RNase free H2O洗脫。

    圖1 二代測序的PCR產(chǎn)物樣品T1-T3Fig.1 PCR products T1-T3 for next-generation sequencing

    將DNA接頭分子稀釋到5 μmol/L后取1 μL到樣品溶液中準(zhǔn)備連接反應(yīng),加6 μL的5×快速連接緩沖液和4 μL的快速DNA連接酶,小心混勻,室溫25℃靜置孵育15 min。結(jié)束后加入試劑盒的USER酶2 μL,在37℃反應(yīng)15 min完成接頭開環(huán)剪切。最后用與樣品等倍體積30 μL的磁珠純化兩遍,檢測濃度。取10 μL(總量在5~10 ng)樣品進(jìn)行PCR反應(yīng),循環(huán)數(shù)減少至15次。最后用與樣品等體積的磁珠純化兩遍,30 μL的ddH2O洗脫,檢測濃度。

    1.3.2 建庫后樣品質(zhì)量檢控 將建庫后樣品稀釋至0.5~2 ng/μL,分別用Fragment Analyzer全自動毛細(xì)管電泳儀進(jìn)行片段分析;用Mx300Ps熒光定量PCR儀對樣品進(jìn)行熒光定量檢測。選用的標(biāo)準(zhǔn)品為已準(zhǔn)確定量并上機(jī)測序后的建庫樣品。用10 mmol/L Tris-HCl pH 8.0,0.05%Tween 20緩沖液按10倍梯度配制成0.002~20 pmol/L五個不同濃度以獲得定量標(biāo)準(zhǔn)曲線。建庫樣品再次稀釋 104倍到 0.1~ 2 pmol/L后取2 μL,2×Brilliant SYBR Green Master Mix 5 μL,10 μmol/L引物P1/P2 Mix 0.4 μL,最后ddH2O補(bǔ)足為10 μL的反應(yīng)體系,陰性對照組NTC(No Template Control)和所有樣品及標(biāo)準(zhǔn)品做三組重復(fù)。

    [Final Conc(nmol/L)]=[QPCR]×340×104×10-6/Fragm

    式中,[QPCR]代表儀器根據(jù)標(biāo)準(zhǔn)曲線檢測出來的平均樣品濃度 (fM);340 bp代表選用的DNA標(biāo)準(zhǔn)品的片段大??;104代表樣品稀釋倍數(shù);Fragm代表建庫后樣品的平均片段大小。

    1.3.3 Illumina Miseq上機(jī)測序 對已建庫樣品T2準(zhǔn)備上機(jī)測序,將T2樣品與Phix基因組標(biāo)準(zhǔn)樣品等摩爾混合,用10 mmol/L Tris-HCl pH 8.0緩沖液稀釋至終濃度2 nmol/L,取10 μL的上機(jī)混合樣品,加入10 μL 2 mol/L的NaOH充分混勻后室溫變性5 min。用雙端測序試劑盒里的雜交緩沖液將樣品稀釋到18 pmol/L,取600 μL上機(jī)測序。Illumina數(shù)據(jù)收集和分析軟件:Miseq Control Software;Real Time Analysis(RTA);Offline Basecaller(OLB);CASAVA軟件用于進(jìn)一步數(shù)據(jù)分析。

    2 結(jié)果與分析

    2.1 建庫后樣品檢測

    PCR產(chǎn)物樣品T1,T2和T3為同一樣本的3個平行生物組重復(fù)。建庫結(jié)束后檢測終質(zhì)量濃度為:22.2、15.3、16.1 ng/μL各30 μL體系。隨機(jī)選取樣品T2的建庫及上機(jī)測序結(jié)果分析。用 Fragment Analyzer分析初始樣品和建庫后樣品片段分布,見圖2-3。樣品建庫后片段的出峰位置出現(xiàn)移動,增加約120 bp,滿足接頭序列成功連接后的片段分布情況。根據(jù)標(biāo)準(zhǔn)品所得的熒光定量標(biāo)準(zhǔn)曲線其相關(guān)系數(shù)為Rsq:99.9%,計(jì)算得到稀釋后的3個建庫樣品濃度12.46、6.78、7.81 nmol/L,滿足上機(jī)測序要求。T2樣品準(zhǔn)備上機(jī)測序。

    圖2 高通量測序樣品T2的DNA片段分析Fig.2 DNA fragment analysis of sample T2

    圖3 測序樣品T2建庫后的DNA片段分析Fig.3 DNA fragment analysis of final library sample T2

    2.2 上機(jī)測序結(jié)果

    Illumina Miseq測序中最主要指標(biāo)參數(shù)是數(shù)據(jù)質(zhì)量Q30值和測序數(shù)據(jù)量統(tǒng)計(jì)結(jié)果見圖4-5。其中百分比≥Q30代表每輪測序堿基的質(zhì)量不低于99.9%所占的百分比。統(tǒng)計(jì)T2樣品測序最終的Q值分布(橫坐標(biāo))和得到的數(shù)據(jù)量(縱坐標(biāo),單位是million)。統(tǒng)計(jì)得到全部測序總數(shù)據(jù)量為5.8G,數(shù)據(jù)質(zhì)量≥Q30的占94.3%。最終統(tǒng)計(jì)得到的T2樣品原始測序數(shù)據(jù)量(Raw Data)為3.725 G。進(jìn)一步生物信息過濾處理得有效可用數(shù)據(jù)(clean data)為3.705 G,數(shù)據(jù)有效率為99.46%,Q20值為94.09%,Q30值為92.85%;GC含量比例為48.19%。結(jié)果滿足后續(xù)對T2測序樣品的生物學(xué)研究分析的基本要求。

    圖4 樣品T2在Miseq上機(jī)測序得到的全部數(shù)據(jù)Q值(Quality Score)分布統(tǒng)計(jì)Fig.4 Statistic of the QScore distribution on illumina Miseq

    圖5 Hiseq2000的測序芯片 (flowcell)中第4條測序通道(lane 4)里所有樣品都只是PCR產(chǎn)物樣品測序時的數(shù)據(jù)Q30分布統(tǒng)計(jì)Fig.5 Statistic of the Qscore distribution without the phix standard mixture on lane 4 at Hiseq 2000

    2.3 上機(jī)測序?qū)嶒?yàn)優(yōu)化分析

    當(dāng)前二代高通量測序正突飛猛進(jìn)地不斷進(jìn)行技術(shù)革新,但以Illumina為平臺的邊合成邊測序的核心技術(shù)原理保持不變[2]。在上機(jī)測序過程中前25個堿基測序循環(huán)(cycle)結(jié)果質(zhì)量的統(tǒng)計(jì)方法無法改變,依舊決定了整個測序的數(shù)據(jù)質(zhì)量。測序中監(jiān)測參數(shù)Cluster PF(Pass Filter)仍是重要的上機(jī)質(zhì)檢標(biāo)準(zhǔn)。當(dāng)測序的前25個堿基低質(zhì)量的數(shù)據(jù)有兩個以上(即PF<60%),測序儀器將判定這條讀長質(zhì)檢不通過。而PCR產(chǎn)物樣品DNA多樣性差,在測序中極易引起兩個堿基的識別可靠性低于60%,從而減損最終測序的Q30值和判定合格測序數(shù)據(jù)總量,見圖6。圖4為Illumina Hiseq2000測序中的第5測序通道(Lane 5)標(biāo)準(zhǔn)基因組樣品測序監(jiān)測的%Base正常分布情況,Cluster PF(%)為95.23±1.15,可獲得的數(shù)據(jù)量約34 G;圖7為同時測序中第4通道(Lane 4)的樣品當(dāng)全為PCR產(chǎn)物時%Base分布情況,Cluster PF(%)為21.24±14.83。使最終得到的合格數(shù)據(jù)量和Q30的統(tǒng)計(jì)值降低,測序通道Lane 4里全部測序數(shù)據(jù)質(zhì)量≥Q30僅為52.9%,最終獲得的測序數(shù)據(jù)量降低到4.5 G。

    作者嘗試將待測的PCR產(chǎn)物樣品T2中混合了等摩爾的Phix基因組標(biāo)準(zhǔn)品后進(jìn)行Miseq測序,監(jiān)測得到的Cluster PF(%)為89.81±0.72。最終目的樣品T2占整個有效測序數(shù)據(jù)量的61.77%,大大改善了PCR混合物測序樣品的數(shù)據(jù)質(zhì)量和有效數(shù)據(jù)總量。對于多樣性差的其它待測序DNA樣品,比如RRBS 樣 品 (Reduced Representation Bisulfite Sequencing,簡化的表觀亞硫酸氫鹽測序),由于含有固定的酶切位點(diǎn)序列降低了樣品的多樣性,但通過混入一定量比例的多樣性好的基因組標(biāo)準(zhǔn)品,實(shí)驗(yàn)結(jié)果顯示得到的測序數(shù)據(jù)質(zhì)量可提高約50%。

    圖6 在Illumina Hiseq2000測序儀器中監(jiān)測的第5條測序通道標(biāo)準(zhǔn)基因組樣品的堿基分布圖Fig.6 Base distribution of genomic DNA sequencing on lane 5 of Illumina Hiseq 2000

    目前還可以嘗試降低待測序樣品上樣量的方法改善低多樣性樣品引起的低數(shù)據(jù)質(zhì)量的情況。但研究表明至少要降低一半的PCR樣品量才可能達(dá)到略微的效果,也大大折損了測序數(shù)據(jù)總量,并非經(jīng)濟(jì)有效的方案。而選擇將多樣性好的基因組樣品摻入到多樣性差的樣品中上機(jī)測序的方法,實(shí)驗(yàn)操作簡單,能改善整個測序通道的樣品多樣性表現(xiàn),從而提高測序質(zhì)量,同時也能保證得到的最終有效的測序數(shù)據(jù)總量損失最少。

    圖7 第4條測序通道為全部PCR產(chǎn)物樣品測序時的堿基分布圖Fig.7 Base distribution of all PCR samples sequencing on lane 4 of Illumina hiseq 2000

    3 結(jié)語

    作者研究探索了一套對普通PCR產(chǎn)物樣品的高通量測序建庫和上機(jī)實(shí)驗(yàn)方法,對低多樣性樣品的高通量測序方法做了探索研究,通過對樣品的上機(jī)測序數(shù)據(jù)監(jiān)測,證實(shí)所建立的方法保證樣品制備的高效經(jīng)濟(jì)和測序質(zhì)量的穩(wěn)定可靠。

    [1]MARDIS E R.A decade’s perspective on DNA sequencing technology[J].Nature,2011,470(7333):198-203.

    [2]METZKER M L.Sequencing technologies-the next generation[J].Nat Rev Genet,2010,11(1):31-46.

    [3]NIEDRINGHAUS T P,Milanova D,Kerby M B,et al.Landscape of next-generation sequencing technologies[J].Anal Chem,2011,83(12):4327-4341.

    [4]王興春,楊致榮,王敏,等.高通量測序技術(shù)及其應(yīng)用[J].中國生物工程雜志,2012,32(1):109-114. WANG Xingchun,YANG Zhirong,WANG Min,et al.High-throughput sequencing technology and its application[J].China Biotechnology,2012,32(1):109-114.(in Chinese)

    [5]MCCANN J C,WICKERSHAM T A,Loor J J,et al.High-throughput methods redefine the rumen microbiome and its relationship with nutrition and metabolism[J].Bioinformatics and Biology Insights,2014,8:109-125.

    [6]MOROZOVA O,MARRA M A.Applications of next-generation sequencing technologies in functional genomics[J].Genomics,2008,92(5):255-264.

    [7]MARDIS E R.Next-generation sequencing platforms[J].Annu Rev Anal Chem(Palo Alto Calif),2013(6):287-303.

    [8]YU C,ZHANG Y,YAO S,et al.A PCR based protocol for detecting indel mutations induced by TALENs and CRISPR/Cas9 in Zebrafish[J].PLoS One,2014,9(6):e98282.

    [9]NELSON M C,MORRISON H G,BENJAMINO J,et al.Analysis,optimization and verification of Illumina-generated 16S rRNAgene amplicon surveys[J].PLoS One,2014,9(4):e94249.

    [10]BOUTIN S,SEVELLEC M,Pavey S A,et al.A fast,highly sensitive double-nested PCR-based method to screen fish immunobiomes[J].Mol Ecol Resour,2012,12(6):1027-1039.

    [11]ZHOU Y,ZHU S,CAI C,et al.High-throughput screening of a CRISPR/Cas9 library for functional genomics in human cells[J]. Nature,2014,509(7501):487-491.

    [12]CAPORASO J G,LAUBER C L,WALTERS W A,et al.Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms[J].ISME J,2012,6(8):1621-1624.

    [13]AMIR A,ZEISEL A,ZUK O,et al.High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions[J].Nucleic Acids Res,2013,41(22):e205.

    [14]BOWMAN S K,SIMON M D,DEATON A M,et al.Multiplexed Illumina sequencing libraries from picogram quantities of DNA [J].BMC Genomics,2013,14:466.

    Methods and Optimization of High Throughput Sequencing Technologies on PCR Products

    LI Guilan1, XU Chuanlai2
    (1.School of Life Sciences,Peking University,Beijing 100871,China;2.State Key Laboratory of Food Science &Technology,Jiangnan University,Wuxi 214122,China)

    The development of the new generation sequencing technology extends the application field of gene sequencing.High throughput sequencing of PCR products has been widely used in functional gene screening,mutation and methylation detection of tumor related genes,etc.In the high throughput sequencing technology,database and computer sequencing experiments are the key factor to decide the quality of DNA data.Based on Illumina sequencing instruments,the sample preparation method and sequencing running for PCR products were optimized,the results make us clearly understand the primary features of high-throughput sequencing techniques,and provide us important references of sequencing methods for low diversity DNA samples to address biological questions of interest.

    PCR products,high-throughput sequencing,sample library preparation,sequencing runs

    Q 523.8

    A

    1673—1689(2016)12—1317—06

    2015-02-03

    國家“十二五”科技支撐計(jì)劃項(xiàng)目(2012BAC01B071)。

    李桂瀾(1980—),女,重慶人,理學(xué)博士,工程師,主要從事DNA測序技術(shù)的應(yīng)用及生物分子相互作用方面的研究。E-mail:liglan@pku.edu.cn

    猜你喜歡
    建庫上機(jī)數(shù)據(jù)量
    周金應(yīng)
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    計(jì)算Lyapunov指數(shù)的模糊C均值聚類小數(shù)據(jù)量法
    標(biāo)準(zhǔn)化護(hù)理程序?qū)w外膜肺氧合術(shù)患者上機(jī)各階段用時的影響
    高刷新率不容易顯示器需求與接口標(biāo)準(zhǔn)帶寬
    寬帶信號采集與大數(shù)據(jù)量傳輸系統(tǒng)設(shè)計(jì)與研究
    電子制作(2019年13期)2020-01-14 03:15:18
    面向建庫與制圖一體化的規(guī)則研究
    中文期刊回溯建庫的實(shí)踐與思考——以貴州省圖書館為例
    建設(shè)用地節(jié)約集約利用評價數(shù)據(jù)庫建庫流程:以西安市為例
    河南科技(2014年5期)2014-02-27 14:08:42
    基于數(shù)據(jù)字典的空間數(shù)據(jù)庫通用建庫技術(shù)
    黄色毛片三级朝国网站| 男人添女人高潮全过程视频| 亚洲精品av麻豆狂野| 各种免费的搞黄视频| 91aial.com中文字幕在线观看| 最近中文字幕2019免费版| 国产在线免费精品| a级片在线免费高清观看视频| 飞空精品影院首页| 狂野欧美激情性bbbbbb| 成人无遮挡网站| 美女福利国产在线| 97在线视频观看| 日韩三级伦理在线观看| 哪个播放器可以免费观看大片| 天天操日日干夜夜撸| 99国产综合亚洲精品| 亚洲国产日韩一区二区| 国产精品一二三区在线看| 午夜老司机福利剧场| freevideosex欧美| 亚洲高清免费不卡视频| 久久久久精品久久久久真实原创| 五月伊人婷婷丁香| 国产精品一国产av| 久久99蜜桃精品久久| 美女国产视频在线观看| 一边摸一边做爽爽视频免费| 欧美xxxx性猛交bbbb| 三级国产精品片| 一边摸一边做爽爽视频免费| 韩国高清视频一区二区三区| 性色av一级| 国产精品国产三级国产专区5o| 最近的中文字幕免费完整| 精品酒店卫生间| 九九爱精品视频在线观看| 久久久久久人妻| av在线老鸭窝| 免费人成在线观看视频色| 国产午夜精品一二区理论片| 热99久久久久精品小说推荐| 国产伦理片在线播放av一区| 久久av网站| 少妇高潮的动态图| 国产成人精品在线电影| 中文字幕另类日韩欧美亚洲嫩草| 天天操日日干夜夜撸| 交换朋友夫妻互换小说| 日韩av不卡免费在线播放| 天天躁夜夜躁狠狠久久av| 如日韩欧美国产精品一区二区三区| 男女啪啪激烈高潮av片| 男女啪啪激烈高潮av片| 欧美精品一区二区免费开放| 欧美日韩精品成人综合77777| 久久这里只有精品19| 女人精品久久久久毛片| 国产成人精品婷婷| 久久影院123| 观看美女的网站| 国产av国产精品国产| 性色avwww在线观看| 男人爽女人下面视频在线观看| 国产精品99久久99久久久不卡 | 国产日韩欧美亚洲二区| 欧美丝袜亚洲另类| 久久久久久久国产电影| 高清在线视频一区二区三区| 久久精品国产亚洲av天美| 久久久国产欧美日韩av| 欧美97在线视频| 欧美激情 高清一区二区三区| 尾随美女入室| 少妇高潮的动态图| 免费观看在线日韩| 亚洲一区二区三区欧美精品| 91国产中文字幕| 亚洲欧洲国产日韩| 黑人高潮一二区| 18禁裸乳无遮挡动漫免费视频| 国产激情久久老熟女| 丰满迷人的少妇在线观看| 精品国产国语对白av| av不卡在线播放| 99热这里只有是精品在线观看| 另类精品久久| 新久久久久国产一级毛片| 巨乳人妻的诱惑在线观看| 热99久久久久精品小说推荐| 五月开心婷婷网| 最近中文字幕2019免费版| 国产免费福利视频在线观看| 日韩av在线免费看完整版不卡| 国产熟女午夜一区二区三区| 亚洲精品,欧美精品| 国产无遮挡羞羞视频在线观看| 亚洲精品456在线播放app| 亚洲精品一区蜜桃| 老女人水多毛片| 国产精品人妻久久久久久| 国产欧美亚洲国产| 高清黄色对白视频在线免费看| 亚洲精品视频女| 一区在线观看完整版| 国产免费福利视频在线观看| 精品福利永久在线观看| 99香蕉大伊视频| 九色成人免费人妻av| 久久99热这里只频精品6学生| 午夜老司机福利剧场| 亚洲欧洲国产日韩| 日韩中文字幕视频在线看片| 免费观看性生交大片5| 国产精品一区www在线观看| 熟妇人妻不卡中文字幕| 秋霞伦理黄片| 女人久久www免费人成看片| 欧美精品一区二区大全| 国产一区二区在线观看日韩| 高清av免费在线| 亚洲综合精品二区| 久久国内精品自在自线图片| 国产精品一二三区在线看| 97在线人人人人妻| 熟女av电影| 香蕉丝袜av| 日本与韩国留学比较| 一区二区av电影网| 大码成人一级视频| 久久 成人 亚洲| 亚洲国产精品成人久久小说| 天天躁夜夜躁狠狠躁躁| 美女中出高潮动态图| 国产白丝娇喘喷水9色精品| 亚洲精品,欧美精品| 爱豆传媒免费全集在线观看| av天堂久久9| 一级毛片我不卡| 久久久久久久国产电影| 免费看不卡的av| 精品国产露脸久久av麻豆| av播播在线观看一区| 日韩中文字幕视频在线看片| 天美传媒精品一区二区| 又黄又粗又硬又大视频| 免费看av在线观看网站| 免费人成在线观看视频色| 性高湖久久久久久久久免费观看| 成年美女黄网站色视频大全免费| 99香蕉大伊视频| 国产亚洲欧美精品永久| 嫩草影院入口| 欧美日韩av久久| 又黄又爽又刺激的免费视频.| 两个人免费观看高清视频| 精品少妇黑人巨大在线播放| 永久网站在线| 亚洲在久久综合| 亚洲av电影在线观看一区二区三区| 国产色爽女视频免费观看| 精品国产露脸久久av麻豆| 狂野欧美激情性xxxx在线观看| 亚洲国产欧美在线一区| 国产亚洲最大av| 一二三四中文在线观看免费高清| 久久久久国产精品人妻一区二区| 国产亚洲av片在线观看秒播厂| 99国产精品免费福利视频| 亚洲一区二区三区欧美精品| 中国三级夫妇交换| 丝袜美足系列| 国产无遮挡羞羞视频在线观看| 日本欧美视频一区| 免费大片黄手机在线观看| 亚洲精品aⅴ在线观看| 国产老妇伦熟女老妇高清| 大话2 男鬼变身卡| 欧美日韩视频精品一区| 久久人人爽人人爽人人片va| 91精品伊人久久大香线蕉| 美女视频免费永久观看网站| 18+在线观看网站| 肉色欧美久久久久久久蜜桃| 亚洲熟女精品中文字幕| 亚洲综合色网址| 亚洲国产最新在线播放| 女的被弄到高潮叫床怎么办| 国产xxxxx性猛交| 亚洲成人av在线免费| 超色免费av| 亚洲内射少妇av| 一本久久精品| 如何舔出高潮| 男女无遮挡免费网站观看| 老司机影院毛片| 视频区图区小说| 黄色怎么调成土黄色| av女优亚洲男人天堂| 国产高清三级在线| 久久久久久久久久久免费av| 欧美亚洲 丝袜 人妻 在线| 五月玫瑰六月丁香| av电影中文网址| 亚洲美女搞黄在线观看| 日本黄色日本黄色录像| 男女啪啪激烈高潮av片| 一级片'在线观看视频| 日韩免费高清中文字幕av| 亚洲美女搞黄在线观看| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看 | 久久久久久久大尺度免费视频| 日韩一区二区视频免费看| 大话2 男鬼变身卡| 亚洲av男天堂| 国产免费一区二区三区四区乱码| 飞空精品影院首页| 久久精品国产亚洲av涩爱| 街头女战士在线观看网站| 久久热在线av| 高清视频免费观看一区二区| 欧美最新免费一区二区三区| 免费观看a级毛片全部| 9色porny在线观看| 日韩一区二区三区影片| 一区二区日韩欧美中文字幕 | 91aial.com中文字幕在线观看| 男女高潮啪啪啪动态图| 久久久久网色| 国产色婷婷99| 国产精品久久久久久久电影| 亚洲欧美成人综合另类久久久| 91在线精品国自产拍蜜月| 啦啦啦啦在线视频资源| 熟女av电影| 只有这里有精品99| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 看非洲黑人一级黄片| 黄片播放在线免费| 麻豆乱淫一区二区| 国产成人午夜福利电影在线观看| 亚洲av欧美aⅴ国产| 亚洲精品乱久久久久久| 精品久久久精品久久久| av线在线观看网站| 亚洲成人av在线免费| 久久久久久伊人网av| 一级毛片我不卡| av在线老鸭窝| 你懂的网址亚洲精品在线观看| 精品久久久久久电影网| 99久久综合免费| 免费黄网站久久成人精品| 伦理电影大哥的女人| 国产欧美日韩综合在线一区二区| 黄色怎么调成土黄色| 如何舔出高潮| 女人久久www免费人成看片| 精品国产一区二区久久| 成人无遮挡网站| 极品少妇高潮喷水抽搐| 国产av精品麻豆| 少妇人妻久久综合中文| 咕卡用的链子| 在线观看一区二区三区激情| 成年美女黄网站色视频大全免费| 欧美精品国产亚洲| av不卡在线播放| 精品国产一区二区三区四区第35| 免费女性裸体啪啪无遮挡网站| 日韩大片免费观看网站| 男女无遮挡免费网站观看| 美女国产高潮福利片在线看| 亚洲av在线观看美女高潮| 亚洲av电影在线进入| 日韩制服骚丝袜av| 大片免费播放器 马上看| a级毛色黄片| 久久久欧美国产精品| 街头女战士在线观看网站| 成人漫画全彩无遮挡| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 精品人妻在线不人妻| 18在线观看网站| 91在线精品国自产拍蜜月| 纯流量卡能插随身wifi吗| videosex国产| 精品午夜福利在线看| 亚洲人与动物交配视频| 人人妻人人爽人人添夜夜欢视频| 一区在线观看完整版| av国产精品久久久久影院| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 精品久久蜜臀av无| 秋霞在线观看毛片| 十八禁网站网址无遮挡| 69精品国产乱码久久久| 99久久综合免费| 精品人妻一区二区三区麻豆| 国产成人精品久久久久久| 一级a做视频免费观看| 超碰97精品在线观看| 成年人免费黄色播放视频| 精品国产一区二区久久| 一本久久精品| 国产精品人妻久久久久久| 国产熟女午夜一区二区三区| 王馨瑶露胸无遮挡在线观看| 成年女人在线观看亚洲视频| 18在线观看网站| 久久青草综合色| 国产色婷婷99| 香蕉国产在线看| 高清欧美精品videossex| 久久久精品免费免费高清| 日本午夜av视频| 国产亚洲精品久久久com| 久久毛片免费看一区二区三区| 69精品国产乱码久久久| 午夜老司机福利剧场| 成年女人在线观看亚洲视频| av在线观看视频网站免费| 精品卡一卡二卡四卡免费| 久久久精品免费免费高清| 国产精品无大码| 热re99久久国产66热| 日本91视频免费播放| 日韩av在线免费看完整版不卡| 久久精品久久精品一区二区三区| 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 夫妻午夜视频| 久久久久国产精品人妻一区二区| 91精品三级在线观看| 国产一区亚洲一区在线观看| 国产精品一区www在线观看| 精品人妻在线不人妻| 久久久久久久久久人人人人人人| 国产欧美另类精品又又久久亚洲欧美| 日本爱情动作片www.在线观看| 免费黄网站久久成人精品| 亚洲欧美一区二区三区黑人 | 亚洲综合色惰| a 毛片基地| 欧美激情极品国产一区二区三区 | 久久久久久久国产电影| kizo精华| 亚洲av.av天堂| 亚洲欧美日韩另类电影网站| 久久鲁丝午夜福利片| 日韩大片免费观看网站| 秋霞在线观看毛片| 免费高清在线观看日韩| 夫妻午夜视频| 亚洲精品国产av蜜桃| 成人毛片60女人毛片免费| 色视频在线一区二区三区| 日日啪夜夜爽| 日韩av在线免费看完整版不卡| 黄色毛片三级朝国网站| 成年动漫av网址| 国产黄频视频在线观看| 春色校园在线视频观看| 日韩在线高清观看一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产黄色视频一区二区在线观看| 九草在线视频观看| 国产激情久久老熟女| 久久久久久久国产电影| 亚洲一码二码三码区别大吗| 日本爱情动作片www.在线观看| 国产亚洲精品第一综合不卡 | 国产精品麻豆人妻色哟哟久久| 日韩 亚洲 欧美在线| 伦理电影大哥的女人| 又粗又硬又长又爽又黄的视频| 久久久精品区二区三区| 日韩成人伦理影院| 色婷婷av一区二区三区视频| 99热6这里只有精品| 国产男女超爽视频在线观看| av天堂久久9| 日本欧美视频一区| 成人黄色视频免费在线看| 久久久久精品性色| 亚洲成人av在线免费| 美女xxoo啪啪120秒动态图| 免费黄频网站在线观看国产| 国产极品天堂在线| 午夜老司机福利剧场| 亚洲国产精品999| 午夜激情久久久久久久| 欧美 日韩 精品 国产| 永久免费av网站大全| 大片免费播放器 马上看| 久久av网站| 亚洲精品色激情综合| 国产无遮挡羞羞视频在线观看| 亚洲av日韩在线播放| 精品少妇内射三级| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 女人被躁到高潮嗷嗷叫费观| 狠狠婷婷综合久久久久久88av| 国产综合精华液| 亚洲精品久久成人aⅴ小说| 日本vs欧美在线观看视频| 欧美3d第一页| 高清黄色对白视频在线免费看| 亚洲欧美日韩卡通动漫| 考比视频在线观看| 色吧在线观看| 久热久热在线精品观看| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 国产伦理片在线播放av一区| 日韩中文字幕视频在线看片| 午夜精品国产一区二区电影| 日韩av不卡免费在线播放| 精品久久久精品久久久| 免费看光身美女| 欧美国产精品一级二级三级| 丁香六月天网| 哪个播放器可以免费观看大片| 欧美精品亚洲一区二区| 久久这里有精品视频免费| 免费不卡的大黄色大毛片视频在线观看| 国产精品蜜桃在线观看| 亚洲精品乱码久久久久久按摩| 老司机影院成人| 国产成人精品在线电影| 久久综合国产亚洲精品| 精品人妻在线不人妻| 一本色道久久久久久精品综合| 久热久热在线精品观看| 亚洲精品日韩在线中文字幕| 日本欧美国产在线视频| 国产成人精品在线电影| 黄色怎么调成土黄色| 国产又色又爽无遮挡免| 日韩av在线免费看完整版不卡| 在线天堂中文资源库| 男女国产视频网站| 日本欧美视频一区| 草草在线视频免费看| 国产一级毛片在线| 午夜日本视频在线| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 中文字幕最新亚洲高清| 久久婷婷青草| 啦啦啦在线观看免费高清www| 亚洲精品一区蜜桃| 婷婷色综合大香蕉| 天堂俺去俺来也www色官网| av网站免费在线观看视频| 国产精品久久久久成人av| videosex国产| av国产精品久久久久影院| 日韩在线高清观看一区二区三区| 2018国产大陆天天弄谢| 国产亚洲av片在线观看秒播厂| 久久精品久久久久久噜噜老黄| 三级国产精品片| 人成视频在线观看免费观看| 我要看黄色一级片免费的| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 欧美+日韩+精品| 久久影院123| 赤兔流量卡办理| 免费av中文字幕在线| 男女无遮挡免费网站观看| 国产精品国产三级国产av玫瑰| 十八禁高潮呻吟视频| 69精品国产乱码久久久| 天堂俺去俺来也www色官网| 亚洲精品视频女| 下体分泌物呈黄色| a级片在线免费高清观看视频| 亚洲伊人久久精品综合| 久久久久久人人人人人| 成年美女黄网站色视频大全免费| 成人漫画全彩无遮挡| 国产毛片在线视频| 亚洲av电影在线观看一区二区三区| 日韩在线高清观看一区二区三区| 国产激情久久老熟女| 国产又爽黄色视频| 亚洲国产欧美在线一区| 新久久久久国产一级毛片| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 九九在线视频观看精品| 热re99久久国产66热| 成人亚洲精品一区在线观看| 国产成人av激情在线播放| 国产成人一区二区在线| 国产精品三级大全| 久久精品aⅴ一区二区三区四区 | 亚洲成人av在线免费| 热re99久久国产66热| av有码第一页| 久久 成人 亚洲| 日韩熟女老妇一区二区性免费视频| 黄色毛片三级朝国网站| 黄色配什么色好看| 啦啦啦视频在线资源免费观看| 亚洲精品一二三| 成人影院久久| 欧美人与善性xxx| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 综合色丁香网| 亚洲精品美女久久久久99蜜臀 | 少妇熟女欧美另类| 汤姆久久久久久久影院中文字幕| 麻豆乱淫一区二区| 街头女战士在线观看网站| 97人妻天天添夜夜摸| 丝袜在线中文字幕| 香蕉国产在线看| 亚洲国产欧美日韩在线播放| 日本av免费视频播放| 九色亚洲精品在线播放| 免费av不卡在线播放| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区综合在线观看 | 日本91视频免费播放| 一本色道久久久久久精品综合| 一个人免费看片子| 各种免费的搞黄视频| 涩涩av久久男人的天堂| 99久久中文字幕三级久久日本| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 国产乱人偷精品视频| 精品久久久久久电影网| 国产亚洲精品第一综合不卡 | 高清在线视频一区二区三区| 少妇人妻久久综合中文| 乱码一卡2卡4卡精品| 日韩中文字幕视频在线看片| 国产xxxxx性猛交| 精品久久国产蜜桃| 国产不卡av网站在线观看| 精品第一国产精品| 欧美人与善性xxx| 亚洲欧美日韩卡通动漫| 国产一区二区激情短视频 | 这个男人来自地球电影免费观看 | 69精品国产乱码久久久| 久久久久久久国产电影| 日韩伦理黄色片| 成年av动漫网址| 搡女人真爽免费视频火全软件| 男女午夜视频在线观看 | 涩涩av久久男人的天堂| 一级毛片 在线播放| 国产亚洲精品第一综合不卡 | 国产一区有黄有色的免费视频| 少妇的逼好多水| 亚洲人成网站在线观看播放| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| xxxhd国产人妻xxx| 黄色一级大片看看| 99九九在线精品视频| 国产白丝娇喘喷水9色精品| 一级毛片我不卡| 国产免费一级a男人的天堂| 欧美精品亚洲一区二区| 两性夫妻黄色片 | 看十八女毛片水多多多| 国产在视频线精品| 亚洲一码二码三码区别大吗| 麻豆乱淫一区二区| 日韩人妻精品一区2区三区| 99热全是精品| 少妇被粗大猛烈的视频| 精品亚洲成a人片在线观看| 一边亲一边摸免费视频| 少妇被粗大的猛进出69影院 | 中文字幕制服av| 亚洲国产av新网站| 日本91视频免费播放| 久久久国产一区二区| 日韩成人伦理影院| 久久 成人 亚洲| 狠狠婷婷综合久久久久久88av| 女性被躁到高潮视频| 免费人妻精品一区二区三区视频| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 日韩一区二区视频免费看| 9色porny在线观看| 90打野战视频偷拍视频| 女性被躁到高潮视频| 欧美亚洲 丝袜 人妻 在线| 内地一区二区视频在线| 在线观看三级黄色| 18禁国产床啪视频网站| 卡戴珊不雅视频在线播放| av.在线天堂| 麻豆乱淫一区二区| 国产精品秋霞免费鲁丝片| 国产又爽黄色视频|