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

    宜昌臘腸與恩施臘腸細(xì)菌多樣性的比較分析

    2020-12-28 03:04周亞澳黃怡王韻博郭壯張振東李楊坤侯強(qiáng)川
    肉類研究 2020年10期
    關(guān)鍵詞:高通量測(cè)序臘腸

    周亞澳 黃怡 王韻博 郭壯 張振東 李楊坤 侯強(qiáng)川

    摘 要:采用Illumina MiSeq高通量測(cè)序技術(shù)對(duì)宜昌臘腸細(xì)菌多樣性進(jìn)行解析,并與之前報(bào)道的恩施臘腸細(xì)菌群落結(jié)構(gòu)進(jìn)行比較分析。結(jié)果表明:宜昌臘腸在門水平主要由厚壁菌門(Firmicutes,64.93%)、變形菌門(Proteobacteria,32.12%)和放線菌門(Actinobacteria,1.14%)組成;在屬水平上,主要由葡萄球菌屬(Staphylococcus,16.57%)、乳酸桿菌屬(Lactobacillus,16.22%)、魏斯氏菌屬(Weissella,14.67%)和明串珠菌屬(Leuconostoc,11.43%)等組成。宜昌臘腸細(xì)菌豐度顯著低于恩施臘腸(P<0.05),且2 個(gè)地區(qū)臘腸的細(xì)菌菌群結(jié)構(gòu)存在較大差異?;蚬δ茴A(yù)測(cè)結(jié)果表明,宜昌臘腸中細(xì)菌在碳水化合物代謝和多糖的生物合成與代謝等方面更加旺盛,這些差異與臘腸中乳酸桿菌的相對(duì)豐度存在顯著相關(guān)性。

    關(guān)鍵詞:臘腸;高通量測(cè)序;微生物群落結(jié)構(gòu);功能預(yù)測(cè)

    Abstract: The bacterial diversity in Yichang sausage was analyzed by Illumina MiSeq high-throughput sequencing technique and compared with that previously reported for Enshi sausage. The results showed that the bacterial community in Yichang sausages were mainly composed of Firmicutes (64.93%), Proteobacteria (32.12%) and Actinobacteria (1.14%) at the phylum level. At the genus level, it was mainly composed of Staphylococcus (16.57%), Lactobacillus (16.22%), Weissella (14.67%), and Leuconostoc (11.43%). The abundance of bacteria in Yichang sausages was significantly lower than that in Enshi sausage (P < 0.05). Besides, there were great differences in the bacterial structure between the two sausages. The results of gene function prediction showed that bacteria in Yichang sausages were more vigorous in carbohydrate metabolism, and polysaccharide biosynthesis and metabolism, and these differences were significantly correlated with the relative abundance of Lactobacillus in sausages.

    Keywords: sausages; high throughput sequencing; microbial community structure; function prediction

    DOI:10.7506/rlyj1001-8123-20200711-169

    中圖分類號(hào):TS251.5? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1001-8123(2020)10-0008-06

    臘腸作為我國(guó)一種傳統(tǒng)的發(fā)酵肉制品,擁有悠久的制作歷史,現(xiàn)廣泛分布于我國(guó)華東、華中和華南等地區(qū)[1]。臘腸的制作主要以豬肉為原料,同時(shí)輔以食鹽、八角、花椒和白酒等調(diào)味料進(jìn)行腌制,灌入腸衣后進(jìn)行晾曬發(fā)酵[2]。因其風(fēng)味獨(dú)特、口感醇厚而深受廣大消費(fèi)者的喜愛(ài)[3]。

    宜昌市位于湖北省西南部,地處長(zhǎng)江上游與中游結(jié)合部,其氣候?yàn)榕D腸的發(fā)酵提供了有利條件。同時(shí)宜昌市是一個(gè)多民族雜居的城市,這在一定程度促進(jìn)了傳統(tǒng)臘腸制作工藝的融合發(fā)展。傳統(tǒng)臘腸的制作環(huán)境較為粗獷,不同的地理環(huán)境和氣候條件可能會(huì)賦予臘腸不同的微生物群系[4]。目前,一些學(xué)者針對(duì)不同品質(zhì)、原料及發(fā)酵和貯藏階段臘腸微生物群落組成和差異開展一些研究,滕安國(guó)等[5]發(fā)現(xiàn)乳桿菌屬(Lactobacillus)和明串珠菌屬(Leuconostoc)為正常臘腸和腸衣中的優(yōu)勢(shì)菌屬,而臘腸發(fā)黏變質(zhì)后微生物組成發(fā)生變化,菌群多樣性降低;劉長(zhǎng)建等[6]從臘腸中分離得到多株具有降膽固醇能力的乳酸菌,其中干酪乳桿菌(Lactobacillus casei)的降膽固醇效果最好;田甜等[7]發(fā)現(xiàn)在廣式臘腸發(fā)酵過(guò)程中添加復(fù)合發(fā)酵劑能加快發(fā)酵進(jìn)程,并降低亞硝酸鹽含量。Fougy等[8]發(fā)現(xiàn)減少臘腸中的食鹽添加量會(huì)導(dǎo)致腐敗菌生長(zhǎng)和細(xì)菌多樣性下降;Li Xinfu等[9]通過(guò)高通量測(cè)序技術(shù)解析真空包裝臘腸貯藏過(guò)程中微生物群落的變化。此外,越來(lái)越多的研究表明,臘腸的發(fā)酵菌群直接影響著發(fā)酵的進(jìn)程和產(chǎn)品品質(zhì),郭壯等[10]研究發(fā)現(xiàn)臘腸中的乳酸菌對(duì)產(chǎn)品發(fā)酵成熟和風(fēng)味品質(zhì)的形成具有至關(guān)重要的作用;黃金枝[11]研究發(fā)酵菌種與廣式臘腸產(chǎn)品品質(zhì)的關(guān)系,發(fā)現(xiàn)戊糖片球菌(Pediococcus pentosaceus)對(duì)廣式臘腸發(fā)酵品質(zhì)的影響大于漢遜德巴利酵母(Saccharomyces)和植物乳桿菌(Lactobacillus plantarum)。上述研究的結(jié)果表明,臘腸品質(zhì)與其細(xì)菌組成之間存在密切聯(lián)系。然而,目前鮮有研究報(bào)道宜昌地區(qū)臘腸細(xì)菌多樣性及其與其他地區(qū)臘腸菌群結(jié)構(gòu)的差異,這在一定程度上限制了宜昌地區(qū)臘腸品質(zhì)的提升。

    近年來(lái),伴隨分子生物學(xué)的發(fā)展,越來(lái)越多的基于非培養(yǎng)的微生物檢測(cè)方法和技術(shù)被廣泛應(yīng)用于發(fā)酵食品微生物檢測(cè)中。其中,Illumina MiSeq高通量測(cè)序技術(shù)因其檢測(cè)速率快、通量大和結(jié)果可信度高等優(yōu)點(diǎn)[12-13],在海洋、土壤、發(fā)酵食品和腸道等環(huán)境樣品微生物多樣性的研究中得到廣泛應(yīng)用[14-15]。鄧風(fēng)等[16]采用聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction,PCR)-變性梯度凝膠電泳與Illumina MiSeq高通量測(cè)序技術(shù)相結(jié)合的方法研究恩施地區(qū)臘腸中的細(xì)菌群落結(jié)構(gòu),發(fā)現(xiàn)恩施臘腸樣品中含有大量的共有細(xì)菌類群,優(yōu)勢(shì)細(xì)菌屬主要包括環(huán)絲菌屬(Brochothrix)、葡萄球菌(Staphylococcus)、嗜冷桿菌(Psychrobacter)等,同時(shí)恩施臘腸所有測(cè)序序列中無(wú)法鑒定到屬水平的僅為13.72%。該研究結(jié)果進(jìn)一步表明Illumina MiSeq高通量測(cè)序可以作為解析臘腸細(xì)菌群落結(jié)構(gòu)的有效技術(shù)手段。

    本研究采用Illumina MiSeq高通量測(cè)序技術(shù)對(duì)采集自宜昌地區(qū)的臘腸細(xì)菌多樣性進(jìn)行解析。此外,通過(guò)整合恩施臘腸測(cè)序數(shù)據(jù)[17],結(jié)合多元統(tǒng)計(jì)學(xué)分析和功能預(yù)測(cè)等手段比較分析2 個(gè)地區(qū)臘腸細(xì)菌菌群結(jié)構(gòu)和功能的差異。以期進(jìn)一步豐富人們對(duì)不同地區(qū)臘腸細(xì)菌群落結(jié)構(gòu)的認(rèn)識(shí)。

    1 材料與方法

    1.1 材料與試劑

    臘腸采集自湖北省宜昌市五峰土家族自治縣菜市場(chǎng),共計(jì)9 份,編號(hào)記為YC1~YC9,所有樣品均為手工制作,且滿足以下標(biāo)準(zhǔn):1)臘腸產(chǎn)品在宜昌地區(qū)生產(chǎn)制作完成;2)臘腸使用豬肉作為制作原料;3)臘腸產(chǎn)品無(wú)霉變、無(wú)異味且無(wú)明顯雜質(zhì)存在。樣品采集時(shí)使用無(wú)菌手術(shù)刀切取約100 g裝入無(wú)菌采樣袋,密封后置于采樣箱中,低溫冷藏運(yùn)回實(shí)驗(yàn)室進(jìn)行后續(xù)操作。

    338F/806R引物 武漢天一輝遠(yuǎn)生物科技有限公司;DNeasy Mericon Food Kit 德國(guó)Qiagen公司;FastPfu Fly DNA Polymerase、5×TransStart? FastPfu緩沖液和dNTPs Mix 北京全式金生物技術(shù)有限公司。

    1.2 儀器與設(shè)備

    Vetiri梯度基因擴(kuò)增儀 美國(guó)AB SCIEX公司;ND-2000C微量紫外分光光度計(jì) 美國(guó)Nano Drop公司;Illumina MiSeq PE300高通量測(cè)序平臺(tái) 美國(guó)Illumina公司;R960機(jī)架式服務(wù)器 美國(guó)Dell公司;UVPCDS8000凝膠成像分析系統(tǒng)、DCode? System 美國(guó)Bio-Rad公司;CR21N型高速離心機(jī) 日本日立金屬株式會(huì)社。

    1.3 方法

    1.3.1 樣品預(yù)處理及微生物宏基因組DNA提取

    將采集的臘腸切碎后,稱取5 g加入到45 mL生理鹽水中,使用拍擊器拍擊5 min,400 r/min離心8 min后取上清液,于12 000 r/min離心8 min,取沉淀備用[16]。使用基因組提取試劑盒提取沉淀物宏基因組DNA。

    1.3.2 PCR擴(kuò)增和Illumina MiSeq測(cè)序

    使用338F/806R引物(帶有7 個(gè)堿基核苷酸標(biāo)簽)對(duì)細(xì)菌16S rRNA V3~V4區(qū)進(jìn)行PCR擴(kuò)增,擴(kuò)增條件參照王玉榮等[17]方法。使用1.00%的瓊脂糖凝膠電泳和紫外分光光度計(jì)對(duì)擴(kuò)增產(chǎn)物的純度和濃度進(jìn)行檢驗(yàn),將合格DNA產(chǎn)物送至上海美吉生物醫(yī)藥科技有限公司進(jìn)行測(cè)序。

    1.3.3 生物信息學(xué)分析

    依照郭壯等[18]研究中序列質(zhì)控所約束的參數(shù)對(duì)下機(jī)序列進(jìn)行質(zhì)量控制,并依照核苷酸標(biāo)簽進(jìn)行序列的分配,將質(zhì)控后的序列分配到不同的樣本中。參照Yang Chengcong等[19]的方法,使用QIIME V1.9.1平臺(tái)對(duì)宜昌和恩施臘腸中的細(xì)菌組成和多樣性進(jìn)行解析。具體的生物學(xué)分析過(guò)程如下:1)采用PyNAST軟件對(duì)序列數(shù)據(jù)進(jìn)行比對(duì)和校準(zhǔn)[20];2)使用UCLUST分別通過(guò)100%和97%相似度對(duì)校準(zhǔn)后的序列進(jìn)行劃分[21],分別得到非冗余數(shù)據(jù)集和操作分類單元(operational taxonomic units,OTU);3)使用ChimeraSlayer軟件識(shí)別OTU中的嵌合體序列[22],并對(duì)相應(yīng)的嵌合體序列進(jìn)行刪除;4)挑選每個(gè)OTU中的代表性序列分別使用Greengene(Version 13.8)[23]、RDP(Relace 11.5)[24]和Sliva(Version 132)數(shù)據(jù)庫(kù)[25]進(jìn)行同源性比對(duì),確定各OTU的分類學(xué)地位,使用內(nèi)部腳本整合3 個(gè)數(shù)據(jù)庫(kù)的注釋結(jié)果,確定各OTU的最終分類學(xué)地位;5)使用FastTree軟件構(gòu)建基于OTU代表性序列的系統(tǒng)發(fā)育樹[26];6)依據(jù)最大序列數(shù)計(jì)算每個(gè)臘腸樣本的α多樣性指數(shù)(Chao 1指數(shù)、Shannon指數(shù)、Observed species指數(shù)和Simpson指數(shù));7)依據(jù)最小序列數(shù)計(jì)算樣品間的β多樣性。

    使用PICRUSt(Phylogenetic Investigation of Communities by Reconstruction of Unobserved States)軟件對(duì)臘腸中細(xì)菌微生物的基因功能進(jìn)行預(yù)測(cè)[27],并參照京都基因與基因組百科全書(Kyoto Encyclopedia of Genes and Genomes,KEGG)數(shù)據(jù)庫(kù)進(jìn)行基因功能注釋。

    1.3.4 多元統(tǒng)計(jì)學(xué)分析

    使用非參數(shù)的Mann-Whitney秩和檢驗(yàn)對(duì)不同地區(qū)臘腸樣品菌群組成進(jìn)行差異顯著性分析;使用LefSe軟件分析并基于LDA Effect Size判別組間豐度有顯著差異的微生物;基于Pearson相關(guān)性分析臘腸中優(yōu)勢(shì)菌屬和差異代謝通路之間的相關(guān)性。使用R軟件(v3.5.0)、STAMP軟件和Origin 2019b軟件進(jìn)行數(shù)據(jù)可視化處理。

    2 結(jié)果與分析

    2.1 宜昌臘腸序列豐富度和多樣性分析

    宜昌臘腸所有樣品共測(cè)得序列427 570 條,每個(gè)樣品的平均測(cè)序量為26 009 條(12 215~40 828,標(biāo)準(zhǔn)偏差3 192),上述序列在97%相似度劃分得到5 330 個(gè)非嵌合體OTU,每個(gè)樣本平均OUT數(shù)為931。

    由表1可知,YC9臘腸的OTU數(shù)和Chao 1指數(shù)最大,而YC2臘腸的Shannon指數(shù)最大。通過(guò)對(duì)OTU數(shù)、Chao 1指數(shù)和Shannon指數(shù)的變異系數(shù)進(jìn)行計(jì)算發(fā)現(xiàn),宜昌臘腸在α多樣性指數(shù)上存在較大的組內(nèi)差異。相關(guān)研究顯示,Chao 1指數(shù)和Shannon指數(shù)常分別用來(lái)評(píng)估環(huán)境中微生物的豐度和多樣性[28]。表明YC9臘腸中微生物豐度最高,而YC2臘腸的微生物多樣性最高,且宜昌臘腸中的微生物豐度和多樣性存在較大差異。

    進(jìn)一步分析宜昌臘腸中的細(xì)菌菌落組成,所有樣品共注釋得到15 個(gè)菌門。由圖1可知,宜昌臘腸中優(yōu)勢(shì)菌門(平均相對(duì)豐度大于1.00%)共有3 個(gè),分別為厚壁菌門(Firmicutes)、變形菌門(Proteobacteria)和放線菌門(Actinobacteria),平均相對(duì)豐度分別為64.93%、32.12%和1.14%??梢?jiàn),宜昌臘腸中的主要菌門為厚壁菌門和變形菌門,這與鄧風(fēng)等[16]研究得到的恩施臘腸細(xì)菌群落結(jié)構(gòu)結(jié)果基本一致。

    本研究進(jìn)一步在屬水平對(duì)宜昌臘腸細(xì)菌群落結(jié)構(gòu)進(jìn)行分析,由圖2可知,宜昌臘腸中優(yōu)勢(shì)菌屬(平均相對(duì)豐度大于1.00%)共有14 個(gè),其中葡萄球菌屬(Staphylococcus,16.57%)、乳酸桿菌屬(Lactobacillus,16.22%)、魏斯氏菌屬(Weissella,14.67%)、明串珠菌屬(Leuconostoc,11.43%)和環(huán)絲菌屬(Brochothrix,2.54%)隸屬于厚壁菌門;拉恩菌屬(Rahnella,7.91%)、假單胞菌屬(Pseudomonas,5.74%)、泛菌屬(Pantoea,3.16%)、克呂沃爾氏菌屬(Kluyvera,2.76%)、腸桿菌屬(Enterobacter,2.21%)、乳球菌屬(Lactococcus,1.67%)、檸檬酸桿菌屬(Citrobacter,1.43%)、拉烏爾菌屬(Raoultella,1.21%)和沙雷菌屬(Serratia,1.14%)隸屬于變形菌門。可見(jiàn),宜昌臘腸中相對(duì)豐度最高的是葡萄球菌屬和乳酸桿菌屬,這與恩施臘腸中的優(yōu)勢(shì)菌屬不盡相同,后者優(yōu)勢(shì)菌屬除葡萄球菌屬和乳酸桿菌屬外,還包括環(huán)絲菌屬,三者在樣品中的平均相對(duì)豐度分別為9.79%、2.80%和38.34%[16]。宜昌臘腸中細(xì)菌主要以乳酸菌為主,其平均相對(duì)豐度接近60%(主要以乳酸桿菌屬、魏斯氏菌屬和明串珠菌屬為主)。乳酸菌作為發(fā)酵肉制品中常見(jiàn)的發(fā)酵菌種,在一定程度上決定著發(fā)酵制品最終品質(zhì)。值得注意的是,本研究在宜昌臘腸中檢測(cè)到大量的機(jī)會(huì)性致病微生物如葡萄球菌屬和沙雷菌屬等,提示需要進(jìn)一步規(guī)范化臘腸制作過(guò)程的安全管理和質(zhì)量控制。

    2.2 宜昌和恩施臘腸中細(xì)菌群落結(jié)構(gòu)的比較分析

    由圖3可知,宜昌臘腸中細(xì)菌的Chao 1指數(shù)極顯著低于恩施臘腸(P<0.01),Observed species指數(shù)顯著低于恩施臘腸(P<0.05),Shannon指數(shù)低于恩施臘腸,且差異不顯著,而Simpson指數(shù)略高于恩施臘腸。Chao 1指數(shù)和Observed species指數(shù)常用于評(píng)估環(huán)境樣本中微生物的豐度,Shannon指數(shù)和Simpson指數(shù)常用于評(píng)估環(huán)境中微生物的多樣性[28]??梢?jiàn),宜昌臘腸中微生物的豐度明顯低于恩施臘腸,而多樣性沒(méi)有明顯差異。

    本研究進(jìn)一步通過(guò)基于非加權(quán)和加權(quán)UniFrac距離的主坐標(biāo)分析比較2 個(gè)地區(qū)臘腸菌群的組間差異。由基于非加權(quán)UniFrac距離的主坐標(biāo)分析結(jié)果(圖4A)可知,宜昌和恩施臘腸的微生物菌群在空間排布上呈現(xiàn)明顯的分離趨勢(shì),宜昌臘腸和恩施臘腸分別位于X軸的正方向和負(fù)方向,多元方差分析結(jié)果也表明二者的微生物群落結(jié)構(gòu)差異顯著(P<0.05);而基于加權(quán)UniFrac距離的分析結(jié)果(圖4B)表明,雖然2 個(gè)地區(qū)的臘腸樣品整體存在分離趨勢(shì),但部分樣品在空間排布上出現(xiàn)一定的重疊。上述結(jié)果表明,宜昌和恩施臘腸中優(yōu)勢(shì)菌群組成較為相似,但二者均含有較多獨(dú)特的低豐度菌群,這部分菌群對(duì)臘腸的整體菌群結(jié)構(gòu)產(chǎn)生較大影響。

    為揭示宜昌和恩施臘腸的菌群結(jié)構(gòu)差異,通過(guò)LefSe分析比較兩地區(qū)臘腸存在顯著差異的菌群。由圖5可知,共有7 個(gè)細(xì)菌類群(判別得分>3)在2 組臘腸樣品中存在顯著差異(P<0.05),其主要為果膠桿菌屬(Pectobacterium)、腸桿菌屬和假單胞菌屬,值得注意的是,甄別出的所有差異菌均在恩施臘腸中的相對(duì)豐度顯著較高(P<0.05)。相關(guān)報(bào)道顯示,腸桿菌屬和假單胞菌屬?gòu)V泛存在于傳統(tǒng)發(fā)酵食品中,能夠有效促進(jìn)發(fā)酵制品的發(fā)酵進(jìn)程,其中假單胞菌屬與發(fā)酵制品的顏色密切相關(guān),且對(duì)酪胺和組胺的產(chǎn)生具有一定的影響[29]。同時(shí),發(fā)酵制品中的菌群結(jié)構(gòu)與地理環(huán)境、制作工藝和原材料等有直接的關(guān)系,不同菌群結(jié)構(gòu)可能是造成不同地區(qū)臘腸色澤、滋味和風(fēng)味差異的主要原因之一。通過(guò)對(duì)不同地區(qū)發(fā)酵制品中微生物多樣性進(jìn)行解析和菌株的收集及保藏,將對(duì)傳統(tǒng)發(fā)酵食品的工業(yè)化和個(gè)性化生產(chǎn)具有重要意思。

    2.3 PICRUSt功能預(yù)測(cè)

    本研究利用PICRUSt軟件對(duì)宜昌和恩施臘腸中細(xì)菌類群的基因功能進(jìn)行預(yù)測(cè),使用KEGG數(shù)據(jù)庫(kù)對(duì)預(yù)測(cè)基因進(jìn)行功能注釋,在此基礎(chǔ)上對(duì)差異代謝通路進(jìn)行分析。Z值是用于衡量樣本均值偏離整體均值的方差倍數(shù),當(dāng)用于差異代謝通路查找時(shí),Z值越大,則相應(yīng)的代謝通路在兩組間的差異越大。本研究選擇Z值>1.6(95%置信區(qū)間)[30],同時(shí)代謝通路中的基因在某一地區(qū)樣品高表達(dá)概率>80%作為判定2 個(gè)地區(qū)差異代謝通路的條件。

    由圖6A可知,宜昌和恩施臘腸中的細(xì)菌菌群在4 條代謝通路上存在顯著差異,且所有差異代謝通路均隸屬于與代謝相關(guān)的一級(jí)功能層。其中與碳水化合物和多糖生物合成及代謝的相關(guān)基因在宜昌臘腸微生物中顯著富集,而與能量代謝和氨基酸代謝的相關(guān)基因在恩施臘腸微生物中顯著富集。由圖6B可知,魏斯氏菌屬與氨基酸代謝呈顯著負(fù)相關(guān)(P<0.05),而乳酸桿菌屬與多糖生物合成及代謝呈顯著正相關(guān)(P<0.05),與碳水化合物代謝呈極顯著正相關(guān)(P<0.01)??梢?jiàn),宜昌地區(qū)臘腸中微生物對(duì)碳水化合物的代謝更為旺盛,并且這與臘腸中的乳酸桿菌密切相關(guān)。乳酸桿菌屬對(duì)碳水化合物的利用效率受到環(huán)境的溫度、濕度等因素的強(qiáng)烈影響[31]。宜昌地區(qū)相較恩施地區(qū)空氣濕度更低,且平均溫度更高,上述環(huán)境因素均有利于臘腸中乳酸桿菌的生長(zhǎng)和代謝,這可能是宜昌臘腸中微生物碳水化合物代謝相關(guān)基因富集的原因之一。亮氨酸和異亮氨酸均為人體的必需氨基酸,由于人體無(wú)法合成,因此必須從飲食中攝取。恩施臘腸中微生物亮氨酸和異亮氨酸生物合成通路更為活躍,提示該地區(qū)生產(chǎn)的臘腸中這2 種必需氨基酸的含量可能會(huì)更為豐富。

    3 結(jié) 論

    本研究采用Illumina MiSeq高通量測(cè)序技術(shù)對(duì)湖北宜昌臘腸中的細(xì)菌多樣性進(jìn)行解析,并和恩施臘腸進(jìn)行比較分析。結(jié)果表明,宜昌臘腸中的細(xì)菌主要隸屬于厚壁菌門、變形菌門和放線菌門,總相對(duì)豐度98.19%;而優(yōu)勢(shì)菌屬主要以葡萄球菌屬、乳酸桿菌屬、魏斯氏菌屬和明串珠菌屬等為主。宜昌臘腸中的微生物豐度顯著低于恩施臘腸,且果膠桿菌屬、腸桿菌屬和假單胞菌屬等菌屬在宜昌地區(qū)的相對(duì)豐度顯著低于恩施地區(qū)。功能預(yù)測(cè)結(jié)果表明,宜昌和恩施地區(qū)臘腸與代謝相關(guān)的通路存在差異顯著,宜昌臘腸中細(xì)菌碳水化合物代謝和多糖生物合成與代謝更加旺盛,這些差異可能與2 個(gè)地區(qū)的氣候不同存在一定關(guān)系。

    參考文獻(xiàn):

    [1] 龍強(qiáng), 聶乾忠, 劉成國(guó). 發(fā)酵肉制品功能性發(fā)酵劑研究現(xiàn)狀[J]. 食品科學(xué), 2016, 37(17): 263-269. DOI:10.7506/spkx1002-6630-201617044.

    [2] 楊江, 楊成聰, 凌霞, 等. 基于電子鼻和氣相色譜-質(zhì)譜聯(lián)用技術(shù)評(píng)價(jià)襄陽(yáng)地區(qū)臘腸風(fēng)味品質(zhì)[J]. 肉類研究, 2018, 32(8): 46-50. DOI:10.7506/rlyj1001-8123-201808008.

    [3] 劉琨毅, 王琪, 王衛(wèi), 等. 茶多酚對(duì)低鹽中式臘腸防腐保鮮的影響[J]. 肉類研究, 2018, 32(3): 34-39. DOI:10.7506/rlyj1001-8123-201803007.

    [4] 杜莎, 李柯, 黃晴, 等. 人工發(fā)酵劑對(duì)湘西臘腸微生物菌群和理化特性的影響[J]. 現(xiàn)代食品科技, 2017, 33(4): 241-247. DOI:10.13982/j.mfst.1673-9078.2017.4.037.

    [5] 滕安國(guó), 齊曉娜, 張芹, 等. 高通量測(cè)序技術(shù)分析天然腸衣香腸的菌群組成[J]. 食品研究與開發(fā), 2016, 37(11): 113-116. DOI:10.3969/j.issn.1005-6521.2016.11.027.

    [6] 劉長(zhǎng)建, 蔣本國(guó), 姜波, 等. 臘腸中降膽固醇乳酸菌的篩選及鑒定[J]. 中國(guó)釀造, 2011, 30(8): 35-37. DOI:10.3969/j.issn.0254-5071.2011.08.011.

    [7] 田甜, 詹銳琪, 張雅琳, 等. 微生物發(fā)酵技術(shù)對(duì)廣味香腸安全性提升的研究[J]. 中國(guó)調(diào)味品, 2019, 44(12): 26-30. DOI:10.3969/j.issn.1000-9973.2019.12.006.

    [8] FOUGY L, DESMONTS M, COEURET G, et al. Reducing salt in raw pork sausages increases spoilage and correlates with reduced bacterial diversity[J]. Applied and Environmental Microbiology, 2016, 82(13): 3928-3939. DOI:10.1128/aem.00323-16.

    [9] LI Xinfu, LI Cong, YE Hua, et al. Changes in the microbial communities in vacuum-packaged smoked bacon during storage[J]. Food Microbiology, 2019, 77(2): 26-37. DOI:10.1016/j.fm.2018.08.007.

    [10] 郭壯, 王玉榮, 葛東穎, 等. 臘腸發(fā)酵過(guò)程中細(xì)菌多樣性評(píng)價(jià)及其對(duì)風(fēng)味的影響[J].食品科學(xué), 2020. DOI:10.7506/spkx1002-6630-20191121-255.

    [11] 黃金枝. 發(fā)酵廣式臘腸加工菌種優(yōu)選及品質(zhì)變化研究[D]. 南昌: 江西農(nóng)業(yè)大學(xué), 2015.

    [12] LIANG Huipeng, YIN Liguo, ZHANG Yahao, et al. Dynamics and diversity of a microbial community during the fermentation of industrialized Qingcai paocai, a traditional Chinese fermented vegetable food, as assessed by Illumina MiSeq sequencing, DGGE and qPCR assay[J]. Annals of Microbiology, 2018, 68(2): 111-122. DOI:10.1007/s13213-017-1321-z.

    [13] SCHIRMER M, IJAZ U Z, DAMORE R, et al. Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform[J]. Nucleic Acids Research, 2015, 43(6): e37. DOI:10.1093/nar/gku1341.

    [14] WALSH A M, FIONA C, KIERAN K, et al. Microbial succession and flavor production in the fermented dairy beverage kefir[J]. mSystems, 2016, 1(5): 2379-5077. DOI:10.1128/mSystems.00052-16.

    [15] HE Guoqing, LIU Tongjie, ASADIQ F, et al. Insights into the microbial diversity and community dynamics of Chinese traditional fermented foods from using high-throughput sequencing approaches[J]. Journal of Zhejiang University-Science B, 2017, 18(4): 289-302. DOI:10.1631/jzus.b1600148.

    [16] 鄧風(fēng), 王玉榮, 尚雪嬌, 等. 恩施地區(qū)臘腸的細(xì)菌多樣性[J]. 肉類研究, 2018, 32(9): 18-22. DOI:10.7506/rlyj1001-8123-201809004.

    [17] 王玉榮, 沈馨, 董蘊(yùn), 等. 鲊廣椒細(xì)菌多樣性評(píng)價(jià)及其對(duì)風(fēng)味的影響[J]. 食品與機(jī)械, 2018, 34(4): 25-30. DOI:10.13652/j.issn.1003-5788.2018.04.005.

    [18] 郭壯, 葛東穎, 尚雪嬌, 等. 退化和正常窖泥微生物多樣性的比較分析[J]. 食品工業(yè)科技, 2018, 39(22): 93-98. DOI:0.13386/j.issn1002-0306.2018.22.018.

    [19] YANG Chengcong, ZHAO Feiyan, HOU Qiangchuan, et al. PacBio sequencing reveals bacterial community diversity in cheeses collected from different regions[J]. Journal of Dairy Science, 2020, 103(2): 1238-1249. DOI:10.3168/jds.2019-17496.

    [20] CAPORASO J G, BITTINGER K, BUSHMAN F D, et al. PyNAST: a flexible tool for aligning sequences to a template alignment[J]. Bioinformatics, 2010, 26(2): 266-267. DOI:10.1093/bioinformatics/btp636.

    [21] EDGAR R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads[J]. Nature Methods, 2013, 10: 996. DOI:10.1038/nmeth.2604.

    [22] EDGAR R C, HAAS B J, CLEMENTE J C, et al. UCHIME improves sensitivity and speed of chimera detection[J]. Bioinformatics, 2011, 27(16): 2194-2200. DOI:10.1093/bioinformatics/btr381.

    [23] DESANTIS T Z, HUGENHOLTZ P, LARSEN N, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB[J]. Applied and Environmental Microbiology, 2006, 72(7): 5069-5072. DOI:10.1128/aem.03006-05.

    [24] MAIDAK B L, COLE J R, LILBURN T G, et al. The RDP (ribosomal database project) continues[J]. Nucleic Acids Research, 2000, 28(1): 173-174. DOI:10.1093/nar/28.1.173.

    [25] MONIKA B, HUSON D H. SILVA, RDP, Greengenes, NCBI and OTT- how do these taxonomies compare?[J]. BMC Genomics, 2017, 18(Suppl 2): 114. DOI:10.1186/s12864-017-3501-4.

    [26] PRICE MORGAN N, DEHAL PARAMVIR S, ARKIN ADAM P. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix[J]. Molecular Biology and Evolution, 2009, 26(7): 1641-1650. DOI:10.1093/molbev/msp077.

    [27] WILKINSON T J, HUWS S A, EDWARDS J E, et al. CowPI: a rumen microbiome focussed version of the PICRUSt functional inference software[J]. Frontiers in Microbiology, 2018, 9(5): 1095. DOI:10.3389/fmicb.2018.01095.

    [28] 姚國(guó)強(qiáng). 傳統(tǒng)發(fā)酵乳中細(xì)菌多樣性及其功能基因研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2017.

    [29] XIE Mengxi, AN Feiyu, ZHAO Yue, et al. Metagenomic analysis of bacterial community structure and functions during the fermentation of da-jiang, a Chinese traditional fermented food[J]. LWT-Food Science and Technology, 2020, 30(4): 109450. DOI:10.1016/j.lwt.2020.109450.

    [30] STADLMAYR A, FENG Qiang, HUBER-SCH?NAUER U, et al. Gut microbiome development along the colorectal adenoma-carcinoma sequence[J]. Nature Communications, 2015, 6: 6528. DOI:10.1038/ncomms7528.

    [31] 張?bào)蓿?孟祥晨. 乳酸菌的脅迫應(yīng)答及其對(duì)碳水化合物代謝的影響[J].?中國(guó)食品學(xué)報(bào), 2017, 17(6): 145-151. DOI:10.16429/j.1009-7848.2017.06.020

    猜你喜歡
    高通量測(cè)序臘腸
    臘腸香味就是年味
    新年的臘腸
    有趣的臘腸樹
    臘腸果化學(xué)成分及其生物活性
    有趣的臘腸樹
    99国产极品粉嫩在线观看| 欧美人与性动交α欧美精品济南到| 女人高潮潮喷娇喘18禁视频| 亚洲人成网站在线播放欧美日韩| 大香蕉久久成人网| 国产三级在线视频| 九色亚洲精品在线播放| 麻豆一二三区av精品| 九色国产91popny在线| 中文字幕av电影在线播放| 亚洲精品国产色婷婷电影| 男人的好看免费观看在线视频 | 欧美成人a在线观看| 久久6这里有精品| 午夜亚洲福利在线播放| 丰满人妻一区二区三区视频av| 3wmmmm亚洲av在线观看| 亚洲美女视频黄频| 国产午夜精品久久久久久一区二区三区 | 99久久精品一区二区三区| 亚洲国产精品成人综合色| 国产不卡一卡二| 赤兔流量卡办理| 亚洲久久久久久中文字幕| 一夜夜www| av福利片在线观看| 国产av一区在线观看免费| 一区二区三区四区激情视频 | a级毛片免费高清观看在线播放| 国产黄片美女视频| 久久久精品大字幕| 麻豆成人午夜福利视频| 精品乱码久久久久久99久播| 91午夜精品亚洲一区二区三区 | 九九热线精品视视频播放| 我的女老师完整版在线观看| 久久久久久久久大av| 国产三级在线视频| 91狼人影院| 蜜桃久久精品国产亚洲av| 成人无遮挡网站| 成年女人永久免费观看视频| 少妇人妻一区二区三区视频| 免费观看在线日韩| 亚洲,欧美,日韩| 一级av片app| 国产亚洲精品久久久com| 少妇猛男粗大的猛烈进出视频 | 99久久精品热视频| 最新中文字幕久久久久| 搡老岳熟女国产| 欧美日韩精品成人综合77777| 日本-黄色视频高清免费观看| 身体一侧抽搐| 久久香蕉精品热| 亚洲中文字幕日韩| 热99re8久久精品国产| 久久精品久久久久久噜噜老黄 | 波多野结衣巨乳人妻| 色视频www国产| 亚洲美女黄片视频| 九九在线视频观看精品| 国内精品一区二区在线观看| 欧洲精品卡2卡3卡4卡5卡区| av在线天堂中文字幕| 日韩精品有码人妻一区| 亚洲av二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 午夜精品在线福利| 狂野欧美激情性xxxx在线观看| 日韩精品有码人妻一区| 国产一区二区三区视频了| 波多野结衣巨乳人妻| 精品久久久久久成人av| 春色校园在线视频观看| 欧美日韩中文字幕国产精品一区二区三区| 婷婷六月久久综合丁香| 亚洲经典国产精华液单| 免费人成在线观看视频色| 久久久久久久久中文| a级一级毛片免费在线观看| 国产色婷婷99| 亚洲欧美日韩高清在线视频| 十八禁国产超污无遮挡网站| 国产在视频线在精品| 免费看日本二区| x7x7x7水蜜桃| 最近在线观看免费完整版| 亚洲黑人精品在线| 久久久久久久午夜电影| 一进一出好大好爽视频| 一级a爱片免费观看的视频| 精品一区二区三区av网在线观看| 成人性生交大片免费视频hd| 欧美成人一区二区免费高清观看| 最近在线观看免费完整版| 美女 人体艺术 gogo| 男人的好看免费观看在线视频| 国产欧美日韩精品亚洲av| 黄色女人牲交| www日本黄色视频网| 亚州av有码| 麻豆国产97在线/欧美| 我要看日韩黄色一级片| 俄罗斯特黄特色一大片| 亚洲熟妇熟女久久| 亚洲精品色激情综合| 精品久久久久久久久亚洲 | 69人妻影院| 少妇猛男粗大的猛烈进出视频 | bbb黄色大片| 特大巨黑吊av在线直播| 99热这里只有精品一区| 亚洲成人久久爱视频| 桃色一区二区三区在线观看| 欧美丝袜亚洲另类 | 日韩在线高清观看一区二区三区 | 午夜久久久久精精品| 搡老妇女老女人老熟妇| 黄色女人牲交| 欧美日韩黄片免| 啦啦啦观看免费观看视频高清| 综合色av麻豆| 两性午夜刺激爽爽歪歪视频在线观看| 久久九九热精品免费| 嫩草影视91久久| 免费看日本二区| 十八禁网站免费在线| 尾随美女入室| 国产精品一区www在线观看 | 中国美女看黄片| 久久精品国产亚洲av香蕉五月| 欧美日韩精品成人综合77777| 久久精品国产自在天天线| 久久这里只有精品中国| 国产亚洲精品久久久久久毛片| 成年女人永久免费观看视频| 精品一区二区免费观看| 男人舔女人下体高潮全视频| 少妇丰满av| 深夜a级毛片| 亚洲 国产 在线| 身体一侧抽搐| 精品人妻偷拍中文字幕| 精品不卡国产一区二区三区| 午夜视频国产福利| 国产午夜精品论理片| 三级男女做爰猛烈吃奶摸视频| 免费在线观看影片大全网站| 国产免费一级a男人的天堂| 国产高清视频在线播放一区| 97碰自拍视频| 日韩欧美精品v在线| 午夜福利欧美成人| 88av欧美| 18禁黄网站禁片免费观看直播| 亚洲美女黄片视频| 长腿黑丝高跟| 亚洲精品影视一区二区三区av| 51国产日韩欧美| 国产成人av教育| 简卡轻食公司| 一边摸一边抽搐一进一小说| 日日干狠狠操夜夜爽| 国产免费av片在线观看野外av| 亚洲中文日韩欧美视频| 国内精品久久久久精免费| 亚洲欧美日韩高清在线视频| av专区在线播放| 国产探花在线观看一区二区| 可以在线观看毛片的网站| 亚洲欧美日韩高清在线视频| 夜夜夜夜夜久久久久| 久9热在线精品视频| 桃红色精品国产亚洲av| 超碰av人人做人人爽久久| 综合色av麻豆| 国产精品日韩av在线免费观看| 99在线人妻在线中文字幕| 一级毛片久久久久久久久女| 亚洲av中文av极速乱 | 国产精品日韩av在线免费观看| 我要看日韩黄色一级片| 嫩草影院入口| 久9热在线精品视频| 国产精品久久久久久av不卡| 一区二区三区免费毛片| 五月玫瑰六月丁香| 成人特级黄色片久久久久久久| 免费看日本二区| 一本一本综合久久| 九色成人免费人妻av| 免费在线观看影片大全网站| 联通29元200g的流量卡| 男人和女人高潮做爰伦理| 天堂√8在线中文| 日日啪夜夜撸| 国产黄片美女视频| 亚洲人成网站在线播| 久久亚洲真实| 少妇的逼好多水| 国产精品福利在线免费观看| .国产精品久久| 国产一区二区三区视频了| 3wmmmm亚洲av在线观看| 日日啪夜夜撸| 国产高清激情床上av| 国产欧美日韩精品亚洲av| 中亚洲国语对白在线视频| 久久久久性生活片| 亚洲欧美日韩高清在线视频| 一个人看视频在线观看www免费| 成人特级黄色片久久久久久久| 精品国产三级普通话版| 色综合站精品国产| or卡值多少钱| 国产高清有码在线观看视频| 熟女人妻精品中文字幕| 精品久久久久久久久久久久久| 欧美绝顶高潮抽搐喷水| 熟女电影av网| 午夜免费成人在线视频| 亚洲av免费在线观看| 18禁在线播放成人免费| 少妇人妻精品综合一区二区 | 乱码一卡2卡4卡精品| 尤物成人国产欧美一区二区三区| 欧美成人免费av一区二区三区| 亚洲久久久久久中文字幕| 亚洲精品456在线播放app | 国产伦一二天堂av在线观看| 国产激情偷乱视频一区二区| 精品久久久噜噜| 中文字幕久久专区| 麻豆一二三区av精品| 黄色一级大片看看| 日韩欧美 国产精品| 久久精品国产亚洲av香蕉五月| 亚洲国产欧洲综合997久久,| 性色avwww在线观看| 黄色丝袜av网址大全| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av五月六月丁香网| 我的老师免费观看完整版| 22中文网久久字幕| av中文乱码字幕在线| 亚洲18禁久久av| 在线观看美女被高潮喷水网站| 欧美成人性av电影在线观看| 99精品久久久久人妻精品| 精品人妻一区二区三区麻豆 | 国产精品女同一区二区软件 | 国产大屁股一区二区在线视频| 亚洲人成网站高清观看| 欧美+日韩+精品| 夜夜看夜夜爽夜夜摸| 精品一区二区三区人妻视频| 国产乱人视频| 国产乱人视频| 午夜a级毛片| 中文字幕av成人在线电影| 九九久久精品国产亚洲av麻豆| 深夜精品福利| 高清日韩中文字幕在线| 深夜精品福利| 国产一区二区激情短视频| 欧美性猛交黑人性爽| 日本 欧美在线| 两人在一起打扑克的视频| 日本撒尿小便嘘嘘汇集6| 桃色一区二区三区在线观看| 最近最新中文字幕大全电影3| 国产av麻豆久久久久久久| 日日摸夜夜添夜夜添av毛片 | 久久精品国产亚洲av香蕉五月| 国产精品女同一区二区软件 | 国产色婷婷99| 日本熟妇午夜| av天堂在线播放| avwww免费| 男女边吃奶边做爰视频| 少妇的逼好多水| 日韩av在线大香蕉| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 在线播放无遮挡| 国产精品亚洲一级av第二区| 91久久精品电影网| 五月玫瑰六月丁香| 午夜爱爱视频在线播放| 日韩大尺度精品在线看网址| 网址你懂的国产日韩在线| 黄片wwwwww| 日本免费a在线| 色视频www国产| 又爽又黄无遮挡网站| 国产美女午夜福利| 亚洲图色成人| 日日摸夜夜添夜夜添小说| 中文字幕人妻熟人妻熟丝袜美| 全区人妻精品视频| 人人妻,人人澡人人爽秒播| 亚洲精品国产成人久久av| 狂野欧美白嫩少妇大欣赏| 少妇的逼水好多| av在线亚洲专区| 日韩,欧美,国产一区二区三区 | 国产美女午夜福利| 中文字幕久久专区| 国产精品1区2区在线观看.| 日韩欧美国产一区二区入口| 精品一区二区免费观看| 我要搜黄色片| 麻豆国产av国片精品| 日本在线视频免费播放| 亚洲午夜理论影院| av.在线天堂| 18禁黄网站禁片免费观看直播| 亚洲熟妇中文字幕五十中出| 一个人观看的视频www高清免费观看| 很黄的视频免费| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| 国产91精品成人一区二区三区| 欧美xxxx黑人xx丫x性爽| 一级av片app| 日韩精品中文字幕看吧| 国产真实伦视频高清在线观看 | 高清在线国产一区| 免费在线观看影片大全网站| 国产精品爽爽va在线观看网站| 中文字幕精品亚洲无线码一区| 淫妇啪啪啪对白视频| 精品一区二区三区视频在线观看免费| 日本色播在线视频| 国产国拍精品亚洲av在线观看| 国产免费av片在线观看野外av| 草草在线视频免费看| 亚洲国产精品成人综合色| 欧美中文日本在线观看视频| 欧美日韩国产亚洲二区| av天堂在线播放| a级毛片a级免费在线| 在线免费观看不下载黄p国产 | 成人三级黄色视频| 黄色视频,在线免费观看| x7x7x7水蜜桃| 91在线观看av| 亚洲 国产 在线| 国产精品三级大全| 99久久精品国产国产毛片| 国产精品永久免费网站| 三级男女做爰猛烈吃奶摸视频| 久久久久久伊人网av| 久久久色成人| 日韩大尺度精品在线看网址| 成人av一区二区三区在线看| 国产国拍精品亚洲av在线观看| 乱系列少妇在线播放| 18+在线观看网站| 久久人人精品亚洲av| 日韩av在线大香蕉| 天堂影院成人在线观看| 欧美区成人在线视频| 国产男人的电影天堂91| 久久天躁狠狠躁夜夜2o2o| 九色国产91popny在线| 中文字幕av成人在线电影| 两个人视频免费观看高清| 久久精品国产鲁丝片午夜精品 | 女人十人毛片免费观看3o分钟| 久久久久久久久久久丰满 | 在线免费观看的www视频| bbb黄色大片| 国语自产精品视频在线第100页| 免费电影在线观看免费观看| 亚洲人与动物交配视频| 午夜日韩欧美国产| 色在线成人网| 中亚洲国语对白在线视频| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 淫秽高清视频在线观看| 精品午夜福利视频在线观看一区| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 一进一出抽搐动态| 免费一级毛片在线播放高清视频| 少妇人妻一区二区三区视频| 色精品久久人妻99蜜桃| 小蜜桃在线观看免费完整版高清| 亚洲美女黄片视频| 午夜福利在线观看吧| 真人做人爱边吃奶动态| 美女大奶头视频| 少妇猛男粗大的猛烈进出视频 | 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看| 十八禁网站免费在线| 97热精品久久久久久| 欧美日本视频| 男人舔女人下体高潮全视频| 在线免费观看不下载黄p国产 | 欧美日韩精品成人综合77777| 国产高清有码在线观看视频| 观看免费一级毛片| 91久久精品电影网| 91久久精品国产一区二区成人| 真人做人爱边吃奶动态| 99久久精品热视频| 亚洲欧美激情综合另类| 国内精品美女久久久久久| 久久6这里有精品| 免费人成在线观看视频色| 丰满的人妻完整版| 精品99又大又爽又粗少妇毛片 | 国产亚洲精品久久久com| 内射极品少妇av片p| 男女边吃奶边做爰视频| 亚洲人成网站高清观看| 国产精品av视频在线免费观看| 成人二区视频| 日韩欧美精品v在线| 午夜福利视频1000在线观看| 国产亚洲精品av在线| 日韩精品有码人妻一区| 99久久精品国产国产毛片| 成人欧美大片| 男女啪啪激烈高潮av片| 欧美一区二区亚洲| 成熟少妇高潮喷水视频| 99视频精品全部免费 在线| 99久久中文字幕三级久久日本| 久9热在线精品视频| 中文字幕av在线有码专区| 亚洲成人精品中文字幕电影| 久久6这里有精品| 国产aⅴ精品一区二区三区波| 一区二区三区高清视频在线| 乱人视频在线观看| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 少妇的逼好多水| 婷婷丁香在线五月| 成年女人看的毛片在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品人妻熟女av久视频| 欧美激情国产日韩精品一区| eeuss影院久久| 嫩草影院入口| 中文字幕久久专区| 91午夜精品亚洲一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 国国产精品蜜臀av免费| 午夜久久久久精精品| 神马国产精品三级电影在线观看| 亚洲精华国产精华精| 人人妻人人看人人澡| 国产高清激情床上av| 久久久国产成人精品二区| 亚洲中文字幕日韩| 深爱激情五月婷婷| 午夜精品一区二区三区免费看| 午夜福利18| 久久久色成人| 精品不卡国产一区二区三区| 99久久无色码亚洲精品果冻| 成人高潮视频无遮挡免费网站| 有码 亚洲区| av黄色大香蕉| 亚洲真实伦在线观看| 给我免费播放毛片高清在线观看| 欧美绝顶高潮抽搐喷水| 亚洲国产欧洲综合997久久,| 特大巨黑吊av在线直播| 国产乱人视频| 午夜福利在线观看吧| 国产亚洲精品av在线| 亚洲一级一片aⅴ在线观看| 日本黄色片子视频| 日本-黄色视频高清免费观看| 女人十人毛片免费观看3o分钟| 18+在线观看网站| 国产私拍福利视频在线观看| 男人狂女人下面高潮的视频| 大型黄色视频在线免费观看| 97热精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 男女边吃奶边做爰视频| 国产色爽女视频免费观看| avwww免费| 成年免费大片在线观看| 99久久久亚洲精品蜜臀av| 18禁裸乳无遮挡免费网站照片| 国产精品99久久久久久久久| 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| 日韩欧美国产一区二区入口| 日本一二三区视频观看| 国内精品一区二区在线观看| 国产男人的电影天堂91| 中文在线观看免费www的网站| 黄片wwwwww| 久久精品国产99精品国产亚洲性色| ponron亚洲| 51国产日韩欧美| 1000部很黄的大片| 一区二区三区免费毛片| 一个人看视频在线观看www免费| 成年女人永久免费观看视频| 午夜久久久久精精品| 婷婷六月久久综合丁香| 美女高潮的动态| 嫩草影视91久久| 成人av在线播放网站| 熟女电影av网| 免费人成视频x8x8入口观看| 久久久久久久久中文| 久久久久久久久久黄片| 色av中文字幕| 亚洲美女黄片视频| bbb黄色大片| 999久久久精品免费观看国产| 日韩av在线大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 天堂影院成人在线观看| 亚洲精品成人久久久久久| 桃色一区二区三区在线观看| 国产精品综合久久久久久久免费| 高清毛片免费观看视频网站| 老熟妇乱子伦视频在线观看| 我要搜黄色片| 亚洲成av人片在线播放无| 欧美性感艳星| 窝窝影院91人妻| 男女边吃奶边做爰视频| 免费av观看视频| 九九热线精品视视频播放| 国产探花在线观看一区二区| 国产男人的电影天堂91| 亚洲av成人精品一区久久| 亚洲18禁久久av| 男女那种视频在线观看| 国产免费av片在线观看野外av| 老司机福利观看| 午夜福利在线在线| 18禁黄网站禁片午夜丰满| 成人特级黄色片久久久久久久| www.www免费av| 亚洲av二区三区四区| 国产高清有码在线观看视频| 99久久九九国产精品国产免费| 狠狠狠狠99中文字幕| 欧美激情国产日韩精品一区| 天美传媒精品一区二区| 搞女人的毛片| 亚洲七黄色美女视频| 国产精品99久久久久久久久| 国产精品精品国产色婷婷| 国产精品99久久久久久久久| 嫩草影院精品99| 日日摸夜夜添夜夜添av毛片 | 乱系列少妇在线播放| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 国产乱人视频| 欧美3d第一页| 国产三级在线视频| 欧美3d第一页| 五月伊人婷婷丁香| 99热这里只有精品一区| 黄片wwwwww| 久久久久久久久久久丰满 | 啪啪无遮挡十八禁网站| 麻豆成人av在线观看| a在线观看视频网站| 中文亚洲av片在线观看爽| 午夜免费激情av| 国产免费一级a男人的天堂| 国产美女午夜福利| 精品人妻1区二区| 久久99热6这里只有精品| 精品人妻视频免费看| 极品教师在线视频| 精品人妻视频免费看| 中出人妻视频一区二区| 久久亚洲精品不卡| 国内久久婷婷六月综合欲色啪| 亚洲av.av天堂| 人妻丰满熟妇av一区二区三区| 欧美性猛交黑人性爽| 韩国av一区二区三区四区| 制服丝袜大香蕉在线| 亚洲性久久影院| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 一个人观看的视频www高清免费观看| 又黄又爽又免费观看的视频| 欧美黑人巨大hd| 国产极品精品免费视频能看的| 好男人在线观看高清免费视频| 国产乱人视频| 女人被狂操c到高潮| 美女免费视频网站| 亚洲不卡免费看| 美女黄网站色视频| 午夜精品一区二区三区免费看| 国产黄a三级三级三级人| 看黄色毛片网站| 一级毛片久久久久久久久女| 国产探花极品一区二区| 免费搜索国产男女视频| 日本-黄色视频高清免费观看| 在线免费十八禁|