江楊陽(yáng)，楊水兵，余海霞，姚潔玉，胡亞芹,*

（1.浙江大學(xué)生物系統(tǒng)工程與食品科學(xué)學(xué)院，馥莉食品研究院，智能食品加工技術(shù)與裝備國(guó)家（地方）聯(lián)合實(shí)驗(yàn)室，農(nóng)業(yè)部農(nóng)產(chǎn)品產(chǎn)后處理重點(diǎn)實(shí)驗(yàn)室，農(nóng)業(yè)部農(nóng)產(chǎn)品營(yíng)養(yǎng)功能評(píng)價(jià)實(shí)驗(yàn)室，浙江省農(nóng)產(chǎn)品加工技術(shù)研究重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 310058；2.浙江大學(xué)舟山海洋研究中心，浙江 舟山 316021）

小龍蝦，學(xué)名克氏原鰲蝦（Procambarus clarkii），作為一種重要的經(jīng)濟(jì)蝦類，廣泛分布在我國(guó)長(zhǎng)江中下游地區(qū)。其肉質(zhì)肥美、鮮嫩，倍受消費(fèi)者青睞，現(xiàn)已逐步形成了集“苗種繁育、健康養(yǎng)殖、加工出口、餐飲物流、節(jié)慶文化”于一體的產(chǎn)業(yè)鏈[1]。與其他水產(chǎn)品類似，在運(yùn)輸、貯藏過(guò)程中，極易受微生物活動(dòng)、化學(xué)氧化和酶反應(yīng)3 種破壞機(jī)制的影響，致其品質(zhì)快速劣變[2]。小龍蝦屬于甲殼類水產(chǎn)品，擁有更多細(xì)菌生長(zhǎng)基質(zhì)的自由氨基酸，故更易遭受微生物破壞[3]。為抑制腐敗微生物活動(dòng)，目前應(yīng)用最為廣泛的方法是低溫冷藏。但即使在冷藏過(guò)程中，小龍蝦貨架期仍受到特定腐敗菌的影響，且這類特定菌隨著貯藏時(shí)間延長(zhǎng)，其數(shù)量逐漸占據(jù)優(yōu)勢(shì)地位，故又將這類特定腐敗菌稱為優(yōu)勢(shì)腐敗菌（specific spoilage organism，SSO）。

培養(yǎng)基鑒定法是較為常用的鑒定優(yōu)勢(shì)腐敗菌的方法，主要是通過(guò)在各種選擇性培養(yǎng)基上培養(yǎng)、篩選不同的菌種后，提取這些單菌落DNA，進(jìn)行聚合酶鏈?zhǔn)椒磻?yīng)（polymerase chain reaction，PCR）擴(kuò)增后測(cè)其基因序列，再對(duì)比確定菌株種類?，F(xiàn)已有報(bào)道將該方法應(yīng)用于商業(yè)海洋魚類[4]、冰鮮銀鯧魚[5]、冷卻豬肉[6]的優(yōu)勢(shì)腐敗菌鑒定。但該方法存在較大缺陷，特別是挑選單菌落時(shí)易受主觀因素影響，導(dǎo)致腐敗菌遺漏；此外，部分環(huán)境微生物對(duì)環(huán)境要求苛刻，用該方法分離出來(lái)的微生物僅占微生物總數(shù)的1%～10%[7]。

高通量測(cè)序法（2 代測(cè)序法）不僅可有效將測(cè)序時(shí)間縮短至幾天甚至幾個(gè)小時(shí)，能同時(shí)測(cè)序幾十萬(wàn)至幾百萬(wàn)條DNA分子，還能檢測(cè)出對(duì)環(huán)境要求苛刻及低豐度的微生物，有效彌補(bǔ)傳統(tǒng)培養(yǎng)基鑒定法的不足[8]。目前已廣泛應(yīng)用分析土壤[9]、腸道[10]等的微生物組成?；诖?，該方法應(yīng)用于鑒定水產(chǎn)品腐敗菌相的報(bào)道也逐年增加。Cao Rong等[11]用高通量法鑒定冷藏超高壓滅菌牡蠣的優(yōu)勢(shì)腐敗菌；江艷華等[12-13]將高通量測(cè)序法用于鑒定冷凍南極磷蝦及蝦夷扇貝柱的腐敗菌相組成；鄧曉影等[14]用高通量法鑒定南美白對(duì)蝦腐敗菌相。但其用于鑒定小龍蝦腐敗菌的研究尚鮮見報(bào)道。本實(shí)驗(yàn)應(yīng)用高通量測(cè)序法分別鑒定冷藏小龍蝦的鮮樣、腐敗中期樣及腐敗末期樣所對(duì)應(yīng)的優(yōu)勢(shì)菌群，并將其腐敗末期鑒定結(jié)果與培養(yǎng)基鑒定法進(jìn)行比較，驗(yàn)證高通量測(cè)序法精準(zhǔn)、高效性的同時(shí)，準(zhǔn)確分析小龍蝦物流貯運(yùn)加工過(guò)程中的優(yōu)勢(shì)腐敗菌，旨在為冷藏小龍蝦的腐敗機(jī)理解析及其保鮮研究提供重要的參考依據(jù)。

1 材料與方法

1.1 材料、試劑與培養(yǎng)基

鮮活的小龍蝦樣品（質(zhì)量約40～50 g），購(gòu)于江蘇盱眙。挑選無(wú)病態(tài)、活力足的小龍蝦若干只，將其進(jìn)行去頭處理后放置4 ℃冰箱冷藏至腐敗。

氯化鈉、硼酸、MgO、鹽酸（均為分析純） 國(guó)藥集團(tuán)化學(xué)試劑有限公司；細(xì)菌基因組DNA提取試劑盒天根生化科技有限公司。

大豆酪蛋白瓊脂（TSA）培養(yǎng)基、BP培養(yǎng)基、假單胞菌選擇性CFC瓊脂培養(yǎng)基、結(jié)晶紫膽鹽葡萄糖中性紅瓊脂（VRBGA）培養(yǎng)基、MR瓊脂培養(yǎng)基、孟加拉紅培養(yǎng)基、營(yíng)養(yǎng)肉湯 山東青島海博生物技術(shù)有限公司。

1.2 儀器與設(shè)備

ZDDN-II自動(dòng)定氮儀 浙江托普有限公司；HE-120核酸電泳儀 上海天能有限公司。

1.3 方法

1.3.1 揮發(fā)性鹽基氮（total volatile basic nitrogen，TVB-N）值測(cè)定

參照GB 5009.228—2016《食品中揮發(fā)性鹽基氮的測(cè)定》。以20 g/L硼酸為吸收液，0.01 mol/L鹽酸滴定，精確稱量5 g研碎蝦肉于750 mL的蒸餾管中，加入0.5 g MgO和少許去離子水，連接到定氮儀上測(cè)定。

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

采用?zogul等[15]的方法并略有改動(dòng)。在無(wú)菌條件下稱取10 g蝦肉（取蝦肉肌肉混合樣），絞碎于90 mL滅菌蒸餾水中。充分振蕩、搖勻，制成樣品菌懸液，并進(jìn)行梯度稀釋。選擇合適的稀釋梯度，取1 mL于無(wú)菌平板中，用平板傾注法倒平板。待平板凝固后，倒置放在37 ℃烘箱培養(yǎng)48 h后數(shù)其菌落數(shù)。

1.3.3 單菌落的分離和純化

將小龍蝦放置在4 ℃冷藏至腐敗，參照1.3.2節(jié)方法制得樣品菌懸液，梯度稀釋后，選取合適的稀釋度吸取0.1 mL于多種選擇性培養(yǎng)基平板上均勻涂布，后放置在培養(yǎng)箱中培養(yǎng)24～48 h。挑選各培養(yǎng)基中典型菌落，平板劃線2～3 次，得到純化的單菌落于－80 ℃冰箱中保存。

1.3.4 無(wú)菌蝦肉的制備

參照唐文靜等[16]的方法略有改動(dòng)。將蝦去殼，于無(wú)菌條件下用無(wú)菌水沖洗，用75%乙醇溶液擦拭蝦肉表面，并放置在紫外條件下晾干。

1.3.5 單菌落優(yōu)勢(shì)腐敗菌的篩選

對(duì)各單菌株進(jìn)行致腐能力測(cè)定分析以確定冷藏小龍蝦的優(yōu)勢(shì)腐敗菌，致腐能力的測(cè)定以TVB-N值和菌落總數(shù)為評(píng)定依據(jù)。

將提取保存的菌株活化后制成菌懸液，通過(guò)平板傾注法確定各菌液的濃度，調(diào)整菌液濃度約為106CFU/mL。將制備好的無(wú)菌蝦肉分別浸入各菌液中10 s后取出，放置無(wú)菌平板中，于4 ℃冷藏。分別測(cè)定其在0、3、6、9、12 d的TVB-N值和菌落總數(shù)，以確定優(yōu)勢(shì)腐敗菌。TVB-N值和菌落總數(shù)的測(cè)定方法參照1.3.1節(jié)和1.3.2節(jié)。

1.3.6 單菌落優(yōu)勢(shì)腐敗菌的分子生物學(xué)鑒定

將各優(yōu)勢(shì)腐敗菌于營(yíng)養(yǎng)肉湯中培養(yǎng)12 h，用細(xì)菌基因組DNA提取試劑盒提取各菌株的DNA，作為16S rDNA序列擴(kuò)增反應(yīng)模板。PCR體系為細(xì)菌16S rDNA序列擴(kuò)增模板2 μL、引物27f和1492r各2 μL，2×Master Mix 25 μL，ddH2O 19 μL。PCR程序：95 ℃預(yù)變性3 min，95 ℃變性20 s，56 ℃退火30 s，72 ℃延伸1 min，循環(huán)35 次后，再在72 ℃徹底延伸5 min，4 ℃保存。取10 μL的PCR擴(kuò)增產(chǎn)物進(jìn)行0.8%的瓊脂糖凝膠電泳。將凝膠電泳中條帶清晰，無(wú)明顯拖尾的PCR擴(kuò)增產(chǎn)物送至生工生物工程（上海）股份有限公司測(cè)序。使用NCBI（National Center for Biotechnology Information）中的BLAST對(duì)各菌株序列分析比對(duì)，選取同源性在98%以上的菌株序列，利用MEGA 7.0.26軟件進(jìn)行多重序列對(duì)比，并構(gòu)建Neighbor-Joining系統(tǒng)發(fā)育樹，以確定優(yōu)勢(shì)腐敗菌種類[17]。

1.3.7 高通量測(cè)序法鑒定優(yōu)勢(shì)腐敗菌

選取貯藏期為0、2、8 d的小龍蝦肉，制成樣品菌懸液，直接提取菌懸液的DNA。將提取得到的DNA進(jìn)行PCR擴(kuò)增后，再進(jìn)行瓊脂糖凝膠電泳。將凝膠電泳中條帶清晰，無(wú)明顯拖尾的DNA樣品送至南京諾唯贊生物科技有限公司進(jìn)行高通量測(cè)定分析。

1.4 數(shù)據(jù)及圖像處理

每組實(shí)驗(yàn)進(jìn)行2～3 次平行。數(shù)據(jù)及圖像分別采用Excel和Origin 8.5進(jìn)行處理。

2 結(jié)果與分析

2.1 冷藏期間小龍蝦TVB-N值及TVC值變化

TVB-N值是用作評(píng)判蝦肉腐敗程度的重要指標(biāo)之一，是三甲胺、二甲胺、氨等揮發(fā)性胺類物質(zhì)的總和。這些胺類物質(zhì)是在貯藏期間，由于微生物和內(nèi)源酶的雙重作用下，降解蝦肉蛋白和非蛋白物質(zhì)產(chǎn)生的[18]。根據(jù)GB 10136—2015《動(dòng)物性水品衛(wèi)生標(biāo)準(zhǔn)》規(guī)定，TVB-N值不高于30.00 mg/100 g為合格品。

圖1 冷藏小龍蝦TVB-N值的變化Fig. 1 Change in TVB-N value in red swamp crayfish during cold storage

如圖1所示，貯藏初期，數(shù)值增長(zhǎng)較為平緩，而當(dāng)達(dá)到30 mg/100 g的安全限值時(shí)，增長(zhǎng)顯著增加。小龍蝦TVB-N初始值為5.6 mg/100 g，貯藏第2天數(shù)值為18.76 mg/100 g，而貯藏第8天后，其數(shù)值已達(dá)到45.36 mg/100 g，伴隨較為強(qiáng)烈的氨臭味，可知已達(dá)到小龍蝦的貨架期終點(diǎn)。

圖2 冷藏小龍蝦菌落總數(shù)的變化Fig. 2 Change in total bacterial count in red swamp crayfish during cold storage

如圖2所示，菌落總數(shù)隨貯藏時(shí)間的延長(zhǎng)而增加，并在第6天達(dá)到6.82（lg（CFU/g）），超過(guò)無(wú)公害水產(chǎn)品安全限值所規(guī)定的6.0（lg（CFU/g）），表示已到達(dá)小龍蝦的貨架期。相比TVB-N值，菌落總數(shù)預(yù)示的貨架期終點(diǎn)稍早，可能是由于細(xì)菌的增長(zhǎng)主導(dǎo)TVB-N的產(chǎn)生，故其增長(zhǎng)滯后于菌落總數(shù)的增長(zhǎng)。這一結(jié)果與符莎露[19]所研究的哈氏仿對(duì)蝦結(jié)果一致。結(jié)合TVB-N值與菌落總數(shù)結(jié)果，可認(rèn)為貯藏0、2、8 d分別為小龍蝦鮮樣期、腐敗中期和腐敗末期，在此時(shí)間節(jié)點(diǎn)提取小龍蝦樣品總DNA，用于高通量測(cè)序。并同時(shí)在第8天進(jìn)行小龍蝦疑似單菌落優(yōu)勢(shì)腐敗菌的提取。

2.2 培養(yǎng)基鑒定法

2.2.1 優(yōu)勢(shì)腐敗菌的分離

蝦肉樣品通過(guò)稀釋涂布平板和平板劃線法，根據(jù)菌落的形態(tài)、顏色差異，在TSA、BP、CFC、MRS、VRBGA、孟加拉紅6 種培養(yǎng)基上，分別分離出4、2、2、1、2、0共11 株菌株。編號(hào)依次為T1、T2、T3、T4、B1、B2、C1、C2、M、V1、V2。

2.2.2 優(yōu)勢(shì)腐敗菌的腐敗特性

將上述11 株菌制成菌懸液，并將菌濃度調(diào)至106CFU/mL后，分別接種無(wú)菌小龍蝦肉，研究冷藏時(shí)間對(duì)接種蝦肉TVB-N值的影響，結(jié)果如圖3所示。實(shí)驗(yàn)組的TVB-N值均高于對(duì)照組，且在第6天時(shí)，絕大部分實(shí)驗(yàn)組的TVB-N值已超出30 mg/100 g，而對(duì)照組在第9天時(shí)，其TVB-N值仍未到達(dá)這一限值。貯藏至第12天時(shí)，對(duì)照組的TVB-N值為39.39 mg/100 g，而實(shí)驗(yàn)組T2組最大，為82.41 mg/100 g，V2組最低，為50.45 mg/100 g。故確定這11 株菌株均為小龍蝦的優(yōu)勢(shì)腐敗菌。

圖3 優(yōu)勢(shì)腐敗菌對(duì)小龍蝦TVB-N值的影響Fig. 3 Change in TVB-N value during storage of crayfish inoculated with suspected dominant spoilage bacteria

圖4 優(yōu)勢(shì)腐敗菌對(duì)小龍蝦菌落總數(shù)的影響Fig. 4 Change in total bacteria count during storage of crayfish inoculated with suspected dominant spoilage bacteria

為驗(yàn)證這些菌株的腐敗能力，將11 株菌株接種到無(wú)菌小龍蝦肉上，研究其對(duì)小龍蝦肉菌落總數(shù)的影響。如圖4所示，當(dāng)貯藏第6天時(shí)，T2、T3、T4、B1、B2、C1、V1組均超過(guò)無(wú)公害水產(chǎn)品安全限值所規(guī)定的6.0（lg（CFU/g）），而對(duì)照組則在第12天時(shí)才超出這一限值。貯藏至第12天時(shí)，菌落總數(shù)從高至低依次為C2、T2、V1、T3、C1、T1、B1、B2、T4、M、V2，分別達(dá)到7.87、7.54、7.52、7.50、7.45、7.36、7.31、7.25、7.12、7.10、7.10（lg（CFU/g）），高于對(duì)照組6.56（lg（CFU/g））。結(jié)合TVB-N值，可知致腐能力最強(qiáng)的腐敗菌為C2、T2、V1、T3、C1，較弱的為T4、M、V2。

2.2.3 優(yōu)勢(shì)腐敗菌的分子鑒定

利用0.8%的瓊脂糖凝膠電泳檢測(cè)冷藏小龍蝦優(yōu)勢(shì)腐敗菌的16S rDNA擴(kuò)增產(chǎn)物，均在1 500 bp左右得到特異性條帶，電泳圖譜如圖5所示。

圖5 冷藏小龍蝦優(yōu)勢(shì)腐敗菌的16S rDNA電泳圖譜Fig. 5 Electropherogram of PCR amplified 16S rDNA genes from four dominant spoilage bacteria in red swamp crayfish

圖6 基于16S rDNA序列的同源性菌株系統(tǒng)發(fā)育樹Fig. 6 Phylogenetic tree of strains based on the 16S rDNA sequences

將測(cè)序所得序列使用NCBI的BLAST分析，選取同源性在98%以上的菌株序列，利用軟件MEGA 7.0.26進(jìn)行對(duì)比，構(gòu)建Neighbor-Joining系統(tǒng)發(fā)育樹（圖6）。菌株T1～T4分別屬于金黃桿菌（Chryseobacterium anthropi）、微球菌屬（Micrococcus endophyticus）、微小桿菌屬（Exiguobacterium himgiriensis）、噬冷菌屬，B1、C1屬于希瓦氏菌屬，B2為格氏乳球菌屬（Lactococcus garvieae），C2為氣單胞菌，M～V2則分別為巨型球菌屬（Macrococcus caseolyticus）、糞腸球菌屬（Enterococcus faecalis）、弗氏檸檬酸桿菌。

結(jié)合致腐能力分析，可知?dú)鈫伟?、微球菌屬、糞腸球菌屬、微小桿菌屬、希瓦氏菌屬致腐能力最強(qiáng)，巨型球菌屬和弗氏檸檬酸桿菌致腐能力較弱。Macé等[20]研究熱帶熟蝦，發(fā)現(xiàn)致腐能力最強(qiáng)的是希瓦氏菌。唐文靜等[16]指出冷藏海鱸魚中致腐能力最強(qiáng)的是草莓假單胞菌。致腐能力的差異可能由物種、菌種差異所致。

2.3 高通量測(cè)序法鑒定優(yōu)勢(shì)腐敗菌

2.3.1 PCR擴(kuò)增結(jié)果及微生物數(shù)量分析

小龍蝦不同腐敗時(shí)期的總DNA提取物，經(jīng)PCR擴(kuò)增V3～V4可變區(qū)后瓊脂糖凝膠電泳如圖7所示。均在1 500 bp處有特異性條帶，純度達(dá)到送檢要求。

圖7 高通量測(cè)序法下V3～V4可變區(qū)DNA電泳圖Fig. 7 Electropherogram of PCR-amplified V3-V4 variable region

采用Uclust分類法，將拼接好的序列聚類為操作分類單元（operational taxonomic units，OTU），根據(jù)貯藏時(shí)間0、2、8 d所對(duì)應(yīng)的OTU繪制韋恩圖。如圖8所示，3 個(gè)不同時(shí)期共有的微生物共有102 種。鮮樣期與腐敗中期共有微生物121 種，低于腐敗中期與末期共有的413 種。說(shuō)明腐敗中期與腐敗末期的物種相似度很高，產(chǎn)生差異的OTU豐度較低。

圖8 不同貯藏時(shí)期小龍蝦OTU韋恩圖Fig. 8 Venn diagram of OTUs in red swamp crayfish at different storage times

2.3.2 物種多樣性分析

Alpha多樣性多用于評(píng)價(jià)微生物菌群的豐富度及其個(gè)體分配均勻性[21]。常用的度量指數(shù)包括Shannon、Simpson、Chao1等。其中，Shannon和Simpson為群落多樣性指數(shù)，數(shù)值越大，表明菌落多樣性越高。而Chao1指數(shù)則用于表示樣品菌群豐富度，Chao1指數(shù)越大，菌落豐富度則越高[22]。

表1 小龍蝦不同貯藏時(shí)期菌群的alpha多樣性比較Table 1 Alpha diversity of bacterial fl ora in crayfish at different storage times

由表1可知，3 個(gè)時(shí)期有效序列范圍為17 034～29 324 條。在相似水平97%上對(duì)其進(jìn)行聚類分析（OTU），得到鮮樣的OTU數(shù)為452，腐敗中期OTU數(shù)為434，腐敗末期OTU數(shù)為327。表明隨著貯藏時(shí)間延長(zhǎng)，小龍蝦的微生物種類呈現(xiàn)降低趨勢(shì)。此外，與鮮樣組比較，蝦肉中期及末期組的Shannon、Simpson、Chao1指數(shù)均減小，說(shuō)明蝦肉樣品中微生物物種多樣性和豐富度均下降。其中，Simpson指數(shù)變化較小，這是由于該數(shù)值主要綜合體現(xiàn)物種的豐富度和均勻度，但不只是考慮種類數(shù)目的變化，還考慮了種類之間的配比變化。數(shù)值變化不大，則說(shuō)明盡管物種種類數(shù)目有一定變化，但整體各物種間的相對(duì)比例變化不明顯。由此可以驗(yàn)證特定腐敗菌的生長(zhǎng)繁殖抑制其他微生物菌群生長(zhǎng)結(jié)論的正確性。

稀釋曲線是以測(cè)序數(shù)量及其對(duì)應(yīng)的OTU數(shù)構(gòu)建的、可用來(lái)評(píng)估樣品中微生物群落物種豐度的曲線[23]。其可直接反映測(cè)序數(shù)量的合理性，并可通過(guò)物種變化趨勢(shì)間接反映樣品物種的豐度[24]。如圖9A所示，當(dāng)樣品測(cè)序數(shù)超過(guò)16 000時(shí)，其稀釋曲線趨于平緩，表明測(cè)序數(shù)量比較合理，再次增加測(cè)序量，只會(huì)產(chǎn)生極少新OTU，故其基本已覆蓋到樣品中所有的物種[25]。Shannon曲線（圖9B）也在此處趨于平緩，說(shuō)明測(cè)序數(shù)量足夠充分合理，已包含所測(cè)樣品中絕大多數(shù)菌群信息。

如圖10所示，描述物種相對(duì)豐度情況。在相對(duì)豐度曲線圖的水平方向（X軸），曲線寬度反映物種豐度，即曲線在X軸上的范圍越大代表物種豐度越高；曲線形狀的平滑度反映物種組成的均勻程度，一般曲線越平滑，物種分布越均勻[26]。從圖10可知，3 組樣品在X軸上的分布長(zhǎng)度由0、2、8 d依次遞減，表明0 d的物種豐度最大，8 d豐度最小。此外，8 d豐度曲線的陡峭程度大于其他兩組，說(shuō)明該組樣品菌群組成分布最不均勻，其內(nèi)部可能存在比列較大的菌群。

圖9 不同貯藏時(shí)間小龍蝦樣品稀釋曲線（A）和Shannon曲線（B）分析Fig. 9 Rarefaction curves (A) and Shannon curves (B) of crayfish at different storage times

圖10 不同貯藏時(shí)間小龍蝦樣品微生物等級(jí)豐度曲線Fig. 10 Rank abundance curves of bacterial species in crayfish at different storage times

2.3.3 菌相組成分析

小龍蝦在冷藏期間的全部菌落按屬分類主要包括動(dòng)球菌屬希瓦氏菌、黃桿菌屬（Flavobacterium sp.）、金黃桿菌屬、不動(dòng)桿菌屬、嗜冷菌屬、巨型球菌屬、布丘氏菌屬（Buttiauxella sp.）、微球菌屬、肉食桿菌、異常球菌屬（Deinococcus sp.）、環(huán)絲菌屬、氣單胞菌屬、叢毛單胞菌屬（Comamonas sp.）、節(jié)桿菌屬（Arthrobacter sp.）、鞘脂桿菌屬（Sphingobacterium sp.）、漫游球菌屬、羅氏菌屬（Rothia sp.）、假單胞菌屬（Pseudomonas sp.）等。

如圖11所示，腐敗末期樣品占據(jù)主要地位的菌屬依次為希瓦氏菌、肉食桿菌、環(huán)絲菌屬、嗜冷菌屬、漫游球菌屬、不動(dòng)桿菌屬，六者之和物種豐度占據(jù)90.57%。與鮮樣菌屬比較，除不動(dòng)桿菌屬由初期的4.76%下降至2.19%，其他5 種比列均上升。希瓦氏菌由10.04%升至63.19%，肉食桿菌、環(huán)絲菌屬均由0.02%分別升至7.33%、6.48%，嗜冷菌屬和漫游球菌屬則分別由1.91%、0.20%升至6.32%和2.93%。吳海虹等[27]比較不同包裝方式下冷藏青蝦菌群差異，發(fā)現(xiàn)空氣冷藏條件下，優(yōu)勢(shì)菌群主要是不動(dòng)桿菌、黃桿菌和希瓦氏菌；黃佳奇[26]發(fā)現(xiàn)空氣條件下，冷藏小黃魚的優(yōu)勢(shì)腐敗菌主要是希瓦氏菌屬、變形桿菌屬、普羅威斯登菌屬和動(dòng)性球菌屬，可見優(yōu)勢(shì)腐敗菌存在品種差異。

與傳統(tǒng)培養(yǎng)基法相比，高通量測(cè)序法檢測(cè)到的菌種類更多。除分離的11 株單菌外，高通量法還檢測(cè)到黃桿菌屬、布丘氏菌屬、微球菌屬、肉食桿菌、異常球菌屬、環(huán)絲菌屬、叢毛單胞菌屬等。故該方法可提高菌種數(shù)量測(cè)定的準(zhǔn)確性，彌補(bǔ)傳統(tǒng)法肉眼識(shí)別的遺漏欠缺，有效提升菌種鑒定水平，且可直接得到對(duì)應(yīng)菌種所占的比列，極大滿足研究需要。但較培養(yǎng)基法而言，其測(cè)序所耗時(shí)間也有所延長(zhǎng)。

圖11 不同貯藏時(shí)間小龍蝦樣品微生物物種組成（屬水平）Fig. 11 Bacterial community composition at genus level in crayfish at different storage times

3 結(jié)論與討論

本實(shí)驗(yàn)利用傳統(tǒng)培養(yǎng)基法和高通量測(cè)序法分析測(cè)定冷藏小龍蝦菌屬變化情況，傳統(tǒng)法提取的11 株菌株基本全包含在高通量法測(cè)序結(jié)果菌株中，可以驗(yàn)證高通量法的高效、精準(zhǔn)性。傳統(tǒng)法共分離11 株菌株，并對(duì)這些菌株進(jìn)行致腐能力測(cè)定，結(jié)果表明氣單胞菌致腐能力最強(qiáng)，弗氏檸檬酸桿菌致腐能力最弱。高通量法則分別比較小龍蝦鮮樣期、腐敗中期及腐敗末期菌屬變化，發(fā)現(xiàn)隨著貯藏時(shí)間延長(zhǎng)，菌屬多樣性下降但某些菌屬的豐度呈現(xiàn)上升趨勢(shì)，尤其在腐敗末期時(shí)，希瓦氏菌、肉食桿菌、環(huán)絲菌屬、嗜冷菌屬、漫游球菌屬、不動(dòng)桿菌屬6 種菌屬共占據(jù)物種豐度的90.57%，證實(shí)特定腐敗菌的生長(zhǎng)繁殖會(huì)抑制其他微生物菌群生長(zhǎng)，并可以確定小龍蝦的優(yōu)勢(shì)腐敗菌為希瓦氏菌、肉食桿菌、環(huán)絲菌屬、嗜冷菌屬、漫游球菌屬、不動(dòng)桿菌屬、氣單胞菌、束毛球菌屬等。

已有研究表明希瓦氏菌屬和假單胞菌屬為海產(chǎn)品的主要優(yōu)勢(shì)腐敗菌，而假單胞菌屬則為淡水產(chǎn)品的優(yōu)勢(shì)腐敗菌[28]。目前已證實(shí)大黃魚在低溫貯藏條件下的優(yōu)勢(shì)腐敗菌為希瓦氏菌[29]，而廣州鯽魚、鳊魚和羅非魚3 種淡水魚的優(yōu)勢(shì)腐敗菌為假單胞菌屬[30]。但本實(shí)驗(yàn)得到淡水小龍蝦的主要優(yōu)勢(shì)腐敗菌希瓦氏菌等，這與已有結(jié)論略有差異?？赡苁怯捎谒a(chǎn)品的種類、生活水域、捕捉方法以及貯藏方式的不同使其特定腐敗菌也不同[31]。此外，希瓦氏菌具有強(qiáng)大的代謝系統(tǒng)，能適應(yīng)多種鹽濃度、溫度及氣體環(huán)境等，是一種嗜冷性及產(chǎn)硫能力極強(qiáng)的菌屬。故其在冷藏至腐敗末期，含量劇增，成為冷藏小龍蝦的優(yōu)勢(shì)腐敗菌屬之一。本實(shí)驗(yàn)結(jié)果可為后期研究抑制優(yōu)勢(shì)腐敗菌生長(zhǎng)、延長(zhǎng)小龍蝦貨架期提供參考基礎(chǔ)。