東北酸菜發(fā)酵過程中細(xì)菌菌群與亞硝酸鹽長消變化關(guān)系

張慶芳 袁文濤 孫溪 魏振勇 周姝靜 遲乃玉 于爽

摘要：為探究東北酸菜發(fā)酵過程中細(xì)菌菌群與亞硝酸鹽長消變化關(guān)系，實(shí)驗(yàn)設(shè)立了兩種封閉發(fā)酵系統(tǒng)，即未添加亞硝酸鹽（自然發(fā)酵組）、添加亞硝酸鹽（對照組），結(jié)合主要理化指標(biāo)（亞硝酸鹽、pH值等）對兩個系統(tǒng)內(nèi)不同發(fā)酵階段的酸菜發(fā)酵液（相對豐度排名前30）中的菌屬進(jìn)行比對分析。實(shí)驗(yàn)研究了不同發(fā)酵階段有、無添加亞硝酸鹽樣品的菌群變化，結(jié)果表明，添加亞硝酸鹽對整個發(fā)酵周期中細(xì)菌的生長都具有影響，亞硝酸鹽使優(yōu)勢乳酸菌菌屬（Lactococcus、Leuconostoc、Lactobacillus）、α-變形菌菌屬（Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、Aureimonas、Brevundimonas、Methylobacterium、Sphingomonas）、β-變形菌菌屬（Acidovorax、Aquabacterium、Comamonas、Delftia、Herbaspirillum、Methylophilus、Methylovorus）、CFB類菌屬（Dyadobacter、Pedobacter、Sphingobacterium）和未培養(yǎng)菌屬（Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast）在發(fā)酵初期相對豐度增加；使致病菌菌屬（Acinetobacter、Alkanindiges、Chryseobacterium、Enterococcus、Serratia、uncultured_bacterium_f_Burkholderiaceae、uncultured_bacterium_f_Enterobacteriaceae、uncultured_bacterium_f_Xanthobacteraceae）相對豐度始終處于低水平狀態(tài)。在S系統(tǒng)內(nèi)，發(fā)酵前期為亞硝酸鹽積累階段，此時(shí)Lactococcus、Leuconostoc為優(yōu)勢乳酸菌菌屬，而達(dá)到發(fā)酵中后期時(shí)Lactobacillus成為絕對優(yōu)勢乳酸菌菌屬；在發(fā)酵初期，Acinetobacter、Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae等致病菌大量繁殖，導(dǎo)致亞硝酸鹽大量積累，而積累的亞硝酸鹽使這些致病菌的相對豐度始終處于低水平狀態(tài)，其中對Serratia的生長影響最顯著。

關(guān)鍵詞：東北酸菜；亞硝酸鹽；細(xì)菌菌群；優(yōu)勢乳酸菌；致病菌

中圖分類號：TS255.53

文獻(xiàn)標(biāo)志碼：A

文章編號：1000-9973（2024）06-0050-14

Relationships Between Bacterial Flora and Growth and Decline of Nitrite

During Fermentation of Northeast Sauerkraut

ZHANG Qing-fang^1，2， YUAN Wen-tao^1，2， SUN Xi^1，2， WEI Zhen-yong^1，2，

ZHOU Shu-jing^1，2， CHI Nai-yu^1，2*， YU Shuang^1，2*

（1.College of Life and Health， Dalian University， Dalian 116622， China; 2.Engineering and Technology

Research Center for Marine Microorganism in Liaoning Province， Dalian 116622， China）

Abstract： To investigate the relationship between bacterial flora and growth and decline of nitrite during fermentation of northeast sauerkraut， in this experiment， two closed fermentation systems are set up， namely， nitrite is not added （natural fermentation group） and nitrite is added （control group）. In combination with the main physicochemical indexes （nitrite， pH value and so on）， the bacterial species in the fermentation broth of sauerkraut at different fermentation stages （relative abundance ranking top 30） in the two systems are compared and analyzed. In this experiment， the changes of the bacterial flora of samples with and without nitrite at different fermentation stages are investigated. The results show that the addition of nitrite has an effect on the growth of bacteria throughout the fermentation cycle， and that nitrite makes the relative abundance of the dominant Lactobacillus genera （Lactococcus， Leuconostoc， Lactobacillus）， α-Proteus （Allorhizobium-Neorhizobium-Parararhizobium-Rhizobium， Aureimonas， Brevundimonas， Methylobacterium， Sphingomonas）， β-Proteus （Acidovorax， Aquabacterium， Comamonas， Delftia， Herbaspirillum， Methylophilus， Methylovorus）， CFB genera （Dyadobacter， Pedobacter， Sphingobacterium） and uncultured genera （Brassica_napus_rape， uncultured_bacterium_f_Chitinophagaceae， uncultured_bacterium_f_env.OPS_17， uncultured_bacterium_o_Chloroplast） increase in the early stages of fermentation; and it makes the relative abundance of pathogenic bacteria genera （Acinetobacter， Alkanindiges， Chryseobacterium， Enterococcus， Serratia， uncultured_bacterium_f_Burkholderiaceae， uncultured_bacterium_f_Enterobacteriaceae， uncultured_bacterium_f_Xanthobacteraceae） always stay in a low-level state.In the S system，the early stage of fermentation is the nitrite accumulation stage. At this time， Lactococcus and Leuconostoc are the dominant Lactobacillus， and Lactobacillus becomes the absolute dominant Lactobacillus at the middle and late stages of fermentation. In the early stage of fermentation， the pathogenic bacteria such as Acinetobacter， Enterococcus， Serratia and uncultured_bacterium_f_Enterobacteriaceae proliferates largely， leading to a large accumulation of nitrite， and the accumulated nitrite keeps the relative abundance of these pathogenic bacteria always at a low-level state， with the most significant effect on the growth of Serratia.

Key words： northeast sauerkraut; nitrite; bacterial flora; dominant Lactobacillus; pathogenic bacteria

收稿日期：2023-11-28

基金項(xiàng)目：國家重點(diǎn)研發(fā)計(jì)劃（2018YFC0311100）

作者簡介：張慶芳（1965—），女，教授，博士，研究方向：微生物及酶工程的基礎(chǔ)理論與應(yīng)用。

*通信作者：遲乃玉（1965—），男，教授，博士，研究方向：微生物及酶工程的基礎(chǔ)理論與應(yīng)用；

于爽（1986—），男，教授，博士，研究方向：微生物及酶工程的基礎(chǔ)理論與應(yīng)用。

東北酸菜（northeast sauerkraut）是一種以大白菜為原材料的發(fā)酵蔬菜，口味咸酸，有開胃健食、潤腸通便、殺菌抗炎、促進(jìn)消化的效果，營養(yǎng)價(jià)值豐富^[1]。在東北酸菜發(fā)酵期間，除乳酸菌之外的一些雜菌生長繁殖產(chǎn)生的硝酸鹽還原酶會將蔬菜本身含有的硝酸鹽還原成亞硝酸鹽^[2-3]，亞硝酸鹽是一種潛在的致癌物質(zhì)，同時(shí)也是一種常用的食品防腐劑。所以，控制發(fā)酵蔬菜食品中的亞硝酸鹽含量對于提高發(fā)酵蔬菜食品的安全性至關(guān)重要。但無法分析東北酸菜在發(fā)酵過程中產(chǎn)生的一定含量的亞硝酸鹽對細(xì)菌菌群的影響，從而不能得出兩者之間的明確關(guān)系。

隨著科學(xué)技術(shù)的不斷發(fā)展，人們已經(jīng)對發(fā)酵蔬菜中微生物多樣性開展了多方位的研究。Wu等^[4]利用DGGE分析東北酸菜發(fā)酵過程中的微生物多樣性，鑒定出14種細(xì)菌。Yang等^[5]對來自3個不同家庭的酸菜進(jìn)行了微生物群落動態(tài)分析，發(fā)現(xiàn)魏斯氏菌屬、乳桿菌屬、梭狀芽孢菌屬、腸桿菌屬是酸菜中菌群結(jié)構(gòu)的主要部分。劉長根^[6]通過高通量測序分析得出東北酸菜微生物中有17個門，35個綱，77個目，145個科，408個屬，283個種和1 414個OTU。Yang等^[7]研究發(fā)現(xiàn)蘿卜泡菜中細(xì)菌多樣性隨著發(fā)酵的進(jìn)行而顯著下降。Chen等^[8]、Gou等^[9]、Huang等^[10]、Liu等^[11]、Ren等^[12]、Shen等^[13]、Yun等^[14]對不同酒曲和原料發(fā)酵的酒類微生物的多樣性進(jìn)行了研究。Shan等^[15]研究表明，云南宣威火腿核心發(fā)酵菌屬為孢囊放線菌屬、乳酸桿菌、棒狀桿菌屬、小單孢菌屬、鏈霉菌。Lin等^[16]和Xiang等^[17]研究發(fā)現(xiàn)Lactobacillus能減少4種苦味氨基酸的含量，具有降低發(fā)酵蔬菜苦味的功效，還能提高發(fā)酵蔬菜的感官品質(zhì)和風(fēng)味。Liang等^[18]和Peng等^[19]也發(fā)現(xiàn)Lactobacillus參與碳水化合物、氨基酸、脂類的代謝，從而形成風(fēng)味物質(zhì)。Guo等^[20]利用宏基因組技術(shù)在白酒窖泥中發(fā)現(xiàn)乳酸菌、芽孢桿菌和梭菌等細(xì)菌與乳酸、乙酸等有機(jī)酸及乙醇、酯類等風(fēng)味化合物產(chǎn)生有關(guān)，而半乳酵母、青霉菌和曲霉菌等真菌微生物在白酒發(fā)酵中可以產(chǎn)生葡萄糖淀粉酶和α-淀粉酶，這對糖酵解過程影響顯著并可將淀粉等大分子物質(zhì)降解為小分子糖。雖然目前已有大量發(fā)酵蔬菜及其他發(fā)酵食品中微生物多樣性和呈味物質(zhì)的相關(guān)研究，但通過高通量測序結(jié)合OD_{600 nm}值和理化指標(biāo)進(jìn)行監(jiān)測進(jìn)而探究亞硝酸鹽對發(fā)酵蔬菜不同發(fā)酵階段細(xì)菌菌群影響的研究尚少。

本研究以東北酸菜為研究對象，采用一種東北酸菜小型發(fā)酵體系^[21]，即使用無菌且統(tǒng)一規(guī)格的礦泉水瓶作為發(fā)酵容器。將發(fā)酵溫度和腌漬方式進(jìn)行優(yōu)化后，確定25 ℃為本研究發(fā)酵體系中的最適發(fā)酵溫度。并基于Illumina NovaSeq測序平臺，采用16S rRNA高通量測序技術(shù)，對生漬法發(fā)酵條件下有、無添加亞硝酸鹽的東北酸菜發(fā)酵過程中不同發(fā)酵時(shí)間點(diǎn)的細(xì)菌菌群進(jìn)行分析；同時(shí)對細(xì)菌菌群生長（OD_{600 nm}）及生理代謝相關(guān)的主要理化指標(biāo)（pH、總酸、亞硝酸鹽）進(jìn)行監(jiān)測，進(jìn)而探究東北酸菜不同發(fā)酵階段亞硝酸鹽與菌群的關(guān)系。本研究旨在探究亞硝酸鹽對東北酸菜不同發(fā)酵階段菌群的影響，從而為東北酸菜及其他發(fā)酵蔬菜的開發(fā)提供更多的基礎(chǔ)理論依據(jù)，并消除了人們對東北酸菜中亞硝酸鹽安全性的疑慮。

1 材料與方法

1.1 材料與試劑

新鮮大白菜；硼砂、亞鐵氰化鉀、乙酸鋅、亞硝酸鈉、對氨基苯磺酸、鹽酸萘乙二胺、氫氧化鈉、酚酞、胰蛋白胨、酵母浸膏、葡萄糖、瓊脂（均為分析純）：上海麥克林生化科技股份有限公司；DP812土壤基因組脫氧核糖核酸（deoxyribonucleic acid，DNA）提取試劑盒：天根生化科技（北京）有限公司；16S rRNA基因V3～V4區(qū)引物合成及建庫測序：由北京百邁客生物科技有限公司完成。

1.2 儀器與設(shè)備

PHS-3E型pH計(jì) 上海儀電科學(xué)儀器股份有限公司；CRY-2112恒溫?fù)u床 上海茸研儀器有限公司；DK-S26電熱恒溫水浴鍋 上海精宏實(shí)驗(yàn)設(shè)備有限公司；Thermo Multiskan 1510酶標(biāo)儀 芬蘭Labsystems公司；AL204電子天平 梅特勒-托利多儀器（上海）有限公司；UV-1200型紫外可見分光光度計(jì) 上海美譜達(dá)儀器有限公司。

1.3 實(shí)驗(yàn)方法

1.3.1 酸菜的制作工藝流程及操作要點(diǎn)

制作工藝流程：挑選新鮮大白菜→清洗→瀝水→切絲混勻→裝瓶→注鹽水→封蓋→恒溫發(fā)酵→酸菜成品。

操作要點(diǎn)：設(shè)置添加亞硝酸鹽（亞硝酸鈉：100 mg/L）組和未添加亞硝酸鹽組，其中未添加亞硝酸鹽組為對照組。挑選新鮮大白菜，清理好外層葉片，將白菜順著菜幫掰成若干片后，用清水清洗干凈，瀝干表面水分后將白菜切成0.5～1 cm的均勻細(xì)絲。注意保證不同部位的白菜細(xì)絲混勻，選用統(tǒng)一規(guī)格的555 mL PET材料制備的東北酸菜發(fā)酵容器若干瓶，每瓶分裝160 g混勻的白菜絲，再注入1.5%鹽水至滿瓶，最終形成封閉發(fā)酵體系，將全部酸菜放置于25 ℃恒溫培養(yǎng)箱中發(fā)酵一定時(shí)間，得到酸菜成品。

1.3.2 酸菜理化指標(biāo)的測定

每瓶酸菜代表一個獨(dú)立的時(shí)間點(diǎn)，從1 h開始取樣測量，然后分別在13 h、25 h、73 h、15 d取樣一次。將添加亞硝酸鹽組中發(fā)酵1 h、13 h、25 h、73 h、15 d的樣本依次命名為YS1、YS2、YS3、YS4、YS5，YS1～YS5的分組命名為Y組；將未添加亞硝酸鹽組中發(fā)酵1 h、13 h、25 h、73 h、15 d的樣本依次命名為S1、S2、S3、S4、S5，S1～S5的分組命名為S組。

將瓶內(nèi)酸菜與發(fā)酵液倒入燒杯中，將其全部研磨后制成勻漿，分別稱取10 g酸菜勻漿作為待測樣品，進(jìn)行理化指標(biāo)的檢測^[22-24]。pH的測定：使用pH計(jì)；總酸含量的測定：參照GB 12456—2021《食品安全國家標(biāo)準(zhǔn) 食品中總酸的測定》中的酸堿滴定法；亞硝酸鹽含量的測定：參照GB 5009.33—2016《食品安全國家標(biāo)準(zhǔn) 食品中亞硝酸鹽與硝酸鹽的測定》中的鹽酸萘乙二胺法；OD_{600 nm}值的測定：采用紫外分光光度計(jì)測定在波長600 nm處的吸光度值。所有實(shí)驗(yàn)重復(fù)測定3次取平均值。

1.3.3 16S rRNA高通量測序技術(shù)

高通量測序技術(shù)（high-throughput sequencing）也稱“下一代”測序技術(shù)（“next-generation” sequencing technology，NGS），能夠一次性完成海量DNA分子的測序任務(wù)。如今，基于高通量測序的擴(kuò)增子測序技術(shù)已在微生物多樣性分析研究中普及，根據(jù)待測樣品的不同又分為16S、ITS、18S擴(kuò)增子測序，細(xì)菌主要是基于16S區(qū)，真菌主要基于18S區(qū)或ITS區(qū)（內(nèi)轉(zhuǎn)錄間隔區(qū)）^[25]。

按照土壤基因組DNA提取試劑盒說明書提取樣本DNA，以其為模板，采用引物對338F（5′-ACTCCTAC-GGGAGGCAGCA-3′）和806R（5′-GGACTACHVGGG-TWTCTAAT-3′）進(jìn)行PCR擴(kuò)增，PCR擴(kuò)增產(chǎn)物委托北京百邁克生物科技有限公司完成建庫測序。

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

首先利用Trimmomatic軟件和Cutadapt軟件對初始序列進(jìn)行質(zhì)量過濾；然后利用Usearch軟件和Uchime軟件對序列進(jìn)行拼接和篩選，得到最終有效數(shù)據(jù)。使用QIIME 2軟件分析樣本中的Alpha多樣性（Alpha diversity）和Beta多樣性（Beta diversity），并比較S組和Y組中菌群豐度的差異^[26-29]。

2 結(jié)果與分析

2.1 酸菜樣本理化指標(biāo)及OD_{600 nm}值的檢測

pH能反映酸菜的發(fā)酵進(jìn)程，總酸（TAN）含量能體現(xiàn)出酸菜發(fā)酵后的酸香味，亞硝酸鹽（NIT）含量反映了發(fā)酵蔬菜的安全性，OD_{600 nm}值在一定程度上可以反映酸菜中的菌群濃度^[30]。

由圖1可知，添加亞硝酸鹽與未添加亞硝酸鹽這兩種發(fā)酵條件下，pH均呈整體下降趨勢，但添加亞硝酸鹽組的pH整體大于未添加亞硝酸鹽組。由此可知，添加亞硝酸鹽會影響整個發(fā)酵過程中pH的變化，導(dǎo)致pH整體偏高。添加亞硝酸鹽組與未添加亞硝酸鹽組的TAN均呈整體上升趨勢，但添加亞硝酸鹽組的TAN含量整體小于未添加亞硝酸鹽組，由此可知，添加亞硝酸鹽會影響整個發(fā)酵過程中TAN含量的變化，導(dǎo)致TAN含量整體偏低，該結(jié)果與pH的變化情況相吻合。發(fā)酵時(shí)間為25 h時(shí)未添加亞硝酸鹽組達(dá)到“亞硝峰”，此時(shí)添加亞硝酸鹽組的NIT含量為33.50 mg/L，當(dāng)發(fā)酵時(shí)間為25～73 h時(shí)，添加亞硝酸鹽組的NIT含量下降趨勢更劇烈。添加亞硝酸鹽與未添加亞硝酸鹽這兩種發(fā)酵條件下OD_{600 nm}值均呈整體上升趨勢，但添加亞硝酸鹽組的OD_{600 nm}值整體小于未添加亞硝酸鹽組，由此可知，添加亞硝酸鹽影響了整個發(fā)酵過程中的菌群濃度，在相同時(shí)間點(diǎn)抑制了發(fā)酵系統(tǒng)內(nèi)菌群的生長。

2.2 酸菜樣本的生物信息學(xué)分析

2.2.1 酸菜樣本細(xì)菌微生物高通量測序結(jié)果及質(zhì)量評估

S1～S5、YS1～YS5樣本細(xì)菌微生物的高通量測序結(jié)果見表1。

由表1可知，從兩個酸菜樣本中共獲得891 634條原始序列，雙端序列質(zhì)控、拼接后共得到888 590條高質(zhì)量序列，最終將高質(zhì)量序列通過拼接、過濾和嵌合后共得到有效序列839 128條，序列平均堿基長度為423 bp，樣本平均GC含量為52.70%。此外，各酸菜樣本的測序質(zhì)量值Q20＞98%、Q30＞95%，有效序列占比＞91%，說明測序結(jié)果數(shù)據(jù)質(zhì)量良好。

2.2.2 酸菜樣本細(xì)菌菌群的Alpha多樣性分析

Alpha多樣性主要的衡量指標(biāo)包括Chao 1指數(shù)、ACE指數(shù)、Shannon指數(shù)、Simpson指數(shù)、覆蓋率^[31]。Chao 1和ACE指數(shù)衡量物種的豐度，指數(shù)越大，豐度越高；Simpson指數(shù)和Shannon指數(shù)反映了物種的多樣性，Simpson指數(shù)越小，Shannon指數(shù)越大，多樣性越高^[32]。兩個酸菜樣本細(xì)菌菌群的Alpha多樣性指數(shù)分析結(jié)果見表2。

由表2可知，樣本S1和S2的物種豐度均高于樣本YS1和YS2，多樣性均低于YS1和YS2，樣本S3、S5的物種豐度和多樣性均高于樣本YS3，樣本S4的物種豐度低于YS4、多樣性高于YS4。兩個樣本的覆蓋率均≥0.999 4，表示測序覆蓋率達(dá)到了99%，樣本中物種幾乎全部被檢測出來。根據(jù)10個樣本的理化指標(biāo)進(jìn)行推斷：第一，在未添加亞硝酸鹽的樣本中，S2達(dá)到了整個發(fā)酵周期菌群多樣性峰值，此時(shí)酸菜中細(xì)菌豐度最大；隨著發(fā)酵進(jìn)程的推進(jìn)，S4中細(xì)菌的多樣性和豐度值逐漸減少至平穩(wěn)。第二，在添加亞硝酸鹽的樣本中，隨著發(fā)酵進(jìn)程的推進(jìn)，YS5達(dá)到了整個發(fā)酵周期菌群多樣性峰值，此時(shí)酸菜中細(xì)菌豐度最大。由此可知，在S組中S2的物種豐度最高，S1的物種多樣性最高；在Y組中YS5的物種豐度最高，YS1的物種多樣性最高。發(fā)現(xiàn)無論是否添加亞硝酸鹽，在發(fā)酵1 h即發(fā)酵初始時(shí)的物種多樣性均最高。S組內(nèi)S1～S3的物種豐度大于Y組內(nèi)YS1～YS3，分析原因?yàn)樘砑觼喯跛猁}抑制了發(fā)酵過程中多種細(xì)菌菌群的生長繁殖。

2.2.3 酸菜樣本細(xì)菌菌群的Beta多樣性分析

本研究采用主坐標(biāo)分析（principal coordinates analysis，PCoA）和非度量多維尺度分析（non-metric multidimensional scaling，NMDS）對S組和Y組間細(xì)菌菌群結(jié)構(gòu)的差異進(jìn)行分析。PCoA將標(biāo)記出的特征值和向量統(tǒng)一排序，在降維的思想中選擇排列在前幾名的主要特征值，基于加權(quán)UniFrac距離和未加權(quán)UniFrac距離的計(jì)算結(jié)果找到距離矩陣中最主要的坐標(biāo)，以觀察兩組間的差異^[33]。NMDS分析是一種非線性模型，其基本特征是將樣本間的相似或相異性數(shù)值看作一種以點(diǎn)之間距離來表示的單調(diào)函數(shù)，在保持初始數(shù)據(jù)排序關(guān)系的前提下，用相同排序的新數(shù)值代替初始數(shù)據(jù)來完成度量多維尺度分析，從而分辨樣本間的差異性^[34]。在兩種分析方法中，樣本在坐標(biāo)圖上距離越近，表明相似性越高，樣本的聚集程度越高，細(xì)菌菌群結(jié)構(gòu)越相似；樣本離散程度越大，細(xì)菌菌群結(jié)構(gòu)差異越大。

由圖2可知，第一主成分PC1和第二主成分PC2的貢獻(xiàn)率分別為39.67%和20.21%，表明PCoA可以高度解釋原始測序信息的差異。S組內(nèi)僅S1分布于PC1的負(fù)坐標(biāo)區(qū)域，組內(nèi)其他樣本均分布于PC1的正坐標(biāo)區(qū)域，而Y組內(nèi)5個樣本均分布于PC1的負(fù)坐標(biāo)區(qū)域。說明發(fā)酵1 h時(shí)細(xì)菌菌群結(jié)構(gòu)與其他樣本的差異最大，分析原因?yàn)榘l(fā)酵環(huán)境內(nèi)原料、水?dāng)y帶的菌群導(dǎo)致此時(shí)菌群結(jié)構(gòu)復(fù)雜多樣。

由圖3可知，Y組整體的離散程度大于S組。在S1與YS1、S2與YS2、S3與YS3、S4與YS4、S5與YS5 5組樣本的間距中，發(fā)現(xiàn)S1與YS1的間距最近，說明隨著發(fā)酵時(shí)間的延長，添加亞硝酸鹽對發(fā)酵體系內(nèi)菌群結(jié)構(gòu)的影響越來越顯著。

2.3 東北酸菜發(fā)酵系統(tǒng)中細(xì)菌菌群與亞硝酸鹽長消變化關(guān)系

2.3.1 發(fā)酵系統(tǒng)中優(yōu)勢乳酸菌菌屬與亞硝酸鹽長消變化關(guān)系

由圖4可知，隨著發(fā)酵進(jìn)程的推進(jìn)，發(fā)酵體系內(nèi)的優(yōu)勢乳酸菌菌屬的相對豐度不斷升高，發(fā)酵過程中乳酸菌不斷產(chǎn)生亞硝酸鹽還原酶和H⁺，導(dǎo)致TAN含量不斷升高，pH不斷降低，使亞硝酸鹽不斷降解至安全值；乳酸菌在發(fā)酵過程中逐漸占據(jù)主導(dǎo)地位，在發(fā)酵中后期成為發(fā)酵體系內(nèi)的優(yōu)勢菌群，導(dǎo)致OD_{600 nm}值迅速升高^[35]。在發(fā)酵時(shí)間點(diǎn)2以后，Y組中優(yōu)勢乳酸菌的相對豐度遠(yuǎn)遠(yuǎn)高于S組，說明添加的亞硝酸鹽影響了Y組中優(yōu)勢乳酸菌的生長，使其相對豐度不斷增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，S組與Y組中優(yōu)勢乳酸菌的相對豐度有所降低但都出現(xiàn)回升現(xiàn)象。由此可知，亞硝酸鹽對優(yōu)勢乳酸菌的生長具有明顯影響，使其在發(fā)酵中后期相對豐度增加。

由圖5可知，在發(fā)酵時(shí)間點(diǎn)2以后，Y組中Lactococcus、Leuconostoc的相對豐度遠(yuǎn)遠(yuǎn)高于S組，說明添加的亞硝酸鹽影響了Y組中Lactococcus、Leuconostoc的生長，使其相對豐度不斷增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，S組與Y組中Lactococcus的相對豐度開始降低，而Y組中Leuconostoc的相對豐度有所降低但又出現(xiàn)回升現(xiàn)象，S組中Leuconostoc的相對豐度產(chǎn)生波動后降低。在發(fā)酵過程中，Y組中Lactobacillus的相對豐度始終低于S組，說明添加的亞硝酸鹽影響了Y組中Lactobacillus的生長，使其相對豐度始終處于低水平狀態(tài)；在S組中，發(fā)酵系統(tǒng)內(nèi)不斷生成的亞硝酸鹽影響了Lactobacillus的生長，使其相對豐度先升高后降低，同時(shí)Lactobacillus迎來了第一個低峰值；當(dāng)?shù)竭_(dá)發(fā)酵時(shí)間點(diǎn)4時(shí)，由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Lactobacillus的相對豐度不斷增加，并將迎來第二個更高的峰值。結(jié)果表明，亞硝酸鹽對Lactococcus、Leuconostoc和Lactobacillus的生長具有明顯影響，使Lactococcus、Leuconostoc在發(fā)酵中后期相對豐度增加；使Lactobacillus的相對豐度始終處于低水平狀態(tài)。

2.3.2 發(fā)酵系統(tǒng)中α-變形菌菌屬與亞硝酸鹽長消變化關(guān)系

由圖6可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中α-變形菌的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中α-變形菌的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，α-變形菌的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對α-變形菌的生長具有明顯影響，使其在發(fā)酵初期相對豐度升高，在發(fā)酵中后期相對豐度不斷降低。

由圖7可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、Aureimonas、Brevundimonas、Methylobacterium、Sphingomonas的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、Aureimonas、Brevundimonas、Methylobacterium、Sphingomonas的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、Aureimonas、

Brevundimonas、Methylobacterium、Sphingomonas的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、

Aureimonas、Brevundimonas、Methylobacterium、Sphingomonas

的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

2.3.3 發(fā)酵系統(tǒng)中β-變形菌菌屬與亞硝酸鹽長消變化關(guān)系

由圖8可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中β-變形菌的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了β-變形菌的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，β-變形菌的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對β-變形菌的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

由圖9可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中Acidovorax、Aquabacterium、Delftia、Herbaspirillum、Methylophilus、Methylovorus的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中Acidovorax、Aquabacterium、Delftia、Herbaspirillum、Methylophilus、Methylovorus的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Acidovorax、Aquabacterium、Delftia、Herbaspirillum、Methylophilus、Methylovorus的相對豐度不斷降低，并處于低水平狀態(tài)。隨著發(fā)酵進(jìn)程的推進(jìn)，S組與Y組中Comamonas的相對豐度呈現(xiàn)基本相同的變化趨勢，但除發(fā)酵時(shí)間點(diǎn)3～4之間，Y組中Comamonas的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中Comamonas的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Comamonas的相對豐度不斷降低，雖有所回升但仍處于低水平狀態(tài)。由此可知，亞硝酸鹽對Acidovorax、Aquabacterium、Delftia、Herbaspirillum、Methylophilus、Methylovorus、Comamonas的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

2.3.4 發(fā)酵系統(tǒng)中CFB類菌屬與亞硝酸鹽長消變化關(guān)系

由圖10可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中CFB類菌屬的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中CFB類菌屬的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，CFB類菌屬的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對CFB類菌屬的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

由圖11可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中Dyadobacter、Pedobacter、Sphingobacterium的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中Dyadobacter、Pedobacter、Sphingobacterium的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Dyadobacter、Pedobacter、Sphingobacterium的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對Dyadobacter、Pedobacter、Sphingobacterium的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

2.3.5 發(fā)酵系統(tǒng)中未培養(yǎng)菌屬與亞硝酸鹽長消變化關(guān)系

由圖12可知，隨著發(fā)酵進(jìn)程的推進(jìn)，S組與Y組中未培養(yǎng)菌屬的相對豐度都呈現(xiàn)下降趨勢，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中未培養(yǎng)菌屬的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中未培養(yǎng)菌屬的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，未培養(yǎng)菌屬的相對豐度不斷降低，并處于低水平狀態(tài)。由此可知，亞硝酸鹽對未培養(yǎng)菌屬的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

由圖13可知，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast的相對豐度明顯高于S組，說明添加的亞硝酸鹽影響了Y組中Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast的相對豐度不斷降低，并處于低水平狀態(tài)。由此可得，亞硝酸鹽對Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast的生長具有明顯影響，使其在發(fā)酵初期相對豐度增加。

2.3.6 發(fā)酵系統(tǒng)中致病菌菌屬與亞硝酸鹽長消變化關(guān)系

由圖14可知，除發(fā)酵初始外，Y組中致病菌的相對豐度始終低于S組，說明添加的亞硝酸鹽影響了Y組中致病菌的生長，使其相對豐度始終處于低水平狀態(tài)。在S組中，發(fā)酵系統(tǒng)內(nèi)不斷生成的亞硝酸鹽影響了致病菌的生長，使其相對豐度先升高后降低，同時(shí)致病菌迎來了第一個低峰值；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，S組中致病菌的相對豐度仍不斷降低。這意味著無論是Y組還是S組，在發(fā)酵后期致病菌菌屬的相對豐度都處于低水平狀態(tài)，此時(shí)東北酸菜不會對人體健康構(gòu)成危害，為可食用狀態(tài)。由此可知，亞硝酸鹽對致病菌的生長具有明顯影響，使其相對豐度始終處于低水平狀態(tài)。

由圖15可知，在發(fā)酵時(shí)間點(diǎn)1～2之間，Y組中Alkanindiges的相對豐度明顯高于S組，在發(fā)酵時(shí)間點(diǎn)1～3之間，Y組中uncultured_bacterium_f_Burkholderiaceae的相對豐度高于S組，說明添加的亞硝酸鹽影響了Y組中Alkanindiges、uncultured_bacterium_f_Burkholderiaceae的生長，使其相對豐度增加；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，Y組中Alkanindiges、uncultured_bacterium_f_Burkholderiaceae的相對豐度一直處于低水平狀態(tài)。Y組中Acinetobacter、Chryseobacterium的相對豐度始終處于低水平狀態(tài)，且在發(fā)酵時(shí)間點(diǎn)1～3之間Acinetobacter、Chryseobacterium的相對豐度遠(yuǎn)遠(yuǎn)低于S組，說明添加的亞硝酸鹽影響了Y組中Acinetobacter、Chryseobacterium的生長，使其相對豐度始終處于低水平狀態(tài)；在S組中，發(fā)酵系統(tǒng)內(nèi)不斷生成的亞硝酸鹽影響了S組中Acinetobacter、Chryseobacterium的生長，使其相對豐度先升高后降低；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，導(dǎo)致乳酸菌大量積累，S組中Acinetobacter、Chryseobacterium的相對豐度仍處于低水平狀態(tài)。在整個發(fā)酵過程中，Y組中Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae的相對豐度處于低水平狀態(tài)且明顯低于S組，說明添加的亞硝酸鹽顯著影響了Y組中Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae的生長，使其相對豐度始終處于低水平狀態(tài)；由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，導(dǎo)致乳酸菌大量積累，Enterococcus的相對豐度仍處于低水平狀態(tài)；在發(fā)酵時(shí)間點(diǎn)1～2之間，分析S組中的Serratia由于具有硝酸鹽還原酶，而促進(jìn)了亞硝酸鹽的生成；在發(fā)酵時(shí)間點(diǎn)2～3之間，該階段表現(xiàn)為亞硝酸鹽影響了Serratia的生長，使其相對豐度降低，由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，乳酸菌大量積累，Serratia的相對豐度有所回升，但其相對豐度始終處于處于低水平狀態(tài)；uncultured_bacterium_f_Enterobacteriaceae是腸桿菌科中具有硝酸鹽還原酶的一種細(xì)菌，可將發(fā)酵體系中存在的硝酸鹽還原成亞硝酸鹽，在S組發(fā)酵前期，發(fā)酵體系內(nèi)含氧量較高，uncultured_bacterium_f_Enterobacteriaceae不斷消耗發(fā)酵系統(tǒng)中的氧氣并隨之不斷增長，同時(shí)將原料所攜帶的硝酸鹽不斷地轉(zhuǎn)化成亞硝酸鹽，導(dǎo)致發(fā)酵體系內(nèi)亞硝酸鹽含量不斷升高，隨著發(fā)酵系統(tǒng)中乳酸菌大量積累、pH不斷降低、TAN不斷升高，發(fā)酵體系內(nèi)呈現(xiàn)出缺氧的酸性環(huán)境，以Lactobacillus為代表的乳酸菌在發(fā)酵體系中逐漸占據(jù)主導(dǎo)地位，使得uncultured_bacterium_f_Enterobacteriaceae的相對豐度在發(fā)酵后期不斷下降，并處于低水平狀態(tài)。在發(fā)酵時(shí)間點(diǎn)1～2之間，Y組uncultured_bacterium_f_Xanthobacteraceae的相對豐度明顯低于S組，說明添加的亞硝酸鹽顯著影響了Y組中uncultured_bacterium_f_Xanthobacteraceae的生長，使其相對豐度降低；在發(fā)酵時(shí)間點(diǎn)2～3之間，發(fā)酵過程中不斷生成亞硝酸鹽，從而使uncultured_bacterium_f_Xanthobacteraceae的相對豐度繼續(xù)降低，同時(shí)其迎來了第一個低峰值，由于發(fā)酵系統(tǒng)內(nèi)亞硝酸鹽不斷降解，TAN含量不斷升高，pH不斷降低，uncultured_bacterium_f_Xanthobacteraceae的相對豐度有所回升，但仍處于低水平狀態(tài)。由此可知，亞硝酸鹽對Alkanindiges、uncultured_bacterium_f_Burkholderiaceae、Acinetobacter、Chryseobacterium、Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae、uncultured_bacterium_f_Xanthobacteraceae的生長具有明顯影響，使Alkanindiges、uncultured_bacterium_f_Burkholderiaceae在發(fā)酵初期相對豐度增加；使Acinetobacter、Chryseobacterium、Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae的相對豐度始終處于低水平狀態(tài)；使uncultured_bacterium_f_Xanthobacteraceae在發(fā)酵前期相對豐度降低。

3 討論

圖1展示了添加亞硝酸鹽組（Y組）與未添加亞硝酸鹽組（S組）中pH、TAN、NIT和OD_{600 nm}值的變化情況，發(fā)現(xiàn)Y組的pH始終大于S組，Y組的NIT含量始終小于S組，該現(xiàn)象與張慶芳等^[36]在亞硝酸鹽對Lactobacillus brevis 4903發(fā)酵中亞硝酸菌對乳酸菌的發(fā)酵有抑制作用的研究結(jié)果一致。Y組發(fā)酵體系內(nèi)產(chǎn)酸量較低，推測原因可能是乳酸菌產(chǎn)生亞硝酸還原酶將亞硝酸鹽還原成氨氣，氨氣與乳酸菌所產(chǎn)的乳酸反應(yīng)，導(dǎo)致發(fā)酵環(huán)境中pH偏高，TAN含量偏低，有利于發(fā)酵早期乳球菌屬的生長與積累。

Alpha多樣性結(jié)果表明，添加亞硝酸鹽對各樣本中共有的菌屬的相對豐度有顯著影響。Beta多樣性結(jié)果表明，發(fā)酵時(shí)間點(diǎn)1時(shí)，S組內(nèi)S1的菌群結(jié)構(gòu)與組內(nèi)其他樣本的差異最大，隨著發(fā)酵時(shí)間的延長，添加亞硝酸鹽對發(fā)酵體系內(nèi)菌群結(jié)構(gòu)的影響越來越顯著。

本文研究表明S組和Y組中共含有Enterococcus、Lactococcus、Leuconostoc和Lactobacillus 4種乳酸菌菌屬，其中Lactococcus、Leuconostoc和Lactobacillus為S組中的優(yōu)勢乳酸菌菌屬。與S組相比，Y組中Lactococcus與Leuconostoc的相對豐度分別增加了18.55%和8.05%，而Lactobacillus的相對豐度減少了14.79%，說明添加亞硝酸鹽使Lactococcus與Leuconostoc的相對豐度升高，但使Lactobacillus的相對豐度降低；隨著發(fā)酵進(jìn)程的推進(jìn)，亞硝酸鹽被大量降解，亞硝酸鹽對Lactobacillus的發(fā)酵影響開始解除。發(fā)酵前期為亞硝酸鹽積累階段，發(fā)酵時(shí)間點(diǎn)2時(shí)亞硝酸鹽含量達(dá)到12.62 mg/L，此時(shí)Lactococcus、Leuconostoc為S組中的優(yōu)勢乳酸菌菌屬，而當(dāng)發(fā)酵25 h時(shí)達(dá)到亞硝峰（32.63 mg/L）后，Lactobacillus成為S組中的絕對優(yōu)勢乳酸菌菌屬；Lactococcus和Leuconostoc為Y組中的優(yōu)勢乳酸菌菌屬，且亞硝酸鹽使Leuconostoc和Lactococcus的相對豐度升高，使Lactobacillus的相對豐度降低，此結(jié)果與Huang等^[37]的研究結(jié)果相同。在發(fā)酵73 h時(shí)，Y組中亞硝酸鹽含量達(dá)到23.33 mg/L，為S組中亞硝酸鹽含量（6.66 mg/L）的3.5倍；在發(fā)酵15 d時(shí)，Y組中亞硝酸鹽含量達(dá)到11.80 mg/L，為S組亞硝酸鹽含量（3.09 mg/L）的3.8倍，由此可知，在發(fā)酵73 h后Y組中亞硝酸鹽對Lactobacillus的生長仍具有影響，而S組中亞硝酸鹽基本不影響Lactobacillus的生長。乳酸菌具有抗菌、抗氧化、免疫、調(diào)節(jié)人體脂代謝、治療高血壓等功效^[38]。Chryseobacterium、Enterococcus、Serratia和uncultured_bacterium_f_Enterobacteriaceae是東北酸菜中比較有代表性的致病菌菌屬，觀察S組和Y組中幾種致病菌菌屬的變化情況，發(fā)現(xiàn)YS2中Chryseobacterium的相對豐度（15.94%）比S2（32.93%）低16.99%，YS2中Enterococcus的相對豐度（0.22%）比S2（0.47%）低0.25%，YS2中uncultured_bacterium_f_Enterobacteriaceae的相對豐度（0.35%）比S2（28.31%）低27.96%，YS2中Serratia的相對豐度（0.01%）比S2（0.28%）低0.27%，其中Serratia在YS3、YS4、YS5中相對豐度均為0，說明亞硝酸鹽對Serratia的生長影響最顯著，使其相對豐度始終處于極低水平狀態(tài)。以上均證明亞硝酸鹽可防止發(fā)酵過程中出現(xiàn)過多雜菌引起腐敗現(xiàn)象，保證了發(fā)酵的順利進(jìn)行。此外，本研究發(fā)現(xiàn)無論是否人為添加亞硝酸鹽，東北酸菜中亞硝酸鹽的殘留量均會隨著時(shí)間的延長降至遠(yuǎn)低于國標(biāo)值，因此在科學(xué)的生產(chǎn)工藝下，東北酸菜不會對人體健康產(chǎn)生危害。

東北酸菜的原料大白菜中含有一定量的硝酸鹽，在東北酸菜的發(fā)酵過程中，發(fā)酵前期存在多種具有硝酸還原酶的細(xì)菌，白菜自身所攜帶的硝酸鹽會被發(fā)酵體系內(nèi)的雜菌還原成亞硝酸鹽，亞硝酸鹽最終又會被乳酸菌逐漸降解至安全值^[39]。由此可知，食用科學(xué)工藝所生產(chǎn)的東北酸菜具有一定的安全性。由本研究中的兩個發(fā)酵系統(tǒng)中各菌屬相對豐度變化情況可知，在東北酸菜的整個發(fā)酵周期中，發(fā)酵前期pH較高，發(fā)酵系統(tǒng)內(nèi)Lactococcus大量增長，其產(chǎn)生的亞硝酸鹽還原酶使發(fā)酵體系內(nèi)的亞硝酸鹽大量降解；發(fā)酵后期Lactobacillus大量增長，其產(chǎn)生大量的乳酸導(dǎo)致發(fā)酵系統(tǒng)內(nèi)pH始終處于較低水平狀態(tài)，此時(shí)發(fā)酵體系內(nèi)的H⁺發(fā)揮作用導(dǎo)致亞硝酸鹽不斷降解。以上結(jié)果與張慶芳等^[35]、王一茜等^[40]的研究結(jié)果相符合。

4 結(jié)論

以25 ℃生漬法發(fā)酵的東北酸菜發(fā)酵液為研究對象，在添加亞硝酸鹽與未添加亞硝酸鹽這兩組中選取的10個樣本進(jìn)行高通量測序，得到以下結(jié)論：10個樣本共注釋得到13個門的142個屬，Y組中YS4的物種豐度大于S4，S組中S1的物種多樣性大于YS1。通過研究添加亞硝酸鹽組與未添加亞硝酸鹽組這兩種發(fā)酵條件下各理化指標(biāo)和微生物指標(biāo)變化情況，由pH、TAN含量、NIT含量和OD_{600 nm}值的變化可知，添加亞硝酸鹽會對酸菜整個發(fā)酵體系內(nèi)的菌群結(jié)構(gòu)產(chǎn)生一定的影響。綜上所述，無論是否添加亞硝酸鹽，兩個發(fā)酵體系中酸菜樣本的細(xì)菌總數(shù)均呈現(xiàn)先上升后下降的趨勢。發(fā)酵開始后，添加亞硝酸鹽對整個發(fā)酵周期中細(xì)菌的生長都具有影響，亞硝酸鹽使優(yōu)勢乳酸菌菌屬（Lactococcus、Leuconostoc、Lactobacillus）、α-變形菌菌屬（Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium、Aureimonas、Brevundimonas、Methylobacterium、Sphingomonas）、β-變形菌菌屬（Acidovorax、Aquabacterium、Comamonas、Delftia、Herbaspirillum、Methylophilus、Methylovorus）、CFB類菌屬（Dyadobacter、Pedobacter、Sphingobacterium）和未培養(yǎng)菌屬（Brassica_napus_rape、uncultured_bacterium_f_Chitinophagaceae、uncultured_bacterium_f_env.OPS_17、uncultured_bacterium_o_Chloroplast）在發(fā)酵初期相對豐度增加；使致病菌菌屬（Acinetobacter、Alkanindiges、Chryseobacterium、Enterococcus、Serratia、uncultured_bacterium_f_Burkholderiaceae、uncultured_bacterium_f_Enterobacteriaceae、uncultured_bacterium_f_Xanthobacteraceae）的相對豐度始終處于低水平狀態(tài)。東北酸菜發(fā)酵過程中乳酸菌逐漸占據(jù)主導(dǎo)地位，在發(fā)酵中后期成為發(fā)酵體系內(nèi)的優(yōu)勢菌群；在S系統(tǒng)內(nèi)，發(fā)酵前期為亞硝酸鹽積累階段，此時(shí)Lactococcus、Leuconostoc為優(yōu)勢乳酸菌菌屬，而達(dá)到發(fā)酵中后期時(shí)Lactobacillus成為絕對優(yōu)勢乳酸菌菌屬；在發(fā)酵初期，Acinetobacter、Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae等致病菌大量繁殖，導(dǎo)致亞硝酸鹽大量積累，而積累的亞硝酸鹽使Acinetobacter、Enterococcus、Serratia、uncultured_bacterium_f_Enterobacteriaceae等致病菌的相對豐度始終處于低水平狀態(tài)，其中對Serratia的生長影響最顯著。

5 展望

發(fā)酵蔬菜因其獨(dú)特的風(fēng)味深受人們喜愛，但其中的亞硝酸鹽含量一直是食品安全問題的熱點(diǎn)。本文對東北酸菜發(fā)酵過程中相對豐度排名前30的菌屬與亞硝酸鹽長消變化關(guān)系進(jìn)行了研究，為后期調(diào)控東北酸菜發(fā)酵環(huán)境中菌群結(jié)構(gòu)、提高東北酸菜發(fā)酵品質(zhì)奠定了基礎(chǔ)，為含有亞硝酸鹽的不同發(fā)酵蔬菜的開發(fā)與研究提供了新的理論依據(jù)和方向。

參考文獻(xiàn)：

[1]FONSECA S， RIVAS I， ROMAGUERA D， et al. Association between consumption of fermented vegetables and COVID-19 mortality at a country level in Europe[J].MedRxiv，2020（7）：

20147025.

[2]肖珍，謝筆鈞，孫智達(dá).蓮房原花青素對紫甘藍(lán)泡菜亞硝酸鹽的抑制作用[J].食品工業(yè)科技，2018，39（11）：67-73.

[3]燕平梅，王炫月，趙文婧.泡菜中亞硝酸鹽形成、還原相關(guān)微生物的研究[J].中國調(diào)味品，2020，45（9）：77-80.

[4]WU R N， YU M L， LIU X Y， et al. Changes in flavour and microbial diversity during natural fermentation of suan-cai， a traditional food made in Northeast China[J].International Journal of Food Microbiology，2015，211：23-31.

[5]YANG X Z， HU W Z， XIU Z L， et al. Microbial dynamics and volatilome profiles during the fermentation of Chinese northeast sauerkraut by Leuconostoc mesenteroides ORC2 and Lactobacillus plantarum HBUAS51041 under different salt concentrations[J].Food Research International，2020，130：108926.

[6]劉長根.我國傳統(tǒng)發(fā)酵蔬菜微生物多樣性比較[D].南昌：南昌大學(xué)，2019.

[7]YANG Y L， FAN Y， LI T， et al. Microbial composition and correlation between microbiota and quality-related physiochemical characteristics in Chongqing radish paocai[J].Food Chemistry，2022，369：130897.

[8]CHEN G M， HUANG Z R， WU L， et al. Microbial diversity and flavor of Chinese rice wine （Huangjiu）： an overview of current research and future prospects[J].Current Opinion in Food Science，2021，77：37-43.

[9]GOU M， WANG H， YUAN H， et al. Characterization of the microbial community in three types of fermentation starters used for Chinese liquor production[J].Journal of the Institute of Brewing，2015，121（4）：620-627.

[10]HUANG Z R， HONG J L， XU J X， et al. Exploring core functional microbiota responsible for the production of volatile flavour during the traditional brewing of Wuyi Hong Qu glutinous rice wine[J].Food Microbiology，2018，76：487-496.

[11]LIU Z B， WANG Z Y， SUN J Y， et al. The dynamics of volatile compounds and their correlation with the microbial succession during the traditional solid-state fermentation of Gutian Hong Qu glutinous rice wine[J].Food Microbiology，2020，86：103347.

[12]REN Q， SUN L， WU H， et al. The changes of microbial community and flavor compound in the fermentation process of Chinese rice wine using Fagopyrum tataricum grain as feedstock[J].Scientific Reports，2019，9（1）：3365.

[13]SHEN X， YAN X， XU J. Fatty acid profiles in muscle of large yellow croakers （Larimichthys crocea） can be used to distinguish between the samples of Dai-qu stock and Min-yuedong stock[J].Biochemical Systematics and Ecology，2018，77：37-43.

[14]YUN J H， KIM J H， LEE J E. Surface film formation in static-fermented rice vinegar： a case study[J].Mycobiology，2019，47（2）：250-255.

[15]SHAN L， LI Y， ZHENG S， et al.Analysis of the bacterial floral structure and diversity of Xuanwei ham by 16S rDNA sequencing[J].Journal of Food Safety，2020，40（4）：12800.

[16]LIN L J， ZENG J， TIAN Q M， et al.Effect of the bacterial community on the volatile flavour profile of a Chinese fermented condiment—red sour soup during fermentation[J].Food Research International，2022，155：111059.

[17]XIANG W， ZHANG N， LU Y， et al. Effect of Weissella cibaria co-inoculation on the quality of Sichuan pickle fermented by Lactobacillus plantarum[J].LWT-Food Science and Technology，2020，121：108975.

[18]LIANG H， CHEN H， JI C， et al. Dynamic and functional characteristics of predominant species in industrial Paocai as revealed by combined DGGE and metagenomic sequencing[J].Frontiers in Microbiology，2018，9：2416.

[19]PENG Q， JIANG S， CHEN J， et al. Unique microbial diversity and metabolic pathway features of fermented vegetables from Hainan， China[J].Frontiers in Microbiology，2018，9：399.

[20]GUO M Y， HOU C J， BIAN M H， et al. Characterization of microbial community profiles associated with quality of Chinese strong-aromatic liquor through metagenomics[J].Journal of Applied Microbiology，2019，127（3）：750-762.

[21]周姝靜，孫全敏，遲乃玉，等.東北酸菜發(fā)酵前后期細(xì)菌菌群多樣性分析[J].中國釀造，2022，41（5）：42-46.

[22]遲雪梅，王一茜，喬慧，等.發(fā)酵蔬菜硝酸鹽、亞硝酸鹽消長變化及其相關(guān)性的研究[J].食品與發(fā)酵工業(yè)，2018，44（1）：25-30.

[23]叢敏.酸菜自然發(fā)酵過程中化學(xué)成分及乳酸菌變化規(guī)律的研究[D].沈陽：沈陽農(nóng)業(yè)大學(xué)，2016.

[24]李欣蔚.東北地區(qū)自然發(fā)酵酸菜中乳酸菌時(shí)空分布規(guī)律的研究[D].沈陽：沈陽農(nóng)業(yè)大學(xué)，2017.

[25]FILIPPIS F D， PARENTE E， ZOTTA T， et al. A comparison of bioinformatic approaches for 16S rRNA gene profiling of food bacterial microbiota[J].International Journal of Food Microbiology，2018，265：9-17.

[26]PARKS D H， TYSON G W， HUGENHOLTZ P， et al. STAMP： statistical analysis of taxonomic and functional profiles[J].Bioinformatics，2014，30（21）：3123-3124.

[27]SUL W J， COLE J R， JESUS E D， et al. Bacterial community comparisons by taxonomy-supervised analysis independent of sequence alignment and clustering[J].Proceedings of the National Academy of Sciences，2011，108（35）：14637-14642.

[28]陳浩，何國慶.乳酸菌對發(fā)酵豆豉品質(zhì)的影響[J].食品安全質(zhì)量檢測學(xué)報(bào)，2017，8（12）：4516-4523.

[29]張慶芳，楊超，于爽，等.黃海海域海洋沉積物細(xì)菌多樣性分析[J].微生物學(xué)通報(bào)，2020，47（2）：370-378.

[30]高鳳.東北酸菜發(fā)酵過程亞硝酸鹽微生物酶降解機(jī)制及應(yīng)用[D].鎮(zhèn)江：江蘇大學(xué)，2019.

[31]王琪琪，田界先，潘宗東，等.基于Illumina MiSeq分析貴州凱里酸湯獨(dú)特風(fēng)味的優(yōu)勢菌群[J].食品與發(fā)酵工業(yè)，2020，46（14）：40-47.

[32]周藝萍，熊智，李選文，等.鹽分對新平酸腌菜主發(fā)酵期細(xì)菌多樣性的影響[J].中國釀造，2021，40（4）：26-32.

[33]孫煒寧.基于高通量測序的瀘酒酒醅和酸菜中細(xì)菌區(qū)系的研究[D].天津：天津大學(xué)，2017.

[34]RIVAS M N， BURTON O T， WISE P， et al. A microbiota signature associated with experimental food allergy promotes allergic sensitization and anaphylaxis[J].Journal of Allergy and Clinical Immunology，2013，131（1）：201-212.

[35]張慶芳，遲乃玉，鄭燕，等.乳酸菌降解亞硝酸鹽機(jī)理的研究[J].食品與發(fā)酵工業(yè)，2002（8）：27-31.

[36]張慶芳，遲乃玉，薛景珍，等.亞硝酸鹽影響Lactobacillus brevis 4903發(fā)酵的研究[J].微生物學(xué)雜志，2005（3）：25-27.

[37]HUANG Y， LIU D， JIA X， et al. Whole genome sequencing of Lactobacillus plantarum DMDL 9010 and its effect on growth phenotype under nitrite stress[J].LWT-Food Science and Technology，2021，149：111778.

[38]孫鈺薇，張?jiān)姮帲瑒⒅炯眩?發(fā)酵食品中乳酸菌的健康功效研究進(jìn)展[J].食品與發(fā)酵工業(yè)，2021，47（23）：280-287.

[39]遲雪梅，張慶芳.食品中常見乳酸菌高效降解NO₂^-發(fā)酵性能評價(jià)[J].食品與發(fā)酵工業(yè)，2017，43（6）：78-84.

[40]王一茜，榮金誠，王曉輝，等.常見乳酸菌降解亞硝酸鹽機(jī)理探討[J].食品與發(fā)酵工業(yè)，2019，45（8）：50-56.