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

    基于宏基因組學(xué)技術(shù)不同地區(qū)蛋源表面微生物的比較

    2021-06-04 02:17:26溫佳奇王秀娟宋雨奇代偉長(zhǎng)王玉華
    食品科學(xué) 2021年10期
    關(guān)鍵詞:吉林安徽群落

    李 俠,溫佳奇,王秀娟,宋雨奇,代偉長(zhǎng),王玉華

    (吉林農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院,吉林 春 130118)

    我國(guó)是生產(chǎn)和消費(fèi)禽蛋大國(guó),禽蛋和蛋制品的質(zhì)量與安全也備受消費(fèi)者關(guān)注[1]。蛋源在母禽生殖道及生產(chǎn)環(huán)境中會(huì)受到污染,表面帶有大量微生物,未經(jīng)處理的禽蛋在打蛋過(guò)程中微生物會(huì)感染到蛋液[2-4];液蛋制品加工時(shí),分離、過(guò)濾及均質(zhì)過(guò)程中蛋液也會(huì)因加工器械等未消毒完全而被污染。雖然灌裝后的液蛋制品會(huì)進(jìn)行滅菌,但前期大量微生物的污染會(huì)導(dǎo)致滅菌難度大,滅菌效果差,影響液蛋制品的質(zhì)量及保質(zhì)期。因此,微生物污染源的控制非常重要[5],在可能產(chǎn)生污染的過(guò)程里,蛋源的微生物污染數(shù)量最多,污染最嚴(yán)重,對(duì)蛋源微生物來(lái)源的控制成為控制液蛋制品微生物污染源重中之重。蛋源表面的微生物種類(lèi)繁多,常見(jiàn)的微生物有大腸桿菌屬、假單胞菌屬、不動(dòng)桿菌屬、變形桿菌屬、產(chǎn)氣單胞菌屬、葡萄球菌屬、沙門(mén)氏菌屬、產(chǎn)堿桿菌屬、副大腸桿菌屬等[6-8]。據(jù)相關(guān)材料報(bào)道,輸卵管帶菌的家禽所產(chǎn)的蛋,70%是帶菌蛋,所帶細(xì)菌如雞白痢沙門(mén)氏菌、雞副傷寒沙門(mén)氏菌及金黃色葡萄球菌等[9-10]。禽蛋腐敗的主要原因是微生物的存在,主要有非致病性細(xì)菌和霉菌。分解蛋白質(zhì)的微生物主要是芽孢桿菌屬、假單胞菌屬、變形桿菌屬、梭菌屬、葡萄球菌屬、無(wú)色桿菌屬、產(chǎn)堿桿菌屬、黃桿菌屬和腸桿菌屬等[11-12];分解糖的微生物有大腸桿菌屬、枯草桿菌屬和丁酸梭狀芽孢桿菌屬等[13];分解脂肪的微生物主要有假單胞菌屬、沙門(mén)氏菌屬等[14]。此外,還可能存在一些不可培養(yǎng)的微生物,特別是致病菌的存在,如霍亂弧菌和大腸埃希氏雖均具有生命活性,但不能培養(yǎng)[15-17],而這些未被培養(yǎng)的種類(lèi)卻代表了巨大的多樣性,此類(lèi)微生物可能存在潛在危害。因此,對(duì)蛋源表面微生物多樣性的探究對(duì)后續(xù)消毒處理方法的選擇具有重要意義,有助于液蛋制品產(chǎn)業(yè)的發(fā)展。

    雞蛋表面微生物種類(lèi)繁多,利用傳統(tǒng)的微生物培養(yǎng)法雖然能檢測(cè)出常見(jiàn)食品致病菌的數(shù)量,但并不能將所有菌相完全高效、準(zhǔn)確地檢測(cè)出來(lái),并且傳統(tǒng)方法耗時(shí)耗力,將高通量測(cè)序宏基因組技術(shù)應(yīng)用于蛋源表面微生物多樣性的檢測(cè),從宏觀角度準(zhǔn)確高效地檢測(cè)出微生物種類(lèi),為后期對(duì)蛋源前處理消毒方法的優(yōu)化研究奠定基礎(chǔ),同時(shí)也拓寬了宏基因組學(xué)技術(shù)在食品領(lǐng)域的探索。本實(shí)驗(yàn)采用高通量測(cè)序的宏基因組學(xué)技術(shù),分析吉林、安徽、北京和河北4 個(gè)不同地區(qū)蛋源表面微生物多樣性,確定蛋源表面微生物菌相分布,并分析不同蛋源之間的差異性,為蛋源消毒前處理方法的優(yōu)化提供一定的參考。

    1 材料與方法

    1.1 材料與試劑

    當(dāng)日產(chǎn)新鮮雞蛋,分別來(lái)自吉林、北京、河北、安徽地區(qū)的大規(guī)模蛋雞養(yǎng)殖企業(yè)(存欄量達(dá)300萬(wàn) 羽以上)。

    營(yíng)養(yǎng)瓊脂培養(yǎng)基、結(jié)晶紫中性紅膽鹽瓊脂培養(yǎng)基、緩沖蛋白胨水、亞硫酸鉍瓊脂培養(yǎng)基 北京奧博星生物技術(shù)有限公司;DNA試劑盒 北京全式金生物技術(shù)有限公司;其他試劑均為國(guó)產(chǎn)分析純。

    1.2 儀器與設(shè)備

    AUM220電子分析天平 日本島津公司;MCV-131超凈工作臺(tái) 上海博訊實(shí)業(yè)有限公司;AHWY-200B超低溫冰箱 日本三洋公司;YX280立式壓力蒸汽滅菌器、YX280立式壓力蒸汽滅菌器 上海申安醫(yī)療有限公司;Infinile M200多功能酶標(biāo)儀 瑞士TECAN公司;DHG-9146電熱恒溫培養(yǎng)箱 上海精宏實(shí)驗(yàn)設(shè)備有限公司。

    1.3 方法

    1.3.1 樣品的采集

    在同一月份,分別從吉林、北京、河北、安徽4 個(gè)地區(qū)的雞場(chǎng)進(jìn)行樣品采集,隨機(jī)采集每個(gè)雞場(chǎng)當(dāng)日產(chǎn)出的新鮮雞蛋,用滅菌生理鹽水棉拭子均勻擦拭蛋殼表面,每組樣品50 枚,每個(gè)地區(qū)采6 組平行樣品。取樣后,-80 ℃冷凍保存,用于后續(xù)蛋源蛋殼表面微生物多樣性分析。

    1.3.2 細(xì)菌菌落總數(shù)測(cè)定

    將擦拭蛋源表面的棉拭子充分溶于無(wú)菌生理鹽水、混勻,獲得樣品菌液,然后進(jìn)行稀釋,分別取稀釋倍數(shù)為10-3、10-4、10-5的樣品,按照GB 4789.2—2016《食品微生物菌落總數(shù)的測(cè)定》進(jìn)行蛋源表面菌落總數(shù)的測(cè)定。

    1.3.3 DNA的提取

    將擦拭蛋源表面的棉拭子充分溶于滅菌生理鹽水,采用EasyPureGenomic DNA試劑盒進(jìn)行DNA提取,嚴(yán)格按照試劑盒說(shuō)明書(shū)進(jìn)行操作。

    1.3.4 微生物多樣性分析

    通過(guò)蛋源表面微生物16S rRNA基因的高度可變區(qū)V3和V4進(jìn)行測(cè)序,基于測(cè)序結(jié)果進(jìn)行生物多樣性分析。實(shí)驗(yàn)流程為:提取的DNA→靶標(biāo)片段PCR擴(kuò)增→擴(kuò)增產(chǎn)物回收和純化→擴(kuò)增產(chǎn)物熒光定量→測(cè)序文庫(kù)制備→高通量測(cè)序,實(shí)驗(yàn)測(cè)序工作由上海派森諾生物科技股份有限公司完成。通過(guò)操作分類(lèi)單元(operational taxonomic unit,OTU)對(duì)獲得的序列進(jìn)行分類(lèi),并且通常使用97%的序列相似度作為OTU分區(qū)閾值,其等于分類(lèi)中的物種水平。根據(jù)不同樣品中OTU的豐度分布,評(píng)價(jià)每個(gè)樣品的多樣性水平,使用稀疏曲線反映測(cè)序深度是否達(dá)到標(biāo)準(zhǔn)。分析不同分類(lèi)水平的各地區(qū)微生物的具體組成,對(duì)吉林、安徽、北京和河北4 個(gè)地區(qū)的蛋源表面微生物菌群進(jìn)行α多樣性、群落結(jié)構(gòu)和β多樣性分析比較。

    1.4 數(shù)據(jù)分析

    實(shí)驗(yàn)數(shù)據(jù)利用Origin 2017、GraphPad Prism 7.0.軟件和Excel 2007對(duì)數(shù)據(jù)進(jìn)行整理、計(jì)算、統(tǒng)計(jì)與分析。P<0.05為差異顯著、P<0.01為差異極顯著、P<0.001為差異高度顯著,具有統(tǒng)計(jì)學(xué)意義。

    2 結(jié)果與分析

    2.1 不同地區(qū)蛋源表面微生物OTU分布

    Venn圖可用于統(tǒng)計(jì)不同樣本中共有和獨(dú)有的OTU數(shù)目,本研究選用相似度為97%的OTU,結(jié)果如圖1所示。吉林地區(qū)蛋源表面微生物OTU總的數(shù)目為2 341,其中獨(dú)有的OTU數(shù)目為1 285,河北地區(qū)蛋源表面微生物OTU總的數(shù)目為1 216,其中獨(dú)有的OTU數(shù)目為328,北京地區(qū)蛋源表面微生物OTU總的數(shù)目為1 708,其中獨(dú)有的OTU數(shù)目為295,安徽地區(qū)蛋源表面微生物OTU總的數(shù)目為1 301,其中獨(dú)有的OTU數(shù)目為135。4 個(gè)地區(qū)共有的OTU數(shù)目為218,吉林和河北地區(qū)共有的OTU數(shù)目為488,吉林和北京地區(qū)共有的OTU數(shù)目為765,吉林和安徽地區(qū)共有的OTU數(shù)目為574,河北和北京地區(qū)共有的OTU數(shù)目為675,河北和安徽地區(qū)共有的OTU數(shù)目為540,北京和安徽地區(qū)共有的OTU數(shù)目為957。由OTU數(shù)目分布Venn圖可以初步得出,吉林地區(qū)微生物豐富度最高,北京和安徽地區(qū)微生物最相似。

    圖1 不同地區(qū)蛋源表面微生物OTU的Venn圖Fig.1 Venn diagram showing unique and shared microbial OTUs among eggs from different regions

    2.2 不同地區(qū)蛋源表面微生物α多樣性分析

    α多樣性即每個(gè)樣本群落的多樣性,可用Shannon指數(shù)曲線和α多樣性指數(shù)等進(jìn)行評(píng)價(jià)。

    2.2.1 不同地區(qū)蛋源表面微生物Shannon指數(shù)曲線

    Shannon指數(shù)曲線是反映樣品中微生物多樣性的指數(shù),利用各樣品的測(cè)序量在不同測(cè)序深度時(shí)的微生物多樣性指數(shù)構(gòu)建曲線,以此反映各樣本在不同測(cè)序數(shù)量時(shí)的微生物多樣性[18]。當(dāng)曲線趨向平坦時(shí),說(shuō)明測(cè)序數(shù)據(jù)量足夠大,可以反映樣品中絕大多數(shù)的微生物信息。可以用來(lái)判斷測(cè)序的有效性[19]。本實(shí)驗(yàn)選用Shannon指數(shù)曲線對(duì)測(cè)序有效性進(jìn)行分析。

    圖2 不同地區(qū)蛋源表面微生物Shannon指數(shù)曲線Fig.2 Shannon index curves showing microbial community diversity on eggs from different regions

    如圖2所示,橫坐標(biāo)表示每個(gè)樣品中隨機(jī)抽取序列的總數(shù);縱坐標(biāo)表示微生物多樣性的Shannon指數(shù)。曲線的長(zhǎng)度反映了樣品測(cè)序的數(shù)量。曲線的漸變程度反映了測(cè)序深度對(duì)觀察樣品多樣性的影響,曲線趨于平坦,表明測(cè)序已趨于飽和,增加測(cè)序數(shù)據(jù)無(wú)法再找到更多的OTU。隨著樣品測(cè)序數(shù)量的不斷增加,當(dāng)其大于5 000時(shí),曲線已趨于平緩,說(shuō)明測(cè)序數(shù)據(jù)量足夠大,足夠反應(yīng)樣本中包含的微生物多樣性。

    2.2.2 不同地區(qū)蛋源表面微生物α多樣性指數(shù)

    對(duì)于微生物群落的α多樣性可以通過(guò)Chao1指數(shù)、ACE指數(shù)、Shannon指數(shù)和Simpson指數(shù)反映。一般而言,Chao1指數(shù)或ACE指數(shù)越大,表明群落的物種種類(lèi)越多、豐富度越高[20]。而Shannon多樣性指數(shù)不僅反映群落的豐富度,還綜合考慮了群落的均勻度。Shannon指數(shù)值越高,表明群落的多樣性越高;Simpson指數(shù)與Shannon指數(shù)類(lèi)似,是群落多樣性的評(píng)價(jià)常用指數(shù)之一,Simpson指數(shù)值越高,說(shuō)明群落的多樣性越高[21]。

    表1 不同地區(qū)蛋源表面菌群微生物多樣性指數(shù)表Table 1 Microbial diversity indexes of surface microbial flora in eggs from different regions

    由表1可知,4 個(gè)地區(qū)都具有較高的菌群多樣性。其中Chao1指數(shù)和ACE指數(shù)表現(xiàn)為JL>BJ>AH>HE,說(shuō)明在菌群豐富度上按照吉林、北京、安徽、河北的順序依次降低;而Shannon指數(shù)和Simpson指數(shù)表現(xiàn)為JL>BJ>HE>AH,表明在群落多樣性上按照吉林、北京、河北、安徽的順序依次降低。由此可知,4 個(gè)地區(qū)相比,吉林地區(qū)蛋源表面微生物的群落豐富度及群落多樣性均最高。

    2.3 不同地區(qū)蛋源表面細(xì)菌群落結(jié)構(gòu)分析

    2.3.1 不同地區(qū)蛋源表面基于門(mén)水平的微生物分析

    根據(jù)OTU劃分和分類(lèi)地位鑒定,獲得了吉林、安徽、北京和河北在門(mén)、綱、目、科和屬分類(lèi)水平的具體微生物分布狀況。在本研究中,選擇門(mén)和屬這2 個(gè)水平,對(duì)4 個(gè)地區(qū)的微生物進(jìn)行分析。4 個(gè)地區(qū)蛋源基于門(mén)水平上的群落組成如圖3所示。

    吉林地區(qū)的蛋源樣品細(xì)菌組成為放線菌門(mén)(Actinobacteria)、厚壁菌門(mén)(Firmicutes)、變形菌門(mén)(Proteobacteria)和擬桿菌門(mén)(Bacteroidetes),相對(duì)豐度分別為49.12%、31.93%、8.41%及4.37%,還包括少量的異常球菌-棲熱菌門(mén)(Deinococcus-Thermus)、藍(lán)藻門(mén)(Cyanobacteria)以及綠彎菌門(mén)(Chloroflex);安徽地區(qū)的蛋源樣品細(xì)菌組成為Proteobacteria、Firmicutes和Actinobacteria,相對(duì)豐度分別為77.03%、18.21%及2.90%,還包括少量的Cyanobacteria;北京地區(qū)的蛋源樣品細(xì)菌組成為Proteobacteria、Firmicutes和Actinobacteria,相對(duì)豐度分別為57.93%、37.68%及2.93%,還包括少量的Cyanobacteria;河北地區(qū)的蛋源樣品細(xì)菌組成為Proteobacteria、Firmicutes和Actinobacteria,相對(duì)豐度分別為84.28%、10.42%及3.12%,還包括少量的Cyanobacteria。安徽、北京和河北3 個(gè)地區(qū)蛋源菌群組成相似,優(yōu)勢(shì)菌門(mén)均為Proteobacteria,而在吉林地區(qū)Proteobacteria豐度僅為8.40%,變形菌門(mén)中包括很多種食品中常見(jiàn)的致病菌,如幽門(mén)螺桿菌、沙門(mén)氏菌和大腸桿菌等,吉林地區(qū)蛋源表面此類(lèi)致病菌可能遠(yuǎn)低于其他3 個(gè)地區(qū);而吉林地區(qū)的優(yōu)勢(shì)菌門(mén)Actinobacteria在其他3 個(gè)地區(qū)豐度較低, Chloroflex僅在吉林地區(qū)發(fā)現(xiàn)。因此,吉林地區(qū)與其他3 個(gè)地區(qū)在蛋源表面微生物菌群有很大差異,但由于其變形菌門(mén)豐度較低,含有常見(jiàn)致病菌的風(fēng)險(xiǎn)可能低于其他3 個(gè)地區(qū),因此更適合作為蛋源生產(chǎn)產(chǎn)區(qū)。

    圖3 基于門(mén)分類(lèi)水平上的不同地區(qū)樣品菌群分布圖Fig.3 Distribution of microflora on eggs from different regions at the phylum level

    2.3.2 不同地區(qū)蛋源表面基于屬水平的微生物分析

    圖4 基于屬分類(lèi)水平上的不同地區(qū)樣品菌群分布圖Fig.4 Distribution of microflora on eggs from different regions at the genus level

    如圖4所示,吉林地區(qū)的蛋源樣品中豐度較高的菌屬為考克氏菌屬(Kocuria)占10.1%,短狀桿菌屬(Brachybacterium)占9.83%,榛狀桿菌屬(Corynebacterium1)占5.2%,uncultured占7.2%,微球菌屬(Micrococcus)占3.95%,羅思氏菌屬(Rothia)占3.80%,嗜冷桿菌屬(Psychrobacter)占0.96%,不動(dòng)桿菌屬(Acinetobacter)占0.93%;安徽地區(qū)蛋源樣品中豐度較高的菌屬為Acinetobacter占75.52%,氣球菌屬(Aerococcus)占8.40%,乳酸菌屬(Lactobacillus)占2.82%,Bromus tectorum占1.55%,葡萄球菌屬(Staphylococcus)占1.20%;北京地區(qū)蛋源樣品中豐度較高的菌屬為Acinetobacter占26.42%,Psychrobacter占20.44%,Lactobacillus占15.12%,Aerococcus占14.36%,Staphylococcus占3.53%,uncultured占2.11%,埃希氏菌志賀氏菌屬(Escherichia-Shigella)占3.15%;河北地區(qū)蛋源樣品中豐度較高的菌屬為Acinetobacter占26.70%,Psychrobacter占13.69%,Lactobacillus占7.32%,假單胞菌屬(Pseudomonas)占16.28%,Massillia占15.39%。安徽、北京和河北3 個(gè)地區(qū)菌屬分布較為相似,優(yōu)勢(shì)菌屬均為Acinetobacter,而吉林地區(qū)Acinetobacter豐度僅為0.89%,Brachybacterium為吉林地區(qū)獨(dú)有的豐度較高的菌屬。由此表明安徽、北京和河北3 個(gè)地區(qū)菌屬分布較為相似,優(yōu)勢(shì)菌屬均為條件致病菌的Acinetobacter,而且存在其他致病菌屬,即安徽地區(qū)蛋源樣品Staphylococcus,北京地區(qū)蛋源樣品Staphylococcus和Escherichia-Shigella,河北地區(qū)蛋源樣品Pseudomonas。可見(jiàn)不同地區(qū)蛋源表面微生物屬水平都存在一定差異性,相對(duì)豐度差異顯著,北方地區(qū)與中部以南地區(qū)差異極顯著(P<0.01),北方地區(qū)更適宜禽蛋生產(chǎn)。同時(shí),為保障液蛋制品的安全性,蛋制品企業(yè)應(yīng)該根據(jù)蛋源微生物的種類(lèi)適當(dāng)調(diào)整消毒處理方法。

    2.3.3 不同地區(qū)蛋源表面微生物多樣性聚類(lèi)熱圖分析

    根據(jù)不同樣本的豐度分布或相似程度所組成的聚類(lèi)分析圖。其目的是通過(guò)聚類(lèi),可以區(qū)分高豐度和低豐度的分類(lèi)單位,并且由顏色組成梯度反映樣本之間的群落組成相似性。

    如圖5所示,4 個(gè)地區(qū)蛋源表面豐度前50 位的菌屬聚類(lèi)熱圖,紅色代表相應(yīng)樣品中豐度較高的屬,綠色代表豐度較低的屬。在吉林地區(qū)豐度較高的29 個(gè)屬在其他地區(qū)豐度較低;而河北地區(qū)豐度較高的屬為Brevundimonas、Pseudomonas、Massillia、Alkanindiges、Cupriavidus、Sphingomonas和Ochrobactrum;安徽和北京地區(qū)蛋源屬的分布最為相似,豐度較高的屬為Acinetobacter、Jeoigalicoccus、Fackiamia、Lysinbacliius、Psychrobacter。

    圖5 不同地區(qū)蛋源表面微生物多樣性聚類(lèi)分析熱圖Fig.5 Heatmap of clustering analysis for microbial community diversity on eggs from different regions

    2.4 不同地區(qū)蛋源表面微生物β多樣性分析

    β多樣性分析的主要目的是評(píng)價(jià)不同樣本之間群落結(jié)構(gòu)的相似性。主要通過(guò)主成分分析(principal component analysis,PCA)、多維尺度分析(multidime-nsional scaling,MDS),對(duì)群落數(shù)據(jù)結(jié)構(gòu)進(jìn)行自然分解并通過(guò)對(duì)樣本排序,從而觀察樣本之間的差異。

    2.4.1 不同地區(qū)蛋源表面菌群PCA

    菌群的PCA主要是通過(guò)線性變換將菌群OTU豐度矩陣組合,并通過(guò)線性變換將其投影到較低維空間坐標(biāo)系中。從而達(dá)到降維的目的,簡(jiǎn)化數(shù)據(jù)結(jié)構(gòu),并顯示樣本的自然分布[22]。PCA可以從原始數(shù)據(jù)中提取樣本之間最重要的差異,并根據(jù)這些差異在新的低維坐標(biāo)系中對(duì)樣本進(jìn)行排序。使樣本接近新坐標(biāo)系可以最大程度地恢復(fù)樣本之間的實(shí)際差異[23]。因此,通過(guò)分析PCA可以知道群落樣本的主要分布特征,從而量化樣本之間的差異和相似性。

    圖6 不同地區(qū)蛋源表面菌群PCA二維排序圖Fig.6 Two-dimensional sequencing diagram of PCA analysis for discrimination of microbial community composition on eggs from different regions

    如圖6所示,橫坐標(biāo)為PC1,其貢獻(xiàn)率為74.26%??v坐標(biāo)為PC2,貢獻(xiàn)率為13.65%。這2 個(gè)PC是樣品微生物群落結(jié)構(gòu)組成差異的主要因子。在PC1方向上,吉林和河北地區(qū)蛋源樣品點(diǎn)與安徽地區(qū)樣品點(diǎn)分布差異顯著(P<0.05),樣本點(diǎn)處于PC1水平的2 個(gè)極端;在PC2水平上,吉林地區(qū)樣品點(diǎn)與河北地區(qū)樣本點(diǎn)分布差異顯著(P<0.05),吉林地區(qū)樣本點(diǎn)分布于Y軸的下方,河北地區(qū)樣本點(diǎn)分布于Y軸的上方。由此可知,吉林地區(qū)蛋源與河北地區(qū)蛋源菌群結(jié)構(gòu)及微生物多樣性差異極顯著(P<0.01)。

    2.4.2 基于UniFrac距離的NMDS非度量多維尺度分析

    NMDS分析是基于樣本距離矩陣的MDS分析方法。因此,在特定的距離尺度下描述了樣本的分布特征。NMDS分析對(duì)樣本距離進(jìn)行排序,以使低維空間中的樣本的順序盡可能彼此接近。因此,NMDS分析不受樣本距離值的影響,僅考慮彼此之間的大小關(guān)系[24]。

    圖7 不同地區(qū)蛋源表面菌群NMDS分析的二維排序圖Fig.7 Two-dimensional sequencing diagram of NMDS analysis for discrimination of microbial community composition on eggs from different regions

    由圖7A可知,僅考慮OTU在樣本中存在與否的情況下,安徽地區(qū)和北京地區(qū)樣本之間的微生物群落結(jié)構(gòu)相似度高,而吉林地區(qū)與河北地區(qū)樣本差異性大,這與PCA中得到的結(jié)論吻合。而考慮由于群落豐度梯度導(dǎo)致差異的情況,如圖7B可以看出,吉林地區(qū)樣本與其他3 個(gè)樣本之間有很大差異。因此,無(wú)論是否考慮OTU的豐度,吉林地區(qū)樣本與其他地區(qū)差異都較大。

    3 討 論

    蛋源微生物污染一直是蛋品領(lǐng)域值得研究的問(wèn)題,采用高通量測(cè)序的宏基因組學(xué)技術(shù)能夠一次并行對(duì)幾十萬(wàn)到幾百萬(wàn)條DNA分子進(jìn)行序列測(cè)定,高效準(zhǔn)確地測(cè)定出蛋源表面所有已知微生物菌群分布,它在分析微生物的群落結(jié)構(gòu)時(shí)有著獨(dú)特的優(yōu)勢(shì),能夠通過(guò)從環(huán)境樣本中直接獲取總DNA進(jìn)行文庫(kù)構(gòu)建并測(cè)序[25-27],用16S rRNA基因的測(cè)序數(shù)據(jù)估計(jì)微生物群落的物種構(gòu)成,更加真實(shí)地揭示原位環(huán)境中微生物群落的復(fù)雜性和多樣性[28-29],對(duì)蛋源微生物有一個(gè)宏觀的了解,可以詳細(xì)地掌握蛋源表面所有基因庫(kù)里已知菌屬的情況;同時(shí),本研究也發(fā)現(xiàn)蛋源表面微生物除了已知的菌屬外,還存在一定量未知菌屬,為今后更加深入研究蛋源表面微生物種類(lèi)奠定良好基礎(chǔ),也為后續(xù)蛋源表面微生物多樣性研究提供理論依據(jù)。本研究選擇吉林、北京、河北、安徽4 個(gè)地區(qū)作為采樣地區(qū),是因?yàn)檫@4 個(gè)地區(qū)為國(guó)內(nèi)主要蛋源產(chǎn)區(qū),安徽以南溫度過(guò)高,不適合大規(guī)模蛋雞養(yǎng)殖,因此沒(méi)有代表性。研究目的在于探究北部、中部和南部蛋源在菌群多樣性上的差異性,從結(jié)果可以看出吉林地區(qū)的微生物多樣性和豐富度與其他3 個(gè)地區(qū)差異性最大,優(yōu)勢(shì)菌門(mén)也不同,但值得關(guān)注的是,擁有較多食品致病菌的變形菌門(mén)在吉林地區(qū)豐度顯著低于其他3 個(gè)地區(qū),蛋液感染此類(lèi)病源菌的風(fēng)險(xiǎn)也隨之減少;而且不同地區(qū)樣品微生物屬水平差異更大,吉林地區(qū)的蛋源樣品中豐度較高的菌屬為Kocuria和Brachybacterium分別為10.1%和9.83%,目前還沒(méi)有關(guān)于這2 個(gè)屬存在致病菌的報(bào)道;安徽、北京和河北3 個(gè)地區(qū)菌屬分布較為相似,優(yōu)勢(shì)菌屬均為條件致病菌的Acinetobacter,另外還存在致病菌屬。Acinetobacter是條件致病菌,當(dāng)機(jī)體抵抗力降低時(shí)易引起機(jī)體感染,是引起醫(yī)院內(nèi)感染的重要機(jī)會(huì)致病菌之一??梢鸷粑栏腥尽⊙Y、腦膜炎、心內(nèi)膜炎、傷口及皮膚感染、泌尿生殖道感染等。重癥者可導(dǎo)致死亡。安徽地區(qū)蛋源樣品致病菌屬為Staphylococcus,北京地區(qū)蛋源樣品致病菌屬為Staphylococcus和Escherichia-Shigella,河北地區(qū)蛋源樣品致病菌屬為Pseudomonas。從相關(guān)文獻(xiàn)報(bào)道可知,目前國(guó)內(nèi)外對(duì)蛋源表面微生物的研究多采用傳統(tǒng)的微生物培養(yǎng)法,集中在細(xì)菌的數(shù)量檢測(cè)上,如Moyle等[30]對(duì)自由放養(yǎng)的2 個(gè)農(nóng)場(chǎng)的雞蛋表面細(xì)菌測(cè)定結(jié)果分別為3.84(lg(CFU/枚))和3.69(lg(CFU/枚)),分離培養(yǎng)得到腸桿菌科細(xì)菌平均水平為0.98(lg(CFU/枚))和0.68(lg(CFU/枚));薛艷蓉等[31]報(bào)道蛋殼表面細(xì)菌數(shù)一般在103CFU/枚。傳統(tǒng)微生物研究采用實(shí)驗(yàn)室分離培養(yǎng),只有少數(shù)的微生物能夠被成功培養(yǎng)鑒定出來(lái),而且用時(shí)長(zhǎng)、篩選菌屬也不夠廣泛,不能進(jìn)一步分析微生物群落組成及其多樣性。宏基因組學(xué)克服了這一限制,以基因組技術(shù)為基礎(chǔ)通過(guò)對(duì)環(huán)境中全部DNA系統(tǒng)全面的研究,從宏觀角度準(zhǔn)確高效地對(duì)微生物種類(lèi)進(jìn)行測(cè)定。本研究中選取安徽、北京、河北、吉林地區(qū)作為取樣點(diǎn),涉及范圍廣泛,利用宏基因組學(xué)技術(shù)經(jīng)過(guò)對(duì)蛋源表面菌群多樣性分析,結(jié)果表明以吉林地區(qū)為代表的北方,更適合作為蛋源的生產(chǎn)區(qū),蛋源表面微生物的菌群多樣性分析結(jié)果為后續(xù)消毒處理方法的優(yōu)化奠定了一定的理論基礎(chǔ)。

    4 結(jié) 論

    采用高通量測(cè)序的宏基因組學(xué)技術(shù)對(duì)4 個(gè)地區(qū)蛋源表面菌群多樣性系統(tǒng)分析,不同地區(qū)蛋源表面主要微生物門(mén)水平和屬水平都存在一定差異性,相對(duì)豐度差異更顯著,北方地區(qū)與中部以南地區(qū)差異明顯,北方地區(qū)更適宜大規(guī)模蛋雞養(yǎng)殖。為保障液蛋制品的安全性,蛋制品企業(yè)應(yīng)該根據(jù)蛋源表面微生物組成情況適當(dāng)調(diào)整消毒處理方法,對(duì)液蛋制品加工行業(yè)有著重要意義。

    猜你喜歡
    吉林安徽群落
    13.吉林卷
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    吉林卷
    吉林卷
    安徽醫(yī)改自我完善主動(dòng)糾錯(cuò)
    安徽藥采如何“三步走”
    安徽 諸多方面走在前列
    安徽為什么選擇帶量采購(gòu)
    春季和夏季巢湖浮游生物群落組成及其動(dòng)態(tài)分析
    免费在线观看影片大全网站| 欧美bdsm另类| 国产欧美日韩一区二区精品| 天美传媒精品一区二区| 亚洲av.av天堂| 91av网一区二区| 嫩草影视91久久| 日本黄大片高清| 午夜两性在线视频| 午夜免费男女啪啪视频观看 | 久久久久精品国产欧美久久久| 成人av在线播放网站| 日本a在线网址| 久久久久久久午夜电影| 色在线成人网| 国产极品精品免费视频能看的| 在线观看舔阴道视频| 成人一区二区视频在线观看| 久久国产精品人妻蜜桃| 国产v大片淫在线免费观看| 夜夜躁狠狠躁天天躁| 久久精品国产99精品国产亚洲性色| 成人av在线播放网站| 给我免费播放毛片高清在线观看| 一个人观看的视频www高清免费观看| 欧美日韩乱码在线| 久久久精品欧美日韩精品| 精品午夜福利在线看| 成人特级黄色片久久久久久久| 国产三级中文精品| 麻豆成人午夜福利视频| 亚洲av成人av| 美女 人体艺术 gogo| 成年女人永久免费观看视频| АⅤ资源中文在线天堂| 美女高潮喷水抽搐中文字幕| 久久99热这里只有精品18| 亚洲精品色激情综合| 色综合站精品国产| 97超视频在线观看视频| 色综合婷婷激情| 国产乱人伦免费视频| 国内少妇人妻偷人精品xxx网站| 波多野结衣高清无吗| 成人三级黄色视频| 亚洲自偷自拍三级| 国产一区二区亚洲精品在线观看| a在线观看视频网站| 精品一区二区三区视频在线| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲av天美| www.熟女人妻精品国产| 人人妻,人人澡人人爽秒播| 欧美最新免费一区二区三区 | 亚洲精品一卡2卡三卡4卡5卡| 在线观看免费视频日本深夜| 久久精品91蜜桃| 国产一级毛片七仙女欲春2| 悠悠久久av| 日本黄大片高清| 成人av一区二区三区在线看| 国产爱豆传媒在线观看| 日韩免费av在线播放| ponron亚洲| 精品一区二区三区人妻视频| 在线天堂最新版资源| 一区二区三区激情视频| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 永久网站在线| 男人舔奶头视频| 国产精品久久久久久亚洲av鲁大| 神马国产精品三级电影在线观看| av天堂在线播放| 国产午夜精品久久久久久一区二区三区 | 99久久成人亚洲精品观看| 内地一区二区视频在线| 国产成人啪精品午夜网站| 色av中文字幕| 高清在线国产一区| 天堂影院成人在线观看| 国产 一区 欧美 日韩| 久久久精品大字幕| 亚洲成av人片在线播放无| 日韩国内少妇激情av| 少妇高潮的动态图| 国产精品综合久久久久久久免费| 麻豆久久精品国产亚洲av| 国产乱人伦免费视频| 日本黄大片高清| 国产爱豆传媒在线观看| 91字幕亚洲| av在线观看视频网站免费| 97人妻精品一区二区三区麻豆| 一级a爱片免费观看的视频| 亚洲欧美日韩无卡精品| 国产野战对白在线观看| 亚洲精品一区av在线观看| 中文字幕人成人乱码亚洲影| 久久精品国产自在天天线| 国产精品av视频在线免费观看| 欧美bdsm另类| 日本 av在线| 国产精品久久久久久精品电影| 91av网一区二区| 欧美午夜高清在线| 亚洲av成人av| 免费高清视频大片| 国产国拍精品亚洲av在线观看| 亚洲自偷自拍三级| 男女床上黄色一级片免费看| 精品久久国产蜜桃| 成人无遮挡网站| 欧美日韩福利视频一区二区| 欧美色欧美亚洲另类二区| 欧美日韩福利视频一区二区| 欧美日韩福利视频一区二区| 男人和女人高潮做爰伦理| 给我免费播放毛片高清在线观看| 最近在线观看免费完整版| x7x7x7水蜜桃| 老司机福利观看| 日本三级黄在线观看| 久久伊人香网站| 国产精品一区二区免费欧美| 欧美最黄视频在线播放免费| a级一级毛片免费在线观看| 亚洲在线自拍视频| 真人做人爱边吃奶动态| 搞女人的毛片| 成年版毛片免费区| 色哟哟哟哟哟哟| 999久久久精品免费观看国产| 我的老师免费观看完整版| 日韩精品青青久久久久久| 亚洲美女黄片视频| 久久精品久久久久久噜噜老黄 | 久久久国产成人精品二区| 91麻豆av在线| 最近最新中文字幕大全电影3| 可以在线观看毛片的网站| 真人做人爱边吃奶动态| 国产伦一二天堂av在线观看| 黄色日韩在线| 99久久久亚洲精品蜜臀av| av天堂在线播放| 热99在线观看视频| 精品久久久久久久久av| 日韩欧美免费精品| 中亚洲国语对白在线视频| 99精品久久久久人妻精品| 女人被狂操c到高潮| 三级毛片av免费| 亚洲不卡免费看| 天堂av国产一区二区熟女人妻| 欧美激情国产日韩精品一区| 午夜两性在线视频| 人人妻,人人澡人人爽秒播| 一区福利在线观看| www.999成人在线观看| 少妇裸体淫交视频免费看高清| 99久久无色码亚洲精品果冻| 成年人黄色毛片网站| 极品教师在线视频| 国产精品野战在线观看| 在线观看66精品国产| 神马国产精品三级电影在线观看| 在线免费观看的www视频| 日日干狠狠操夜夜爽| 99riav亚洲国产免费| aaaaa片日本免费| 男女那种视频在线观看| 一本综合久久免费| 国产高潮美女av| 亚洲精品一区av在线观看| 精品不卡国产一区二区三区| 1000部很黄的大片| av天堂中文字幕网| 欧美日本亚洲视频在线播放| 免费电影在线观看免费观看| 欧美性感艳星| 久久九九热精品免费| 成人亚洲精品av一区二区| 丁香六月欧美| 全区人妻精品视频| 国内揄拍国产精品人妻在线| 偷拍熟女少妇极品色| 亚洲av免费高清在线观看| 内地一区二区视频在线| a级一级毛片免费在线观看| 精品不卡国产一区二区三区| www日本黄色视频网| 国产蜜桃级精品一区二区三区| 国产成人福利小说| 日本一二三区视频观看| 国产单亲对白刺激| 高潮久久久久久久久久久不卡| 成人国产综合亚洲| 欧美+亚洲+日韩+国产| 丰满的人妻完整版| 国产精品女同一区二区软件 | 午夜精品在线福利| 内地一区二区视频在线| 国产熟女xx| 99国产极品粉嫩在线观看| 久久久久久久久久黄片| 亚洲七黄色美女视频| 亚洲成人久久爱视频| www.熟女人妻精品国产| 成人永久免费在线观看视频| 免费在线观看影片大全网站| 国产真实伦视频高清在线观看 | 99热这里只有精品一区| 日本三级黄在线观看| 小蜜桃在线观看免费完整版高清| 中文资源天堂在线| 成年版毛片免费区| 天堂影院成人在线观看| 成人高潮视频无遮挡免费网站| 在线a可以看的网站| 成人毛片a级毛片在线播放| 伊人久久精品亚洲午夜| 国产视频一区二区在线看| 97碰自拍视频| 一a级毛片在线观看| 久久人人爽人人爽人人片va | or卡值多少钱| 欧美黑人巨大hd| 最近在线观看免费完整版| 国产精品久久久久久久久免 | 天堂av国产一区二区熟女人妻| netflix在线观看网站| 国产精品99久久久久久久久| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线观看免费| 人人妻人人澡欧美一区二区| 最新中文字幕久久久久| 国产欧美日韩一区二区三| 大型黄色视频在线免费观看| 少妇丰满av| 久久人人爽人人爽人人片va | 天天躁日日操中文字幕| 性欧美人与动物交配| 欧美日本视频| 人妻夜夜爽99麻豆av| 亚洲av美国av| 精品欧美国产一区二区三| 精品一区二区三区av网在线观看| 又紧又爽又黄一区二区| 久久久成人免费电影| 欧美xxxx黑人xx丫x性爽| www.熟女人妻精品国产| 国产精品乱码一区二三区的特点| 久久人人精品亚洲av| 亚洲无线在线观看| 久久精品国产99精品国产亚洲性色| 亚洲国产精品合色在线| 在线观看一区二区三区| 久久6这里有精品| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼| 人人妻,人人澡人人爽秒播| 嫩草影视91久久| 国产欧美日韩一区二区三| 97热精品久久久久久| 成人亚洲精品av一区二区| 又黄又爽又免费观看的视频| 亚洲乱码一区二区免费版| 天堂√8在线中文| 在现免费观看毛片| 欧美潮喷喷水| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产亚洲av香蕉五月| 69av精品久久久久久| 特大巨黑吊av在线直播| 久久久久久久久久成人| 日本与韩国留学比较| 日韩av在线大香蕉| 亚洲av电影不卡..在线观看| 日韩欧美在线乱码| 国产乱人视频| 丰满乱子伦码专区| 一进一出抽搐动态| 免费一级毛片在线播放高清视频| 在线观看av片永久免费下载| 色哟哟·www| 少妇人妻精品综合一区二区 | 日韩欧美精品v在线| 九九热线精品视视频播放| 久9热在线精品视频| 久久国产精品人妻蜜桃| 国产一级毛片七仙女欲春2| 日本精品一区二区三区蜜桃| 欧美在线一区亚洲| 免费黄网站久久成人精品 | 露出奶头的视频| 日韩高清综合在线| 亚洲天堂国产精品一区在线| 精品不卡国产一区二区三区| 五月伊人婷婷丁香| 国产高清有码在线观看视频| 亚洲人成伊人成综合网2020| 波多野结衣巨乳人妻| 波野结衣二区三区在线| 少妇被粗大猛烈的视频| 久久国产乱子伦精品免费另类| 亚洲av日韩精品久久久久久密| 午夜免费激情av| 观看免费一级毛片| 99久久久亚洲精品蜜臀av| 久久久精品欧美日韩精品| 国产成人欧美在线观看| 亚洲经典国产精华液单 | 99热6这里只有精品| 在线十欧美十亚洲十日本专区| 久久这里只有精品中国| 亚洲精品456在线播放app | 国产主播在线观看一区二区| 国产美女午夜福利| 日韩免费av在线播放| 免费大片18禁| 欧美高清成人免费视频www| 九九久久精品国产亚洲av麻豆| 久久草成人影院| 免费观看精品视频网站| 精品一区二区免费观看| 欧美丝袜亚洲另类 | 国产精品一区二区三区四区免费观看 | 午夜激情欧美在线| 婷婷精品国产亚洲av| 国产一区二区激情短视频| 成人美女网站在线观看视频| 最新中文字幕久久久久| 亚洲久久久久久中文字幕| 美女大奶头视频| 美女cb高潮喷水在线观看| 国产精品98久久久久久宅男小说| 一级毛片久久久久久久久女| 在线播放无遮挡| 美女高潮喷水抽搐中文字幕| 又爽又黄a免费视频| 日韩欧美精品v在线| 免费在线观看亚洲国产| 国产精品影院久久| 国产日本99.免费观看| 日韩精品青青久久久久久| 99精品久久久久人妻精品| 精品一区二区三区人妻视频| 久久久国产成人免费| 熟女电影av网| 亚洲最大成人手机在线| 亚洲精品在线美女| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 一个人观看的视频www高清免费观看| 757午夜福利合集在线观看| 九色成人免费人妻av| 99热精品在线国产| 国产一区二区在线av高清观看| 中文字幕久久专区| 床上黄色一级片| 别揉我奶头 嗯啊视频| 最新中文字幕久久久久| 日韩欧美在线二视频| 日韩欧美精品v在线| 国产精品亚洲av一区麻豆| 成人特级黄色片久久久久久久| a在线观看视频网站| 91狼人影院| 俄罗斯特黄特色一大片| 日日摸夜夜添夜夜添av毛片 | 天天躁日日操中文字幕| 国产成+人综合+亚洲专区| 欧美一区二区亚洲| 午夜免费男女啪啪视频观看 | 久久热精品热| 午夜激情福利司机影院| 动漫黄色视频在线观看| 精品熟女少妇八av免费久了| 日韩欧美国产一区二区入口| av女优亚洲男人天堂| 我的老师免费观看完整版| av女优亚洲男人天堂| 变态另类成人亚洲欧美熟女| 国产亚洲av嫩草精品影院| 在线观看舔阴道视频| 中文字幕熟女人妻在线| 日韩国内少妇激情av| 成年女人永久免费观看视频| 中文资源天堂在线| 性色av乱码一区二区三区2| 成人美女网站在线观看视频| 久久精品综合一区二区三区| 午夜精品一区二区三区免费看| 长腿黑丝高跟| av专区在线播放| 免费电影在线观看免费观看| 精品国内亚洲2022精品成人| 国产成年人精品一区二区| 欧美国产日韩亚洲一区| 最新中文字幕久久久久| 国产大屁股一区二区在线视频| 亚洲av美国av| 看片在线看免费视频| 色吧在线观看| 夜夜躁狠狠躁天天躁| 床上黄色一级片| 中文字幕av在线有码专区| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 别揉我奶头 嗯啊视频| 国产主播在线观看一区二区| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 日韩亚洲欧美综合| 精品一区二区三区视频在线观看免费| 国产一区二区三区在线臀色熟女| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图 | 偷拍熟女少妇极品色| 首页视频小说图片口味搜索| 国产伦精品一区二区三区视频9| 午夜久久久久精精品| av国产免费在线观看| 亚洲国产精品久久男人天堂| 97人妻精品一区二区三区麻豆| 99riav亚洲国产免费| 91在线观看av| 亚洲最大成人av| 欧美最黄视频在线播放免费| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 观看美女的网站| 久久久久免费精品人妻一区二区| 亚洲性夜色夜夜综合| 少妇熟女aⅴ在线视频| 深爱激情五月婷婷| 波野结衣二区三区在线| 欧美xxxx黑人xx丫x性爽| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 成年女人看的毛片在线观看| 日韩 亚洲 欧美在线| 国产免费一级a男人的天堂| 国产中年淑女户外野战色| 麻豆久久精品国产亚洲av| 成年人黄色毛片网站| 久久精品人妻少妇| 欧美xxxx性猛交bbbb| 国产真实伦视频高清在线观看 | 国内久久婷婷六月综合欲色啪| 欧美一区二区精品小视频在线| 精品乱码久久久久久99久播| 中文字幕熟女人妻在线| 美女高潮的动态| 日日摸夜夜添夜夜添av毛片 | 夜夜看夜夜爽夜夜摸| 免费人成视频x8x8入口观看| 麻豆一二三区av精品| 久久久久久久午夜电影| 天堂网av新在线| 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 免费无遮挡裸体视频| 美女xxoo啪啪120秒动态图 | 久久精品国产亚洲av香蕉五月| 91麻豆av在线| 久久国产乱子伦精品免费另类| 久久久久性生活片| 欧美+亚洲+日韩+国产| 国产真实乱freesex| 亚洲第一区二区三区不卡| 国产成人影院久久av| av福利片在线观看| 色播亚洲综合网| 国产毛片a区久久久久| 偷拍熟女少妇极品色| 亚洲人成网站高清观看| 亚洲精品色激情综合| 中文在线观看免费www的网站| 国产av在哪里看| 在线免费观看不下载黄p国产 | 国产高清激情床上av| 久久久久久久亚洲中文字幕 | 国产成人啪精品午夜网站| 国产精品一区二区三区四区久久| 淫妇啪啪啪对白视频| 亚洲狠狠婷婷综合久久图片| 国产精品98久久久久久宅男小说| 亚洲人成电影免费在线| 中文字幕av成人在线电影| 精品日产1卡2卡| 亚洲精品色激情综合| 色综合欧美亚洲国产小说| 久久久久久久久久黄片| 色5月婷婷丁香| 国模一区二区三区四区视频| 中文亚洲av片在线观看爽| 无人区码免费观看不卡| 国产主播在线观看一区二区| 乱码一卡2卡4卡精品| or卡值多少钱| 又紧又爽又黄一区二区| 欧美xxxx性猛交bbbb| av在线观看视频网站免费| 在线播放国产精品三级| 淫秽高清视频在线观看| 国产精品99久久久久久久久| 免费观看的影片在线观看| 日本五十路高清| 亚洲精品乱码久久久v下载方式| 亚洲综合色惰| 午夜免费男女啪啪视频观看 | 国产高清视频在线观看网站| 欧美日韩国产亚洲二区| 国产高清有码在线观看视频| av在线蜜桃| 久久久国产成人精品二区| 不卡一级毛片| 老司机午夜十八禁免费视频| 免费av毛片视频| 久久久久久久久中文| 听说在线观看完整版免费高清| 女同久久另类99精品国产91| 老司机午夜十八禁免费视频| 精品一区二区免费观看| 免费在线观看日本一区| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区性色av| 51国产日韩欧美| 黄色女人牲交| 看黄色毛片网站| 毛片一级片免费看久久久久 | 国产男靠女视频免费网站| 国产单亲对白刺激| 亚洲熟妇中文字幕五十中出| 久久久久国产精品人妻aⅴ院| 青草久久国产| 又粗又爽又猛毛片免费看| 色尼玛亚洲综合影院| 淫秽高清视频在线观看| 欧美高清成人免费视频www| 色哟哟·www| 国产成人a区在线观看| 波多野结衣高清作品| 国产探花在线观看一区二区| 亚洲精品日韩av片在线观看| 成年女人毛片免费观看观看9| 亚洲 国产 在线| 1024手机看黄色片| 午夜福利高清视频| 在线看三级毛片| 丁香欧美五月| 国产精品久久电影中文字幕| 精品国产三级普通话版| 国产在视频线在精品| 亚洲18禁久久av| 国产精品野战在线观看| 久久久精品欧美日韩精品| 欧美激情国产日韩精品一区| 69人妻影院| 国产精品亚洲美女久久久| 国产高清有码在线观看视频| 欧美不卡视频在线免费观看| 欧美性感艳星| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 伊人久久精品亚洲午夜| 男人的好看免费观看在线视频| 日韩成人在线观看一区二区三区| 色尼玛亚洲综合影院| 听说在线观看完整版免费高清| 成年女人看的毛片在线观看| 亚洲精品456在线播放app | 亚洲精品影视一区二区三区av| 亚洲精品在线观看二区| 99久国产av精品| 久久久久久国产a免费观看| 在线a可以看的网站| 成人鲁丝片一二三区免费| 男女床上黄色一级片免费看| 久久久久免费精品人妻一区二区| 男人和女人高潮做爰伦理| 男女下面进入的视频免费午夜| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩卡通动漫| 男人舔奶头视频| 亚洲人成网站高清观看| 两个人视频免费观看高清| 999久久久精品免费观看国产| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利免费观看在线| 少妇高潮的动态图| 日本撒尿小便嘘嘘汇集6| 天堂动漫精品| 亚洲av中文字字幕乱码综合| 免费搜索国产男女视频| 国产男靠女视频免费网站| 亚洲av电影不卡..在线观看| 国产一区二区三区在线臀色熟女| 成年女人看的毛片在线观看| 中文字幕精品亚洲无线码一区| 久久精品国产亚洲av香蕉五月| 欧美乱色亚洲激情| 日韩av在线大香蕉| 欧美一区二区精品小视频在线| 级片在线观看| 国产高清有码在线观看视频| 午夜两性在线视频| 成年免费大片在线观看| 久久久久久久久大av| 啦啦啦观看免费观看视频高清| 欧美成人性av电影在线观看| 97碰自拍视频|