高通量测序技术在探索芽菜微生物群落结构中的应用

　　摘要：文章主要通过高通量测序技术探索市场上不同品牌宜宾芽菜之间细菌菌群结构的差异性。采用Alpha多样性分析和Beta分析，从门水平上对细菌菌群结构分析，发现优势菌门是厚壁菌门(Firmicutes)，在3个样本中占比分别为43.77%、50.73%、45.59%;在种的水平上进行分析发现，优势菌种是鼠乳杆菌(Lactobacillusmurinus)，在3个样本中占比分别为10.51%、9.87%、11.92%。此外，还在芽菜中检测出脑膜炎嗜冷杆菌(Psychrobactermeningitidis)和荧光假单胞菌(Pseudomonasfluorescens)等腐败致病菌。结果表明，3种品牌宜宾芽菜之间细菌菌群结构有一定差异，并且市场上流通芽菜的微生物安全需要改善。

　　关键词：芽菜;高通量测序技术;群落结构;优势菌种

　　芽菜是在我国四川南部较为盛行的一道传统酱腌菜，由新鲜的芥菜茎部通过清洗、划丝、晾晒、腌制、发酵等工艺制作而成[1]。芽菜发酵是芽菜制作过程中重要的环节，因为芽菜的风味主要来源于发酵环境中微生物代谢活动[2-3]。芽菜自然发酵过程中会有微生物群落结构变化，这也是影响芽菜风味和食用安全性的关键点。但是，目前对于芽菜中的微生物菌群结构和多样性分析并未完全了解。在过去的研究中，主要有依赖微生物的培养特点和形态学的传统培养分离方法来探索发酵食品中菌群。

　　马长路等[4]通过传统培养分离的方法从东北酸菜中分离出植物乳杆菌(Lactobacillusplantarum);尹曦等[5]用传统培养分离方法发现宜宾芽菜中主要有戊糖片球菌(Pediococcuspentosaceus)、植物乳杆菌。但依靠形态学来鉴定微生物常常不够全面，因为能通过培养基进行培养的菌群很少，所以传统的方式不能全面地认识发酵环境中的菌群结构。随后基于 PCR 的检 测手段出现，主要有 PCR-DGGE 技术、EIRC-PCR 技术和 RFLP技术等。

　　张先琴[6]通过 PCR-DGGE 技术分析了芽菜中的微生物群落结构，但无法有一个准确的定量;荆雪 娇[7]通 过 PCR-DGGE 技术发现短乳杆菌为泡菜和酸菜的优势菌群，弯曲乳杆菌在泡菜和浙江榨菜中为优势菌群，消化乳杆菌和食用糖乳杆菌分别为酸菜和酱菜的优势菌群，但基于 PCR的技术手段往往因为要凝胶电泳，图谱复杂且存在工作量大、灵敏度有限等缺点。近几年，高通量测序技术的诞生为了解微生物菌群结构提供了新的解决方案，因其具有无需克隆、测序量 大、速 度 快 等 特 点[8-10]而 广 泛应 用 于 检测植物病毒、畜禽、食源性病原体[11]等。

　　杨柳等[12]利用高通量测序发现在辣白菜中 Weissella 菌属与菌落计数相距最近，同时也与总酸含量的呈味物质呈正相关，其贡献率最好，其他菌属呈负相关，说明 Weissella菌属对主要呈味物质有重要作用。贾 晶 晶 等[13]利 用高通量测序技术对12%食盐浓度发酵的腌渍白菜 液分析得出共有6个菌门。最优势菌门为变形菌门，其相对含量高达67.01%，丙型变形菌纲相对含量为99.17%，为其主要优势菌纲，其中普罗威登斯菌属为最优势菌属。由此可以看出高通量测序技术对微生物群落多样分析是较准确且完备的。本研究重点在于利用高通量测序技术从门、种水平上研究市场上不同品牌的宜宾芽菜之间微生物群落结构差异性及是否存在有害菌，在检测市场产品的安全性的同时也为厂家在发酵过程中规避一些不必要的有害微生物繁殖起到提醒作用。

　　1 材料与方法

　　1.1 验材宜宾芽菜：购于宜宾叙州区当地贸易市场，购买3种不同厂家的袋装芽菜，分别编号 C1、C2、C3。所有样品均为无菌采样，使用无菌塑封袋包装后置于-80 ℃冰箱中保存待测。

　　1.2 验试剂与仪器E.Z.N.ATM Mag-BindSoilDNA提取试剂盒 美国Omega公司;通用型恒压恒流电泳仪;AppliedBiosystemsⓇGeneAmpⓇ PCR System 9700 PCR 仪;QIAquick GelExtractionKit凝胶回收试剂盒 Qiagen公司;Qubit2.0微量荧光核酸定量仪 ThermoFisher公司;X1R高速冷冻离心机 美国 Thermo公司;HiSeq测序仪、HiSeqRapidSBSKitv2(FC-402-4023500Cycle)测序试剂盒 美国Illumina公司;其他所需试剂均购自成都市科隆化学品有限公司。

　　1.3 DNA 提取和 PCR 、扩增使用特定 的 DNA 提取 试剂 盒 进 行 基 因 组 DNA提取后，使用0.8%琼脂糖凝胶电泳检测 DNA。引物引 用：515F(5′-GTGCCAGCMGCCGCGGTAA-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′)，以稀释后的基因组 DNA为模板，根据测序区域的选择，使用带Barcode的特异引物进行PCR。每一个25μL的体系包括1×PCRbuffer、1.5mmol/LMgCl2、0.4μmol/LdNTPs、正向和反向 引 物 各1.0μmol/L、0.5 U KOD-Plus-Neo酶(Toyobo公司)和10ng模板。

　　PCR 程序包括起始94℃1min，然后30次循环(变性94℃20s，退火54℃30s和延伸72℃30s)，最后72℃5min。每个样本进行3个PCR技术重复。PCR与1/6体积的6loadingbuffer混合，使用2%琼脂糖凝胶电泳检测。取目的条带用来回收，回收使用 QIAquickGelExtractionKit(Qiagen公司)。使用 Qubit@2.0Fluorometer(ThermoScientific公司)定量。最后等摩尔量混合。建库使用 TruSeqDNAPCR-FreeSamplePrepKit，构建好的文库经过定量和文库检测合格后，使用 HiSeq2500平台 PE250模式测序。

　　1.4 据分首先使用FLASH[14]拼接双端序列，截去Barcode序列得到 原 始 数 据，然 后 使 用 Trimmomatic进 行 质控。在 Usearch软 件[15]基 础 上，使用 UPARSE 算 法在97%水平上进行 OTU 聚类，再用 R 语言进行各种数据转换，最 后 用 ggplot2包 作 图。对于 Alpha多样性分析，PD 指 数 的 计 算 使 用 Picante包，其 他 指数 使用 Vegan包。对于 Beta多样性分 析，使 用 GUniFrac包计算 UniFrac距离。

　　2 结果与分析

　　2.1 芽菜中细菌群落 Alpha多样性分析Alpha多样性指数主要有细菌菌群丰度指数 Chao1指数和 ACE指数，细菌菌群多样性指数 Simpson指数和Shannon指数。C3样本中细菌菌群丰度最高，而 C2样本中细菌菌群的丰度最低。从芽菜中细菌菌群多样性指数的结 果 发 现，C1 样本中的细菌菌群多样性最高，C2样本中的细菌菌群多样性最低。3种芽菜在细菌菌群多样性和丰度上具有一定的差异，这可能是由于原料的选择、制作工艺以及灭菌条件的不同造成的。

　　2.2 菌群落韦恩图分析将细菌序列在97%相似水平上定义为一个 OTU，对3种宜宾芽菜中细菌的 OTU 进行分析得到 OTU 韦恩图。3个样 品 共 测 得1948个 OTUs。C1中测得1351个 OTUs，C2中测得1315个 OTUs，C3中测得1444个OTUs，3种宜宾芽菜中共有727个OTUs，结果表明，3种芽菜细菌群落结构之间存在一定差异，但相比 3 种 芽 菜 共 有 菌 落，每 种 芽 菜 独 有 菌 落 占 比更小。

　　2.3 生物相对丰度曲线根据 OTU 丰度，使用Rank-Abundance曲线对样本多样性进行评估，得到相对丰度曲线。

　　曲线的平滑程度反映了样本中物种的均匀度，可以看出随着测序加深，序列增加，各个样品的曲线斜率逐渐减小，呈现一个趋于平缓的态势，说明物种分布较均匀。而曲线在水平横轴上的范围越大，说明物种丰度越高，可 知 C3样 品 的 物 种 丰 度 最 高，C1次之，C2最低。

　　2.4 细菌群落在门水平上的结构差异分析在门的水平上对3种宜宾芽菜中细菌菌群结构进行分析。

　　3种 芽菜中优势菌门均是厚壁菌门(Firmicutes)，C1样本中厚壁菌门占比为43.77%，在C2样本中占比为50.73%，在C3样本中占比为45.59%。厚壁菌门会参与膳食植物多糖代谢[16];其次是拟杆菌门(Bacteroidetes)，拟杆菌门在 C1、C2和 C3样本中占比分别是22.53%、21.98%和23.77%。Ley等[17]发现厚壁菌门与拟杆菌门的比例降低对减轻体重有直接的影响，在肥胖受试者中通常厚壁菌门含量较高，拟杆菌门含量较低，因此在样本中若拟杆菌门比例较高，可能起到抑制人体脂肪沉积的积极作用。在样本中还检测出变形菌门(Proteobacteria)，少数的奇古菌门(Thaumarchaeota)、绿弯菌 门(Chloroflexi)、广 古 菌门(Euryarchaeota)、放线菌门 (Actinobacteria)、醋 杆 菌 门 (Acidobacteria)、螺旋体门(Spirochaetes)等。其中，变形菌门属于革兰氏阴性细菌，它们具有分泌脂多糖的能力，而脂多糖可以作为内毒 素，容 易 引 起 肠 道 炎 症[18]。对 此，我 们日常饮食可以注意多食用高纤维或碳水化合物改善肠道菌群来降低炎症发生的概率。

　　2.5 细菌群落在种水平上的结构差异分析在种的水平上对3种宜宾芽菜中细菌结构进行分析。

　　在3组样品中丰度前三的优势菌种一致，分别 为 鼠 乳 杆 菌(Lactobacillusmurinus)、禽 乳 杆菌(Lactobacillusaviariu)、格氏乳杆菌(Lactobacillusgasseri)，共属于乳酸杆菌属(Lactobacillusspp.)，乳 酸 杆菌 具有嗜酸的特点，因此广泛分布于传统酱腌菜中。

　　优势菌种鼠乳杆菌在3组样品中占比分别为10.51%、9.87%、11.92%。鼠乳杆菌在肠道中起着调节肠道菌群，增强肠道屏障的功能[19-20]，Huang等[21]研究发现，鼠乳杆菌可以作为一种功能性的益生菌来对抗食物过敏。禽乳杆菌属于革兰氏阳性菌，具有产生细菌素和有机酸的能力，可以抑制病原体和炎症反应[22]，它在3种芽菜样品中占比分别为4.70%、9.82%、5.05%，其中禽乳杆菌在C2中占比最大，而在 C1中占比最小。

　　格氏乳杆菌在3组样品中占比分别为2.56%、2.94%、2.42%，Kyoung等[23]研究发现格氏乳杆菌会影响体内代谢因子，与Cudraniatricuspidata 协 同 作用 调 节 代 谢 功 能，发 现 了 一 种 新的作为肥胖与消瘦的生物标志物。此外，还发现了弯曲乳杆菌(Lactobacilluscurvatus)，有研究表明弯曲乳杆菌具有较高的稳定性，对单核增生细菌有显著的抑制作用，从而起到改善发酵香肠安全的作用[24-25]。

　　同时，3个样品中均含有罗伊乳杆菌(Lactobacillusreuteri)，它常用于婴幼儿保健食品益生菌使用，可以促进肠道蠕动，增加人体免疫力[26]。特别注意的是有腐败致病菌被检测出，其中一个是脑膜炎嗜冷杆菌(Psychrobactermeningitidis)，它可以产生褐藻酸降解菌快速吸收水分从而导致紫菜出现病斑进而腐败[27]，因此在芽菜中可能存在加速腐败的效 果;此外还检测出荧光假单胞菌 (Pseudomonasfluorescens)，有研究表明荧 光 假 单 胞 菌 常 附 着 于 酱 腌菜表面，存在还原硝酸盐生成亚硝酸盐的情况，因此可能对食品安全有不良影响[28]。

　　另外，我们还检测到了产酸拟杆菌(Bacteroidesacidifaciens)，MarquesFZ等[29]发现Bacteroidesacidifaciens等肠道菌群数量会增加高血压小鼠 发 生 心 力 衰 竭 概 率，KangChil-sung等[30]发 现Akkermansiamuciniphila、Bacteroidesacidifaciens的减少与 DSS诱导的结肠炎相关。

　　总体上，在种水平上的有益功能菌主要有鼠乳杆菌(Lactobacillusmurinus)、禽 乳 杆 菌 (Lactobacillusaviariu)、格氏乳杆菌(Lactobacillusgasseri)、罗伊乳杆菌(Lactobacillusreuteri)，腐败致病菌主要有脑膜炎嗜冷 杆 菌(Psychrobactermeningitidis)、荧 光 假单 胞 菌(Pseudomonasfluorescens)、产 酸 拟 杆 菌 (Bacteroidesacidifaciens)。在传统工艺下生产的芽菜，发酵菌群主要来源于原辅料和生产环境中，如何减少腐败致病菌和扩大优势菌种是提高芽菜食用安全性和风味的关键点。

　　3 结论

　　本文采用高通量测序技术对3种不同宜宾芽菜中的细菌菌群 结 构 进 行 分 析，结 果 表 明，3 种 宜 宾 芽 菜中，C3样本中细菌菌群丰度最高，而在 C1样本中，细菌菌群多样性最高。此外，共检测到3种宜宾芽菜中共有727个 OTUs，3种芽菜具有一定相似性。在门的水平上，对3种宜宾芽菜中细菌菌群结构分析，发现厚壁菌门是优势菌门，在C1、C2和C3中占比分别为43.77%、50.73%和45.59%。

　　在种的水平上对3种宜宾芽菜中细菌菌群结构进行分析，发现鼠乳杆菌是3种宜宾芽菜中的优势菌种，它在C1、C2、C3中占比分别为10.51%、9.87%、11.92%。其次分别是 禽 乳 杆 菌(Lactobacillusaviariu)、格 氏 乳杆 菌(Lactobacillusgasseri)。禽乳 杆菌 在 C1、C2、C3中 占 比分别为 4.70%、9.82%、5.05%，格 氏 乳 杆 菌 在 C1、C2、C3中占比分别为2.56%、2.94%、2.42%。

　　通过高通量测序技 术，能快速、准确地分析芽菜发酵环境中的细菌菌群结构，本研究发现市面售卖的宜宾芽菜中主要功能菌是鼠乳杆菌(Lactobacillusmurinus)、禽乳杆 菌(Lactobacillusaviariu)、格 氏 乳杆 菌(Lactobacillusgasseri)、罗伊乳杆菌(Lactobacillusreuteri)，但也存在部分腐败致病菌，如脑膜炎嗜冷杆菌(Psychrobactermeningitidis)、荧光 假 单 胞 菌 (Pseudomonasfluorescens)、产 酸 拟 杆 菌(Bacteroidesacidifaciens)。本研究为芽菜的生产、包装、运输等过程中的微生物控制提供了参考，同时也为后期对芽菜的深入研究提供了理论基础。

　　参考文献 ：

　　[1]刘达玉，吴枚枚，肖龙泉，等.腌制芽菜加工现状与发展技术路线[J].中国调味品，2019，44(9)：180-182.

　　[2]左勇，鞠帅，谢晖，等.宜宾芽菜生产技术展望[J].中 国 调 味品，2011，36(5)：12-14.

　　[3]杜玫，蔡玥，贾利蓉，等.不同前处理方式对泡青菜挥发性风味成分的影响[J].中国调味品，2021，46(4)：38-42，49.

　　[4]马长路，焦梦丽，罗红霞，等.传统东北酸菜中乳杆菌的分离鉴定及其降胆固醇特性研究[J].食品研究与开发，2019，40(21)：209-213.

　　[5]尹曦，蒲彪，陈安均，等.四川宜宾芽菜中乳酸菌分离筛选及菌相分析[J].食品科学，2013，34(3)：163-168.

　　[6]张先琴.PCR-DGGE技术分析传统发酵蔬菜中微生物群落结构[D].成都：四川农业大学，2012

　　作者：王新惠