鲊广椒真菌中滋味品质检测方案(感官智能分析)

食品科学※生物工程2018，Vol.39，No.18173 鲊广椒真菌多样性及其对滋味品质影响的评价 王玉荣，代凯文文，沈馨，董童蕴，周书楠'，郭 壮* (1.湖北文理学院化学工程与食品科学学院，鄂西北传统发酵食品研究所，湖北襄阳 441053：2.当阳市食品药品监督管理局，湖北当阳 444100) 摘 要： 采用MiSeq高通量测序技术对鲊广椒真菌微生物的群落结构进行解析的基础上，使用电子舌对不同样品各滋味指标的相对强度进行测定，最终对真菌多样性与鲊广椒滋味品质进行关联性分析。结果发现，在门水平上，Ascomycota(子囊菌门)的平均相对含量高达99.97%。在属水平上，相对含量大于1.0%的真菌属分别为Candida(念珠菌菌)、Eurotium(曲霉菌属)、Pichia(毕赤酵母属)、Fusarium(镰刀霉属)、Galactomyces(地霉属)、Cladosporium (分枝子菌属)、Debaryomyces(德巴利氏酵母属)和Guehomyces(久浩酵母属)，其平均相对含量分别为56.63%、10.12%、11.19%、2.67%、2.48%、2.21%、1.50%和1.01%。尔皮尔森相关性分析发现优势真菌属与鲊广椒的苦味、涩味、后味A和后味B呈正相关。由此可见，真菌含量过高不利于鲊广椒滋味品质的形成。 关键词：鲊广椒；高通量测序；电子舌；多样性：品质 Characterization of Fungal Microflora and Its Influence on Taste Quality of Zhaguangjiao， a Chinese TraditionalFermented Chili Product WANG Yurong'， DAI Kaiwen’， SHEN Xin’， DONG Yun’， ZHOU Shunan， GUO Zhuang* (1. Northwest Hubei Research Institute of Traditionally Fermented Food， School of Chemical Engineering and Food Science，Hubei University of Arts and Science， Xiangyang 441053， China；2. Dangyang Municipal Food and Drug Administration， Dangyang 444100， China) Abstract： The fungal community structure in Zhaguangjiao was analyzed by MiSeq high throughput sequencing and therelative intensity of various taste attributes of Zhaguangjiao samples was determined using an electronic tongue. Furthermore，the correlation between the fungal diversity and taste quality of was evaluated. The results showed that Ascomycota was themost dominant fungal phyla with a relative abundance of 99.97%. Candida， Eurotium， Pichia， Fusarium， Galactomyces，Cladosporium， Guehomyces， and Debaryomyces each accounted for more than 1.0% of total sequences at the genus level，with a relative abundance of 56.63%，10.12%，11.19%，2.67%，2.48%，2.21%， 1.50% and 1.01%， respectively. The resultsof Pearson correlation analysis indicated that there were positive correlations among the relative abundance of the dominantfungal genera and the defective indexes of taste quality. Thus， we can consider that the fungal microbiota played inactiveroles in the taste quality formation of Zhaguangjiao. Keywords： Zhaguangjiao；high throughput sequencing； electronic tongue； diversity； quality DOI：10.7506/spkx1002-6630-201818027 中图分类号：TS201.4 文献标志码：A 文章编号：1002-6630(2018)18-0173-06 引文格式： 王玉荣，代凯文，沈馨，等.鲊广椒真菌多样性及其对滋味品质影响的评价[J].食品科学，2018，39(18)：173-178.DOI：10.7506/spkx1002-6630-201818027.http：//www.spkx.net.cnWANG Yurong， DAI Kaiwen， SHEN Xin， et al. Characterization of fungal microflora and its influence on taste qualityof Zhaguangjiao， a Chinese traditional fermented chili product[J]. Food Science， 2018，39(18)： 173-178. (in Chinese withEnglish abstract) DOI：10.7506/spkx1002-6630-201818027. http：//www.spkx.net.cn ( 收稿日期：20 1 7-08- 2 0 ) ( 基金项目 ：湖 北文 理 学院教师科研能力培育基金项目 (2017k ypy051) ) ( 第 一 作 者 简 介：王玉 荣 (1993一 )，女 ，硕士 ， 研 究 方向为食 品 生物技术。E- m ail： w a n g y urong1993@163. co m ) ( * 通 信作者简介：郭 壮 (198 4 一) ， 男 ，副教授，博士，研究方向为食 品 生物技术。E-mail： guoz huang1984@163 . com ) 作为我国特色传统发酵食品，湖北宜昌地区制作的鲊广椒以大米面、鲜红辣椒和食盐为主要原料，室温下发酵15~20d而成，具有酸辣适口、滋味鲜香的特点。鲊广椒的发酵通常在倒扣于水盆的坛子中完成，同时坛口会塞入半坛稻草以防止原料浸入水中""。虽然发酵过程均在密封环境中完成，但发酵初始阶段坛中存在的一定量空气，为专性好氧真菌的生长提供了必要的氧气。此外，鲊广椒的制作环境相对开放，加之环境中存在少量兼性厌氧和微好氧真菌，因而发酵好的鲊广椒产品中存在一定量的真菌，而该部分真菌的存在可能会对产品品质及使用安全性产生一定的影响，所以对鲊广椒中真菌多样性进行评价是极为必要的。 以Illumina MiSeq测序平台为代表的第2代测序技术具有通量高和检测速度快的特点，通过为每个样品在测序前加上一段被称作样品特异性条形码或者标签的核苷酸序列，实现了多样本的平行实验，解决了传统指纹图谱技术不能进行关联分析的缺点。近年来，该技术被广泛应用于酸奶、腊肠、葡萄酒、窖泥和泡菜等发酵食品微生物多样性研究和特异菌种的检测中，为全面和深入研究不同发酵食品基质中微生物分子生态学提供了新的技术手段。电子舌实现了食品酸味、苦味、涩味、鲜味、咸味、甜味、苦味回味、涩味回味和鲜味回味的数字化评价，具有结果稳定且受外界影响小的优点，目前已经广泛应用于米酒1121、啤酒131、水产品4和醋等发酵食品的滋味品质评价中。 本研究从湖北省宜昌市当阳地区农户家中采集了10个鲊广椒样品，首先采用MiSeq高通量测序技术对其真菌微生物多样性进行评价，继而使用电子舌系统对各滋味指标相对强度进行测定，最终对真菌多样性与鲊广椒滋味品质进行了关联性分析，以期对后续鲊广椒产品品质改良及安全性评价提供理论支持。 1 材料与方法 1.1 材料与试剂 鲊广椒采集自湖北省当阳市草埠湖镇和玉泉办事处10个农户家。 QIAGEN DNeasy mericon Food Kit食品DNA基因组提取试剂盒 德国QIAGEN公公；5XTransStartFastPfu Buffer、FastPfu Fly DNA Polymerase、dNTPs Mix北京全式金金物技术有限公司；DNA 1000试剂盒美国Agilent公司；味觉标准溶液、阴离子和阳离子溶液、内部溶液及参比溶液 日本Insent公司。 1.2 仪器与设备 MiSeq高通量测序平台 美国Illumina公司； R920机架式服务器 美国Dell公司； SA 402B电子舌(配备AAE、CT0、CA0、AE1、C00和GL1测试传感器) 日本Insent公司；2100芯片生物分析仪 美国Agilent公司；5810R台式高速冷冻离心机 德国Eppendorf公司；DYY-12电泳仪 北京六一仪器厂；ND-2000C微量紫外分光光度计 美国Nano Drop公司：UVPCDS8000凝胶成像分析系统 美国Bio-Rad公司；vetiri梯度基因扩增仪 美国AB公司。 1.3 方法 1.3.1 样品微生物宏基因组DNA提取 取1.0g鲊广椒， 参照QIAGEN DNeasy mericon FoodKit试剂盒说明书进行DNA提取，检测合格后的DNA样品置-20℃暂存备用。 1.3.2 真菌18S rRNA V~V，区聚合酶链式反应 (polymerasechain reaction， PCR)扩增 正向引物为SSU0817F (5'-TTAGCATGGAATAATRRAATAGGA-3')，反向引物为SSU1196R(5'-TCTGGACCTGGTGAGTTTCC-3')116，其中在正向引物中加入7个核苷酸标签(barcode)。PCR扩增体系为：4pL5×PCR缓冲液， 2 pL2.5 mmol/L dNTPs mix，0.8 pL 5 umol/L正向引物， 0.8 pL 5 pmol/L反向引物，0.4 pL 5 U/uL DNA聚合酶， 10 ng DNA模板， 体系用ddH，O补充至20 uL。PCR扩增条件为：95℃3 min；95 ℃ 30 s， 55 ℃ 30s， 72 ℃ 45 s， 30个循环；72℃10 min。 1.3.3 样品平衡及MiSeq高通量测序 将10个经1.0%琼唐糖凝胶电泳检验合格的扩增产物按照100 nmol/L浓度进行稀匀混匀后，送至上海美吉生物医药科技有限公司使用MiSeq高通量测序平台进行测序。 1.3.4 序列质控 根据成对序列之间的重叠(overlap) 关系，将下机的双端序列数据拼接(merge)成一条序列。在序列拼接过程中根据以下标准进行质控： overlap区的碱基数需大于10 bp；最大错配比率需小于等于0.2；拼接好的序列需无barcode碱基错配；引物碱基错配数需小于等于2 bp。 序列拼接后，依照barcode信息将下机序列划分到各样品，对序列方向进行校正后切掉barcode和引物，进而得到高质量的序列。若切掉barcode和引物后的序列其碱基数小于50bp， 则该序列亦予以切除。 1.3.5 生物信息学分析 质控合格后的序列，采用QIME(v1.70)平台进行真菌微生物物种分析和多样性评价4。主要的处理流程为：1)采用PyNAST校准并把序列排齐：2)首先在100%相似性下进行UCLUST归并，建立无重复的单一的18S rRNAV~V，区序列集，继而在97.0%相似性下进行序列归并，构建分类操作单元(operational taxonomicunits， OTU) 矩阵；3)应用ChimeraSlayer去除含有嵌合体序列的OTU序列201；4)使用SILVA数据库进行序列 同源性比对，在门、纲、目、科和属水平上对其分类学地位进行明确；5)使用FastTree软件，绘制基于OTU代表序列的系统发育进化树22；6)采用发现物种数和香农指数(Shannon-Wiener index) 等a多样性指标分别对单个样品真菌微生物的丰富度和多样性进行评价，同时，采用稀疏曲线和香农指数曲线对测序深度是否满足后续分子生物学分析进行评估231；7)分别基于UniFrac距离P241进行主坐标分析和非加斤组平均(unweighted pair-groupmethod with arithmetic means， UPGMA)法聚类分析， 进而完成不同样品间真菌微生物的β多样性分析。 1.3.6 核酸登录号 本研究中所有序列数据已提交至MG-RAST数据库，登录号为mgp81949。 1.3.7 基于电子舌技术鲊广椒滋味品质的评价 称取10g鲊广椒样品于190mL蒸馏水中，3000 r/min离心10 min后取上清液过滤，按照下述步骤对鲊广椒样品酸、苦、涩、咸、鲜等基本味及涩、苦和鲜味回味进行测定251：1)将C00和AE1传感器浸泡在阳离子溶液中，CA0、CT0和AAE传感器浸泡在阴离子溶液中，5个传感器浸泡时间均为90 s， 以去除传感器上的吸附物；2)5个传感器在参比溶液1和2中分别洗涤120s后， 在参比溶液3中浸泡30 s， 测得参比溶液的电势值V：3)5个传感器在某鲊广椒样品中浸泡30s， 测得电势值V，通过计算V一V值即可得到样品酸、苦、涩、咸和鲜味5个基本味的相对强度值；4)传感器C00、AE1和AAE于参比溶液4和5中分别洗涤3s后，于参比溶液6中浸泡30 s， 测得电势V'， 通过计算V'一V，值即可得到样品后味A(涩的回味)、后味B(苦的回味)和丰度(鲜的回味)的相对强度值。每个样品重复测定4次，选后3次纳入数据分析。 GL1传感器在阴离子溶液中浸泡90s后，参照上述程序单独对鲊广椒样品的甜味进行测定。每个样品重复测定5次，选后3次纳入数据分析。 1.3.8 多元统计学分析 使用皮尔森相关性分析法对平均相对含量大于1.0%的真菌属之间的相关性进行计算，并采用热图对结果进行展示；亦使用皮尔森相关性分析法对鲊广椒中优势真菌属和电子舌各滋味指标相对强度之间的相关性进行计算，同时选取相关系数绝对值大于0.5的指标，采用Cytoscape软件(v3.5.1)进行相关性网络图绘制。使用Origin 8.5软件 (OriginLab Corp， MA， USA) 和R软件(v3.3.2)进行作图。 2 结果与分析 2.1 序列丰富度和多样性分析 在提取鲊广椒样品宏基因组DNA的基础上，本研究采用MiSeq高通量测序技术，对鲊广椒样品中真菌微生物多样性进行了评价。样品18SrRNA测序情况及各分类地位数量如表1所示。 表1 样品18SrRNA测序情况及各分类地位数量 Table 1 18S rRNA read counts and number of identifiable units atdifferent taxonomical levels 样品 序列 OTU 门/个纲/个 目/个科/个属/个 发现物香农 编号 数/条 数/个 种数 指数 DY1130970 913 4 912 4.89 DY2 42007 603 4 475 4.15 DY3 43 032 857 10 31 702 3.66 DY4436020 681 4 11 21 615 1.13 DY5 443 514 664 5 12 21 501 2.01 DY638271 559 5 11 470 1.26 DY7 37 697 440 7 364 0.31 DY8413911 1280 5 14 43 1 072 5.07 DY99331206 318 2 10 11 12 314 1.16 DY10 32 148 44D 4 6 428 1.75 注：计算每个样品发现物种数和香农指数时，样品的测序量均为30810条序列。 由表1可知，10个鲊广椒样品共产生了376256条高质量的18SrRNA序列，平均每个样品产生37626条。本研究采用两步UCLUST法进行OTU的划分，根据序列100%相似性聚类分析后，得到69496条代表性序列，继而经过序列97.0%相似性聚类分析后，得到5457个OTU， 经嵌合体检查后去除14个OTU， 还剩余5443个OTU进行后续分析。本研究有18.45%的OTU和3.91%的序列(范围为0%~14.11%，标准差为5.15%)不能鉴定到属水平。 采用稀疏曲线和香农指数曲线，本研究进一步对现有测序深度是否满足后续分子生物学分析进行了评估，其结果如图1所示。 图1 稀疏曲线图(A)和香农指数图(B) Fig.1 Rarefaction analysis (A) and Shannon diversity estimates (B) 由图1可知，随着测序深度的增加，鲊广椒样品中发现真菌物种数的数量亦随之增大且稀疏曲线尚未进入平台期，而香农指数曲线已完全达到平衡状态。由此可见，随着测序深度的增加虽然有新的真菌种系型可能被发现，但本研究已经完全捕获到了鲊广椒样品中真菌类群的多样性信息，因而本研究的测序深度可以满足后续生物信息学分析要求201。 2.2 基于各分类学地位相对含量的鲊广椒真菌微生物构成分析 在门水平上，鲊广椒样品中微生物主要隶属于7个真菌门， 其中Ascomycota(子囊菌门)的平均相对含量高达99.97%。在属水平上，共发现135个真菌属，鲊广椒中优势真菌属相对含含的比较分析如图2所示。 图2 鲊广椒中优势真菌属相对含量的比较分析 Fig.2 Relative abundances of the major fungal genera inZhaguangjiao samples 由图2可知，鲊广椒中隶属于Ascomycota(子囊菌门)且平均相对含量大于1.0%的真菌属及其相对含量分别为：Candida(念珠菌属， 56.63%)、Eurotium (曲霉菌属，10.12%)、Pichia(毕赤酵母属， 11.19%)、Fusarium(镰刀霉属，2.67%)、Galactomyces(地霉属，2.48%)、Cladosporium(分枝孢子菌属，2.21%)和Debaryomyces(德巴利氏酵母属，，1.01%)。此外，还有隶属于Basidiomycota(担子菌门)的Guehomyces(久浩酵母属)，其平均相对含量为1.50%。由此可见，当阳地区鲊广椒主要是由若干个隶属于Ascomycota(子囊菌门)已知的优势真菌属组成，其平均累计含量为86.32%。Candida(念珠菌属)和Pichia (毕赤酵母属)在10个样品中均存在，其平均累计含量分别为67.82%。平均相对含量大于1.0%的真菌属相关性的热图如图3所示。 由图3可知， Galactomyces(地霉属)与Candida(念珠菌属)呈显著负相关(R=-0.673，P<0.05)，与Eurotium (曲霉菌属)呈极显著正相关(R=0.893，P<0.01)， Guehomyces(久浩酵母属)与Fusarium(镰刀霉属)呈极显著正相关(R=0.867， P<0.01)，Debaryomyces (德巴利氏酵母属)与Fusarium(镰刀霉属， R=0.667)和Guehomyces(久浩酵母属，R=0.719) 均呈显著正相关(P<0.05)。本研究进一步在OTU水平上对鲊广椒真菌微生物的多样性进行分析， OTU在10个样品中出现次数统计如图4所示。 *.显著相关(P<0.05)；**.极显著相关(P<0.01). 图3 平均相对含量大于1.0%的真菌属相关性的热图 Fig. 3 Heat map of correlation among the fungal genera with relativeabundances more than 1.0% 图4 OTU在10个样品中出现次数统计 Fig.4 Distribution of OTU as a function of their prevalence in the10 samples 由图4可知，纳入本研究的样品共产生5443个OTU， 然而在10个样品中仅出现1次的OTU有4715个，占OTU总数的86.63%，所包含序列数为13070条，仅占所有质控后合格序列数的3.47%。10个样品没有共有的OTU， 仅有5个OTU存在9个样品中，其包含序列数为203208条，高达所有质控后合格序列数的54.01%。由此可见，在OTU水平上，虽然若干鲊广椒样品存在大量的核心真菌菌群，但某些样品间微生物群落结构的构成存在较大差异。 2.3 基于多元统计学分析的鲊广椒真菌微生物构成分析 通过上述分析，本研究发现虽然多数鲊广椒样品存在大量的核心真菌菌群，但某些样品间微生物群落结构的构成存在较大差异，因而本研究进一步采用基于OTU水平的加权UniFrac距离主坐标分析和UPGMA聚类分析，对10个鲊广椒样品的β多样性进行了分析，基于分类操作单元的加权UniFrac离离的主坐标分析如图5所示。 图5 基于分类操作单元的加权UniFrac距离的主坐标分析 Fig.5 Principal coordinate analysis of OTUs based on weightedUniFrac distance 由图5可知，10个鲊广椒样品呈现明显的分离趋势，其中DY1、DY2、DY5和DY8， DY4、DY6、DY7和DY9， DY3和DY10分别呈现出聚类趋势。基于分类操作单元的加权UniFrac距离的UPGMA聚类分析如图6所示。 图6 基于分类操作单元的加权UniFrac距离的UPGMA聚类分析 Fig.6 UPGMA clustering analysis of OTUs based on weightedUniFrac distance 由图6可知，10个样品整体上可以分为3个聚类，其中聚类I由样品DY4、DY6、DY7和DY9构成，聚类II由样品DY1、DY2、DY5和DY8构成，聚类Ⅲ由DY3和DY10构成，这与基于分类操作单元的加权UniFrac距离的主坐标分析结果一致，亦进一步证实了虽然若干鲊广椒样品存在大量的核心真菌菌群，但某些样品间其微生物群落结构构成存在较大的差异。通过Krushkal-Wallis检验发现，隶属于不同聚类的样品除Eurotium(曲霉菌属)和Guehomyce(久浩酵母属)外，其他优势真菌属差异均不显著(P>0.05)。 2.4 鲊广椒优势真菌属及滋味品质的关联性分析 Fig.7 Box plot of relative intensity of each taste index (n=10) 在采用MiSeq高通量测序技术对真菌微生物多样性进行揭示的基础上，使用电子舌对鲊广椒样品的滋味品质进行评价，在此基础上构建了优势真菌属与产品滋味品质相关性的网络图。鲊广椒各滋味指标相对强度值的箱形图如图7所示。 由图7可知，鲊广椒样品在酸味上的差异最大，其次为咸味、甜味和鲜味，而在苦味、涩味、后味A(涩的回味)、后味B(苦的回味)和丰度(鲜的回味)指标上的差异较小。鲊广椒优势真菌属和滋味物质相关性的网络图如图8所示。 实线表示正相关，虚线表示负相关。 图8 鲊广椒优势真菌属和滋味物质相关性的网络图 Fig. 8 Correlation network diagram of major fungal genera and tasteattributes in Zhaguangjiao samples 由图8可知， Pichia (毕赤酵母属)与鲊广椒苦味呈正相关； Debaryomyces(德巴利氏酵母属)和Guehomyce(久浩酵母属)与鲊广椒涩味呈正相关； Debaryomyces(德巴利氏酵母属)、Galactomyces(地霉属) 和Eurotium(曲霉菌属)与鲊广椒后味A(涩味的回味)呈正相关； Galactomyces (地霉属)与后味B(苦味的回味)呈正相关。苦味、涩味、后味A和后味B均为鲊广椒的缺陷型指标，由此可见，优势真菌属多与鲊广椒品质的缺陷性指标相关，真菌含量过高可能不利于产品滋味品质的形成。 3 讨寸1论 近年来出现的以Illumina MiSeq为代表的第2代高通量测序技术可以全面、客观、无偏地了解某一微生态系统中微生物群落结构，具有通量高、检测速度快和检测结果准确的特点271。本研究采用该技术对当阳地区鲊广椒中真菌微生物多样性进行了评价，结果发现念珠菌属、曲霉菌属和毕赤酵母属为其优势真菌属，相对含量占到了真菌数的近八成。念珠菌是真菌中最常见的条件致病菌，曲霉菌属亦是引起多种物质霉腐的主要微生物之一，且鲊广椒样品中亦存在隶属于镰刀菌属、地霉属和分枝孢子菌属的腐败真菌281。制作鲊广椒的发酵方式 为厌氧发酵，然而发酵初始阶段坛中存在的空气及发酵后期居民常开坛取食鲊广椒，均为真菌的生长提供了必需的氧气，使鲊广椒受腐败真菌污染的机率大大增加。此外，制作鲊广椒的坛子常倒扣于水盆中并置于墙角或房檐下，用于密封的水亦常因家禽饮用或雨水溅入而发臭，这也增加了鲊广椒被腐败真菌污染的机会。由此可见，改善鲊广椒加工环境进而抑制腐败霉菌的生长，对保障鲊广椒的食用安全性具有积极的意义。 在采用MiSeq高通量测序技术对真菌微生物多样性进行评价的基础上，本研究使用电子舌系统对各滋味指标相对强度进行测定，结果发现鲊广椒样品在酸味上的差异最大，这可能与样品中乳酸菌含量不同有关1291，因为除真菌外鲊广椒中亦存在大量细菌。此外，通过将真菌多样性与鲊广椒滋味品质进行关联性分析发现，优势慎菌属多与鲊广椒品质的缺陷性指标相关，因而综上所述，当阳地区鲊广椒中真菌微生物主要隶属于念珠菌属、曲霉菌属和毕赤酵母属，且其含量过高可能不利于产品滋味品质的形成。此外， MiSeq高通量测序技术存在读长短的缺陷，因而二代测序只能在属的水平上研究样本中菌群的组成和变化，难以精确到种水平上301，因而在后续研究中积极引入新的测序技术并实现鲊广椒样品中活的真菌的绝对定量，对鲊广椒产品品质改良及安全性评价将具有积极的意义。 ( 参考文献： ) ( [1 ] 王巧 碧， 王丹 ，赵欠， 等 . SD E 和SP ME法对鲊海椒发酵 中 香气组 成的比较分析[J].食 品 科 学 ，2016，37 ( 4 ) ：10 8 - 1 1 4 . DO I ：10. 7506 / sp k x10 0 2-6 6 3 0 -2 0160 402 0 . ) ( [2 ] 葛平珍， 王 丹，周才琼 . 不同淀粉源对鲊海椒发酵过程中功能成分 的影响[J]. 食 品 科学，2015，3 6 (21)：1 91 -19 5 .DOI ： 1 0. 75 0 6 / spkx10 0 2- 6 6 30 -2 0 152 1036 . ) ( [3 ] C AP ORAS O J G， L AU B E R C L ， W A L T ER S W A， et a l . U l t r a -h i gh- th ro u g h p u t m icr o bia l c o m mun it y a n alysis on t h e Il lu mi na H i S eq a nd M iS e q platf orms[ J ] . The I S ME Jo u r n a l， 2 0 1 2 ， 6( 8 )：1 6 2 1 - 1 62 4 .DOI：10 . 103 8/i s m e j.2012.8. ) ( [4 ] ROH S W， K I M K H ， NA M Y D， e t a l . In v es ti ga t io n of archaeala n d bact e rial d iver si t y i n f e r me n t e d s ea fo od u sing b arco d ed p yr os equen c i n g[ J] . Th e ISME J o ur nal， 2 010 ， 4 ( 1) ： 1-1 6. DOI： 10. 10 3 8/i s me j.2 00 9.83. ) ( [5 ] KING J L ， L A RUE B L， N O V R OSKI N M ， et al. H igh -quali t ya n d h i g h-t h rou ghp u t m assiv e ly p ar all el sequencing of th e hu ma n mito ch o ndr ial g enom e u si n g t h e I llum i na M iSeq [ ]. F o r e n s ic S c i e n ce I n t e rn at i onal： Genetics ， 2 0 14 ， 12 ( 10) ： 128 - 1 35. DO I ： 10 . 1 0 1 6 / j . f si gen.20 1 4. 06 . 0 0 1 . ) ( [ 6 ] 智楠楠. I l lu m ina Miseq平台深度测定酸奶中微生物多样性[J] . 食品工 业 科技，201 6 ，3 7 (24 ) ： 78- 8 2. D O 1 ： 1 0 . 13386/ i . issn1 0 0 2- 0 30 6 .2 016 .24 . 0 0 7. ) ( [7 ] J U S T Y NA P， AN N A LI S A R， V INCEN Z A P ， e t a l . Bac t e ria l di ve rsi t y in t ypi c al I t a l ian sala m i a t d iff e r ent r i pe ni n g s t ages as r e v e a l e d b y h ig h-t hr oughp u t o f 1 6S r D NA a m p l i c ons[J] . F o o d M i c rob io lo g y ，2 0 1 5 ，46 ( 4) ： 342-3 56 .D OI：1 0 .10 1 6/ j .fm. 2 0 1 4 . 0 8.0 23 . ) ( [ 8] BOYNT O N P J ， GR E I G D. F ungal d ive r sity a nd e cosys t e m f un ctionda t a f ro m w in e f e r m e n ta t ion vats an d mi c r o c o s m s[J ]. D a ta i n Brie f，2 016， 8(12)：225 -2 2 9. DO I：10 . 1 0 16/ j. d ib.2 016 .05 . 03 8. ) ( [ 9] HU X ， D U H ， R E N C . et al， Il l um i na t ing anaero b ic m i c r ob i alc om munity an d c o o cc ur re n c e p at t e r n s ac r o ss a q u a lity gradien t in C h inese l i qu or f erm en ta t i o n p it mu ds[ J ]. Ap pl i e d an d E nvi r onment a lM ic r o b i o l o gy ，2 016 ，8 2 (8)：25 06 -25 1 5.D OI：1 0.1 12 8 / A EM . 034 0 9 - 15. ) ( [10 ] YA N G H ， WU H. G AO L ， et a l . E f f e c t s of L a c tobacillus curvatusan d L e uco nos to c m e se nt er oid e s on Sua n C a i f e rm e nt at i o n in n orthea st ) ( Ch ina[J ]. J o urna l o f Mic r o bi ol ogy an d B i ot e c hnol o gy ， 2 01 6 ， 2 6( 12 )： 2 148- 2 158 .D OI ： 1 0 .4 01 4 / jm b. 1607 . 070 1 0. ) ( [11 ] KHAN M R R， K H ALIL I AN A ， K A N G S W. A hi gh s ensiti v ity I D C - elec t r o n ic t o n gu e u s ing d i e l e ct r ic/sensi n g m e m bra nes w i th so lv atoc hr omic d ye s[ J]. Se ns o rs， 2 0 16 ， 16( 5) ： 66 8- 6 8 7. D O I： 1 0 . 3 390 / s16 050668 . ) ( [12] 王玉荣，张俊英，胡欣洁，等.湖北孝感和四川成都地区来源的酒曲 对米酒滋味品质影响的评价[].食品科 学， 2015，36( 16 ) ： 2 0 7- 21 0 . D OI：1 0. 7 50 6/ s p kx1 002 -6630 - 2 015 1 603 8 . ) ( [ 1 3] GU T I ERREZ J M， HA D DI Z， A M A RI A . et a l . H y b r id el ectr onic tongue ba s ed o n m u lti sen so r dat a f u si o n f o r discri m i n a t ion o f b e e r s[J ] . S e ns o r s an d A c t uat o r s B： Che m ica l， 2 0 1 3 ， 1 7 7(2 ) ：98 9- 996 .D O I：10 . 1 0 1 6 /j . s n b.2 012 .1 1. 1 10 . ) ( [ 1 4] 张晶晶，顾赛麒，丁玉庭，等.电子舌在中华绒螯蟹产地鉴别及等 级评定的应用[J].食 品科 学 ，2015，3 6 (4)；141-1 46. DOI： 1 0.7506 / s p kx1 002 -6 63 0-20 15 0 4027 . ) ( [1 5] L I U M ， W AN G M， W ANG J ， et a l . C o m p a ri so n o f ra n d om f o re s t ，s u pport v e c tor mac hin e a n d b a ck p ropag a tio n n e ur a l n etw o r k f orelect ro n ic t ongue d ata c la ss if ic atio n ： a p pl i catio n to t h e r ec o gni t io n of o r a nge be ve ra g e and C h i n es e vin eg ar[ J ] . S ens or s an d A c tu a tors B ： C he mica l， 2 01 3， 1 77(2 ) ：97 0 -98 0. D OI：10 .1 016/ j.sn b.2 0 1 2. 1 1.07 1 . ) ( [16] W U H . Z H ANG S， M A Y， e t a l. C om p ari son of mi c r obi a l commun i ti e s i n t he fermen t a ti o n s t a rt e r use d t o br e w X i aoqu li qu o r[J ]. J o u rn a l o f t h e In s tit u t e o f Br e wi ng， 201 7， 1 23 ( 1 )： 1 13 -1 2 0. D O I： 1 0 . 1 00 2/ j i b . 388. ) ( [ 1 71 C A PORASO J G， K UCZYNSK I J， S TOMBAU G H J， e t a l . QlME a l l ows a nalysi s of hig h-th roughp ut comm u n i ty seq uen c in g d ata[ J] .N a t ure Me t ho ds， 201 0， 7 (5) ： 335- 3 36 . D OI：1 0 .10 3 8 /nmeth .f .3 0 3. ) ( [1 8] C A P OR A S O J G， B I T TINGER K， B U SH M AN F D ， e t al . P yNA S T： a fle xib le t o o l f or al ign i ng seq u enc e s t o a t e m pla te a lig n men t[ J ].Bioi nfo r ma t ics，2 0 1 0， 26 ( 2)：2 6 6- 26 7 . DO I：10. 1 09 3/ b i o inf o rmati cs/btp63 6. ) ( [ 1 9 ] E DGAR R C . S e a r ch a nd c lus te ri ng order s of m ag nit ude f a s t e r t ha nBLA S T[J ] . B i oi nf o rma t ic s ， 2 0 1 0，2 6 (1 9)： 2 4 6 0 - 2 4 6 1. D OI ： 1 0. 1 0 93/ b io i n forma ti cs/ bt q 4 61 . ) ( [ 20 1 HA AS B J ， GE V E RS D. EAR L A M. e t al . C him eri c 1 6 S rRNA s e que nce f o r m a t i o n a nd d e t e cti o n i n S a n ger and 454 -p yr os eq u encedP CR a m pl icon s[ J] . G en o me R e sear c h ， 201 1，21 (3) ： 49 4 -504 . D OI：1 0. 1101 /gr .1 1273 0 . 110 . ) ( [2 11 Q UA ST C ， PR UES S E E ， Y I LMAZ P ， e t al. Th e SI LVA r i b o somal RN A g e n e databas e pro j e c t： i mpr oved d a ta p r o c e ssi ng an d web -b a sedto ol s [J] . Nu clei c Aci d s R e s e a rch ， 2 0 12 ， 41 ( 1 )：5 90 -596. DOI：10. 1 0 9 3 /nar/gks 1 219. ) ( [ 22 ] P R IC E M N . D EH A L P S ， A R K I N A P . Fas t t re e： com pu t i n g la r g e mini mu m e v o luti o n tr e e s w i t h pr o f i les ins tea d o f a d i sta n ce m a tr i x[ J] . Mo l e cula r Bio lo g y a n d E v o l u ti o n，20 09 ， 2 6( 7) ：16 4 1 -1 6 50 .D O I：1 0 .10 93/molbev/msp077. ) ( [ 23] Z HAN G J， G U O Z ， X U E Z . e t a l . A p h y lo-f u n cti on a l cor e of gu t m ic r o bio ta i n h ealt h y yo u n g C h i ne s e c ohor t s acr os s life s tyl e s，geogr a p h y an d eth nic ities[ J ]. T he ISME J ou r n a l， 2 0 1 5， 9 ( 9 ) ：1- 1 2 . D O I：1 0. 1 038 / ism e j . 2 0 15. 11. ) ( [2 4] LO ZU PONE C ， KN I G HT R. Un i Fr a c： a n ew p hy l o genet i cmethod fo r c om par ing m i c r ob i a l c o mm unitie s [ J]. App li ed a n dE n vi ro n m e n t a l Mic rob i o l ogy， 2 0 05，7 1 ( 12)： 822 8 -82 3 5. D O I： 1 0 . 11 2 8/ A EM.7 1. 12 .8 2 2 8- 82 3 5 . 20 0 5 . ) ( [2 5] 郭壮，汤尚文， 王 玉荣，等.基于电子舌技术的襄阳市售米酒滋味 品质评价[].食品 工 业科技，2015，36( 1 5 ) ：28 9- 2 9 3 . D OI ：10.133 86 / i .is s n10 0 2-0 3 06 .2 0 1 5 .1 5. 05 2. ) ( [26 ] WA NG T， C AI G， QI U Y ， e t a l . S t r uct u r al se g r ega tio n of gu t m i c r obio ta b etw e e n c olor e ct a l ca ncer pa ti en t s an d health y v o l u nte er s J] . Th e ISMEJ ou r na l ， 2 0 1 2 ，6 (2 ) ：32 0- 32 9. DO I ： 10 . 1 0 3 8/ i s m ej . 2 0 1 1 .1 09. ) ( [ 2 7 ) Q U A I L M A， S M ITH M， C OUP L AN D P ， et a l . A t al e o f th re e n ex t g e n e r ation s e q ue n cing p l a t fo rms ： c om p arison of I o n To r rent， P a c ificB i o s c ien ces an d Il l umina MiSeq se q ue n cers [ J]. BMC G e n om i c s ， 201 2 ，13( 1) ：34 1 . D OI：1 0.1 18 6/14 7 1-216 4-1 3 -34 1 . ) ( [ 28] Q U A N L S ， T S A O M . B io c on t r o l o f s p o i l age y ea st s an d mo u ld s b y W i lliop s is s a turnus v a r . s a turnus i n y o ghu r t[J] . N u tri t i o n & F o o d S c i e nce， 2 01 0.4 0 (2)： 1 66 - 175 .DOI ： 10.11 08 / 003 46 6 51 01 1 0 2 919 2 . ) ( [ 2 9] C O RS E T T I A ， S E TT AN N I L . L ac t ob a cilli i n so u rdoug h f e rm en tati o n [J ].F o od Re search In t e r na tiona l， 2007 ， 40( 5 )： 5 3 9 -5 5 8. D OI ：10 . 10 16 / j.f o od re s.2 006.11.0 0 1. ) ( [ 30 ] KO Z I C H JJ， W EST C OT T S L ， B A XTER N T， e t al. D ev e l o pme nt o f a d u a l - i ndex se q u en cing s trategy a n d cur at i on p ipel in e f o r an a lyzing a mp l ico n se que n c e d a ta o n th e M iS eq I l l u m in a s eq u e ncin g p la t f o rm [ J ]. A pp l ie d a nd E n v ir on men ta l M ic rob i olog y ， 2 01 3 ， 7 9 (17 )：511 2 - 51 20 . DOI： 10 . 11 28 /A EM.0 10 4 3-13. ) 鲊广椒真菌多样性及其对滋味品质影响的评价《食品科学》王玉荣1代凯文2沈馨1董蕴1周书楠1郭壮11. 湖北文理学院化学工程与食品科学学院鄂西北传统发酵食品研究所2. 当阳市食品药品监督管理局 摘  要：采用Mi Seq高通量测序技术对鲊广椒真菌微生物的群落结构进行解析的基础上，使用电子舌对不同样品各滋味指标的相对强度进行测定，最终对真菌多样性与鲊广椒滋味品质进行关联性分析。结果发现，在门水平上，Ascomycota(子囊菌门)的平均相对含量高达99.97%。在属水平上，相对含量大于1.0%的真菌属分别为Candida(念珠菌属)、Eurotium(曲霉菌属)、Pichia(毕赤酵母属)、Fusarium(镰刀霉属)、Galactomyces(地霉属)、Cladosporium(分枝孢子菌属)、Debaryomyces(德巴利氏酵母属)和Guehomyces(久浩酵母属)，其平均相对含量分别为56.63%、10.12%、11.19%、2.67%、2.48%、2.21%、1.50%和1.01%。经皮尔森相关性分析发现优势真菌属与鲊广椒的苦味、涩味、后味A和后味B呈正相关。由此可见，真菌含量过高不利于鲊广椒滋味品质的形成。 关键词：鲊广椒； 高通量测序； 电子舌； 多样性； 品质；