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

    全球副溶血弧菌環(huán)境菌株的基于堿基序列和肽鏈序列的多位點(diǎn)序列分型

    2015-12-26 09:43:39徐嘉良杜小莉
    食品科學(xué) 2015年21期
    關(guān)鍵詞:復(fù)合體弧菌分型

    徐嘉良,杜小莉,盧 昕,*

    (1.北京工商大學(xué)食品學(xué)院,食品添加劑與配料北京高校工程研究中心,北京 100048;2.中國(guó)疾病預(yù)防控制中心傳染病預(yù)防控制所,傳染病預(yù)防控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102206)

    全球副溶血弧菌環(huán)境菌株的基于堿基序列和肽鏈序列的多位點(diǎn)序列分型

    徐嘉良1,杜小莉2,盧 昕2,*

    (1.北京工商大學(xué)食品學(xué)院,食品添加劑與配料北京高校工程研究中心,北京 100048;2.中國(guó)疾病預(yù)防控制中心傳染病預(yù)防控制所,傳染病預(yù)防控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102206)

    為了了解來(lái)自全球的副溶血弧菌環(huán)境菌株的群組結(jié)構(gòu)和克隆復(fù)合體構(gòu)成,利用PubMLST公共數(shù)據(jù)的數(shù)據(jù),篩選其中具有完整序列分型(sequence type,ST)及肽序列型(peptide sequence type,pST)的來(lái)自全球的副溶血弧菌環(huán)境菌株數(shù)據(jù),進(jìn)行亞群分析和克隆復(fù)合體分析,并分別構(gòu)建了基于ST型和pST型的最小生成樹。結(jié)果表明,從數(shù)據(jù)庫(kù)中共篩選到具有ST型及pST型的來(lái)自全球的副溶血弧菌環(huán)境菌株數(shù)據(jù)886 條,共含有663 個(gè)ST型,以ST3型為最多,但也僅含27 條菌株信息。其中,有436 條菌株數(shù)據(jù)來(lái)自中國(guó)(包括大陸、香港和臺(tái)灣),覆蓋了410 個(gè)ST型,每個(gè)ST型含1~4 條菌株數(shù)據(jù),這436 條菌株數(shù)據(jù)又可分為128 個(gè)pST型,其中pST1型為最多,達(dá)116 條菌株數(shù)據(jù)。利用eBURST軟件對(duì)數(shù)據(jù)進(jìn)行分析,共發(fā)現(xiàn)了73 個(gè)組,483 個(gè)單體。STRUCTURE軟件分析顯示來(lái)自全球的副溶血弧菌環(huán)境菌株的適宜亞群數(shù)為4,各亞群內(nèi)的樣品平均距離為0.981 7。綜上結(jié)果表明,來(lái)自全球的副溶血弧菌環(huán)境菌株具有高度多態(tài)性,可更細(xì)分為4 個(gè)亞群,MLST分型時(shí)以ST3型為多,AA-MLST分型時(shí)以pST1型為主。

    副溶血弧菌;基于堿基序列的多位點(diǎn)序列分型;基于肽鏈的多位點(diǎn)序列分型;克隆復(fù)合體

    副溶血弧菌(Vibrio parahaemolyticus)屬于弧菌科弧菌屬,菌體呈弧狀、桿狀等形態(tài),革蘭染色呈陰性。1950年,F(xiàn)ujino在日本的一次食物中毒中首次發(fā)現(xiàn)并成功分離培養(yǎng)到副溶血弧菌[1]。食用含有副溶血弧菌的食物可能引發(fā)食源性胃腸道炎癥,少數(shù)情況下還可引發(fā)敗血癥,危及患者生命。該菌主要來(lái)自于海產(chǎn)品,常見者為蟹類、烏賊、海蜇、魚、黃泥螺等,其次為蛋品、肉類或蔬菜,已在全球頻繁引起食物中毒事件[2]。從遺傳角度和血清型角度看,副溶血弧菌都存在高度多態(tài)性。目前,副溶血弧菌包括13 種不同的O抗原和71 個(gè)不同的K抗原[3-4],但與霍亂弧菌(Vibrio cholera)不同,目前的研究顯示副溶血弧菌的致病性并不與特定的血清型相關(guān),目前有71 種血清型的副溶血弧菌均可以引起胃腸炎[5]。自1996年起,副溶血弧菌在全球地區(qū)開始暴發(fā),其中以O(shè)3:K6血清型的副溶血弧菌為主要的暴發(fā)菌株。之后,開始在全球迅速傳播。

    近年來(lái),隨著分子生物學(xué)和基因組學(xué)的發(fā)展,一些分子分型技術(shù),如擴(kuò)增片段長(zhǎng)度多態(tài)性分析技術(shù)(amplified fragment length polymorphism,AFLP)[6]、隨機(jī)擴(kuò)增多態(tài)性DNA片段(random amplified polymorphic DNA,RAPD)[7]、核糖體基因分型[8]、脈沖場(chǎng)凝膠電泳(pulsed-field gel electrophoresis,PFGE)[9]、多位點(diǎn)序列分型(multilocus sequence typing,MLST)[10]和多位點(diǎn)可變數(shù)目串聯(lián)重復(fù)序列分析(multilocus variable number tandem repeat analysis,MLVA)[11]等,已廣泛應(yīng)用于細(xì)菌鑒定、分型溯源、種群結(jié)構(gòu)、遺傳進(jìn)化分析及分子流行病學(xué)[12-13]的調(diào)查和食品安全溯源等方面研究[14]。MLST方法因其豐富的數(shù)據(jù)庫(kù)可研究不同地域不同來(lái)源的菌株遺傳變異情況,已被用于區(qū)分副溶血弧菌[15]。副溶血弧菌的MLST在2008年由González-Escalona等[16]提出,由7 個(gè)基因片段獲得菌株的序列分型(sequence type,ST)。隨后2014年,Urmersbach等[17]在副溶血弧菌的群體結(jié)構(gòu)分析中又提出了基于肽鏈序列的多位點(diǎn)分型(amino acid multilocus sequence typing,AA-MLST)方案。由于副溶血弧菌遺傳上的多態(tài)性,尤其是環(huán)境菌株呈現(xiàn)出高度多態(tài)性,因此應(yīng)用AA-MLST進(jìn)行菌株遺傳關(guān)系研究可能會(huì)得到更為可信的結(jié)果[17]。本研究將采用MLST方案和AA-MLST方案分別研究來(lái)自PubMLST數(shù)據(jù)庫(kù)中的來(lái)自全球地區(qū)的副溶血弧菌環(huán)境菌株的群體結(jié)構(gòu),為我國(guó)副溶血弧菌的預(yù)防控制提供參考,通過(guò)采取針對(duì)性的手段,減少致病性的環(huán)境菌株對(duì)食品的污染,為保障我國(guó)的食品安全提供理論依據(jù)。

    1 材料與方法

    1.1 數(shù)據(jù)

    全部數(shù)據(jù)均來(lái)自于截止到2015年2月10日之前的PubMLST公共數(shù)據(jù)庫(kù)(http://pubmlst.org/ vparahaemolyticus/),共有來(lái)自于全球副溶血弧菌的1 689 條數(shù)據(jù),其中包括臨床菌株(分離自病人)和環(huán)境菌株(分離自外環(huán)境水體和水產(chǎn)品中)。篩選其中的環(huán)境菌株,且具備ST型信息及肽序列型(peptide suquence type,pST)的數(shù)據(jù)條目,共篩選出886 條,用于后續(xù)分析。

    1.2 克隆群結(jié)構(gòu)分析

    采用eBURST version 3(http://eburst/mlst.net)對(duì)數(shù)據(jù)進(jìn)行克隆復(fù)合體的分析,分析中采用了1 000次Bootstrap重排。本研究中,限定7 個(gè)管家基因中必須有6 個(gè)及以上相同的歸為同分型組,若該組中含4 個(gè)及以上的ST型,該組即為一個(gè)克隆復(fù)合體(clonal complex,CC),并確定其中一個(gè)ST型在進(jìn)化關(guān)系上可能是祖先ST型,定為起源。不屬于任何1 個(gè)克隆復(fù)合體的ST型被稱為單體。

    1.3 STRUCTURE軟件分群

    應(yīng)用STRUCTURE2.3.4軟件,計(jì)算來(lái)自全球的副溶血弧菌環(huán)境菌株歸屬第k亞群的概率Q值。在確定亞群數(shù)時(shí),選擇獨(dú)立等位基因頻率模型,依次設(shè)定亞群數(shù)目K為1~11,將馬爾科夫鏈蒙特卡洛(Markov chain Monte Carlo,MCMC)迭代設(shè)為110 000 次,去掉前10 000 次,每個(gè)K重復(fù)獨(dú)立運(yùn)行4 次。

    1.4 Bionumerics分析

    采用Bionumerics軟件分別對(duì)MLST數(shù)據(jù)及AA-MLST數(shù)據(jù)進(jìn)行分析,構(gòu)建最小生成樹,以菌株來(lái)源國(guó)家或菌株來(lái)源地區(qū)(亞洲、拉丁美洲、歐洲、美洲)為顏色區(qū)分,圓圈大小代表菌株數(shù)量,研究不同國(guó)家來(lái)源菌株群體結(jié)構(gòu)和遺傳進(jìn)化關(guān)系,比較MLST方案及AA-MLST方案對(duì)數(shù)據(jù)的分析結(jié)果。

    2 結(jié)果與分析

    2.1 MLST數(shù)據(jù)基本分析

    從數(shù)據(jù)庫(kù)中共篩選到具有ST型及pST型的來(lái)自全球的副溶血弧菌環(huán)境菌株數(shù)據(jù)886 條,用于后續(xù)分析,國(guó)別(菌株數(shù)據(jù)量)分別為:孟加拉(7)、智利(19)、中國(guó)(包括大陸、香港和臺(tái)灣)(436)、丹麥(8)、厄瓜多爾(30)、德國(guó)(27)、洪都拉斯(2)、印度(3)、印度尼西亞(1)、日本(1)、韓國(guó)(1)、馬來(lái)西亞(1)、墨西哥(1)、荷蘭(3)、挪威(17)、西班牙(4)、斯里蘭卡(43)、泰國(guó)(116)、英國(guó)(1)、美國(guó)(161)、越南(4)。在886條數(shù)據(jù)中,共含有663 個(gè)ST型,其中ST3型(3、4、19、4、29、4、22)所含數(shù)據(jù)數(shù)最多,但也僅為27 條。在886 條菌株數(shù)據(jù)中,有436 條菌株數(shù)據(jù)來(lái)自中國(guó)(包括大陸、香港和臺(tái)灣),覆蓋了410 個(gè)ST型,每個(gè)ST型含1~4 條菌株數(shù)據(jù),這436 條菌株數(shù)據(jù)又可分為128 個(gè)pST型,其中pST1型為最多,達(dá)116 條菌株數(shù)據(jù)。

    2.2 eBURST克隆復(fù)合體分析

    采用eBURST分析,886 條菌株數(shù)據(jù)中,可以分為663 個(gè)ST型,包含了73 個(gè)組,483 個(gè)單體。其中CC34群為含ST型別最多的克隆復(fù)合體,包括了來(lái)自7 個(gè)ST型(ST34、ST334、ST35、ST77、ST324、ST95、ST33)的24 株菌株數(shù)據(jù),這其中ST34型被認(rèn)為是CC34群的祖先型。CC3群為含菌株數(shù)最多的克隆復(fù)合體,包括了來(lái)自4 個(gè)ST型(ST3、ST266、ST220、ST2)的30 株菌株數(shù)據(jù),其中ST3型被認(rèn)為是CC3群的祖先型。

    2.3 STRUCTURE對(duì)亞群的劃分

    STRUCTURE軟件分析結(jié)果顯示,LnP(D)和α1在亞群數(shù)為4時(shí)有明顯的拐點(diǎn),因而確定來(lái)自全球的副溶血弧菌環(huán)境菌株群體的適宜亞群數(shù)為4(圖1)。各亞群之間的F值變化范圍是0.010 1~0.012 2,平均值為0.011 1。STRUCTURE軟件分群各亞群內(nèi)樣品的平均距離為0.980 6~0.984 0,均值為0.981 7。

    圖1 全球副溶血弧菌環(huán)境菌株群體歸屬于4 個(gè)亞群的后驗(yàn)概率圖Fig.1 A posteriori probability graph of worldwide environmental Vibrio parahaemolyticus strains belonging to 4 subgroups

    2.4 AA-MLST數(shù)據(jù)基本分析

    從數(shù)據(jù)庫(kù)中共篩選到具有pST型的來(lái)自全球的副溶血弧菌環(huán)境菌株數(shù)據(jù)886 條,共含有187 個(gè)pST型,其中pST1型(1、1、1、1、1、1、1)所含數(shù)據(jù)數(shù)最多,可達(dá)228 條。pST2型(2、1、1、1、1、1、1)所含數(shù)據(jù)數(shù)次之,可達(dá)188 條。根據(jù)pST型分析克隆群結(jié)構(gòu),顯示pST75、pST92、pST221、pST178、pST153、pST79、pST148、pST125、pST142、pST147和pST175與其他pST型有多于1 個(gè)等位基因的差異。

    2.5 最小生成樹的構(gòu)建

    Bionumerics軟件應(yīng)用Prim算法,根據(jù)各等位基因型之間的差異對(duì)不同ST型進(jìn)行分組聚類,找到各ST型進(jìn)化關(guān)系上共同的祖先,最終生成最小生成樹(minimum spanning tree,MST)。在MST圖中(圖2、3),每個(gè)圓代表一種ST型,圓圈大小與菌株數(shù)目成正比。比較eBURST軟件與Bionumerics軟件的分析結(jié)果發(fā)現(xiàn),兩種算法對(duì)克隆復(fù)合體及祖先ST型的認(rèn)定存在一定差異,而且Bionumerics軟件對(duì)基于MLST數(shù)據(jù)和AA-MLST數(shù)據(jù)的克隆復(fù)合體分析也存在差異。并不是同屬于一個(gè)ST型克隆復(fù)合體的菌株也會(huì)劃歸于同一pST型克隆復(fù)合體中。世界各地報(bào)道的副溶血弧菌環(huán)境菌株呈現(xiàn)出高度多態(tài)性,并沒(méi)有特別優(yōu)勢(shì)的克隆群存在,且副溶血弧菌的血清型分型與MLST分型及AA-MLST分型結(jié)果均表現(xiàn)出不同,基于ST型的分析找到比基于pST型數(shù)據(jù)分析更多的克隆復(fù)合體。

    圖2 全球副溶血弧菌886 株環(huán)境菌株基于MLST分型方案數(shù)據(jù)的最小生成樹Fig.2 Minimum spanning tree based on MLST of 886 worldwide environmental Vibrio parahaemolyticus strains

    由圖3可知,有些ST型的分布比較局限,僅分布于一個(gè)國(guó)家或一個(gè)大洲,而有些ST型則廣泛分布于多地,如ST3型在美洲、拉丁美洲、亞洲均有分布;ST711型分布于拉丁美洲和亞洲;ST539型分布于拉丁美洲和美洲;ST471型分布于歐洲和亞洲;ST424型分布于亞洲和拉丁美洲;ST246分布于拉丁美洲和亞洲;ST141型也歐洲和美洲均被發(fā)現(xiàn);ST121型目前在歐洲和亞洲均被檢測(cè)到;ST114型為美洲和亞洲分布;ST56型為拉丁美洲和亞洲分布;ST28型分布于拉丁美洲和歐洲;ST12型分布于美洲和亞洲;ST6型分布于拉丁美洲和歐洲。亞洲有507 個(gè)ST型,歐洲有38 個(gè)ST型,美洲有94 個(gè)ST型別,拉丁美洲有36 個(gè)ST型別。由此可見,亞洲的副溶血弧菌環(huán)境分離株的型別和多樣性要顯著高于其他洲。

    圖3 全球副溶血弧菌886 株環(huán)境菌株基于AA-MLST分型方案數(shù)據(jù)的最小生成樹Fig.3 Minimum spanning tree based on AA-MLST of 886 worldwide environmental Vibrio parahaemolyticus strains

    3 討 論

    國(guó)內(nèi)外副溶血弧菌食物中毒暴發(fā)主要原因與生食海鮮和海產(chǎn)品污染有關(guān)[18]。也包括肉制品、腌制品和一些涼菜等[19],由副溶血弧菌引起的食物中毒事件居細(xì)菌性食物中毒的首位[20-21],開展全球副溶血弧菌環(huán)境菌株的群體遺傳學(xué)特征研究對(duì)中國(guó)的副溶血弧菌的預(yù)防控制具有重要意義。在副溶血弧菌中,血清學(xué)在追溯流行性克隆復(fù)合體傳播上是不可靠的[16]。而MLST是以選定的管家基因序列分析為基礎(chǔ)的,可更好地揭示同一個(gè)屬中菌株間的遺傳親緣關(guān)系[22]。2014年,有研究者對(duì)來(lái)自國(guó)際上17 個(gè)沿海城市的副溶血弧菌臨床菌株進(jìn)行了MLST分析,發(fā)現(xiàn)490 株菌可分為161 個(gè)ST型,8 個(gè)克隆群復(fù)合體,其中以CC3型為主[23]。本研究采用來(lái)自PubMLST公共數(shù)據(jù)庫(kù)中的來(lái)自全球地區(qū)的副溶血弧菌環(huán)境菌株,進(jìn)行了MLST及AA-MLST分型方案分析,研究全球副溶血弧菌環(huán)境菌株的群體結(jié)構(gòu)和克隆群復(fù)合體組成。從數(shù)據(jù)庫(kù)中共篩選到具有ST型及pST型的來(lái)自全球的副溶血弧菌環(huán)境菌株數(shù)據(jù)886 條,共含有663 個(gè)ST型,這說(shuō)明來(lái)自全球地區(qū)的副溶血弧菌環(huán)境菌株的高度多態(tài)性。其中有436 條菌株數(shù)據(jù)來(lái)自中國(guó)(包括大陸、香港和臺(tái)灣),覆蓋了410 個(gè)ST型。由此可以看出,副溶血弧菌的環(huán)境菌株比臨床菌株體現(xiàn)出更高的遺傳多樣性。因?yàn)楦比苎【h(huán)境菌株的遺傳高度多態(tài)性,因此確定菌株之間遺傳相關(guān)可能缺乏可信度,因此基于肽鏈順序的AA-MLST方案可能更適合副溶血弧菌的群組結(jié)構(gòu)分析[24],于是在本研究中也采用了基于pST型的分析。從結(jié)果中,可以發(fā)現(xiàn)有些型別是跨洲分布的,這提示這些型別可能更適合外界壞境,有著更廣泛的分布,人們?cè)谕猸h(huán)境監(jiān)測(cè)時(shí)需格外注意這些型別的發(fā)現(xiàn)。亞洲的副溶血弧菌環(huán)境分離株的型別和多樣性要顯著高于其他洲,這提示作為亞洲國(guó)家,應(yīng)該密切關(guān)注副溶血弧菌在環(huán)境中的進(jìn)化變異。與其他研究結(jié)果一致[21],本結(jié)果也顯示副溶血弧菌的血清型分型對(duì)其流行病學(xué)研究?jī)r(jià)值有限。因此,本研究認(rèn)為在副溶血弧菌的監(jiān)測(cè)中,應(yīng)將分子分型納入到菌株分型中,在全球副溶血弧菌環(huán)境菌株整體呈現(xiàn)高度多態(tài)性的大環(huán)境下,也應(yīng)該引起我國(guó)專業(yè)技術(shù)人員和食品安全監(jiān)督部門的警示,對(duì)于環(huán)境中的副溶血弧菌可能引起的養(yǎng)殖和食品加工過(guò)程中的污染,尤其是海產(chǎn)品,比如副溶血弧菌傳播的主要媒介:蝦、蟹、牡蠣等海產(chǎn)品污染食品中的副溶血弧菌要加大監(jiān)測(cè)力度,此外,在監(jiān)測(cè)工作中的重點(diǎn)是關(guān)注不同克隆群的副溶血弧菌,加強(qiáng)對(duì)環(huán)境污染的監(jiān)測(cè)和管理,這樣可及時(shí)掌控分離自環(huán)境中的副溶血弧菌的群體遺傳構(gòu)成、優(yōu)勢(shì)菌群變化,預(yù)警致病類似菌株在海產(chǎn)品中的分布變化以及可能的危險(xiǎn)食品,為暴發(fā)疫情和人群中流行菌株的溯源提供基礎(chǔ)數(shù)據(jù),采取必要的措施減少食源性疾病暴發(fā)和流行[25]。

    [1] FUJINO K. Discovery of Vibrio parahaemolyticus[J]. Nihon Saikingaku Zasshi, 1951, 21(8): 373-388.

    [2] ALAM M, ChOwDhURy w B, BhUIyAN N A, et al. Serogroup, virulence, and genetic traits of Vibrio parahaemolyticus in the estuarine ecosystem of Bangladesh[J]. Applied and Environmental Microbiology, 2009, 75(19): 6268-6274.

    [3] 李維彬, 郭辰虹, 王玉炳. 病原微生物快速檢測(cè)方法及研究進(jìn)展[J].生物技術(shù)通報(bào), 2006(2): 67-70.

    [4] hARAKUDO y, DEPAOLA A. Pathogenic vibrios and food safety[M]. New York: Nova Science Publishers Inc., 2012: 179-210.

    [5] NAIR G B, HORMAZáBAL J C. The Vibrio parahaemolyticus pandemic[J]. Revista Chilena de Infectologia, 2005, 22(2): 125-130.

    [6] VANDENBERGHE J, VERDONCK L, ROBLES-AROZARENA R, et al. Vibrios associated with Litopenaeusvannamei larvae, postlarvae, broodstock, and hatchery probionts[J]. Applied and Environmental Microbiology, 1999, 65(6): 2592-2597.

    [7] SADOK K, MEJDI S, NOURhEN S, et al. Phenotypic characterization and RAPD fingerprinting of Vibrio parahaemolyticus and Vibrio alginolyticus isolated during Tunisian fish farm outbreaks[J]. Folia Microbiology (Praha), 2013, 58(1): 17-26.

    [8] CHEN Wanyi, XIE Yanping, XU Jingye, et al. Molecular typing of Vibrio parahaemolyticus isolates from the middle-east coastline of China[J]. International Journal of Food Microbiology, 2012, 153(3): 402-412.

    [9] hAENDIGES J, ROCK M, MyERS R A, et al. Pandemic Vibrio parahaemolyticus, Maryland, USA, 2012[J]. Emerging Infectious Diseases, 2014, 20(4): 718-720.

    [10] HAN Dongsheng, TANG Hui, REN Chuanli, et al. Prevalence and genetic diversity of clinical Vibrio parahaemolyticus isolates from China, revealed by multilocus sequence typing scheme[J]. Front Microbiology, 2015, 6: 291. doi:10.3389/fmicb.2015.00291.

    [11] KIMURA B, SEKINE y, TAKAhAShI h, et al. Multiplelocus variable-number of tandem-repeats analysis distinguishes Vibrio parahaemolyticuspandemic O3:K6 strains[J]. Journal of Microbiological Methods, 2008, 72(3): 313-320.

    [12] 張蔚, 潘勁草, 陳坤. 副溶血性弧菌的分子流行病學(xué)研究進(jìn)展[J]. 國(guó)外醫(yī)學(xué): 流行病學(xué)傳染病學(xué)分冊(cè), 2004, 31(3): 182-184.

    [13] 魯開國(guó), 史英輝. 副溶血弧菌的分子生物學(xué)研究進(jìn)展[J]. 中國(guó)衛(wèi)生檢驗(yàn)雜志, 1994, 4(5): 314-316.

    [14] 劉丁. 細(xì)菌基因分型技術(shù)方法及其評(píng)價(jià)[J]. 國(guó)外醫(yī)學(xué): 臨床生物化學(xué)與檢驗(yàn)學(xué)分冊(cè), 2002(4): 220-222.

    [15] TURNER J, PARANJPyE R N, LANDIS E D, et al. Population structure of clinical and environmental Vibrio parahaemolyticus from the Pacific Northwest coast of the United States[J]. PLoS ONE, 2013, 8(2): e55726. doi: 10.1371/journal.pone.0055726.

    [16] GONZáLEZ-ESCALONA N, MARTINEZ-URTAZA J, ROMERO J, et al. Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing[J]. Journal of Bacteriology, 2008, 190(8): 2831-2840.

    [17] URMERSBACh S, ALTER T, KORALAGE M S, et al. Population analysis of Vibrio parahaemolyticus originating from different geographical regions demonstrates a high genetic diversity[J]. BMC Microbiology, 2014, 14(3): 106-111.

    [18] 劉秀梅, 陳艷, 郭云昌, 等. 2005年中國(guó)食源性疾病暴發(fā)事件監(jiān)測(cè)資料分析[J]. 中國(guó)食品衛(wèi)生雜志, 2008, 20(6): 506-509.

    [19] 林祥田. 94起副溶血弧菌食物中毒流行病學(xué)特征分析[J]. 上海預(yù)防醫(yī)學(xué), 2006, 18(3): 141-143.

    [20] 李薇薇, 梅玲玲, 唐震, 等. 2007—2009年中國(guó)副溶血性弧菌臨床分離株分子特征分析[J]. 中華預(yù)防醫(yī)學(xué)雜志, 2014, 48(1): 44-52.

    [21] 李婧, 邱少富, 宋宏彬, 等. 沈陽(yáng)地區(qū)副溶血弧菌毒力基因檢測(cè)及脈沖場(chǎng)凝膠電泳分析[J]. 中國(guó)公共衛(wèi)生, 2012, 28(5): 696-698.

    [22] ELLIS C N, SChUSTER B M, STRIPLIN M J, et al. Influence of seasonality on the genetic diversity of Vibrio parahaemolyticus in New Hampshire shellfish waters as determined by multilocus sequence analysis[J]. Applied and Environmental Microbiology, 2012, 78(10): 3778-3782.

    [23] HAN Dongsheng, TANG Hui, LU Jingye, et al. Population structure of clinical Vibrio parahaemolyticus from 17 coastal countries, determined through multilocus sequence analysis[J]. PLoS ONE, 2014, 9(9): e107371. doi:10.1371/journal.pone.0107371.

    [24] YAN Yanfeng, CUI Yujun, HAN Haihong, et al. Extended MLST-based population genetics and phylogeny of Vibrio parahaemolyticus with high levels of recombination[J]. International Journal of Food Microbiology, 2011, 145(1): 106-112.

    [25] 樊曉琳, 張德福, 付緒磊, 等. 錦州筆架山海域海產(chǎn)品中副溶血弧菌的分子流行病學(xué)調(diào)查[J]. 食品工業(yè)科技, 2014, 35(16): 49-51.

    Multilocus Sequence Typing of Global Environmental Strains of Vibrio parahaemolyticus Based on Nucleotide and Peptide Sequences

    XU Jialiang1, DU Xiaoli2, LU Xin2,*
    (1. Beijing Higher Institution Engineering Research Center of Food Additives and Ingredients, School of Food and Chemical Engineering, Beijing Technology and Business University, Beijing 100048, China; 2. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 102206, China)

    Objective: To understand the group structure and clonal complex of worldwide environmental Vibrio parahaemolyticus strains. Methods: Worldwide environmental V. parahaemolyticus strains with complete sequence types (ST) and peptide sequence types (pST) were screened from the PubMLST public database, the subgroup analysis and complex analysis was conducted, and the minimum spanning tree was established based on ST and pST types. Results: From the database, 886 items were screened for analyzing the ST and pST types of worldwide environmental V. parahaemolyticus strains, which belonged to 663 ST types with ST3 type being the most dominant, but only containing information on 27 strains. Out of these, 436 items came from China (including mainland, Hong Kong and Taiwan), including 410 ST types, each of which contained information on 1-4 strains and these 436 strains could be divided into 128 pST types, with pST1 type being the most dominant type, containing information on 116 strains. The eBURST software was used to analyze the data and 73 groups and 483 singletons were observed. Analysis with the software STRUCTURE showed that appropriate subset number of global V. parahaemolyticus strains was 4, and the average distance between samples in each subgroup was 0.981 7. Conclusion: Worldwide environmental V. parahaemolyticus strains in Asia are characterized by high polymorphism and can be subdivided into 4 subgroups. ST3 type is the most dominant strain when MLST typing is used while pST1 type is the most dominant strain by AA-MLST typing.

    Vibrio parahaemolyticus; multilocus sequence typing; peptide based multilocus sequence typing; clonal complex

    Q939.99

    A

    1002-6630(2015)21-0095-05

    10.7506/spkx1002-6630-201521019

    2015-05-18

    國(guó)家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(31400161);北京工商大學(xué)青年教師科研啟動(dòng)基金項(xiàng)目(QNJJ2014-25)

    徐嘉良(1984—),男,講師,博士,研究方向?yàn)槭称焚|(zhì)量安全控制。E-mail:xujialiang@btbu.edu.cn

    *通信作者:盧昕(1981—),女,副研究員,博士,研究方向?yàn)槟c道致病微生物的進(jìn)化分析。E-mail:luxin@icdc.cn

    猜你喜歡
    復(fù)合體弧菌分型
    銷量增長(zhǎng)200倍!“弧菌克星”風(fēng)靡行業(yè),3天殺滅98%弧菌
    副溶血弧菌檢測(cè)方法的研究進(jìn)展
    失眠可調(diào)養(yǎng),食補(bǔ)需分型
    如何有效防控對(duì)蝦養(yǎng)殖中的弧菌病
    便秘有多種 治療須分型
    副溶血弧菌噬菌體微膠囊的制備及在餌料中的應(yīng)用
    CoFe2O4/空心微球復(fù)合體的制備與吸波性能
    基于分型線驅(qū)動(dòng)的分型面設(shè)計(jì)研究
    頸椎病中醫(yī)辨證分型與影像探討
    3種多糖復(fù)合體外抗腫瘤協(xié)同增效作用
    国产片特级美女逼逼视频| 99热网站在线观看| 成人毛片a级毛片在线播放| 波野结衣二区三区在线| 91精品一卡2卡3卡4卡| 熟妇人妻久久中文字幕3abv| 级片在线观看| 免费看av在线观看网站| 99久久精品国产国产毛片| 日韩成人av中文字幕在线观看| 欧美精品国产亚洲| 色综合亚洲欧美另类图片| 22中文网久久字幕| 99视频精品全部免费 在线| 欧美日韩国产亚洲二区| 国产高清三级在线| 国产乱人偷精品视频| 性欧美人与动物交配| 能在线免费观看的黄片| 男插女下体视频免费在线播放| 天美传媒精品一区二区| 一级毛片电影观看 | 久久亚洲国产成人精品v| 最近最新中文字幕大全电影3| 精品99又大又爽又粗少妇毛片| 久久人人爽人人片av| 亚洲欧洲日产国产| 亚洲综合色惰| eeuss影院久久| 又粗又爽又猛毛片免费看| 欧美日韩综合久久久久久| 国产精品国产高清国产av| 国产精品人妻久久久久久| 99热这里只有是精品50| 亚洲不卡免费看| 麻豆av噜噜一区二区三区| 午夜爱爱视频在线播放| 在线观看午夜福利视频| av视频在线观看入口| 丰满乱子伦码专区| 亚洲精品粉嫩美女一区| 爱豆传媒免费全集在线观看| 男女下面进入的视频免费午夜| 久久99蜜桃精品久久| 日韩一本色道免费dvd| 波野结衣二区三区在线| 欧美日韩国产亚洲二区| 亚洲欧美日韩高清在线视频| 免费黄网站久久成人精品| 亚洲成人久久爱视频| 老司机影院成人| а√天堂www在线а√下载| 欧美一区二区国产精品久久精品| 97超视频在线观看视频| 久久精品久久久久久噜噜老黄 | 亚洲性久久影院| 91狼人影院| 岛国毛片在线播放| 日本黄色视频三级网站网址| 一个人看的www免费观看视频| 国产成人一区二区在线| 一个人看视频在线观看www免费| 国内少妇人妻偷人精品xxx网站| av又黄又爽大尺度在线免费看 | 国产av不卡久久| 插阴视频在线观看视频| 成人亚洲精品av一区二区| 亚洲无线观看免费| 哪个播放器可以免费观看大片| 日韩强制内射视频| 久久99热6这里只有精品| 成人美女网站在线观看视频| 高清午夜精品一区二区三区 | 国产v大片淫在线免费观看| 欧美日本视频| 久久精品久久久久久噜噜老黄 | 欧美高清性xxxxhd video| 欧美不卡视频在线免费观看| 日韩三级伦理在线观看| 久久精品国产亚洲网站| 尾随美女入室| 精品久久久久久久久亚洲| av视频在线观看入口| 成年av动漫网址| 看十八女毛片水多多多| 波多野结衣高清作品| 久久久色成人| 啦啦啦啦在线视频资源| 精品99又大又爽又粗少妇毛片| 美女xxoo啪啪120秒动态图| 欧美性感艳星| 禁无遮挡网站| 91精品国产九色| 极品教师在线视频| 国产日韩欧美在线精品| 大型黄色视频在线免费观看| 一夜夜www| 如何舔出高潮| 91狼人影院| 亚洲av中文av极速乱| 亚洲18禁久久av| 一本久久中文字幕| 亚洲久久久久久中文字幕| 人人妻人人看人人澡| 91精品一卡2卡3卡4卡| 九九在线视频观看精品| 久久久午夜欧美精品| 亚洲精品色激情综合| 我要看日韩黄色一级片| 精品欧美国产一区二区三| 亚洲久久久久久中文字幕| 亚洲国产日韩欧美精品在线观看| 国产熟女欧美一区二区| 男女下面进入的视频免费午夜| 黄色视频,在线免费观看| 乱系列少妇在线播放| 国产在视频线在精品| 亚洲在线自拍视频| 久久久久久久久中文| 日韩国内少妇激情av| 国产免费一级a男人的天堂| 永久网站在线| 欧美日韩在线观看h| 久久精品国产亚洲av香蕉五月| 亚洲18禁久久av| 99久久无色码亚洲精品果冻| 欧美极品一区二区三区四区| 啦啦啦观看免费观看视频高清| 精品少妇黑人巨大在线播放 | 日韩一区二区视频免费看| 两性午夜刺激爽爽歪歪视频在线观看| 国内精品美女久久久久久| 免费av观看视频| 成年免费大片在线观看| 欧美xxxx黑人xx丫x性爽| 成年av动漫网址| 亚洲在线观看片| 99久久九九国产精品国产免费| 成年女人看的毛片在线观看| 日本五十路高清| av国产免费在线观看| 国国产精品蜜臀av免费| 给我免费播放毛片高清在线观看| 欧美高清性xxxxhd video| 天堂av国产一区二区熟女人妻| 国产一级毛片七仙女欲春2| 久久久久久久久大av| 欧美xxxx黑人xx丫x性爽| 蜜桃亚洲精品一区二区三区| 看片在线看免费视频| 日本免费一区二区三区高清不卡| 日韩三级伦理在线观看| 午夜福利视频1000在线观看| 免费黄网站久久成人精品| 尾随美女入室| 成人高潮视频无遮挡免费网站| 22中文网久久字幕| 久久久久久久久久黄片| 亚洲精品国产av成人精品| 国内久久婷婷六月综合欲色啪| 搞女人的毛片| 亚洲经典国产精华液单| 国产成人91sexporn| 一级毛片电影观看 | 久久精品国产清高在天天线| 嫩草影院入口| 最近的中文字幕免费完整| 国产v大片淫在线免费观看| 国产又黄又爽又无遮挡在线| 久久这里有精品视频免费| 久久鲁丝午夜福利片| 国产爱豆传媒在线观看| 一夜夜www| 级片在线观看| 深夜a级毛片| videossex国产| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 欧美激情在线99| 欧美成人免费av一区二区三区| 尾随美女入室| 一本久久中文字幕| 好男人在线观看高清免费视频| 看片在线看免费视频| 五月伊人婷婷丁香| 毛片女人毛片| 久久久久国产网址| 国产一区二区在线av高清观看| 一级av片app| 欧美日韩综合久久久久久| 成人特级av手机在线观看| 国产精品蜜桃在线观看 | 免费观看人在逋| 乱码一卡2卡4卡精品| 大型黄色视频在线免费观看| 国产精品日韩av在线免费观看| 一级毛片电影观看 | 色尼玛亚洲综合影院| 美女高潮的动态| 直男gayav资源| 97超碰精品成人国产| 99视频精品全部免费 在线| 少妇丰满av| 欧美在线一区亚洲| 国产精品电影一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 美女高潮的动态| 欧美区成人在线视频| 综合色av麻豆| 观看免费一级毛片| 免费看a级黄色片| 插阴视频在线观看视频| 黄色配什么色好看| 我要看日韩黄色一级片| 99热精品在线国产| 亚洲av.av天堂| 久久精品夜色国产| 国产精品人妻久久久久久| 欧美激情久久久久久爽电影| 蜜桃亚洲精品一区二区三区| 男女那种视频在线观看| 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 一夜夜www| 成人美女网站在线观看视频| 亚洲精品成人久久久久久| 亚洲四区av| 两性午夜刺激爽爽歪歪视频在线观看| av又黄又爽大尺度在线免费看 | 在线免费观看不下载黄p国产| 不卡视频在线观看欧美| 综合色av麻豆| 啦啦啦韩国在线观看视频| 婷婷色综合大香蕉| 男人的好看免费观看在线视频| 亚洲国产日韩欧美精品在线观看| 国产成人福利小说| 国产在线精品亚洲第一网站| 亚洲欧洲国产日韩| 欧美丝袜亚洲另类| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 91精品一卡2卡3卡4卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产白丝娇喘喷水9色精品| 成年女人永久免费观看视频| 国产精品麻豆人妻色哟哟久久 | 久久久成人免费电影| 非洲黑人性xxxx精品又粗又长| 亚洲经典国产精华液单| 国产黄片视频在线免费观看| 免费av观看视频| 日韩欧美三级三区| 国产精品永久免费网站| 此物有八面人人有两片| 国产老妇伦熟女老妇高清| 看免费成人av毛片| 听说在线观看完整版免费高清| 高清毛片免费看| 小说图片视频综合网站| 蜜臀久久99精品久久宅男| 国产精品一区二区三区四区免费观看| 1000部很黄的大片| 国产精品国产三级国产av玫瑰| 国产精品久久久久久av不卡| 色视频www国产| 18+在线观看网站| 黄片无遮挡物在线观看| 欧美色视频一区免费| 黄色一级大片看看| 我要搜黄色片| 啦啦啦观看免费观看视频高清| 国产高清激情床上av| 欧美色欧美亚洲另类二区| av.在线天堂| 中国国产av一级| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 婷婷色av中文字幕| 国内少妇人妻偷人精品xxx网站| 成人无遮挡网站| 天天一区二区日本电影三级| 青青草视频在线视频观看| 黄色欧美视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲在久久综合| 国产一区二区在线av高清观看| 一本精品99久久精品77| 日韩精品有码人妻一区| 国产日本99.免费观看| 成人毛片60女人毛片免费| 毛片一级片免费看久久久久| 亚洲国产精品合色在线| 嫩草影院精品99| 久久久精品大字幕| 日韩三级伦理在线观看| 真实男女啪啪啪动态图| 91麻豆精品激情在线观看国产| 精品日产1卡2卡| 日日啪夜夜撸| 日韩亚洲欧美综合| 欧美日韩一区二区视频在线观看视频在线 | 少妇人妻精品综合一区二区 | 欧美zozozo另类| 小蜜桃在线观看免费完整版高清| 成年女人永久免费观看视频| 可以在线观看毛片的网站| 在线免费观看的www视频| 99热这里只有精品一区| 精品久久久久久久久亚洲| 人妻系列 视频| 久久久久久九九精品二区国产| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看| 成人午夜高清在线视频| 久久鲁丝午夜福利片| 成人漫画全彩无遮挡| 美女被艹到高潮喷水动态| 国产私拍福利视频在线观看| 老师上课跳d突然被开到最大视频| 男女下面进入的视频免费午夜| 性欧美人与动物交配| 午夜福利在线在线| 精品人妻熟女av久视频| 午夜福利在线观看免费完整高清在 | 国产欧美日韩精品一区二区| 中文精品一卡2卡3卡4更新| 级片在线观看| 色吧在线观看| 国产色爽女视频免费观看| 色综合亚洲欧美另类图片| 精品99又大又爽又粗少妇毛片| 少妇人妻一区二区三区视频| 91精品国产九色| 亚洲欧美精品专区久久| 免费在线观看成人毛片| 亚洲真实伦在线观看| 能在线免费观看的黄片| 久久精品人妻少妇| 欧美xxxx性猛交bbbb| 精品一区二区三区人妻视频| 插阴视频在线观看视频| 久久精品91蜜桃| 午夜福利在线观看吧| 白带黄色成豆腐渣| 久久韩国三级中文字幕| 亚洲av二区三区四区| 亚洲七黄色美女视频| 欧美日韩国产亚洲二区| 国产黄色小视频在线观看| 久久亚洲国产成人精品v| 综合色av麻豆| 亚洲av中文字字幕乱码综合| 国产精品av视频在线免费观看| 免费在线观看成人毛片| 久久亚洲国产成人精品v| 亚洲欧美日韩东京热| 九九久久精品国产亚洲av麻豆| 一本一本综合久久| 久久鲁丝午夜福利片| 国产高清视频在线观看网站| 夜夜夜夜夜久久久久| 女同久久另类99精品国产91| 在线观看美女被高潮喷水网站| 真实男女啪啪啪动态图| 免费看a级黄色片| 一本一本综合久久| 九草在线视频观看| 欧美性感艳星| 精品久久久噜噜| 久久精品国产自在天天线| 久久精品夜色国产| 欧美最新免费一区二区三区| 久久久久久大精品| 日本黄色片子视频| 国内精品美女久久久久久| 国产黄色小视频在线观看| 久久精品综合一区二区三区| 国产精品久久久久久av不卡| 成人午夜精彩视频在线观看| 男人舔女人下体高潮全视频| 欧美xxxx性猛交bbbb| 国产精品1区2区在线观看.| 中文精品一卡2卡3卡4更新| 婷婷色av中文字幕| 99视频精品全部免费 在线| 狂野欧美激情性xxxx在线观看| 亚洲欧美精品综合久久99| 日本-黄色视频高清免费观看| 欧美性猛交╳xxx乱大交人| 国产亚洲av片在线观看秒播厂 | 啦啦啦观看免费观看视频高清| 又粗又硬又长又爽又黄的视频 | 美女国产视频在线观看| 老女人水多毛片| 美女黄网站色视频| 黄色欧美视频在线观看| 久久久精品94久久精品| 国产伦在线观看视频一区| 一级二级三级毛片免费看| 亚洲无线观看免费| 欧美激情在线99| 久久久久久久久久久免费av| 一本久久中文字幕| 边亲边吃奶的免费视频| 麻豆乱淫一区二区| 久久国产乱子免费精品| 免费黄网站久久成人精品| 国产精品一区二区在线观看99 | 日本五十路高清| 丝袜美腿在线中文| 麻豆一二三区av精品| 精品国产三级普通话版| 国产成人aa在线观看| 波野结衣二区三区在线| 国产精品精品国产色婷婷| 国产色爽女视频免费观看| 搡女人真爽免费视频火全软件| 国内久久婷婷六月综合欲色啪| 国产私拍福利视频在线观看| 中文字幕久久专区| 嫩草影院精品99| 国产精品人妻久久久久久| 美女xxoo啪啪120秒动态图| 一个人看视频在线观看www免费| 99国产精品一区二区蜜桃av| 三级男女做爰猛烈吃奶摸视频| 日本一二三区视频观看| 亚洲成人久久性| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 高清毛片免费看| 久久中文看片网| 联通29元200g的流量卡| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 夜夜爽天天搞| 一级黄色大片毛片| 亚洲国产精品成人久久小说 | 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 中文字幕久久专区| 一级二级三级毛片免费看| 人妻制服诱惑在线中文字幕| 国产毛片a区久久久久| 日韩av在线大香蕉| 天堂中文最新版在线下载 | 白带黄色成豆腐渣| 99久国产av精品国产电影| 国产精品综合久久久久久久免费| 日本黄色片子视频| 国产精品av视频在线免费观看| 十八禁国产超污无遮挡网站| 国产精品日韩av在线免费观看| 成人永久免费在线观看视频| 午夜福利在线观看吧| 日本黄大片高清| 美女 人体艺术 gogo| 深夜a级毛片| 亚洲不卡免费看| 99热只有精品国产| 久久久久久久久中文| 日本av手机在线免费观看| 看十八女毛片水多多多| 欧美色视频一区免费| 欧美又色又爽又黄视频| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 99精品在免费线老司机午夜| 美女xxoo啪啪120秒动态图| 欧美不卡视频在线免费观看| 亚洲四区av| 久久欧美精品欧美久久欧美| 99久国产av精品| 婷婷六月久久综合丁香| 欧美潮喷喷水| 最近中文字幕高清免费大全6| 国产高清不卡午夜福利| 国产美女午夜福利| 日本色播在线视频| 国产高清三级在线| 青春草亚洲视频在线观看| 岛国毛片在线播放| 久久久色成人| 亚洲成人精品中文字幕电影| 在线播放国产精品三级| 久久久久久久午夜电影| 天天一区二区日本电影三级| 春色校园在线视频观看| 国产亚洲5aaaaa淫片| 亚洲av电影不卡..在线观看| 亚洲国产欧美在线一区| 国产伦一二天堂av在线观看| 日韩高清综合在线| 亚洲精品亚洲一区二区| 亚洲久久久久久中文字幕| 亚洲国产欧洲综合997久久,| 欧美人与善性xxx| 在线观看午夜福利视频| 此物有八面人人有两片| 亚洲精品乱码久久久v下载方式| 精品午夜福利在线看| 免费av观看视频| 卡戴珊不雅视频在线播放| 国产一级毛片七仙女欲春2| 久久99热6这里只有精品| 久久99热这里只有精品18| 成人一区二区视频在线观看| 午夜亚洲福利在线播放| 亚洲第一区二区三区不卡| 国产日本99.免费观看| 婷婷六月久久综合丁香| 精品久久久久久久人妻蜜臀av| 亚洲av熟女| 色哟哟哟哟哟哟| 婷婷色av中文字幕| 深夜精品福利| 色5月婷婷丁香| 日韩高清综合在线| 日产精品乱码卡一卡2卡三| 亚洲图色成人| 91久久精品国产一区二区三区| 免费无遮挡裸体视频| 亚洲欧美中文字幕日韩二区| 免费观看a级毛片全部| 久久久久久久亚洲中文字幕| 欧美性猛交╳xxx乱大交人| 小说图片视频综合网站| 日韩一区二区三区影片| 久久久久久久久大av| 欧美变态另类bdsm刘玥| av卡一久久| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看| 色综合站精品国产| 黄色日韩在线| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 国产在线男女| 亚洲五月天丁香| 22中文网久久字幕| 亚洲美女搞黄在线观看| 久久久久久大精品| 一级黄色大片毛片| 国产精品一二三区在线看| 精品一区二区免费观看| 国产亚洲精品av在线| 国产成人影院久久av| 日本与韩国留学比较| 又爽又黄a免费视频| 中国国产av一级| 在线天堂最新版资源| 国产亚洲欧美98| 97超视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品欧美国产一区二区三| av国产免费在线观看| 国产av一区在线观看免费| 美女黄网站色视频| 联通29元200g的流量卡| 噜噜噜噜噜久久久久久91| 亚洲性久久影院| 国产精品久久久久久av不卡| 国产v大片淫在线免费观看| 熟妇人妻久久中文字幕3abv| 性插视频无遮挡在线免费观看| 身体一侧抽搐| 国产成人aa在线观看| 2022亚洲国产成人精品| 国产视频内射| 婷婷色综合大香蕉| 国产久久久一区二区三区| av又黄又爽大尺度在线免费看 | 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区| 99热这里只有是精品在线观看| 日本-黄色视频高清免费观看| 午夜精品在线福利| 午夜福利高清视频| 国产精品一二三区在线看| 3wmmmm亚洲av在线观看| 国产久久久一区二区三区| 亚洲18禁久久av| 日韩在线高清观看一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲国产精品合色在线| 一级毛片我不卡| 边亲边吃奶的免费视频| 日本撒尿小便嘘嘘汇集6| 黄片wwwwww| 乱码一卡2卡4卡精品| 国产黄a三级三级三级人| 97超碰精品成人国产| 免费人成视频x8x8入口观看| 久久欧美精品欧美久久欧美| 国产精品麻豆人妻色哟哟久久 | 精品久久国产蜜桃| 在线a可以看的网站| 免费看av在线观看网站| 九色成人免费人妻av| 在线观看66精品国产| 亚洲一区高清亚洲精品| 亚洲熟妇中文字幕五十中出| 久久久久久久久久久免费av| 日本撒尿小便嘘嘘汇集6| 亚洲熟妇中文字幕五十中出| 国产黄片美女视频| 秋霞在线观看毛片| 亚州av有码| 亚洲av成人av| 国产综合懂色| 亚洲在久久综合| 一个人观看的视频www高清免费观看| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影小说 |