基于Illumina Miseq测序技术分析剁辣椒产地预加工过程中菌群变化

陈梦娟; 蒋立文; 刘洋; 王蓉蓉; 周辉; 赵玲艳; 覃业优

doi:10.13386/j.issn1002-0306.2021110259

基于Illumina Miseq测序技术分析剁辣椒产地预加工过程中菌群变化

陈梦娟^1,,
蒋立文^{1, 2, ,},
刘洋¹,
王蓉蓉¹,
周辉^{1, 2},
赵玲艳¹,
覃业优³

1.
湖南农业大学食品科学技术学院，湖南长沙 410128
2.
食品科学与生物技术湖南省重点实验室，湖南长沙 410128
3.
湖南坛坛香生物科技有限公司，湖南长沙 410128

基金项目: 湖南省重点研发计划项目（2020NK2027）；湖南省研究生科研创新项目（CX20200662）。

详细信息

作者简介:
陈梦娟（1992−）（ORCID：0000-0001-9801-1215），女，博士研究生，研究方向：园艺产品采后科学与技术，E-mail：715657977@qq.com

通讯作者:
蒋立文（1968−）（ORCID:0000-0002-8536-8289），男，博士，教授，研究方向：食品生物技术，E-mail：hnndjlw@163.com

中图分类号: TS201.3
计量
- 文章访问数: 172
- HTML全文浏览量: 36
- PDF下载量: 13
出版历程
- 收稿日期: 2021-11-21
- 网络出版日期: 2022-06-12
- 刊出日期: 2022-08-14

Analysis of Microorganism Changes in Origin Pre-processing of Chopped Pepper Based on Illumina Miseq Sequencing Technology

1.
College of Food Science and Technology, Hunan Agricultural University, Changsha 410128, China
2.
Hunan Key Laboratory of Food Science and Biotechnology, Changsha 410128, China
3.
Hunan Tantanxiang Biotechnology Co., Ltd., Changsha 410128, China

摘要

摘要: 基于Illumina Miseq测序技术，对山东、山西的朝天椒（Capsicum annuum L.）和线椒（Capsicum frutescens L.）产地加工过程中细菌16S rDNA和真菌ITS进行测序，解析其细菌及真菌群落结构多样性变化，并探讨不同产地加工工艺的合理性和科学性。结果表明：两种辣椒产地原料上附着的细菌和真菌菌群差异显著，经盐渍加工后，两种辣椒中细菌菌群仍存在显著差异（P<0.05），而真菌菌群差异性则随着加工的进行逐步减小。根据可培养微生物计数结果发现，经两次清洗后，线椒样品中可培养细菌和真菌数量分别增至2.6×10⁶和1.6×10⁶ CFU/g，而朝天椒的可培养细菌和真菌数量分别降至9.0×10⁴和6.6×10³ CFU/g，说明相同设备用于清洗线椒时的适用性低，导致清洗过程中微生物积累。拌盐后的朝天椒和线椒样品中，可培养细菌和真菌数量降至10²~10⁴ CFU/g。该研究揭示了盐渍辣椒加工中微生物的群落结构和多样性变化规律，为提升产地加工水平、设备优化及品质安全与控制提供科学依据。
- 剁辣椒 /
- 预加工 /
- Illumina Miseq测序技术 /
- 细菌多样性 /
- 真菌多样性
Abstract: Based on the Illumina Miseq sequencing technology, the bacterial and fungal flora of the fresh Chaotian pepper (Capsicum annuum L.) and the line pepper (Capsicum frutescens L.) in Shandong and Shanxi province were sequenced during origin pre-processing. The diversity of bacterial and fungal community structure was analyzed, and the rationality and scientific of processing technology in different producing areas were also discussed. The results showed that the bacterial and fungal flora attached to the raw materials of the two red peppers were significantly different (P<0.05). After salting process, there were still significant differences in bacterial flora between the two kinds of pepper, while the differences in fungal flora gradually decreased with the progress of processing. According to the results of the count of cultivable microorganisms, the number of cultivable bacteria and fungi in the line pepper samples increased to 2.6×10⁶ and 1.6×10⁶ CFU/g after two-stage cleaning. However, in Chaotian pepper samples, the number of culturable bacteria and fungi after cleaning decreased to 9.0×10⁴ and 6.6×10³ CFU/g, respectively. It indicated that the same equipment was used to clean the line pepper with low applicability, resulting in the accumulation of microorganisms during the cleaning process. After salt mixing, the number of culturable bacteria and fungi decreased to 10²~10⁴ CFU/g. This study would reveal the community structure and diversity of microorganism during the primary processing of fresh red pepper, and provide a scientific basis for improving the processing level, the optimization of equipment, and quality control.
- chopped pepper /
- pre-processing /
- Illumina Miseq sequencing technology /
- bacterial community diversity /
- fungal community diversity

HTML全文

辣椒（Capsicum annuum L.）为木兰纲茄科辣椒属草本植物果实，原产于南美洲的热带地区。目前，辣椒已成为世界上最重要的农业经济作物之一，尤其是在中国和墨西哥^[1]。而中国是世界上第一大辣椒生产国和消费国^[2]。辣椒自明朝传入我国以来，其接受程度、加工方式和食用方法逐步形成了地域性差异。剁辣椒是湖南的特色辣椒加工产品，传统加工一般是以新鲜红辣椒为原料，经清洗、晾干、去梗、剁碎，拌料及密封盐渍发酵而成的产品，以香辣咸鲜而著名。

近年来，我国剁辣椒产业的发展趋势稳定，产业产值稳步增加，但其产业的发展却受到许多限制。首先，我国辣椒的主产区集中于山东、河北、山西等地北方地区，但南方地区才是主流的消费地区^[3-5]。由于新鲜辣椒果实的含水量较高，采收后辣椒易霉烂变质，且在运输过程中还存在自然蒸发脱水引起的果实失重、失鲜、机械损伤以及生理性病害等问题^[6]，损耗大，进而导致剁辣椒的生产成本升高。其次，由于辣椒为季节性蔬菜，因此生产时还需要考虑原料的季节性问题。此外，剁辣椒产业的发展还受到制作工艺与加工机械的限制。为了产业链对接紧密，一般剁辣椒加工采取异地收购辣椒并产地腌渍处理，再输送到生产企业加工等待发酵和加工，盐渍辣椒半成品更便于运输及贮藏，满足维持工厂全年生产的需求。因此，产地处理效果好坏直接影响到产品品质。

目前，剁辣椒厂家在辣椒产地预加工方式是以18%~25%的盐度盐渍新鲜红辣椒制成大量的盐渍辣椒半成品，装袋密封备用。但绝大部分厂家在产地加工过程中存在生产工艺落后、经验性加工、原料未清洗、粗放管理、添加剂使用不规范等问题^[7]，直接影响盐渍辣椒半成品的质量与产品的后续生产。此外，高盐盐渍辣椒在后续生产时需进行脱盐处理，产生大量高盐废水，给环保带来极大压力。随着人们对食品安全和环境保护的高度重视，产地加工的卫生加工条件引起生产企业的极大关注。并且根据国家对农产品加工环节减损增效的要求^[8]，辣椒产地加工必将成为未来的发展趋势。因此，对辣椒产地加工工艺优化、建立标准化管理、辣椒保脆保鲜方法研究和配套加工设备优化等方面的研究显得尤为重要。剁辣椒生产中最重要的盐渍辣椒半成品的产地预加工环节缺乏加工方法考察、配套设备适用性考察、技术规范、食品安全性以及产品品质变化等多方面的研究。

目前，规范企业在产地加工标准化、连续化方面做了大量的工作，也开发非标准设备包括：分选、清洗、振动脱水、分流斩拌、分级、拌料、包装等连续性较强的半自动化生产线，但具体控制要求尚不明晰。本文针对不同产地盐渍辣椒半成品产地预加工过程中影响品质和安全的生物学问题开展研究，采用Illumina MiSeq高通量测序技术对山东、山西不同产地及不同品种辣椒样品在预加工为盐渍辣椒过程中的细菌和真菌菌群进行分析，考虑辣椒产地预加工中关键的工艺步骤对微生物的影响，探讨加工过程中的工艺合理性和科学性，为提升产地加工水平以及工艺及设备改进提供科学依据。

1. 材料与方法

1.1 材料与仪器

新鲜辣椒　源于山东曹县和山西忻州辣椒种植加工基地，品种分别是线椒和朝天椒。样品采集时间为2020年9月中旬，两种红辣椒均为新鲜、成熟状态。样品采集见表1；平板计数琼脂培养基（PCA培养基）、孟加拉红琼脂培养基　广东环凯微生物科技有限公司；食品级焦亚硫酸钠　山东凯龙化工科技发展有限公司；异抗坏血酸钠　湖南明瑞祥盛贸易有限公司；柠檬酸　山东顾诚化工科技有限公司；食盐　中盐皓龙盐化有限责任公司；Pusion Hot start flex 2X Master Mix　上海仪涛生物仪器有限公司；DL2000 DNA Maker　宝日医生物技术（北京）有限公司；Gene colour核酸染料　北京金博益生物技术有限公司；Qubit dsDNA HS Assay Kit　赛默飞世尔科技（中国）有限公司；Biowest Agarose G-10　西班牙BIOWEST；50×TAE Buffer　上海生工生物工程股份有限公司；AxyPrep PCR Cleanup Kit AXYGEN　美国Bio-Rad Laboratories Inc.；Soil DNA Kit　美国OMEGA bio-tek Inc.。

表 1 产地加工辣椒样品采集信息

Table 1. Collection information of processed pepper samples from origin

种类	产地	品种	样品	样品名称
朝天椒	山东省菏泽市曹县	艳红	原料	CYL
			一级清洗	CQX1
			二级清洗	CQX2
			分切过筛	CFQ
			拌盐	CBY
线椒	山西省忻州市	线椒3号	原料	XYL
			一级清洗	XQX1
			二级清洗	XQX2
			分切过筛	XFQ
			拌盐	XBY

下载: 导出CSV

| 显示表格

5424型常温离心机　Eppendorf Centrifuge；Microfuge 22R型冷冻离心机　美国Beckman Coulter公司；WH-861型旋涡震荡仪　太仓市华利达实验设备有限公司；A200型PCR仪　杭州朗基科学仪器有限公司；MiSeq pe300测序仪　Illumina公司；Tanon-2500型凝胶成像仪　上海天能科技有限公司；JY300型电泳仪　北京君意东方电泳设备有限公司；DW-HL388型超低温冷冻储存箱　中科美菱低温科技有限责任公司。

1.2 实验方法

1.2.1 样品的处理和采集

新鲜红辣椒经采收后送至当地工厂采用相同加工流程与机械，具体操作见图1。取样时间点分别为原料采收、一级清洗、二级清洗、切分过筛及拌料（盐渍半成品）步骤，完成后立即取样，每个环节随机取样4袋，每袋1 kg，低温保藏运输至实验室进行试验。

图 1 辣椒产地预加工流程图

Figure 1. Flow chart of pepper origin pre-processing

下载: 全尺寸图片幻灯片

1.2.2 细菌及真菌计数

霉菌和酵母菌的菌落计数根据GB/T 4789.15-2016进行操作^[9]；菌落总数的测定则根据GB 4789.2-2016进行操作^[10]。进行微生物计数时，每个样品设置4个重复。

1.2.3 细菌16S rDNA序列及ITS 2序列分析

使用OMEGA的土壤DNA提取试剂盒对各关键工序的辣椒样本的总DNA进行提取，提取步骤均按照试剂盒说明书完成，然后将总DNA加入到50 μL缓冲液中并储存在−80 ℃备用。

扩增细菌16S rDNA V3~V4区域。引物为：338F（5'-ACTCCTACGGGAGGCAGCAG-3'）和806R（5'-GGACTACHVGGGTWTCTAAT-3'）。PCR扩增体系：25 μL包含DNA模板25 ng、上下游引物各为2.5 μL、12.5 μL PCR预混合物。PCR扩增条件：98 ℃预变性30 s；98 ℃变性10 s，54 ℃退火30 s，72 ℃延伸45 s，循环35次；最后72 ℃延伸10 min。

扩增真菌ITS 2区所选用的引物为：ITS1FI2（5'-GTGARTCATCGAATCTTTG-3'）和ITS2（5'-TCCTCCGCTTATTGATATGC-3'）。PCR体系为25 μL包含DNA模板25 ng、上下游引物各为2.5 μL、12.5 μL PCR预混合物。PCR条件：98 ℃预变性30 s；98 ℃变性10 s，52 ℃退火30 s，72 ℃延伸45 s，循环35次；最后72 ℃延伸10 min。

取PCR产物10 μL用2%琼脂糖凝胶电泳进行检测。将检测合格的样品送至杭州联川生物技术有限公司进行高通量测序。按照生物学样品处理基本要求，所有样品均有4个技术重复样品。

1.2.4 测序数据分析

对下机数据进行拆分，使用PEAR合并成对端读数，对原始标签进行质量过滤，以获得符合FastQC（V 0.10.1）的高质量tags^[11]。通过Verseach（V2.3.4）将具有97%以上相似性的序列分配给相同的可操作分类单元（operational taxonomic unit，OUT）^[12-13]。利用PyNAST软件对不同类群优势种群的差异进行多序列比对。以Chao 1、Shannon、Simpson和Observed species为指标，用QIME（1.8.0版）计算所有样品的各项指标并进行聚类，来分析各个样品中的Alpha多样性，以及进行主成分分析和主坐标分析分析^[14-16]。同时，根据每个OTU代表序列与RDP数据库和Unite数据库的比对，得到各个样本所有OTU的物种注释^[17]，经分类和统计后以R语言绘制各分类水平物种丰度表及热图。

1.3 数据处理

采用SPSS 22.0的Duncan多区间检验方法进行单因素方差分析。数据绘图使用了Origin 2018。利用OmicStudio工具进行生物信息学分析（https://www.omicstudio.cn/tool）。平均值和标准差通过使用Excel 2019从3个独立实验获得的数据并进行计算。

2. 结果与分析

2.1 线椒和朝天椒加工过程中细菌群落结构的变化

2.1.1 加工辣椒样品细菌α-多样性

本次16S测序共获得的3035040条有效数据，进行OTU聚类分析后两种产地加工辣椒中的细菌群落多样性指数结果见表2。由表2可知，Goods coverage指数均≥0.99，测序结果可反映样品中细菌群落的真实情况。由α-多样性结果可知，线椒样品中的Observed Species显著高于朝天椒中的（P<0.05）。在朝天椒样品中，根据Observed Species、Shannon以及Simpson指数可以看出，细菌的物种数量、丰富度及均匀度在清洗两次的过程中呈现显著降低趋势（P<0.05）。但在分切过筛的过程中Observed Species又呈现升高的趋势，可能与相应设备中残留的细菌有关。而在线椒样品中，在原料清洗经第二次后，Observed Species和Shannon指数均显著上升（P<0.05），说明线椒的第二次清洗没有达到清洗的效果。但随着线椒在分切过筛和拌盐过程，Shannon和Simpson逐渐降低，表明在此过程中细菌菌群的丰富度和均匀度是逐渐降低的，且变化规律与朝天椒中的相似。

表 2 加工辣椒样品中细菌群落多样性指数统计

Table 2. Alpha diversity index of bacterial community in processed pepper samples

Simple	Observed Species	Shannon	Simpson	Chao1
CYL	199±28^bc	4.73±0.25^a	0.93±0.01^ab	237±13^bc
CQX1	153±16^d	3.30±0.28^d	0.80±0.04^c	154±16^d
CQX2	151±27^d	4.07±0.30^bc	0.87±0.04^abc	151±27^d
CFQ	195±8^c	4.89±0.09^a	0.93±0.01^a	195±8^c
CBY	207±10^bc	4.67±0.37^ab	0.91±0.04^ab	207±9^bc
XYL	336±34^a	5.15±0.17^a	0.94±0.01^a	337±34^a
XQX1	183±50^cd	3.81±1.00^cd	0.82±0.18^bc	183±50^cd
XQX2	320±28^a	4.57±0.35^ab	0.87±0.05^abc	320±28^a
XFQ	311±19^a	4.95±0.21^a	0.93±0.01^a	311±19^a
XBY	237±13^b	4.56±0.14^ab	0.92±0.01^ab	237±13^b
注：表中数值表示为4次重复的平均值±标准差；同一列中不同字母表示存在显著性差异，P<0.05；表3同。

下载: 导出CSV

| 显示表格

2.1.2 加工辣椒样品细菌β-多样性

细菌群落的主成分分析（Principal Component Analysis, PCA）见图2。由图2可知，朝天椒样品点与线椒的样品点的距离较远，说明在两种辣椒在相同的加工流程中细菌菌群结构差异大，这种差异主要与朝天椒（CYL）和线椒（XYL）原料所附着细菌菌群的差异相关。此外，随着加工的进行，朝天椒样品中的细菌菌群结构的变化明显大于线椒样品中的。

图 2 加工辣椒样品中细菌群落的主成分分析

Figure 2. Principal component analysis of bacterial community in processed pepper samples

下载: 全尺寸图片幻灯片

2.1.3 辣椒样品中细菌物种注释及分类

图3A和图3B分别为产地加工辣椒样品中细菌在门水平和属水平上的物种注释结果。在朝天椒和线椒的样品中门水平上共注释到17个菌门，主要优势菌均为厚壁菌门（Firmicutes）和变形菌门（Proteobacteria），在所有样品中这两门类相对丰度之和均超过了95%。而厚壁菌门（Firmicutes）和变形菌门（Proteobacteria）也经常作为优势门类出现在发酵辣椒制品中^[18-19]。此外，还有拟杆菌门（Bacteroidetes）、放线菌门（Actinobacteria）、蓝细菌（Cyanobacteria）等。

图 3 门和属分类水平上细菌菌群相对丰度

Figure 3. Relative abundance of bacterial community at the phylum and genus classification level

下载: 全尺寸图片幻灯片

由图3B可知，在属水平上原料上优势菌属种类差异明显，朝天椒原料附着的优势菌为乳球菌属（Lactococcus，相对丰度为14.21%）、假单胞菌属（Pseudomonas，13.30%）、芽孢杆菌属（Bacillus，12.39%）、不动杆菌属（Acinetobacter，10.37%）以及肠球菌属（Enterococcus，11.02%）。而线椒原料表面附着的优势菌为假单胞菌属（Pseudomonas，18.48%）、泛菌属（Pantoea，13.94%）、肠球菌属（Enterococcus，13.72%）和沙雷氏菌属（Serratia，10.36%）。

其次，在经过两次清洗后，朝天椒样品中乳球菌属（Lactococcus）和芽孢杆菌属（Bacillus）以及线椒样品中的肠球菌属（Enterococcus）的相对丰度显著增加（P<0.05）。但在对比CQX2和CFQ后发现，朝天椒样品中乳球菌属（Lactococcus）和芽孢杆菌属（Bacillus）的相对丰度分别降低至8.19%和2.93%。线椒样品中肠球菌属（Enterococcus）的相对丰度也由二次清洗后的29.91%降低至14.26%，这可能与二次清洗后的振动脱水工序有关。

在朝天椒的半成品（CBY）中，丛毛单胞菌属（Comamonas）、肠球菌属（Enterococcus）和假单胞菌属（Pseudomonas）3种菌属相对丰度最高；而线椒半成品（XBY）中假单胞菌属（Pseudomonas，23.17%）、寡养单胞菌属（Stenotrophomonas，23.20%）和肠球菌属（Enterococcus，20.12%）为优势菌，均与原料中优势菌种类和结果有较大差异。本研究中XBY的细菌注释结果与赵玲艳等^[20]采用454焦磷酸测序对同品种同产地的盐渍线椒中细菌菌群部分相似。此外，两种辣椒样品前30的菌属中，常见的乳酸菌菌属有乳球菌属（Lactococcus）和明串珠菌属（Leuconostoc），但盐渍半成品中乳酸菌的相对丰度均低于原料。

2.2 线椒和朝天椒产地加工过程中真菌群落结构的变化

2.2.1 产地加工辣椒样品真菌α-多样性

本次ITS测序共获得的3337111条有效数据，进行分析后两种产地加工辣椒中的真菌群落多样性指数结果见表3。由真菌α-多样性的结果可知，线椒与朝天椒原料上的真菌种类数量上无明显差异。并且两种辣椒中真菌的Observed Species均呈现出先降低后增高的规律，与16S的结果呈现出的规律一致。根据Observed Species、Shannon以及Simpson指数可以看出，朝天椒和线椒样品中微生物的多样性及均匀度在第二次清洗后显著降低（P<0.05）。在分切过筛时，CFQ和XFQ样品中真菌Observed Species显著高于CQX2和XQX2的（P<0.05）。

表 3 加工辣椒样品中真菌群落多样性指数统计

Table 3. Alpha diversity index of fungi community in processed pepper samples

Simple	Observed Species	Shannon	Simpson	Chao1
CYL	152±6^a	3.61±0.31^bc	0.81±0.02^abc	154±7^a
CQX1	136±8^ab	4.11±0.05^ab	0.86±0.00^ab	137±7^ab
CQX2	66±1^e	0.43±0.17^e	0.08±0.03^e	70±10^d
CFQ	137±12^ab	3.96±0.23^ab	0.85±0.02^ab	140±12^ab
CBY	121±17^bc	3.13±0.80^cd	0.73±0.12^cd	125±19^bc
XYL	138±7^ab	4.44±0.39^a	0.89±0.04^a	140±9^ab
XQX1	74±4^e	3.77±0.07^b	0.87±0.01^ab	74±4^d
XQX2	82±11^de	2.98±0.37^d	0.69±0.09^d	83±11^d
XFQ	106±17^c	3.93±0.31^ab	0.86±0.05^ab	109±17^c
XBY	104±35^cd	2.95±0.44^d	0.79±0.05^bc	107±37^c

下载: 导出CSV

| 显示表格

2.2.2 产地加工辣椒样品真菌β-多样性

在真菌群落的PCA分析中（见图4），两坐标轴的贡献率之和为70.55%。由图4可知，朝天椒（CYL）和线椒（XYL）原料附着的真菌也具有显著差异（P<0.01），但随着加工的进行，两种辣椒中的真菌菌群结构的差异性在逐渐的减小。

图 4 加工辣椒样品中真菌群落的主成分分析

Figure 4. Principal component analysis of fungi community in processed pepper samples

下载: 全尺寸图片幻灯片

2.2.3 产地加工辣椒样品中真菌物种注释及分类

图5A和图5B分别为产地加工辣椒样品中真菌在门水平和属水平上的物种注释结果。在门水平上共注释到6个菌门，子囊菌门（Ascomycota）和担子菌门（Basidiomycota）为绝对优势菌门，所有样品中二者相对丰度的总和均超过99.7%。此外，还检测到接合菌门（Zygomycota）、被孢霉菌门（Mortierellomycota）、壶菌门（Chytridiomycota）以及Fungi unclassified。

图 5 门和属分类水平上真菌菌群相对丰度

Figure 5. Relative abundance of fungi community at the phylum and genus classification level

下载: 全尺寸图片幻灯片

由图5B可知，朝天椒原料上真菌的优势菌属为拟酵母属（Pseudozyma，相对丰度为36.24%）、Gibellulopsis（相对丰度为13.44%）以及链格孢属（Alternaria，12.13%）。而线椒原料上的优势菌为链格孢属（Alternaria，相对丰度29.39%）、织球壳菌属（Plectosphaerella，19.43%）、腥掷孢菌属（Tilletiopsis，11.63%）。

在朝天椒的清洗二次（CQX2）样品中，链格孢属（Alternaria）的相对丰度高达96.05%，并且在线椒清洗二次（XQX2）的样品中，该菌属的相对丰度也增高至55.51%。链格孢属（Alternaria）是植物上常见的腐生菌，易引果实细胞壁坍塌且色泽变暗，对辣椒保脆保鲜不利^[21-22]。在朝天椒和线椒的半成品中，链格孢属（Alternaria）的相对丰度仍然较高，分别为29.62%和24.20%。此外，还存在织球壳菌属（Plectosphaerella）、赤霉菌属（Gibberella）以及镰刀菌属（Fusarium）等植物致病菌或致腐真菌属^[23-25]。

2.3 可培养微生物的变化

从图6可以看出，在加工过程中细菌与霉菌数量级在10²~10⁶ CFU/g之间。在盐渍朝天椒半成品中可培养细菌数量由原料的2.2×10⁵降至1.2×10⁴ CFU/g，可培养真菌则由原料中8.3×10³降至8.4×10² CFU/g，但第一次清洗后样品中的在由8.3×10³升高至1.6×10⁴ CFU/g，第二次清洗后降至6.6×10³ CFU/g。而在线椒原料在两次清洗后，可培养细菌数量由2.6×10⁴升高至2.6×10⁶ CFU/g，并且真菌数量也由1.4×10⁴上升至1.6×10⁶ CFU/g。在拌盐后可培养细菌和真菌的数量才分别降低至1.4×10³和2.3×10³ CFU/g。结合菌落计数和测序结果发现，现有机械的清洗效果较差，尤其是在红线椒加工中，这可能是由于两次清洗均使用自来水冲洗且进水和出水的速率不一致导致在清洗环节中辣椒表明微生物出现积累。其次，线椒果实较长且表面积更大，因此同样的清洗工艺无法满足不同品种原料的清洁要求，需要采取延长清洗工序时间或改善清洗过程中换水的频率才有利于微生物降低。

图 6 加工辣椒样品中可培养微生物的变化

注：图中不同小写字母表示差异显著，P<0.05。

Figure 6. Changes of culturable microorganisms in samples of processed pepper

下载: 全尺寸图片幻灯片

3. 讨论

本研究利用高通量测序探讨山东、山西的朝天椒和线椒在盐渍辣椒加工过程中微生物的变化规律，结果表明两种辣椒产地原料上附着的细菌和真菌菌群差异显著，经盐渍加工后，两种辣椒中细菌菌群仍存在显著差异（P<0.05），而真菌菌群差异性则随着加工是进行逐步减小，且细菌和真菌的丰富度、均匀度均极显著降低（P<0.01）。根据计数结果，在盐渍朝天椒半成品中可培养细菌数量由原料的2.2×10⁵降至1.2×10⁴ CFU/g，可培养真菌则由原料中8.3×10³降至8.4×10² CFU/g；而盐渍线椒样品中可培养细菌和真菌的数量分别降至1.4×10³和2.3×10³ CFU/g。预加工过程中两种辣椒样品中微生物数量呈现下降趋势，产地预加工效果较好。

剁辣椒的最大特点就是脆辣酸爽，产地加工是保持辣椒原料脆度是重要前提，除品种差异外，微生物分泌果胶酶和内源酶对辣椒果实软化也是影响脆度的重要因素，因而，在产地加工过程中，降低微生物数量是必要的，同时也降低了条件致病菌带来的风险。本研究结果表明，相较于朝天椒样品中菌群和计数结果的变化，线椒样品在两次清洗的效果较差，清洗后反而出现微生积累现象，需采取延长清洗工序时间或改善清洗过程中换水的频率等方法来减少微生物数量。所以，现有的产地加工工艺是可行的，但应结合辣椒原料的特点进行差异化设计，并关注清洗阶段的效果可能对后期加工的影响，进行下一步研究和探讨。若产地预加工清洗处理可有效降低微生物污染，则可适当降低拌入的食盐含量，对后期加工及环境保护具有重要的影响。从产地预加工环节控制辣椒的质量，充分研究辣椒加工过程的内在品质变化，可为稳定产品品质提供支撑。本文主要针对产地预加工环节中微生物多样性及数量上进行研究，不同盐度及产地预处理对产品品质的影响则需要进一步研究、验证。

图 1 辣椒产地预加工流程图

Figure 1. Flow chart of pepper origin pre-processing

下载: 全尺寸图片幻灯片

图 2 加工辣椒样品中细菌群落的主成分分析

Figure 2. Principal component analysis of bacterial community in processed pepper samples

下载: 全尺寸图片幻灯片

图 3 门和属分类水平上细菌菌群相对丰度

Figure 3. Relative abundance of bacterial community at the phylum and genus classification level

下载: 全尺寸图片幻灯片

图 4 加工辣椒样品中真菌群落的主成分分析

Figure 4. Principal component analysis of fungi community in processed pepper samples

下载: 全尺寸图片幻灯片

图 5 门和属分类水平上真菌菌群相对丰度

Figure 5. Relative abundance of fungi community at the phylum and genus classification level

下载: 全尺寸图片幻灯片

图 6 加工辣椒样品中可培养微生物的变化

注：图中不同小写字母表示差异显著，P<0.05。

Figure 6. Changes of culturable microorganisms in samples of processed pepper

下载: 全尺寸图片幻灯片

表 1 产地加工辣椒样品采集信息

Table 1 Collection information of processed pepper samples from origin

种类	产地	品种	样品	样品名称
朝天椒	山东省菏泽市曹县	艳红	原料	CYL
			一级清洗	CQX1
			二级清洗	CQX2
			分切过筛	CFQ
			拌盐	CBY
线椒	山西省忻州市	线椒3号	原料	XYL
			一级清洗	XQX1
			二级清洗	XQX2
			分切过筛	XFQ
			拌盐	XBY

下载: 导出CSV

表 2 加工辣椒样品中细菌群落多样性指数统计

Table 2 Alpha diversity index of bacterial community in processed pepper samples

Simple	Observed Species	Shannon	Simpson	Chao1
CYL	199±28^bc	4.73±0.25^a	0.93±0.01^ab	237±13^bc
CQX1	153±16^d	3.30±0.28^d	0.80±0.04^c	154±16^d
CQX2	151±27^d	4.07±0.30^bc	0.87±0.04^abc	151±27^d
CFQ	195±8^c	4.89±0.09^a	0.93±0.01^a	195±8^c
CBY	207±10^bc	4.67±0.37^ab	0.91±0.04^ab	207±9^bc
XYL	336±34^a	5.15±0.17^a	0.94±0.01^a	337±34^a
XQX1	183±50^cd	3.81±1.00^cd	0.82±0.18^bc	183±50^cd
XQX2	320±28^a	4.57±0.35^ab	0.87±0.05^abc	320±28^a
XFQ	311±19^a	4.95±0.21^a	0.93±0.01^a	311±19^a
XBY	237±13^b	4.56±0.14^ab	0.92±0.01^ab	237±13^b
注：表中数值表示为4次重复的平均值±标准差；同一列中不同字母表示存在显著性差异，P<0.05；表3同。

下载: 导出CSV

表 3 加工辣椒样品中真菌群落多样性指数统计

Table 3 Alpha diversity index of fungi community in processed pepper samples

Simple	Observed Species	Shannon	Simpson	Chao1
CYL	152±6^a	3.61±0.31^bc	0.81±0.02^abc	154±7^a
CQX1	136±8^ab	4.11±0.05^ab	0.86±0.00^ab	137±7^ab
CQX2	66±1^e	0.43±0.17^e	0.08±0.03^e	70±10^d
CFQ	137±12^ab	3.96±0.23^ab	0.85±0.02^ab	140±12^ab
CBY	121±17^bc	3.13±0.80^cd	0.73±0.12^cd	125±19^bc
XYL	138±7^ab	4.44±0.39^a	0.89±0.04^a	140±9^ab
XQX1	74±4^e	3.77±0.07^b	0.87±0.01^ab	74±4^d
XQX2	82±11^de	2.98±0.37^d	0.69±0.09^d	83±11^d
XFQ	106±17^c	3.93±0.31^ab	0.86±0.05^ab	109±17^c
XBY	104±35^cd	2.95±0.44^d	0.79±0.05^bc	107±37^c

下载: 导出CSV

参考文献(25)

[1]	FAOSTAT (Food and Agriculture Organization of the United Nation)[DB/OL]. (2019)[2021-01-23]. http://www.fao.org/faostat/en/#data/QC.
[2]	李洪. 印度魔鬼椒果实转录组测序及分析[D]. 吉林: 吉林大学, 2020 LI H. Transcriptome sequencing and analysis of Capsicum chinense[D]. Jilin: Jilin University, 2020.
[3]	甘增鹏, 蒋立文, 吴章华, 等. 剁辣椒加工的研究进展[J]. 食品安全质量检测学报,2019,10(13):4278−4283. [GAN Z P, JIANG L W, WU Z H, et al. Research progress on the chopped chilli processing[J]. Journal of Food Safety & Quality,2019,10(13):4278−4283. doi: 10.3969/j.issn.2095-0381.2019.13.042 GAN Z P, JIANG L W, WU Z H, et al. Research progress on the chopped chilli processing[J]. Journal of Food Safety & Quality, 2019, 10(13): 4278-4283. doi: 10.3969/j.issn.2095-0381.2019.13.042
[4]	赵泽阳, 赵邦宏, 乔立娟. 利用集中度分析研究河北辣椒产业发展现状[J]. 中国商贸,2019,7:211−212. [ZHAO Z Y, ZHAO B H, QIAO L J. Using concentration analysis to study the development status of Hebei pepper industry[J]. China Journal of Commerce,2019,7:211−212. ZHAO Z Y, ZHAO B H, QIAO L J. Using concentration analysis to study the development status of Hebei pepper industry[J]. China Journal of Commerce, 2019, 7: 211-212.
[5]	徐青, 田朝辉, 申爱民, 等. 柘城县辣椒产业发展现状, 问题及建议[J]. 辣椒杂志,2020(1):1−5. [XU Q, TIAN Z H, SHEN A M, et al. Current status, problems and suggestions of the hot pepper industry in Zhecheng County[J]. Journal of China Capsicum,2020(1):1−5. doi: 10.3969/j.issn.1672-4542.2020.01.001 XU Q, TIAN Z H, SHEN A M, et al. Current status, problems and suggestions of the hot pepper industry in Zhecheng County[J]. Journal of China Capsicum, 2020(1): 1-5. doi: 10.3969/j.issn.1672-4542.2020.01.001
[6]	SANTOS S, CARDOSO R, BORGES S, et al. Post-harvest losses of fruits and vegetables in supply centers in Salvador, Brazil: Analysis of determinants, volumes and reduction strategies[J]. Waste Management,2020,101:161−170. doi: 10.1016/j.wasman.2019.10.007
[7]	徐浩. 发酵脆鲜风味剁辣椒发酵剂的制备及应用[D]. 长沙: 湖南农业大学, 2013 XU H. Preparation and application of crisp and fresh flavor chopped pepper fermentation agent[D]. Changsha: Hunan Agricultural University, 2013.
[8]	农业农村部《关于促进农产品加工环节减损增效的指导意见》[J]. 中国食品, 2021(2): 149−150 Guiding opinions on promoting loss reduction and efficiency increase in agricultural product processing links issued by the Ministry of Agriculture and Rural Areas[J]. Chinese Food, 2021(2): 149−150.
[9]	国家卫生和计划生育委员会. GB 4789.15-2016 食品微生物学检验霉菌和酵母菌计数[S]. 北京: 中国标准出版社, 2016 State Health and Family Planning Commission. GB 4789.15-2016 Microbiological examination of food mould and yeast count[S]. Beijing: China Standard Press, 2016.
[10]	国家食品药品监督管理总局, 国家卫生和计划生育委员会. GB 4789.2-2016 食品微生物学检验菌落总数测定[S]. 北京: 中国标准出版社, 2016 State Food and Drug Administration, State Health and Family Planning Commission. GB 4789.2-2016 Microbiological examination of food Determination of total bacterial count[S]. Beijing: China Standard Press, 2016.
[11]	BLAXTER M, MANN J, CHAPMAN T, et al. Defining operational taxonomic units using DNA barcode data[J]. Philosophical Transactions Biological Sciences,2005,360(1462):1935−1943. doi: 10.1098/rstb.2005.1725
[12]	ROGENS T, FLOURI T, NICHOLS B, et al. VSEARCH: A versatile open source tool for metagenomics[J]. Peerj,2016,4(10):e2584.
[13]	CAPORASO J G, KUCZYNSKI J, STOMBAUGH J, et al. QIIME allows analysis of high-throughput community sequencing data[J]. Nature Methods,2010,7(5):335−336. doi: 10.1038/nmeth.f.303
[14]	LOZUPONE C A, HAMADY M, KELLEY S T, et al. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities[J]. Applied & Environmental Microbiology,2007,73(5):1576−1585.
[15]	RAMETTE A. Multivariate analyses in microbial ecology[J]. FEMS Microbiology Ecology,2010,62(2):142−160.
[16]	KÕLJALG U, LARSSON K H, ABARENKOV K, et al. UNITE: A database providing web-based methods for the molecular identification of ectomycorrhizal fungi[J]. New Phytologist,2005,166(3):1063−1068. doi: 10.1111/j.1469-8137.2005.01376.x
[17]	COLE J R, CHAI B, FARRIS R J, et al. The Ribosomal Database Project (RDP-II): Sequences and tools for high-throughput rRNA analysis[J]. Nucleic Acids Research, 2005, 33(Database issue): 294−296.
[18]	刘洋, 丁悦, 孙劲松, 等. 基于高通量测序分析四川辣椒酱自然发酵过程中细菌群落结构演替规律[J]. 食品工业科技,2020,41(24):81−86. [LIU Y, DING Y, SUN J S, et al. Succession of bacterial community structure during Sichuan chili sauce natural fermentation analyzed by high-throughput sequencing[J]. Science and Technology of Food Industry,2020,41(24):81−86. LIU Y, DING Y, SUN J S, et al. Succession of bacterial community structure during Sichuan chili sauce natural fermentation analyzed by high-throughput sequencing[J]. Science and Technology of Food Industry, 2020, 41(24): 81-86.
[19]	WANG J, WANG R, XIAO Q, et al. Analysis of bacterial diversity during fermentation of Chinese traditional fermented chopped pepper[J]. Letters in Applied Microbiology,2019,69(5):346−352. doi: 10.1111/lam.13212
[20]	赵玲艳, 黄嘉欣, 杨剑, 等. 盐渍辣椒细菌多样性分析[J]. 食品科学,2020,41(20):105−111. [ZHAO L Y, HUANG J X, YANG J, et al. Analysis of bacterial diversity of salted peppers (Capsicum annuum)[J]. Food Science,2020,41(20):105−111. doi: 10.7506/spkx1002-6630-20190829-318 ZHAO L Y, HUANG J X, YANG J, et al. Analysis of bacterial diversity of salted peppers (Capsicum annuum)[J]. Food Science, 2020, 41(20): 105-111. doi: 10.7506/spkx1002-6630-20190829-318
[21]	LORENZINI M, ZAPPARPLI G. Characterization and pathogenicity of Alternaria spp. strains associated with grape bunch rot during post-harvest withering[J]. International Journal of Food Microbiology,2014,186:1−5. doi: 10.1016/j.ijfoodmicro.2014.06.008
[22]	郑鄢燕, 赵力卉, 王宇滨, 等. 链格孢, 粉红单端孢通过降解果皮细胞壁导致哈密瓜病害[J]. 食品与发酵工业,2020,46(10):124−131. [ZHENG Y Y, ZHAO L H, WANG Y B, et al. Alternaria alternata and Trichothecium roseum cause disease of Hami melon by degrading the cell wall in the pericarp[J]. Food and Fermentation Industries,2020,46(10):124−131. ZHENG Y Y, ZHAO L H, WANG Y B, et al. Alternaria alternata and Trichothecium roseum cause disease of Hami melon by degrading the cell wall in the pericarp[J]. Food and Fermentation Industries, 2020, 46(10): 124-131.
[23]	陈婷, 谷羚毓, 韩晓云. 腌蛋盐水中致腐菌的分离鉴定[J]. 中国农学通报,2018,34(16):151−156. [CHEN T, GU L Y, HAN X Y. Separation and identification of rot pathogens in salted egg brine[J]. Chinese Agricultural Science Bulletin,2018,34(16):151−156. doi: 10.11924/j.issn.1000-6850.casb17040109 CHEN T, GU L Y, HAN X Y. Separation and identification of rot pathogens in salted egg brine[J]. Chinese Agricultural Science Bulletin, 2018, 34(16): 151-156. doi: 10.11924/j.issn.1000-6850.casb17040109
[24]	MURAD N A, MOHAMED N N, SHOHAIMI S, et al. Genetic diversity and pathogenicity of Fusarium species associated with fruit rot disease in banana across Peninsular Malaysia[J]. Journal of Applied Microbiology,2017,123:1533−1546. doi: 10.1111/jam.13582
[25]	LEE J H, OH S, KIM D S, et al. First report of leaf blight disease caused by Fusarium verticillioides (Teleomorph Gibberella moniliformis) in Kentucky Bluegrass[J]. Journal of Phytopathology,2014,162(5):345−347. doi: 10.1111/jph.12193

施引文献

资源附件(0)

图(6) / 表(3)

计量

文章访问数: 172
HTML全文浏览量: 36
PDF下载量: 13
被引次数: 0

1. 材料与方法
1.1 材料与仪器
1.2 实验方法
1.2.1 样品的处理和采集
1.2.2 细菌及真菌计数
1.2.3 细菌16S rDNA序列及ITS 2序列分析
1.2.4 测序数据分析
1.3 数据处理
2. 结果与分析
2.1 线椒和朝天椒加工过程中细菌群落结构的变化
2.1.1 加工辣椒样品细菌α-多样性
2.1.2 加工辣椒样品细菌β-多样性
2.1.3 辣椒样品中细菌物种注释及分类
2.2 线椒和朝天椒产地加工过程中真菌群落结构的变化
2.2.1 产地加工辣椒样品真菌α-多样性
2.2.2 产地加工辣椒样品真菌β-多样性
2.2.3 产地加工辣椒样品中真菌物种注释及分类
2.3 可培养微生物的变化
3. 讨论

1. 材料与方法
1.1 材料与仪器
1.2 实验方法
1.2.1 样品的处理和采集
1.2.2 细菌及真菌计数
1.2.3 细菌16S rDNA序列及ITS 2序列分析
1.2.4 测序数据分析
1.3 数据处理
2. 结果与分析
2.1 线椒和朝天椒加工过程中细菌群落结构的变化
2.1.1 加工辣椒样品细菌α-多样性
2.1.2 加工辣椒样品细菌β-多样性
2.1.3 辣椒样品中细菌物种注释及分类
2.2 线椒和朝天椒产地加工过程中真菌群落结构的变化
2.2.1 产地加工辣椒样品真菌α-多样性
2.2.2 产地加工辣椒样品真菌β-多样性
2.2.3 产地加工辣椒样品中真菌物种注释及分类
2.3 可培养微生物的变化
3. 讨论

参考文献(25)

施引文献

资源附件(0)

基于Illumina Miseq测序技术分析剁辣椒产地预加工过程中菌群变化

作者简介: 陈梦娟（1992−）（ORCID：0000-0001-9801-1215），女，博士研究生，研究方向：园艺产品采后科学与技术，E-mail：715657977@qq.com

通讯作者: 蒋立文（1968−）（ORCID:0000-0002-8536-8289），男，博士，教授，研究方向：食品生物技术，E-mail：hnndjlw@163.com

计量

出版历程