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

    不同貯藏溫度下調(diào)理牛排菌群動態(tài)變化和相互作用

    2025-04-16 00:00:00史娜邢超宋麗萍毛婷杜建萍李龍耿健強孫曉冬楊麗梅李黎
    肉類研究 2025年4期

    摘 要:為了解不同貯藏溫度下調(diào)理牛排微生物多樣性和細菌之間的相互作用關(guān)系,比較分析冷藏(4 ℃)和脫冷貯藏(15 ℃)牛排菌群結(jié)構(gòu)。將新鮮調(diào)理牛排進行冷藏(4 ℃)和脫冷貯藏(15 ℃),基于傳統(tǒng)的微生物培養(yǎng)計數(shù)和高通量測序技術(shù)分析不同貯藏溫度下過保質(zhì)期牛排菌群動態(tài)變化和相互作用。結(jié)果表明,牛排4 ℃貯藏下的菌落總數(shù)(7.4×108 CFU/g)極顯著低于15 ℃(5.2×109 CFU/g)(P<0.001)。菌群動態(tài)變化結(jié)果表明,冷藏或脫冷貯藏不影響牛排菌群組成,但優(yōu)勢菌屬發(fā)生更替。肉桿菌屬(Carnobacterium)、彎曲乳桿菌屬(Latilactobacillus)及假單胞菌屬(Pseudomonas)是2 種貯藏溫度下的優(yōu)勢腐敗菌屬。索絲菌屬(Brochothrix)僅為4 ℃貯藏下的優(yōu)勢菌屬,沙雷氏菌屬(Serratia)僅為15 ℃貯藏下的優(yōu)勢菌屬。菌屬豐度相關(guān)性網(wǎng)絡(luò)中PageRank>0.4且出現(xiàn)頻率最高的菌屬分別為不動桿菌屬(Acinetobacter)和索絲菌屬。不動桿菌屬與沙雷氏菌屬、假單胞菌屬、索絲菌屬,沙雷氏菌屬與索絲菌屬,索絲菌屬與假單胞菌屬相對豐度在特定貯藏溫度和時間下存在相關(guān)性,這表明牛排腐敗菌之間存在潛在種間相互作用。索絲菌屬、彎曲乳桿菌屬及假單胞菌屬可能是牛排優(yōu)勢腐敗菌。

    關(guān)鍵詞:貯藏溫度;高通量測序技術(shù);微生物多樣性;種間相互作用

    Dynamic Changes and Interactions of Bacterial Communities in Processed Beef Steak at Different Storage Temperatures

    SHI Na1, XING Chao1, SONG Liping1, MAO Ting2, DU Jianping1, LI Long3, GENG Jianqiang1,*, SUN Xiaodong1, YANG Limei1, LI Li1

    (1. Beijing Institute of Food Inspection and Research (Beijing Municipal Center for Food Safety Monitoring and Risk Assessment), Laboratory of Key Technologies of Major Comprehensive Guarantee of Food Safety for State Administration for Market Regulation, Beijing 100094, China; 2. Beijing Inspection and Testing Certification Center, Beijing 101300, China;

    3. Beijing Institute of Metrology and Testing Science, Beijing 100029, China)

    Abstract: In this study, we aimed to understand the bacterial diversity and interactions in processed beef steak at different storage temperatures. Comparative analysis of the bacterial community structure in chilled (4 ℃) and non-chilled (15 ℃) steak was performed, and the dynamic changes and interactions of bacterial communities in spoiled meat during storage at different temperatures after shelf life were analyzed using traditional microbial culture method and high throughput sequencing technology. The results showed that the total number of colonies was significantly lower at 4 ℃ than that at 15 ℃ (7.4 × 108 CFU/g versus 5.2 × 109 CFU/g) (P lt; 0.001). The microbiota in steak was unaffected by storage temperature, but the dominant bacterial genera varied. Carnobacterium, Latilactobacillus and Pseudomonas were the dominant spoilage genera at both storage temperatures, Brochothrix being the only dominant genus at 4 ℃ and Serratia at 15 ℃. The results of microbial community abundance correlation network analysis showed that the genera with the highest frequency of occurrence at

    PageRank gt; 0.4 were Acinetobacter and Brochothrix. There was a correlation between the relative abundance of Acinetobacter and that of Serratia, Pseudomonas, and Brochothrix, between the relative abundance of Serratia and that of Brochothrix, and between the relative abundance of Brochothrix and that of Pseudomonas at a certain storage temperature and time, indicating a potential interaction between steak spoilage bacteria. Brochothrix, Latilactobacillus and Pseudomonas may be dominant spoilage bacteria in steak.

    Keywords: storage temperature; high throughput sequencing technology; bacterial diversity; inter-species interactions

    DOI:10.7506/rlyj1001-8123-20241016-264

    中圖分類號:TS251.52 " " " " " " " " " " " " " " " " " " 文獻標志碼:A 文章編號:1001-8123(2025)04-0030-10

    引文格式:

    史娜, 邢超, 宋麗萍, 等. 不同貯藏溫度下調(diào)理牛排菌群動態(tài)變化和相互作用[J]. 肉類研究, 2025, 39(4): 30-39. DOI:10.7506/rlyj1001-8123-20241016-264. " "http://www.rlyj.net.cn

    SHI Na, XING Chao, SONG Liping, et al. Dynamic changes and interactions of bacterial communities in processed beef steak at different storage temperatures[J]. Meat Research, 2025, 39(4): 30-39. (in Chinese with English abstract) DOI:10.7506/rlyj1001-8123-20241016-264. " "http://www.rlyj.net.cn

    2022年,中國牛肉市場消費量為9 86.1萬 t,位居全球第二。在牛排、生冷鮮牛肉、牛肉卷、牛肉丸等眾多牛肉產(chǎn)品中,國民對牛排的消費熱情持續(xù)升高。據(jù)《2022年牛排消費趨勢洞察》報告顯示,2021年線上牛排消費增速持續(xù)上漲,同比增長7 倍。牛肉中的水分和營養(yǎng)物質(zhì)能夠為微生物提供生長條件,如假單胞菌屬(Pseudomonas)和不動桿菌屬(Acinetobacter)等,這些微生物可導(dǎo)致牛肉腐敗[1-4]。因此,腐敗微生物的生長和代謝活動引起的牛肉腐敗成為肉類工業(yè)的關(guān)注熱點。

    貯藏溫度是影響鮮肉貨架期的主要因素[5]。目前,冷藏與冷凍是最常用的2 種肉類貯藏方式,其對應(yīng)溫度范圍分別為0~10、-18~-20 ℃[6]。盡管大多數(shù)國家已經(jīng)規(guī)定了食品最高貯藏溫度限值,但據(jù)報道,由于配送系統(tǒng)、運輸距離和時間的不穩(wěn)定性,運輸過程中的溫度可能達到10 ℃甚至更高。這種不受控制的溫度條件可能會顯著縮短肉類貨架期[7-8]。而且,一些鄉(xiāng)村市場的肉類仍然在室溫下露天保存,這無疑會增加肉類變質(zhì)風險[9]。因此,有必要更深入地了解不同貯藏溫度下肉品中的微生物組成動態(tài)變化及其相互作用。

    高通量測序技術(shù)和生物信息學(xué)分析技術(shù)為食品中可培養(yǎng)和不可培養(yǎng)細菌的研究提供了高效的分析工具,使全面描述菌群的動態(tài)變化、分析微生物代謝活動與肉類腐敗之間的相關(guān)性成為可能[10-11]。目前,已有關(guān)于采用高通量測序技術(shù)表征肉類腐敗中微生物組成的報道[12],然而肉類菌群組成在不同貯藏溫度下的差異較大。目前,大多數(shù)研究集中在不同冷藏或過冷溫度(包括4、-1、-2、-3 ℃)或不同包裝方式(如氣調(diào)包裝)對菌群的影響[6,13-14],比較脫冷貯藏下肉類菌群組成和種間相互作用的研究較為少見。

    本研究將新鮮調(diào)理牛排分別進行冷藏(4 ℃)和脫冷貯藏(15 ℃),評估不同貯藏溫度與時間對牛排菌群組成的影響,旨在深入了解貯藏溫度對牛排菌群及菌種之間潛在相互作用的影響,確定牛排的關(guān)鍵致腐微生物。

    1 材料與方法

    1.1 材料

    在北京大型商超采買3 個品牌的調(diào)理牛排(F1、F2、F3),共63 個樣品(每種品牌各21 個),樣品出廠日期均為當天,保質(zhì)期2 d。

    1.2 儀器與設(shè)備

    BagMixer 400稀釋用電子秤、DiluFow基本型均質(zhì)器 " 法國Interscience公司;HPX-25085-III恒溫恒濕箱 上海新苗醫(yī)療器械制造有限公司。

    1.3 方法

    1.3.1 樣品制備

    因從保質(zhì)期內(nèi)牛排中鑒定的高豐度細菌無法與腐敗建立有效聯(lián)系,因此本研究以不同溫度下長時間貯藏至過保質(zhì)期的調(diào)理牛排為樣品,分別在冷藏(4 ℃)、脫冷(15 ℃)條件下貯藏,在保質(zhì)期結(jié)束當天、過保質(zhì)期3 d及過保質(zhì)期6 d取樣,進行菌落總數(shù)測定和16S rRNA基因高通量測序。

    1.3.2 菌落總數(shù)測定

    參考GB 4789.2—2022《食品安全國家標準 食品微生物學(xué)檢驗 菌落總數(shù)測定》測定菌落總數(shù)。

    1.3.3 16S rRNA基因高通量測序

    1.3.3.1 PacBio文庫制備和測序

    提取樣品總DNA后,根據(jù)全長引物序列合成帶有Barcode的特異性引物,進行聚合酶鏈式反應(yīng)(polymerase chain reaction,PCR)擴增并對其產(chǎn)物進行純化、定量和均一化,形成的測序文庫(SMRT Bell)質(zhì)檢合格后,采用PacBio Sequel II進行測序。

    1.3.3.2 數(shù)據(jù)預(yù)處理

    將PacBio下機數(shù)據(jù)導(dǎo)出為CCS文件后,使用Lima v1.7.0軟件,通過Barcode對CCS進行識別,得到Raw-CCS序列數(shù)據(jù)。使用Cutadapt 1.9.1軟件對引物序列進行識別與去除,并進行長度過濾,得到不包含引物序列的Clean-CCS序列。使用UCHME v4.2軟件鑒定并去除嵌合體序列,得到Effective-CCS序列。

    1.3.3.3 α-多樣性、主坐標分析(principal coordinate analysis,PCoA)及屬水平網(wǎng)絡(luò)分析α-多樣性、PCoA及屬水平分布分析在BMKCloud中完成。其中,α-多樣性分類水平選擇屬,組間差異檢驗方法選擇Student’s t-test;PCoA分類水平選擇屬;屬水平分布分析選擇相對豐度前10。屬水平豐度相關(guān)性網(wǎng)絡(luò)圖基于R包igraph實現(xiàn)[15],物種數(shù)量選擇前30,相關(guān)性類型選擇Spearman,相關(guān)系數(shù)閾值選擇|0.6|,顯著性水平選擇0.05。采用Gephi v.0.9.2對構(gòu)建網(wǎng)絡(luò)圖的圖密度和PageRank值進行計算,圖密度反映菌屬之間相互作用程度,圖密度越大,表示菌屬之間相互作用越頻繁[16]。

    1.4 數(shù)據(jù)處理

    數(shù)據(jù)表示為平均值±標準差,多組間數(shù)據(jù)比較使用SPSS 24.0中方差分析的鄧肯多重范圍檢驗。使用R v.4.0.2軟件的Hmisc包對屬水平相對豐度的Pearson相關(guān)性進行分析。

    2 結(jié)果與分析

    2.1 調(diào)理牛排菌落總數(shù)

    由圖1a可知,4 ℃貯藏下,菌落總數(shù)為7.4×108 CFU/g,顯著低于15 ℃貯藏下菌落總數(shù)(5.2×109 CFU/g),脫冷貯藏下調(diào)理牛排的菌落總數(shù)極顯著升高(P<0.001),加速了其腐敗變質(zhì)進程。

    由圖1b~d可知,F(xiàn)1在出廠當天菌落總數(shù)為4.3×106 CFU/g。4 ℃和15 ℃貯藏下,在保質(zhì)期當天,F(xiàn)1菌落總數(shù)分別為4.3×108 CFU/g和1.1×109 CFU/g;在過保質(zhì)期3 d,F(xiàn)1菌落總數(shù)分別為1.5×109 CFU/g和1.0×109 CFU/g;在過保質(zhì)期6 d,F(xiàn)1菌落總數(shù)分別為6.1×108 CFU/g和1.3×109 CFU/g。F2在出廠當天菌落總數(shù)為1.0×106 CFU/g。4 ℃和15 ℃貯藏下,在保質(zhì)期當天,F(xiàn)2菌落總數(shù)分別為5.0×106 CFU/g和2.5×108 CFU/g;

    在過保質(zhì)期3 d,F(xiàn)2菌落總數(shù)分別為7.8×106 CFU/g和4.9×109 CFU/g;在過保質(zhì)期6 d,F(xiàn)2菌落總數(shù)分別為8.5×108 CFU/g和7.9×109 CFU/g。F3在出廠當天菌落總數(shù)為2.3×107 CFU/g。4 ℃和15 ℃貯藏下,在保質(zhì)期當天,F(xiàn)3菌落總數(shù)分別為2.3×108 CFU/g和1.2×1010 CFU/g;

    在過保質(zhì)期3 d,F(xiàn)3菌落總數(shù)分別為1.0×109 CFU/g和9.1×109 CFU/g;在過保質(zhì)期6 d,F(xiàn)3菌落總數(shù)分別為1.7×109 CFU/g和5.9×109 CFU/g。目前,我國尚未針對生肉制定統(tǒng)一的菌落總數(shù)限量指標。歐盟、美國及日本規(guī)定每克牛肉需氧平板計數(shù)約在5.0×106 CFU/g以下,視為合格。由上述結(jié)果可知,本研究所取調(diào)理牛排F1和F2菌落總數(shù)在出廠當天均小于該限值,F(xiàn)3菌落總數(shù)則大于該限值。在保質(zhì)期當天,僅4 ℃貯藏的F2菌落總數(shù)符合該限值標準,F(xiàn)1和F3菌落總數(shù)均大于該限值。15 ℃貯藏下的牛排細菌生長更快,可加速牛排腐敗變質(zhì)。

    2.2 調(diào)理牛排菌群α-多樣性

    α-多樣性反映牛排菌群物種豐富度和多樣性[11]。Shannon指數(shù)用以描述菌群多樣性。Chao1指數(shù)用以描述菌群物種豐富度,其值越高,菌群物種豐富度越高。由表1、2可知,不同貯藏溫度和時間下,F(xiàn)1、F2、F3菌群多樣性均未發(fā)生顯著改變,這表明初始菌群中的部分細菌參與了牛排的腐敗變質(zhì)過程,多種適應(yīng)性強的細菌在牛排腐敗過程依然保持生長與繁殖。值得注意的是,F(xiàn)2菌群豐富度在一定程度上受到貯藏時間和溫度的影響,但與對照組相比差異不顯著。

    2.3 調(diào)理牛排菌群屬水平相對豐度

    本研究選擇相對豐度排名前10的菌屬(占總菌屬豐度的90%以上,為優(yōu)勢菌屬)進行相對豐度分析。由圖2a可知,在出廠當天,F(xiàn)1中相對豐度排名前3的優(yōu)勢菌屬分別為嗜冷桿菌屬(Psychrobacter)(52.4%)、肉桿菌屬(Carnobacterium)(15.3%)、索絲菌屬(Brochothrix)(13.5%);F2中相對豐度排名前3的優(yōu)勢菌屬分別為彎曲乳桿菌屬(Latilactobacillus)(34.8%)、肉桿菌屬(29.0%)、索絲菌屬(10.0%);F3中相對豐度排名前3的優(yōu)勢菌屬分別為肉桿菌屬(46.2%)、彎曲乳桿菌屬(20.1%)、索絲菌屬(12.6%)。其中,嗜冷桿菌屬、彎曲乳桿菌屬及肉桿菌屬分別為3 種牛排出廠時相對豐度最高的優(yōu)勢菌屬。肉桿菌屬和索絲菌屬在3 種牛排中均為優(yōu)勢菌屬。Liang Ce等[6]對不同貯藏溫度下新鮮羊肉菌群變化的研究顯示,冷藏(4 ℃)和過冷(-1.5 ℃)貯藏下索絲菌屬為優(yōu)勢腐敗菌屬,嗜冷桿菌屬則為超冷(-4 ℃)和超冷凍

    (-9 ℃)貯藏下的優(yōu)勢腐敗菌屬。Fang Jinyu等[17]研究表明,Brochothrix thermosphata BT27為4 ℃好氧貯藏條件下牛肉的優(yōu)勢腐敗菌。由此可見,本研究中出廠當天牛排已存在報道的腐敗菌屬。

    a.出廠當天;b、c.分別為4 ℃和15 ℃貯藏下保質(zhì)期當天;d、e.分別為4 ℃和15 ℃貯藏下過保質(zhì)期3 d;f、g.分別為4 ℃和15 ℃下過保質(zhì)期6 d。圖3、4、6、7同。

    圖 2 不同貯藏溫度和時間下調(diào)理牛排優(yōu)勢菌屬分布

    Fig. 2 Distribution of dominant bacterial genera in processed beef steak at different temperatures and times

    由圖2b、c可知,在保質(zhì)期當天,與4 ℃貯藏相比,脫冷貯藏并不影響F1、F2與F3的菌群組成,但不同貯藏溫度下牛排優(yōu)勢菌屬存在差異。具體地,4 ℃貯藏F1中相對豐度最高的優(yōu)勢菌屬為索絲菌屬(34.1%),F(xiàn)2、F3中相對豐度最高的優(yōu)勢菌屬均為彎曲乳桿菌屬(F2:56.5%,F(xiàn)3:61.1%);15 ℃貯藏下F1、F3中相對豐度最高的優(yōu)勢菌屬均為彎曲乳桿菌屬(F1:41.8%,F(xiàn)3:74.9%),F(xiàn)2中相對豐度最高的優(yōu)勢菌屬為肉桿菌屬(63.2%)。

    由圖2d、e可知,在過保質(zhì)期3 d時,4 ℃和15 ℃貯藏下牛排均已徹底腐敗。4 ℃貯藏F1、F2、F3中相對豐度最高的優(yōu)勢菌屬分別為索絲菌屬(51.5%)、肉桿菌屬(67.0%)、彎曲乳桿菌屬(43.8%),這與出廠當天、保質(zhì)期當天牛排的優(yōu)勢菌屬類似。15 ℃貯藏F1、F2、F3中相對豐度最高的優(yōu)勢菌屬分別為彎曲乳桿菌屬(40.2%)、假單胞菌屬(78.3%)、沙雷氏菌屬(Serratia)(30.2%),表明過保質(zhì)期后15 ℃貯藏牛排的優(yōu)勢腐敗菌發(fā)生改變。

    由圖2f、g可知,過保質(zhì)期6 d時,隨著腐敗程度的增加,4 ℃貯藏F1、F3中相對豐度最高的優(yōu)勢菌屬均為索絲菌屬(F1:54.8%,F(xiàn)3:32.1%),F(xiàn)2中相對豐度最高的優(yōu)勢菌屬為假單胞菌屬(86.8%);15 ℃貯藏F1、F3中相對豐度最高的則為彎曲乳桿菌屬(F1:41.3%,F(xiàn)3:33.9%),F(xiàn)2中相對豐度最高的優(yōu)勢菌屬仍為假單胞菌屬(45.9%)。

    對比可知,4 ℃貯藏時,F(xiàn)1中的索絲菌屬豐度逐漸升高,而15 ℃貯藏時,F(xiàn)1中的彎曲乳桿菌屬則始終為相對豐度最高的優(yōu)勢菌屬;F2中相對豐度最高的優(yōu)勢菌屬的出現(xiàn)受冷藏影響呈現(xiàn)滯后性,即假單胞菌屬在過保質(zhì)期6 d時成為相對豐度最高的優(yōu)勢菌屬,而在15 ℃貯藏時這一現(xiàn)象過保質(zhì)期3 d便已出現(xiàn)。F3中優(yōu)勢菌屬在2 種貯藏條件下呈現(xiàn)類似的動態(tài)演替,即在保質(zhì)期當天出現(xiàn)相對豐度最高的優(yōu)勢菌屬,隨著貯藏時間的延長,更多細菌出現(xiàn)相對豐度增高的現(xiàn)象。肉類涉及分切、運輸、銷售等諸多環(huán)節(jié),可能導(dǎo)致樣品中不同細菌豐度存在差異[18-20]。然而,本研究所抽取的F1和F3中的菌群結(jié)構(gòu)在4 ℃和15 ℃貯藏下比較穩(wěn)定,優(yōu)勢菌屬的相對豐度變化明顯,未引入新的優(yōu)勢菌屬。這可能是由于從生產(chǎn)到銷售過程中F1和F3處于較為穩(wěn)定的環(huán)境,而F2在貯藏過程中,優(yōu)勢菌屬發(fā)生改變,可能是由細菌種間的相互作用引起的生態(tài)位競爭所導(dǎo)致。

    2.4 調(diào)理牛排菌群PCoA

    由圖3可知,在出廠當天,3 個品牌牛排樣品的置信區(qū)間相交,表明此時各樣品菌群組成相似。隨著貯藏時間的延長,4 ℃和15 ℃下,F(xiàn)1、F2、F3牛排樣品各自聚類,表明3 個品牌的牛排菌群組成出現(xiàn)差異。

    2.5 調(diào)理牛排菌群屬水平豐度相關(guān)性網(wǎng)絡(luò)

    為探究腐敗過程中低豐度菌屬與高豐度細菌屬之間的相關(guān)性,對F1、F2、F3菌屬相對豐度進行相關(guān)性分析,以相對豐度前30的菌屬構(gòu)建豐度相關(guān)性網(wǎng)絡(luò),見圖4。

    由圖5可知,保質(zhì)期當天、過保質(zhì)期3 d與過保質(zhì)期6 d牛排在4 ℃和15 ℃下的菌屬豐度相關(guān)性網(wǎng)絡(luò)圖的圖密度均高于出廠當天,表明隨著貯藏時間的延長,細菌之間豐度變化更為復(fù)雜。此外,4 ℃貯藏下,相關(guān)性網(wǎng)絡(luò)圖的圖密度明顯大于15 ℃,這表明4 ℃貯藏牛排中細菌之間的相互作用更加頻繁與復(fù)雜。

    性[21-23]。由圖6可知,PageRank值排名前5的菌屬中,不動桿菌屬和索絲菌屬出現(xiàn)頻率最高,共在5 個網(wǎng)絡(luò)圖中出現(xiàn)。不動桿菌屬在4 ℃下的2 個網(wǎng)絡(luò)圖中出現(xiàn)。雖然不動桿菌屬并不是優(yōu)勢菌屬,但相關(guān)性網(wǎng)絡(luò)分析結(jié)果提示其重要性,強調(diào)非優(yōu)勢菌屬在牛排腐敗進程中也具有關(guān)鍵作用。已有研究[24-27]表明,不動桿菌屬能夠造成冷藏期間南美白對蝦、金槍魚、冰鮮雞肉及水煮小龍蝦等多種食品基質(zhì)腐敗。

    微生物群落并非由眾多微生物簡單拼湊而成,而是不同物種間通過資源競爭、代謝互補、信號交流和基因水平轉(zhuǎn)移等相互作用形成的高度復(fù)雜的生態(tài)系統(tǒng)。微生物群落在地球化學(xué)元素馴化、農(nóng)業(yè)生產(chǎn)和人類健康等領(lǐng)域承擔著重要功能[28-29]。在食品領(lǐng)域,微生物平衡體系與食品腐敗密切相關(guān)。Andreevskaya等[30]探究3 種嗜冷乳酸菌(Leuconostoc gelidum、Lactococcus piscium和Lactobacillus oligofermentans)在肉類腐敗中的相互作用發(fā)現(xiàn),3 種嗜冷乳酸菌的種間相互作用強烈影響碳水化合物分解代謝/運輸、發(fā)酵、能量產(chǎn)生和翻譯相關(guān)基因表達。Joishy等[31]通過構(gòu)建生凝乳(由生乳制備的凝乳)、煮沸乳與凝乳的核心細菌共生網(wǎng)絡(luò)發(fā)現(xiàn),凝乳中Lactobacillus豐度具有優(yōu)勢,且與Enterococcus、Streptococcus及Leuconostoc豐度呈負相關(guān),推測與其內(nèi)部存在生態(tài)位競爭有關(guān)。對與相關(guān)性網(wǎng)絡(luò)圖中關(guān)鍵節(jié)點菌屬相對豐度呈現(xiàn)正/負相關(guān)的菌屬進行統(tǒng)計,由圖7可知,優(yōu)勢腐敗菌屬相對豐度和關(guān)鍵節(jié)點菌屬相對豐度之間存在相關(guān)性。例如,出廠當天牛排菌群相關(guān)性網(wǎng)絡(luò)中,不動桿菌屬與沙雷氏菌屬相對豐度相關(guān)性為0.72,保質(zhì)期當天(4 ℃)的不動桿菌屬與假單胞菌屬、索絲菌屬相對豐度相關(guān)性分別為0.75、0.78,沙雷氏菌屬和索絲菌屬相v對豐度相關(guān)性為0.75。4 ℃貯藏下,過保質(zhì)期3 d時不動桿菌屬與索絲菌屬相對豐度相關(guān)性為0.83,過保質(zhì)期6 d時不動桿菌屬與假單胞菌屬、索絲菌屬的相關(guān)性分別為-0.87、0.80;15 ℃貯藏下,過保質(zhì)期3 d時索絲菌屬與假單胞菌屬相對豐度相關(guān)性為-0.83。這表明致腐細菌在牛排腐敗過程中可能存在相互作用。值得注意的是,4 ℃貯藏下,保質(zhì)期當天不動桿菌和假單胞菌相對豐度呈現(xiàn)正相關(guān)性,而在過保質(zhì)期后,不動桿菌和假單胞菌相對豐度呈現(xiàn)負相關(guān),提示致腐菌之間可能也存在相互競爭。

    3 結(jié) 論

    索絲菌屬、假單胞菌屬、沙雷氏菌屬、肉桿菌屬及彎曲乳桿菌屬是牛排在不同貯藏溫度與時間下的優(yōu)勢菌屬。其中,肉桿菌屬、彎曲乳桿菌屬及假單胞菌屬為4、15 ℃下的優(yōu)勢菌屬。而索絲菌屬僅在4 ℃下成為優(yōu)勢菌屬、沙雷氏菌屬僅在15 ℃下成為優(yōu)勢菌屬,表明不同貯藏溫度下牛排腐敗菌群組成存在差異。不動桿菌屬與沙雷氏菌屬、假單胞菌屬、索絲菌屬,沙雷氏菌屬與索絲菌屬,索絲菌屬與假單胞菌屬相對豐度在特定貯藏溫度和時間下存在相關(guān)性,表明牛排腐敗菌之間存在潛在相互作用,這一結(jié)果提供了肉類腐敗菌種間相互作用的高通量測序證據(jù)。本研究可為牛排質(zhì)量控制中腐敗菌評估體系的完善提供理論基礎(chǔ)。

    參考文獻:

    [1] Ercolini D, Russo F, Nasi A, et al. Mesophilic and psychrotrophic bacteria from meat and their spoilage potential in vitro and in beef[J]. Applied and Environmental Microbiology, 2009, 75(7): 1990-2001. DOI:10.1128/AEM.02762-08.

    [2] ZHU Y L, WANG W, LI M, et al. Microbial diversity of meat products under spoilage and its controlling approaches[J]. Frontiers in Nutrition, 2022, 9: 1078201. DOI:10.3389/fnut.2022.1078201.

    [3] 方金玉. 冷鮮牛肉中假單胞菌和熱死環(huán)絲菌致腐機制及其光動力殺菌控制研究[D]. 杭州: 浙江工商大學(xué), 2022. DOI:10.27462/d.cnki.ghzhc.2022.000876.

    [4] 鄧森榮, 向蓉, 魏琦麟, 等. 冷鮮肉微生物種類及微生物預(yù)測模型的研究進展[J]. 廣東農(nóng)業(yè)科學(xué), 2021, 48(3): 133-141. DOI:10.16768/j.issn.1004-874X.2021.03.016.

    [5] NYCHAS G J, SKANDAMIS P N, TASSOU C C, et al. Meat spoilage during distribution[J]. Meat Science, 2008, 78(1/2): 77-89. DOI:10.1016/j.meatsci.2007.06.020.

    [6] LIANG C, ZHANG D Q, ZHENG X C, et al. Effects of different storage temperatures on the physicochemical properties and bacterial community structure of fresh lamb meat[J]. Food Science of Animal Resources, 2021, 41(3): 509-526. DOI:10.5851/kosfa.2021.e15.

    [7] DOULGERAKI A I, ERCOLINI D, VILLANI F, et al. Spoilage microbiota associated to the storage of raw meat in different conditions[J]. International Journal of Food Microbiology, 2012, 157(2): 130-141. DOI:10.1016/j.ijfoodmicro.2012.05.020.

    [8] GIARRATANA F, NALBONE L, ZIINO G, et al. Temperature fluctuations along food supply chain: a dynamic and stochastic predictive approach to establish the best temperature value in challenge tests for Listeria monocytogenes[J]. Italian Journal of Food Safety, 2022, 11(1): 9981. DOI:10.4081/ijfs.2022.9981.

    [9] MIAO Z M, LI S, QIN K, et al. Prevalence and antimicrobial resistance of Salmonella isolates recovered from retail pork in major village markets in Tai’an region, China[J]. Journal of Food Protection, 2017, 80(10): 1635-1640. DOI:10.4315/0362-028X.JFP-17-019.

    [10] JOHNSON J S, SPAKOWICZ D J, HONG B Y, et al. Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis[J]. Nature Communications, 2019, 10(1): 5029. DOI:10.1038/s41467-019-13036-1.

    [11] ZHANG T Y, LI H Z, MA S L, et al. The newest oxford nanopore r10.4.1 full-length 16S rRNA sequencing enables the accurate resolution of species-level microbial community profiling[J]. Appiled and Environmental Microbiology, 2023, 89(10): e0060523. DOI:10.1128/aem.00605-23.

    [12] YANG X Y, ZHU L X, ZHANG Y M, et al. Microbial community dynamics analysis by high-throughput sequencing in chilled beef longissimus steaks packaged under modified atmospheres[J]. Meat Science, 2018, 141: 94-102. DOI:10.1016/j.meatsci.2018.03.010.

    [13] YANG J, YANG X Y, LIANG R R, et al. The response of bacterial communities to carbon dioxide in high-oxygen modified atmosphere packaged beef steaks during chilled storage[J]. Food Research International, 2022, 151: 110872. DOI:10.1016/j.foodres.2021.110872.

    [14] TOOMIK E, ROOD L, BOWMAN J P, et al. Microbial spoilage mechanisms of vacuum-packed lamb meat: a review[J]. International Journal of Food Mcirobiology, 2023, 387: 110056. DOI:10.1016/j.ijfoodmicro.2022.110056.

    [15] LIU D H, ZHANG R T, S CHEN, et al. Analysis of gastric microbiome reveals three distinctive microbial communities associated with the occurrence of gastric cancer[J]. BMC Microbiology, 2022, 22: 184. DOI:10.1186/s12866-022-02594-y.

    [16] JACOMY M, VENTURINI T, HEYMANN S, et al. ForceAtlas2,

    a continuous graph layout algorithm for handy network visualization designed for the Gephi software[J]. PLoS ONE, 2014, 9: e98679. DOI:10.1371/journal.pone.0098679.

    [17] FANG J Y, FENG L F, LU H X, et al. Metabolomics reveals spoilage characteristics and interaction of Pseudomonas lundensis and Brochothrix thermosphacta in refrigerated beef[J]. Food Research International, 2022, 156: 111139. DOI:10.1371/journal.pone.0098679.

    [18] VAN BA H, SEO H W, SEONG P N, et al. The fates of microbial populations on pig carcasses during slaughtering process, on retail cuts after slaughter, and intervention efficiency of lactic acid spraying[J]. International Journal of Food Microbiology, 2019, 294: 10-17. DOI:10.1016/j.ijfoodmicro.2019.01.015.

    [19] PARK J, BAE D, KIM S A. Microbial trace investigation throughout the entire chicken supply chain based on metagenomic high-throughput sequencing[J]. Food Research International, 2023, 169: 112775. DOI:10.1016/j.foodres.2023.112775.

    [20] CAMPOS CALERO G, CABALLERO GOMEZ N, et al. In silico mapping of microbial communities and stress responses in a porcine slaughterhouse and pork products through its production chain, and the efficacy of HLE disinfectant[J]. Food Research International, 2020, 136: 109486. DOI:10.1016/j.foodres.2020.109486.

    [21] TANG H, SHI R X, HE T S, et al. Tensorflow solver for quantum pagerank in large-scale networks[J]. Science Bulletin, 2021, 66(2): 120-126. DOI:10.1016/j.scib.2020.09.009.

    [22] GONG Y H. Hedgerank: heterogeneity-aware, energy-efficient partitioning of personalized PageRank at the edge[J]. Micromachines, 2023, 14(9): 1714. DOI:10.3390/mi14091714.

    [23] MA K, LI J W, ZHAO M Y, et al. PPRTGI: a personalized pagerank graph neural network for tf-target gene interaction detection[J]. IEEE/ACM Transactions on Computational Biology and Bioinformatics, 2024, 21(3): 480-491. DOI:10.1109/TCBB.2024.3374430.

    [24] XIA J Y, JIANG N, ZHANG B, et al. Bacterial changes in boiled crayfish between different storage periods and characterizations of the specific spoilage bacteria[J]. Foods, 2023, 12(16): 3006. DOI:10.3390/foods12163006.

    [25] ZHU S Q, WU H H, ZHANG C L, et al. Spoilage of refrigerated Litopenaeus vannamei: eavesdropping on Acinetobacter acyl-homoserine lactones promotes the spoilage potential of Shewanella baltica[J]. Journal of Food Science and Technology, 2018, 55(5): 1903-1912. DOI:10.1007/s13197-018-3108-z.

    [26] WANG X Y, XIE J. Growth kinetics and spoilage potential of co-culturing Acinetobacter johnsonii and Pseudomonas fluorescens from bigeye tuna (Thunnus obesus) during refrigerated storage[J]. Current Microbiology, 2020, 77(8): 1637-1646. DOI:10.1007/s00284-020-01978-5.

    [27] ZHANG T, DING H, CHEN L, et al. Characterization of chilled chicken spoilage using an integrated microbiome and metabolomics analysis[J]. Food Research International, 2021, 144: 110328. DOI:10.1016/j.foodres.2021.110328.

    [28] FALKOWSKI P G, FENCHEL T, DELONG E F. The microbial engines that drive earth’s biogeochemical cycles[J]. Science, 2008, 320: 1034-1039. DOI:10.1126/science.1153213.

    [29] CHO I, BLASER M J. The human microbiome: at the interface of health and disease[J]. Nature Reviews Genetics, 2012, 13: 260-270. DOI:10.1038/nrg3182.

    [30] ANDREEVSKAYA M, JAASKELAINEN E, et al. Food spoilage-associated Leuconostoc, Lactococcus and Lactobacillus species display different survival strategies in response to competition[J]. Applied and Environmental Microbiology, 2018, 84: e00554-18. DOI:10.1128/AEM.00554-18.

    [31] JOISHY T K, DEHINGIA M, KHAN M R. Bacterial diversity and metabolite profiles of curd prepared by natural fermentation of raw milk and back sloping of boiled milk[J]. World Journal of Microbiology amp; Biotechnology, 2019, 35: 102. DOI:10.1007/s11274-019-2677-y.

    精品第一国产精品| 欧美老熟妇乱子伦牲交| 动漫黄色视频在线观看| 久久香蕉精品热| 超碰成人久久| 脱女人内裤的视频| 亚洲 欧美一区二区三区| 淫妇啪啪啪对白视频| 窝窝影院91人妻| 99国产精品99久久久久| 欧美人与性动交α欧美软件| 亚洲色图av天堂| 日本黄色视频三级网站网址 | 久久天躁狠狠躁夜夜2o2o| 欧美成狂野欧美在线观看| 无人区码免费观看不卡| 精品久久蜜臀av无| 多毛熟女@视频| 精品欧美一区二区三区在线| 又大又爽又粗| 十八禁人妻一区二区| 18在线观看网站| 久久久国产成人免费| 久久久久久亚洲精品国产蜜桃av| 国产精品永久免费网站| 国产精品电影一区二区三区 | 老司机影院毛片| cao死你这个sao货| 国产欧美日韩一区二区精品| bbb黄色大片| 99久久综合精品五月天人人| 黄色a级毛片大全视频| 黄色毛片三级朝国网站| 看免费av毛片| 日韩欧美在线二视频 | 中文字幕精品免费在线观看视频| 香蕉久久夜色| 亚洲精品一卡2卡三卡4卡5卡| 午夜视频精品福利| 一区二区日韩欧美中文字幕| 国产精品免费一区二区三区在线 | 视频区图区小说| 欧美另类亚洲清纯唯美| 国产野战对白在线观看| 国产不卡av网站在线观看| 老司机影院毛片| 青草久久国产| av福利片在线| 怎么达到女性高潮| 香蕉久久夜色| 99riav亚洲国产免费| 免费看a级黄色片| 久久国产精品影院| 又黄又爽又免费观看的视频| 日日爽夜夜爽网站| 亚洲中文日韩欧美视频| 热99国产精品久久久久久7| 国产亚洲精品久久久久久毛片 | 久久人妻熟女aⅴ| 91国产中文字幕| 中文字幕最新亚洲高清| 老汉色∧v一级毛片| 一本一本久久a久久精品综合妖精| 两性夫妻黄色片| 两个人看的免费小视频| 91字幕亚洲| 99re6热这里在线精品视频| 国产精品成人在线| 美女福利国产在线| 日韩成人在线观看一区二区三区| 最近最新免费中文字幕在线| 人人妻人人添人人爽欧美一区卜| www.自偷自拍.com| 日本a在线网址| 91麻豆av在线| 亚洲三区欧美一区| 大香蕉久久网| 久久精品国产99精品国产亚洲性色 | 男女床上黄色一级片免费看| 三级毛片av免费| 午夜日韩欧美国产| ponron亚洲| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 黄色成人免费大全| 国产精品国产av在线观看| 久久精品91无色码中文字幕| 精品少妇久久久久久888优播| 老司机靠b影院| 黄色视频不卡| 久久久国产成人精品二区 | 国产xxxxx性猛交| 久久亚洲真实| 久久人妻av系列| 婷婷成人精品国产| 精品国产一区二区三区久久久樱花| 最新在线观看一区二区三区| 在线观看66精品国产| 欧美激情 高清一区二区三区| 中文欧美无线码| 日日爽夜夜爽网站| 一a级毛片在线观看| 日韩人妻精品一区2区三区| 国内久久婷婷六月综合欲色啪| 麻豆乱淫一区二区| 宅男免费午夜| 欧美日韩福利视频一区二区| 国产精品乱码一区二三区的特点 | 极品少妇高潮喷水抽搐| 黄片小视频在线播放| 国产蜜桃级精品一区二区三区 | 欧美国产精品va在线观看不卡| 国产精品电影一区二区三区 | 国产亚洲精品第一综合不卡| 免费黄频网站在线观看国产| 欧美av亚洲av综合av国产av| 国产高清videossex| av天堂久久9| 中文字幕人妻熟女乱码| 99精品久久久久人妻精品| 99精品久久久久人妻精品| bbb黄色大片| 亚洲精品国产色婷婷电影| 亚洲欧美一区二区三区久久| 丝袜在线中文字幕| 一级a爱片免费观看的视频| 满18在线观看网站| 精品久久久久久久久久免费视频 | 又紧又爽又黄一区二区| 国产精品.久久久| e午夜精品久久久久久久| 欧美精品啪啪一区二区三区| 99精品欧美一区二区三区四区| 99热国产这里只有精品6| 亚洲欧美一区二区三区久久| 不卡一级毛片| 精品少妇一区二区三区视频日本电影| 亚洲国产精品一区二区三区在线| 极品教师在线免费播放| 男女下面插进去视频免费观看| 超色免费av| 色尼玛亚洲综合影院| 一本大道久久a久久精品| 天天躁夜夜躁狠狠躁躁| 久久精品国产清高在天天线| 日韩有码中文字幕| 亚洲专区字幕在线| 热99久久久久精品小说推荐| x7x7x7水蜜桃| 黑人巨大精品欧美一区二区mp4| 天堂动漫精品| 大片电影免费在线观看免费| 国产一区二区三区视频了| 国产成人精品在线电影| 91麻豆av在线| 日韩中文字幕欧美一区二区| 国产有黄有色有爽视频| 多毛熟女@视频| 久99久视频精品免费| 91精品国产国语对白视频| 无遮挡黄片免费观看| 亚洲五月天丁香| 黄色a级毛片大全视频| av天堂在线播放| 久久久久久久精品吃奶| 一区二区三区国产精品乱码| 他把我摸到了高潮在线观看| 色婷婷av一区二区三区视频| 日本黄色视频三级网站网址 | 欧美最黄视频在线播放免费 | 欧美av亚洲av综合av国产av| 久9热在线精品视频| 天天躁日日躁夜夜躁夜夜| 欧美黑人精品巨大| 操美女的视频在线观看| 中文亚洲av片在线观看爽 | 又黄又爽又免费观看的视频| av有码第一页| 国产精品 欧美亚洲| 久久久精品区二区三区| 18禁美女被吸乳视频| 成人18禁高潮啪啪吃奶动态图| 91字幕亚洲| 天天操日日干夜夜撸| 夜夜躁狠狠躁天天躁| 午夜免费成人在线视频| 一区二区三区精品91| 黑人操中国人逼视频| 精品视频人人做人人爽| 午夜精品久久久久久毛片777| 午夜福利影视在线免费观看| 99re6热这里在线精品视频| 日本wwww免费看| 亚洲色图综合在线观看| 777久久人妻少妇嫩草av网站| 欧美亚洲日本最大视频资源| 露出奶头的视频| av网站在线播放免费| 999久久久精品免费观看国产| 成在线人永久免费视频| 亚洲在线自拍视频| 日韩精品免费视频一区二区三区| 色婷婷av一区二区三区视频| 免费在线观看影片大全网站| 最近最新免费中文字幕在线| 午夜91福利影院| 国产蜜桃级精品一区二区三区 | 精品视频人人做人人爽| 亚洲欧美激情综合另类| 热re99久久精品国产66热6| 久久中文字幕人妻熟女| 日日爽夜夜爽网站| a级片在线免费高清观看视频| 极品人妻少妇av视频| 一级a爱视频在线免费观看| 国产午夜精品久久久久久| 人妻久久中文字幕网| 黄色 视频免费看| 18禁美女被吸乳视频| 精品少妇久久久久久888优播| 一级a爱片免费观看的视频| 后天国语完整版免费观看| 国产免费男女视频| 成人亚洲精品一区在线观看| 丰满迷人的少妇在线观看| 村上凉子中文字幕在线| 欧美日韩一级在线毛片| 午夜两性在线视频| 精品免费久久久久久久清纯 | 免费观看人在逋| 9色porny在线观看| 一区福利在线观看| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 中文字幕另类日韩欧美亚洲嫩草| 一区二区三区国产精品乱码| 亚洲av成人av| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 操出白浆在线播放| 亚洲精品在线美女| 精品福利永久在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲熟妇熟女久久| 欧美中文综合在线视频| 精品亚洲成国产av| 女警被强在线播放| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 99精品在免费线老司机午夜| a级片在线免费高清观看视频| 亚洲一区高清亚洲精品| av线在线观看网站| 成人手机av| 久久精品亚洲熟妇少妇任你| 成在线人永久免费视频| 在线观看免费高清a一片| 另类亚洲欧美激情| 在线av久久热| 国产精品影院久久| 日本一区二区免费在线视频| 老汉色∧v一级毛片| 一级片免费观看大全| 91九色精品人成在线观看| 99热只有精品国产| 在线播放国产精品三级| avwww免费| 免费不卡黄色视频| 亚洲全国av大片| 又大又爽又粗| 亚洲色图综合在线观看| 国产视频一区二区在线看| 国产亚洲欧美98| 国产成人免费观看mmmm| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 日韩视频一区二区在线观看| av片东京热男人的天堂| 午夜91福利影院| 不卡一级毛片| 免费在线观看完整版高清| 亚洲精品乱久久久久久| 丝袜在线中文字幕| 免费高清在线观看日韩| 最新在线观看一区二区三区| 又紧又爽又黄一区二区| 一a级毛片在线观看| 18禁美女被吸乳视频| 高清av免费在线| 美女国产高潮福利片在线看| 国产精品免费一区二区三区在线 | www.熟女人妻精品国产| 日韩免费av在线播放| 香蕉丝袜av| 国产精品美女特级片免费视频播放器 | 色婷婷av一区二区三区视频| 在线十欧美十亚洲十日本专区| 午夜福利乱码中文字幕| 午夜精品在线福利| 精品卡一卡二卡四卡免费| 一级a爱视频在线免费观看| 天堂动漫精品| 久久人人爽av亚洲精品天堂| 欧美午夜高清在线| 免费av中文字幕在线| 在线观看一区二区三区激情| 在线av久久热| 女人高潮潮喷娇喘18禁视频| 老司机午夜十八禁免费视频| 丝瓜视频免费看黄片| 日本a在线网址| 日韩欧美在线二视频 | 男女之事视频高清在线观看| 国产精品久久电影中文字幕 | 国产亚洲av高清不卡| 国产一区二区三区综合在线观看| 操美女的视频在线观看| 午夜激情av网站| 国产精品九九99| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| 免费在线观看亚洲国产| av天堂在线播放| 国产成人av教育| 12—13女人毛片做爰片一| 欧美日韩成人在线一区二区| 亚洲五月色婷婷综合| 久久99一区二区三区| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 少妇猛男粗大的猛烈进出视频| 欧美成人午夜精品| 露出奶头的视频| 亚洲欧洲精品一区二区精品久久久| av天堂在线播放| 国产黄色免费在线视频| 中文字幕高清在线视频| 国产一卡二卡三卡精品| 首页视频小说图片口味搜索| 国产一区二区三区视频了| 国产成人精品无人区| 亚洲中文日韩欧美视频| av在线播放免费不卡| 久久精品亚洲av国产电影网| 真人做人爱边吃奶动态| 成年动漫av网址| 一夜夜www| 久久草成人影院| 国产视频一区二区在线看| 一区福利在线观看| 在线观看免费日韩欧美大片| 亚洲一区二区三区不卡视频| 国产日韩欧美亚洲二区| 一区二区三区精品91| 黑人操中国人逼视频| av网站在线播放免费| 免费在线观看视频国产中文字幕亚洲| 日本精品一区二区三区蜜桃| 99riav亚洲国产免费| 欧美精品av麻豆av| 国产亚洲av高清不卡| 另类亚洲欧美激情| 精品无人区乱码1区二区| 老司机在亚洲福利影院| 美女高潮到喷水免费观看| 亚洲avbb在线观看| 成人av一区二区三区在线看| 一个人免费在线观看的高清视频| 99精品在免费线老司机午夜| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利欧美成人| 国产精品欧美亚洲77777| 两个人免费观看高清视频| 欧洲精品卡2卡3卡4卡5卡区| 一区福利在线观看| 中文字幕最新亚洲高清| 天堂动漫精品| 真人做人爱边吃奶动态| 俄罗斯特黄特色一大片| 老司机深夜福利视频在线观看| 性色av乱码一区二区三区2| 无限看片的www在线观看| 午夜亚洲福利在线播放| 激情在线观看视频在线高清 | 12—13女人毛片做爰片一| 久久久国产成人免费| 国产高清视频在线播放一区| 在线观看日韩欧美| 亚洲伊人色综图| 免费在线观看亚洲国产| 国产在线观看jvid| 午夜久久久在线观看| 久久精品国产亚洲av香蕉五月 | 国产精品九九99| 真人做人爱边吃奶动态| 久久中文字幕一级| 亚洲午夜精品一区,二区,三区| 久久久久久久久久久久大奶| 日本黄色视频三级网站网址 | 日韩成人在线观看一区二区三区| 校园春色视频在线观看| 他把我摸到了高潮在线观看| 人人妻人人添人人爽欧美一区卜| 女人久久www免费人成看片| 一级毛片高清免费大全| 国产免费av片在线观看野外av| 99久久精品国产亚洲精品| a级毛片黄视频| 一区福利在线观看| 青草久久国产| 热99国产精品久久久久久7| 黄色视频不卡| 欧美人与性动交α欧美精品济南到| 香蕉久久夜色| 国产一区二区三区综合在线观看| 欧美中文综合在线视频| 国产野战对白在线观看| 99热国产这里只有精品6| 侵犯人妻中文字幕一二三四区| 久久中文字幕一级| 久久青草综合色| av欧美777| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 亚洲精品中文字幕在线视频| 99国产精品免费福利视频| 青草久久国产| 麻豆成人av在线观看| 露出奶头的视频| 亚洲国产欧美网| 国产亚洲欧美在线一区二区| 国产人伦9x9x在线观看| 国产真人三级小视频在线观看| 日本一区二区免费在线视频| 国产成人欧美| 国产色视频综合| 一级黄色大片毛片| 国产精品久久久久久人妻精品电影| 午夜91福利影院| 天天添夜夜摸| 亚洲午夜理论影院| 视频区欧美日本亚洲| 啦啦啦视频在线资源免费观看| 黄色视频不卡| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频| 亚洲 欧美一区二区三区| 国产不卡av网站在线观看| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 大型av网站在线播放| 伦理电影免费视频| 久久狼人影院| 在线av久久热| 嫩草影视91久久| 成人手机av| 侵犯人妻中文字幕一二三四区| 国产成人精品无人区| 一本大道久久a久久精品| 久久久国产精品麻豆| 久久青草综合色| www.999成人在线观看| 国内久久婷婷六月综合欲色啪| av欧美777| 国产精品一区二区免费欧美| videosex国产| 欧美激情久久久久久爽电影 | 人人妻人人爽人人添夜夜欢视频| 丝袜美足系列| 亚洲av欧美aⅴ国产| 久久久久久免费高清国产稀缺| 午夜老司机福利片| 国产淫语在线视频| xxxhd国产人妻xxx| 欧美日韩乱码在线| 国产精品亚洲一级av第二区| 国产成人av教育| 99re在线观看精品视频| 我的亚洲天堂| 免费观看人在逋| 亚洲精品成人av观看孕妇| 69精品国产乱码久久久| 亚洲欧美色中文字幕在线| 黑人猛操日本美女一级片| www.999成人在线观看| 精品久久久精品久久久| 国产精品久久视频播放| 999久久久国产精品视频| 久久久久国产一级毛片高清牌| 亚洲精品一二三| 中文字幕av电影在线播放| 一夜夜www| 亚洲中文字幕日韩| 在线观看日韩欧美| 午夜福利欧美成人| 看免费av毛片| 狠狠狠狠99中文字幕| 午夜两性在线视频| 99热网站在线观看| 欧美午夜高清在线| 精品国产国语对白av| av超薄肉色丝袜交足视频| 999精品在线视频| 最新在线观看一区二区三区| 1024香蕉在线观看| 精品国产国语对白av| 亚洲国产欧美一区二区综合| 999久久久国产精品视频| 成人三级做爰电影| 99re6热这里在线精品视频| 美国免费a级毛片| 最近最新免费中文字幕在线| 成人国产一区最新在线观看| 国产精品久久电影中文字幕 | 国产高清videossex| 视频区欧美日本亚洲| 欧美不卡视频在线免费观看 | 91精品国产国语对白视频| 国产成人精品久久二区二区91| 国产精品久久久久久精品古装| 国产精品亚洲av一区麻豆| 国产一卡二卡三卡精品| 王馨瑶露胸无遮挡在线观看| 中亚洲国语对白在线视频| 天天躁夜夜躁狠狠躁躁| 免费久久久久久久精品成人欧美视频| 日本精品一区二区三区蜜桃| 高清欧美精品videossex| 91国产中文字幕| 精品亚洲成国产av| 久久国产精品人妻蜜桃| 99精国产麻豆久久婷婷| 亚洲欧美日韩另类电影网站| 黄色片一级片一级黄色片| 一区二区日韩欧美中文字幕| 91精品国产国语对白视频| 90打野战视频偷拍视频| 91九色精品人成在线观看| 黄片小视频在线播放| 91av网站免费观看| 美女福利国产在线| 亚洲人成77777在线视频| 欧美日韩国产mv在线观看视频| 美女高潮到喷水免费观看| 夜夜躁狠狠躁天天躁| 亚洲精品国产精品久久久不卡| 国产99久久九九免费精品| 国产精品二区激情视频| av国产精品久久久久影院| 亚洲精品美女久久久久99蜜臀| 亚洲色图综合在线观看| 国产免费现黄频在线看| 91九色精品人成在线观看| 午夜亚洲福利在线播放| 色老头精品视频在线观看| 1024视频免费在线观看| 丰满饥渴人妻一区二区三| 欧美人与性动交α欧美精品济南到| 欧美日韩亚洲高清精品| 首页视频小说图片口味搜索| 亚洲久久久国产精品| 亚洲国产精品合色在线| 久99久视频精品免费| 亚洲片人在线观看| 亚洲av片天天在线观看| 嫁个100分男人电影在线观看| 免费黄频网站在线观看国产| 精品国产乱子伦一区二区三区| 国产成人精品无人区| 亚洲全国av大片| 国产高清videossex| 国产精品av久久久久免费| 亚洲三区欧美一区| 亚洲第一欧美日韩一区二区三区| 亚洲精品美女久久av网站| 午夜激情av网站| aaaaa片日本免费| 人人妻人人爽人人添夜夜欢视频| 在线十欧美十亚洲十日本专区| 亚洲人成电影免费在线| 在线永久观看黄色视频| 美国免费a级毛片| 国产欧美日韩一区二区精品| 精品一品国产午夜福利视频| 国产成人精品久久二区二区免费| 欧美日韩成人在线一区二区| 免费少妇av软件| 国产1区2区3区精品| 最新的欧美精品一区二区| 国产一区有黄有色的免费视频| 国产精品电影一区二区三区 | 欧美日韩国产mv在线观看视频| 高清黄色对白视频在线免费看| 一边摸一边做爽爽视频免费| 可以免费在线观看a视频的电影网站| 国产欧美亚洲国产| 老司机深夜福利视频在线观看| 国产黄色免费在线视频| 精品电影一区二区在线| 久久ye,这里只有精品| 一夜夜www| 成人国产一区最新在线观看| 精品第一国产精品| 交换朋友夫妻互换小说| 女人高潮潮喷娇喘18禁视频| 精品少妇久久久久久888优播| 国产精品 国内视频| 亚洲少妇的诱惑av| 欧美最黄视频在线播放免费 | 国产成人av教育| 欧美黄色片欧美黄色片| 老熟妇乱子伦视频在线观看| 色婷婷av一区二区三区视频| 日本黄色日本黄色录像| 少妇 在线观看|