|
基于高通量测序技术分析蜂粮微生物多样性基于高通量测序技术分析蜂粮微生物多样性 刘玥佳1,2，韩业君2，彭文君1，牛庆生3，方小明1，赵亚周1，田文礼1,* （1.中国农业科学院蜜蜂研究所，北京 100093；2.中国科学院过程工程研究所，北京 100190；3.吉林省养蜂科学研究所，吉林 吉林 132013） 摘 要：为分析蜂粮发酵过程中微生物多样性和菌群动态变化，探究不同发酵时间对蜂粮微生物菌群组成的影响，以发酵0.5、4、7、12 d的卡尼鄂拉蜂（Apis mellifera carnica）蜂粮为研究对象，采用高通量测序技术分析不同发酵时间蜂粮微生物的16S rDNA序列，比较各组菌群构成和丰度信息，通过α多样性分析和蜂粮菌群组成分析，考察蜂粮中的优势菌群。在4 个蜂粮样品中共鉴定出27 个门、112 个属的细菌，其中4 个样品中有25 个共有菌属。结果表明：测序深度有效地覆盖了微生物种类，随着发酵时间的延长，蜂粮微生物群落的物种丰富度和多样性先升高后降低。基于门水平，生氧光细菌门（Oxyphotobacteria）和变形菌门（Proteobacteria）为优势菌群；基于属水平，unidentified Oxyphotobacteria为优势微生物。发酵过程中，unidentified Oxyphotobacteria表现出先增加后减少的趋势。 关键词： 蜂粮；微生物群落多样性；高通量测序 蜂粮是蜜蜂将蜂花粉、蜂蜜、蜜蜂分泌物等充分混合，经特定微生物充分发酵后，贮存于巢房中供蜜蜂食用的物质。其营养丰富、风味独特，是一种极具开发潜力的发酵食品[1]。由于从蜂群中直接获取量少且非常繁琐，理想的方式是实现离群、工业化生产，而制约工业化生产的因素众多，其中关键是发酵微生物的确定。对于研究食品中微生物的群落组成及变化，常用的传统培养分离方法费时、费力，且无法真正解析微生物的组成、丰度及其变化情况。高通量测序技术可省去培养过程，直接从样品中提取微生物基因组DNA后进行测序分析，结合生物信息学工具，能够检测出传统方法无法分离培养的微生物[2-3]。蜂粮发酵过程中的菌群一直处于动态变化中，因此，利用高通量测序技术研究和分析不同发酵周期的蜂粮菌群构成及其演替规律，具有重要的理论意义和潜在的实用价值。 目前，在食品领域中已有大量将高通量测序技术用于微生物群落分析的研究[4-6]，例如对干酪、韩国泡菜[7]等发酵食品中微生物群落结构的分析，揭示了发酵食品酿造过程中微生物群落结构、丰度变化及其与代谢产物之间的联系。蜂粮中菌群丰富，包括乳酸菌、芽孢杆菌、酵母菌、霉菌等，但关于蜂粮的菌群结构研究还比较少。刘赛[8]从蜂粮中分离纯化出39 株菌株，并从这些菌株中筛选出了3 株产酸性能较好的乳酸菌，经鉴定分别为：Lactobacillus plantarum、L. casei ssp. casei、L. debrueckii ssp. delbrueckii。Gilliam[9]从人工采集花粉、蜜蜂采集蜂花粉以及贮存在蜂箱中1、2、6 周的（蜂粮）花粉中分离出了6 个属，113 种酵母菌。此外，花粉中大部分的酵母菌种在蜂粮中并未发现，这说明蜂粮酿制过程中酵母菌群发生了变化。张言周等[10]从蜂粮样品中分离出产芽孢的细菌，经鉴定分别隶属于Lysinibacillus fusiformis、Bacillus tequilensis、Fictibacillus nanhaiensis、B. aerophilus、B. invictus、B. sonorensis、L. xylanilyticus、B. methylotrophicus等芽孢杆菌；崔学沛[11]从意蜂蜂粮中筛选得到了B. amyloliquefaciens和B. licheniformis。Gilliam等[12]从蜂粮样品中分离出了曲霉菌、青霉菌和毛霉菌等霉菌。苏松坤[13]从不同酿制时期的蜂粮中亦分离出了曲霉菌、青霉菌和毛霉菌等霉菌，由于花粉源植物和地区的不同，其含量存在一定差异。 目前，通常把蜂粮的发酵过程分为4 个阶段：1）异源微生物生长阶段；2）厌氧发酵起始阶段；3）酸化阶段；4）稳定阶段[14]。开始的0.5 d，花粉中营养物质丰富，乳酸菌和酵母菌大量繁殖。约在发酵开始7 d以后，基质内pH值达到4.0～4.5，高浓度的乳酸抑制乳酸菌和酵母菌等的生长，某些乳酸菌和酵母菌开始消失，蜂粮酿制完成[15]。乳酸菌在发酵第0.5天时丰度最高，发酵第4天时丰度下降，7 d以后乳酸菌含量下降明显。 卡尼鄂拉蜂（Apis mellifera carnica）采蜜能力强，范围广，故本研究以卡尼鄂拉蜂的蜂粮为研究对象，采用Illunina MiSeq高通量测序平台对细菌的16S rDNA V4-V5区进行测序，分析蜂粮的微生物群落多样性及各类微生物的构成及变化，明确蜂粮发酵过程中的优势，为人工蜂粮的生产提供理论依据，以期更好地提高蜂粮品质，推动其产业化发展。 1 材料与方法1.1 材料与试剂实验样品采集于吉林省养蜂科学研究所。采集方法：首先，取样的第1天，蜂箱中只保留一张大幼虫脾外，其他脾全部抽出，加一张育过5～6 代蜜蜂的空脾（事先用75%乙醇溶液喷洒消毒）；用隔王栅将蜂王限制在子脾上，使工蜂给空脾装花粉。5 h后将花粉脾抽出，均匀选取位于巢脾不同位置的蜂粮样品收集于离心管（已灭菌）中，样品采集完后置于-80 ℃冰箱冻存。10 h后取下多用器盖板，第2天早上再盖回。同时在寄放花粉脾的蜂群内加一张大花粉脾供蜂群食用[16]。按照上述方法依次取得发酵0.5、4、7、12 d的蜂粮样品。 首先，优势特色产业快速发展，奠定了农村一、二、三产业融合发展的良好基础。近年来，各级政府大力培植农业优势主导产业，逐步形成了优质专用小麦、优质稻米、食用菌、蔬菜、奶牛、大闸蟹、冬枣等一批专业化生产基地及优势产业带。 PowerSoil DNA Isolation Kit （Cat# 12888-50）试剂盒 上海斯信生物科技有限公司；胶回收试剂盒德国Qiagen公司。 卡夫卡的婚恋显然并非一般意义的，他更像是在找寻一位依赖对象，找寻自己灵魂的出口。某种程度上，卡夫卡试图通过婚姻来缓和自己同世界的关系，以此来完成对自我的拯救。但很显然，在这几段婚恋关系中，无论是从过程还是结果上来说，卡夫卡“在精神上显然没有能力结婚”。[4]461-501 1.2 仪器与设备ME104E精密天平 梅特勒-托利多国际贸易（上海）有限公司；Haier磁力搅拌器 常州迈科诺仪器有限公司。 1.3 方法1.3.1 蜂粮样品总DNA提取 由于每个巢房内可获取的蜂粮样品较少，将若干个巢房的蜂粮样品混匀进行DNA提取以提高DNA浓度，同时使蜂粮样品更具代表性。参照DNeasy PowerSoil Kit（Cat# 12888-50）试剂盒说明，分别提取不同发酵时间蜂粮细菌基因组DNA。所得DNA采用0.8%琼脂糖凝胶电泳检测DNA完整性。 本文结合仿真工具Matlab，实现高阶累积量算法识别6种通信调制信号。图2为识别2ASK、4ASK、2FSK、4FSK、2PSK、4PSK信号的识别率仿真图。图中每一种调制信号识别率曲线都用不同形状表示。 1.3.2 蜂粮16S rRNA基因V4-V5区的扩增 利用琼脂糖凝胶电泳对基因组DNA进行初步检测，即取适量的样本DNA于离心管中，使用无菌水稀释样本至1 ng/μL。以稀释后的基因组DNA为模板进行聚合酶链式以应（polymerase chain reaction，PCR）测定。 重金属Cd污染已成为贵州喀斯特高原环境下发展无公害农业土地利用的一个重要障碍。环境中Cd超标对人体的骨骼和肾脏造成不利影响，严重时能导致骨痛病，人体过量Cd的摄入，会产生拮抗作用，阻碍微量元素如锌的吸收[6]。Cd在土壤固液相之间的分配主要受吸附解吸反应影响，且在一定程度上决定了各相态Cd浓度及其生物有效性[7-10]。鉴于此，采用批量平衡法，探究镉在贵州不同地区农田土壤中的吸附解吸规律，以期为镉在贵州土壤中的迁移规律和生态治理提供科学依据。 1.4 蜂粮高通量测序和数据分析将PCR扩增产物进行琼脂糖凝胶电泳检测；切下目的条带并用胶回收试剂盒回收。委托北京诺禾致源生物信息科技有限公司基于Illumina MiSeq技术测序平台构建文库，进行高通量双末端测序。 得到下机数据后利用Barcode序列和PCR扩增引物序列进行拆分，截去Barcode和引物序列后，使用FLASH[17]对每个样本的reads进行拼接，得到的原始Tags数据需要经过严格的过滤处理[18]得到高质量的Tags数据。参照QIIME（Version 1.9.1）[19]的Tags质量控制流程，得到最终的有效数据。 利用UPARSE软件（UPARSE v7.0.1001）[20]对所有样本的全部 Effective Tags进行聚类，默认以97%的一致性将序列聚类成为可操作分类单元（operational taxonomic units，OTUs），筛选出现频率最高的OTUs序列作为代表序列。用Mothur方法与SILVA132的SSU rRNA方法将筛选到的OTUs序列进行物种注释并根据其分类学信息在各分类水平上进行分类统计，分析各个样本的群落组成[21-22]。 使用QIIME软件（Version 1.9.1）计算Observed-OTUs、Chao1指数、Shannon指数、Simpson指数、ACE指数、Coverage指数，使用R软件（Version 2.15.3）绘制稀释曲线，并进行α多样性指数组间差异分析。 2 结果与分析2.1 蜂粮样品总DNA的提取提取蜂粮基因组总DNA，用琼脂糖凝胶电泳检测其DNA完整性，结果如图1所示，可进行后续测序。  图1 蜂粮细菌基因组DNA提取
Fig. 1 Agarose gel electrophoresis of bacterial genomic DNA samples extracted from bee bread
2.2 蜂粮发酵过程中细菌群落的α多样性分析利用稀释曲线评价蜂粮菌群的测序深度和物种丰富度。从图2可以看出，随着测序深度增加，稀释曲线基本达到饱和，说明样本的测序数据量足以以应蜂粮中的微生物多样性。 每个样品的Coverage指数、Chao1指数和Shannon指数分别用于评估物种的测序深度、丰度和多样性。Coverage指数以映了样本的真实情况，Coverage指数越高，则样本中序列没有被测出的概率越低；Chao1指数是群落丰度指数，Chao1指数越大，表明样品微生物群落丰度越高；Shannon指数是群落分布多样性指数，Shannon指数越大，表明样品微生物多样性越高；OTUs数量亦可以代表样品物种的丰度[23]。由表1可知，不同发酵周期的蜂粮中微生物Coverage指数都大于0.99，说明样品文库中序列基本上都被测出，即样本测序结果可以以映样品的真实情况。随着发酵时间的延长，Chao1指数、ACE指数和Shannon指数先升高后降低，发酵中期，物种丰富度达到最高。由图3可知，4 个样品共得到429 个OTUs，不同发酵周期共有OTUs数为99，各个时期独有的OTUs数相对较少。 治疗及预后 轻度肾损伤患者监测肾功能指标，中度肾损伤患者按照损伤程度调整替诺福韦的剂量，重度肾损伤 （包括FS）患者将替诺福韦更换为阿巴卡韦。电解质紊乱补充电解质，重度代谢性酸中毒给予碳酸氢钠中和等对症治疗，严重肾损伤给予人工肾治疗。发生肾损伤的98例患者中，82例 （84%）经过治疗1～2个月后电解质紊乱、尿糖、尿蛋白和高氯性代谢性酸中毒好转；16例2个月后电解质、尿糖、尿蛋白、肌酐等无明显缓解，其中4例病情出现加重，经过1年的跟踪和随访，6例 （6%）无明显缓解，大部分与处理不及时有关。  图2 蜂粮测序样品稀释曲线
Fig. 2 Rarefaction curves of bee bread
表1 α多样性指数
Table 1 α Diversity indices of different samples  蜂粮发酵时间/d OTUs数/个Shannon指数Chao1指数Simpson指数ACE指数Coverage指数0.5 205 2.725 254.214 0.702 247.385 0.999 4 279 2.749 294.122 0.667 300.518 0.999 7 283 3.251 626.462 0.749 388.924 0.998 12 215 2.524 225.973 0.651 230.902 0.999
 图3 OTUs分布Venn图Fig. 3 Venn diagram showing the distribution of OTUs
2.3 蜂粮菌群多样性2.3.1 基于门分类水平的蜂粮菌群多样性分析  图4 蜂粮菌群基于门分类水平的分布
Fig. 4 Bacterial distribution at the phylum level
如图4所示，高通量测序分析表明蜂粮菌群主要由软壁菌门（Tenericutes）、广古菌门（Euryarchaeota）、（Acetothermia）、奇古菌门（Thaumarchaeota）、放线菌门（Actinobacteria）、厚壁菌门（Firmicutes）、拟杆菌门（Bacteroidetes）、变形菌门（Proteobacteria）、生氧光细菌门（Oxyphotobacteria）、螺旋体门（Spirochaetes）等组成。随着发酵时间的延长，基于门水平的微生物多样性发生了变化。 不同发酵周期的蜂粮在门分类水平上，优势菌门所占比例相似，生氧光细菌和变形菌的含量最高，为优势菌群。发酵时间为0.5、4、7、12 d的蜂粮中生氧光细菌和变形菌的含量分别为：90.2%、89.0%、85.6%、92.2%；1.7%、4.2%、5.6%、1.5%。 衣原体门（C hla m y dia e）、梭杆菌门（Fusobacteria）、泉古菌门（Crenarchaeota）在发酵初期含量比较少，而在发酵4、7、12 d的蜂粮样品中均未检测出，可能是来自于蜂花粉，随着发酵的进行，加上蜜蜂的转化，这几种菌逐渐消失。Calditrichaeota、硝化螺旋菌门（Nitrospirae）在发酵0.5、4、7 d的蜂粮中均未检测出，而在发酵12 d的蜂粮中少量存在，说明其产生于蜂粮的发酵后期。 2.3.2 基于纲分类水平的蜂粮菌群多样性分析 现代翻译学中等值或者说对等理论是一重要的思想。诸多西方翻译理论家都在探讨翻译等值问题。当提及等值之时，并不是指本体论上的绝对等同，而是认识论上的等值，也就是说，在认识事物的实践中是等值的。  图5 蜂粮菌群基于纲分类水平的分布
Fig. 5 Bacterial distribution at the class level
如图5所示，高通量测序分析表明蜂粮菌群主要由unidentified Oxyphotobacteria、α-变形杆菌纲（Alphaproteobacteria）、γ-变形杆菌纲（Gammaproteobacteria）、拟杆菌纲（Bacteroidia）、梭菌纲（Clostridia）、δ-变形菌纲（Deltaproteobacteria）、杆菌纲（Bacilli）、Nitrososphaeria、unidentified Actinobacteria、酸微菌纲（Acidimicrobiia）等组成。随着发酵时间的延长，基于纲水平的微生物群落结构发生了变化。Nitrososphaeria和Acidimicrobiia在发酵第0.5天和12天时均为0，仅在发酵中期存在。 蜂粮中主要的菌群为生氧光细菌门Oxyphotobacteria，在蜂粮中占到了80%以上，但其只鉴定到了门分类水平，无法具体鉴定到门以下水平，可能是现有数据库中关于此类菌群的信息较少，所以序列没有得到高相似度的16S DNA比对结果。 2.3.3 基于目水平的蜂粮菌群多样性分析  图6 蜂粮菌群基于目分类水平的分布
Fig. 6 Bacterial distribution at the order level
如图6所示，测序结果表明蜂粮菌群主要由unidentified Oxyphotobacteria、unidentified Alphaproteobacteria、拟杆菌目（Bacteroidales）、假单胞菌目（Pseudomonadales）、梭菌目（Clostridiales）、肠杆菌目（Enterobacteriales）、unidentified Gammaproteobacteria、立克次氏体目（Rickettsiales）、互营杆菌目（Syntrophobacterales）、黄杆菌目（Flavobacteriales）等组成。这些菌均出现于整个发酵周期。 2.3.4 基于科水平的蜂粮菌群多样性分析 1996年，比尔·盖茨在意大利的佛罗伦萨大学演讲时，发现有一个男学生，高举着一块写着自己名字的白色硬纸板。  图7 蜂粮菌群基于科分类水平的分布
Fig. 7 Bacterial distribution at the family level
如图7所示，测序结果表明蜂粮菌群主要由unidentified Oxyphotobacteria、醋杆菌科（Acetobacteraceae）、unidentified Alphaproteobacteria、肠杆菌科（Enterobacteriaceae）、莫拉氏菌科（Moraxellaceae）、（Muribaculaceae）、普雷沃氏菌科（Prevotellaceae）、unidentified Rickettsiales、unidentified Gammaproteobacteria、互营杆菌科（Syntrophobacteraceae）等组成。 2.3.5 基于属水平的蜂粮菌群多样性分析  图8 蜂粮菌群基于属分类水平的分布
Fig. 8 Bacterial distribution at the genus level
如图8所示，共鉴定出112 个属：在发酵0.5 d的蜂粮中共鉴定出46 个属，4 d的蜂粮中鉴定出个69 个属，7 d的蜂粮中鉴定出77 个属，12 d的蜂粮中鉴定出47 个属，其中4 个样品中有25 个共有菌属。 基于属水平，蜂粮最优势微生物为unidentified Oxyphotobacteria，随着发酵的进行，unidentified Oxyphotobacteria受到抑制，发酵中期减少，发酵后期又增加。但unidentified Oxyphotobacteria始终是优势菌。此外，玫瑰单胞菌属（Roseomonas）、unidentified Alphaproteobacteria、不动杆菌属（Acinetobacter）、热硫还原杆菌属（Thermodesulforhabdus）、unidentified Gammaproteobacteria、短波单胞菌属（Brevundimonas）均呈先增加后减少的趋势，与unidentified Oxyphotobacteria基本一致。unidentified Alphaproteobacteria与unidentified Gammaproteobacteria在发酵第1天时，均为0，发酵中期才存在，很可能来自蜜蜂的肠道菌群。Gammaproteobacteria正是蜜蜂肠道的八大菌群；Stenotrophomonas、Sphingobacterium、Anaerobiospirillum、Odoribacter、Rhodococcus、Bacteroides、Geothermobacter、Calorithrix、Anaeromyxobacter在发酵第1天时，均为0，发酵后期才存在，说明这些菌产生于蜂粮的发酵后期；还有一些菌，如unidentified Gammaproteobacteria、念珠菌（Candidatus）、Nitrosopelagicus、SUP05 cluster、Candidatus Actinomarina、unidentified Bacteria、unidentified Rhodospirillales、unidentified Thermoplasmata、弧形菌属（Vibrio）、Puniceispirillum、Aquibacter、嗜热弯曲甲烷热杆菌（Methanothermobacter）、Limnobacter、瘤胃梭菌属（Ruminiclostridium）等在发酵初期和末期均不存在，只存在于发酵中期。 3 讨 论由于每个巢房内可挖取的蜂粮样品较少，小于0.1 g，进行DNA提取的浓度太低，无法建库进行后续测序，故参考Vigneron等[24]的实验方法，每次取样时均匀选取位于巢脾不同位置的蜂粮样品，将若干个巢房的蜂粮样品混匀，进行DNA的提取以提高DNA浓度，同时可以使蜂粮样品更具代表性。 对不同发酵时间的蜂粮中的细菌种群结构和多样性变化的研究表明，随着发酵的进行，细菌物种多样性和丰富度先升高后降低。但主要优势微生物并没有改变，始终为生氧光细菌和变形菌，生氧光细菌是蜂粮中丰度最高的菌群。unidentified Oxyphotobacteria在发酵过程中丰度一直最高，为发酵过程的主要优势菌。 近年来，科研人员基于16S rRNA序列和宏基因组分析技术，对不同地区的西方蜜蜂的肠道菌群进行了检测，结果显示有95%以上都属于8大类群，分别是Gilliamella、Gammaproteobacteria、Snodgrassella、Bartonella、Acetobacteraceae、Lactobacillus、Lactobacillus acidophilus、Bifidobacterium[25-33]。本研究中发现的unidentified Alphaproteobacteria与unidentified Gammaproteobacteria在发酵第1天时均未被检出，发酵中期才存在，很可能来自蜜蜂的肠道菌群。 在项目水资源论证报告书审查过程中，专家组详细分析了各项目区水资源开发条件、各水源供水基础设施数量及供水能力、机电井分布及供水状况，明确要求水源确定时要优先考虑地表水，合理开发地下水，禁止开采深层承压水。地表水资源相对丰富、具备供水条件和供水能力的项目区，应进一步压缩地下水开采量；项目实施必须开采地下水资源的，要实行地下水总量与水位双控制，科学论证地下水取水规模、新建和改造井数及适宜井深；严格控制生态脆弱区取用地下水，禁止在地下水超采区新增地下水开采，确保项目实施后地下水的采补平衡。 Gilliamella、Lactobacillus、Bartonella、Bifidobacterium等几种菌在蜂粮样品中均被检出。Gilliamella与Bifidobacterium均在第7天的样品中含量最高；Lactobacillus在发酵初期含量最高，而后一直呈下降趋势；Bartonella在发酵末期含量最高。Engel等[34]对西方蜜蜂肠道菌群宏基因组的分析发现：其肠道菌可能有多种功能，例如与宿主相互作用，形成生物膜，分解碳水化合物等，对病原菌的防御和免疫起着至关重要的作用。 在电力系统安全生产的过程中，员工的安全保障意识还有待提升，许多安全事故的发生大多是由于工作人员没有足够的安全生产意识。这不仅为工作人员带来了许多安全隐患，也不利于电力公司的发展。所以，我国的电力公司必须建立一个科学合理、安全性强的安全生产制度，并且加强对员工进行安全教育和安全培训。同时还要对电力系统中的生产工具进行严格的检修与管理，避免因工具有缺陷而造成的安全生产事故、 Engel等[35]通过对Bartonella进一步分析发现：Bartonella含有VB12的合成系统，可以合成VB12供给蜜蜂；此外，Bifidobacterium、Lactobacillus和Gilliamella还具有合成果胶降解酶、糖苷水解酶和多糖水解酶等功能。而果胶纤维素是花粉内壁的主要成分，果胶降解酶的合成表明这些菌群有助于蜂花粉释放营养物质。这些潜在的功能也从基因水平上以映了乳酸菌可能参与了蜜蜂体内碳水化合物的代谢。 4 结 论本实验通过对卡尼鄂拉蜂不同发酵时间的蜂粮进行高通量测序，分析和研究了蜂粮发酵过程中细菌群落的组成及其丰度变化，考察不同发酵时间对蜂粮发酵过程中菌群的影响，共27 个门、112 个属的细菌被鉴定出。基于门水平，Oxyphotobacteria和Proteobacteria是优势菌群。基于属水平，蜂粮最优势微生物为unidentified Oxyphotobacteria，随着发酵的进行，unidentified Oxyphotobacteria受到抑制，发酵中期减少，发酵后期又增加，但unidentified Oxyphotobacteria始终是优势菌。生氧光细菌与光合作用有关[36]，从而推断本实验的优势菌unidentified Oxyphotobacteria可能与植物的光合作用有关，说明这种菌有可能来自于花粉源植物。此外，一些菌如unidentified Alphaproteobacteria与unidentified Gammaproteobacteria等，在发酵第1天时，均为0，发酵中期才存在，很可能来自蜜蜂的肠道菌群。进一步验证了Mattila等[37]认为蜂粮中的菌群一部分来自于蜜蜂肠道菌群的研究结果。 无花果树叶绿油油的，蓄势待发，那架势，有点逼人。再加上侯大同的放任，无花果简直不把女主人放在眼里。平时侧一下身子也就过了，一到春天，连侯大同自己都觉着不好意思，赶紧剪了东边的枝杈，留下一条进出主房的通道。汤翠嫁过来的时候倒不觉得院子多逼仄，那时候，东边还没有盖偏房，无花果跟未发育的小孩子一样，还小，还没有完全长开。没几年，家里添了女儿，偏房也建起来，院子就显得狭小。现在汤翠不愁了，无花果即使堵着门也无所谓了，反正早晚要拆了。 蜂粮中主要的菌群为生氧光细菌门Oxyphotobacteria，在蜂粮中占到了80%以上，但其只鉴定到了门分类水平，无法具体鉴定到门以下水平，可能是现有数据库中关于此类菌群的信息较少，所以序列没有得到高相似度的16S DNA比对结果。研究结果对于蜂粮中该类微生物的分离鉴定具有重要意义。本研究为实现蜂粮工业化生产过程中发酵菌株的确定及发酵条件确定、进一步筛选优势发酵菌株提供一定的参考依据。 李白的《东武吟》改为叙写自己的经历，回顾三年志得意满的翰林生活及离开长安后凄凉的境况，失意中潜藏着一份用舍由时、行藏在我的豁达，昂扬的情感成了诗歌的主旋律。 参考文献： [1] 段传人, 暴佳芳, 石轶松, 等. 枇杷花和油菜花蜂粮中微生物的分离与鉴定[J]. 食品工业科技, 2015, 36(5): 175-180. DOI:10.13386/j.issn1002-0306.2015.05.028. [2] HANDELSMAN J, RONDON M R, BRADY S F, et al. Molecular biological access to the chemistry of unknown soil microbes: a new frontier for natural products[J]. Chemistry Biology, 1998, 5(10): 245-249. DOI:10.1016/s1074-5521(98)90108-9. [3] HANDELSMAN J. Metagenomics application of genomics touncultured in microorganisms[J]. Microbiology and Molecular Biology Reviews, 2004, 68(4): 669-685. DOI:10.1128/MMBR.68.4.669-685.2004. [4] WOLFE B E, DUTTON R J. Fermented foods as experimentally tractable microbial ecosystems[J]. Cell, 2015, 161(1): 49-55.DOI:10.1016/j.cell.2015.02.034. [5] KEEGOURLAY G, TAMINIAU B, DAUBE G, et al. Metagenomic insights into the dynamics of microbial communities in food[J].International Journal of Food Microbiology, 2015, 213: 31-39.DOI:10.1016/j.ijfoodmicro.2015.09.010. [6] VELASCO G L, WELBAUM G E, BOYER R R, et al. Changes in spinach phylloepiphytic bacteria communities following minimal processing and refrigerated storage described using pyrosequencing of 16S rRNA amplicons[J]. Journal of Applied Microbiology, 2011,110(5): 1203-1214. DOI:10.1111/j.1365-2672.2011.04969.x. [7] TAPIA A A, RAMIREZ M C E, TAMPLIN M L, et al. Highthroughput sequencing of microbial communities in Poro cheese, an artisanal Mexican cheese[J]. Food Microbiology, 2014, 44: 136-141.DOI:10.1016/j.fm.2014.05.022. [8] 刘赛. 蜂粮中乳酸菌的分离鉴定及其在蜂花粉发酵中的应用[D].福建: 福建农林大学, 2012: 184-193. [9] GILLIAM M. Microbiology of pollen and bee bread: the yeasts[J].Apidologie, 1979, 10(1): 43-53. DOI:10.1051/apido:19790106. [10] 张言周, 王爽, 李韵雅, 等. 蜂粮中芽孢杆菌的分离鉴定及其耐酸耐高糖特性的分析[J]. 食品与生物技术学报, 2017, 36(4): 410-415.DOI:10.3969/j.issn.1673-1689. [11] 崔学沛. 意大利天然成熟蜂粮菌种鉴定及人工蜂粮发酵影响因素的研究[D]. 泰安: 山东农业大学, 2015. [12] GILLIAM M, PREST D B, LORENZ B J. Microbiology of pollen and bee bread: taxonomy and enzymology of molds[J]. Apidologie, 1989,20(1): 53-68. DOI:10.1051/apido:19890106. [13] 苏松坤. 蜂粮的酿制机理和营养价值的研究[D]. 杭州: 浙江大学,2000. [14] 刘玥佳, 韩业君, 彭文君, 等. 蜂粮的营养成分和微生物发酵过程研究进展[J]. 中国酿造, 2018, 37(3): 23-27. DOI:10.11882/j.issn.0254-5071.2018.03.006. [15] 江武军, 曾志将. 蜂粮研究进展[J]. 黑龙江畜牧兽医, 2015, 1(2):118-119. DOI:10.13881/j.cnki.hljxmsy.2015.013. [16] 苏松坤, 陈盛禄, 刘华, 等. 天然蜂粮生产方法[J]. 养蜂科技,2006(6): 32-33. [17] MAGOC T, SALZBERG S L. FLASH: fast length adjustment of short reads to improve genome assemblies[J]. Bioinformatics, 2011, 27(21):2957-2963. DOI:10.1093/bioinformatics/btr507. [18] BOKULICH N A, SUBRAMANIAN S, FAITH J J, et al. Qualityfi ltering vastly improves diversity estimates from Illumina amplicon sequencing[J]. Nature Methods, 2013, 10(1): 57-60. DOI:10.1038/nmeth.2276. [19] 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. [20] HAAS B J, GEVERS D, EARL A M, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons[J]. Genome Research, 2011, 21(3): 494-504.DOI:10.1101/gr.112730.110. [21] EDGAR R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads[J]. Nature Methods, 2013, 10(10): 996-998.DOI:10.1038/NMETH.2604. [22] 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. [23] KE L Q, WANG L L, LI H B, et al. Molecular identification of lactic acid bacteria in Chinese rice wine using species-specific multiplex PCR[J]. European Food Research and Technology, 2014, 239(1): 59-65. DOI:10.1007/s00217-014-2193-0. [24] VIGNERON A, ALSOP E B, CHAMBERS B, et al. Complementary microorganisms in highly corrosive biofilms from an offshore oil production facility[J]. Apply Environment Microbiology, 2016, 82(8):2545-2554. DOI:10.1128/AEM.03842-15. [25] 徐龙龙, 吴杰, 李继莲. 蜜蜂共生菌研究进展[J]. 中国农业科技导报,2013, 15(6): 107-112. DOI:10.3969/j.issn.1008-0864.2013.06.16. [26] FOSTER D L C, CONLAN S, HOLMES E C, et al. A metagenomic survey of microbes in honey bee colony collapse disorder[J]. Science,2007, 318: 283-287. DOI:10.1126/science.1146498. [27] ENGEL P, MARTINSON V G, MORAN N A. Functional diversity within the simple gut microbiota of the honey bee[J]. PNAS, 2012,109(27): 11002-11007. DOI:10.1073/pnas.1202970109. [28] BABENDREIER D, JOLLER D, ROMEIS J, et al. Bacterial community structures in honeybee intestines and their response to two insecticidal proteins[J]. FEMS Microbiol Ecology, 2007, 59(3): 600-610. DOI:10.1111/j.1574-6941.2006.00249.x. [29] JEYAPRAKASH A, HOY M A, ALLSOPP M H. Bacterial diversity in worker adults of Apis mellifera capensis and Apis mellifera scutellata (Insecta: Hymenoptera) assessed using 16S rRNA sequences[J]. Journal of Invertebrate Pathology, 2003, 84(2): 96-103.DOI:10.1016/j.jip.2003.08.007. [30] MARTINSON V G, DANFORTH B N, MINCKLEY R L, et al.A simple and distinctive microbiota associated with honey bees and bumble bees[J]. Molecular Ecology, 2011, 20(3): 619-628.DOI:10.1111/j.1365-294x.2010.04959.x. [31] MOHR K I, TEBBE C C. Diversity and phylotype consistency of bacteria in the guts of three bee species (Apoidea) at an oilseed rape field[J]. Environmental Microbiology, 2006, 8(2): 258-272.DOI:10.1111/j.1462-2920.2005.00893.x. [32] MORAN N A, HANSEN A K, POWELL J E, et al. Distinctive gut microbiota of honey bees assessed using deep sampling from individual worker bees[J]. PLoS ONE, 2012, 7(4): e36393.DOI:10.1371/journal.pone.0036393. [33] OLOFSSON T C, VASQUEZ A. Detection and identification of a novel lactic acid bacterial flora with in the honey stomach of the honeybee Apis mellifera[J]. Current Microbiology, 2008, 57(4): 356-363. DOI:10.1007/s00284-008-9202-0. [34] ENGEL P, MARTINSON V G, MORAN N A. Functional diversity within the simple gut microbiota of the honey bee[J]. PNAS, 2012,109(27): 11002-11007. DOI:10.1073/pnas.1202970109. [35] ENGEL P, MORAN N A. Functional and evolutionary insights into the simple yet specific gut microbiota of the honey bee from metagenomic analysis[J]. Gut Microbes, 2013, 4(1): 60-65.DOI:10.4161/gmic.22517. [36] BEKKER A, HOLLAND H D, WANG P L, et al. Dating the rise of atmospheric oxygen[J]. Nature, 2004, 427: 117-120. DOI:10.1038/nature02260. [37] MATTILA H R, RIOS D, WALKER-SPERLING V E, et al.Characterization of the active microbiotas associated with honey bees reveals healthier and broader communities when colonies are genetically diverse[J]. PLoS ONE, 2012, 7(3): e36962.
Analysis of Microbial Community Diversity in Bee Bread by High-throughput Sequencing LIU Yuejia1,2, HAN Yejun2, PENG Wenjun1, NIU Qingsheng3, FANG Xiaoming1, ZHAO Yazhou1, TIAN Wenli1,*
(1. Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China;2. Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China；3. Jilin Province Institute of Apicultural Science, Jilin 132013, China) Abstract: Bee bread is a fermented food with great development potential. The objective of this study was to analyze the dynamic changes in microbial communities of bee bread during fermentation, and to evaluate the effect of fermentation time on the microflora compositions. Bee bread samples of Apis mellifera carnica fermented for 0.5, 4, 7 and 12 d were collected,and their microbial 16S rDNA sequences were evaluated by high-throughput sequencing. The microbial community diversity and composition in these samples were comparatively analyzed. A total of 27 phyla and 112 genera of bacteria were identified in the 4 samples. Besides, there were 25 genera common to the 4 samples. Moreover, the sequencing depth in this study effectively covered all the microbial species, along with the fermentation process, microbial richness and diversity were increased fi rst and then decreased. Oxyphotobacteria and Proteobacteria were dominant at the phylum level. At the genus level, unidentified Oxyphotobacteria was the dominant microorganism in bee bread. Along with the fermentation process, the unidentified Oxyphotobacteria increased fi rst and then decreased. Keywords: bee bread; microbial community diversity; high-throughput sequencing
收稿日期：2019-04-02 基金项目：国家现代农业（蜜蜂）产业技术体系建设专项（NCYTI-43-KXJ17）；中国农业科学院科技创新工程项目（CAAS-ASTIP-2015-IAR）；国家自然科学基金面上项目（31871861） *通信作者简介：田文礼（1980—）（ORCID: 0000-0003-3803-9616），男，研究员，博士，研究方向为功能食品与生物活性物质。E-mail: leroytian@126.comDOI:10.7506/spkx1002-6630-20190402-023 中图分类号：Q939.9 文献标志码：A 文章编号：1002-6630（2020）10-0094-07 引文格式： 刘玥佳, 韩业君, 彭文君, 等. 基于高通量测序技术分析蜂粮微生物多样性[J]. 食品科学, 2020, 41(10): 94-100.DOI:10.7506/spkx1002-6630-20190402-023. http://www.spkx.net.cnLIU Yuejia, HAN Yejun, PENG Wenjun, et al. Analysis of microbial community diversity in bee bread by high-throughput sequencing[J]. Food Science, 2020, 41(10): 94-100. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20190402-023. http://www.spkx.net.cn
| 
|Archiver|手机版|小黑屋|www.openhelp100.com
( 冀ICP备19026749号-1 )
发表于 2021-2-10 10:44:10
