基于低壓變頻電場技術對帶魚(Trichiurus lepturus)保鮮過程中微生物群落影響分析*

史恬恬 謝 超① 張家瑋 李遠慧 朱亞猛 吳玉婷 周卓穎 虞 舟

基于低壓變頻電場技術對帶魚()保鮮過程中微生物群落影響分析*

史恬恬1謝 超1①張家瑋1李遠慧1朱亞猛1吳玉婷1周卓穎1虞 舟2

(1. 浙江海洋大學食品與藥學學院 浙江舟山 316022; 2. 舟山匯豐冷藏物流發(fā)展有限公司 浙江舟山 316002)

帶魚()在低溫條件運輸和貯藏過程中, 易受微生物和內(nèi)源性蛋白酶的影響而產(chǎn)生胺類和醛類物質導致腐敗。為了研究低壓變頻電場對微凍帶魚在貯藏過程中的微生物群落特征, 采用Miseq基因組測序技術, 對貯藏至40 d的舟山帶魚的微生物群落進行堿基信息測序。將產(chǎn)生的堿基序列通過OTU聚類、多樣性分析、相關性分析等手段分析了影響帶魚保鮮效果的微生物群落組成和發(fā)育信息。在此基礎上, 研究了不同強度電場(0, 2, 2.5, 3 kV/m)處理對帶魚微生物群落組成的影響。根據(jù)菌屬相對豐度結果, 推測希瓦氏菌屬()、假單胞菌屬()、紫色桿菌屬()是帶魚微生物腐敗主要致腐菌。結果表明, 40 d后, 經(jīng)3 kV/m處理的樣品其OTU6嗜冷菌屬和OTU13假單胞菌屬c的相對豐度均小于1%, 而且OTU7、OTU14、OTU8、OTU9等菌屬未被OTU聚類。結果表明, 經(jīng)3 kV/m電場強度的低壓靜電場處理的抑菌保鮮效果最好, 表現(xiàn)為在貯藏過程中雜菌種類最少, 樣本微生物群落物種豐度和多樣性均最低, 腐敗程度最低。

帶魚(); 低壓變頻電場; 微生物; 高通量測序

帶魚()又名牙帶魚, 其風味鮮美、營養(yǎng)物質種類豐富, 蛋白質高, 易被人體吸收, 深受大眾喜愛(楊勝平, 2010)。帶魚肉質銀白細嫩, 富含人體所必需的多種營養(yǎng)元素, 而且還可以預防血壓較高、動脈狹窄、癌癥等疾病的發(fā)生(賀瑩, 2019)。通常來說帶魚采用低溫條件進行運輸和貯藏, 并且在這期間易受微生物和內(nèi)源性蛋白酶的影響, 易產(chǎn)生一些胺類和醛類物質, 導致魚肉發(fā)生腐敗, 使帶魚體表顏色發(fā)生變化, 引起品質劣變, 而導致魚肉腐敗的一部分原因就是由于微生物的作用(沈妮, 2019)。

在微生物腐敗菌的分離培養(yǎng)鑒定等方法中如實時定量熒光聚合酶鏈式反應、多重聚合酶鏈式反應等在檢測過程中會受到多重因素的干擾, 容易造成錯誤判定, 而高通量測序技術能夠更完美、準確反映樣品中微生物菌群分布, 可以同時完成對測定微生物基因分子的序列測定, 更精準、全面地鑒定樣品中微生物的單一或全面基因組(Di Bella, 2013)。除此之外, Miseq技術在測序過程中無需大量PCR擴增, 實現(xiàn)準確、快速、全面的測序及多樣性分析(徐鵬昊, 2016; 張國林等, 2021)。而一些優(yōu)勢腐敗菌如希瓦氏菌屬、嗜冷菌屬等是傳統(tǒng)分離鑒定方法難以對群落結構進行全面分析的。

低壓靜電場保鮮是一種物理保鮮方法, 通過控制水分子偶極振動減少水的相變時間, 冰晶粗糙度均勻并重新排列形成的小冰晶, 使蛋白質結構得以完整保留, 食物的口感和咀嚼性得到提升, 其次就是降低細胞呼吸速率(Xie, 2021), 外部電場還可以改變水分子與活性酶的結合態(tài), 導致酶的活性降低甚至沒有活性, 在電場反應下, 空氣被電離產(chǎn)生臭氧和負離子, 臭氧能直接氧化并破壞細菌的細胞壁和細胞質膜, 然后進入細胞并作用于其DNA (丹陽等, 2004; Xanthakis, 2013)。

因此, 本研究基于高通量測序技術, 研究低壓靜電場對貯藏至40 d帶魚保鮮過程中微生物組成差異情況, 旨在為保鮮中菌群組成及低壓靜電保鮮技術提供基礎信息(毛俊龍等, 2020)。

1 材料與方法

1.1 試驗材料

新鮮舟山帶魚()采購自舟山沈家門漁港, 體長約0.8～1.0 m, 厚約2～3 cm, 帶魚通體銀白錚亮, 魚眼小且透明, 魚皮濕滑, 無破損品質完好, 經(jīng)過處理后進行包裝袋包裝, 設定溫度為-4 °C, 分別進行電場保存。

1.2 試驗試劑

本試驗所用試劑包括AS132 DNA提取試劑盒、Qubit3.0 DNA試劑盒、Hieff NGS? DNA分選磁株等, 具體見表1所示。

1.3 儀器與設備

本試驗所用儀器與設備包括恒溫振蕩器、電場板、PCR檢測儀等, 具體見表2所示。

1.4 方法

1.4.1 樣品處理 取同一批次的新鮮帶魚為試驗原料, 運送至實驗室后進行清洗、分裝。

對照組(CK): 將帶魚放入無電場的環(huán)境下, 微凍–4 °C進行保鮮試驗, 標記樣品為DY1。

表1 實驗試劑

表2 儀器與設備

實驗組(LVEF): 將帶魚樣品放入低壓靜電場, 微凍–4 °C下進行持續(xù)三種電場強度(2, 2.5, 3 kV/m)處理, 標記樣品為(DY2, DY3, DY4), 電場頻率為50 Hz。

1.4.2 帶魚微生物的提取 將貯藏至40 d的四組帶魚樣品進行研磨攪碎后取10 g帶魚樣品, 加入9倍體積的生理鹽水, 之后進行離心分離(12 500 r/min, 10 min), 取菌體沉淀物進行微生物菌群研究, 每組樣品進行三次平行實驗(Hu, 2015; Li, 2020)。

1.4.3 總DNA的提取 參考廖娟等(2020)的方法。使用無菌槍頭刮取表面菌落, 轉移適量樣品至2 mL樣品管中進行提取。Qubit定量檢測DNA樣本濃度。DY1、DY2、DY3、DY4點樣濃度分別為48.6、51.2、53.6、66.2 ng/μL, 上樣量400 ng, DNA提取電泳如圖1所示。

1.4.4 PCR擴增 將PCR管置于PCR儀中進行擴大范圍遞增, 第一輪PCR作用條件: 94 °C預變性3 min; (94 °C變性30 s → 45 °C退火20 s → 65 °C延伸30 s) 5個循環(huán); (94 °C預變性20 s → 55 °C退火20 s →72 °C延伸30 s) 20個循環(huán); 72 °C再延伸5 min, 10 °C保存。第二輪引入橋式PCR兼容引物, 27F、1492R為通用引物, 將PCR管置于PCR儀中進行擴大, PCR反應條件: 95 °C預變性3 min; (94 °C變性20 s→ 55 °C退火20 s → 72 °C延伸30 s) 5個循環(huán); 72 °C再延伸5 min, 10 °C保存。PCR擴增電泳如圖2所示(張皖君, 2018)。

圖1 DNA提取電泳圖

注: 數(shù)據(jù)單位為bp; 從左至右分別為DY1、DY2、DY3、DY4

圖2 PCR擴增電泳圖

注: 從左至右分別為DY1、DY2、DY3、DY4

1.4.5 MiSeq測序 為了得到均勻的長簇效果和高質量的測序數(shù)據(jù), 使用Qubit 4.0熒光定量系統(tǒng)進行文庫濃度的定量測定。文庫經(jīng)完整性查驗后送往上海生工生物公司進行Illumina MiSeq測序(郭倩倩等, 2019)。

1.4.6 菌群生物信息學分析 采用cutadapt、PEAR、PRINSEQ、RDP classifier、GraPhlAn等軟件分析微生物群落生物信息。

2 結果與討論

2.1 OTU聚類

MiSeq測序過程中會產(chǎn)生數(shù)以萬計的16S序列, 逐一進行物種標注耗時太長, 且第二代測序在經(jīng)PCR擴增后會發(fā)生堿基錯配的概率會大大提升, 故引入OTU聚類, 且后續(xù)微生物信息學分析都是OTU聚類結果展開的。OTU (Operational Taxonomic Units)是在遺傳生物學和統(tǒng)計系統(tǒng)學研究中, 為方便分析而人為對微生物菌群分類添加的標志。通常所有樣本去除冗余序列合并, 去除沒有重復的單序列, 相似度95%以下認為不同屬, 相似度97%以下可以認為測序微生物不同種, 一般采用Usearch作OTU聚類(Edgar, 2013)。

由表3可知, 帶魚樣本的有效序列信息, DY1、DY2、DY3、DY4的堿基信息分別是CGTATC、AGCAGT、GAGGAA、CGAGTC, 有效序列條數(shù)分別為88 714、84 447、83 544、74 711, 堿基數(shù)分別為38 137 408、36 304 277、35 919 079、32 124 163, 序列平均長度大致相同。差異說明電場強度與檢出有效序列條數(shù)、堿基數(shù)呈負相關。

表4列出了各電場所有優(yōu)化序列參數(shù)對應的OTU代表的序列閾值數(shù), 選出與對應代表序列相似性在97%以上的序列OTU名稱, 可以看出各個組間樣品OTU類別和序列數(shù)不盡相同, 需要結合表5的OTU物種注釋分析。OTU物種注釋表明, 帶魚貯藏過程中體內(nèi)微生物群落全部由細菌群組成, 并沒有產(chǎn)生真菌群落。在門水平上, 變形菌門(Proteobacteria)為優(yōu)勢菌門, 僅有OTU3聚類為厚壁菌門(Firmicutes)。在綱水平上, 丙型變形菌綱(Gammaproteobacteria)為優(yōu)勢菌綱, 丙型變形菌綱分布廣泛, 各組樣品均有OTU聚類。在目水平上, 假單胞菌目(Pseudomonadales)聚類水平較高。在科水平表現(xiàn)出了較高的豐度, 聚類結果為希瓦氏菌科(Shewanellaceae)、假單胞菌科(Pseudomonadaceae)、微球菌科(Micrococcaceae)、莫拉氏菌科(Moraxellaceae)、李斯特菌科(Listeriaceae)、柄桿菌科(Caulobacteraceae)、草酸桿菌科(Oxalobacteraceae)、鞘脂桿菌科(Sphingobacteriaceae)。屬水平是Miseq測序OTU分類信息學Usearch分析的最小單位(Caporaso, 2012), 能較客觀地描述該菌的生理特征和作用, 可以用于腐敗原因分析, 本次Miseq測序的OTU聚類屬主要有希瓦氏菌屬()、節(jié)桿菌屬()、柄桿菌屬()、工地桿菌屬()、假單胞菌屬()、索絲菌屬()、嗜冷菌屬()、紫色桿菌屬()、不動桿菌屬() (曹榮等, 2016)。

表3 樣本有效序列統(tǒng)計

表4 OTU序列數(shù)統(tǒng)計表

表5 OTU分類學信息表

2.2 多樣性分析

Miseq測序后的OTU聚類能夠表征微生物群落的具體信息, 但僅靠OTU聚類的分析仍不夠全面, 還需要結合Mothur和R軟件對OTU菌屬做多樣性分析(Lu, 2014)。

表6列出了DY1、DY2、DY3、DY4的Alpha多樣性指數(shù)。在97%的OTU聚類水平上, 微生物豐度指數(shù)Chao和Ace越大, 樣本豐度指數(shù)越高(楊春敏等, 2019); 微生物菌群的多樣性指數(shù)Shannon越大, 樣本多樣性越高, Simple指數(shù)趨勢則與Shannon相反(Gong, 2019)。可以看出四組樣品Chao和Ace豐度指數(shù)大小排序為DY1>DY2>DY3>DY4, Shannon多樣性指數(shù)大小排序與Chao和Ace相同, Simple則相反。樣本文庫覆蓋率均為1, 各樣本測序可信度高。

表6 Alpha多樣性指數(shù)統(tǒng)計表

圖3列出了OTU聚類微生物的屬級相對豐度, 不同顏色區(qū)塊代表各菌屬所占樣本的百分比。DY1、DY2、DY3、DY4中相對豐度最高的是假單胞菌屬a, 分別占了OTU總量的61%、76%、87%、97%。對照組DY1顏色區(qū)塊最多, 有多個獨有OTU聚類菌屬, 其中希瓦氏菌屬(31%)、假單胞菌屬的惡臭假單胞菌(4%)、紫色桿菌屬(1%)是水產(chǎn)品主要致腐菌(曹榮等, 2019)。電場處理組DY2、DY3、DY4也有多個獨有菌屬, 經(jīng)分析均無強致腐性和致病性, 說明電場能夠有效抑制帶魚貯藏過程中腐敗菌的生長, 這與段偉文(2019)的結論相同。DY4微生物群落豐度最小, 表現(xiàn)為嗜冷菌屬(0.1%)、假單胞菌屬c (0%)的相對豐度較少。DY4在40 d時OTU6嗜冷菌屬和OTU13假單胞菌屬c的相對豐度均小于1%, 還有OTU7、OTU14、OTU8、OTU9等菌屬未被OTU聚類。說明3 kV/m電場強度的低壓靜電場抑菌保鮮效果最好, DY2、DY3次之, DY1效果最差。印證了前文中貯藏理化特性、蛋白質特性和水分遷移情況的數(shù)據(jù)。

圖3 菌屬相對豐度柱狀圖

Alpha稀釋性曲線圖變化情況如圖4所示, 橫軸代表多樣性分析抽取的數(shù)據(jù)量, 縱軸代表Alpha指數(shù), 根據(jù)曲線是否能達到平緩來判斷Miseq測序數(shù)據(jù)量的充分性。從曲線趨平程度來看, Alpha稀釋性曲線中DY1>DY2>DY3>DY4, DY4最先趨平, DY3、DY2次之, DY1直至最后仍未趨平。說明DY1貯藏過程中雜菌種類多, 樣本微生物群落物種豐度和多樣性均高于電場組, 腐敗程度也最大。電場處理組DY2、DY3、DY4也有多個獨有菌屬, 但均無致腐性和致病性。

圖4 Alpha指數(shù)稀釋性曲線圖

樣本共線性關系如圖5所示。以門水平為標準, 右弧表征微生物群落物種組成情況, 左弧表征微生物群落的分布比例情況, 一個顏色代表一個物種, 長度代表豐度比例, 第三圈的彩色區(qū)帶一角連接樣本, 另一角連接物種, 區(qū)帶寬度代表具體物種豐度。樣本門中96.82%為變形菌門(Proteobacteria), 僅有2.81%厚壁菌門(Firmicutes), 其他菌門占比均小于1%。

圖5 共線性關系圖

2.3 相關性分析

組間相關性是檢驗實驗可靠性和樣本合理性的重要標準,2越接近1, 表示關聯(lián)性越好。由圖6可知, DY4與DY3、DY2、DY1相關系數(shù)分別為0.72、0.67、0.49, 相關程度逐漸縮小。DY3與DY2、DY1相關系數(shù)分別為0.45、–0.07, 相關性不強。DY2與DY1相關性僅為0.32。DY4、DY3、DY2與對照組DY1主成分差異顯著(<0.05), DY4與DY3樣本最相似(>0.05)。表明了在四組電場中DY4與DY3的關聯(lián)密切程度最好, 而與DY1呈極度不相關。

圖7為PCA主成分分析結果, PCA可以通過數(shù)據(jù)降維運用ANOVA方差分析不同樣本組間差異和距離, 其軸和軸并無實際意義, 百分比表示差異解釋度值, 點越接近, 其菌群組成就越相似。數(shù)據(jù)表明除樣本本身外, 電場處理的DY4、DY3、DY2與未經(jīng)電場處理DY1樣本間點的距離遠, 表明主成分差異顯著(<0.05), DY4與DY3樣本最相似(>0.05), 表明了經(jīng)電場設備處理后對帶魚微生物菌群組成結構有相當顯著的影響。

3 結論

采用Miseq基因組測序技術, 對貯藏至40 d的帶魚微生物群落堿基信息測序, 產(chǎn)生的堿基序列通過OTU聚類、多樣性分析、相關性分析等手段分析影響帶魚保鮮效果的微生物群落組成和發(fā)育信息, 研究電場強度(0, 2, 2.5, 3 kV/m)對帶魚微生物群落組成的影響。OTU聚類屬主要有希瓦氏菌屬()、工地桿菌屬()、節(jié)桿菌屬()等菌屬, Chao和Ace豐度指數(shù)大小排序為DY1>DY2> DY3>DY4。對照組DY1的OTU多樣性高, 有多個獨有OTU聚類菌屬, 電場處理組DY2、DY3、DY4也有多個獨有菌屬, 但均無致腐性和致病性。說明3 kV/m電場強度的低壓靜電場抑菌保鮮效果最好, DY2、DY3次之, DY1效果最差; Alpha稀釋性曲線中DY4最先趨平, DY3、DY2次之, DY1直至最后仍未趨平。DY4與DY3、DY2、DY1相關系數(shù)分別為0.72、0.67、0.49, 相關程度逐漸縮小。綜上說明, 電場能夠有效抑制帶魚貯藏過程中腐敗菌的繁衍生長, 3 kV/m電場強度的抑菌保鮮效果最好, 2、2.5 kV/m次之, 未施加電場效果最不理想。研究結果為低壓靜電保鮮技術提供了基礎信息。

圖6 相關性熱圖

圖7 PCA主成分分析

毛俊龍, 祁雪兒, 張賓, 等, 2020. 基于高通量測序的大黃魚表面微生物群落組成及變化分析[J]. 現(xiàn)代食品科技, 36(12): 27-35, 167.

丹陽, 李里特, 2004. 高壓靜電場(HVEF)臭氧產(chǎn)生能力以及所產(chǎn)生臭氧對毛霉菌的抑制作用[J]. 食品工業(yè)科技, 25(1): 49-51.

楊春敏, 黃建初, 查麗艷, 2019. 海南山蘭米酒酒曲微生物多樣性分析[J]. 食品工程(4): 44-48, 59.

楊勝平, 2010. 帶魚生物保鮮劑及氣調包裝保鮮技術的研究[D]. 上海: 上海海洋大學: 1-64.

沈妮, 2019. 帶魚低溫貯藏蛋白氧化、組織蛋白酶活性及魚肉質地結構的變化規(guī)律[D]. 杭州: 浙江大學: 1-91.

張國林, 景榮先, 劉昆梅, 等, 2021. 新一代測序技術進展及其在藥品質量控制中應用和發(fā)展趨勢[J]. 藥物分析雜志, 41(1): 1-12.

張皖君, 2018. 流化冰處理對鱸魚流通期間品質及微生物菌群變化影響[D]. 上海: 上海海洋大學: 1-71.

段偉文, 2019. 應用低壓靜電場與氣調包裝進行凡納濱對蝦保鮮的研究[D]. 湛江: 廣東海洋大學: 1-63.

賀瑩, 2019. 紫外線殺菌結合氣調包裝技術對帶魚品質的影響[J].肉類研究, 33(1): 37-41.

徐鵬昊, 2016. 高通量測序中拼接問題的研究現(xiàn)狀[J]. 山東農(nóng)業(yè)工程學院學報, 33(1): 42-44.

郭倩倩, 盧彪, 2019. 基于高通量測序不同生產(chǎn)工藝晴隆酸菜細菌多樣性分析[J]. 中國釀造, 38(8): 73-76.

曹榮, 劉淇, 趙玲, 等, 2016. 基于高通量測序的牡蠣冷藏過程中微生物群落分析[J]. 農(nóng)業(yè)工程學報, 32(20): 275-280.

曹榮, 孟輝輝, 劉淇, 等, 2019. 超高壓處理對牡蠣冷藏過程腐敗菌群結構的影響及典型菌株——腐敗希瓦氏菌致死機理研究[J]. 中國食品學報, 19(2): 192-198.

廖娟, 王鋼, 喻世剛, 等, 2020. 利用Illumina MiSeq測序平臺分析沐川烏骨黑雞空腸微生物多樣性[J]. 江蘇農(nóng)業(yè)科學, 48(24): 178-182.

CAPORASO J G, LAUBER C L, WALTERS W A,, 2012. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms [J]. The ISME Journal, 6(8): 1621-1624.

DI BELLA J M, BAO Y G, GLOOR G B,, 2013. High throughput sequencing methods and analysis for microbiome research [J]. Journal of Microbiological Methods, 95(3): 401-414.

EDGAR R C, 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads [J]. Nature Methods, 10(10): 996-998.

GONG X W, LIU C J, LI J,, 2019. Responses of rhizosphere soil properties, enzyme activities and microbial diversity to intercropping patterns on the Loess Plateau of China [J]. Soil and Tillage Research, 195: 104355.

HU D X, PAN L Q, ZHAO Q,, 2015. Transcriptomic response to low salinity stress in gills of the Pacific white shrimp,[J]. Marine Genomics, 24: 297-304.

LI H M, MENG Q L, CHEN T,, 2020. Transcriptomic response to low pH stress in gills of the pacific white shrimp,[J]. Aquaculture Research, 51(1): 175-186.

LU X M, LU P Z, 2014. Characterization of bacterial communities in sediments receiving various wastewater effluents with high-throughput sequencing analysis [J]. Microbial Ecology, 67(3): 612-623.

XANTHAKIS E, HAVET M, CHEVALLIER S,, 2013. Effect of static electric field on ice crystal size reduction during freezing of pork meat [J]. Innovative Food Science & Emerging Technologies, 20: 115-120.

XIE Y, ZHOU K, CHEN B,, 2021. Applying low voltage electrostatic field in the freezing process of beef steak reduced the loss of juiciness and textural properties [J]. Innovative Food Science & Emerging Technologies, 68: 102600.

EFFECT OF LOW-VOLTAGE ALTERNATING FREQUENCY ELECTRIC FIELD ON MICROBIAL COMMUNITY OF HAIRTAILDURING PRESERVATION

SHI Tian-Tian1, XIE Chao1, ZHANG Jia-Wei1, LI Yuan-Hui1, ZHU Ya-Meng1, WU Yu-Ting1, ZHOU Zhuo-Ying1, YU Zhou2

(1. Zhejiang Ocean University College of Food and Medicine, Zhoushan 316022, China; 2. Zhoushan HSBC Cold Storage Logistics Development, Co. Ltd., Zhoushan 316002, China)

To study the effect of low-voltage variable-frequency electric field on the microbial community characteristics of microfrozen hairtail () during storage, Miseq genome sequencing technology was used to sequence the base information of the microbial community of Zhoushan hairtail stored for 40 days. The microbial community composition and development information affecting the preservation effect of hairtail were analyzed by OTU clustering, diversity analysis, and correlation analysis. On this basis, the effects of electric field in different intensities (0, 2, 2.5, 3 kV/m) on the microbial community composition of hairtail were studied. According to the results of relative abundance of bacteria, we speculated that,, andwere the main spoilage bacteria of hairtail microbial spoilage. The results show that after 40-day storage, the relative abundances of OTU6 psychrophilic bacteria and OTU13C were less than 1%, and OTU7, OTU14, OTU8, OTU9 and other bacterial genera were not clustered by OTU. Specifically, the antibacterial and fresh-keeping effect of low-voltage electrostatic field in 3 kV/m intensity was the best, under which the species of miscellaneous bacteria were the least, the species abundance and diversity of microbial community were the lowest, and the degree of corruption was the lowest.

hairtail; low-voltage electrostatic field; microorganism; high-throughput sequencing

TS254

10.11693/hyhz20210800189

*“十三五”國家重點研發(fā)計劃重點專項, 2019YFD0901604號; 舟山市科技計劃項目, 2020C12016號; 舟山市普陀區(qū)科技計劃項目, 2020JH102號。史恬恬, 碩士研究生, E-mail: 2020033759@qq.com

謝 超, 博士, 副教授, E-mail: xc750205@163.com

2021-08-27,

2021-09-25