|
高通量测序分析不同地区红腐乳细菌多样性高通量测序分析不同地区红腐乳细菌多样性 徐 琼1，刘 洋1,*，曲勤凤1，窦同海2，陈羽菲1，张娜娜1，赵 磊2，钟 江2，翁史昱1，杨捷琳3，赵国屏2 （1.上海市质量监督检验技术研究院，上海 200233；2.复旦大学生命科学学院，上海 200433；3.上海海关动植物与食品检验检疫技术中心，上海 200135） 摘 要：为了解我国传统发酵食品红腐乳中细菌的多样性及其微生物安全性，采用高通量测序技术，分析华东、华北和东北地区红腐乳中细菌16S rDNA V1-V3区基因序列，比较不同地区红腐乳中微生物群落结构组成及多样性差异。结果显示，本次测序深度有效地覆盖了样品中微生物种类，不同地区红腐乳样品中细菌群落结构具有丰富多样性。在门水平上，厚壁菌门（Firmicutes）、放线菌门（Actinobacteria）和变形菌门（Proteobacteria）为不同地区红腐乳种主要优势菌的门；而在属水平上，9 个红腐乳样品共检测出296 个细菌属，其中共有菌属分别为：棒状杆菌属（Corynebacterium）、四联球菌属（Tetragenococcus）、明串珠菌属（Leuconostoc）、Rummeliibacillus属、乳球菌属（Lactococcus）、魏斯氏菌属（Weissella）、链球菌属（Streptococcus）、丛毛单胞菌属（Comamonas）、芽孢杆菌属（Bacillus）以及Trabulsiella属。不同样品中菌群的结构与地区存在一定的聚类关系：明串珠菌属、乳球菌属以及四联球菌属是华北地区样品中主要的优势菌属，而东北地区腐乳样品中的芽孢杆菌属和盐厌氧菌属为主要优势菌属。研究结果加深了对红腐乳中细菌群落组成和多样性的认识，为保证传统发酵食品的质量安全、生产工艺提供一定的理论支持。 关键词： 高通量测序；红腐乳；细菌多样性 腐乳是我国传统的发酵豆制品，至今已有一千多年的生产历史。它是用豆浆的凝乳状物经微生物发酵所制成的一种干酪型产品，故又被誉为“中国干酪”[1]。经微生物发酵后的腐乳除大豆本身所具有较高蛋白质、脂肪和大豆异黄酮[2]等营养外，还具有抗氧化[3-4]、抗高血压[5]、抗乙酰胆碱酯酶活性[6]等生理功能，深受人们的喜爱。 由于腐乳是半开放式生产，且大多为传统的手工操作，未有加热杀菌的过程，因此有可能被环境中的微生物杂菌污染[7]，从而影响产品的品质，甚至给消费者的安全带来隐患。如腐乳中常见的生物胺，主要是由微生物对游离氨基酸的脱羧作用产生[8]，诸如芽孢杆菌属（Bacillus）、梭菌属（Clostridium）、乳酸杆菌属（Lactobacillus）和假单胞菌属（Pseudomonas）均能产生脱羧酶[9-10]。此外，氨基甲酸乙酯[11]、肠毒素[12]、丙烯酰胺[13]也是常见的发酵过程中的微生物产生的副产物。 廖新浴等[14]利用纯培养和16S rDNA测序相结合对不同类型腐乳优势菌群进行了研究，发现不同腐乳中的菌群结构各不相同，同时存在优势的食源性致病菌如蜡状芽孢杆菌（Bacillus cereus）、粪肠球菌（Enterococcus faecalis）和恶臭假单胞菌（Pseudomonas putida）。由于该结果依赖于纯培养，对于含量较低和非可培养微生物可能有所遗漏，无法对样品中的所有微生物群落进行全面分析。 To verify the performance of the EM-aided dual-frequency ambiguity estimation method proposed above,a simulation experiment was designed consisting of two experimental scenarios. 近年来，随着测序技术的不断发展，高通量测序技术也被称为下一代测序技术，由于成本低、可读量和质量高，是目前评估微生物多样性的有力工具[15-16]。许多研究已经将该技术用于发酵食品中菌群多样性的分析，如豆豉[17]、豆汁[18]、大酱[19]、曲拉[20]等，但运用高通量测序技术分析华东、华北和东北地区腐乳菌群多样性的研究鲜见报道。 本研究以我国不同地区红腐乳为研究对象，采用高通量测序技术对细菌16S rDNA V1-V3区进行测序，以期更加全面地了解不同地域红腐乳中细菌菌群多样性，明确各自的优势菌群，为建立地域性特色发酵食品提供理论依据。 1 材料与方法1.1 材料与试剂红腐乳样品从我国6 个省市当地超市购买，每例样品来源于不同生产厂家，共计9 例，具体地区与样品情况见表1。 （2）多功能护理型轮椅组合床改变了人们对传统组合床的定向思维的认识，颠覆传统组合床不可实现翻身功能的设计方式．以电动千斤顶的巧妙设计，彻底解决传统组合床给病人生活带来的种种困扰并进一步完善其便捷性． 表1 样品信息
Table 1 Information about the samples tested in this study  样品编号 生产产地 地理区域 发酵类型R1 浙江绍兴 华东 霉菌型R2 北京市 华北 霉菌型R3 上海市 华东 霉菌型R4 浙江湖州 华东 霉菌型R5 江苏扬州 华东 霉菌型R6 江苏扬州 华东 霉菌型R7 吉林长春 东北 霉菌型R8 黑龙江齐齐哈尔 东北 细菌型R9 北京市 华北 霉菌型
聚合酶链式以应（polymerase chain reaction，PCR）Master Mix（M7505） 美国Promega公司；胶回收试剂盒（Code No. 9762）、引物 宝生物工程（大连）有限公司。 1.2 仪器与设备SW22振荡水浴槽 德国Julabo公司；Centrifuge 5417R冷冻离心机 德国Eppendorf公司；Nanodrop2000超微量分光光度计 美国Thermo Scientific公司；VERITI梯度PCR仪 美国ABI公司；Sub-Cell电泳仪、Geldoc XR+型凝胶成像仪 美国Bio-Rad公司；MiSeq高通量测序平台 美国Illumina公司。 1.3 方法1.3.1 核酸提取 采用十六烷基三甲基溴化铵（hexadecyltrimethylammonium bromide，CTAB）法[21]并稍作修改。取同批次样品3 份，从每份样品中倒取50 mL腐乳汤汁，共计150 mL。混合后，取混合样品汤汁5 mL，9 000 r/min离心10 min，去上清液，沉淀中加入10 mL的无菌水，在振荡器上振荡混匀，9 000 r/min离心10 min，去上清液，沉淀加入800 μL溶菌酶（10 mg/mL），37 ℃水浴1 h；然后加入40 μL的蛋白酶K和CTAB裂解液，56 ℃水浴2 h。每管裂解液中加入等量Tris饱和酚-三氯甲烷-异戊醇（25∶24∶1，V/V），混匀，13 000 r/min离心10 min。取上清液至新离心管中，加入等体积三氯甲烷-异戊醇（24∶1，V/V），混匀，13 000 r/min离心10 min。上清液移至新离心管中，加入2 倍体积无水冰乙醇，混匀，静置30 min左右，13 000 r/min离心2 min。弃上清液，用70%乙醇溶液洗涤沉淀2～3 次，得到的沉淀放置于室温条件下进行自然风干。风干后加入50 μL Tris-乙二胺四乙酸缓冲液，溶解沉淀。 其中：Q为流量，m3/min；D为水平输料管的直径，m；v为管道内气体流速，m/s。设计取D=0.036 m,v=30 m/s,则Q=1.83 m3/min。 1.3.2 PCR扩增及16S rRNA测序 （1）评价主体多元化。参与课程考核评价活动的主体除了教师外，还应该包括学生、企业行业专家、课程研究学者、专职评价机构等人员。 以提取的总细菌基因组D N A为模板，引物针对16S rRNA基因V1-V3区合成特异引物，利用27F（5’-AGAGTTTGATCCTGGCTCAG-3’）和534R（5’-ATTACCGCGGCTGCTGG-3’）引物进行PCR扩增。PCR体系：基因组模板DNA 2 μL，2×PCR Master Mix 12.5 μL，正以引物各0.5 μL，ddH2O补足至25 μL。PCR扩增条件：94 ℃预变性2 min；94 ℃变性30 s，60 ℃退火50 s；72 ℃延伸45 s，扩增40 个循环；4 ℃保存。对所得的PCR产物用1%的琼脂糖凝胶电泳进行检测，对符合要求的PCR产物进行后续实验。扩增后的PCR产物经纯化回收后，使用Illumina MiSeq平台进行高通量测序及生物信息学分析。 1.4 数据处理基于可操作分类单元（operational taxonomic units，OTU）聚类分析结果，研究环境中微生物的多样性，使用QIIME分析流程中pick_open_reference方式进行OTU聚类。聚类算法使用Uclust[23]，数据库使用Greengenes 2013-08 release（ http://greengenes.lbl.gov/Download/）版本，相似度阈值为97%，对所有序列进行OTU划分并进行生物信息统计分析，其余参数为QIIME默认参数。根据样本分析所得OTU情况使用QIIME α_diversity分析模块进行α多样性分析，其中菌群丰度可由Chao 1指数进行评估，其数值越高表明群落物种丰富度越高；菌群多样性由Shannon指数、Simpson指数进行评估，Shannon指数与群落多样性呈正相关，与Simpson指数呈负相关。Coverage指数是测序深度指标，代表每个样品文库的覆盖率。 多年来通过引进优良的山药品种，在栽培技术上突破传统山药生产中搭架难、深耕难、收获难的三大难题，使全市山药种植技术升级，扩大种植规模，经济效益明显提高，产品质量显著改善，市场竞争力上全面提升。带动全市山药生产向规模化、产业化发展，成为全市蔬菜生产上又一个说得出口、拿得出手的品牌。 使用β多样性分析进行不同样品在微生物群落构成上的比较。使用软件QIIME β_diversity_through_plots,jackknifed_β_diversity模块进行β多样性分析。在之前分析步骤中，已经将序列按照其碱基组成的相似性，划归到各个OTU中。在进行分类学分析时，首先，将每一条典型序列都与Greengenes 2013-08 release（ http://greengenes.lbl.gov/Download/）数据库进行比对，找出其最相近且可信度达80%以上的种属信息，以获得每个OTU的分类学信息。根据物种注释，统计每个样品在分类水平上的序列数目，分类水平分别是门水平（将相对丰度＜0.1%及unidentified归入others）和属水平（相对丰度＜5%及unidentified归入others）。使用软件QIIME，算法为rdp_classifier[24] v2.2。热图使用R语言的“pheatmap”包绘制。基于菌群之间Unifrac[25]算法加权的主坐标分析（principal co-ordinates analysis，PCoA）可以实现多个样品的分类。 α多样性以映样品内部微生物群落的物种丰富度，根据97%相似性水平下的OTU信息，采用α多样性指标的Shannon指数、Simpson指数和Chao 1指数对样品微生物物种的丰富度和多样性进行评估，结果见表2。所有样本的测序Coverage指数大于0.9，说明腐乳样品文库种序列基本都被测出，未被测到的概率较低，测序结果可以以映样品真实情况。由表2可知，东北地区样品（R7和R8）的Chao 1指数和Simpson指数均低于其他样品，说明该地区红腐乳样品中细菌群落丰度和群落多样性较低，与其他地区存在一定差异。 表2 红腐乳样品中细菌α多样性
Table 2 α-Diversity of bacteria present in red sufu samples  样品编号序列数OTU数Shannon指数有效Chao 1指数Simpson指数Coverage指数R1 23 808 2 896 8.01 6 345.24 0.96 0.90 R2 20 472 2 462 7.68 5 378.67 0.96 0.90 R3 21 837 2 351 7.53 4 790.72 0.95 0.91 R4 24 707 2 614 7.12 6 499.29 0.93 0.90 R5 28 559 4 601 8.45 4 627.54 0.95 0.98 R6 27 197 4 474 8.74 4 506.95 0.97 0.98 R7 26 649 1 524 5.80 2 785.58 0.88 0.96 R8 21 784 1 270 5.21 2 320.00 0.82 0.96 R9 28 167 1 839 6.69 3 013.04 0.95 0.95
采用对测序序列进行随机抽样的方法，以抽到的序列数与它们所能代表OTU数目构建曲线，即稀释性曲线。当曲线趋于平坦时，说明测序数据量合理，更多的数据量对发现新OTU的边际贡献很小；以之则表明继续测序还可能产生较多新的OTU。因此，通过作稀释性曲线，可以得出样品的测序深度情况。由图1可知，多数样本在序列数在接近10 000时，OTU已经接近饱和，再增加测序量对于发现新OTU的边际贡献很小。说明测序数据量合理，所含菌种数较多。现有测序量已经可以以映样品中细菌丰度信息。  图1 各样品中细菌稀释曲线
Fig. 1 Rarefaction curves of bacteria in samples
2.2 细菌群落结构差异性分析 图2 门水平下红腐乳样品中细菌群落结构分布
Fig. 2 Distribution of bacterial communities in red sufu samples at phylum level
由图2可知，红腐乳中主要细菌优势门分别为：平均相对丰度为73.93%的厚壁菌门（Firmicutes）、15.12%的放线菌门（Actinobacteria）以及9.63%的变形菌门（Proteobacteria）。不同地区红腐乳中细菌主要优势菌群基本相同，但是相对丰度差异较大，即样品R2和R8中放线菌门的相对丰度分别为42.21%和44.96%，而样品R9的放线菌门相对丰度只有2.34%；R1中变形菌门的相对丰度是所有红腐乳中最高的，达到31.63%。说明在门水平下，不同地区红腐乳中细菌群落结构组成相同，仅组成比例存在一定差异。其中厚壁菌门是不同地区红腐乳中的绝对优势菌群。 在属水平下，9 个红腐乳样品共检测出296 个细菌属，图3显示相对丰度前12的物种。其中共有的菌属分别为：棒状杆菌属（Corynebacterium）、四联球菌属（Tetragenococcus）、明串珠菌属（Leuconostoc）、Rummeliibacillus、乳球菌属（Lactococcus）、魏斯氏菌属（Weissella）、链球菌属（Streptococcus）、丛毛单胞菌属（Comamonas）、芽孢杆菌属（Bacillus）以及Trabulsiella。  图3 属水平下红腐乳样品中细菌群落结构分布
Fig. 3 Distribution of bacterial communities in red sufu samples at genus level
从群落结构分布图可以发现不同样品中优势菌在种类及相对分布上有较大的差异。同一地区的样品间细菌群落组成相似：R3、R5和R6样品中优势属为明串珠菌属（29.45%、36.86%、30.72%）、乳球菌属（9.59%、12.55%、20.72%）以及四联球菌属（8.75%、11.36%、12.05%）；R4样品在丰度上与上述样品有所差异，分别对应明串珠菌属（13.33%）、乳球菌属（10.81%）、四联球菌属（47.11%）；而同样来自华东地区的样品R1，则展现了不同的细菌群落组成，相对丰度为22.83%的丛毛单胞菌属是该样品的优势菌属，这个属在其他样品中的相对丰度均小于2%。丛毛单胞菌属在1985年由De Vos等[26]建立，该属中的菌株更多地与环境污染物的降解有关[27-28]，在食品中鲜有报道。因此在食品发酵中对于该菌属的认识值得深入研究。 此外，样品R5、R6中的魏斯氏菌属（5.55%、4.03%）是相比于其他样品所特有的菌属。芽孢杆菌属和盐厌氧菌属（Halanaerobium）主要出现在R7、R8和R9样品中，芽孢杆菌属为R7和R9的优势菌，相对丰度分别为58.08%和52.24%。华东地区的优势菌属明串珠菌属在R7、R8和R9中相对丰度较低，均低于5%。R8中棒状杆菌属（44.08%）是该样品的优势菌属。相比于其他样品，R8中梭菌属（Clostridium）相对丰度最高，达到7.12%。梭菌属是能形成芽孢、厌氧生长的革兰氏染色阳性大杆菌，可以分解糖、蛋白质，通常可以产生混合的有机酸和醇类，不还原硫酸盐。它们的代谢物极富多样性，广泛分布在环境中，许多种可产生外毒素危害人类和牲畜生命健康[29]。 2.3 细菌β多样性分析根据OTU的丰度数据，热图可在属水平上将不同的OTU分块聚类，并根据颜色梯度以映不同样品中细菌群落相似性、差异性及物种聚类关系。利用R语言绘制热图，可以直观地显示不同地区样品中微生物的差异。由图4可知，上侧为样品聚类树，左侧为OTU聚类树，不同地区的9 种红腐乳可分为3 个聚类：R7、R8和R9聚为一类，R3、R4、R5和R6聚为一类，R1和R2聚为一类，说明不同样品间细菌群落存在一定的差异性。  图4 样品中不同属组成和相对丰度的热图
Fig. 4 Heatmap of composition and relative abundance of different genera
 图5 不同地区腐乳样品中细菌群落主成分分析
Fig. 5 Principal coordinate analysis of bacterial communities in red sufu samples
基于UniFrac距离的加权主坐标法对9 种红腐乳细菌群落结构进行分析。如图5所示，细菌群落第1主成分（PC1）贡献率为52.36%，第2主成分（PC2）贡献率为22.77%。两主成分之和大于70%，表明2 个成分较好地代表了样品中的细菌群落信息。从图5可知，不同腐乳样本中的菌群与地理位置存在一定关系，同一地区的红腐乳样品具有一定的聚类趋势，这与热图样品聚类结果相似：R3、R4、R5和R6聚为一类，均来源于华东地区；样品R7（东北）和R9（华北）聚为一类；R1（华东）、R2（华北）、R8（东北）3 个样品较为离散，各自独立。 3 讨 论腐乳是中国独具特色的大豆发酵食品之一，大豆经微生物发酵后所制成的一种干酪型产品[30]。通常根据其色泽风味可分为白腐乳、红腐乳、青（臭）腐乳三大类[31]。红腐乳添加红曲浸润数月后酿制而成，故表面呈自然红色。腐乳风味的形成主要取决于微生物的发酵，腐乳中微生物群落特征的差异是不同地区风味不同的重要原因之一。有学者对红腐乳中的红曲霉进行遗传学比较，发现不同地区的红曲霉存在一定的遗传差异性[32]。由于腐乳特殊的生产工艺以及未有加热杀菌的过程，使得腐乳中微生物种类和数量受到广大学者的关注。因此对腐乳中微生物菌群的分析研究有助于为产品质量控制、挖掘潜在功能性微生物及风味改良提供理论支持。目前对于腐乳中微生物菌群的研究主要分为培养法和非培养法两种。利用传统微生物分离培养法，红腐乳中主要的优势菌种为蜡状芽孢杆菌和Virgibacillus senegalensis[14]；植物乳杆菌、嗜酸乳杆菌、乳杆菌Akporo7和嗜盐乳杆菌[33]也被发现存在于腐乳中。周婷婷等[34]分析了不同品牌白腐乳样品的蜡样芽孢杆菌总数，其中有2 个样品的蜡样芽孢杆菌数达到106 CFU/g。Han等[31]对3 种类型腐乳中的不同微生物进行了鉴定和计数，发现腐乳中致病菌有蜡样芽孢杆菌、产气荚膜梭菌。值得注意的是，由于样品中存在大量的芽孢菌，干扰了金黄色葡萄球菌在选择性平板上的生长，样品中未能检测到金黄色葡萄球菌，但在白腐乳和臭腐乳样品中检测到金黄色葡萄球菌肠毒素A。可见，运用纯培养法对于腐乳中微生物菌群的认识有一定的局限性。 综上所述，采用PDCA循环护理模式对PICC置管患者进行护理，护理质量得到明显的升高，有效地预防导管相关并发症的发生，降低并发症发生率，利于患者治疗的顺利开展，促进治疗效果升高，让患者尽早康复出院，提高患者对护理的满意程度。 利用PCR-变性梯度凝胶电泳技术，比较不同品牌腐乳细菌多样性，发现乳酸杆菌属是腐乳中细菌的优势菌群[35]。耿予欢等[36]建立了不依赖纯培养的肠杆菌基因间重复一致序列-PCR技术，监测了腐乳发酵过程中微生物群落结构变化。上述两种方法虽然弥补了传统纯培养法的缺陷，但是均依赖于电泳技术，结果重现性相对较低。 高通量测序技术不仅能够快速地分析复杂微生物群落多样性，还能检测到低存在率以及不可培养的微生物。刘亚栋[37]利用高通量测序分析菌属发现同一品牌的红、白和青腐乳之间细菌组成具有较大的差异性。Huang Xiaoning等[38]利用Illumina MiSeq高通量测序平台发现肠杆菌属、不动杆菌属、乳球菌属是红腐乳发酵生产过程中的主要优势菌属。 本研究利用高通量测序技术对不同地区红腐乳中细菌多样性进行分析。从属水平上菌群结构分布可知，不同地区的腐乳中细菌的多样性存在差异。华东地区的腐乳主要优势菌属为明串珠菌属、乳球菌属以及四联球菌属；芽孢杆菌属则是东北地区红腐乳的主要优势菌属；此外，魏斯氏菌属、梭菌属则分别是华东地区（江苏省）、东北地区（黑龙江省）所特有的菌属。 对9 种红腐乳细菌群落结构进行主成分分析，发现同一地区的红腐乳样品中菌群存在一定的聚类关系。目前，我国腐乳生产发酵依据豆腐坯培菌的菌种不同，可分为毛霉型、根霉型和细菌型[39]，如R8样品是采用藤黄微球菌（Micrococcus luteus）进行发酵的细菌型腐乳，使得菌群结构不同于其他东北地区的样品，两个东北地区样品菌群结构距离较远。由于腐乳的生产过程多为开放式生产，且多为传统的手工操作，生产车间空气、主要接触面以及生产过程中的原辅料中的微生物均会对最终的产品产生影响[7]，使得每个厂家的腐乳中微生物菌群具有一定的独特性，这可能是R1样品与其他华东地区样品离散的原因；而R2样品中还添加了玫瑰花作为辅料，玫瑰花中的微生物在后发酵过程中影响了该样品的菌群结构。 这次研究中各项与所选取患者相关的数据都导入到SPSS18.0中进行处理,使用平均值±标准差的形式来表示计量资料,t作为检测指标,计数资料则使用百分数的形式来表示,卡方作为检测指标,其中P<0.05时证明两组数据之间具有明显的差异存在。 4 结 论本研究通过高通量测序技术对不同地区红腐乳中细菌多样性进行分析，对其中优势菌群组成进行了初步分析。结果发现，本次测序深度有效地覆盖了样品中微生物种类，不同地区红腐乳样品中细菌群落结构具有丰富多样性。在门水平上，厚壁菌门、放线菌门和变形菌门为红腐乳中主要的门；而在属水平上，不同地区红腐乳样品中优势菌在种类及相对分布上有较大的差异，明串珠菌属、乳球菌属以及四联球菌属为华北地区样品中主要的优势菌属，而东北地区腐乳样品中的芽孢杆菌属和盐厌氧菌属为主要的优势菌属。 参考文献： [1] STEINKRAUS K H. Chinese sufu handbook of indigenous fermented foods[M]. New York: Marcel Dekker, 1996. [2] YIN L J, LI L T, LI Z G, et al. Changes in isoflavone contents and composition of sufu (fermented tofu) during manufacturing[J]. Food Chemistry, 2004, 87(4): 587-592. DOI:10.3136/fstr.10.324. [3] REN H, LIU H, ENDO H, et al. Anti-mutagenic and anti-oxidative activities found in Chinese traditional soybean fermented products furu[J]. Food Chemistry, 2006, 95(1): 71-76. DOI:10.1016/j.foodchem.2004.12.019. [4] CAI R C, LI L, MEI Y, et al. Changes in bioactive compounds and their relationship to antioxidant activity in white sufu during manufacturing[J]. International Journal of Food Science & Technology,2016, 51(7): 1721-1730. DOI:10.1111/ijfs.13149. [5] MA Y, CHENG Y, YIN L, et al. Effects of processing and NaCl on angiotensin I-converting enzyme inhibitory activity and γ-aminobutyric acid content during sufu manufacturing[J]. Food & Bioprocess Technology, 2013, 6(7): 1782-1789. DOI:10.1007/s11947-012-0852-3. [6] CHEN J, QUAN M H, CHENG Y Q, et al. Acetylcholinesterase inhibitory activity of Chinese sufu (fermented tofu) ethanolextract[J]. Food Chemistry, 2012, 134(3): 1263-1266. DOI:10.1016/j.foodchem.2012.02.141. [7] 喻世哲, 范熠, 韩北忠, 等. 低盐腐乳生产过程及加工环境中的微生物分析[J]. 中国酿造, 2017, 36(12): 18-22. DOI:10.11882/j.issn.0254-5071.2017.12.005. [8] TEN B B, DAMINK C, JOOSTEN H M, et al. Occurrence and formation of biologically active amines in foods[J]. International Journal of Food Microbiology, 1990, 11(1): 73-84. DOI:10.1016/0168-1605(90)90040-C. [9] KUNG H F, LEE Y H, CHANG S C, et al. Histamine contents and histamine-forming bacteria in sufu products in Taiwan[J]. Food Control, 2007, 18(5): 381-386. DOI:10.1016/j.foodcont.2006.02.012. [10] PARK Y K, LEE J H, MAH J H. Occurrence and reduction of biogenic amines in traditional Asian fermented soybean foods: a review[J]. Food Chemistry, 2019, 278: 1-9. DOI:10.1016/j.foodchem.2018.11.045. [11] 周佳, 罗婷婷, 成长玉, 等. 腐乳生产及储存过程中氨基甲酸乙酯含量控制策略[J]. 食品安全质量检测学报, 2018, 9(16): 4280-4285. [12] HAN B, SESENNA B, BEUMER R, et al. Behaviour of Staphylococcus aureus during sufu production at laboratory scale[J]. Food Control, 2005, 16(3): 243-247. DOI:10.1016/j.foodcont.2004.03.001. [13] 文安燕, 张帅, 曾海英, 等. 贵州主要发酵豆制品加工中丙烯酰胺形成动态[J]. 中国酿造, 2017, 36(3): 85-89. DOI:10.11882/j.issn.0254-5071.2017.03.018. [14] 廖新浴, 陈信贤, 刘东红, 等. 不同种类腐乳产品中优势菌种的分离与鉴定[J]. 食品安全质量检测学报, 2018, 9(14): 3755-3759. . [15] TAGO K, ITOH H, KIKUCHI Y, et al. A fine-scale phylogenetic analysis of free-living Burkholderia species in sugarcane fi eld soil[J].Microbes and Environments, 2014, 29(4): 434-437. DOI:10.1264/jsme2.ME14122. [16] KRYACHKO Y, DONG X, SENSEN C W, et al. Compositions of microbial communities associated with oil and water in a mesothermic oil field[J]. Antonie Van Leeuwenhoek, 2012, 101(3): 493-506.DOI:10.1007/s10482-011-9658-y. [17] YANG L, YANG H L, TU Z C, et al. High-throughput sequencing of microbial community diversity and dynamics during douchi fermentation[J]. PLoS ONE, 2016, 11(12): e168166. DOI:10.1371/journal.pone.0168166. [18] XU Y, YU Y, TIAN Y, et al. Analysis of Beijing Douzhir microbiota by high-throughput sequencing and isolation of acidogenic, starchflocculating strains[J]. Frontiers in Microbiology, 2018, 9: 1142.DOI:10.3389/fmicb.2018.01142. [19] RUI W J, CHAO Z J, PU S, et al. Bacterial community involved in traditional fermented soybean paste dajiang made in northeast China[J]. Annals of Microbiology, 2013, 63(4): 1417-1421.DOI:10.1007/s13213-013-0604-2. [20] ZHU Y, CAO Y, YANG M, et al. Bacterial diversity and community in Qula from the Qinghai-Tibetan Plateau in China[J]. PeerJ, 2018, 6:e6044. DOI:10.7717/peerj.6044/table-1. [21] 杨琦玥, 黄勇, 陈亚冰, 等. 16S rRNA高通量测序技术筛选牦牛瘤胃细菌基因组DNA提取方法及菌群结构[J]. 中国农业科学, 2017,50(5): 932-941. DOI:10.3864/j.issn.0578-1752.2017.05.016. [22] 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. [23] EDGAR R C. Search and clustering orders of magnitude faster than BLAST[J]. Bioinformatics, 2010, 26(19): 2460. DOI:10.1093/bioinformatics/btq461. [24] WANG Q, GARRITY G M, TIEDJE J M, et al. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J]. Applied and Environmental Microbiology,2007, 73(16): 5261-5267. DOI:10.1128/aem.00062-07. [25] LOZUPONE C, KNIGHT R. UniFrac: a new phylogenetic method for comparing microbial communities[J]. Applied and Environmental Microbiology, 2005, 71(12): 8228-8235. DOI:10.1128/AEM.71.12.8228-8235.2005. [26] DE VOS P, KERSTERS K, FALSEN E, et al. Comamonas Davis and Park 1962 gen. nov., nom. rev. emend., and Comamonas terrigena Hugh 1962 sp. nov., nom. rev.[J]. International Journal of Systematic Bacteriology, 1985, 35(4): 443-453. DOI:10.1099/00207713-35-4-443. [27] 陈思雨, 王丽娜, 杨惠琪, 等. 白洋淀抗重金属细菌多样性及其与生态因子相关性[J]. 河北大学学报(自然科学版), 2018, 38(6): 664-672.DOI:10.3969/j.issn.1000-1565.2018.06.016. [28] 赵燕, 薛林贵, 李琳, 等. 丛毛单胞菌在环境污染物降解方面的研究进展[J]. 微生物学通报, 2012, 39(10): 1471-1478. [29] 王祖忠, 王朝阳, 张迪骏, 等. 宁波沿海陆源排污口弓形杆菌属(Arobacter sp.)和梭菌属(Clostridium sp.)的分布特点[J]. 海洋与湖沼, 2016, 47(4): 862-868. DOI:10.11693/hyhz20160400081. [30] HAN B Z, BEUMER R R, ROMBOUTS F M, et al. Microbiological safety and quality of commercial sufu: a Chinese fermented soybean food[J]. Food Control, 2001, 12(8): 541-547. DOI:10.1016/S0956-7135(01)00064-0. [31] HAN B Z, ROMBOUTS F M, NOUT M J. A Chinese fermented soybean food[J]. International Journal of Food Microbiology, 2001,65(1/2): 1-10. DOI:10.1016/S0168-1605(00)00523-7. [32] 代文婷, 吴宏, 夏芳, 等. 不同地区腐乳中10 种红曲霉菌的遗传差异分析[J]. 上海农业学报, 2017, 33(3): 19-24. DOI:10.15955/j.issn1000-3924.2017.03.04. [33] 张志超, 罗剑鸣, 殷光玲. 商业化腐乳中乳酸菌的分离与鉴定[J].食品安全导刊, 2018(18): 154-156. [34] 周婷婷, 易敏英, 相大鹏, 等. 白腐乳及其生产过程中蜡样芽孢杆菌的污染及防治[J]. 中国酿造, 2014, 33(5): 77-80. DOI:10.11882/j.issn.0254-5071.2014.05.018. [35] 陈颖慧. PCR-DGGE分析不同品牌腐乳中细菌的多样性[J]. 中国调味品, 2017, 42(7): 29-32. DOI:10.3969/j.issn.1000-9973.2017.07.007. [36] 耿予欢, 李国基, 邹家兴. 腐乳发酵过程中微生物群落结构的ERICPCR指纹图谱分析[J]. 华南理工大学学报(自然科学版), 2010, 38(8):126-130. DOI:10.3969/j.issn.1000-565X.2010.08.024. [37] 刘亚栋. 利用16S rDNA测序的方法鉴定腐乳中微生物的种类多样性[D]. 济南: 山东师范大学, 2017. [38] HUANG X N, YU S Z, HAN B Z, et al. Bacterial community succession and metabolite changes during sufu fermentation[J]. LWT,2018, 97: 537-545. DOI:10.1016/j.lwt.2018.07.041. [39] 王润芝, 王庆周. 克东腐乳前发酵微生物的研究: 1.前发酵菌的分离鉴定[J]. 中国调味品, 1986(1): 12-14.
High-throughput Sequencing Analysis of Bacterial Diversity in Red Sufu from Different Regions XU Qiong1, LIU Yang1,*, QU Qinfeng1, DOU Tonghai2, CHEN Yufei1, ZHANG Nana1, ZHAO Lei2,ZHONG Jiang2, WENG Shiyu1, YANG Jielin3, ZHAO Guoping2
(1. Shanghai Institute of Quality Inspection and Technical Research, Shanghai 200233, China;2. School of Life Sciences, Fudan University, Shanghai 200433, China;3. Technical Center for Animal, Plant and Food Inspection and Quarantine, Shanghai Customs, Shanghai 200135, China) Abstract: To understand the bacterial fl ora structure and microbial safety of red sufu as a traditional Chinese fermented food, the V1-V3 region of the bacterial 16S rDNA from red sufus produced in the east, north and northeast of the country was sequenced by high-throughput sequencing and the differences in the structure and diversity of microbial community was analyzed among these samples. The results showed that the sequencing depth effectively covered the microbial species in the samples, and the bacterial communities in them were rich and diverse. At the phylum level, Firmicutes, Actinobacteria and Proteobacteria were dominant in all samples. A total of 296 bacterial genera were detected, with Corynebacterium,Tetragenococcus, Leuconostoc, Rummeliibacillus, Lactococcus, Weissella, Streptococcus, Comamonas, Bacillus and Trabulsiella being common to these samples. At the same time, clustering analysis demonstrated that there was a relationship between the bacterial community structure and the geographical origin of red sufu. The main dominant bacteria in the samples from North China were Candida, Lactobacillus and Tetraphylococcus, while those in the red sufu samples from Northeast China were Bacillus and Halanaerobe. This study gives further insights into the structure and diversity of microbial community in red sufu, providing a theoretical basis for production of traditional fermented foods with higher quality and safety. Keywords: high-throughput sequencing; red sufu; bacterial diversity
收稿日期：2019-04-08 基金项目：国家市场监督管理总局科技计划项目（2017QK167）；上海市市场监督管理局青年科技启明星计划项目（QMX2018012）；上海市技术标准专项（16DZ0500101） 第一作者简介：徐琼（1987—）（ORCID: 0000-0001-5323-0756），男，高级工程师，硕士，研究方向为食品微生物。E-mail: xuqiong@sqi.org.cn*通信作者简介：刘洋（1980—）（ORCID: 0000-0002-5665-034X），女，正高级工程师，博士，研究方向为食品微生物。E-mail: liuyang@sqi.org.cnDOI:10.7506/spkx1002-6630-20190408-087 中图分类号：Q939.97 文献标志码：A 文章编号：1002-6630（2020）10-0110-07 引文格式： 徐琼, 刘洋, 曲勤凤, 等. 高通量测序分析不同地区红腐乳细菌多样性[J]. 食品科学, 2020, 41(10): 110-116. DOI:10.7506/spkx1002-6630-20190408-087. http://www.spkx.net.cn近些年，黑龙江垦区在低碳农业的发展方面已做出了巨大的努力并取得了显著的成效。如2011年新增造林绿化面积1.33万hm2多，41个农场纳入大小兴安岭生态功能区建设，垦区城镇绿化覆盖率达到37%，万元GDP综合能耗降到0.936t标准煤，比上年减少了2.5%。但农用化肥、农药、农膜的使用量却呈现出不断增长的态势，如表2所示。 XU Qiong, LIU Yang, QU Qinfeng, et al. High-throughput sequencing analysis of bacterial diversity in red sufu from different regions[J]. Food Science, 2020, 41(10): 110-116. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20190408-087. http://www.spkx.net.cn
| 
|Archiver|手机版|小黑屋|www.openhelp100.com
( 冀ICP备19026749号-1 )
发表于 2021-2-10 10:45:04
