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

    基于染色質(zhì)免疫共沉淀的高通量測序數(shù)據(jù)集的順式調(diào)控模體發(fā)現(xiàn)算法

    2018-08-28 08:53:08馮艷霞張志紅張少強(qiáng)
    計(jì)算機(jī)應(yīng)用 2018年6期
    關(guān)鍵詞:模體堿基賦權(quán)

    馮艷霞,張志紅,張少強(qiáng)

    (天津師范大學(xué)計(jì)算機(jī)與信息工程學(xué)院,天津300387)

    (*通信作者電子郵箱zhangshaoqiang@tjnu.edu.cn)

    0 引言

    基于新一代基因測序的高通量技術(shù),特別是染色質(zhì)免疫共沉淀的高通量測序(ChromatinImmunoprecipitation Sequencing,ChIP-Seq)技術(shù)[1],在全基因組范圍內(nèi)的使用,很大程度上改變了生物學(xué)家定位大規(guī)模真核基因組順式調(diào)控元件和模體的方式[2]。順式調(diào)控模體是同一轉(zhuǎn)錄因子結(jié)合位點(diǎn)的序列集合。模體發(fā)現(xiàn)問題一直是生物信息學(xué)比較重要的研究課題之一。模體所蘊(yùn)含的遺傳信號(hào),對(duì)于在全基因組范圍內(nèi)研究基因的轉(zhuǎn)錄與調(diào)控機(jī)制具有重大意義。在ChIP-Seq技術(shù)之前的模體發(fā)現(xiàn),主要是通過比較同源基因或共調(diào)控基因非編碼區(qū)域的保守性來預(yù)測調(diào)控元件模體,其主要方法包括:通過多序列比對(duì)尋找保守序列片段;通過位置賦權(quán)矩陣更新來遍歷調(diào)控元件;運(yùn)用統(tǒng)計(jì)取樣的估計(jì)方法等。這些方法均取得了較好的預(yù)測結(jié)果,并應(yīng)運(yùn)而生了大量模體發(fā)現(xiàn)算法和工具,例如多重期望最大模體提取(Multiple Expectation maximization for Motif Elicitation,MEME)算法[3]、MotifSampler算法[4]、Weeder算法[5]等,為生物學(xué)家對(duì)真核生物全基因組[6]的解讀和研究提供了有利保障,推動(dòng)著生物信息學(xué)的蓬勃發(fā)展。模體發(fā)現(xiàn)算法的工作流程如圖1所示。

    由于ChIP-Seq單次實(shí)驗(yàn)產(chǎn)生了海量的堿基讀序(reads),通過序列比對(duì)拼接及結(jié)合峰預(yù)測[7]后仍有成千上萬條序列,這比傳統(tǒng)的同源基因或共調(diào)控基因序列規(guī)模大很多。為了處理如此大規(guī)模的數(shù)據(jù),近來的算法主要采用“有區(qū)別的”(discriminative)模型[8]策略,主要包括:一是選擇數(shù)據(jù)集的子集作為實(shí)驗(yàn)的輸入;二是在消耗較大數(shù)據(jù)存儲(chǔ)空間的情況下,運(yùn)用新型數(shù)據(jù)結(jié)構(gòu)來節(jié)省運(yùn)行時(shí)間;三是將實(shí)驗(yàn)數(shù)據(jù)集分割成兩個(gè)不等的部分,即實(shí)驗(yàn)組和對(duì)照組,僅將實(shí)驗(yàn)組作為模體發(fā)現(xiàn)算法的核心運(yùn)行數(shù)據(jù),然后根據(jù)實(shí)驗(yàn)組發(fā)現(xiàn)的模體,與對(duì)照組進(jìn)行比對(duì)。上述第一種方法數(shù)據(jù)損失較多,要預(yù)測的結(jié)合位點(diǎn)模體信息在子集中可能不保守,從而大大影響算法準(zhǔn)確性,例如獨(dú)立成分分析(Independent Component Analysis,ICA)和套層抽樣相結(jié)合的的模體發(fā)現(xiàn)算法——NestedMICA(Nested saMpling and ICA)[9];第二種方法的主要問題是隨著數(shù)據(jù)規(guī)模增加,會(huì)成倍增加內(nèi)存需求和時(shí)間復(fù)雜性,從而降低了算法運(yùn)行速度,例如有區(qū)別的正則表達(dá)式模體提取(Discriminative Regular Expression Motif Elicitation,DREME)算法[10];相比較而言,第三種方法更注重序列數(shù)據(jù)集的分隔方法,但不恰當(dāng)?shù)臄?shù)據(jù)分隔可能導(dǎo)致錯(cuò)誤的結(jié)果,例如Amadeus[11]、Trawler[12]。特別是,隨著 ChIP-Seq 數(shù)據(jù)的海量增長,其噪聲也越來越多,從而導(dǎo)致預(yù)測結(jié)果具有較高的錯(cuò)誤發(fā)生率(False Discovery Rate,F(xiàn)DR)。此外,模體發(fā)現(xiàn)算法在追求時(shí)間與空間的高效外,還應(yīng)準(zhǔn)確識(shí)別出模體的真實(shí)長度。

    為了解決上述問題,本文設(shè)計(jì)了一個(gè)基于ChIP-Seq數(shù)據(jù)集的順式調(diào)控模體發(fā)現(xiàn)算法——FisherNet,該算法既能提高模體預(yù)測的準(zhǔn)確率,又能提高算法運(yùn)行速度。

    圖1 模體發(fā)現(xiàn)流程Fig.1 Flow chart of motif finding

    1 數(shù)據(jù)來源與預(yù)處理

    為了檢驗(yàn)算法的優(yōu)越性,F(xiàn)isherNet算法使用與DREME[10]算法相同的數(shù)據(jù)集,包括13個(gè)小鼠胚胎干細(xì)胞(Mouse Embryonic Stem Cell,mESC)ChIP-Seq 數(shù)據(jù)集[13]、3 個(gè)小鼠紅細(xì)胞數(shù)據(jù)集和1個(gè)人類淋巴母細(xì)胞系數(shù)據(jù)集。mESC數(shù)據(jù)集為12個(gè)轉(zhuǎn)錄因子,它是維持多能性和CTCF的關(guān)鍵。小鼠紅細(xì)胞的數(shù)據(jù)集為 Gata1[14]、Klf1[15]和 Tal1,是紅細(xì)胞生成的關(guān)鍵調(diào)節(jié)者。人體淋巴細(xì)胞數(shù)據(jù)集是維生素D受體包括在使用骨化三醇的刺激前和之后的ChIP-Seq數(shù)據(jù)。另外,基于人的 ENCODE(Encyclopedia of DNA Elements)[16]數(shù)據(jù)庫的ChIP-Seq數(shù)據(jù)集規(guī)模超大,本文選取其中的10個(gè)有已知模體的轉(zhuǎn)錄因子(ATF1、ATF3、CTCF、E2F4、EBF1、ECF1、POL2、ZNF274、USF2、YY1)數(shù)據(jù)集進(jìn)行比較。這些數(shù)據(jù)集源文件格式為FASTQ格式。

    2 算法的主要步驟

    本文算法的流程如圖2所示。

    2.1 生成背景數(shù)據(jù)集

    FisherNet算法采用了三階馬爾可夫鏈[17]生成背景數(shù)據(jù),即遍歷實(shí)驗(yàn)組ChIP-Seq數(shù)據(jù)集統(tǒng)計(jì)每三個(gè)堿基短序列(共43=64種)后出現(xiàn)各個(gè)堿基的頻率得到概率轉(zhuǎn)移矩陣。然后用該轉(zhuǎn)移矩陣生成與實(shí)驗(yàn)組數(shù)據(jù)中序列長度和數(shù)目相同的背景序列數(shù)據(jù)集。

    2.2 遍歷短序生成哈希表

    為了統(tǒng)計(jì)每個(gè)長度為k的短序(k-mer)[18]在各個(gè)序列是否存在,F(xiàn)isherNet分別遍歷實(shí)驗(yàn)數(shù)據(jù)集和背景數(shù)據(jù)集中的每條序列并用兩個(gè)哈希表分別存儲(chǔ)所有的k-mer信息。哈希表的鍵為k-mer,鍵值為該k-mer所在序列編號(hào)集合。因?yàn)檎婧松锏慕Y(jié)合位點(diǎn)長度主要介于4~8個(gè)堿基長。FisherNet算法默認(rèn)k從4到8個(gè)堿基為窗口長度分別進(jìn)行遍歷。在遍歷過程中,同時(shí)考慮k-mer在序列的互補(bǔ)鏈?zhǔn)欠翊嬖?,若存在,則將該序列的編號(hào)也放到該k-mer鍵值對(duì)應(yīng)的集合中。遍歷過程如圖3所示。

    圖2 本文算法流程Fig.2 Flow chart of proposed algorithm

    2.3 根據(jù)費(fèi)舍爾精確檢驗(yàn)計(jì)算P值

    FisherNet算法采用費(fèi)舍爾精確檢驗(yàn)分析每個(gè)k-mer的富集情況[19]。對(duì)于每個(gè)k-mer,通過哈希表查詢記錄下它在實(shí)驗(yàn)數(shù)據(jù)集和背景數(shù)據(jù)集中出現(xiàn)的序列條數(shù)。如表1所示,a代表某k-mer在實(shí)驗(yàn)組數(shù)據(jù)集中出現(xiàn)過的序列數(shù),b代表該短序在背景數(shù)據(jù)集中出現(xiàn)過的序列數(shù),相對(duì)地,c和d分別代表該短序在兩組實(shí)驗(yàn)數(shù)據(jù)集中未出現(xiàn)的序列數(shù)。

    根據(jù)表1通過超幾何分布公式計(jì)算P值[20]為:

    若P值大于意義臨界值0.05(假設(shè)檢驗(yàn)常規(guī)設(shè)P值的閾值為0.05,在算法中該閾值作為輸入?yún)?shù)可由使用者自主調(diào)整),說明該條DNA短序列存在基因信息與背景差異較小,則對(duì)其不再考慮。上述計(jì)算P值的式(1),對(duì)于小規(guī)模數(shù)據(jù)集效果良好,對(duì)于大型ChIP-Seq數(shù)據(jù)集,則階乘n!會(huì)是超大的數(shù),甚至超出計(jì)算機(jī)整數(shù)計(jì)數(shù)范圍。對(duì)此,F(xiàn)isherNet引入了斯特林(Stirling)公式,對(duì)組合數(shù)的計(jì)算先取其對(duì)數(shù)ln,再取指數(shù)。斯特林如式(2)所示,與ln C推導(dǎo)如式(3)、(4)所示:

    則P值公式就等價(jià)于:

    對(duì)式(1)兩邊取對(duì)數(shù)得:

    對(duì)于每個(gè)k-mer,用式(6)計(jì)算其P值,若P值小于0.05,以此k-mer為鍵,其P值為鍵值保存到一個(gè)新的哈希表。通過比較各k-mer的鍵值,對(duì)k-mer進(jìn)行排序得到候選k-mer集合。

    表1 k-mer(模體)序列在實(shí)驗(yàn)組和背景組的數(shù)據(jù)統(tǒng)計(jì)Tab.1 Statistics of k-mer(motif)sequence in experimental group and background group

    2.4 構(gòu)建初始模體

    在候選k-mer集合中選取P值最小的k-mer作為種子,依次向下遍歷與種子長度相同且只有一個(gè)位置不同的k-mer,將其合并形成初始模體。每合并一個(gè)k-mer就要通過計(jì)算出現(xiàn)模體和未出現(xiàn)模體的序列數(shù)來更新費(fèi)舍爾精確檢驗(yàn)表并計(jì)算合并后模體的P值,若P值 >0.05,則停止k-mer合并。統(tǒng)計(jì)初始序列上各堿基的出現(xiàn)頻率得到位置頻率矩陣P=(P(b,i))4×k[21],其中 P(b,i)表示堿基 b 在位置 i出現(xiàn)的頻率;然后再根據(jù)式(7)計(jì)算出對(duì)應(yīng)的位置賦權(quán)矩陣W =(W(b,i))4×k[22],其中:P(b)是堿基 b 在實(shí)驗(yàn)數(shù)據(jù)集上出現(xiàn)的頻率(為了能夠計(jì)算對(duì)數(shù)值,當(dāng)P(b,i)=0,令P(b,i)=10-5)。構(gòu)建初始模體的流程如圖4中第① ~③步所示。

    2.5 生成最終模體

    用上述生成的位置賦權(quán)矩陣掃描候選k-mer集合中的每一個(gè)相同長度的k-mer={b1,b2,…,bk},為了衡量一個(gè)k-mer是否屬于模體,本文用式(8)(公式中的Smax、Smin分別為位置賦權(quán)矩陣中每列元素的最大值的和與最小值的和)計(jì)算k-mer的相對(duì)分?jǐn)?shù)值R并只保留R≥0.80(根據(jù)已知模體統(tǒng)計(jì)有80% 的位置具有保守性)的k-mer,以此生成最終模體。如圖4中第④到⑤步所示。

    3 實(shí)驗(yàn)與結(jié)果分析

    FisherNet算法運(yùn)用mESC的ChIP-Seq數(shù)據(jù)集和一個(gè)意義臨界值對(duì)結(jié)果進(jìn)行具體分析。如表2為FisherNet算法預(yù)測出的模體發(fā)現(xiàn)情況。將算法發(fā)現(xiàn)的模體序列與Tomtom數(shù)據(jù)庫[23]中的真實(shí)模體進(jìn)行比較,針對(duì)Tomtom標(biāo)準(zhǔn)數(shù)據(jù)庫中已有的13個(gè)重要模體,F(xiàn)isherNet算法發(fā)現(xiàn)了其中的10個(gè)。

    此外,針對(duì)mESC ChIP-Seq數(shù)據(jù)集,本文分別運(yùn)行了其他六種被廣泛應(yīng)用且能發(fā)現(xiàn)輔調(diào)控因子模體的模體發(fā)現(xiàn)算法,即 MEME[3]、Weeder[5]、NestedMICA[9]、DREME[10]、Amadeus[11]、Trawler[12],并將這些算法的預(yù)測結(jié)果與 FisherNet順式調(diào)控模體算法相比較,統(tǒng)計(jì)結(jié)果如表3所示。其中:N為各算法發(fā)現(xiàn)的模體總數(shù);S為與Tomtom標(biāo)準(zhǔn)數(shù)據(jù)庫中相匹配的模體數(shù);C為算法發(fā)現(xiàn)的輔調(diào)控因子模體數(shù)。在發(fā)現(xiàn)轉(zhuǎn)錄因子模體上FisherNet與 DREME、Weeder、NestedMICA 和 MEME 均達(dá)到最高,但在發(fā)現(xiàn)輔調(diào)控因子模體上FisherNet算法要優(yōu)于其他算法。這說明該算法在保證發(fā)現(xiàn)模體的精確度的同時(shí),能發(fā)現(xiàn)更多的輔調(diào)控因子模體,F(xiàn)isherNet總體能夠達(dá)到80%(即(6.1+10)/20)的精確度。表3最后一列為各算法對(duì)13個(gè)mESC數(shù)據(jù)集的平均運(yùn)行時(shí)間,可以看出,F(xiàn)isherNet是所有比對(duì)算法中平均運(yùn)行時(shí)間最短的,而NestedMICA和Weeder是運(yùn)行時(shí)間最長的兩個(gè)算法。

    表2 小鼠的ES細(xì)胞ChIP-Seq數(shù)據(jù)集發(fā)現(xiàn)的模體Tab.2 Found motif in mouse ES cells ChIP-Seq dataset

    表3 不同模體發(fā)現(xiàn)算法結(jié)果統(tǒng)計(jì)Tab.3 Results statistics of different motif finding algorithms

    由于DREME是MEME開發(fā)者設(shè)計(jì)的改進(jìn)算法,為了進(jìn)一步檢驗(yàn)除了MEME、Weeder和NestedMICA之外的四種算法的效率,分別取序列規(guī)模為1兆、2兆、3兆和4兆堿基對(duì)的數(shù)據(jù)集進(jìn)行運(yùn)算,圖5為四種算法的運(yùn)行時(shí)間評(píng)估折線圖。由圖5可以看出,與其他算法相比,隨著數(shù)據(jù)集序列堿基規(guī)模的增加,F(xiàn)isherNet算法運(yùn)行時(shí)間最短且增加幅度不大,充分說明了該算法在運(yùn)行效率上的優(yōu)越性。將上述算法應(yīng)用到選取的ENCODE[16]數(shù)據(jù)集中,發(fā)現(xiàn)每一個(gè)真實(shí)模體均存在于FisherNet算法輸出的前10個(gè)模體中,而其余算法的前10個(gè)模體中含有真實(shí)模體的準(zhǔn)確度均小于80%。

    圖5 不同算法運(yùn)行時(shí)間比較Fig.5 Running time comparison of different motif finding algorithms

    4 結(jié)語

    本文基于ChIP-Seq數(shù)據(jù)的模體發(fā)現(xiàn),以ChIP-Seq結(jié)合峰數(shù)據(jù)集為算法輸入數(shù)據(jù),提出一個(gè)全新的順式調(diào)控模體發(fā)現(xiàn)算法(本文的FisherNet由 Perl語言編寫,代碼下載地址為http://www.escience.cn/people/sqzhang/program.html.)。該算法運(yùn)用費(fèi)舍爾假設(shè)檢驗(yàn)來尋找區(qū)分度高的k-mer并通過位置賦權(quán)矩陣進(jìn)行優(yōu)化篩選。與其他多種模體發(fā)現(xiàn)算法進(jìn)行比較,驗(yàn)證了該算法大大提高了模體發(fā)現(xiàn)的數(shù)目和準(zhǔn)確度,成功預(yù)測到許多核心模體及輔調(diào)控因子模體。該算法目前只對(duì)幾種代表性的ChIP-Seq數(shù)據(jù)進(jìn)行了驗(yàn)證,對(duì)其他物種的數(shù)據(jù)是否具有魯棒性有待進(jìn)一步驗(yàn)證。

    猜你喜歡
    模體堿基賦權(quán)
    論鄉(xiāng)村治理的有效賦權(quán)——以A縣扶貧項(xiàng)目為例
    中國西部(2022年2期)2022-05-23 13:28:20
    基于Matrix Profile的時(shí)間序列變長模體挖掘
    企業(yè)數(shù)據(jù)賦權(quán)保護(hù)的反思與求解
    應(yīng)用思維進(jìn)階構(gòu)建模型 例談培養(yǎng)學(xué)生創(chuàng)造性思維
    中國科學(xué)家創(chuàng)建出新型糖基化酶堿基編輯器
    試論新媒體賦權(quán)
    活力(2019年15期)2019-09-25 07:22:12
    生命“字母表”迎來4名新成員
    生命“字母表”迎來4名新成員
    植入(l, d)模體發(fā)現(xiàn)若干算法的實(shí)現(xiàn)與比較
    基于改進(jìn)AHP熵博弈賦權(quán)的輸變電工程評(píng)價(jià)
    男女边吃奶边做爰视频| 性欧美人与动物交配| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在 | 国产蜜桃级精品一区二区三区| 欧美日韩国产亚洲二区| 久99久视频精品免费| 国产伦精品一区二区三区视频9| 国产成人午夜福利电影在线观看| 亚洲四区av| av又黄又爽大尺度在线免费看 | 一边摸一边抽搐一进一小说| 最新中文字幕久久久久| 色综合色国产| 精品日产1卡2卡| 欧美zozozo另类| 村上凉子中文字幕在线| 久久久a久久爽久久v久久| 日本av手机在线免费观看| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久久久久久久| 亚洲av一区综合| 在线观看av片永久免费下载| 成人美女网站在线观看视频| 亚洲最大成人手机在线| av又黄又爽大尺度在线免费看 | 麻豆国产av国片精品| 在线观看一区二区三区| 91久久精品电影网| a级一级毛片免费在线观看| 51国产日韩欧美| 亚洲精品粉嫩美女一区| 亚洲成人久久性| 熟女电影av网| 久久久久久久久久久免费av| 99久国产av精品| 色播亚洲综合网| 久久精品国产清高在天天线| 寂寞人妻少妇视频99o| 亚洲人成网站在线观看播放| 韩国av在线不卡| 午夜激情欧美在线| 免费看av在线观看网站| 国产单亲对白刺激| 免费在线观看成人毛片| 99riav亚洲国产免费| 国产片特级美女逼逼视频| 中国美女看黄片| 男女做爰动态图高潮gif福利片| 两个人的视频大全免费| 一个人看的www免费观看视频| av在线亚洲专区| 精品久久久久久久久久免费视频| 成人二区视频| 两个人视频免费观看高清| 国产男人的电影天堂91| 国产高潮美女av| 99热精品在线国产| 国产日本99.免费观看| 欧美性猛交╳xxx乱大交人| 久久人人爽人人爽人人片va| 18禁黄网站禁片免费观看直播| 亚洲成人久久性| 成人亚洲欧美一区二区av| 在线播放国产精品三级| 99久久精品国产国产毛片| 成年女人永久免费观看视频| 村上凉子中文字幕在线| 人人妻人人澡欧美一区二区| 黄色日韩在线| 日韩在线高清观看一区二区三区| 久久人人爽人人片av| or卡值多少钱| 国产精品无大码| 日本熟妇午夜| 九九爱精品视频在线观看| 国产精品伦人一区二区| 97在线视频观看| 午夜福利在线观看免费完整高清在 | 日韩高清综合在线| 自拍偷自拍亚洲精品老妇| 国产伦理片在线播放av一区 | 久久精品久久久久久久性| 中国国产av一级| 男人舔女人下体高潮全视频| 亚洲第一区二区三区不卡| 国产精品乱码一区二三区的特点| 又爽又黄a免费视频| videossex国产| 一区二区三区高清视频在线| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 国产伦理片在线播放av一区 | 国产一区二区激情短视频| 国产黄片美女视频| 国产乱人偷精品视频| 啦啦啦啦在线视频资源| 国产亚洲精品久久久com| 一级毛片我不卡| 高清毛片免费观看视频网站| 黑人高潮一二区| 国产精品福利在线免费观看| 国产伦理片在线播放av一区 | 精品少妇黑人巨大在线播放 | 内射极品少妇av片p| 国产成人一区二区在线| 欧美激情国产日韩精品一区| 夫妻性生交免费视频一级片| 看免费成人av毛片| 人妻制服诱惑在线中文字幕| 欧美日韩综合久久久久久| 国产免费一级a男人的天堂| 1000部很黄的大片| 久久99蜜桃精品久久| av在线观看视频网站免费| 黄片wwwwww| 草草在线视频免费看| 日韩制服骚丝袜av| 国产精品,欧美在线| 国产精品国产三级国产av玫瑰| 国产精品免费一区二区三区在线| 亚洲美女视频黄频| 久久99蜜桃精品久久| 六月丁香七月| 丰满的人妻完整版| 看十八女毛片水多多多| 色尼玛亚洲综合影院| 日本av手机在线免费观看| 嫩草影院精品99| 日日撸夜夜添| 97人妻精品一区二区三区麻豆| 免费av不卡在线播放| 卡戴珊不雅视频在线播放| 亚洲色图av天堂| 黄色视频,在线免费观看| 日韩大尺度精品在线看网址| 日韩欧美一区二区三区在线观看| 老司机福利观看| 男女做爰动态图高潮gif福利片| 久久久精品大字幕| 国产亚洲精品久久久com| 婷婷六月久久综合丁香| 黄色日韩在线| 亚洲美女搞黄在线观看| a级毛片免费高清观看在线播放| 精品不卡国产一区二区三区| 可以在线观看毛片的网站| 中文字幕精品亚洲无线码一区| 国产精品一二三区在线看| 亚洲av不卡在线观看| 国产综合懂色| 国产日韩欧美在线精品| 亚洲在线自拍视频| 丝袜喷水一区| 丝袜喷水一区| 九九爱精品视频在线观看| 在线观看av片永久免费下载| 伦理电影大哥的女人| 亚洲精华国产精华液的使用体验 | 久久精品久久久久久噜噜老黄 | 亚洲综合色惰| 久久国产乱子免费精品| 国产女主播在线喷水免费视频网站 | 天天一区二区日本电影三级| 亚洲人成网站高清观看| 高清毛片免费看| av天堂中文字幕网| 亚洲无线观看免费| 青春草视频在线免费观看| 久久精品国产鲁丝片午夜精品| 国产一级毛片在线| 99热全是精品| 麻豆成人午夜福利视频| 午夜老司机福利剧场| 狠狠狠狠99中文字幕| 欧美性猛交╳xxx乱大交人| 男女边吃奶边做爰视频| 天堂√8在线中文| 国产精品人妻久久久影院| 久久这里只有精品中国| 3wmmmm亚洲av在线观看| 国产真实伦视频高清在线观看| 此物有八面人人有两片| 日本免费一区二区三区高清不卡| 在线天堂最新版资源| 亚洲色图av天堂| 亚洲性久久影院| 成人综合一区亚洲| 日本一本二区三区精品| 高清毛片免费看| 欧美成人免费av一区二区三区| 夜夜夜夜夜久久久久| 51国产日韩欧美| 熟女人妻精品中文字幕| 久久久久久大精品| 久久欧美精品欧美久久欧美| 国产精品一区二区在线观看99 | 麻豆国产97在线/欧美| 亚洲国产高清在线一区二区三| 亚洲人成网站在线播| 毛片一级片免费看久久久久| 国产午夜精品久久久久久一区二区三区| 91狼人影院| 久久精品久久久久久久性| 国产成人aa在线观看| 国产精品美女特级片免费视频播放器| 在线免费十八禁| 99在线视频只有这里精品首页| 免费观看的影片在线观看| 毛片女人毛片| 禁无遮挡网站| 精品欧美国产一区二区三| 一本一本综合久久| 最近的中文字幕免费完整| 非洲黑人性xxxx精品又粗又长| 日本撒尿小便嘘嘘汇集6| 少妇熟女欧美另类| 欧美成人免费av一区二区三区| 亚洲av电影不卡..在线观看| 欧美极品一区二区三区四区| 午夜精品国产一区二区电影 | 精品久久久久久久久av| 国产av不卡久久| 久久久久久国产a免费观看| 久久久久久国产a免费观看| 天美传媒精品一区二区| 嫩草影院新地址| 欧美成人一区二区免费高清观看| 久久久久网色| 久99久视频精品免费| 女的被弄到高潮叫床怎么办| 亚洲人成网站在线播| 天堂网av新在线| 国产午夜精品一二区理论片| 亚洲在久久综合| 在线观看免费视频日本深夜| 身体一侧抽搐| 亚洲真实伦在线观看| 少妇被粗大猛烈的视频| 久久久成人免费电影| 亚洲18禁久久av| 乱系列少妇在线播放| 最近2019中文字幕mv第一页| 欧美丝袜亚洲另类| 中文字幕制服av| 级片在线观看| 成人综合一区亚洲| 久久精品国产亚洲av涩爱 | 色综合站精品国产| 中文字幕熟女人妻在线| 99国产精品一区二区蜜桃av| 美女cb高潮喷水在线观看| 可以在线观看的亚洲视频| 久久精品夜色国产| 国产一级毛片在线| 精品人妻一区二区三区麻豆| 国产在线精品亚洲第一网站| 久久久欧美国产精品| 三级男女做爰猛烈吃奶摸视频| 边亲边吃奶的免费视频| 黄片无遮挡物在线观看| 免费电影在线观看免费观看| АⅤ资源中文在线天堂| 国产黄色小视频在线观看| 99riav亚洲国产免费| 欧美xxxx性猛交bbbb| 久久草成人影院| 国产色爽女视频免费观看| 国产亚洲精品av在线| 内射极品少妇av片p| 夫妻性生交免费视频一级片| 男的添女的下面高潮视频| 免费观看a级毛片全部| 黑人高潮一二区| 国产精品三级大全| 女人十人毛片免费观看3o分钟| 男人的好看免费观看在线视频| 1024手机看黄色片| 色视频www国产| 在线播放国产精品三级| 欧美潮喷喷水| 白带黄色成豆腐渣| 麻豆久久精品国产亚洲av| 日韩欧美精品v在线| 亚洲自偷自拍三级| а√天堂www在线а√下载| 亚洲欧美日韩高清专用| 日韩 亚洲 欧美在线| 亚洲人成网站在线播| av专区在线播放| 婷婷精品国产亚洲av| 日韩精品青青久久久久久| 搡女人真爽免费视频火全软件| 亚洲熟妇中文字幕五十中出| 成人漫画全彩无遮挡| 国产在线精品亚洲第一网站| 免费黄网站久久成人精品| 天堂网av新在线| 国产白丝娇喘喷水9色精品| 久久欧美精品欧美久久欧美| 国产精品av视频在线免费观看| 国产一区二区激情短视频| 免费人成视频x8x8入口观看| 日韩 亚洲 欧美在线| 老女人水多毛片| 中国美白少妇内射xxxbb| 成年版毛片免费区| 国产成人精品久久久久久| 日产精品乱码卡一卡2卡三| 午夜a级毛片| av又黄又爽大尺度在线免费看 | 国产成人一区二区在线| 国产一级毛片七仙女欲春2| 久久国产乱子免费精品| 亚洲精品色激情综合| а√天堂www在线а√下载| 非洲黑人性xxxx精品又粗又长| 成人美女网站在线观看视频| 成人午夜精彩视频在线观看| 欧美激情久久久久久爽电影| 亚洲第一区二区三区不卡| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 乱人视频在线观看| 国产亚洲精品久久久com| 午夜福利成人在线免费观看| 91精品一卡2卡3卡4卡| 成年女人看的毛片在线观看| 国产真实伦视频高清在线观看| 国产高清视频在线观看网站| 丝袜美腿在线中文| 精品无人区乱码1区二区| 中文字幕人妻熟人妻熟丝袜美| 91精品一卡2卡3卡4卡| 免费观看在线日韩| 国产 一区精品| 亚洲一级一片aⅴ在线观看| 两个人视频免费观看高清| 亚洲国产精品成人久久小说 | 男女啪啪激烈高潮av片| 丰满人妻一区二区三区视频av| 成人午夜高清在线视频| 热99re8久久精品国产| 午夜激情欧美在线| 亚洲欧美精品综合久久99| 国产在线男女| 成人特级黄色片久久久久久久| 只有这里有精品99| 国产av不卡久久| 夜夜爽天天搞| 老司机福利观看| 中国美白少妇内射xxxbb| 亚洲中文字幕一区二区三区有码在线看| 成人三级黄色视频| 国产av不卡久久| 18禁在线无遮挡免费观看视频| 桃色一区二区三区在线观看| 欧美最黄视频在线播放免费| 欧美日韩综合久久久久久| 亚洲成人精品中文字幕电影| 亚洲内射少妇av| 青春草国产在线视频 | 欧美激情国产日韩精品一区| 在线观看av片永久免费下载| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产女主播在线喷水免费视频网站 | 精品免费久久久久久久清纯| 乱系列少妇在线播放| 可以在线观看毛片的网站| 亚洲中文字幕日韩| 99久久精品热视频| 亚洲av一区综合| 美女 人体艺术 gogo| 亚洲欧美日韩高清在线视频| 女同久久另类99精品国产91| 婷婷色综合大香蕉| 最近视频中文字幕2019在线8| 一本久久中文字幕| 国产69精品久久久久777片| 99久久精品一区二区三区| 国内精品美女久久久久久| 久久久精品欧美日韩精品| 青春草国产在线视频 | 日本撒尿小便嘘嘘汇集6| 成年av动漫网址| 一个人看视频在线观看www免费| 你懂的网址亚洲精品在线观看 | 国产精品一区二区在线观看99 | 性色avwww在线观看| 免费黄网站久久成人精品| 成人亚洲欧美一区二区av| 精品一区二区免费观看| 免费观看的影片在线观看| 18+在线观看网站| 蜜臀久久99精品久久宅男| 搡老妇女老女人老熟妇| 成人亚洲精品av一区二区| 国产极品精品免费视频能看的| 免费不卡的大黄色大毛片视频在线观看 | 日本一二三区视频观看| 国产久久久一区二区三区| 中文字幕av在线有码专区| 欧美人与善性xxx| 99久久久亚洲精品蜜臀av| 亚洲综合色惰| 舔av片在线| 午夜精品在线福利| 日韩一区二区三区影片| 国产色爽女视频免费观看| 国内精品一区二区在线观看| 永久网站在线| 美女黄网站色视频| 国产精品久久久久久av不卡| 亚洲欧美精品综合久久99| 黄色一级大片看看| 国产精品久久久久久精品电影| 亚洲国产色片| 97超碰精品成人国产| 日韩中字成人| 超碰av人人做人人爽久久| 简卡轻食公司| a级一级毛片免费在线观看| 神马国产精品三级电影在线观看| 观看美女的网站| 亚洲国产精品久久男人天堂| 你懂的网址亚洲精品在线观看 | 嫩草影院精品99| 级片在线观看| 久久人人精品亚洲av| 成人午夜高清在线视频| 精品久久久久久久久久久久久| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站| 久久久精品大字幕| 精华霜和精华液先用哪个| 在现免费观看毛片| 伦精品一区二区三区| 黑人高潮一二区| 精品一区二区三区视频在线| 亚洲欧美日韩高清在线视频| 久久精品久久久久久噜噜老黄 | 欧美极品一区二区三区四区| 久久精品国产鲁丝片午夜精品| 变态另类成人亚洲欧美熟女| 亚洲国产欧美在线一区| 久久九九热精品免费| 免费观看精品视频网站| 在线观看av片永久免费下载| 丰满乱子伦码专区| 免费在线观看成人毛片| 国产精品乱码一区二三区的特点| 日本色播在线视频| 国产精品久久久久久av不卡| 丝袜喷水一区| 两个人的视频大全免费| 久久久久久大精品| 一级av片app| 岛国毛片在线播放| 国产精品国产三级国产av玫瑰| 国产老妇伦熟女老妇高清| 美女被艹到高潮喷水动态| 2021天堂中文幕一二区在线观| 日本爱情动作片www.在线观看| 亚洲av不卡在线观看| 国产熟女欧美一区二区| 一边亲一边摸免费视频| 非洲黑人性xxxx精品又粗又长| 国产精品麻豆人妻色哟哟久久 | 国产精品人妻久久久影院| 在线观看美女被高潮喷水网站| 韩国av在线不卡| 国产精品野战在线观看| 三级男女做爰猛烈吃奶摸视频| 长腿黑丝高跟| 韩国av在线不卡| 久久久久国产网址| 久久中文看片网| 国产真实乱freesex| 在线国产一区二区在线| 国产精品无大码| 免费人成视频x8x8入口观看| 亚洲国产精品国产精品| 两个人的视频大全免费| 欧美性猛交╳xxx乱大交人| 国产91av在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美中文字幕日韩二区| 久久久欧美国产精品| 色噜噜av男人的天堂激情| 亚洲欧美清纯卡通| 久久久久久久久中文| 少妇被粗大猛烈的视频| 国产精华一区二区三区| kizo精华| 婷婷六月久久综合丁香| 亚洲欧美中文字幕日韩二区| 伦精品一区二区三区| 99久久精品热视频| 国产成人aa在线观看| 黄色配什么色好看| 亚洲18禁久久av| 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 乱人视频在线观看| 又黄又爽又刺激的免费视频.| 国产精品久久久久久亚洲av鲁大| 国产精品久久视频播放| 亚洲,欧美,日韩| 国产精品一区二区三区四区久久| 国产精品久久久久久久久免| 精品熟女少妇av免费看| 一夜夜www| 成人美女网站在线观看视频| 日韩欧美精品v在线| 美女cb高潮喷水在线观看| 亚洲欧美成人精品一区二区| 国产蜜桃级精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 色综合色国产| 国产精品三级大全| 高清毛片免费看| 大型黄色视频在线免费观看| 国产精品久久视频播放| 18+在线观看网站| 亚洲综合色惰| 亚洲精品456在线播放app| 日本-黄色视频高清免费观看| 日本免费a在线| 国产成人精品一,二区 | 麻豆成人午夜福利视频| www.av在线官网国产| 亚洲精品粉嫩美女一区| 精品99又大又爽又粗少妇毛片| 精品一区二区三区人妻视频| а√天堂www在线а√下载| 国产色婷婷99| 亚洲国产精品sss在线观看| 一区二区三区四区激情视频 | 国产成人freesex在线| 日韩欧美精品v在线| 日本av手机在线免费观看| 国产单亲对白刺激| 少妇熟女aⅴ在线视频| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 老司机福利观看| 一个人看的www免费观看视频| 99久久精品热视频| 最后的刺客免费高清国语| av在线老鸭窝| 亚洲一区二区三区色噜噜| 亚洲精品粉嫩美女一区| 青春草国产在线视频 | 啦啦啦韩国在线观看视频| 麻豆久久精品国产亚洲av| 长腿黑丝高跟| 熟女人妻精品中文字幕| 国产av麻豆久久久久久久| 日本黄色视频三级网站网址| 草草在线视频免费看| 最后的刺客免费高清国语| 国产在线精品亚洲第一网站| 国产免费一级a男人的天堂| 色哟哟·www| 99热这里只有是精品在线观看| 欧美日韩国产亚洲二区| 国产真实伦视频高清在线观看| 日本爱情动作片www.在线观看| 亚洲av不卡在线观看| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 99久久人妻综合| 美女大奶头视频| 三级毛片av免费| 亚洲中文字幕一区二区三区有码在线看| 一夜夜www| 伦理电影大哥的女人| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 亚州av有码| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 精品无人区乱码1区二区| 国产一级毛片七仙女欲春2| 嫩草影院入口| 国产毛片a区久久久久| 一个人免费在线观看电影| 人人妻人人看人人澡| 精品午夜福利在线看| 国产精品永久免费网站| 91av网一区二区| 亚洲成a人片在线一区二区| 国产三级中文精品| 亚洲国产精品成人久久小说 | 中国国产av一级| 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 色噜噜av男人的天堂激情| 亚洲精品456在线播放app| 深夜a级毛片| 99热这里只有是精品在线观看| 国产爱豆传媒在线观看| 亚洲av免费高清在线观看| 美女脱内裤让男人舔精品视频 | 乱系列少妇在线播放| 国产极品精品免费视频能看的| 日韩高清综合在线| h日本视频在线播放| 亚洲七黄色美女视频| 国产探花极品一区二区| 一级二级三级毛片免费看| 日韩成人伦理影院|