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

    層次聚類在進(jìn)化樹構(gòu)建中的應(yīng)用

    2014-09-10 09:15:56李國(guó)寶
    關(guān)鍵詞:進(jìn)化樹線粒體聚類

    李國(guó)寶,業(yè) 寧

    (南京林業(yè)大學(xué) 信息科學(xué)與技術(shù)學(xué)院,南京 210037)

    0 引言

    從達(dá)爾文時(shí)代開始,人類對(duì)生物起源的研究加以重視,各種生物進(jìn)化學(xué)說相繼被提出。由于生物進(jìn)化歷史是沒有文字記載的,后人只能通過史前生物的化石等片面信息來盡可能準(zhǔn)確的模擬生物進(jìn)化的順序,這就可能會(huì)形成錯(cuò)誤的生物進(jìn)化推斷歷史。隨著20世紀(jì)中期生物遺傳信息研究取得突破進(jìn)展,人類通過生物的遺傳物質(zhì)來研究其進(jìn)化歷史成為可能。

    現(xiàn)代生物學(xué)用生物進(jìn)化樹來描述生物之間進(jìn)化關(guān)系,兩種(或者多種)生物如果在同一層節(jié)點(diǎn),則表明該組生物進(jìn)化距離較近(即從同一祖先進(jìn)化而來的可能性較大)[1];反之,表明這些生物之間的生物差異性較大。

    生物進(jìn)化樹可以根據(jù)其是否按照進(jìn)化距離構(gòu)建來分類,這樣就有基于進(jìn)化距離構(gòu)建的方法和基于統(tǒng)計(jì)特征或者生物離散特征構(gòu)建的方法?;谶M(jìn)化距離的構(gòu)建方法主要有最近鄰法[2],UPGMA法等;基于統(tǒng)計(jì)方法的構(gòu)建主要有最大似然法(Maximum likelihood)[3];基于生物離散特征的構(gòu)建方法主要是最大簡(jiǎn)約法(Maximum Parsimony)[4]。

    UPGMA方法是基于距離的進(jìn)化樹構(gòu)建方法,該方法思想是:將兩個(gè)進(jìn)化距離最近的物種合成到一個(gè)復(fù)合物種組中,然后將新的距離矩陣中距離最小的兩個(gè)物種再次合成一個(gè)復(fù)合物種組,如此反復(fù),直到所有的物種都被聚為一棵進(jìn)化樹[5]。UPGMA方法的使用有一個(gè)前提,即假設(shè)一棵進(jìn)化樹中所有物種的進(jìn)化速率[6]是相同的。

    NJ(Neighbor Joining Method, NJ)方法是距離法建樹中比較有實(shí)用價(jià)值的方法。與UPGMA相比,NJ方法不用假設(shè)進(jìn)化樹中所有物種的進(jìn)化速率相同,因此在大多數(shù)情況下比較令人信服。該方法思想是:通過確定距離最近的成對(duì)分類物種組來使進(jìn)化樹的進(jìn)化距離之和達(dá)到最小。在進(jìn)行序列合并時(shí),不僅要滿足待合并序列進(jìn)化距離的相近,同時(shí)也要求待合并的序列與其它序列的近似距離較遠(yuǎn)。本文實(shí)驗(yàn)一中生成的進(jìn)化樹將會(huì)與NJ方法生成的進(jìn)化樹進(jìn)行比較。

    ML(Maximum Likelihood)方法于1981年被提出,該方法構(gòu)建思想基于統(tǒng)計(jì)學(xué)。在預(yù)先選擇的進(jìn)化模型下計(jì)算每一種進(jìn)化樹生成的可能性,選擇最大可能性的進(jìn)化樹即為最大似然樹。最大似然法在構(gòu)建進(jìn)化樹的準(zhǔn)確度方面很高,但是在處理大數(shù)據(jù)量時(shí)效率比較低,并且對(duì)模型的依賴比較嚴(yán)重。

    MP(Maximum Parsimony)方法依據(jù)各個(gè)位置上由一條生物序列突變成另一條生物序列所需最小數(shù)量突變來進(jìn)行比較分析和聚類樹生成的,最終的進(jìn)化樹是基于整條序列所需的突變總數(shù)的。

    1 研究思想

    層次聚類的基本思想是:通過迭代分類,把相似的樣本放在一層,直到樣本都被歸到某一層中[7]。具體的層次聚類算法分為兩種,從頂層到底層的方法和從底層到頂層的方法。

    從頂層到底層的算法思想是:(1)先把樣本當(dāng)作一個(gè)集合,從這個(gè)集合中取出2個(gè)(或者多個(gè))最為相似的樣本,形成一顆二叉(或者多叉)形狀的樹,此時(shí),集合剩下的樣本作為樹根節(jié)點(diǎn),2個(gè)(或者多個(gè))相似的節(jié)點(diǎn)作為第一層葉子節(jié)點(diǎn);(2)再?gòu)臉涓?jié)點(diǎn)集合中選取一類相似的樣本作為第二層葉子節(jié)點(diǎn),繼續(xù)擴(kuò)展剛才生成的樹;(3)重復(fù)步驟2,直到根節(jié)點(diǎn)集合中樣本數(shù)為0;(4)按照2叉樹的生成規(guī)則整理生成的進(jìn)化樹,保證結(jié)果是一棵2叉樹。

    從底層到頂層的層次聚類方法與前者描述的由頂層至底層算法相似,只是在進(jìn)化樹形成順序上是從葉子節(jié)點(diǎn)到根節(jié)點(diǎn)。

    先通過序列比對(duì)計(jì)算生物之間的進(jìn)化距離,然后運(yùn)用層次聚類方法[8]對(duì)生物樣本進(jìn)行分類,最終生成層次聚類結(jié)果以二叉樹形式表現(xiàn)。與基于距離的其他方法如NJ、UPGMA等進(jìn)化樹構(gòu)建方法以及基于統(tǒng)計(jì)的ML方法,基于生物表現(xiàn)特征的MP方法等生成的進(jìn)化樹比較后驗(yàn)證層次聚類方法的可行性,以及如何提高進(jìn)化樹的準(zhǔn)確率。

    2 序列比對(duì)計(jì)算生物進(jìn)化距離

    進(jìn)化距離的計(jì)算是通過比較DNA序列得到的。序列的比對(duì)有兩兩比對(duì)和多重比對(duì)之分。比對(duì)算法有blast,clustal,fasta等,本實(shí)驗(yàn)采用的是多序列比對(duì)中的clustalw方法。

    實(shí)驗(yàn)一:8個(gè)物種線粒體DNA距離矩陣

    由于線粒體DNA變異速率很慢,它們以每一百萬(wàn)年百分之二點(diǎn)二的速度變異,因此選擇線粒體DNA作為研究對(duì)象是比較合理的。

    本實(shí)驗(yàn)將甘藍(lán)型油菜種子萌發(fā)出現(xiàn)胚根第1天的幼苗轉(zhuǎn)移至含有eBL(epi-Brassinolide)或BRZ的固體1/2 MS培養(yǎng)基上(方形皿)。eBL生理處理實(shí)驗(yàn)在光照條件下進(jìn)行,6種處理濃度分別為0、1 × 10-9、1 × 10-8、1 × 10-7、1 × 10-6和 1 × 10-5 mol/L,處理4 d后觀察不同濃度eBL對(duì)幼苗發(fā)育的影響。BRZ生理處理實(shí)驗(yàn)在黑暗下進(jìn)行,分為0和1 × 10-6 mol/L 兩個(gè)處理濃度,處理3 d后觀察BRZ對(duì)幼苗黑暗條件下發(fā)育的影響。

    序列來自GenBank,詳細(xì)信息如表1所示。

    表1 8 條線粒體 DNA 物種 GenBank 相關(guān)信息

    由表1構(gòu)建的距離矩陣見表2。

    表2 8 個(gè)物種所對(duì)應(yīng)的相似距離矩陣

    實(shí)驗(yàn)二:十一種脊椎動(dòng)物序列數(shù)據(jù)的實(shí)驗(yàn)

    生物進(jìn)化樹的好壞,主要看與它與生物真實(shí)的進(jìn)化歷史差別多大。生物的真實(shí)進(jìn)化史往往難以知曉,但是也有例外。表3給出的生物序列信息的真實(shí)進(jìn)化樹可以用古生物學(xué)和形態(tài)學(xué)方面的數(shù)據(jù)來構(gòu)建[9]。11種脊椎動(dòng)物線粒體全DNA信息如表3。

    表3 11 種脊椎動(dòng)物的名稱及其線粒體全基因組 GenBank 編號(hào)

    由表3構(gòu)建的距離矩陣如表4所示。

    表4 11種脊椎動(dòng)物線粒體DNA的距離矩陣

    續(xù)表4 11種脊椎動(dòng)物線粒體DNA的距離矩陣

    3 層次聚類生成生物進(jìn)化樹

    層次聚類的步驟:

    (1)從n個(gè)序列樣本中選擇2個(gè)距離最小的放到集合s中,此時(shí)集合s和剩下的n-2個(gè)樣本作為n-1個(gè)候選樣本;

    (2)從剩下的n-1個(gè)候選樣本中選擇2個(gè)距離最小的放到集合s1中,類似步驟1,此時(shí)n-2個(gè)候選樣本;

    (3)重復(fù)步驟2,直到候選樣本樹為0;

    (4)此時(shí)n個(gè)序列樣本層次聚類完成,通過層次號(hào)來逐層構(gòu)建二叉樹。

    實(shí)驗(yàn)一通過層次聚類方法構(gòu)建的進(jìn)化樹見圖1。

    圖1 8個(gè)物種層次聚類結(jié)果

    對(duì)生成的聚類結(jié)果統(tǒng)一注釋后得到的進(jìn)化樹見圖2。

    圖2 實(shí)驗(yàn)一的層次聚類生物進(jìn)化樹

    實(shí)驗(yàn)二通過層次聚類的結(jié)果見圖3。

    圖3 實(shí)驗(yàn)二11種脊椎動(dòng)物線粒體

    實(shí)驗(yàn)二層次聚類結(jié)果經(jīng)過解釋后的生物進(jìn)化樹見圖4。

    圖4 實(shí)驗(yàn)二11種脊椎動(dòng)物的生物進(jìn)化樹

    4 對(duì)層次聚類生物進(jìn)化樹的評(píng)價(jià)

    由于生物的真實(shí)進(jìn)化歷史無(wú)法得知,后人只能通過生物化石等信息來推斷生物最有可能形成的進(jìn)化歷史。因此,對(duì)于構(gòu)建的生物進(jìn)化樹,不存在唯一的評(píng)價(jià)標(biāo)準(zhǔn)。不過,可以在確保數(shù)據(jù)準(zhǔn)確的前提下,對(duì)多種方法構(gòu)建的進(jìn)化樹進(jìn)行比較分析,如果多種方法構(gòu)建的進(jìn)化樹拓?fù)浣Y(jié)構(gòu)一致,則可以側(cè)面驗(yàn)證構(gòu)建方法的有效性。

    實(shí)驗(yàn)一中層次聚類構(gòu)建的進(jìn)化樹與最大似然法、NJ法構(gòu)建的進(jìn)化樹(見圖5)結(jié)構(gòu)完全一致,因此可以驗(yàn)證層次聚類在構(gòu)建這8種生物線粒體DNA進(jìn)化樹中的有效性。

    圖5 最大似然法和最近鄰法構(gòu)建的8種生物進(jìn)化樹

    實(shí)驗(yàn)二中層次聚類構(gòu)建的11種脊椎動(dòng)物的進(jìn)化樹和生物真實(shí)的進(jìn)化樹拓?fù)浣Y(jié)構(gòu)有細(xì)微差別,見圖6。

    圖6 11種脊椎動(dòng)物的真實(shí)系統(tǒng)樹

    經(jīng)過查閱相關(guān)資料,作者發(fā)現(xiàn)這種進(jìn)化樹結(jié)構(gòu)差別源于所選取的序列長(zhǎng)度較長(zhǎng)。在計(jì)算相似距離時(shí),序列越長(zhǎng),出現(xiàn)程序性錯(cuò)誤的可能性就越大。

    因此,本文的結(jié)論:在層次聚類產(chǎn)生進(jìn)化樹時(shí),盡可能選取長(zhǎng)度較短的能夠代替生物進(jìn)化信息的DNA序列,這樣可以生成較準(zhǔn)確的進(jìn)化樹。

    參考文獻(xiàn):

    [1] SNEATH P, SOKAL R. Numerical taxonomy—the principles and practice of numerical classification[M].San Francisco:W. H. Freeman and Company, 1973.

    [2] SAITOU N, NEI M. The neighbor-joining method: a new method for reconstructing phylogenetic trees[J].Mol Biol Evol, 1987(4):406-425.

    [3] Felsentein J. Evolutionary trees from DNA sequences: a Maximum likelihood approach[J].J Mol Evol, 1981,17:368-376.

    [4] Kimura M. Evolutionary rate at the molecular level[J].Nature,1968,217:624-626.

    [5] 劉金桂.分?jǐn)?shù)階超混沌系統(tǒng)的自適應(yīng)函數(shù)投影同步[J].淮陰工學(xué)院學(xué)報(bào),2012,21(1):1-4.

    [6] Zhaxybayeva, O. and W.F. Doolittle, Lateral gene transfer[J].Current Biology, 2011,21(7):242-246.

    [7] 孫亂,陸祖宏,謝建明.生物信息學(xué)基礎(chǔ)[M].北京:清華大學(xué)出版社,2005.

    [8] 丁淑妍.進(jìn)化分析與結(jié)構(gòu)預(yù)測(cè)中的若干問題研究[D].大連:大連理工大學(xué),2012.

    [9] 孫士保,秦克云.改進(jìn)的k-平均聚類算法研究[J].計(jì)算機(jī)工程,2007,33(13):200-201.

    猜你喜歡
    進(jìn)化樹線粒體聚類
    基于心理旋轉(zhuǎn)的小學(xué)生物進(jìn)化樹教學(xué)實(shí)驗(yàn)報(bào)告
    棘皮動(dòng)物線粒體基因組研究進(jìn)展
    常見的進(jìn)化樹錯(cuò)誤概念及其辨析*
    線粒體自噬與帕金森病的研究進(jìn)展
    基于DBSACN聚類算法的XML文檔聚類
    福州2009—2014年甲型H1N1流感病毒株HA基因進(jìn)化分析
    艾草白粉病的病原菌鑒定
    基于改進(jìn)的遺傳算法的模糊聚類算法
    一種層次初始的聚類個(gè)數(shù)自適應(yīng)的聚類方法研究
    NF-κB介導(dǎo)線粒體依賴的神經(jīng)細(xì)胞凋亡途徑
    国产一区二区在线av高清观看| 国产私拍福利视频在线观看| 一本精品99久久精品77| 女人十人毛片免费观看3o分钟| 久久久久免费精品人妻一区二区| 成人无遮挡网站| 日韩欧美精品v在线| 国产色婷婷99| 欧美xxxx性猛交bbbb| 深夜精品福利| 亚洲成人av在线免费| 日韩大尺度精品在线看网址| 一本一本综合久久| 国产黄a三级三级三级人| 国语自产精品视频在线第100页| 韩国av在线不卡| 九色成人免费人妻av| 男人舔奶头视频| 欧美绝顶高潮抽搐喷水| 亚洲成a人片在线一区二区| 久久久精品欧美日韩精品| 蜜臀久久99精品久久宅男| 在线国产一区二区在线| 亚洲婷婷狠狠爱综合网| 色吧在线观看| 亚洲国产高清在线一区二区三| 国产 一区 欧美 日韩| 男女视频在线观看网站免费| 床上黄色一级片| 日韩精品青青久久久久久| 黄色日韩在线| 亚洲欧美日韩高清在线视频| 麻豆乱淫一区二区| 亚洲aⅴ乱码一区二区在线播放| 黄片wwwwww| 精品福利观看| 男人狂女人下面高潮的视频| 亚洲精品久久国产高清桃花| 日韩精品青青久久久久久| 亚洲av成人精品一区久久| 午夜老司机福利剧场| 久久精品久久久久久噜噜老黄 | 精品久久久久久久久久免费视频| 亚洲av免费高清在线观看| 亚洲美女视频黄频| 色哟哟哟哟哟哟| 一进一出抽搐gif免费好疼| 悠悠久久av| 日韩在线高清观看一区二区三区| 99热全是精品| 午夜精品在线福利| 欧美日韩国产亚洲二区| 亚洲最大成人手机在线| 亚洲av第一区精品v没综合| 亚洲自偷自拍三级| 欧美一区二区精品小视频在线| 精品福利观看| 97超碰精品成人国产| 国产高清三级在线| 99热这里只有精品一区| 99热全是精品| 在现免费观看毛片| 一个人观看的视频www高清免费观看| 久久精品影院6| eeuss影院久久| av在线蜜桃| 亚洲国产精品合色在线| 成人欧美大片| 亚洲天堂国产精品一区在线| 国产日本99.免费观看| 少妇猛男粗大的猛烈进出视频 | 亚洲成a人片在线一区二区| 精品欧美国产一区二区三| 欧美在线一区亚洲| 久久久精品94久久精品| 99久久中文字幕三级久久日本| 免费无遮挡裸体视频| 亚洲一区二区三区色噜噜| 亚洲人成网站在线观看播放| av在线亚洲专区| 国产男人的电影天堂91| 夜夜看夜夜爽夜夜摸| 欧美激情国产日韩精品一区| 中文字幕av在线有码专区| 日韩大尺度精品在线看网址| 丝袜喷水一区| 啦啦啦啦在线视频资源| 国产美女午夜福利| 最近在线观看免费完整版| 91在线精品国自产拍蜜月| 日韩三级伦理在线观看| 日产精品乱码卡一卡2卡三| 永久网站在线| 亚洲,欧美,日韩| 日韩在线高清观看一区二区三区| 亚洲在线自拍视频| 欧美中文日本在线观看视频| 久久国内精品自在自线图片| 丝袜美腿在线中文| av在线天堂中文字幕| 欧美一区二区国产精品久久精品| 在线免费观看的www视频| 日本 av在线| 少妇丰满av| 91av网一区二区| 神马国产精品三级电影在线观看| .国产精品久久| 精品少妇黑人巨大在线播放 | 婷婷精品国产亚洲av在线| 全区人妻精品视频| 日本熟妇午夜| 中文字幕av在线有码专区| 国产色婷婷99| 精品久久久噜噜| 色尼玛亚洲综合影院| 91精品国产九色| 国产片特级美女逼逼视频| 毛片女人毛片| 亚洲五月天丁香| 99久久无色码亚洲精品果冻| 日日摸夜夜添夜夜爱| 精品免费久久久久久久清纯| 成年女人永久免费观看视频| 99久久中文字幕三级久久日本| 最新中文字幕久久久久| eeuss影院久久| 日韩欧美精品v在线| 国产精品无大码| 欧美+亚洲+日韩+国产| 99久久久亚洲精品蜜臀av| 亚洲欧美精品自产自拍| 婷婷精品国产亚洲av在线| www日本黄色视频网| 色在线成人网| 亚洲一区二区三区色噜噜| 国产精品一区二区三区四区免费观看 | 亚洲中文字幕一区二区三区有码在线看| 日本一本二区三区精品| 久久精品国产99精品国产亚洲性色| 97在线视频观看| 在线观看一区二区三区| 日韩欧美 国产精品| 日韩精品中文字幕看吧| 国产在线男女| 国产色爽女视频免费观看| 亚洲精品久久国产高清桃花| 国产av不卡久久| 日本与韩国留学比较| 国语自产精品视频在线第100页| 久久久久久久亚洲中文字幕| 欧美色视频一区免费| 欧美性猛交╳xxx乱大交人| 一级a爱片免费观看的视频| 国产女主播在线喷水免费视频网站 | 黑人高潮一二区| 男女下面进入的视频免费午夜| 悠悠久久av| 国产极品精品免费视频能看的| 少妇熟女aⅴ在线视频| 黄色日韩在线| 成人高潮视频无遮挡免费网站| 搞女人的毛片| 日韩,欧美,国产一区二区三区 | 91久久精品电影网| 欧美三级亚洲精品| 成人亚洲欧美一区二区av| 欧美最新免费一区二区三区| 免费看光身美女| 伦精品一区二区三区| 成人永久免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品99久久久久久久久| 伦精品一区二区三区| 婷婷亚洲欧美| 国产免费一级a男人的天堂| 久久热精品热| 在线观看av片永久免费下载| 精品人妻偷拍中文字幕| 此物有八面人人有两片| 永久网站在线| 亚洲第一区二区三区不卡| 国产精品一区二区性色av| 两个人的视频大全免费| 国产爱豆传媒在线观看| 97碰自拍视频| 自拍偷自拍亚洲精品老妇| 国产精品无大码| 18禁黄网站禁片免费观看直播| ponron亚洲| 国产高潮美女av| 熟女电影av网| 两个人的视频大全免费| 小说图片视频综合网站| 婷婷精品国产亚洲av| 欧美区成人在线视频| 中文亚洲av片在线观看爽| 网址你懂的国产日韩在线| 久久久久免费精品人妻一区二区| 国产欧美日韩精品一区二区| 免费观看精品视频网站| 男女那种视频在线观看| 亚洲在线自拍视频| 小蜜桃在线观看免费完整版高清| 伦精品一区二区三区| 三级经典国产精品| 99视频精品全部免费 在线| 国语自产精品视频在线第100页| 精品福利观看| av在线蜜桃| 日本色播在线视频| 久久99热这里只有精品18| 黑人高潮一二区| 成年女人永久免费观看视频| 最近中文字幕高清免费大全6| 亚洲婷婷狠狠爱综合网| 毛片女人毛片| 麻豆乱淫一区二区| 国产精品一二三区在线看| 国产久久久一区二区三区| 亚洲在线自拍视频| 国产伦在线观看视频一区| 不卡视频在线观看欧美| 大型黄色视频在线免费观看| 日本精品一区二区三区蜜桃| АⅤ资源中文在线天堂| 变态另类丝袜制服| 日日摸夜夜添夜夜爱| 免费观看的影片在线观看| 少妇人妻精品综合一区二区 | 亚洲成人精品中文字幕电影| 国产精品久久久久久亚洲av鲁大| 最好的美女福利视频网| 欧美日韩在线观看h| 搡老岳熟女国产| 免费不卡的大黄色大毛片视频在线观看 | 91在线观看av| 国产精品99久久久久久久久| 国产精品亚洲一级av第二区| 欧美一区二区精品小视频在线| 成人三级黄色视频| 久久久久久久久大av| 亚洲国产欧洲综合997久久,| 日本爱情动作片www.在线观看 | 久久鲁丝午夜福利片| 欧美日韩一区二区视频在线观看视频在线 | 级片在线观看| 亚洲人与动物交配视频| 欧美xxxx性猛交bbbb| 免费大片18禁| 在线观看午夜福利视频| 一进一出抽搐动态| 国产色爽女视频免费观看| 久久久欧美国产精品| 国产白丝娇喘喷水9色精品| 国产高清视频在线播放一区| 久久99热这里只有精品18| 久久久精品94久久精品| 男女视频在线观看网站免费| 能在线免费观看的黄片| 国产久久久一区二区三区| 午夜视频国产福利| 成人亚洲欧美一区二区av| 亚洲精品日韩在线中文字幕 | 久久精品人妻少妇| 日韩中字成人| 亚洲精品粉嫩美女一区| 亚洲精品影视一区二区三区av| 99热全是精品| 丝袜喷水一区| 久久久久国产精品人妻aⅴ院| 久久久久久九九精品二区国产| 99久久精品热视频| 一个人免费在线观看电影| 国产精品嫩草影院av在线观看| 欧美3d第一页| 啦啦啦观看免费观看视频高清| 免费无遮挡裸体视频| 亚洲精品亚洲一区二区| 99久久精品热视频| 嫩草影院精品99| 少妇猛男粗大的猛烈进出视频 | 伊人久久精品亚洲午夜| 亚洲av成人av| 亚洲精品一卡2卡三卡4卡5卡| 国产老妇女一区| 亚洲av中文av极速乱| 久久亚洲精品不卡| 国产一区二区在线av高清观看| 日本免费a在线| 午夜精品在线福利| 少妇熟女欧美另类| 老女人水多毛片| 亚洲成av人片在线播放无| 亚洲精品一卡2卡三卡4卡5卡| 日本免费一区二区三区高清不卡| 一级a爱片免费观看的视频| 国产精品,欧美在线| 观看美女的网站| 国产亚洲av嫩草精品影院| 人妻久久中文字幕网| 亚洲av成人精品一区久久| 黄色配什么色好看| 三级国产精品欧美在线观看| 日本 av在线| 美女xxoo啪啪120秒动态图| 给我免费播放毛片高清在线观看| 啦啦啦啦在线视频资源| 国产精品一区二区免费欧美| 一个人看视频在线观看www免费| 在线天堂最新版资源| 久久人人精品亚洲av| 欧美精品国产亚洲| 中国美女看黄片| 免费看av在线观看网站| 国产色爽女视频免费观看| 欧美激情久久久久久爽电影| 午夜福利18| 国产欧美日韩精品亚洲av| 国产成人影院久久av| 国产精品久久电影中文字幕| 免费无遮挡裸体视频| 国产精品久久久久久精品电影| 舔av片在线| 男女视频在线观看网站免费| av女优亚洲男人天堂| 国产日本99.免费观看| 十八禁国产超污无遮挡网站| 欧美人与善性xxx| 久久精品国产亚洲网站| 国产在视频线在精品| 国产午夜精品论理片| 国产亚洲精品久久久久久毛片| 精品人妻视频免费看| 综合色av麻豆| 亚洲国产色片| 乱码一卡2卡4卡精品| 伦理电影大哥的女人| 亚洲自拍偷在线| 五月玫瑰六月丁香| 免费搜索国产男女视频| 成人鲁丝片一二三区免费| 一级毛片aaaaaa免费看小| 欧美最新免费一区二区三区| 日产精品乱码卡一卡2卡三| 在线免费观看的www视频| 欧美xxxx黑人xx丫x性爽| 别揉我奶头 嗯啊视频| 大香蕉久久网| 最新中文字幕久久久久| 亚洲av第一区精品v没综合| 国产精品一二三区在线看| 91狼人影院| 搡老岳熟女国产| 99国产精品一区二区蜜桃av| 欧美极品一区二区三区四区| 国产免费男女视频| 亚洲自拍偷在线| 国产欧美日韩精品亚洲av| 狠狠狠狠99中文字幕| 中文字幕精品亚洲无线码一区| 女的被弄到高潮叫床怎么办| 伊人久久精品亚洲午夜| www日本黄色视频网| 亚洲七黄色美女视频| 成熟少妇高潮喷水视频| 国产精品亚洲美女久久久| 美女免费视频网站| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 女同久久另类99精品国产91| 国产蜜桃级精品一区二区三区| 99久久无色码亚洲精品果冻| 久久精品国产亚洲av天美| 免费无遮挡裸体视频| 一进一出好大好爽视频| 男人的好看免费观看在线视频| 老女人水多毛片| 12—13女人毛片做爰片一| 亚洲最大成人av| 99热这里只有是精品在线观看| 国产精品亚洲美女久久久| 3wmmmm亚洲av在线观看| 欧美成人一区二区免费高清观看| 国产成人a区在线观看| 国产乱人视频| 美女xxoo啪啪120秒动态图| 精品午夜福利视频在线观看一区| 看十八女毛片水多多多| 日日撸夜夜添| 俺也久久电影网| 亚洲三级黄色毛片| 男女啪啪激烈高潮av片| 深爱激情五月婷婷| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 有码 亚洲区| 国产日本99.免费观看| 又黄又爽又免费观看的视频| 国产伦精品一区二区三区视频9| 亚洲四区av| 欧美绝顶高潮抽搐喷水| 午夜a级毛片| 欧美成人一区二区免费高清观看| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 亚洲欧美成人综合另类久久久 | 啦啦啦韩国在线观看视频| 久久精品人妻少妇| 综合色丁香网| 国产91av在线免费观看| 免费无遮挡裸体视频| 亚洲欧美成人综合另类久久久 | 欧美日本视频| 国产精品无大码| 美女内射精品一级片tv| 精品一区二区三区av网在线观看| 日本免费一区二区三区高清不卡| 18禁在线播放成人免费| 国产精品久久久久久久久免| av福利片在线观看| 国产欧美日韩精品亚洲av| 亚洲国产色片| av在线亚洲专区| 精品久久久久久久久亚洲| 高清毛片免费观看视频网站| 亚洲av中文字字幕乱码综合| 久久国产乱子免费精品| 日韩高清综合在线| 天堂av国产一区二区熟女人妻| 美女免费视频网站| 亚洲四区av| 在线国产一区二区在线| 色噜噜av男人的天堂激情| 男插女下体视频免费在线播放| 联通29元200g的流量卡| ponron亚洲| 超碰av人人做人人爽久久| 性欧美人与动物交配| 波野结衣二区三区在线| 女生性感内裤真人,穿戴方法视频| 久久精品国产99精品国产亚洲性色| 免费看光身美女| 久久精品久久久久久噜噜老黄 | 两性午夜刺激爽爽歪歪视频在线观看| 搡女人真爽免费视频火全软件 | 麻豆av噜噜一区二区三区| 色综合色国产| 欧美性猛交黑人性爽| 午夜精品一区二区三区免费看| 欧美+亚洲+日韩+国产| 69av精品久久久久久| 日本在线视频免费播放| 国内精品美女久久久久久| or卡值多少钱| 日本三级黄在线观看| 日韩欧美 国产精品| 搡老妇女老女人老熟妇| 亚洲av一区综合| 国产精品伦人一区二区| 成人av一区二区三区在线看| 精品少妇黑人巨大在线播放 | 精品午夜福利视频在线观看一区| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 啦啦啦啦在线视频资源| 欧美日韩一区二区视频在线观看视频在线 | 日韩中字成人| 老司机影院成人| av黄色大香蕉| 在线观看免费视频日本深夜| 亚洲精品色激情综合| 禁无遮挡网站| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 精品不卡国产一区二区三区| 免费看美女性在线毛片视频| 乱码一卡2卡4卡精品| 国产精品无大码| 亚洲欧美日韩高清专用| 黄色一级大片看看| 哪里可以看免费的av片| 一级av片app| 你懂的网址亚洲精品在线观看 | 午夜视频国产福利| 色综合亚洲欧美另类图片| 色视频www国产| 精品一区二区三区人妻视频| av视频在线观看入口| 赤兔流量卡办理| 99久久中文字幕三级久久日本| 最新在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 国产成人a∨麻豆精品| 亚洲成av人片在线播放无| 亚洲18禁久久av| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 一级毛片我不卡| 伊人久久精品亚洲午夜| 精品久久国产蜜桃| 精品久久久久久久久久久久久| 自拍偷自拍亚洲精品老妇| 精品欧美国产一区二区三| 亚洲色图av天堂| 久久久精品94久久精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲自拍偷在线| 欧美色视频一区免费| 亚洲经典国产精华液单| 18禁在线播放成人免费| av专区在线播放| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 午夜精品在线福利| 在线a可以看的网站| 观看美女的网站| 少妇的逼水好多| 亚洲自偷自拍三级| 高清毛片免费看| 国产真实伦视频高清在线观看| 成人性生交大片免费视频hd| 99视频精品全部免费 在线| 91在线精品国自产拍蜜月| 少妇的逼水好多| 日本精品一区二区三区蜜桃| 九九爱精品视频在线观看| 欧美日本亚洲视频在线播放| 亚洲国产色片| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 免费人成在线观看视频色| 免费av毛片视频| 中出人妻视频一区二区| 日产精品乱码卡一卡2卡三| 国产淫片久久久久久久久| 精品不卡国产一区二区三区| 岛国在线免费视频观看| 欧美人与善性xxx| 露出奶头的视频| 国产色婷婷99| 久久中文看片网| 亚洲国产精品久久男人天堂| 精华霜和精华液先用哪个| 美女高潮的动态| 午夜福利在线在线| 男人和女人高潮做爰伦理| 亚洲av二区三区四区| 欧美激情在线99| 成人高潮视频无遮挡免费网站| 中文字幕熟女人妻在线| 十八禁国产超污无遮挡网站| 免费av毛片视频| 午夜福利成人在线免费观看| 亚洲欧美日韩无卡精品| 日韩 亚洲 欧美在线| 国产女主播在线喷水免费视频网站 | 日韩一本色道免费dvd| 1024手机看黄色片| av.在线天堂| 久久天躁狠狠躁夜夜2o2o| 99久久九九国产精品国产免费| 亚洲av一区综合| 91狼人影院| 国产精品女同一区二区软件| 免费av观看视频| 色综合站精品国产| 看片在线看免费视频| 久久人人爽人人爽人人片va| 成人综合一区亚洲| 久久久久久久久久黄片| 精品久久久久久久久久久久久| 国产高清激情床上av| 亚洲人成网站在线播| 欧美日韩精品成人综合77777| 男人狂女人下面高潮的视频| 欧美日韩在线观看h| 日本黄色片子视频| 免费看a级黄色片| a级毛色黄片| 最好的美女福利视频网| 国产麻豆成人av免费视频| 久久久午夜欧美精品| 99久国产av精品| 国产精品无大码| 麻豆精品久久久久久蜜桃| 麻豆久久精品国产亚洲av| 99久国产av精品国产电影| 日本黄大片高清| 亚洲国产精品成人久久小说 | 人妻制服诱惑在线中文字幕| 俄罗斯特黄特色一大片| 成人欧美大片| 成人特级黄色片久久久久久久| 色在线成人网| 在现免费观看毛片| 欧美绝顶高潮抽搐喷水| 久久婷婷人人爽人人干人人爱| 国产精品美女特级片免费视频播放器| 日韩成人伦理影院| 久久久久久伊人网av| 亚洲国产色片| 欧美潮喷喷水| 午夜福利在线在线| 亚洲av免费高清在线观看| 寂寞人妻少妇视频99o| 简卡轻食公司| 国产精品久久久久久亚洲av鲁大| 99热6这里只有精品| 波多野结衣高清无吗| 美女高潮的动态| 18禁黄网站禁片免费观看直播| 午夜精品国产一区二区电影 | 俺也久久电影网|