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

    微生物二元網(wǎng)絡(luò)作用關(guān)系研究

    2015-07-28 12:40高紅艷
    現(xiàn)代電子技術(shù) 2015年11期
    關(guān)鍵詞:信息論生物信息學(xué)

    高紅艷

    摘 要: 基因測(cè)序技術(shù)的快速發(fā)展產(chǎn)生了大量海洋微生物數(shù)據(jù),使難培養(yǎng)的海洋微生物種群研究成為可能。很多學(xué)者在海洋微生物種群作用模式多樣性、組成多樣性及功能多樣性方面展開(kāi)了相應(yīng)研究,有助于大家認(rèn)識(shí)海洋微生物的分布特征及其在海洋生態(tài)系統(tǒng)中的功能和所扮演的角色。然而,海洋微生物在生長(zhǎng)過(guò)程中不僅受到其他微生物的影響,也受到海洋環(huán)境(溫度,鹽度等)的影響,而在這方面的研究進(jìn)展很小。通過(guò)信息論的方法構(gòu)建微生物與環(huán)境因子的二元網(wǎng)絡(luò),并對(duì)其進(jìn)行社團(tuán)挖掘,從環(huán)境的角度去發(fā)現(xiàn)海洋微生物種群作用模式的多樣性。

    關(guān)鍵詞: 生物信息學(xué); 網(wǎng)絡(luò)構(gòu)建; 信息論; 社團(tuán)挖掘

    中圖分類(lèi)號(hào): TN99?34 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1004?373X(2015)11?0091?04

    Research on interaction relationship in microorganism binary networks

    GAO Hong?yan

    (Institute of Physics and Optoelectronics Technology, Baoji University of Arts and Science, Baoji 721016, China)

    Abstract: With the rapid development of gene sequencing technology, a large number of marine microorganism data is produced, which makes the research of more hard cultivating marine microbial population possible. The corresponding researched in diversity of interaction mode, constitute, function for marine microbial population were conducted by many scholars, which is helpful to understand the distribution features of marine microorganism, and functions and characters in marine ecosystem. The growth process of marine microorganism is not only influenced by other microbes, but also influenced by marine environment (temperature, salinity, etc.). However, the research progress in this field is little. The binary networks of microorganism and environment factors are constructed by information theory, the community detection for binary networks is proceeded. The diversity of interaction mode in marine microbial population is discovered in the view of environment.

    Keywords: bioinformatics; network construction; information theory; community detection

    0 引 言

    海洋微生物是海洋生態(tài)系統(tǒng)中的重要組成部分,對(duì)海洋微生物種群作用模式多樣性、組成多樣性及功能多樣性開(kāi)展研究,將有助于人們認(rèn)識(shí)海洋微生物的分布特征及其在海洋生態(tài)系統(tǒng)中的功能和所扮演的角色,對(duì)深入開(kāi)展海洋生態(tài)環(huán)境研究具有重要的意義[1?2]?,F(xiàn)代微生物學(xué)研究表明,99%以上的海洋微生物目前不能(或很難)被純培養(yǎng),也就是說(shuō)目前對(duì)海洋微生物的認(rèn)知只占其總量[1]的l%。近幾年,隨著微生物基因組學(xué)、PCR和高通量、低成本測(cè)序技術(shù)的發(fā)展,環(huán)境基因組學(xué)推動(dòng)了海洋微生物多樣性研究的快速發(fā)展。運(yùn)用環(huán)境基因組學(xué)技術(shù)(如16S/18S rRNA基因序列分析)研究海洋微生物生態(tài)系統(tǒng),不僅可以獲得有關(guān)海洋微生物生理多樣性和生物功能的詳細(xì)信息,還有助于了解微生物如何響應(yīng)環(huán)境脅迫,并將其應(yīng)用于海洋病害防治。

    16SrRNA基因是細(xì)菌染色體上編碼rRNA相對(duì)應(yīng)的DNA序列,存在于所有細(xì)菌的染色體基因組中,由于其高信息量、高保守性及與大多數(shù)生理、遺傳標(biāo)記一致,16SrRNA基因檢測(cè)技術(shù)已經(jīng)成為微生物種群多樣性及作用模式研究的有力工具。如Hewson等在澳大利亞莫頓灣的八個(gè)采樣點(diǎn)上取樣本[3],采用Pyrosequencing技術(shù)產(chǎn)生細(xì)菌浮游生物16S/23S rRNA序列,在不同序列相似度閾值水平上生成微生物操作分類(lèi)單元(Operation Tax Unit,OTU),通過(guò)統(tǒng)計(jì)OTU在不同采樣點(diǎn)的微生物豐度值,研究不同地點(diǎn)的海洋微生物種群分布特性;Huhrman和Martiny等人采用遺傳指紋圖譜技術(shù)[4],基于判別分析、時(shí)間序列分析、多元回歸分析等多種統(tǒng)計(jì)方法對(duì)美國(guó)加州南部海域微生物的時(shí)間和空間分布模式進(jìn)行了研究;Andersson等在波羅的海海域一個(gè)采樣位點(diǎn)上進(jìn)行8次表層海水采樣[5],時(shí)間跨度長(zhǎng)達(dá)兩年,采用Pyrosequencing技術(shù)產(chǎn)生269 420條細(xì)菌浮游生物16SrRNA序列,在97%序列相似性水平下聚類(lèi)生成4 624個(gè)OTU,基于Spearman秩相關(guān)性分析了OTU與環(huán)境因子的作用模式;Gilbert等人在2003~2008年間從西英海峽觀(guān)測(cè)站采樣點(diǎn)采集76個(gè)表層海水樣本[6],采用Pyrosequencing技術(shù)產(chǎn)生968 140條16SrRNA序列,采用α?多樣性計(jì)算、判別分析等統(tǒng)計(jì)方法分析了不同季節(jié)下的海洋微生物分布特性,發(fā)現(xiàn)了豐度值最高、分布最廣泛的微生物種群,并且得出了微生物種群在不同季節(jié)下的分布差異較大,而不同年份間分布差異較小的結(jié)論。上述這些研究為分析海洋微生物種群組成多樣性、物種間復(fù)雜聯(lián)系多樣性及生物功能提供了基礎(chǔ)。然而這些研究大多基于統(tǒng)計(jì)學(xué)方法,且沒(méi)有考慮海洋環(huán)境因子的影響,無(wú)法深程度地挖掘隱含在這些大量數(shù)據(jù)中的海洋微生物作用模式。

    一般地,同一個(gè)生態(tài)系統(tǒng)下的各個(gè)物種與環(huán)境因子(Environment,ENV)之間也存在相關(guān)關(guān)系,可以通過(guò)種群數(shù)量隨時(shí)間的變化趨勢(shì)來(lái)判斷,當(dāng)物種和環(huán)境之間的數(shù)量呈同步變化時(shí),它們可能存在依賴(lài)關(guān)系;也有可能存在非線(xiàn)性關(guān)系,有時(shí)從曲線(xiàn)變化趨勢(shì)中很難得到。本文以季節(jié)為時(shí)間單位,分別構(gòu)建春、夏、秋、冬四個(gè)季節(jié)下的OTU?ENV二元作用網(wǎng)絡(luò),通過(guò)二元網(wǎng)絡(luò)挖掘方法研究海洋微生物的作用模式。

    1 數(shù)據(jù)集和理論方法

    1.1 數(shù)據(jù)集

    本文16SrRNA數(shù)據(jù)和環(huán)境因子來(lái)自文獻(xiàn)[6],在2003年1月—2008年12月共采集76個(gè)樣本,通過(guò)多種子策略的啟發(fā)式454序列聚類(lèi)方法,獲得969 400條16SrRN序列,這些序列可從VAMPS(http://vamps.mbl.edu/index.php)上下載。其春(4~6月)、夏(7~9月)、秋(10~12月)、冬(1~3月)四個(gè)季節(jié)的采樣數(shù)分別為24、21、15和16,四個(gè)季節(jié)相應(yīng)的16SrRNA序列數(shù)分別為276 932、247 907、212 921和231 640條。環(huán)境因子共有18個(gè),也是76個(gè)時(shí)間點(diǎn),與16SrRNA的采樣時(shí)間點(diǎn)一一對(duì)應(yīng),但有的時(shí)間點(diǎn)有缺失值,本文通過(guò)取近鄰平均值來(lái)補(bǔ)全缺失值,得到18個(gè)ENV。

    生物學(xué)研究表明,序列相似性越高的物種在生物進(jìn)化樹(shù)上越相近,一般認(rèn)為相似度大于99%的序列在微生物分類(lèi)學(xué)上屬于同一物種。本文采用多種子策略基于啟發(fā)式的序列聚類(lèi)算法[7],該算法具有良好的魯棒性能,與傳統(tǒng)的啟發(fā)式序列聚類(lèi)算法相比,該算法能降低OTU過(guò)估計(jì)問(wèn)題,提高聚類(lèi)精度,有效地進(jìn)行操作分類(lèi)單元計(jì)算。在99%序列相似性閾值下對(duì)968 140條16SrRNA序列進(jìn)行聚類(lèi),可生成6 793個(gè)OTUs。

    1.2 理論方法

    信息論:相對(duì)于皮爾遜相關(guān)系數(shù)而言,互信息考慮了向量間的非線(xiàn)性關(guān)系[8]。在介紹互信息之前,先介紹一下信息熵。一個(gè)隨機(jī)向量[X={x1,x2,…,xN},][D=][max(X)-min(X)]。將[D]分成[M]個(gè)區(qū)間,得到一個(gè)區(qū)間集[{Dj}, j=1,2,…,M。]滿(mǎn)足[?j{Dj}=D]和[Dj?Dk=?,]如果[j≠k,][M]是將[D]分成的區(qū)間個(gè)數(shù)。定義[δ]函數(shù):

    [δ(Δxi,Dj)=1,if Δxi∈Dj0,else, i=1,2,…,N;j=1,2,…,M]

    區(qū)間[Dj]的概率為:

    [pX(Dj)=1Ni=1Nδ(Xi,Dj), j=1,2,…,M] (1)

    根據(jù)香農(nóng)理論,變量[X]的信息熵是:

    [H(X)=-j=1MpX(Dj)logpX(Dj)] (2)

    根據(jù)上文提到的[δ]函數(shù),現(xiàn)有兩個(gè)隨機(jī)向量[X={x1,x2,…,xN}]和[Y={y1,y2,…,yN}。]定義聯(lián)合[δ]函數(shù):

    [δ(Δxi,Δyi,Dj,Dk)=1, if Δxi∈Dj and Δyi∈Dk0, else, i=1,2,…,N;j,k=1,2,…,M]

    聯(lián)合區(qū)間[{Dj,Dk}]的概率是:

    [pX,Y(Dj,Dk)=1Ni=1Nδ(Δxi,Δyi,Dj,Dk), j,k=1,2,…,M]

    類(lèi)似的,可定義聯(lián)合熵[H(X,Y)]:

    [H(X,Y)=-j=1Mk=1MpX,Y(Dj,Dk)logpX,Y(Dj,Dk)]

    聯(lián)合熵不大于各自的熵之和,即:

    [H(X,Y)≤H(X)+H(Y)] (3)

    如果[X]和[Y]在統(tǒng)計(jì)上是獨(dú)立的,則有:

    [H(X,Y)=H(X)+H(Y)] (4)

    兩個(gè)向量間的互信息定義為:

    [MI(X,Y)=H(X)+H(Y)-H(X,Y)] (5)

    根據(jù)等式(3),有:[MI(X,Y)≥0。]根據(jù)等式(4)可知:如果[X]和[Y]在統(tǒng)計(jì)上是獨(dú)立的,則它們的互信息是0。綜上所述,由公式(5)得到的互信息可以用來(lái)度量?jī)蓚€(gè)OTU之間或OTU與ENV之間的依賴(lài)關(guān)系,不管是線(xiàn)性的還是非線(xiàn)性的。

    LP(Label Propagation)算法,目前,二元網(wǎng)絡(luò)的研究通常有兩種思路:第一種把二元網(wǎng)絡(luò)投影到單頂點(diǎn)網(wǎng)絡(luò),然后進(jìn)行網(wǎng)絡(luò)分析。單頂點(diǎn)網(wǎng)絡(luò)無(wú)法完整描述出原始二分網(wǎng)絡(luò)的全部信息,即造成了原始二分網(wǎng)絡(luò)信息的丟失。第二種是直接基于原始二分網(wǎng)絡(luò)進(jìn)行分析,最大程度上保留了原始二分網(wǎng)絡(luò)的信息。Newman將科學(xué)家合作網(wǎng)投影到單頂點(diǎn)網(wǎng)絡(luò)進(jìn)行了研究[9],網(wǎng)絡(luò)的連邊是帶權(quán)重的。Lambiotte和Ausloos提出了一種聚類(lèi)方法——PIB(percolation Idea?Based),并在聽(tīng)眾與歌曲網(wǎng)中加以實(shí)現(xiàn)[10],得到聽(tīng)眾節(jié)點(diǎn)的社團(tuán)結(jié)構(gòu)和歌曲的社團(tuán)結(jié)構(gòu),但無(wú)法得出兩類(lèi)節(jié)點(diǎn)之間的聯(lián)系。Barber在二元網(wǎng)絡(luò)模塊度[Q]基礎(chǔ)上去挖掘二元網(wǎng)絡(luò)[11],當(dāng)[Q]值最大時(shí)得到最佳社團(tuán)結(jié)構(gòu),但其方法得需先知道分團(tuán)數(shù)。本文采用Raghavan等人提出的LP算法[12]去分析海洋微生物網(wǎng)絡(luò),LP算法思想簡(jiǎn)單,復(fù)雜度小,準(zhǔn)確度也高。LP算法思想如下:

    (1) 初始化,二元網(wǎng)絡(luò)中每個(gè)節(jié)點(diǎn)都有一個(gè)單獨(dú)的標(biāo)簽,即開(kāi)始時(shí)都屬于一個(gè)單獨(dú)的團(tuán)。

    (2) 更新每個(gè)節(jié)點(diǎn)的標(biāo)簽。使每個(gè)節(jié)點(diǎn)的標(biāo)簽等于其鄰居節(jié)點(diǎn)標(biāo)簽個(gè)數(shù)最大的標(biāo)簽。如果其鄰居節(jié)點(diǎn)標(biāo)簽個(gè)數(shù)最大的標(biāo)簽不止一個(gè),則隨機(jī)選取。

    (3) 重復(fù)步驟(2),直到每個(gè)節(jié)點(diǎn)的標(biāo)簽都等于其鄰居節(jié)點(diǎn)標(biāo)簽個(gè)數(shù)最大的標(biāo)簽。

    2 實(shí)驗(yàn)結(jié)果仿真分析

    為了研究四個(gè)季節(jié)微生物與環(huán)境因子之間的相互作用,用互信息構(gòu)建海洋微生物相關(guān)作用網(wǎng)絡(luò),得到四個(gè)季節(jié)下的OTU?ENV作用網(wǎng)絡(luò)。春季網(wǎng)絡(luò)有29個(gè)節(jié)點(diǎn),22條邊;夏季網(wǎng)絡(luò)有45個(gè)節(jié)點(diǎn),45條邊;秋季網(wǎng)絡(luò)有106個(gè)節(jié)點(diǎn),159條邊;冬季網(wǎng)絡(luò)有66個(gè)節(jié)點(diǎn),64條邊。其中秋季的作用網(wǎng)絡(luò)如圖1所示,圓點(diǎn)表示OTU,三角形表示ENV。由此可知,春、夏、冬三個(gè)季節(jié)的網(wǎng)絡(luò)圖比較簡(jiǎn)單,有的就是一個(gè)單獨(dú)的團(tuán),秋季網(wǎng)絡(luò)圖節(jié)點(diǎn)比較多,相對(duì)來(lái)說(shuō)復(fù)雜一點(diǎn)。

    用LP算法直接在四個(gè)季節(jié)網(wǎng)絡(luò)中進(jìn)行挖掘,秋季網(wǎng)絡(luò)的挖掘結(jié)果如圖2所示。每個(gè)季節(jié)挖掘出的團(tuán)數(shù)如表1所示。實(shí)驗(yàn)結(jié)果顯示,環(huán)境因子E3、E5和E12在四個(gè)季節(jié)里均出現(xiàn),說(shuō)明這三個(gè)環(huán)境因子在微生物的生長(zhǎng)階段有一定的影響。E3、E5和E12分別是DX1、PAR和NO2+NO3。

    文獻(xiàn)[6]也構(gòu)建了微生物和環(huán)境因子的網(wǎng)絡(luò),可以看出E12(NO2+NO3)在微生物的四季生長(zhǎng)中起著重要作用。本文通過(guò)二元網(wǎng)絡(luò)的挖掘,也發(fā)現(xiàn)E12對(duì)海洋微生物的重要影響,所以重點(diǎn)分析E12。生物分類(lèi)學(xué)家用域、界、門(mén)、綱、目、科、屬、種七個(gè)類(lèi)對(duì)生物進(jìn)行分類(lèi),本文參照此分類(lèi)對(duì)各個(gè)季節(jié)的網(wǎng)絡(luò)進(jìn)行了分析。其中秋季網(wǎng)絡(luò)中有7個(gè)OTUs和E12相連,7個(gè)OTUs全是變形菌門(mén)(Proteobacteria)。在綱級(jí)下,5個(gè)OTUs是α?變形菌綱(Alphaproteobacteria),1個(gè)OTU是γ?變形菌綱(Gammaproteobacteria),1個(gè)OTU是δ?變形菌綱(Deltaproteobacteria)。在目級(jí)下,有3個(gè)OTUs是立克次體目(Rickettsiales),1個(gè)OTU是紅螺菌目(Rhodospirillales),1個(gè)OTU是交替單胞菌目(Alteromonadales),1個(gè)OTU是脫硫桿菌目(Desulfobacterales),1個(gè)OTU是海洋螺菌目(Oceanospirillales)。在科級(jí)下,2個(gè)OTUs是SAR11群,1個(gè)OTU是紅螺菌科(Rhodospirillaceae),1個(gè)OTU是交替單胞菌科(Alteromonadaceae),1個(gè)OTU是Nitrospinaceae科,其他2個(gè)OTU未知。

    通過(guò)四個(gè)季節(jié)對(duì)E12(NO2+NO3)所在團(tuán)的分析,得出E12(NO2+NO3)環(huán)境因子在變形菌的生長(zhǎng)過(guò)程中起著一定的作用,而且在不同的季節(jié)里,變形菌和擬桿菌共同出現(xiàn)的可能性很大,說(shuō)明這兩種菌類(lèi)在生長(zhǎng)過(guò)程中相互影響,究竟是共生、寄生、捕食和競(jìng)爭(zhēng)哪種關(guān)系,需要進(jìn)一步的研究。

    3 結(jié) 語(yǔ)

    挖掘海洋微生物的多樣性對(duì)海洋資源的開(kāi)發(fā)起著關(guān)鍵的作用。運(yùn)用宏基因組技術(shù)研究海洋微生物,不僅可以獲得有關(guān)海洋微生物的作用模式、生理多樣性和生物功能的詳細(xì)信息,還有助于了解微生物如何響應(yīng)環(huán)境脅迫。本文采用多種子策略基于啟發(fā)式的序列聚類(lèi)算法將16SrRNA基因序列聚類(lèi)成OTU,用微生物豐度向量表示OTU單元,基于互信息定義OTU單元與環(huán)境因子間是否存在網(wǎng)絡(luò)連接邊,分別構(gòu)建春、夏、秋、冬四季微生物二元作用網(wǎng)絡(luò),進(jìn)而采用二元網(wǎng)絡(luò)LP算法進(jìn)行挖掘。挖掘結(jié)果分析表明:春、夏、秋、冬海洋微生物與環(huán)境因子間作用網(wǎng)絡(luò)滿(mǎn)足復(fù)雜網(wǎng)絡(luò)“小世界”和“無(wú)尺度”特性;四個(gè)季節(jié)下的海洋微生物與環(huán)境因子的作用模式存在一定的差異,且四個(gè)季節(jié)下海洋微生物間的作用模式也存在一定差異。

    參考文獻(xiàn)

    [1] KARL D M. Microbial oceanography: paradigms, processes and promise [J]. Nature Reviews Microbiology, 2007, 5(10): 759?767.

    [2] FALKOWSKI P G, VARGAS C D. Shotgun sequencing in the sea: a blast from the past? [J]. Science, 2004, 304(5667): 58?60.

    [3] HEWSON I, FUHRMAN J A. Richness and diversity of bacterioplankton species along an estuarine gradient in Moreton Bay, Australia [J]. Applied and Environmental Microbiology, 2004, 70(6): 3425?3433.

    [4] FUHRMAN J A, STEELE J A, HEWSON I, et al. A latitudinal diversity gradient in planktonic marine bacteria [J]. Proceedings of the National Academy of Sciences of the United States of America, 2008, 105(22): 7774?7778.

    [5] ANDERSON A F, RIEMANN L, BERTILSSON S. Pyrosequencing reveals contrasting seasonal dynamics of taxa within Baltic Sea bacterioplankton communities [J]. The ISME Journal, 2009, 4(2): 171?181.

    [6] GILBERT J A, STEELE J A, CAPORASO J G, et al. Defining seasonal marine microbial community dynamics [J]. The ISME journal, 2011, 6(2): 298?308.

    [7] CHEN W, CHENG Y, ZHANG C, et al. MSClust: A multi?seeds based clustering algorithm for microbiome profiling using 16S rRNA sequence [J]. Journal of Microbiological Methods, 2013, 94(3): 347?355.

    [8] COVER T M, THOMAS J A. Elements of information theory [M]. New York: Wiley, 1991.

    [9] NEWMAN M E. Scientific collaboration networks. I. Network construction and fundamental results [J]. Physical Review E, 2001, 64(1): 543?551.

    [10] LAMBIOTTE R, AUSLOOS M. Uncovering collective listening habits and music genres in bipartite networks [J]. Physical Review E, 2005, 72(6): 1102?1108.

    [11] BARBER M J. Modularity and community detection in bipartite networks [J]. Physical Review E, 2007, 76(6): 326?331.

    [12] RAGHAVAN U N, ALBERT R, KUMARA S. Near linear time algorithm to detect community structures in large?scale networks [J]. Physical Review E, 2007, 76(3): 426?432.

    猜你喜歡
    信息論生物信息學(xué)
    信息論與應(yīng)用數(shù)學(xué)關(guān)聯(lián)教學(xué)法研究
    基于超像素和信息論的SAR圖像目標(biāo)檢測(cè)研究
    從信息論角度探析漢譯英過(guò)程中文化信息的冗余平衡
    淺談醫(yī)學(xué)院校生物信息學(xué)專(zhuān)業(yè)青年教師規(guī)范培訓(xùn)模式的建立
    “PBL+E—learning”教學(xué)模式探索
    移動(dòng)教學(xué)在生物信息學(xué)課程改革中的應(yīng)用
    中醫(yī)大數(shù)據(jù)下生物信息學(xué)的發(fā)展及教育模式淺析
    生物信息學(xué)課堂危機(jī)及對(duì)策研究
    安全通論(11)——《信息論》、《博弈論》與《安全通論》的融合:刷新您的通信觀(guān)念
    信息論翻譯的多維度探索
    国产亚洲精品av在线| 国产在线男女| 国语自产精品视频在线第100页| 国产老妇女一区| 亚洲欧美精品专区久久| 婷婷色综合大香蕉| 一级二级三级毛片免费看| 精品人妻熟女av久视频| 免费人成在线观看视频色| av福利片在线观看| 在线观看一区二区三区| 精品免费久久久久久久清纯| 国产精品一二三区在线看| 乱码一卡2卡4卡精品| 久久亚洲国产成人精品v| 男的添女的下面高潮视频| 久久婷婷人人爽人人干人人爱| 国内揄拍国产精品人妻在线| 波多野结衣高清作品| 夜夜爽天天搞| 日韩 亚洲 欧美在线| 午夜激情福利司机影院| 国产乱人偷精品视频| 中文在线观看免费www的网站| 黄色配什么色好看| 亚洲精品成人久久久久久| 亚洲最大成人av| 亚洲国产精品成人综合色| 99热6这里只有精品| 国产成人精品一,二区 | 长腿黑丝高跟| av女优亚洲男人天堂| 中文在线观看免费www的网站| 亚洲色图av天堂| 91av网一区二区| 欧美性猛交╳xxx乱大交人| 最近的中文字幕免费完整| 成人美女网站在线观看视频| 免费观看在线日韩| 黄片wwwwww| 一个人免费在线观看电影| 国产精品一区www在线观看| 久久久久久久亚洲中文字幕| 亚洲国产高清在线一区二区三| 色噜噜av男人的天堂激情| 中国国产av一级| 18禁在线无遮挡免费观看视频| 一级av片app| 午夜精品一区二区三区免费看| 你懂的网址亚洲精品在线观看 | 日本爱情动作片www.在线观看| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| ponron亚洲| 亚洲av成人av| 国产淫片久久久久久久久| 国产又黄又爽又无遮挡在线| 国产亚洲精品久久久久久毛片| 日产精品乱码卡一卡2卡三| 狂野欧美激情性xxxx在线观看| 午夜爱爱视频在线播放| 美女被艹到高潮喷水动态| 国语自产精品视频在线第100页| 国产亚洲av片在线观看秒播厂 | 身体一侧抽搐| 国产精品一区二区三区四区免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av电影不卡..在线观看| 成年av动漫网址| 亚洲欧美精品综合久久99| 性欧美人与动物交配| 精品久久久久久成人av| 联通29元200g的流量卡| 日韩,欧美,国产一区二区三区 | 波多野结衣巨乳人妻| 欧美一区二区精品小视频在线| 国产一区亚洲一区在线观看| 日韩亚洲欧美综合| 精品久久久噜噜| 国产午夜福利久久久久久| av在线蜜桃| 亚洲色图av天堂| 亚洲三级黄色毛片| 十八禁国产超污无遮挡网站| 亚洲久久久久久中文字幕| 久久99热6这里只有精品| 能在线免费看毛片的网站| 国产亚洲av嫩草精品影院| 能在线免费观看的黄片| 秋霞在线观看毛片| 亚洲四区av| 在线观看66精品国产| 亚洲av中文字字幕乱码综合| 亚洲国产精品合色在线| 国产熟女欧美一区二区| 日韩一本色道免费dvd| 国产精品av视频在线免费观看| 久久精品综合一区二区三区| 久久综合国产亚洲精品| 国产精品嫩草影院av在线观看| 午夜激情欧美在线| 国产亚洲欧美98| 日本一本二区三区精品| 精品99又大又爽又粗少妇毛片| 99热这里只有是精品在线观看| 能在线免费看毛片的网站| 日本黄色片子视频| 国内精品久久久久精免费| 日韩国内少妇激情av| 人人妻人人澡人人爽人人夜夜 | 只有这里有精品99| 精品熟女少妇av免费看| 久久热精品热| 国产精品综合久久久久久久免费| 看非洲黑人一级黄片| 日韩中字成人| 亚洲av男天堂| 色综合亚洲欧美另类图片| 一级毛片电影观看 | 欧美人与善性xxx| 亚洲色图av天堂| 久久国产乱子免费精品| 超碰av人人做人人爽久久| 久久国内精品自在自线图片| 国产69精品久久久久777片| 精品久久久久久久久久免费视频| 国产精品.久久久| 亚洲成人久久性| 边亲边吃奶的免费视频| 91aial.com中文字幕在线观看| АⅤ资源中文在线天堂| 欧美色视频一区免费| 成熟少妇高潮喷水视频| 亚洲在久久综合| 日产精品乱码卡一卡2卡三| 亚洲自偷自拍三级| 六月丁香七月| 久久午夜亚洲精品久久| 国产高清三级在线| 亚洲经典国产精华液单| 成人毛片a级毛片在线播放| 国产精品无大码| 久久精品影院6| 亚洲四区av| 国产大屁股一区二区在线视频| 国产欧美日韩精品一区二区| av.在线天堂| 免费观看a级毛片全部| 久久精品国产亚洲av香蕉五月| 我要搜黄色片| 国产黄片视频在线免费观看| 在线免费观看不下载黄p国产| 国产探花在线观看一区二区| 在线a可以看的网站| 亚洲欧美精品专区久久| 狂野欧美激情性xxxx在线观看| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 麻豆成人av视频| 午夜福利高清视频| 人人妻人人澡欧美一区二区| 少妇人妻一区二区三区视频| 色哟哟·www| 免费搜索国产男女视频| 一级av片app| 久久精品国产亚洲av涩爱 | 亚洲自偷自拍三级| 一个人观看的视频www高清免费观看| 亚洲人与动物交配视频| 赤兔流量卡办理| 精品一区二区三区人妻视频| 久久99热这里只有精品18| or卡值多少钱| 成人av在线播放网站| 久久久久久久亚洲中文字幕| 美女脱内裤让男人舔精品视频 | 3wmmmm亚洲av在线观看| 狠狠狠狠99中文字幕| 中出人妻视频一区二区| 精品少妇黑人巨大在线播放 | 欧美日韩国产亚洲二区| 成年av动漫网址| 99久久精品一区二区三区| 欧美区成人在线视频| 久久99蜜桃精品久久| 中文精品一卡2卡3卡4更新| 国产精品精品国产色婷婷| 性插视频无遮挡在线免费观看| 国产成人a∨麻豆精品| 校园人妻丝袜中文字幕| 亚洲内射少妇av| 天天一区二区日本电影三级| 亚洲在久久综合| 亚洲国产高清在线一区二区三| 欧美日韩一区二区视频在线观看视频在线 | 欧美精品一区二区大全| 精品少妇黑人巨大在线播放 | 久久久色成人| 两个人视频免费观看高清| av女优亚洲男人天堂| 熟女人妻精品中文字幕| 国产老妇伦熟女老妇高清| 又粗又爽又猛毛片免费看| 精品免费久久久久久久清纯| 国产精品综合久久久久久久免费| 亚洲国产日韩欧美精品在线观看| 丝袜喷水一区| 欧美成人免费av一区二区三区| 国产亚洲欧美98| 男人的好看免费观看在线视频| 天天躁日日操中文字幕| 久久精品综合一区二区三区| 国产蜜桃级精品一区二区三区| 国产一级毛片七仙女欲春2| 国产精品人妻久久久久久| 欧美成人免费av一区二区三区| 久久久国产成人免费| 久久久色成人| 嫩草影院入口| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 99久国产av精品| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产色片| 午夜老司机福利剧场| 美女 人体艺术 gogo| 国产精品一二三区在线看| 国产乱人偷精品视频| 中文字幕av在线有码专区| 亚洲人成网站在线播放欧美日韩| 欧美最新免费一区二区三区| 男女啪啪激烈高潮av片| 少妇人妻一区二区三区视频| 成人国产麻豆网| 嫩草影院入口| 99热这里只有是精品在线观看| 国产日韩欧美在线精品| 麻豆成人av视频| 成年女人看的毛片在线观看| 亚洲欧美精品自产自拍| 夜夜看夜夜爽夜夜摸| 22中文网久久字幕| 婷婷六月久久综合丁香| 午夜免费男女啪啪视频观看| 少妇的逼水好多| а√天堂www在线а√下载| 欧美成人精品欧美一级黄| 国产91av在线免费观看| 国产伦在线观看视频一区| 亚洲中文字幕日韩| 日本五十路高清| 免费不卡的大黄色大毛片视频在线观看 | 男人的好看免费观看在线视频| 桃色一区二区三区在线观看| 啦啦啦韩国在线观看视频| 国产精品1区2区在线观看.| 欧美又色又爽又黄视频| 深爱激情五月婷婷| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 国产午夜精品久久久久久一区二区三区| av视频在线观看入口| 天堂av国产一区二区熟女人妻| 亚洲精品成人久久久久久| a级毛片免费高清观看在线播放| 亚洲国产精品国产精品| 我的老师免费观看完整版| 床上黄色一级片| 一本久久中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 少妇熟女欧美另类| 91精品国产九色| 日韩大尺度精品在线看网址| 国产精品一区二区性色av| eeuss影院久久| 深爱激情五月婷婷| 一进一出抽搐动态| 麻豆av噜噜一区二区三区| 少妇的逼水好多| 国产高潮美女av| www.色视频.com| 久久午夜福利片| 午夜精品国产一区二区电影 | 26uuu在线亚洲综合色| 国产蜜桃级精品一区二区三区| 日韩欧美三级三区| 在线免费观看的www视频| 最好的美女福利视频网| 夜夜看夜夜爽夜夜摸| 亚洲成人久久性| 精品一区二区三区人妻视频| 亚洲中文字幕日韩| 变态另类成人亚洲欧美熟女| 又粗又爽又猛毛片免费看| 日韩精品有码人妻一区| 日日摸夜夜添夜夜爱| 少妇熟女欧美另类| .国产精品久久| 亚洲欧洲国产日韩| 亚洲精品影视一区二区三区av| 精品久久久久久久久久久久久| 岛国毛片在线播放| 美女脱内裤让男人舔精品视频 | 最近手机中文字幕大全| 99久久无色码亚洲精品果冻| 99久国产av精品国产电影| 欧美一级a爱片免费观看看| 一级黄色大片毛片| 村上凉子中文字幕在线| 中文亚洲av片在线观看爽| 久久精品人妻少妇| 深夜精品福利| 国产成人aa在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产在线男女| 成人特级黄色片久久久久久久| 国产成人精品婷婷| 久久九九热精品免费| 国产亚洲av片在线观看秒播厂 | 久久久久久久亚洲中文字幕| 成人亚洲欧美一区二区av| 国产精品福利在线免费观看| 免费无遮挡裸体视频| 国产白丝娇喘喷水9色精品| av福利片在线观看| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 五月伊人婷婷丁香| 男女视频在线观看网站免费| 国产单亲对白刺激| 一本精品99久久精品77| 免费av不卡在线播放| 只有这里有精品99| 成年女人永久免费观看视频| 禁无遮挡网站| 乱码一卡2卡4卡精品| 国产熟女欧美一区二区| 久久中文看片网| 51国产日韩欧美| 美女脱内裤让男人舔精品视频 | 青青草视频在线视频观看| 久久久久久久久中文| 男女下面进入的视频免费午夜| 老司机影院成人| 亚洲国产高清在线一区二区三| 免费大片18禁| 老师上课跳d突然被开到最大视频| 天堂网av新在线| 国产高清有码在线观看视频| 观看美女的网站| 非洲黑人性xxxx精品又粗又长| 中文字幕久久专区| 欧美一区二区精品小视频在线| 亚洲最大成人av| 婷婷色av中文字幕| 特级一级黄色大片| 亚洲国产高清在线一区二区三| av又黄又爽大尺度在线免费看 | 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 91在线精品国自产拍蜜月| av专区在线播放| 亚洲美女视频黄频| 欧美3d第一页| 日韩,欧美,国产一区二区三区 | 国产伦精品一区二区三区视频9| 深爱激情五月婷婷| 亚洲欧美中文字幕日韩二区| 精品久久久久久久久久免费视频| 人妻夜夜爽99麻豆av| 成人午夜高清在线视频| 麻豆av噜噜一区二区三区| 欧美+日韩+精品| 亚洲精品乱码久久久久久按摩| 亚洲国产精品成人综合色| 老师上课跳d突然被开到最大视频| 欧美日韩在线观看h| 麻豆av噜噜一区二区三区| 人妻夜夜爽99麻豆av| 亚洲国产精品成人综合色| 熟妇人妻久久中文字幕3abv| 亚洲欧美日韩高清在线视频| 国产精品,欧美在线| 亚洲一区高清亚洲精品| 国产真实伦视频高清在线观看| 一进一出抽搐动态| 日韩av不卡免费在线播放| 男人和女人高潮做爰伦理| 网址你懂的国产日韩在线| 日韩欧美一区二区三区在线观看| 久久久a久久爽久久v久久| 成人国产麻豆网| 99热这里只有是精品在线观看| 网址你懂的国产日韩在线| 搡老妇女老女人老熟妇| 白带黄色成豆腐渣| 97热精品久久久久久| 亚洲中文字幕日韩| 麻豆成人午夜福利视频| 蜜臀久久99精品久久宅男| 日本撒尿小便嘘嘘汇集6| 国产精品福利在线免费观看| 人妻系列 视频| 悠悠久久av| 日韩高清综合在线| 国产精品一区二区在线观看99 | 亚洲人成网站在线观看播放| 亚洲精品久久久久久婷婷小说 | 婷婷精品国产亚洲av| av国产免费在线观看| 卡戴珊不雅视频在线播放| 别揉我奶头 嗯啊视频| 老熟妇乱子伦视频在线观看| 亚洲精品亚洲一区二区| 免费观看a级毛片全部| 日本爱情动作片www.在线观看| 一级av片app| 国产精品一区二区在线观看99 | 热99在线观看视频| 麻豆成人午夜福利视频| 国产午夜福利久久久久久| 国产淫片久久久久久久久| 亚洲人成网站在线播放欧美日韩| 女人十人毛片免费观看3o分钟| 午夜免费男女啪啪视频观看| 美女xxoo啪啪120秒动态图| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 尾随美女入室| 听说在线观看完整版免费高清| 麻豆精品久久久久久蜜桃| 1000部很黄的大片| 狂野欧美白嫩少妇大欣赏| 身体一侧抽搐| 性色avwww在线观看| 99久久精品一区二区三区| 日韩三级伦理在线观看| 国产乱人偷精品视频| 国产高清激情床上av| 国产色爽女视频免费观看| 永久网站在线| 99久久精品一区二区三区| 国产精品久久久久久亚洲av鲁大| а√天堂www在线а√下载| 别揉我奶头 嗯啊视频| 日本色播在线视频| 男插女下体视频免费在线播放| 国产欧美日韩精品一区二区| 国产私拍福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 寂寞人妻少妇视频99o| 日日摸夜夜添夜夜爱| 看黄色毛片网站| 有码 亚洲区| 成人一区二区视频在线观看| 欧美一区二区亚洲| 日韩国内少妇激情av| 亚洲最大成人av| 国产午夜精品论理片| or卡值多少钱| 白带黄色成豆腐渣| 国产一区二区三区在线臀色熟女| 成人永久免费在线观看视频| 国产精品女同一区二区软件| 免费av不卡在线播放| 免费搜索国产男女视频| 国产一级毛片七仙女欲春2| 国产亚洲av嫩草精品影院| 不卡一级毛片| av在线播放精品| 成年av动漫网址| 久久久久久久久久久丰满| 日本黄色视频三级网站网址| 悠悠久久av| 亚洲av免费高清在线观看| .国产精品久久| а√天堂www在线а√下载| 亚洲欧美日韩东京热| 亚洲内射少妇av| 国产色婷婷99| 爱豆传媒免费全集在线观看| 嘟嘟电影网在线观看| 亚洲精品乱码久久久久久按摩| 69人妻影院| 麻豆成人av视频| 搞女人的毛片| 国产亚洲5aaaaa淫片| 久久久国产成人免费| 精品熟女少妇av免费看| 国产 一区精品| 亚洲综合色惰| 少妇熟女欧美另类| www.av在线官网国产| 99久久九九国产精品国产免费| 精品人妻视频免费看| 99久国产av精品国产电影| 国产老妇女一区| 久久久午夜欧美精品| 自拍偷自拍亚洲精品老妇| 91久久精品国产一区二区成人| 成人综合一区亚洲| 美女大奶头视频| 免费看美女性在线毛片视频| 人人妻人人看人人澡| 三级男女做爰猛烈吃奶摸视频| 欧美bdsm另类| 亚洲精品日韩av片在线观看| 亚洲av一区综合| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲| 爱豆传媒免费全集在线观看| 嫩草影院入口| 亚洲av成人精品一区久久| 99久久中文字幕三级久久日本| 成年免费大片在线观看| 日本黄大片高清| 久久韩国三级中文字幕| 欧美潮喷喷水| 成人鲁丝片一二三区免费| 欧美不卡视频在线免费观看| 亚洲av免费在线观看| 简卡轻食公司| 美女被艹到高潮喷水动态| 精品久久久久久久久av| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av天美| 极品教师在线视频| 亚洲av中文av极速乱| 毛片女人毛片| 国产精品国产高清国产av| 97热精品久久久久久| 国产黄片美女视频| 欧美日韩乱码在线| 麻豆成人av视频| 久久久精品欧美日韩精品| 欧美日本视频| 国产精品免费一区二区三区在线| 国产一级毛片在线| 波多野结衣高清无吗| 一级毛片电影观看 | 亚洲天堂国产精品一区在线| av天堂中文字幕网| 亚洲熟妇中文字幕五十中出| 国产在视频线在精品| 日日啪夜夜撸| 我要看日韩黄色一级片| 可以在线观看毛片的网站| 日本与韩国留学比较| 搡老妇女老女人老熟妇| 99热精品在线国产| 日韩av不卡免费在线播放| 国产一区二区三区av在线 | 国产亚洲精品久久久com| 亚洲国产精品久久男人天堂| 日韩国内少妇激情av| ponron亚洲| 午夜福利在线观看吧| 天堂影院成人在线观看| 不卡一级毛片| 欧美一区二区国产精品久久精品| 成年av动漫网址| 中文资源天堂在线| 国产精品精品国产色婷婷| 国产伦一二天堂av在线观看| 69人妻影院| h日本视频在线播放| 成人无遮挡网站| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器| 麻豆久久精品国产亚洲av| 美女 人体艺术 gogo| 欧美日韩精品成人综合77777| 亚洲欧美精品综合久久99| 99国产极品粉嫩在线观看| 国产精品久久久久久久久免| 亚洲欧美精品自产自拍| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 国产高清有码在线观看视频| 亚洲,欧美,日韩| 男女边吃奶边做爰视频| 中文资源天堂在线| 国产一区二区激情短视频| 直男gayav资源| kizo精华| 国产精品三级大全| 国产精品久久电影中文字幕| 午夜a级毛片| 久久精品久久久久久久性| 亚洲真实伦在线观看| 国产精品日韩av在线免费观看| 亚洲久久久久久中文字幕| 日韩欧美精品v在线| 国产精品美女特级片免费视频播放器| 国产日本99.免费观看| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 美女被艹到高潮喷水动态| а√天堂www在线а√下载| 国产熟女欧美一区二区| 久久午夜亚洲精品久久| 久久精品人妻少妇| 久久久国产成人精品二区| 国产精品,欧美在线| 久久这里有精品视频免费| 亚洲av第一区精品v没综合| 在线观看66精品国产| 可以在线观看的亚洲视频| 精品人妻视频免费看| 欧美另类亚洲清纯唯美| 亚洲18禁久久av| 青春草视频在线免费观看| 最近视频中文字幕2019在线8| 少妇熟女aⅴ在线视频| 日韩欧美精品免费久久|