浓香型白酒酿造过程中酒醅微生物区系变化特征分析.pdf

DOI：10.13746/j.njkj.2022273作者简介：陈进。浓香型白酒酿造过程中酒醅微生物区系变化特征分析陈进1，3，陈萍2，谭光迅2，杨博2，龚大春1，3，4（1.湖北省生物酵素工程技术研究中心(三峡大学)，湖北 宜昌 430072；2.湖北稻花香酒业股份有限公司，湖北 宜昌 443112；3.三峡大学生物与制药学院，湖北 宜昌 443002；4.中国轻工业功能酵母重点实验室(三峡大学)，湖北 宜昌 430072）摘要：采用高通量测序技术研究浓香型白酒发酵过程中窖池内不同深度糟醅的细菌群落组成，探讨理化因子的变化规律及与细菌群落结构变化的相关性，为提高浓香型白酒的品质提供生物学理论依据。结果表明，不同深度糟醅的群落结构无明显差异，不同发酵阶段糟醅的微生物菌群差异相对较大。入窖时乳酸杆菌属（Lactobacillus）、假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）、葡萄球菌属（Staphylococcus）是主要的优势菌属，发酵前期醋酸菌（Acetobacter）是主要的菌属（52.74%），乳酸杆菌属（Lactobacillus）在中后期占绝对优势，最高相对丰度达96.42%。冗余分析（RDA）表明，还原糖、乙醇、含水量、乳酸、乳酸乙酯、乙酸乙酯和丁酸乙酯是与菌群相关性较高的理化因子。上述研究揭示了浓香型白酒发酵过程中微生物区系的特征和变化规律，以及与理化性质的相关性，可为科学认识浓香型白酒发酵过程与机制提供理论依据。关键词：浓香型白酒；酒醅；微生物菌群；动态演变中图分类号：TS262.3；TS261.1文献标识码：A文章编号：1001-9286（2023）08-0017-07Changes of Microbial Flora in Fermented Grainsof Nongxiang Baijiu during FermentationCHEN Jin1,3,CHEN Ping2,TAN Guangxun2,YANG Bo2and GONG Dachun1,3,4(1.Hubei Research Center of Bioenzyme Engineering Technology,China Three Gorges University,Yichang,Hubei 430072;2.Daohuaxiang Distillery Co.Ltd.,Yichang,Hubei 443112;3.College of Biological and Pharmaceutical Sciences,China Three Gorges University,Yichang,Hubei 443002;4.Key Laboratory of Functional Yeastof China Light Industry,China Three Gorges University,Yichang,Hubei 430072,China)Abstract:High-throughput sequencing technology was used to study the bacterial community in the fermented grains of NongxiangBaijiu at different depths in the pit during fermentation,and the changes of physicochemical factors in the pit and their correlationwith the changes of bacterial community structure were explored.The results showed that there was no significant difference in bacte-rial community structure at different depths,while there was obvious difference in the microbial flora in the fermented grains at differ-ent fermentation stages.Lactobacillus,Pseudomonas,Bacillus and Staphylococcus were the dominant bacteria when the fermentedgrains were put into the pit.Acetobacter was dominant(52.74%)at the early fermentation stage,while Lactobacillus was dominant(96.42%)at the middle and late fermentation stages.Redundancy analysis(RDA)showed that reducing sugar,ethanol,water content,lactic acid,ethyl lactate,ethyl acetate and ethyl butyrate were the physicochemical factors with high correlation with bacterial commu-nity.The study has revealed the characteristics and change rules of microbial flora in the fermented grains of Nongxiang Baijiu duringfermentation,as well as their correlation with the physicochemical properties,thus providing theoretical basis for scientific under-standing of the fermentation process and mechanism of Nongxiang Baijiu.Key words:Nongxiang Baijiu;fermented grains;microbial flora;dynamic evolution酿酒科技2023年第8期（总第350期）LIQUOR-MAKING SCIENCE&TECHNOLOGY2023 No.8(Tol.350)1717酿酒科技2023年第8期（总第350期）LIQUOR-MAKING SCIENCE&TECHNOLOGY2023 No.8(Tol.350)中国白酒作为世界六大蒸馏酒之一，广受世界人民的喜爱。根据主体香型的不同，中国白酒大体上可分为三大类：浓香型、酱香型和清香型。浓香型白酒以其“窖香浓郁、绵甜爽净、香味协调、余味悠长”的特点深受消费者喜爱1。浓香型白酒的生产以泥窖窖池为基础，发酵过程是栖息在窖池糟醅、窖泥中的庞大微生物区系在糟醅固、液、气三相界面的复杂的物质能量代谢过程2，以富含淀粉的粮谷类为原料，以酒曲为糖化剂，采用续糟拌和、混蒸混烧、凉醅下曲、分层分批入窖、固态发酵、蒸馏出酒、精心勾兑等工艺进行生产3。浓香型白酒的酿造过程是传统的封闭式固态发酵，复杂的微生物群落除了受到环境微生物、原料、环境温度、环境湿度等外源因素的影响，还受到群落自身演替、代谢产物变化、生物热变化、含水量变化等内源发酵条件的影响。目前国内主要对浓香型白酒的酒曲、窖泥以及短期发酵过程的微生物菌群进行了跟踪研究，对于长期的窖池发酵以及空间上菌群的结构差异研究不足。对发酵过程中窖池内微生物区系变化特征的研究就是对在窖池环境下微生物生态学的研究。微生物生态学的研究内容包括确定微生物的数量、群体结构和活性，对于保护和改善生态环境、构建人工生态系统、提高废水生物处理的效率有着重要的科学意义和实用价值4。传统的微生物鉴定方法操作简单，可以直接获取目标微生物，一直沿用至今。自然界中超过 99%的微生物都是不可直接培养的，绝大多数微生物不能通过纯培养的方式获得。传统分离培养的方法工作量大，不能完整反映微生物复杂的群落组成5。施安辉6采用稀释平板涂布的方法，发现了酒醅中的细菌、酵母菌和霉菌的数量变化特征；吴衍庸等7考查了窖池中不同位置窖泥中六大类厌氧菌（甲烷菌、梭状芽孢杆菌属、乳酸菌、硫酸盐还原菌、硝酸盐还原菌及厌氧异氧菌）的差异性；岳元媛等8采用厌氧培养方法对窖泥中的细菌进行筛选分离，共筛选到8个属的细菌，绝大多数为兼性厌氧细菌，梭菌属细菌仅占1%左右，这与采用未培养技术得到的结果大相径庭，采用传统培养技术对环境微生态的研究不足以全面准确的发现微生物的菌群结构和演替特征。随着科学的进步，分子生物学手段逐渐进入了生态学研究的视野并成为生态学研究的重要技术，“高通量测序技术（High-Throughput Sequencing）”以其测序量大、灵敏、高速、高信息量、低成本等特点成为微生态学研究中应用的热门技术9。高通量测序技术是对第一代测序技术的一次技术革新，可以同时对千万条DNA序列进行测序，所以高通量测序技术也被称为“第二代测序技术”“下一代测序技术（Next Generation Sequencing，NGS）”和“深度测序”10。高通量测序为研究人员和科学家提供了一种革命性的工具，可以快速、经济地获取 DNA 序列。与传统的Sanger测序相比，使用高通量测序获得的每单位DNA序列（即碱基对）显著更具成本和时间效益。雷振河11应用高通量测序技术对清香型白酒酿造用大曲和酒醅中微生物构成进行了分析，发现大曲中优势真核微生物主要包括曲霉菌属、热子囊菌属、根霉菌属、嗜热真菌属、假丝酵母，醋酸杆菌科和乳杆菌科是清香型发酵酒醅中优势原核微生物，占原核微生物70%以上；邓杰等12利用高通量测序技术研究浓香型白酒窖池窖泥微生物群落结构，发现Clostridiales Family XI.IncertaeSedis随着窖龄的增加呈现减少的趋势，Actinobac-teria类群和Synergistetes类群随着窖龄的增加呈现增加的趋势，Haloplasmataceae类群和Clostridiaceae类群只在30年窖龄的窖池中有检测出，Lactobacil-laceae类群和Porphyromonadaceae类群只在30年窖龄的一口窖池中是优势细菌类群；李欣等13通过高通量测序技术对酱香型白酒福矛窖酒发酵的不同时期酒醅微生物进行多样性分析，发现堆积发酵酒醅的优势微生物为乳球菌属、醋酸杆菌属、芽孢杆菌属、梭形杆菌属、假丝酵母属、曲霉属，窖内发酵酒醅的优势微生物为乳杆菌属、醋酸杆菌属、青霉菌属、毕赤酵母属、曲霉属、嗜热子囊菌属、假丝酵母属；李斌等14基于 Illumina Miseq 测序方法研究了浓香型白酒中高温酒曲和芝麻香型白酒高温酒曲的微生物菌群结构，结果表明浓香型白酒中高温酒曲的优势类群主要有9个属，包括乳杆菌属（Lac-tobacillus）、魏斯氏菌属（Weissella）和片球菌属1818（Pediococcus）等，芝麻香型白酒高温酒曲的优势类群主要有芽孢杆菌属（Bacillus）和别样芽孢杆菌属（Allobacillus）；栗连会等15利用16S rRNA基因高通量测序技术比较了大曲、酒醅和窖泥中的乳酸菌群落结构，并解析了酒醅乳酸菌在发酵过程中的动态变化，结果表明，发酵过程酒醅中存在 49 种乳酸菌，其中16种来源于大曲，3种来源于窖泥，酒醅乳酸菌群落结构在发酵的第一周出现剧烈变动，发酵2 d后Lactobacillus acetotolerans取代Weissella con-fusa成为优势微生物并保持至发酵结束。本研究拟利用高通量测序技术，探究发酵过程中微生物结构以及主要的功能微生物，结合环境因子，分析其对浓香型白酒发酵过程微生态结构的影响。了解浓香型白酒酿造过程中微生物资源状况，进一步认识其生态功能，加强对窖池微生物资源的利用。有利于改善白酒的生产工艺，为进一步提高浓香型白酒的品质、人工改造窖池微生物菌群等提供生物学依据。1材料与方法1.1材料、试剂样品：取样于湖北某酒厂生产一车间中的某一窖池，选取发酵1 d、3 d、9 d、15 d、30 d、70 d、90 d为发酵节点，采用四点取样法对窖池四个角共取30 g左右样品混合，分别取窖池每个节点的上层（距离窖顶40 cm）、中层（距离窖顶1 m）、底层（距离窖底40 cm）三层糟醅的样品，共计21个样本（样本编号见表1）。混合均匀后分成两份样品，一份-20 保存并及时进行酒醅微生物基因组DNA提取，另一份进行理化性质和有机酸含量检测。试剂及耗材：葡萄糖、3,5-二硝基水杨酸、氢氧化钠、酒石酸钾钠、重蒸苯酚、亚硫酸钠、四氢呋喃、乙酸乙酯、乙醇、乙酸、丁酸、己酸、乳酸（均为分析纯），天津科密欧化学试剂有限公司；Fast DNATM-SPIN Kit for Soil，美国MP Biomedicals公司。1.2试验方法1.2.1酒醅相对酸度的测定取5 g酒醅，加入45 mL蒸馏水，置于漩涡振荡器上混合均匀。静置5 min，测定上清液 pH值。1.2.2酒醅还原糖的测定采用 3,5-二硝基水杨酸(简称DNS)分光光度计法于540 nm波长处测定。样品预处理：取酒醅样品 1.0 g 于 25 mL 比色管中，加入 10 mL 蒸馏水 100 水浴 20 min，冷却后，取上清液 1 mL 于1.5 mL EP管中8000 r/min离心10 min取上清液。1.2.3高通量测序与数据分析高通量测序由生工生物工程(上海)股份有限公司完成。高通量测序针对细菌 16S rRNA 基因V3V4 区域，扩增引物为 338F(5-ACTCCTAC-GGGAGGCAGCAG-3)和 806R(5-GGACTACH-VGGGTWTCTAAT-3)。测 序 平 台 为 IlluminaMiSeq。去掉 Barcode 序列和引物序列，得到有效的序列文件，在此基础上对测序数据的质量进行控制，去掉测序质量不好的序列，保留测序长度大于400 bp的序列，并去除嵌合体序列，得到合格的有效数据。使用Usearch 将序列按照97%相似度聚类，进 行 OTU(operational taxonomic units，操作分类单元)划分，提取代表序列，得到OTU表。将细菌利用 Mothur 方法与 SILVA 的 SSU rRNA 数据库进行物种注释分析。使用 Mothur 软件计算样品的 Chao1、Ace、Shannon、Simpson 和 Goods-cover-age指数。1.2.4数据分析酒醅理化指标的试验结果用平均值标准差表示；利用 R 语言（3.6.3版本）进行物种多样性聚类分析和相关性分析。2结果与分析2.1发酵过程中酒醅pH值、还原糖的变化规律浓香型白酒发酵过程中pH值、还原糖的变化如图1所示。三层糟醅的pH值呈先上升后下降的趋势，09 d，酒醅pH值由4.80下降至3.30左右，中层下降幅度大；990 d，pH值总体呈上升趋势，表1糟醅样品编号层次上层中层底层0 dZ01UZ01MZ01D3 dZ03UZ03MZ03D9 dZ09UZ09MZ09D15 dZ15UZ15MZ15D30 dZ30UZ30MZ30D70 dZ70UZ70MZ70D90 dZ90UZ90MZ90D陈进，陈萍，谭光迅，杨博，龚大春 浓香型白酒酿造过程中酒醅微生物区系变化特征分析1919酿酒科技2023年第8期（总第350期）LIQUOR-MAKING SCIENCE&TECHNOLOGY2023 No.8(Tol.350)上升幅度较小。三层糟醅还原糖含量总体呈下降趋势，由发酵初始的 100 mg/g 干糟醅下降到10 mg/g干糟醅。2.2酒醅细菌群落-多样性分析本研究利用高通量测序技术对浓香型白酒发酵过程酒醅中细菌群落结构进行分析。21个样品测序深度为 5442884789 条，经过对测序获得数据的拼接和过滤处理，21个样品平均每个样品有54049条有效序列，613个OTU。每个样品中有效序列数为 14891317，OTU 数量在 261058 之间，每 个 样 品 有 效 序 列 的 覆 盖 率 为 98.8%99.8%，表明基于有效序列的-及-多样性分析能较好地反映酒醅微生物群落信息。通 过-多 样 性 分 析 可 知，Chao1 指 数 在89.331618.53之间，Ace指数在155.692237.38，Shannon 指数在 0.263.28 之间，Simpson 指数在0.090.93 之间。从原核微生物群落演替性来看（图2），随着发酵的进行，Chao1和Ace指数表征菌群的丰富度（richness），Chao1 和 Ace 指数越大，菌群的丰富度越高。两个指数总体上呈先上升后下降的趋势，表明糟醅中的原核微生物丰富度先升高后降低，可能是窖泥和环境中的原核微生物菌群在糟醅中大量生长繁殖，后期因窖池环境恶化导致原核微生物菌群多样性降低。底层糟醅在90 d时原核微生物丰富度有明显的下降，底层有机酸大量聚集，较上层和中层酸度和渗透压更低。Shannon及Simpson 指数表征物种的均匀度，Shannon 指数越高、Simpson指数越小，菌群的均匀度越高。这两个指数的变化趋势在时间上波动性较大，在空间上也有明显的区别，表明菌群的均匀度变化较大。将这几个指标与发酵时间做差异性分析，Chao1、Ace和Shannon 指数与发酵时间的差异性指数为 0.017、0.016、0.046、0.162，均小于 0.05，表明不同发酵时间，菌群多样性有显著性差异。不同层糟醅指数的差异性均大于0.05，表明在空间上，不同层糟醅的多样性不存在显著差异。2.3糟醅细菌群落组成根据高通量测序分析结果，选取在属分类水平上相对丰度排名前十的菌属，生成物种相对丰度分布柱状图，见图3。10个优势菌分别为乳酸杆菌属（Lactobacillus）、醋酸杆菌属（Acetobacter）、假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）、葡萄球菌属（Staphylococcus）、Sphingopyxis、Rummeliiba-cillus、魏斯氏菌属（Weissella）、片球菌属（Pediococ-cus）、根瘤菌（Rhizobium）。发酵初始菌群与其他发酵节点菌群有较大差异，主要有乳酸杆菌属（Lacto-bacillus）、假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）、葡萄球菌属（Staphylococcus），在三层糟醅中平均占比分别为 19.11%、17.16%、12.84%、9.95%。发酵第3天，醋酸杆菌属（Acetobacter）占比较大，大约为52.74%，在13 d大量增殖后丰度逐渐降低，发酵30 d时丰度小于1%。醋酸杆菌是好氧菌，13 d 内窖池内存在一定量的溶解氧含量，丰度迅速升高，随后溶解氧含量降低，其丰度逐渐下降。乳酸杆菌属（Lactobacillus）在发酵130 d 丰度急剧升高，由 19.11%上升到 96.42%，图1糟醅理化指标(pH值及还原糖)发酵变化20207090 d乳酸杆菌属（Lactobacillus）相对丰度有所下降，仍占较高优势。乳酸杆菌是兼性厌氧菌，在氧气含量较低时也能生长繁殖。醋酸杆菌的主要发酵产物是醋酸，是丁酸、己酸生成的前体物质，乳酸菌的主要发酵产物有乳酸、乙醇、乙酸等，是乳酸乙酯和乙酸乙酯生成的前体物质，乳酸菌的产酸特性能维持窖池内的低 pH 值环境，抑制杂菌的生长。因此，醋酸杆菌和乳酸杆菌对白酒酿造过程有重要影响，是各种风味产物形成的基础。2.4酒醅细菌群落聚类分析每个样品中含量0.5%的属定义为优势属并将其用于多样性分析，以 R 语言 vegan 包计算Bray-Curtis指数，绘制聚类树图。从图4可看出，21个样品可以聚类为3簇，第一簇包含了1 d和3 d的样品，第二簇包含了9 d和15 d的样品，第三簇包含了30 d、70 d和90 d的样品。表明糟醅样品可按发酵时间分为三个阶段，即发酵前期（13 d）、发酵中期（315 d）、发酵后期（1570 d）。PCA分析图与聚类树图的结果一致。另外，在窖池不同层面上，同一发酵时间三层酒醅能够聚类在一起，表明空间上菌群的差异较小且弱于时间的差异，中层与底层的菌群相似度更高，与上层糟醅菌群存在一定差异。Adonis多因素方差分析结果显示，上述分组对样品差异的解释度为0.60184，P值为0.001，组间差异性显著。充分说明浓香型白酒发酵每个阶段的细菌菌群组成发生了明显的变化，且将发酵过程图2糟醅发酵过程中CHAO1、ACE、SHANNON及SIMPSON多样性指数分析图3在属水平浓香型白酒糟醅细菌群落组成陈进，陈萍，谭光迅，杨博，龚大春 浓香型白酒酿造过程中酒醅微生物区系变化特征分析2121酿酒科技2023年第8期（总第350期）LIQUOR-MAKING SCIENCE&TECHNOLOGY2023 No.8(Tol.350)分为3个发酵阶段是合理的。2.5细菌群落演替与环境因素变化的相关性采用 RDA 分析浓香型白酒发酵过程中细菌群落结构与 11个理化因子之间的关系，结果发现（图5），前2个排序轴分别解释群落结构与理化因子变化率的 82.45%和15.22%。还原糖与pH值之间夹角为锐角，二者呈正相关性，且均与其他理化指标呈负相关性。还原糖、乙醇、含水量、乳酸、乳酸乙酯、乙酸乙酯和丁酸乙酯是对菌群结构影响较大的理化因子。RDA结果还揭示了主要的细菌类群与注：a-聚类树图；b-PCA分析图图4不同发酵时间及深度糟醅样品多样性分析聚类树图及PCA图图5糟醅细菌菌群与环境因子的RDA分析2222理化因子的关系，具体而言，Bacillus、Pseudomo-nas、Staphylococcus、Weissella 与 pH 值呈正相关，Lactobacillus与乙醇、含水量、丁酸乙酯、乙酸乙酯、乳酸乙酯、己酸乙酯和乳酸呈正相关。3结论本研究以浓香型白酒发酵过程中糟醅为研究对象，通过高通量测序技术探究了浓香型白酒发酵过程中的细菌多样性。假单胞菌属（Pseudomo-nas）、芽孢杆菌属（Bacillus）、葡萄球菌属（Staphylo-coccus）、醋酸杆菌属（Acetobacter）为发酵前期的主要菌属，在发酵前期（19 d），群落物种丰富，随着发酵继续进行，细菌群落组成发生显著改变，最终乳酸杆菌属成为相对优势菌属。对浓香型白酒发酵过程中的理化特性变化情况进行了分析，总体而言，随着发酵的进行，pH值和还原糖含量降低。由RDA分析可知，糟醅中微生物的种群结构与糟醅的理化指标存在明显关联，且还原糖含量与微生物物种分布相关程度最大，乳酸杆菌属对糟醅理化指标的影响较大。本研究探究了浓香型白酒糟醅在不同发酵阶段中细菌群落及理化特性的变化规律及相互关系，有助于构建关于浓香型白酒糟醅微生物多样性的系统认识体系，后续可通过分析发酵过程中微生物具体代谢途径和机理探究其对糟醅及浓香型白酒的风味和品质的作用，以期为浓香型白酒工业化发展提供理论基础。参考文献：1吴飞.浓香型白酒窖池发酵过程中酵母类微生物的分析D.成都：四川大学,2006.2郭来虎.中国第一窖M.北京：中国工人出版社,1998：28-29.3胡晓龙,王康丽,余苗等.浓香型酒醅微生物菌群演替规律及其空间异质性J.食品与发酵工业,2020,46(10)：66-73.4谢冰,徐亚同.环境微生物生态研究方法进展C/中国生态学学会.生态安全与生态建设中国科协2002年学术年会论文集.2002：123-128.5严超,红枣白兰地发酵过程中细菌群落结构及与风味物质的相关性研究D.保定：河北农业大学,2018.6施安辉,浓香型白酒发酵过程窖中微生物区系的分析J.酿酒,1986(4)：24-29.7WU YY,XUE T R,CHEN Z R,et al.The distributionand role of anaerobic bacteria in older pit mud ofWuliangyeJ.Acta microbiologica sinica,1991,31(4)：299-307.8YUE YY,ZHANG W X,LIU X,et al.Isolation andidentification of facultative anaerobes in the pit mud ofChinese Luzhou-flavor liquorJ.Microbiology,2007,34(2)：251-255.9PORTILLO M D C,MASA.Analysis of microbialdiversity and dynamics during wine fermentation ofGrenache grape variety by high-throughput barcodingsequencingJ.LWT food science and technology,2016,72：317-321.10MARDIS E R.The impact of next-generationsequencing technology on geneticsJ.Trends in genetics,2008,24(3)：133-141.11雷振河.采用高通量测序技术分析清香型白酒酿造微生物J.食品与发酵工业,2015,41(9)：164-167.12邓杰,黄治国,卫春会,等.基于高通量测序的浓香型白酒窖池细菌群落结构分析J.现代食品科技,2015,31(7)：50-55.13李欣,王彦华,林静怡,等.高通量测序技术分析酱香型白酒酒醅的微生物多样性J.福建师范大学学报(自然科学版),2017,33(1)：51-59.14李斌,闫志鹏,李慧星,等.基于高通量测序技术的浓香型和芝麻香型白酒酒曲细菌群落结构分析J.中国酿造,2018,37(8)：148-152.15栗连会,肖辰,陆震鸣,等.泸型酒发酵酒醅中乳酸菌群落的来源、演替规律及功能预测J.食品与生物技术学报,2018,37(12)：1242-1247.陈进，陈萍，谭光迅，杨博，龚大春 浓香型白酒酿造过程中酒醅微生物区系变化特征分析更 正 说 明我刊2023年第7期“目次”中 不同窖龄窖泥产酒发酵前后的理化指标、风味物质含量和细菌菌群结构差异分析 一文的第三作者“韩怀武”应为“韩怀斌”。特此更正。酿酒科技杂志社2023年8月18日2323