• <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成人精品| 91精品伊人久久大香线蕉| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区久久| 亚洲中文日韩欧美视频| 精品国产国语对白av| 99九九在线精品视频| 一级毛片我不卡| 国产麻豆69| 国产欧美日韩一区二区三 | 黄色视频在线播放观看不卡| 国产精品 欧美亚洲| 国产伦人伦偷精品视频| 人人妻,人人澡人人爽秒播 | 亚洲,一卡二卡三卡| 极品人妻少妇av视频| 极品少妇高潮喷水抽搐| 高清不卡的av网站| 久久久久久久久久久久大奶| 国产亚洲av高清不卡| 亚洲精品日韩在线中文字幕| 久久久国产欧美日韩av| 久久精品亚洲熟妇少妇任你| 精品国产一区二区久久| 国产欧美日韩一区二区三 | 亚洲欧洲日产国产| 欧美日韩福利视频一区二区| 亚洲激情五月婷婷啪啪| 午夜福利,免费看| avwww免费| 久久 成人 亚洲| 久久人人爽av亚洲精品天堂| 国产男女内射视频| 日韩制服骚丝袜av| 超色免费av| 你懂的网址亚洲精品在线观看| 狠狠婷婷综合久久久久久88av| 老汉色av国产亚洲站长工具| 久久久欧美国产精品| 丁香六月天网| 在线观看人妻少妇| 精品亚洲成a人片在线观看| tube8黄色片| 欧美日韩黄片免| 麻豆乱淫一区二区| 下体分泌物呈黄色| 成年人午夜在线观看视频| 欧美日韩av久久| 中文字幕精品免费在线观看视频| 婷婷色综合www| 亚洲国产最新在线播放| 久久精品国产综合久久久| 精品卡一卡二卡四卡免费| 国产在视频线精品| 亚洲欧洲精品一区二区精品久久久| 一级毛片 在线播放| 精品国产乱码久久久久久小说| 国产免费福利视频在线观看| 国产精品一区二区在线观看99| 久久鲁丝午夜福利片| 久久精品人人爽人人爽视色| 欧美日韩综合久久久久久| 性高湖久久久久久久久免费观看| 午夜福利,免费看| 尾随美女入室| 在线精品无人区一区二区三| 啦啦啦 在线观看视频| 国产精品九九99| 老司机影院毛片| 纵有疾风起免费观看全集完整版| 日韩一本色道免费dvd| 视频在线观看一区二区三区| 亚洲一码二码三码区别大吗| 亚洲av美国av| 国产爽快片一区二区三区| 国精品久久久久久国模美| 国产欧美日韩综合在线一区二区| 狠狠婷婷综合久久久久久88av| 免费少妇av软件| 国产免费一区二区三区四区乱码| 午夜免费观看性视频| 成年女人毛片免费观看观看9 | 亚洲精品日本国产第一区| 亚洲欧美日韩高清在线视频 | 97精品久久久久久久久久精品| 色播在线永久视频| 天堂8中文在线网| 看免费av毛片| 国产精品欧美亚洲77777| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品亚洲av一区麻豆| 亚洲欧洲日产国产| av天堂久久9| 啦啦啦在线免费观看视频4| 亚洲中文日韩欧美视频| 久久精品国产综合久久久| 我要看黄色一级片免费的| 亚洲专区国产一区二区| 一二三四在线观看免费中文在| 在线观看www视频免费| 欧美成人午夜精品| 欧美日韩av久久| 国产精品久久久久久精品古装| 免费日韩欧美在线观看| av电影中文网址| 久久久久久久精品精品| av天堂久久9| 成年女人毛片免费观看观看9 | 在线观看人妻少妇| 国产一级毛片在线| 在线观看www视频免费| 麻豆乱淫一区二区| 亚洲专区国产一区二区| 亚洲欧美中文字幕日韩二区| 99国产精品一区二区三区| 中文字幕精品免费在线观看视频| 国精品久久久久久国模美| 精品一区二区三区av网在线观看 | 精品一区在线观看国产| 国产成人免费观看mmmm| 久久午夜综合久久蜜桃| 99国产精品99久久久久| 亚洲欧美中文字幕日韩二区| 国产一级毛片在线| 免费黄频网站在线观看国产| 国产精品香港三级国产av潘金莲 | 精品一区在线观看国产| 国产一区二区 视频在线| 国产高清videossex| 夫妻午夜视频| av又黄又爽大尺度在线免费看| 国产视频首页在线观看| 黄色毛片三级朝国网站| 久久久久网色| 晚上一个人看的免费电影| 伊人久久大香线蕉亚洲五| 久久久久久久国产电影| 亚洲av成人不卡在线观看播放网 | 国产爽快片一区二区三区| 日韩电影二区| 国产一区二区激情短视频 | 亚洲av日韩精品久久久久久密 | 亚洲免费av在线视频| 2018国产大陆天天弄谢| 色播在线永久视频| 亚洲中文字幕日韩| 国产亚洲精品第一综合不卡| 亚洲精品久久成人aⅴ小说| 亚洲av美国av| 久久国产精品男人的天堂亚洲| 国产一级毛片在线| 热99国产精品久久久久久7| 国产亚洲av高清不卡| 高清不卡的av网站| 午夜免费鲁丝| 国产一区二区在线观看av| 国产精品亚洲av一区麻豆| 免费高清在线观看视频在线观看| 亚洲av日韩精品久久久久久密 | 中国美女看黄片| 交换朋友夫妻互换小说| 汤姆久久久久久久影院中文字幕| 18禁黄网站禁片午夜丰满| 一级毛片电影观看| 亚洲国产最新在线播放| 性色av一级| 日韩精品免费视频一区二区三区| 欧美av亚洲av综合av国产av| 极品少妇高潮喷水抽搐| 考比视频在线观看| 男女床上黄色一级片免费看| 男女无遮挡免费网站观看| 欧美97在线视频| 亚洲国产毛片av蜜桃av| 日本欧美视频一区| 久热这里只有精品99| 国产91精品成人一区二区三区 | 精品一区在线观看国产| netflix在线观看网站| 国产成人精品久久二区二区免费| 免费一级毛片在线播放高清视频 | 亚洲美女黄色视频免费看| 欧美亚洲 丝袜 人妻 在线| 视频在线观看一区二区三区| 老司机午夜十八禁免费视频| 亚洲欧洲国产日韩| 少妇人妻 视频| 桃花免费在线播放| 成人国语在线视频| 如日韩欧美国产精品一区二区三区| 国产精品一国产av| 日韩制服骚丝袜av| 亚洲av电影在线观看一区二区三区| 五月天丁香电影| 成年av动漫网址| 1024香蕉在线观看| 国产精品久久久久久人妻精品电影 | 午夜福利视频精品| 99国产精品一区二区三区| 亚洲精品日本国产第一区| 国产老妇伦熟女老妇高清| 免费日韩欧美在线观看| 黄色 视频免费看| 日韩 欧美 亚洲 中文字幕| 午夜福利免费观看在线| 国产成人影院久久av| av网站在线播放免费| 好男人视频免费观看在线| 国产精品久久久久久精品古装| 国产日韩欧美亚洲二区| 国产精品国产av在线观看| 精品国产一区二区三区久久久樱花| 中文字幕人妻熟女乱码| svipshipincom国产片| 人妻人人澡人人爽人人| 久久久久久人人人人人| 18禁国产床啪视频网站| 母亲3免费完整高清在线观看| 国产97色在线日韩免费| av在线播放精品| 免费看av在线观看网站| 成年动漫av网址| 人体艺术视频欧美日本| 在线观看免费日韩欧美大片| 极品人妻少妇av视频| 丁香六月欧美| 秋霞在线观看毛片| 国产精品香港三级国产av潘金莲 | 亚洲人成电影观看| 久9热在线精品视频| 日韩av在线免费看完整版不卡| 桃花免费在线播放| 波多野结衣一区麻豆| 国产高清国产精品国产三级| 两人在一起打扑克的视频| 脱女人内裤的视频| 久久国产精品人妻蜜桃| 真人做人爱边吃奶动态| 国产99久久九九免费精品| 国产日韩一区二区三区精品不卡| 看免费av毛片| 9热在线视频观看99| 国产爽快片一区二区三区| 欧美老熟妇乱子伦牲交| 久久久亚洲精品成人影院| 亚洲成国产人片在线观看| 午夜日韩欧美国产| 久久久久久人人人人人| 亚洲精品一二三| 国产成人精品无人区| 一级毛片 在线播放| 大片免费播放器 马上看| 另类精品久久| 国产不卡av网站在线观看| 日本a在线网址| 中文字幕人妻熟女乱码| 久久午夜综合久久蜜桃| 9热在线视频观看99| 成人国产一区最新在线观看 | 1024视频免费在线观看| 午夜av观看不卡| 免费不卡黄色视频| 大码成人一级视频| 国产一区二区 视频在线| 国产免费视频播放在线视频| 亚洲午夜精品一区,二区,三区| 欧美日韩亚洲高清精品| 久久狼人影院| 91成人精品电影| 蜜桃在线观看..| 中文字幕人妻丝袜制服| 男女下面插进去视频免费观看| 欧美激情高清一区二区三区| 黄色一级大片看看| 久久精品国产亚洲av高清一级| 国产成人影院久久av| 亚洲免费av在线视频| 50天的宝宝边吃奶边哭怎么回事| 男男h啪啪无遮挡| 亚洲精品国产区一区二| 久久av网站| 精品国产乱码久久久久久男人| 一级片'在线观看视频| 日韩 亚洲 欧美在线| 欧美精品高潮呻吟av久久| 99精国产麻豆久久婷婷| 亚洲欧洲国产日韩| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲高清精品| 19禁男女啪啪无遮挡网站| 欧美日韩视频精品一区| 两个人看的免费小视频| 日韩中文字幕欧美一区二区 | 国产免费福利视频在线观看| 欧美精品一区二区免费开放| 性少妇av在线| 我的亚洲天堂| 老司机靠b影院| 精品少妇黑人巨大在线播放| 免费高清在线观看视频在线观看| 高清av免费在线| 丁香六月欧美| 另类精品久久| 国产成人一区二区在线| 日韩大码丰满熟妇| 多毛熟女@视频| 妹子高潮喷水视频| 老汉色av国产亚洲站长工具| 日韩一卡2卡3卡4卡2021年| 国产成人a∨麻豆精品| 国产精品久久久久久精品古装| 一本一本久久a久久精品综合妖精| 免费人妻精品一区二区三区视频| 在线观看国产h片| 久久久久久免费高清国产稀缺| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 91精品国产国语对白视频| 亚洲国产精品国产精品| 男女之事视频高清在线观看 | 亚洲欧美激情在线| 国产日韩欧美亚洲二区| 香蕉丝袜av| 天天添夜夜摸| 欧美精品一区二区大全| 在线观看免费午夜福利视频| 欧美黄色淫秽网站| av一本久久久久| 亚洲欧美一区二区三区久久| 欧美国产精品va在线观看不卡| 亚洲精品一区蜜桃| 男人添女人高潮全过程视频| 99精国产麻豆久久婷婷| 51午夜福利影视在线观看| 18禁黄网站禁片午夜丰满| 国产高清国产精品国产三级| 美女国产高潮福利片在线看| 亚洲国产毛片av蜜桃av| 99精国产麻豆久久婷婷| 自拍欧美九色日韩亚洲蝌蚪91| 国产片特级美女逼逼视频| 91麻豆精品激情在线观看国产 | 一级a爱视频在线免费观看| 亚洲精品久久成人aⅴ小说| 伊人亚洲综合成人网| www.熟女人妻精品国产| 国产亚洲一区二区精品| 国产一卡二卡三卡精品| 国产麻豆69| www日本在线高清视频| 丝袜喷水一区| 久久女婷五月综合色啪小说| 美女主播在线视频| 黄色a级毛片大全视频| 亚洲欧美一区二区三区国产| 欧美av亚洲av综合av国产av| 国产成人精品在线电影| 免费不卡黄色视频| 日韩免费高清中文字幕av| 久久久久久久国产电影| 国产不卡av网站在线观看| 永久免费av网站大全| 国产一区亚洲一区在线观看| 亚洲成人手机| 中文字幕亚洲精品专区| 亚洲精品一二三| 国产成人精品无人区| 爱豆传媒免费全集在线观看| netflix在线观看网站| 免费在线观看影片大全网站 | 久久久欧美国产精品| 欧美精品一区二区大全| 各种免费的搞黄视频| 国产成人精品久久二区二区免费| 少妇 在线观看| 亚洲成人免费av在线播放| 在线精品无人区一区二区三| 久久中文字幕一级| 汤姆久久久久久久影院中文字幕| 亚洲五月婷婷丁香| 午夜激情久久久久久久| 一边摸一边抽搐一进一出视频| 一区福利在线观看| 男女边摸边吃奶| 欧美少妇被猛烈插入视频| 妹子高潮喷水视频| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 91九色精品人成在线观看| 91字幕亚洲| 肉色欧美久久久久久久蜜桃| 亚洲国产中文字幕在线视频| 色视频在线一区二区三区| 亚洲欧美中文字幕日韩二区| 18禁观看日本| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 久久精品久久精品一区二区三区| 纯流量卡能插随身wifi吗| 成年人免费黄色播放视频| 操美女的视频在线观看| 老司机影院成人| 国产欧美日韩一区二区三 | 免费久久久久久久精品成人欧美视频| 久久久久网色| 久久久精品免费免费高清| 最新的欧美精品一区二区| 国产亚洲午夜精品一区二区久久| 欧美成人午夜精品| 国产一区二区激情短视频 | 国产成人精品久久二区二区91| 国产精品欧美亚洲77777| 日韩中文字幕欧美一区二区 | 亚洲国产欧美在线一区| 侵犯人妻中文字幕一二三四区| 日韩,欧美,国产一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品熟女久久久久浪| h视频一区二区三区| 久久九九热精品免费| 91精品国产国语对白视频| 国产爽快片一区二区三区| 成人午夜精彩视频在线观看| 黄频高清免费视频| 国产成人一区二区三区免费视频网站 | 在线观看一区二区三区激情| 国产福利在线免费观看视频| 操出白浆在线播放| 欧美人与性动交α欧美软件| 欧美+亚洲+日韩+国产| 久久久亚洲精品成人影院| 精品一区二区三区av网在线观看 | 久久久久久久国产电影| 免费高清在线观看日韩| 可以免费在线观看a视频的电影网站| 一级a爱视频在线免费观看| 好男人电影高清在线观看| 最近中文字幕2019免费版| 欧美精品人与动牲交sv欧美| 欧美精品av麻豆av| 男女午夜视频在线观看| 一区二区三区精品91| av不卡在线播放| 久久久国产精品麻豆| 搡老乐熟女国产| 欧美av亚洲av综合av国产av| 欧美精品一区二区免费开放| 王馨瑶露胸无遮挡在线观看| 欧美黑人精品巨大| 99香蕉大伊视频| 欧美 亚洲 国产 日韩一| 人妻一区二区av| 亚洲av男天堂| 亚洲精品成人av观看孕妇| 黄色一级大片看看| 各种免费的搞黄视频| 亚洲精品第二区| 女人久久www免费人成看片| 中文字幕色久视频| 一边摸一边做爽爽视频免费| 国产成人精品在线电影| 欧美少妇被猛烈插入视频| 成人国产一区最新在线观看 | 波多野结衣av一区二区av| 欧美xxⅹ黑人| 国产精品一区二区免费欧美 | 老司机靠b影院| 99国产综合亚洲精品| 国产亚洲午夜精品一区二区久久| 久久人人97超碰香蕉20202| 欧美在线一区亚洲| 中文字幕av电影在线播放| 精品一区二区三卡| 国产成人免费观看mmmm| 丰满饥渴人妻一区二区三| 国产免费福利视频在线观看| 男男h啪啪无遮挡| 日本黄色日本黄色录像| svipshipincom国产片| 亚洲成人免费av在线播放| 日本一区二区免费在线视频| 国产精品国产三级国产专区5o| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 亚洲久久久国产精品| 国产亚洲午夜精品一区二区久久| 99久久人妻综合| 欧美日韩视频精品一区| 亚洲专区国产一区二区| 免费观看人在逋| 啦啦啦中文免费视频观看日本| avwww免费| 在线观看免费日韩欧美大片| 国产片内射在线| 97精品久久久久久久久久精品| 美女中出高潮动态图| 国产精品av久久久久免费| 成年女人毛片免费观看观看9 | 欧美亚洲日本最大视频资源| a级片在线免费高清观看视频| 这个男人来自地球电影免费观看| 国产片内射在线| 在线观看www视频免费| 女性生殖器流出的白浆| 极品人妻少妇av视频| 2021少妇久久久久久久久久久| 久久狼人影院| 免费观看a级毛片全部| 好男人视频免费观看在线| 无遮挡黄片免费观看| 天堂中文最新版在线下载| √禁漫天堂资源中文www| 亚洲中文字幕日韩| 日本黄色日本黄色录像| 精品一区二区三区四区五区乱码 | 9191精品国产免费久久| 天天躁狠狠躁夜夜躁狠狠躁| 欧美精品一区二区大全| 日日夜夜操网爽| 国产淫语在线视频| 99精品久久久久人妻精品| 咕卡用的链子| 看免费成人av毛片| 色视频在线一区二区三区| 国产主播在线观看一区二区 | 乱人伦中国视频| 99精品久久久久人妻精品| 日韩一区二区三区影片| 久久女婷五月综合色啪小说| 一级毛片电影观看| 午夜免费男女啪啪视频观看| 亚洲av成人精品一二三区| 丝袜脚勾引网站| 国产成人精品无人区| 久久久久久久精品精品| 日韩视频在线欧美| 色综合欧美亚洲国产小说| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 一本大道久久a久久精品| 亚洲精品美女久久av网站| 高清av免费在线| 男女午夜视频在线观看| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区 | 侵犯人妻中文字幕一二三四区| 男女之事视频高清在线观看 | 免费黄频网站在线观看国产| 极品人妻少妇av视频| av在线老鸭窝| 成人三级做爰电影| 成人午夜精彩视频在线观看| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| 久久午夜综合久久蜜桃| 久久久欧美国产精品| 狂野欧美激情性xxxx| 免费黄频网站在线观看国产| 免费不卡黄色视频| 极品少妇高潮喷水抽搐| 国产成人91sexporn| av在线app专区| 一区二区av电影网| 久久久久久久国产电影| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av综合色区一区| 香蕉国产在线看| 国产精品国产三级国产专区5o| 18在线观看网站| 亚洲中文字幕日韩| 欧美亚洲 丝袜 人妻 在线| 一边摸一边做爽爽视频免费| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 国产欧美亚洲国产| 丁香六月欧美| 99re6热这里在线精品视频| 亚洲情色 制服丝袜| 宅男免费午夜| 两个人看的免费小视频| 国产一卡二卡三卡精品| 制服人妻中文乱码| 另类亚洲欧美激情| av一本久久久久| 国产午夜精品一二区理论片| 欧美中文综合在线视频| 午夜免费观看性视频| 欧美日韩亚洲国产一区二区在线观看 | videosex国产| 啦啦啦视频在线资源免费观看| 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 国产不卡av网站在线观看| 激情五月婷婷亚洲| 国产精品国产av在线观看| 国产高清国产精品国产三级| 国产91精品成人一区二区三区 | 青春草视频在线免费观看| 一级片'在线观看视频| 久久亚洲精品不卡| av一本久久久久| 1024香蕉在线观看| 中文字幕最新亚洲高清| 亚洲成人国产一区在线观看 | 精品亚洲成a人片在线观看| 99久久精品国产亚洲精品| tube8黄色片|