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

    2009至2016年中國甲型H1N1流感病毒神經(jīng)氨酸酶基因進化分析及流行趨勢預測

    2017-08-17 10:04:02宋玥顏文娟方坤李小杉李偉衛(wèi)平民
    臨床檢驗雜志 2017年7期
    關(guān)鍵詞:分析

    宋玥,顏文娟,方坤,李小杉,李偉,衛(wèi)平民

    (1.東南大學公共衛(wèi)生學院,南京 210009;2.南京市疾病預防控制中心急性傳染病防制與學校衛(wèi)生科,南京 210009)

    ·調(diào)查研究·

    2009至2016年中國甲型H1N1流感病毒神經(jīng)氨酸酶基因進化分析及流行趨勢預測

    宋玥1,顏文娟1,方坤1,李小杉1,李偉2,衛(wèi)平民1

    (1.東南大學公共衛(wèi)生學院,南京 210009;2.南京市疾病預防控制中心急性傳染病防制與學校衛(wèi)生科,南京 210009)

    目的 分析2009至2016年中國甲型H1N1流感病毒神經(jīng)氨酸酶基因(NA)進化特征,預測其流行趨勢,為流感疫苗評價提供依據(jù)。方法 從全球共享禽流感數(shù)據(jù)倡議組織(GISAID)數(shù)據(jù)庫及美國國家生物技術(shù)中心(NCBI)數(shù)據(jù)庫下載甲型H1N1流感病毒NA基因序列1 141條。利用鄰接法進行系統(tǒng)進化分析和遺傳變異分析,利用Bayesian skyline Plot預測流行趨勢。結(jié)果 2009至2016年中國甲型H1N1毒株NA基因與參考毒株的同源性逐年降低。遺傳進化分析顯示同一年份的毒株在系統(tǒng)進化樹上基本呈現(xiàn)集中分布。除2012年,其余各年份的毒株均通過不同模型得到正向壓力選擇位點。動力學分析顯示2017年甲型H1N1可能達到一個流行高峰。結(jié)論 甲型H1N1流感病毒NA基因所編碼的氨基酸逐年變異,需進一步加強流感監(jiān)測。

    甲型H1N1流感;神經(jīng)氨酸酶;分子特征;進化分析

    甲型H1N1流感病毒[influenza A viruses(H1N1 09pdm)]是由禽流感、豬流感和人流感病毒3者基因重配而成[1-2],自感染暴發(fā)以來,H1N1(09pdm)在人群中一直引起較高的發(fā)病率,也因此受到國內(nèi)外研究的廣泛關(guān)注。

    流感病毒基因組由8個片段構(gòu)成,其中神經(jīng)氨酸酶(neuraminidase,NA)是流感病毒在傳播和擴散中最關(guān)鍵的酶。NA附著在流感病毒顆粒表面[3],具有水解唾液酸的活性,從而促使病毒的釋放和擴散[4],且NA抑制劑在合成抗流感病毒藥物中占主要地位。NA的抗原位點、糖基化位點和耐藥位點等關(guān)鍵的氨基酸突變是病毒抗原性發(fā)生變異的重要基礎(chǔ)[5-6],其基因的改變對流感疫情的發(fā)生與流行有著密切的影響。本研究通過分析2009年后甲型H1N1流感病毒在中國地區(qū)的流行態(tài)勢,對NA核苷酸序列進行系統(tǒng)進化及氨基酸位點變異分析,為我國甲型H1N1流感的疾病防控及預測提供理論依據(jù)。

    1 材料與方法

    1.1 材料 收集2009年至2016年國家流感中心(http://www.cnic.org.cn/chn/)流感周報中國流感流行的數(shù)據(jù),從美國國家生物技術(shù)中心流感病毒(National Center for Biotechnology Information,NCBI)數(shù)據(jù)庫(http://www.ncbi.nlm.nih.gov/genomes/FLU/FLU.html)及全球共享禽流感數(shù)據(jù)倡議組織(The Global Initiative on Sharing All Influenza Data,GISAID)數(shù)據(jù)庫(http://platform.gisaid.org/epi3/frontend)中下載中國地區(qū)甲型H1N1流感病毒NA基因序列,共1 141條。用Molecular Evolutionary Genetics Analysis (MEGA) version 6.0軟件進行序列比對,并用鄰接法繪制進化樹,根據(jù)樹中終末分支同源性較高的分支情況,刪除部分同地區(qū)同年份的序列。刪除重復的序列和部分同年份同地區(qū)同源性較高的序列,得到354條來自安徽、北京、重慶、福建、甘肅、廣東、廣西、貴州、海南、河北、黑龍江、河南、香港、湖北、湖南、江蘇、吉林、遼寧、內(nèi)蒙古、寧夏、青海、山西、陜西、山東、上海、四川、天津、西安、新疆、云南、浙江各省份的序列,包括2009年87條、2010年55條、2011年26條、2012年4條、2013年34條、2014年24條、2015年54條和2016年70條。

    1.2 方法 采用MEGA version 6.0軟件進行序列比對、氨基酸位點的變異分析,采用鄰接法構(gòu)建甲型H1N1病毒NA基因進化樹(Bootstrap檢驗重復值設(shè)為1 000)。通過Data monkey在線軟件的Single likelihood ancestor counting(SLAC)、Internal fixed effects likelihood(IFEL)和Random effects method(REL)模型進行選擇性壓力分析。利用Figtree v1.4.2(http://beast.bio.ed.ac.uk)讀取系統(tǒng)進化樹,以bootstrap值≥0.9判定為流行枝(Lineage)。采用BEAST v1.8.0軟件的Bayesian skyline Plot(BSP)法構(gòu)建Bayesian skyline reconstruction,進行動力學分析,模型為GTR+Relaxed clock(uncorrelated)+Bayesian skyline。MCMC分析鏈長為200 000 000,迭代頻次為1 000,收斂評價標準為有效樣本數(shù)(effective sampling size, ESS)>200。

    2 結(jié)果

    2.1 流行病學分布特征 從全國流感監(jiān)測網(wǎng)絡實驗室監(jiān)測數(shù)據(jù)顯示,2009年至2016年全國甲型H1N1流感共發(fā)生4次流行高峰,分別是2009年5月至2010年2月、2011年1至3月、2012年12月至2014年3月、2015年末至2016年又呈現(xiàn)上升趨勢。見圖1。

    圖1 中國地區(qū)2009至2016年流感流行情況

    2.2 同源性分析 將不同年份甲型H1N1流感毒株(354條)與參考毒株A/California/07/2009(H1N1)的NA基因序列進行比對,其中2009至2016年的NA基因與參考毒株的同源性分別為99.2%~99.7%、99.2%~99.6%、98.8%~99.5%、98.8%~98.9%、98.4%~99.6%、98.2%~98.9%、97.9%~99.6%和97.5%~98.9%。2009至2016年甲型H1N1流感毒株與參考毒株的NA基因核苷酸和氨基酸遺傳距離及各年內(nèi)核苷酸的平均遺傳距離見表1。

    2.3 系統(tǒng)進化分析 甲型H1N1亞型NA基因系統(tǒng)發(fā)育進化樹,如圖2所示(圖中的進化樹只包含部分毒株,剔除了部分位于進化樹上相同終末分支并且相同年份、同源性較高的序列,其中標為紅色分支的是參考序列),最終選取了有代表性的71支毒株建樹,不同年份的毒株由不同的顏色表示。結(jié)果顯示,同一年份的甲型H1N1流感毒株在系統(tǒng)進化樹上基本呈集中分布的態(tài)勢。其中2009年與2010年的流感毒株形成一個分支;2011年的流感毒株另外形成2個小分支,可能是由于V241I和N369K的變異;2012年毒株N200S突變,獨立形成一個小分支;2013年至2016年的毒株聚集形成一個大支。毒株間存在明顯的年代差異,隨著年代相隔越久,序列差異越大。

    注:標尺值0.002代表毒株間的進化距離。

    圖2 中國2009-2016年甲型H1N1流感毒株NA基因系統(tǒng)進化樹

    2.4 NA分子特征分析

    2.4.1 潛在抗原位點的氨基酸序列分析 NA蛋白N端第140至157位氨基酸很有可能是甲型H1N1流感病毒抗原決定簇所在[4]。中國地區(qū)的354條毒株序列與疫苗株A/California/07/2009(H1N1)相比,2010年潛在抗原位點氨基酸發(fā)生變異(G147E);2015年潛在抗原位點發(fā)生變異(I145N),并且第151位氨基酸出現(xiàn)缺失情況;2016年潛在抗原位點第151位氨基酸也出現(xiàn)缺失情況,潛在抗原位點發(fā)生變異(Y155H)。

    2.4.2 酶活性中心及周圍相關(guān)位點的氨基酸序列分析 NA蛋白的酶催化活性位點包括催化部位的第118、151、152、224、276、292、371和406位等8個位點,除此還有酶活性中心周圍的輔助位點11個,分別是第119、156、178、179、198、222、227、274、277、294和425位[4]。通過軟件比對,發(fā)現(xiàn)2009年至2016年354條甲型H1N1流感病毒的NA基因所編碼氨基酸在所有A型流感病毒的NA蛋白中都是保守的,未發(fā)現(xiàn)氨基酸位點突變。

    2.4.3 糖基化位點的氨基酸序列分析 甲型H1N1流感病毒NA蛋白分子糖基化位點有8個,分別是第50、58、63、68、88、146、235和386位[3-4]。通過軟件比對,發(fā)現(xiàn)2009至2016年354條甲型H1N1流感病毒的NA基因所編碼的氨基酸在第50、58、63、68、88、146、235位未發(fā)生變異,保持穩(wěn)定,但在第386位逐年變異,直至2016年該位點消失。

    2.4.4 耐藥性分析 NA某些氨基酸殘基的突變會使病毒產(chǎn)生對這類抗流感病毒藥物的耐藥性,包括E119V、R292K、 R293K、N294S和H274Y位點[7-8],2009至2016年354條毒株序列尚未發(fā)現(xiàn)上述位點的突變, 表明對NA抑制劑仍敏感。

    2.5 選擇性壓力分析 對2009至2016年各年份的NA序列分別采用SLAC、IFEL和REL模型進行選擇性壓力分析。各年份通過SLAC模型均未發(fā)現(xiàn)正向選擇壓力位點。IFEL模型發(fā)現(xiàn)2009年的序列有正向選擇壓力位點64;2010年有正向選擇壓力位點82和264;2013年有正向選擇壓力位點34、275和321;2016年有正向選擇壓力位點275,其余年份未發(fā)現(xiàn)。除2012年的序列外,通過REL模型均得到正向選擇壓力位點,結(jié)果見表2。

    表2 中國2009至2016年甲型H1N1流感毒株正向選擇壓力位點分布

    年份SLACIFELREL2009N/D6464、3912010N/D82、26462、82、3862011N/DN/D145、1792012N/DN/DN/D2013N/D34、275、32134、44、286、3092014N/DN/D144、2892015N/DN/D34、40、314、3462016N/D27534、40、93、96、257、275、314、396

    注:N/D表示未發(fā)現(xiàn)。

    2.6 進化動力學分析 Bayesian skyline reconstruction建立的模型表明,流感病毒遺傳多樣性在2009年維持穩(wěn)定狀態(tài),并且在2009年和2010年達到高峰,隨后分別在2011年及2012年有2次下降趨勢,2012年6月呈現(xiàn)上升趨勢,于2013年和2014年再次達到流行高峰,此后至2015年下降并維持于相對穩(wěn)定狀態(tài),2016年開始又呈現(xiàn)上升趨勢,預測2017年可能又將達到一個流行高峰。見圖3。

    注:黑色實線為估計中位數(shù),藍色區(qū)域為95%置信區(qū)間。

    圖3 甲型H1N1流感NA基因進化動力學

    3 討論

    根據(jù)國家流感中心的監(jiān)測結(jié)果,2009至2016年中國流感流行主要以甲型H1N1、B型和季節(jié)性H3N2為主,其中甲型H1N1亞型作為我國流感流行的重要亞型,共發(fā)生4次流行高峰,主要在冬春季。2009至2016年甲型H1N1流感病毒NA基因與參考毒株A/California/07/2009(H1N1)的同源性逐年降低,類似研究也提示甲型H1N1流感病毒逐漸進化變異[9]。本研究發(fā)現(xiàn)NA有16個氨基酸位點改變,分別為V13I、I34V、N44S、V67I、N200S、V241I、N248D、V264I、N270K、I314M、I321V、I329N、N369K、T381I、N386K、K432E,其他變異位點呈散發(fā)型。核苷酸與氨基酸的遺傳距離隨著年份的增加而不斷增大,遺傳進化速度加快。流感病毒NA所編碼的氨基酸關(guān)鍵位點突變增多??乖稽c的突變、糖基化位點的消失、酶活性位點的變異等都可能會影響流感病毒的抗原性和致病力[10],從而使病毒適應宿主感染環(huán)境,在人群中更廣泛地傳播。NA蛋白的酶活性中心位于NA的頭部[11],可以直接或間接參與酶的催化作用。2009至2016年甲型H1N1流感病毒NA蛋白分子酶活性中心及周圍相關(guān)位點相對保守,有利于穩(wěn)定酶的活性及病毒擴散能力。NA催化活性位點可以直接或間接參與NA的催化功能,如有氨基酸變異將會影響NA的活性[12],NA抑制劑設(shè)計的原則就是基于NA高度保守的酶活性位點[13-14]。本研究并未發(fā)現(xiàn)酶活性位點或耐藥位點的變異,表明相關(guān)位點仍可作為抗流感藥物研發(fā)的基因靶點[11,15]。糖基化位點是維持NA蛋白分子結(jié)構(gòu)和正常生物功能的重要影響因素,糖基化位點的變異對穩(wěn)定病毒蛋白結(jié)構(gòu)和生物特性有很大影響[16]。2009年新型甲型H1N1流感病毒具有8個糖基化位點,其中N386K隨著年份的遞增,直至徹底消失。研究表明,糖基化位點的非糖基化以及潛在糖基化位點數(shù)量的變化都可能會影響蛋白質(zhì)折疊和運輸?shù)忍匦院凸δ?,進而導致流感病毒結(jié)合受體的能力和病毒致病性降低[17]。病毒在進化過程中會獲得新的糖基化位點,新的糖基化位點可能會掩蓋抗原位點,成為病毒免疫逃逸的一種方式,需要加強監(jiān)測。正向選擇位點與病毒受體結(jié)合能力有關(guān),因此這些正向位點的選擇可能有助于H1N1流感的宿主從豬到人的適應。正向選擇位點常發(fā)生于免疫水平較低的人群,會影響疫苗的近期效果評價[18]。

    Ebranati等[19]認為病毒遺傳多樣性的波動與其所引起的暴發(fā)流行相對應。進化動力學趨勢顯示,甲型H1N1流感的流行趨勢呈階段性,每流行2年會間隔一年,與類似研究結(jié)果一致[20]。其中趨勢圖下降可能由于流感暴發(fā)后,國家采取了及時的公共衛(wèi)生措施,從而控制了流感疫情,趨勢圖的上升可能由于病毒在自然界的不斷變異,從而逃避宿主免疫引起流感的流行。

    綜上所述,隨著流感病毒在人群中的進一步適應和持續(xù)存在,病毒分子特征將發(fā)生變異,預示該亞型在未來的時間里很有可能繼續(xù)成為優(yōu)勢毒株。本研究追蹤甲型H1N1流感病毒的進化趨勢,特別關(guān)注基因變化對病毒的傳播力和致病性的影響,從而及時發(fā)現(xiàn)具有流行病學意義的變異株和耐藥株,這對疾病防控及預測具有重要的意義。

    [1]韓一芳, 謝佳新, 殷建華, 等. 2009年新型甲型H1N1流感病毒基質(zhì)蛋白及核蛋白基因進化分析[J]. 第二軍醫(yī)大學學報, 2009, 30(6): 622-627.

    [2]方斌, 劉琳琳, 李翔, 等. 湖北省大流行期和流行后期新甲型H1N1流感病毒基因特性分析 [J]. 中華疾病控制雜志, 2017, 21(1): 13-18.

    [3]趙紅玲, 章文, 黃濤. 甲型H1N1流感病毒基因特征及耐藥性研究[J]. 中國病原生物學雜志, 2014,9(9): 811-815.

    [4]蘇彤, 李淑華, 常文軍, 等. 2009年新型甲型H1N1流感病毒神經(jīng)氨酸酶基因進化分析[J]. 第二軍醫(yī)大學學報, 2009, 30(6): 618-621.

    [5]Koel B F, Mogling R, Chutinimitkul S,etal. Identification of amino acid substitutions supporting antigenic change of influenza A(H1N1)pdm09 viruses [J]. J Virol, 2015, 89(7): 3763-3775.

    [6]王楊,崔大偉,黃明珠,等.新型甲型H7N9禽流感病毒HA與NA的進化分析[J].臨床檢驗雜志,2013,31(7):499-502.

    [7]盧桂蘭, 彭曉旻, 劉醫(yī)萌, 等. 北京市2010年-2011年兒童甲型H3N2流感病毒神經(jīng)氨酸酶基因進化和耐藥性位點分析[J]. 中國衛(wèi)生檢驗雜志, 2012,22(6): 1232-1236.

    [8]祁賢, 湯奮揚, 李亮, 等. 新甲型H1N1(2009)流感病毒的早期分子特征[J]. 微生物學報, 2010,50(1): 81-90.

    [9]顏文娟, 郭艷, 李小杉, 等. 2009-2015年中國甲型H1N1流感病毒血凝素基因進化分析[J]. 微生物學通報, 2017, 44(2): 420-427.

    [10]Khandaker I, Suzuki A, Kamigaki T,etal. Molecular evolution of the hemagglutinin and neuraminidase genes of pandemic (H1N1) 2009 influenza viruses in Sendai, Japan, during 2009-2011 [J]. Virus Genes, 2013,47(3): 456-466.

    [11]張建國, 張繼寧, 黃勛娟. 神經(jīng)氨酸酶抑制劑的研究進展[J]. 化學與生物工程, 2014, 31(4): 1-5,23.

    [12]吳英萍,崔大偉,秦志梅,等. 2009~2012年杭州地區(qū)甲型H1N1流感病毒HA與NA基因的進化分析[J]. 臨床檢驗雜志,2013,31(4):272-275.

    [13]Richard M, Deleage C, Barthelemy M,etal. Impact of influenza A virus neuraminidase mutations on the stability, activity, and sensibility of the neuraminidase to neuraminidase inhibitors [J]. J Clin Virol, 2008, 41(1): 20-24.

    [14]陳寅,茅海燕,周敏,等. 浙江省近年甲型H1N1流感病毒神經(jīng)氨酸酶基因序列分析[J].疾病監(jiān)測,2013,28(10):811-814.

    [15]楊先知,崔大偉,謝國良,等.杭州市2013-2014年流感暴發(fā)流行期甲型流感病毒H3N2亞型基因特性分析[J].臨床檢驗雜志,2015,33(8):586-590.

    [16]謝劍鋒, 沈曉娜, 王美愛, 等. 福建省甲型H1N1流感病毒基因特征研究[J]. 病毒學報, 2014, 30(1): 37-43.

    [17]李向忠,方芳,陳則. 流感病毒神經(jīng)氨酸酶不同區(qū)域的作用[J]. 生命科學研究, 2005, 9(2): 59-65.

    [18]Smith DJ, Lapedes AS, De Jong JC,etal. Mapping the antigenic and genetic evolution of influenza virus[J]. Science, 2004, 305(5682): 371-376.

    [19]Ebranati E, Pariani E, Piralla A,etal. Reconstruction of the evolutionary dynamics of A(H3N2) influenza viruses circulating in Italy from 2004 to 2012[J]. PLoS One, 2015, 10(9): e0137099.

    [20]劉丹, 王建紅, 陳洪永, 等. 2009年-2015年唐山地區(qū)流感病毒病原學特征[J]. 中國衛(wèi)生檢驗雜志, 2017, 27(2): 185-187.

    (本文編輯:劉群)

    Evolution analysis and trend prediction of neuraminidase gene of influenza A/H1N1(09pdm) virus in China from 2009 to 2016

    SONGYue1,YANWen-juan1,FANGKun1,LIXiao-shan1,LIWei2,WEIPing-min1

    (1.SchoolofPublicHealth,South-EastUniversity,Nanjing210009,Jiangsu; 2.InfectiousDiseaseControlandSchoolHealth,NanjingCenterforDiseaseControlandPrevention,Nanjing210009,Jiangsu,China)

    Objective To analyze the evolutionary characteristics and predict the variation trend of neuraminidase (NA) gene of influenza A/H1N1 (09pdm) virus in China from 2009 to 2016 in order to provide the basis of assessment for flu vaccines. Methods A total of 1 141 sequences of NA gene of influenza A/H1N1 (09pdm) virus were screened out from the Global Initiative on Sharing All Influenza Data and National Center for Biotechnology Information. The phylogenetic trees and the mutations of amino acids sequences were constructed and analyzed by biological softwares. The prediction for epidemic trend of influenza was analyzed by Bayesian skyline Plot. Results Compared with the sequence of reference strain, the homology of nucleic acid sequence of NA gene decreased year by year from 2009 to 2016. The phylogenetic analysis showed that NA gene clustered nearly on the identical phylogenetic tree in one year. The positive selection pressure site of NA strain was observed by different models in each year except 2012. The dynamics analysis showed that the popularity of influenza A/H1N1 virus may continue to increase to a peak in 2017. Conclusion The amino acid encoded by NA gene of influenza A/H1N1 virus is varying gradually, so the importance of surveillance for influenza virus should be reinforced for every year.

    influenza A/H1N1(09pdm); neuraminidase; molecular characteristics; evolution analysis

    10.13602/j.cnki.jcls.2017.07.19

    江蘇省重點研發(fā)計劃(社會發(fā)展)項目(BE2015679);江蘇省“六大人才高峰”第十一批高層次人才選拔培養(yǎng)資助項目(2014-WSN-40)。

    宋玥,1992年生,女,碩士研究生,研究方向:流行病與衛(wèi)生統(tǒng)計學。

    衛(wèi)平民,教授,博士,E-mail: mpw1963@126.com。

    R373.1+3

    A

    2017-03-29)

    猜你喜歡
    分析
    禽大腸桿菌病的分析、診斷和防治
    隱蔽失效適航要求符合性驗證分析
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    經(jīng)濟危機下的均衡與非均衡分析
    對計劃生育必要性以及其貫徹實施的分析
    GB/T 7714-2015 與GB/T 7714-2005對比分析
    出版與印刷(2016年3期)2016-02-02 01:20:11
    網(wǎng)購中不良現(xiàn)象分析與應對
    中西醫(yī)結(jié)合治療抑郁癥100例分析
    偽造有價證券罪立法比較分析
    丝袜美腿诱惑在线| 99国产精品99久久久久| 亚洲成色77777| 少妇粗大呻吟视频| 精品一区二区三区av网在线观看 | 一级毛片女人18水好多 | 性色av乱码一区二区三区2| 精品人妻熟女毛片av久久网站| 激情视频va一区二区三区| 看免费av毛片| 久久九九热精品免费| 少妇精品久久久久久久| 又紧又爽又黄一区二区| xxxhd国产人妻xxx| 最新在线观看一区二区三区 | 精品亚洲乱码少妇综合久久| 高清欧美精品videossex| 后天国语完整版免费观看| 视频区欧美日本亚洲| 日本五十路高清| 下体分泌物呈黄色| 日本一区二区免费在线视频| 欧美日韩亚洲综合一区二区三区_| 男女无遮挡免费网站观看| 亚洲精品第二区| 国产免费视频播放在线视频| 亚洲av成人精品一二三区| 成人午夜精彩视频在线观看| a 毛片基地| 久久影院123| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品第一综合不卡| 黑人欧美特级aaaaaa片| 国产黄频视频在线观看| 久久热在线av| 国产高清不卡午夜福利| 黄色视频在线播放观看不卡| 日韩大码丰满熟妇| 两个人免费观看高清视频| 纯流量卡能插随身wifi吗| 亚洲中文日韩欧美视频| 大片电影免费在线观看免费| av有码第一页| 亚洲五月色婷婷综合| 黄片播放在线免费| 亚洲美女黄色视频免费看| 国产精品.久久久| 亚洲少妇的诱惑av| 精品久久久精品久久久| 啦啦啦在线观看免费高清www| 国产成人精品久久二区二区91| 最近最新中文字幕大全免费视频 | tube8黄色片| 久久久久网色| 男女高潮啪啪啪动态图| 在线观看免费日韩欧美大片| www.999成人在线观看| 国产日韩欧美视频二区| 亚洲国产中文字幕在线视频| 老鸭窝网址在线观看| 尾随美女入室| 妹子高潮喷水视频| 国产99久久九九免费精品| 亚洲美女黄色视频免费看| av有码第一页| 亚洲国产毛片av蜜桃av| 国产99久久九九免费精品| 亚洲欧美成人综合另类久久久| 精品人妻一区二区三区麻豆| 国产日韩一区二区三区精品不卡| 香蕉国产在线看| 精品久久久久久久毛片微露脸 | 视频区图区小说| 美女脱内裤让男人舔精品视频| 成年动漫av网址| 欧美精品啪啪一区二区三区 | 午夜久久久在线观看| 久久精品亚洲av国产电影网| 女性生殖器流出的白浆| 国产欧美日韩一区二区三区在线| 2021少妇久久久久久久久久久| 国产精品秋霞免费鲁丝片| 精品国产乱码久久久久久男人| 亚洲国产毛片av蜜桃av| 黄色视频在线播放观看不卡| 精品少妇久久久久久888优播| 男女下面插进去视频免费观看| 欧美黑人欧美精品刺激| 亚洲成人免费av在线播放| 国产视频一区二区在线看| 91九色精品人成在线观看| 欧美大码av| 极品人妻少妇av视频| 亚洲视频免费观看视频| 久久人人97超碰香蕉20202| 无限看片的www在线观看| 中文字幕最新亚洲高清| 女警被强在线播放| 亚洲av欧美aⅴ国产| 最近中文字幕2019免费版| 精品视频人人做人人爽| 国产精品久久久久久精品古装| av电影中文网址| 亚洲,一卡二卡三卡| 两人在一起打扑克的视频| 国产成人精品无人区| av不卡在线播放| av片东京热男人的天堂| 大香蕉久久成人网| 亚洲av日韩在线播放| 两个人免费观看高清视频| 悠悠久久av| 在线 av 中文字幕| 精品少妇久久久久久888优播| 一边亲一边摸免费视频| 国产精品秋霞免费鲁丝片| 色精品久久人妻99蜜桃| 亚洲精品美女久久av网站| h视频一区二区三区| 免费日韩欧美在线观看| 99久久99久久久精品蜜桃| 又大又爽又粗| 免费观看人在逋| 国产一区二区三区综合在线观看| 国产免费又黄又爽又色| 国产成人欧美| 一级黄片播放器| 最近手机中文字幕大全| 亚洲欧美一区二区三区久久| 日本vs欧美在线观看视频| av福利片在线| 大香蕉久久网| 国产精品国产av在线观看| 国产野战对白在线观看| 飞空精品影院首页| 久久青草综合色| 国产精品一国产av| 汤姆久久久久久久影院中文字幕| 美女大奶头黄色视频| 香蕉丝袜av| 女人久久www免费人成看片| 免费在线观看完整版高清| 看免费成人av毛片| 亚洲av电影在线进入| 久久精品国产亚洲av高清一级| 男女边摸边吃奶| 亚洲天堂av无毛| 亚洲伊人久久精品综合| 大码成人一级视频| 成人影院久久| 久久99热这里只频精品6学生| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区久久| 美女国产高潮福利片在线看| 日韩中文字幕欧美一区二区 | 午夜免费观看性视频| 尾随美女入室| 80岁老熟妇乱子伦牲交| 人人妻,人人澡人人爽秒播 | 夜夜骑夜夜射夜夜干| 国产成人啪精品午夜网站| 成人国产av品久久久| 男男h啪啪无遮挡| 女人爽到高潮嗷嗷叫在线视频| a级毛片黄视频| 老司机靠b影院| 狠狠精品人妻久久久久久综合| 91成人精品电影| 国产精品香港三级国产av潘金莲 | 亚洲精品国产一区二区精华液| 国产高清videossex| 亚洲精品久久久久久婷婷小说| 一级片免费观看大全| 国产欧美日韩精品亚洲av| av天堂在线播放| 久久久久久久国产电影| 亚洲国产精品国产精品| 热re99久久精品国产66热6| 精品一区二区三区av网在线观看 | 老鸭窝网址在线观看| 国产女主播在线喷水免费视频网站| 九草在线视频观看| 精品一区二区三卡| av天堂在线播放| 欧美 日韩 精品 国产| 亚洲av电影在线观看一区二区三区| 99国产综合亚洲精品| 国产精品.久久久| 日韩电影二区| 韩国精品一区二区三区| 中文字幕av电影在线播放| 可以免费在线观看a视频的电影网站| 国产极品粉嫩免费观看在线| 大型av网站在线播放| 久久精品国产亚洲av涩爱| 悠悠久久av| 脱女人内裤的视频| 久久青草综合色| 日本av手机在线免费观看| 一本综合久久免费| 一区二区av电影网| 日本av免费视频播放| 制服诱惑二区| 麻豆国产av国片精品| 欧美日韩福利视频一区二区| 青春草视频在线免费观看| 9热在线视频观看99| 真人做人爱边吃奶动态| 中文字幕高清在线视频| 精品卡一卡二卡四卡免费| 成人三级做爰电影| 国产成人啪精品午夜网站| 大片电影免费在线观看免费| 在线观看免费午夜福利视频| 国产一区二区三区综合在线观看| 亚洲中文av在线| 国产又爽黄色视频| 9191精品国产免费久久| 2021少妇久久久久久久久久久| 我要看黄色一级片免费的| 99九九在线精品视频| 国产成人a∨麻豆精品| 亚洲精品av麻豆狂野| 国产成人精品久久久久久| 色94色欧美一区二区| 国产主播在线观看一区二区 | 久久久久网色| 一本久久精品| 精品第一国产精品| 欧美大码av| xxx大片免费视频| 美女主播在线视频| 精品熟女少妇八av免费久了| 又大又黄又爽视频免费| 亚洲国产av影院在线观看| e午夜精品久久久久久久| 大型av网站在线播放| 叶爱在线成人免费视频播放| 纵有疾风起免费观看全集完整版| 亚洲欧美色中文字幕在线| avwww免费| 亚洲国产精品国产精品| 亚洲专区中文字幕在线| 国产欧美日韩一区二区三区在线| 一二三四在线观看免费中文在| 午夜福利影视在线免费观看| 高潮久久久久久久久久久不卡| 久久久久久亚洲精品国产蜜桃av| 久久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 如日韩欧美国产精品一区二区三区| 亚洲七黄色美女视频| 欧美在线一区亚洲| 欧美精品人与动牲交sv欧美| 精品人妻熟女毛片av久久网站| 久热爱精品视频在线9| 亚洲欧美日韩另类电影网站| 王馨瑶露胸无遮挡在线观看| 免费在线观看视频国产中文字幕亚洲 | av视频免费观看在线观看| 男女床上黄色一级片免费看| 国产三级黄色录像| 国产在线免费精品| 丝袜美腿诱惑在线| 电影成人av| 一级片免费观看大全| 最新的欧美精品一区二区| 69精品国产乱码久久久| 亚洲天堂av无毛| 亚洲伊人久久精品综合| 色综合欧美亚洲国产小说| 国产精品一区二区在线不卡| 人体艺术视频欧美日本| 中国美女看黄片| 国产高清不卡午夜福利| 熟女av电影| 免费看不卡的av| 91字幕亚洲| 搡老岳熟女国产| 免费在线观看完整版高清| 考比视频在线观看| 欧美97在线视频| 欧美性长视频在线观看| 国产成人系列免费观看| 久久 成人 亚洲| 一级毛片我不卡| 美女国产高潮福利片在线看| 久久人人爽人人片av| 亚洲精品国产区一区二| 18禁裸乳无遮挡动漫免费视频| 亚洲伊人久久精品综合| 一本色道久久久久久精品综合| xxx大片免费视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品国产av蜜桃| 日韩av在线免费看完整版不卡| 国产成人精品无人区| 美女主播在线视频| 久久久久久久久免费视频了| 国产片特级美女逼逼视频| 国产在线一区二区三区精| 午夜福利免费观看在线| 69精品国产乱码久久久| 我的亚洲天堂| 精品国产超薄肉色丝袜足j| 亚洲精品久久成人aⅴ小说| 性色av乱码一区二区三区2| 最新在线观看一区二区三区 | 人人妻人人添人人爽欧美一区卜| 久久国产精品影院| 精品一区二区三卡| 亚洲七黄色美女视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 国产成人精品久久久久久| 国产精品人妻久久久影院| 久久久久网色| 人妻 亚洲 视频| 交换朋友夫妻互换小说| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| av视频免费观看在线观看| 一级片免费观看大全| 黑人猛操日本美女一级片| 九色亚洲精品在线播放| 久久久精品免费免费高清| 免费看十八禁软件| 久久热在线av| 国产精品久久久久久精品古装| 日韩,欧美,国产一区二区三区| 美女主播在线视频| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 中文字幕另类日韩欧美亚洲嫩草| 日韩大片免费观看网站| 欧美 亚洲 国产 日韩一| 无遮挡黄片免费观看| 性色av乱码一区二区三区2| 午夜激情久久久久久久| 欧美日韩av久久| kizo精华| 波多野结衣一区麻豆| 久久国产亚洲av麻豆专区| 亚洲av成人不卡在线观看播放网 | 亚洲久久久国产精品| 一区二区三区精品91| 国产成人av激情在线播放| 午夜视频精品福利| 国产成人精品在线电影| 日韩大片免费观看网站| 中文乱码字字幕精品一区二区三区| 久热这里只有精品99| 日韩av不卡免费在线播放| 国产日韩一区二区三区精品不卡| 日韩中文字幕欧美一区二区 | 日本色播在线视频| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美中文字幕日韩二区| 亚洲国产精品国产精品| 午夜福利在线免费观看网站| 可以免费在线观看a视频的电影网站| 一本大道久久a久久精品| 亚洲欧美激情在线| 另类亚洲欧美激情| 国产色视频综合| 欧美精品一区二区大全| 国产视频首页在线观看| 亚洲欧美清纯卡通| 91成人精品电影| 老司机影院毛片| 精品久久久久久电影网| 天天躁日日躁夜夜躁夜夜| 国产高清视频在线播放一区 | 亚洲精品第二区| 久久精品aⅴ一区二区三区四区| 日本a在线网址| 午夜av观看不卡| 一本色道久久久久久精品综合| 真人做人爱边吃奶动态| 欧美成人精品欧美一级黄| videos熟女内射| 老司机影院毛片| 丝袜人妻中文字幕| 日本av手机在线免费观看| 一边摸一边做爽爽视频免费| 国产亚洲av片在线观看秒播厂| 国产精品二区激情视频| 欧美亚洲 丝袜 人妻 在线| 久久综合国产亚洲精品| 在现免费观看毛片| 久久精品国产a三级三级三级| 女人久久www免费人成看片| 免费人妻精品一区二区三区视频| 日韩av免费高清视频| 午夜福利影视在线免费观看| 免费少妇av软件| www.av在线官网国产| 国产成人一区二区三区免费视频网站 | 宅男免费午夜| 在线观看免费日韩欧美大片| www.熟女人妻精品国产| 国产日韩欧美视频二区| 18禁国产床啪视频网站| 一级毛片黄色毛片免费观看视频| 国产熟女午夜一区二区三区| 欧美日韩视频高清一区二区三区二| 精品少妇一区二区三区视频日本电影| 欧美精品高潮呻吟av久久| 少妇裸体淫交视频免费看高清 | 在线 av 中文字幕| 亚洲人成网站在线观看播放| 婷婷色麻豆天堂久久| 国精品久久久久久国模美| 在线亚洲精品国产二区图片欧美| 成人免费观看视频高清| 国产亚洲欧美精品永久| 侵犯人妻中文字幕一二三四区| 人人妻,人人澡人人爽秒播 | 亚洲一码二码三码区别大吗| 久久久精品94久久精品| 一级黄片播放器| 国产欧美日韩精品亚洲av| 国产成人91sexporn| 午夜福利影视在线免费观看| 亚洲精品久久成人aⅴ小说| 男女边摸边吃奶| 人人澡人人妻人| 亚洲精品国产区一区二| 国产欧美日韩精品亚洲av| 日韩大码丰满熟妇| 啦啦啦在线免费观看视频4| 中文字幕亚洲精品专区| 一级毛片黄色毛片免费观看视频| 日韩中文字幕视频在线看片| 在现免费观看毛片| av视频免费观看在线观看| 国产免费福利视频在线观看| 精品第一国产精品| 婷婷成人精品国产| 国产一区二区激情短视频 | 亚洲国产成人一精品久久久| 捣出白浆h1v1| 日韩一卡2卡3卡4卡2021年| 免费一级毛片在线播放高清视频 | 999精品在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 成年美女黄网站色视频大全免费| xxx大片免费视频| 成年人午夜在线观看视频| 少妇粗大呻吟视频| 久久亚洲国产成人精品v| 欧美精品一区二区免费开放| 青青草视频在线视频观看| 美女国产高潮福利片在线看| 国产亚洲av片在线观看秒播厂| 午夜福利免费观看在线| av福利片在线| 欧美日韩国产mv在线观看视频| 国产在线观看jvid| 狠狠婷婷综合久久久久久88av| 久久综合国产亚洲精品| 欧美日韩亚洲高清精品| 女性生殖器流出的白浆| 2018国产大陆天天弄谢| 久久久久国产精品人妻一区二区| 日本欧美视频一区| 这个男人来自地球电影免费观看| 一级毛片女人18水好多 | 日本a在线网址| 黑人猛操日本美女一级片| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| 精品第一国产精品| 悠悠久久av| 美女主播在线视频| 国产精品偷伦视频观看了| 9色porny在线观看| 青春草视频在线免费观看| a级毛片在线看网站| 日日摸夜夜添夜夜爱| 国产精品久久久av美女十八| 只有这里有精品99| 一本久久精品| 日本欧美视频一区| 亚洲美女黄色视频免费看| 五月天丁香电影| 国产一卡二卡三卡精品| 丝袜人妻中文字幕| 欧美乱码精品一区二区三区| 91成人精品电影| 男人操女人黄网站| 久久精品亚洲熟妇少妇任你| 国产三级黄色录像| 日韩一卡2卡3卡4卡2021年| 黄色视频在线播放观看不卡| 香蕉丝袜av| 美女扒开内裤让男人捅视频| 国产日韩欧美视频二区| 色网站视频免费| 两个人看的免费小视频| 欧美在线黄色| 香蕉国产在线看| 99国产精品一区二区三区| 亚洲中文av在线| 一边亲一边摸免费视频| 热re99久久精品国产66热6| 久久精品国产亚洲av高清一级| 欧美国产精品一级二级三级| 999久久久国产精品视频| 欧美亚洲 丝袜 人妻 在线| 国产免费视频播放在线视频| 丝袜在线中文字幕| 老鸭窝网址在线观看| xxxhd国产人妻xxx| 国产爽快片一区二区三区| 国产成人影院久久av| 色播在线永久视频| 精品国产一区二区三区四区第35| 日韩av在线免费看完整版不卡| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区| 波多野结衣一区麻豆| 99国产精品一区二区三区| 亚洲,欧美精品.| 亚洲人成电影免费在线| 国产免费一区二区三区四区乱码| 考比视频在线观看| 这个男人来自地球电影免费观看| 热re99久久国产66热| 亚洲视频免费观看视频| 性少妇av在线| 97精品久久久久久久久久精品| 久久天堂一区二区三区四区| 精品熟女少妇八av免费久了| 久久精品亚洲熟妇少妇任你| 色94色欧美一区二区| 一级片'在线观看视频| 欧美成人午夜精品| 成人免费观看视频高清| 国产精品三级大全| 亚洲精品国产av蜜桃| 男的添女的下面高潮视频| 亚洲 欧美一区二区三区| av天堂在线播放| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品| 亚洲一卡2卡3卡4卡5卡精品中文| 免费日韩欧美在线观看| 亚洲国产精品999| 午夜老司机福利片| a 毛片基地| 欧美人与善性xxx| av一本久久久久| av天堂在线播放| 麻豆av在线久日| 老司机深夜福利视频在线观看 | 老司机深夜福利视频在线观看 | 久久午夜综合久久蜜桃| 中文字幕亚洲精品专区| 亚洲成色77777| www.精华液| 亚洲国产成人一精品久久久| 肉色欧美久久久久久久蜜桃| 在线观看国产h片| 18禁观看日本| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区av网在线观看 | 欧美黄色淫秽网站| 久久久久国产精品人妻一区二区| 精品国产乱码久久久久久男人| 黑人猛操日本美女一级片| 一级a爱视频在线免费观看| av有码第一页| 日韩大片免费观看网站| 天堂俺去俺来也www色官网| 久久99一区二区三区| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 香蕉丝袜av| 9色porny在线观看| 精品人妻在线不人妻| 日本一区二区免费在线视频| 一二三四在线观看免费中文在| 成人亚洲欧美一区二区av| 国产高清国产精品国产三级| 美女中出高潮动态图| 国产国语露脸激情在线看| 建设人人有责人人尽责人人享有的| 色精品久久人妻99蜜桃| netflix在线观看网站| 亚洲av美国av| 国产97色在线日韩免费| 国产精品麻豆人妻色哟哟久久| 18禁裸乳无遮挡动漫免费视频| 婷婷成人精品国产| 纵有疾风起免费观看全集完整版| 久久久久网色| 国产精品av久久久久免费| 久久青草综合色| 91精品三级在线观看| 夫妻午夜视频| 国产深夜福利视频在线观看| 精品久久久久久电影网| 超碰97精品在线观看| 天堂中文最新版在线下载| 99九九在线精品视频| 黄色视频不卡| 一级毛片女人18水好多 | 日本av免费视频播放| 午夜激情av网站| 亚洲国产欧美一区二区综合| 老司机深夜福利视频在线观看 | 亚洲国产毛片av蜜桃av|