• <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)分析
    久久久久久久午夜电影| 天天躁日日操中文字幕| 在线播放无遮挡| 成年人黄色毛片网站| 国产高清激情床上av| 久久国产精品影院| a级一级毛片免费在线观看| 在现免费观看毛片| 成人特级av手机在线观看| 国产精品久久久久久久电影| 99在线视频只有这里精品首页| 精品久久久久久久久亚洲 | 中出人妻视频一区二区| 亚洲五月天丁香| 国产高清视频在线观看网站| 亚洲成人免费电影在线观看| 大型黄色视频在线免费观看| 亚洲熟妇熟女久久| 欧美成狂野欧美在线观看| 精品久久久久久成人av| 天堂av国产一区二区熟女人妻| 在线观看免费视频日本深夜| 丰满的人妻完整版| 91久久精品国产一区二区成人| 97超级碰碰碰精品色视频在线观看| 男女那种视频在线观看| 全区人妻精品视频| 精品一区二区免费观看| 毛片一级片免费看久久久久 | 1024手机看黄色片| 夜夜躁狠狠躁天天躁| 亚洲最大成人中文| 97人妻精品一区二区三区麻豆| 欧美日韩中文字幕国产精品一区二区三区| 18禁黄网站禁片午夜丰满| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 神马国产精品三级电影在线观看| 精品久久国产蜜桃| 国产美女午夜福利| 啦啦啦韩国在线观看视频| 精品不卡国产一区二区三区| 亚洲激情在线av| 亚洲欧美日韩卡通动漫| 内地一区二区视频在线| 亚洲精品一卡2卡三卡4卡5卡| 18禁黄网站禁片免费观看直播| 噜噜噜噜噜久久久久久91| 男人和女人高潮做爰伦理| 国产精品永久免费网站| 久久香蕉精品热| 久久精品夜夜夜夜夜久久蜜豆| 日本免费a在线| 非洲黑人性xxxx精品又粗又长| 啦啦啦观看免费观看视频高清| 九九在线视频观看精品| 国产欧美日韩一区二区三| 欧美成人a在线观看| 国产v大片淫在线免费观看| 在线播放无遮挡| 此物有八面人人有两片| 听说在线观看完整版免费高清| 精品人妻熟女av久视频| 在线播放无遮挡| 久久精品人妻少妇| 最近视频中文字幕2019在线8| 精品人妻熟女av久视频| 搡老熟女国产l中国老女人| 成人国产一区最新在线观看| 最近视频中文字幕2019在线8| 在线免费观看不下载黄p国产 | 99久久久亚洲精品蜜臀av| 高清日韩中文字幕在线| a在线观看视频网站| 亚洲第一电影网av| 国产一级毛片七仙女欲春2| 丰满人妻一区二区三区视频av| 久久亚洲真实| 亚洲天堂国产精品一区在线| 动漫黄色视频在线观看| 亚洲熟妇熟女久久| 夜夜夜夜夜久久久久| 18+在线观看网站| 国产爱豆传媒在线观看| 亚洲狠狠婷婷综合久久图片| 欧美丝袜亚洲另类 | 91久久精品电影网| 90打野战视频偷拍视频| 狂野欧美白嫩少妇大欣赏| 黄色丝袜av网址大全| 网址你懂的国产日韩在线| 国内揄拍国产精品人妻在线| 亚洲欧美清纯卡通| 国产大屁股一区二区在线视频| 黄色日韩在线| 97超视频在线观看视频| 久99久视频精品免费| 天堂av国产一区二区熟女人妻| 日韩人妻高清精品专区| 嫩草影院精品99| 能在线免费观看的黄片| 国产伦一二天堂av在线观看| 一a级毛片在线观看| 99国产精品一区二区蜜桃av| 天堂影院成人在线观看| 成人毛片a级毛片在线播放| 变态另类丝袜制服| 桃色一区二区三区在线观看| 黄色一级大片看看| 美女黄网站色视频| 少妇人妻精品综合一区二区 | 日本精品一区二区三区蜜桃| 无人区码免费观看不卡| 久久精品国产亚洲av涩爱 | 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 国产一区二区激情短视频| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 精品人妻1区二区| 久久精品综合一区二区三区| 一级av片app| 国产探花极品一区二区| 精品一区二区免费观看| 国产成人aa在线观看| av专区在线播放| 国产三级中文精品| 两个人视频免费观看高清| 精品99又大又爽又粗少妇毛片 | 亚洲美女视频黄频| 国产真实乱freesex| 亚洲中文日韩欧美视频| 精品久久久久久久久亚洲 | 又爽又黄a免费视频| 成人特级黄色片久久久久久久| 在线免费观看的www视频| 禁无遮挡网站| 赤兔流量卡办理| 亚洲av成人不卡在线观看播放网| 精品乱码久久久久久99久播| 国产精品嫩草影院av在线观看 | 欧美中文日本在线观看视频| 99视频精品全部免费 在线| 黄片小视频在线播放| 免费人成视频x8x8入口观看| 久久久久久久亚洲中文字幕 | 欧美黄色淫秽网站| 亚洲精品一卡2卡三卡4卡5卡| 国产精品免费一区二区三区在线| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 18禁黄网站禁片免费观看直播| 最后的刺客免费高清国语| 啦啦啦韩国在线观看视频| 国产男靠女视频免费网站| 91久久精品国产一区二区成人| 十八禁国产超污无遮挡网站| 又粗又爽又猛毛片免费看| 少妇裸体淫交视频免费看高清| 亚洲人成网站在线播| 熟女电影av网| 深爱激情五月婷婷| 国内精品久久久久精免费| 日韩人妻高清精品专区| 69人妻影院| av在线老鸭窝| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 亚洲国产精品合色在线| 精品久久久久久久久av| 国产v大片淫在线免费观看| 午夜两性在线视频| 99国产综合亚洲精品| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 动漫黄色视频在线观看| 久久久久久久精品吃奶| 丁香六月欧美| 3wmmmm亚洲av在线观看| 久久久国产成人免费| 好看av亚洲va欧美ⅴa在| 国产蜜桃级精品一区二区三区| 最后的刺客免费高清国语| 亚洲天堂国产精品一区在线| 深夜精品福利| 欧美一区二区亚洲| 精品一区二区三区视频在线| 国内久久婷婷六月综合欲色啪| 午夜久久久久精精品| 久久久久免费精品人妻一区二区| 在线天堂最新版资源| 亚洲国产精品sss在线观看| 午夜福利免费观看在线| 午夜激情福利司机影院| 欧美最黄视频在线播放免费| .国产精品久久| 精品国产三级普通话版| 我要看日韩黄色一级片| 在线看三级毛片| 搡老熟女国产l中国老女人| 欧美三级亚洲精品| 日本 av在线| 久久天躁狠狠躁夜夜2o2o| 亚洲av电影不卡..在线观看| 一本综合久久免费| 国产av在哪里看| 黄片小视频在线播放| 麻豆国产av国片精品| 最近最新中文字幕大全电影3| 久久香蕉精品热| 久久精品影院6| 亚洲精品久久国产高清桃花| 天堂av国产一区二区熟女人妻| av欧美777| 简卡轻食公司| 亚洲av二区三区四区| 丰满乱子伦码专区| 精品久久久久久久久久免费视频| 最近中文字幕高清免费大全6 | 99久国产av精品| 两个人的视频大全免费| 99riav亚洲国产免费| 亚洲无线观看免费| 亚洲国产欧洲综合997久久,| 高清在线国产一区| 五月伊人婷婷丁香| 精品无人区乱码1区二区| 国产老妇女一区| 亚洲狠狠婷婷综合久久图片| 日日夜夜操网爽| 亚洲精品成人久久久久久| 国产亚洲精品久久久com| 在线十欧美十亚洲十日本专区| 国产成年人精品一区二区| 国产单亲对白刺激| 亚洲精品在线观看二区| 简卡轻食公司| 999久久久精品免费观看国产| 最近中文字幕高清免费大全6 | 日韩国内少妇激情av| 宅男免费午夜| 老司机午夜福利在线观看视频| 三级国产精品欧美在线观看| 中亚洲国语对白在线视频| 高潮久久久久久久久久久不卡| 国内少妇人妻偷人精品xxx网站| 亚洲精品影视一区二区三区av| av专区在线播放| 精品一区二区免费观看| 中文字幕免费在线视频6| 日韩欧美 国产精品| 最近最新免费中文字幕在线| 亚洲一区二区三区不卡视频| 欧美精品国产亚洲| 在线播放国产精品三级| 在线观看一区二区三区| 欧美色欧美亚洲另类二区| 色哟哟哟哟哟哟| 欧美黑人巨大hd| 免费av观看视频| 国产精品免费一区二区三区在线| 直男gayav资源| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 一个人看视频在线观看www免费| 麻豆国产97在线/欧美| 亚洲精品乱码久久久v下载方式| 此物有八面人人有两片| 国产成人aa在线观看| 全区人妻精品视频| 97热精品久久久久久| 中文资源天堂在线| 国产亚洲欧美在线一区二区| 99热这里只有是精品50| 精品国产三级普通话版| 99久国产av精品| 国产高清视频在线播放一区| 桃红色精品国产亚洲av| 欧美区成人在线视频| 波多野结衣高清作品| 十八禁人妻一区二区| 特级一级黄色大片| www.色视频.com| av女优亚洲男人天堂| 国产毛片a区久久久久| 成人特级av手机在线观看| 精品久久久久久久久久免费视频| 亚洲美女视频黄频| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 首页视频小说图片口味搜索| 成人毛片a级毛片在线播放| 岛国在线免费视频观看| 国产精品野战在线观看| 麻豆国产97在线/欧美| 国产色爽女视频免费观看| 亚洲国产精品sss在线观看| 精品人妻熟女av久视频| 日韩免费av在线播放| 嫩草影院入口| 在线十欧美十亚洲十日本专区| 久久精品综合一区二区三区| a级一级毛片免费在线观看| 深夜精品福利| 九色国产91popny在线| 午夜福利在线在线| 中文字幕精品亚洲无线码一区| 黄色日韩在线| 欧美xxxx黑人xx丫x性爽| 内射极品少妇av片p| 哪里可以看免费的av片| 窝窝影院91人妻| 久久久久久久久久黄片| 99热这里只有是精品50| 精品乱码久久久久久99久播| 又黄又爽又刺激的免费视频.| 一区福利在线观看| 欧美高清性xxxxhd video| av国产免费在线观看| 免费在线观看成人毛片| 亚洲真实伦在线观看| 狂野欧美白嫩少妇大欣赏| 免费在线观看影片大全网站| 久久性视频一级片| 亚洲国产高清在线一区二区三| 性色av乱码一区二区三区2| 一本精品99久久精品77| 国产伦在线观看视频一区| 午夜影院日韩av| 乱人视频在线观看| 欧美中文日本在线观看视频| 三级男女做爰猛烈吃奶摸视频| 天堂√8在线中文| 欧美黄色淫秽网站| 亚洲中文字幕日韩| 舔av片在线| 国产在视频线在精品| 国内精品美女久久久久久| 亚洲熟妇熟女久久| 成人av在线播放网站| 午夜免费男女啪啪视频观看 | 高清在线国产一区| 日韩精品青青久久久久久| 日本 欧美在线| 国产精品人妻久久久久久| 国内揄拍国产精品人妻在线| 亚洲人成电影免费在线| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 成人鲁丝片一二三区免费| 亚洲电影在线观看av| 观看美女的网站| 十八禁国产超污无遮挡网站| 日韩成人在线观看一区二区三区| 亚洲成av人片免费观看| 色5月婷婷丁香| 亚洲av免费在线观看| 中出人妻视频一区二区| 成人毛片a级毛片在线播放| 男女视频在线观看网站免费| 窝窝影院91人妻| 久久精品影院6| 看免费av毛片| 久久精品综合一区二区三区| 色综合婷婷激情| a级毛片a级免费在线| 在线国产一区二区在线| 国产一区二区三区在线臀色熟女| 99久久99久久久精品蜜桃| 国产精品日韩av在线免费观看| 看免费av毛片| 午夜福利免费观看在线| 国产精品亚洲美女久久久| 欧美日韩乱码在线| 少妇熟女aⅴ在线视频| 一个人免费在线观看电影| 日韩国内少妇激情av| 69人妻影院| 国产毛片a区久久久久| 91久久精品国产一区二区成人| 看十八女毛片水多多多| 久久久久九九精品影院| 日韩欧美 国产精品| 欧美日韩福利视频一区二区| 久久人人精品亚洲av| 免费在线观看成人毛片| 丁香欧美五月| 中文字幕免费在线视频6| 久久这里只有精品中国| 成人亚洲精品av一区二区| 最新中文字幕久久久久| 别揉我奶头 嗯啊视频| 色噜噜av男人的天堂激情| 内地一区二区视频在线| 在线免费观看的www视频| 亚洲在线自拍视频| 欧美性猛交黑人性爽| 国产亚洲精品综合一区在线观看| 色综合站精品国产| 国产亚洲精品av在线| 亚洲av一区综合| 757午夜福利合集在线观看| 欧美成人性av电影在线观看| av福利片在线观看| 在现免费观看毛片| 神马国产精品三级电影在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲美女视频黄频| 精品一区二区三区av网在线观看| 亚洲av日韩精品久久久久久密| 国产精品不卡视频一区二区 | 一个人看的www免费观看视频| 成人无遮挡网站| 亚洲第一区二区三区不卡| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 免费观看的影片在线观看| 精品国产亚洲在线| 成人av在线播放网站| 亚洲av成人不卡在线观看播放网| 老司机午夜福利在线观看视频| 最新中文字幕久久久久| 国产精品一区二区免费欧美| 麻豆国产97在线/欧美| 9191精品国产免费久久| 丝袜美腿在线中文| 亚洲第一区二区三区不卡| 日韩有码中文字幕| 精品人妻视频免费看| 亚洲人成网站在线播| 热99在线观看视频| 久久草成人影院| 亚洲人成网站在线播放欧美日韩| 又黄又爽又免费观看的视频| 久久中文看片网| 国产日本99.免费观看| 国产午夜福利久久久久久| 草草在线视频免费看| 欧美又色又爽又黄视频| 国模一区二区三区四区视频| 搞女人的毛片| 九九热线精品视视频播放| 国产色爽女视频免费观看| 久久这里只有精品中国| 亚洲第一欧美日韩一区二区三区| 国产三级中文精品| 好男人在线观看高清免费视频| 亚洲经典国产精华液单 | 99久久成人亚洲精品观看| www日本黄色视频网| 麻豆av噜噜一区二区三区| 久久精品人妻少妇| 午夜福利视频1000在线观看| 精品福利观看| 99国产精品一区二区蜜桃av| 亚洲最大成人中文| 国产黄色小视频在线观看| 日韩欧美在线乱码| 亚洲精品粉嫩美女一区| 脱女人内裤的视频| 久久久精品欧美日韩精品| 国产伦一二天堂av在线观看| 国产精品日韩av在线免费观看| 成人欧美大片| 每晚都被弄得嗷嗷叫到高潮| 久久久久亚洲av毛片大全| 成年人黄色毛片网站| 内射极品少妇av片p| eeuss影院久久| aaaaa片日本免费| 每晚都被弄得嗷嗷叫到高潮| 欧美乱妇无乱码| 国产毛片a区久久久久| 亚洲国产欧洲综合997久久,| 久久久久精品国产欧美久久久| avwww免费| 99热这里只有精品一区| 搡女人真爽免费视频火全软件 | 亚洲精品一卡2卡三卡4卡5卡| 男女之事视频高清在线观看| 18+在线观看网站| 婷婷六月久久综合丁香| 久久伊人香网站| 精品一区二区免费观看| av在线老鸭窝| 亚洲激情在线av| 国内精品久久久久久久电影| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 免费av不卡在线播放| 女人十人毛片免费观看3o分钟| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 欧美丝袜亚洲另类 | 成人三级黄色视频| 欧美中文日本在线观看视频| 永久网站在线| 丝袜美腿在线中文| h日本视频在线播放| 亚洲五月天丁香| 一个人看的www免费观看视频| 夜夜爽天天搞| 亚洲av五月六月丁香网| 国产精品久久视频播放| 欧美最新免费一区二区三区 | 免费人成在线观看视频色| 变态另类丝袜制服| 色综合婷婷激情| 18禁在线播放成人免费| 国内精品一区二区在线观看| 91久久精品国产一区二区成人| 免费无遮挡裸体视频| 51午夜福利影视在线观看| 熟妇人妻久久中文字幕3abv| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看| 俄罗斯特黄特色一大片| 精品国内亚洲2022精品成人| 免费高清视频大片| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 亚洲男人的天堂狠狠| 亚洲综合色惰| 天堂影院成人在线观看| 久久久成人免费电影| 99热这里只有是精品50| 长腿黑丝高跟| 国产伦精品一区二区三区四那| 九九热线精品视视频播放| 日韩欧美在线乱码| 欧美日韩中文字幕国产精品一区二区三区| 免费av观看视频| 中文字幕熟女人妻在线| 国产成+人综合+亚洲专区| 超碰av人人做人人爽久久| 国产精品亚洲美女久久久| 91久久精品国产一区二区成人| 少妇高潮的动态图| 最近最新免费中文字幕在线| 网址你懂的国产日韩在线| 亚洲精品久久国产高清桃花| 亚洲 欧美 日韩 在线 免费| 午夜福利在线观看吧| 在线观看66精品国产| 午夜免费激情av| 欧美最新免费一区二区三区 | 亚洲欧美日韩高清在线视频| 亚洲av免费在线观看| 日日摸夜夜添夜夜添小说| 国产精品女同一区二区软件 | av天堂中文字幕网| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 我要搜黄色片| 在线国产一区二区在线| 美女xxoo啪啪120秒动态图 | 日日干狠狠操夜夜爽| 少妇高潮的动态图| 国产成人影院久久av| 12—13女人毛片做爰片一| 国产中年淑女户外野战色| 嫩草影视91久久| 亚洲欧美激情综合另类| 国产野战对白在线观看| 特大巨黑吊av在线直播| 男人舔女人下体高潮全视频| 国产色爽女视频免费观看| 深夜精品福利| 亚洲中文日韩欧美视频| 中文字幕熟女人妻在线| 伊人久久精品亚洲午夜| 日韩av在线大香蕉| 免费av不卡在线播放| 一本久久中文字幕| 三级国产精品欧美在线观看| 97超级碰碰碰精品色视频在线观看| 成人三级黄色视频| 国产成人啪精品午夜网站| 国产精品,欧美在线| 国产大屁股一区二区在线视频| 免费高清视频大片| av专区在线播放| 欧美性猛交黑人性爽| 内射极品少妇av片p| 最近最新免费中文字幕在线| 两个人的视频大全免费| x7x7x7水蜜桃| 亚洲美女黄片视频| 亚洲精品乱码久久久v下载方式| 十八禁网站免费在线| 国产精品不卡视频一区二区 | 国产毛片a区久久久久| 男人的好看免费观看在线视频| 91久久精品电影网| 麻豆av噜噜一区二区三区| eeuss影院久久| 亚洲第一电影网av| 久久久久久久午夜电影| 国产亚洲精品av在线| 一个人观看的视频www高清免费观看| 精品不卡国产一区二区三区| 18禁黄网站禁片午夜丰满| 美女大奶头视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 69人妻影院| 在线观看66精品国产| bbb黄色大片| 美女大奶头视频| 亚洲第一电影网av| 亚洲熟妇中文字幕五十中出| 看黄色毛片网站| 两个人视频免费观看高清| 小说图片视频综合网站|