宋思家,黃子彧,林瑩瑩,郭慧媛*
1(中國(guó)農(nóng)業(yè)大學(xué) 食品科學(xué)與營(yíng)養(yǎng)工程學(xué)院,精準(zhǔn)營(yíng)養(yǎng)與食品質(zhì)量重點(diǎn)實(shí)驗(yàn)室,教育部功能乳品重點(diǎn)實(shí)驗(yàn)室,北京,100083) 2(中國(guó)農(nóng)業(yè)大學(xué) 食品質(zhì)量安全北京實(shí)驗(yàn)室,北京,100083)
雞蛋干是中國(guó)市面上常見的休閑蛋制品,因其營(yíng)養(yǎng)豐富、風(fēng)味獨(dú)特而深受消費(fèi)者喜愛。依據(jù)GB 2749—2015《食品安全國(guó)家標(biāo)準(zhǔn) 蛋與蛋制品》內(nèi)容,在保質(zhì)期內(nèi)的雞蛋干在感官方面應(yīng)具有產(chǎn)品正常的色澤、滋味、氣味,無異味,且具有產(chǎn)品正常的形狀、形態(tài)、無酸敗、霉變等;在微生物限量方面,其菌落總數(shù)應(yīng)<104CFU/g,大腸菌群應(yīng)<10 CFU/g[1]。目前在實(shí)際生產(chǎn)過程中,雞蛋干的保質(zhì)期為9個(gè)月。雞蛋干的生產(chǎn)原料以全蛋液或蛋清液為主,在加工制作過程中極易受到微生物污染而導(dǎo)致產(chǎn)品腐敗變質(zhì)[2]。相關(guān)調(diào)查研究表明,預(yù)包裝雞蛋干即便經(jīng)過高壓蒸汽滅菌后也可能出現(xiàn)脹包、湯汁渾濁、有異味的變質(zhì)現(xiàn)象,尤其在夏秋高溫季節(jié),雞蛋干發(fā)生腐敗變質(zhì)的概率更高[3]。最近數(shù)十年,雞蛋干的市場(chǎng)在不斷擴(kuò)大,因此對(duì)雞蛋干的防腐保鮮開展研究愈顯重要。而明確雞蛋干的優(yōu)勢(shì)致腐微生物,是防腐的首要關(guān)鍵。然而目前對(duì)于雞蛋干的腐敗微生物群落多樣性的報(bào)道十分有限,賀燕等[3]從腐敗雞蛋干中分離鑒定得到2株芽孢桿菌,袁先鈴等[4]在腐敗雞蛋干中檢測(cè)到彭氏變形菌和奇異變形菌。因此,開展雞蛋干中腐敗微生物群落多樣性的研究具有重要意義。
對(duì)于微生物群落多樣性的研究,長(zhǎng)期以來主要是基于《伯杰細(xì)菌鑒定手冊(cè)》,對(duì)微生物菌落進(jìn)行分離、純化,進(jìn)而通過形態(tài)學(xué)觀察和一系列復(fù)雜的生理生化實(shí)驗(yàn)進(jìn)行鑒定。這種傳統(tǒng)方法不但耗費(fèi)時(shí)間和精力,而且無法對(duì)微生物進(jìn)行精確鑒定,不能檢測(cè)到微生物群落多樣性的真正概貌[5-6]。因?yàn)橛行┪⑸锏呐囵B(yǎng)條件非??量?,甚至不可培養(yǎng),常規(guī)方法難以檢測(cè)到不可培養(yǎng)的微生物。近年來,在微生物群落多樣性研究中已經(jīng)越來越傾向于采用Illumina MiSeq高通量測(cè)序技術(shù)[7-10]。高通量測(cè)序技術(shù)也稱“下一代測(cè)序技術(shù)”,可同時(shí)測(cè)序數(shù)百萬個(gè)核酸分子,根據(jù)16S rDNA基因序列分析特定樣品中微生物群落的構(gòu)成情況[11-12]。目前高通量測(cè)序技術(shù)已成為監(jiān)測(cè)食品貯藏過程中微生物群落演替、分析發(fā)酵及腐敗食品微生物群落結(jié)構(gòu)等食品微生物多樣性研究的首選方法[11]。
本研究以保質(zhì)期不同的腐敗變質(zhì)雞蛋干為研究對(duì)象,基于Illumina Miseq測(cè)序平臺(tái)對(duì)其16S rDNA V3~V4進(jìn)行擴(kuò)增子測(cè)序分析,研究其微生物群落結(jié)構(gòu)組成和相對(duì)豐度,探究造成雞蛋干腐敗變質(zhì)的核心優(yōu)勢(shì)微生物以及造成雞蛋干保質(zhì)期差異的微生物種類,為靶向抑制雞蛋干的腐敗變質(zhì)提供理論依據(jù)。
雞蛋,北京德青源農(nóng)業(yè)科技股份有限公司;細(xì)菌基因組DNA提取試劑盒,美國(guó)Omega科技公司;Axy Prep DNA凝膠回收試劑盒,美國(guó)Axygen公司;2×Taq PCR MasterMix、DNA Marker,天根生化科技(北京)有限公司。
5418R高速離心機(jī),德國(guó)Eppendorf公司;拍擊式均質(zhì)機(jī),法國(guó)Interscience公司;RS-400型真空包裝機(jī),北京日上公司;T100 Thermal Cycler 聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction,PCR)儀、PowerPac Basic型電泳儀,美國(guó)Bio-Rad公司;ChampGel5000增強(qiáng)型全自動(dòng)凝膠成像儀,北京賽智創(chuàng)業(yè)科技有限公司;Illumina MiSeq測(cè)序儀,美國(guó)Illumina公司。
1.3.1 樣品制備
雞蛋干的制作工藝參照參考文獻(xiàn)[13],具體為取預(yù)處理蛋液于模具中,水浴加熱使其中心溫度達(dá)到80 ℃,并維持20 min,待其冷卻收縮后脫模。將雞蛋干樣品真空包裝后,進(jìn)行90 ℃、30 min滅菌,于25 ℃保存直至變質(zhì)。預(yù)實(shí)驗(yàn)發(fā)現(xiàn),對(duì)于使用不同批次原料、于不同的時(shí)間點(diǎn)制作的雞蛋干樣品,即使其加工工藝完全一致,樣品的保質(zhì)期也會(huì)有較大差異。為探究不同保質(zhì)期的雞蛋干樣品中菌群結(jié)構(gòu)的差異,采用上述加工工藝,分別使用不同批次原料,于不同的時(shí)間點(diǎn)制作了多批樣品,從其中選取貯藏15 d、22 d和29 d腐敗變質(zhì)的3批雞蛋干產(chǎn)品,取樣保存于-80 ℃冰箱待測(cè)。將3批雞蛋干樣品分別編號(hào)為15 d、22 d和29 d,每組樣品設(shè)置3個(gè)平行。
1.3.2 細(xì)菌總DNA的提取
參考高乾坤等[14]的方法,對(duì)樣品進(jìn)行處理后,采用美國(guó)Omega公司的細(xì)菌基因組DNA提取試劑盒,提取細(xì)菌總DNA,提取步驟參照試劑盒產(chǎn)品說明書,并進(jìn)行1%(質(zhì)量分?jǐn)?shù))的瓊脂糖凝膠電泳,以檢測(cè)提取效果。
1.3.3 PCR擴(kuò)增及高通量測(cè)序
以檢驗(yàn)合格的基因組DNA為模板,擴(kuò)增細(xì)菌16S rDNA基因V3~V4區(qū)。正向引物為338F(5′-ACTCCTACGGGAGGCAGCA-3′),反向引物為806R(5′-GGACTACHVGGGTWTCTAAT-3′)。所得PCR產(chǎn)物經(jīng)2%(質(zhì)量分?jǐn)?shù))瓊脂糖凝膠電泳檢測(cè)擴(kuò)增效果,然后將樣品進(jìn)行純化并定量檢測(cè)。將檢測(cè)合格的樣品送至上海美吉生物醫(yī)藥科技有限公司進(jìn)行Illumina MiSeq高通量測(cè)序。
對(duì)MiSeq測(cè)序得到的雙端序列數(shù)據(jù)進(jìn)行處理,首先將成對(duì)的reads拼接成一條序列,同時(shí)對(duì)reads的質(zhì)量和拼接的效果進(jìn)行質(zhì)控過濾、得到有效序列,并校正序列方向[15]。
基于操作分類單元(operational taxonomic units, OTU)聚類分析結(jié)果,研究腐敗變質(zhì)雞蛋干樣品中微生物的多樣性,通過Alpha多樣性檢測(cè)單樣本測(cè)序深度,并分析微生物群落的豐度和多樣性[16]?;诜诸悓W(xué)信息,在各個(gè)分類水平上進(jìn)行群落結(jié)構(gòu)統(tǒng)計(jì),綜合討論保質(zhì)期不同的雞蛋干樣品中微生物群落結(jié)構(gòu)組成和相對(duì)豐度。
采用Illumina Miseq高通量測(cè)序獲得了3種不同保質(zhì)期雞蛋干樣品的原始序列數(shù)據(jù)后,經(jīng)Fastp及FLASH軟件優(yōu)化數(shù)據(jù)后得到了9個(gè)樣本共529 384條有效序列,序列平均長(zhǎng)度為409.70~429.92 bp,主要集中在401~440 bp(圖1-a)。稀釋曲線可反映在高通量測(cè)序過程中的樣品取樣深度,因此是評(píng)價(jià)測(cè)序數(shù)據(jù)量是否全面覆蓋微生物群落的常用指標(biāo)[14]。為了探究雞蛋干樣品的測(cè)序數(shù)據(jù)量是否科學(xué),從測(cè)序獲得的所有原始數(shù)據(jù)中隨機(jī)抽取數(shù)據(jù),并統(tǒng)計(jì)其物種多樣性指數(shù),得到稀釋曲線(圖1-b和圖1-c)。結(jié)果表明,隨著數(shù)據(jù)量的增加,OTU曲線和Shannon-Wiener曲線的趨勢(shì)已經(jīng)逐漸平緩,即再增加數(shù)據(jù)量對(duì)于發(fā)現(xiàn)新物種也無顯著影響。樣本的OTU覆蓋度已經(jīng)基本飽和,更多的數(shù)據(jù)量對(duì)發(fā)現(xiàn)新OTU的邊際貢獻(xiàn)很小,這個(gè)結(jié)果說明測(cè)序數(shù)據(jù)量足夠大,可以反映樣本中絕大多數(shù)的微生物多樣性信息[17]。
a-樣品序列長(zhǎng)度分布圖;b-OTU曲線;c-Shannon-Wiener曲線圖1 樣品序列長(zhǎng)度分布圖、OTU曲線和Shannon-Wiener曲線Fig.1 Length distribution profile of trimed sequences, OTU curves and Shannon-Wiener curves
根據(jù)97%相似性水平下的OTU信息,采用Alpha多樣性指標(biāo)的Shannon、Simpson、ACE、Chao1及Coverage指數(shù)對(duì)雞蛋干樣品微生物物種的豐富度和多樣性進(jìn)行評(píng)估[14]。如表1所示,3組樣品的Coverage指數(shù)都達(dá)到了0.999以上,即樣品文庫的覆蓋率足夠高,有未被檢出序列的可能性較小,表明本次試驗(yàn)所建文庫能夠較準(zhǔn)確地反映出待測(cè)雞蛋干的微生物多樣性情況。Chao1、ACE指數(shù)用于計(jì)算菌群豐度,二者的值越大,說明樣品的群落豐富度越高;Shannon、Simpson指數(shù)用于計(jì)算菌群多樣性指數(shù),Shannon的值越大,說明樣品群落多樣性越高,而Simpson指數(shù)值越大,則說明其群落多樣性越低[18]。22 d組的雞蛋干樣品ACE、Chao1及Shannon指數(shù)最高,Simpson值最低,這說明該組樣品的細(xì)菌群落豐富度及多樣性高于其他2組樣品,可能由于該組樣品的原料中及加工過程中污染的細(xì)菌種類最多;29 d組雞蛋干樣品的ACE、Chao1及Shannon指數(shù)均為3組中最小,這說明了該組樣品中細(xì)菌群落豐度及多樣性最低,可能由于該組樣品的原料中及加工過程中受到的污染程度較低。綜上,經(jīng)Alpha多樣性指數(shù)統(tǒng)計(jì)分析可知,雞蛋干樣品的物種多樣性由高到低依次為22 d組、15 d組和29 d組。此結(jié)果顯示,相比于污染的微生物物種的多樣性,雞蛋干所污染微生物的具體種類可能與雞蛋干的保質(zhì)期相關(guān)性更大,這與孫同輝等[11]的報(bào)道相一致。
基于樣品總OTU數(shù)作Veen圖,利用Venn圖可以展示多組樣品間共有和每組樣品特有的OTU數(shù)目,并直觀展示樣品間OTU的重疊情況[19]。如圖2所示,15 d組樣品與22 d組樣品共有的OTU數(shù)107個(gè),15 d組與29 d組共有的OTU數(shù)131個(gè),22 d組與29 d組共有的OTU數(shù)92個(gè),3組樣本共有的OTU數(shù)為81個(gè),占各樣本OTU數(shù)比例較小,而各樣本特有的OTU數(shù)在113~259,表明各樣本之間的菌群結(jié)構(gòu)差異性較大。
表1 Alpha多樣性指數(shù)Table 1 Alpha diversity indexes
圖2 不同保質(zhì)期樣品中微生物物種間OUT分析圖Fig.2 OUT analysis diagram of microbial species in samples with different shelf life
2.4.1 基于門水平的細(xì)菌菌群結(jié)構(gòu)分析
在門水平上,3組雞蛋干樣品的群落組成如圖3所示。15 d組的雞蛋干樣品總細(xì)菌組成以厚壁菌門(Firmicutes)和變形菌門(Proteobacteria)為主,相對(duì)豐度分別為96.73%和1.96%。除此之外還包括少量的放線菌門(Actinobacteria)和擬桿菌門(Bacteroidetes)微生物。而22 d組的雞蛋干樣品菌群結(jié)構(gòu)以厚壁菌門、變形菌門和擬桿菌門為主,相對(duì)豐度分別為87.55%、9.41%和1.32%,同時(shí)還含有少量的放線菌門微生物。29 d組的雞蛋干樣品中微生物以厚壁菌門(70.94%)、放線菌門(27.81%)和變形菌門(1.09%)為主,還含有少量的擬桿菌門微生物。以上結(jié)果表明,在3組不同保質(zhì)期的樣品中,厚壁菌門的豐度均為最高,可見厚壁菌門是腐敗變質(zhì)雞蛋干中的優(yōu)勢(shì)微生物。近年來許多研究表明,厚壁菌門微生物是高蛋白食品中的優(yōu)勢(shì)腐敗菌,在乳制品、肉制品及水產(chǎn)品中都十分常見[18-21],這與本研究的結(jié)果一致。
圖3 基于門分類水平上的樣品菌群分布圖Fig.3 Bacterial distribution pattern at the phylum level
2.4.2 基于科水平的細(xì)菌菌群結(jié)構(gòu)分析
基于科分類水平的3組雞蛋干樣品的微生物群落結(jié)構(gòu)如圖4所示。測(cè)序結(jié)果表明,15 d組的雞蛋干樣品中類芽孢桿菌科(Paenibacillaceae)所占比例最大,為47.43%,其次為芽孢桿菌科(Bacillaceae)和梭菌科(Clostridiaceae_1),所占比例分別為33.14%和14.96%,還有乳桿菌科(Lactobacillaceae)、伯克氏菌科(Burkholderiaceae)、鏈球菌科(Streptococcaceae)、假單胞菌科(Pseudomonadaceae)等相對(duì)豐度較低的菌種,說明雞蛋干在制作和包裝過程中污染微生物的種類復(fù)雜,導(dǎo)致樣品菌落結(jié)構(gòu)呈現(xiàn)多樣性。22 d組的雞蛋干樣品菌群結(jié)構(gòu)以芽孢桿菌科、乳桿菌科、鏈球菌科、伯克氏菌科為主,相對(duì)豐度分別為47.39%、29.11%和3.76%、3.29%,除此之外還包括少量的紅桿菌科(Rhodobacteraceae)、假單胞菌科、類芽孢桿菌科微生物。29 d組雞蛋干樣品中的菌群在科水平上主要有芽孢桿菌科(68.85%)和小單孢菌科(Micromonosporaceae,27.60%),還含有少量的假單胞菌科和伯克氏菌科微生物。以上結(jié)果表明,不同保質(zhì)期雞蛋干樣品所含微生物種類和數(shù)量有明顯差異,如相比于其他2組,29 d組樣品中卻未檢出乳桿菌科微生物,可能由于該組樣品的原料中及加工過程中污染的乳桿菌科較少。29 d組雞蛋干樣品的保質(zhì)期最長(zhǎng),其微生物多樣性低于其余2組,說明樣品中污染的微生物種類和數(shù)量與雞蛋干的保質(zhì)期密切相關(guān)。22 d組樣品微生物多樣性最高,但此組保質(zhì)期不是最短,可能由于樣品中含有豐度較高的乳桿菌科微生物[22-23],對(duì)其他微生物產(chǎn)生拮抗作用、抑制其生長(zhǎng)繁殖,相對(duì)延長(zhǎng)了雞蛋干保質(zhì)期。此外,在15 d組樣品中檢測(cè)出梭菌科微生物,而其他2組中均未檢出,梭菌科中有些菌種為厭氧菌,可在真空包裝雞蛋干中迅速生長(zhǎng)繁殖,這可能是本組雞蛋干保質(zhì)期遠(yuǎn)小于其他2組的原因。3組不同保質(zhì)期的樣品中,均含有芽孢桿菌科微生物,且豐度較高,表明芽孢桿菌科是腐敗變質(zhì)雞蛋干中的核心優(yōu)勢(shì)微生物。
圖4 基于科分類水平上的樣品菌群分布圖Fig.4 Bacterial distribution pattern at the family level
2.4.3 基于屬水平的細(xì)菌菌群結(jié)構(gòu)分析
選取在雞蛋干樣品豐度排名前15的屬,根據(jù)其在每組雞蛋干樣品中的豐度信息聚類繪制成熱圖(圖5-a)。熱圖可通過顏色變化直觀地反映出各組雞蛋干樣品群落組成的相似性和差異性,更便于發(fā)現(xiàn)在各組樣品中分布較多或較少的菌屬[10]。圖中顏色越偏磚紅色,表明樣品中所含該菌屬比例越高,顏色越偏深藍(lán)色,比例則越低。如圖5-a所示,在屬水平,22 d組和29 d組聚為一類,15 d組單獨(dú)聚為一類,說明22 d組和29 d組菌群整體結(jié)構(gòu)相似度較高,與15 d組相似度較低,表明雞蛋干樣品的保質(zhì)期與其菌群結(jié)構(gòu)相關(guān)性較大。
從圖5-b各樣品菌群柱狀分布圖中可以看出,在屬水平上,15 d組優(yōu)勢(shì)菌群主要是芽孢桿菌屬 (Bacillus,33.14%)、短芽孢桿菌屬(Brevibacillus,30.96%)、類芽孢桿菌屬(Paenibacillus,16.46%)和狹義梭菌屬(Clostridium_sensu_stricto_13,14.86%)。在優(yōu)勢(shì)菌群中,芽孢桿菌屬所占比例最大,它屬于需氧或兼性厭氧菌,此屬中大部分細(xì)菌為腐生菌,廣泛分布于自然環(huán)境中。該菌屬可形成芽孢,對(duì)不良環(huán)境條件有很強(qiáng)的抵抗力,不易殺滅。賀燕等[3]在腐敗雞蛋干中也有檢測(cè)到芽孢桿菌屬微生物。短芽孢桿菌屬為需氧或兼性厭氧菌,近年來相關(guān)報(bào)道多集中于該屬工程菌株的研究,鮮見有研究發(fā)現(xiàn)其在食品腐敗中的作用[24-25]。類芽孢桿菌屬為需氧或兼性厭氧菌,近年來許多研究表明,其為多種乳制品和即食食品的腐敗微生物[26-27]。狹義梭菌屬?gòu)V泛分布于土壤、下水污泥、人和其他哺乳動(dòng)物的腸道等環(huán)境中,是乳制品和肉制品中常見的腐敗微生物[20, 28-29]。
22 d組優(yōu)勢(shì)菌群主要是芽孢桿菌屬(47.38%)、乳桿菌屬(Lactobacillus,28.26%)和鏈球菌屬(Streptococcus,3.57%)。乳桿菌屬為兼性厭氧菌,廣泛分布于自然界中,污染乳制品、肉制品、水產(chǎn)品等高蛋白食品后可引起腐敗變質(zhì)。鏈球菌屬為兼性厭氧菌,該屬中有少數(shù)種為腐生菌,污染食品后可引起腐敗變質(zhì)[29]。
a-樣品菌群相對(duì)豐度熱圖;b-樣品菌群相對(duì)豐度柱形圖圖5 基于屬分類水平上的樣品菌群相對(duì)豐度熱圖 和柱形圖Fig.5 Heat map and stacked-column of relative abundances of bacterial communities at the genus level
29 d組優(yōu)勢(shì)菌屬主要是芽孢桿菌屬(68.85%)和小單孢菌屬(Micromonospora,27.60%)。小單孢菌屬為好氧或微好氧的腐生菌,是土壤和水域的常見微生物群體[30],其在真空包裝雞蛋干中占比達(dá)到27.60%,可能由于雞蛋干中初始污染小單孢菌屬較多,且微好氧微生物能在較低的氧分壓下生長(zhǎng),雞蛋干雖然經(jīng)過真空包裝,但仍無法將樣品中氧氣完全排出,因此小單孢菌屬中部分微好氧菌種可以生長(zhǎng)。
以上結(jié)果表明,在3組不同保質(zhì)期的樣品中,芽孢桿菌屬的相對(duì)豐度均為最高,可見芽孢桿菌屬是雞蛋干的核心致腐微生物。保質(zhì)期最長(zhǎng)的29 d組樣品中優(yōu)勢(shì)菌屬為芽孢桿菌屬和小單孢菌屬,其中小單孢菌屬為好氧或微好氧菌,在真空包裝雞蛋干中生長(zhǎng)緩慢,而芽孢桿菌屬中部分菌種為兼性厭氧菌,可在真空包裝條件下良好生長(zhǎng),因此芽孢桿菌屬應(yīng)為限制此組樣品保質(zhì)期的微生物。而相比于29 d組,22 d組樣品中還含有較高豐度的乳桿菌屬,這可能是造成本組雞蛋干保質(zhì)期<29 d的原因。此外,狹義梭菌屬僅在保質(zhì)期為15 d的雞蛋干樣品中檢測(cè)到,它屬于厭氧菌,因此在真空包裝雞蛋干中的生長(zhǎng)速度大于非厭氧菌,可快速分解雞蛋干中營(yíng)養(yǎng)物質(zhì)造成腐敗變質(zhì),這可能是15 d組雞蛋干保質(zhì)期遠(yuǎn)低于其他2組的原因。
本研究以保質(zhì)期不同的3組腐敗變質(zhì)雞蛋干為研究對(duì)象,采用Illumina Miseq高通量測(cè)序技術(shù)分析樣品中菌群組成情況。結(jié)果表明,3組腐敗變質(zhì)雞蛋干的微生物的多樣性以及豐度存在較大差異,而總體來看,厚壁菌門為腐敗變質(zhì)雞蛋干的優(yōu)勢(shì)菌門,芽孢桿菌科為優(yōu)勢(shì)菌科。在屬水平,芽孢桿菌屬在3組樣品中的相對(duì)豐度均為最高,是雞蛋干的核心致腐微生物。此外,造成雞蛋干保質(zhì)期差異的菌屬包括乳桿菌屬和狹義梭菌屬,它們可加速雞蛋干的腐敗變質(zhì),是雞蛋干保質(zhì)期的限制因素,應(yīng)在殺菌防腐環(huán)節(jié)重點(diǎn)關(guān)注。高通量測(cè)序技術(shù)檢測(cè)不同保質(zhì)期雞蛋干樣品的微生物組成,不僅可以分析雞蛋干腐敗過程中的核心優(yōu)勢(shì)菌群,還可以高效快速地分析樣品中微生物種類的差異,為今后雞蛋干質(zhì)量安全監(jiān)測(cè)及靶向研究雞蛋干防腐保鮮技術(shù)提供理論依據(jù)。