浓香型酒醅中梭菌群落演替及其与理化性质的相关性分析_胡晓龙.pdf

食品与发酵工业FOOD AND FEMENTATION INDUSTIES942023 Vol.49 No.5(Total 473)DOI:10 13995/j cnki 11 1802/ts 031918引用格式:胡晓龙，冯大鸿，张勇，等 浓香型酒醅中梭菌群落演替及其与理化性质的相关性分析J 食品与发酵工业，2023，49(5):94 100 HU Xiaolong，FENG Dahong，ZHANG Yong，et al Dynamic clostridial community succession and its corre-lation with physicochemical factors in Jiupei for strong flavor Baijiu productionJ Food and Fermentation Industries，2023，49(5):94 100浓香型酒醅中梭菌群落演替及其与理化性质的相关性分析胡晓龙1，冯大鸿1，张勇1，迟雷1，郑燕1，韩素娜2，沈祥坤3，魏涛1*1(郑州轻工业大学 食品与生物工程学院，河南 郑州，450001)2(河南仰韶酒业有限公司博士后科研工作站，河南 三门峡，47240)3(河南省食品工业科学研究所有限公司，河南 郑州，450053)摘要采用高通量测序和实时荧光定量 PC 技术解析了浓香型酒醅发酵过程中梭菌纲细菌群落多样性、组成及演替规律，并对酒醅梭菌群落及其理化因子进行了冗余分析(redundancy analysis，DA)和 Pearson 相关性分析。结果表明，梭菌的多样性和相对丰度在发酵第 7 天最低，45 d 达到最高。供试样品中梭菌群落由 1 个目、13 个科和 39 个属组成，其中 Caproiciproducens(己酸菌属)和 Clostridium_sensu_stricto_1(狭义梭菌属 1)是主要的菌属，均呈现先减少后增加再减少的趋势。DA 结果表明，理化因子的贡献率依次为温度(45 5%)含水量(17 8%)铵态氮(11 0%)还原糖(10 3%)乙醇(8 0%)pH(5 2%)总酯(1 5%)总酸(0 6%)，温度与梭菌群落结构呈极显著关系(P 0 01)。Caproiciproducens、Fonticella、Anaerofilu 与总酸呈显著正相关(0.01 P 0 05)，与含水量呈极显著正相关(0 001 P 0 01)。此外，含水量与 Butyrivibrio_2(丁酸弧菌属)具有更强的相关性(0 000 1 P 0 001)，温度与相对丰度前 10 种梭菌均呈负相关。该研究丰富了浓香型白酒发酵过程中对梭菌群落及其影响因素的认识，为发酵过程中梭菌群落的调控提供了理论支撑。关键词浓香型白酒;酒醅;梭菌;理化性质;相关性第一作者:博士，副教授(魏涛教授为通信作者，E-mail:weit8008 zzuli edu cn)基金项目:国家自然科学基金项目(31801535);河南省重点研发与推广专项项目(222102110214);中原产业创新领军人才项目(214100510009);中原科技创新领军人才项目(224200510017);三门峡市引进高层次创新创业人才(团队)产业化项目(2020111)收稿日期:2022-04-11，改回日期:2022-05-17中国白酒根据酿造工艺、原料和风味特征等差异可将其分为浓香型、酱香型、清香型和米香型等多种香型1。其中浓香型白酒因其具有窖香浓郁、绵柔醇厚、香味协调、回味悠长的风格特点而深受广大消费者的青睐2。酒醅是浓香型白酒酿造微生物的主要载体，微生物群落多样性及组成是影响浓香型白酒独特风味和品质的主要因素之一3。近年来，高通量测序技术的广泛应用极大地丰富了人们对白酒酿造微生物的认识。胡晓龙等3 发现酒醅中微生物主要为厚壁菌门(Firmicutes)、拟杆菌门(Bacteroidetes)和变形菌门(Proteobacteria)等，其中厚壁菌门在酒醅和窖泥微生物中占绝对优势，其相对丰度能达到 50%以上。梭菌纲细菌(梭菌)是浓香型白酒酿造过程中一类重要的功能菌，隶属于厚壁菌门，具有多种代谢活动，包括将淀粉转化为己酸(己酸乙酯的前体物质)、丁酸、醇等4。早在 20 世纪 60年代梭菌就被认为是浓香型白酒重要风味物质己酸乙酯的重要贡献菌5。近年来已经从不同酒企的窖泥中分别分离得到克鲁氏梭菌(Clostridium kluyve-ri)6、耳蜗形梭菌(C luticellarii)7 和丁酸梭菌(Cbutyrieum)8 等，并证明和产己酸、丁酸密切相关。尽管梭菌是酒醅中的优势微生物之一，但是目前关于梭菌的研究主要集中在窖泥样品及梭菌的代谢物分析中9 11，而对酒醅样品中梭菌群落多样性及其演替规律研究鲜有报道。因此，本研究采用高通量测序技术结合实时荧光定量方法对酒醅梭菌群落的多样性、组成及其含量动态变化进行解析。同时，通过结合酒醅理化性质的变化，探索了梭菌群落和酒醅理化性质的关联性，有助于进一步揭示梭菌群落在整个酿造周期内的演替规律及其影响因素。1材料与方法1 1材料与试剂样品采集:随机选取正常发酵的 3 口平行窖池，对每口窖池中发酵第 3、7、15、25、45、60 天的酒醅进研究报告2023 年第49 卷第5 期(总第473 期)95行跟踪取样。对每口窖池下层(距窖底约 50 cm)平面上任一对角线的 1/2、1/4 及其距窖壁 20 cm 处 3个不同位置的酒醅分别采集并混匀作为该窖池的代表性酒醅样品。因此，本研究共选取了 18 个样品，采集完成后装入已标记编号的无菌袋中并保存于30 备用，其编号分别为 B3-1、B3-2、B3-3、B7-1、B7-2、B7-3、B15-1、B15-2、B15-3、B25-1、B25-2、B25-3、B45-1、B45-2、B45-3、B60-1、B60-2、B60-3。以编号B3-1、B3-2、B3-3 为例，其中 B 代表下层酒醅，第一个数字 3 代表发酵时间为 3 d，-1、-2 和-3 代表3 个平行样品。主要试剂:浓硫酸、酒石酸钾钠、苯酚、KCl、亚硝基铁氰化钠、NaClO、NH4Cl，天津市科密欧化学试剂有限公司;K2Cr2O7，天津市安吉瑞化工有限公司;酚酞、邻苯二甲酸氢钾、(NH4)2SO3，上海生工生物工程有限公司;DNS，上海麦克林生化科技有限公司。1 2仪器与设备MP200A 精密电子天平，上海良丰仪器仪表有限公司;C1000 温度梯度 PC 仪，伯乐生命医学产品(上海)有限公司;TGL-20M 高速离心机，上海卢湘仪离心机仪器有限公司;DYY-8C 电泳仪，北京市六一仪器厂;Mixer4K 微型涡旋混合仪，生工生物工程(上海)股份有限公司;PHS-3C 型 pH 计，雷磁上海仪电科学仪器股份有限公司;Multiskan FC 酶标仪，赛默飞世尔科技公司;电热鼓风干燥箱，上海一恒科学仪器有限公司。1 3实验方法1 3 1酒醅理化测定含水量测定:采用 105 干燥箱烘干恒重法测定。pH 检测:将酒醅按料液比 1 3(g mL)与超纯水混合均匀，然后用 pH 计检测。还原糖、乙醇、铵态氮、总酸、总酯均参照文献 12的方法进行测定。1 3 2酒醅样品 DNA 提取按照 PowersoilDNA Isolation Kit 试剂盒说明书提取酒醅 DNA，然后利用 1%琼脂糖凝胶电泳检测提取 DNA 质量。1 3 3酒醅 DNA 检测使用引物 338F:5-ACTCCTACGGGAGGCAGCA-3，806:5-GGACTACHVGGGTWTCTAAT-3 扩增合格总 DNA 样品，并对其产物进行纯化、定量和均一化形成测序文库，建好的文库先进行文库质检，质检合格的文库用 Illumina HiSeq 2500 进行测序。1 3 4生物信息分析首先使用 Trimmomatic v0 33 软件对测序得到的aw eads 进行过滤，然后使用 cutadapt 1 9 1 软件进行引物序列的识别与去除，得到不包含引物序列的高质量 eads;使用 FLASH v1 2 7 软件，通过 overlap对每个样品高质量的 eads 进行拼接，得到的拼接序列即 Clean eads;使用 UCHIME v4 2 软件，鉴定并去除嵌合体序列，得到最终有效数据(effective reads)。其中相似性97%的序列归为 1 个操作分类单元(operational taxonomic unit，OTU)，将 OTU 代表序列与 Silva 数据库进行比对和物种注释。挑选梭菌纲细菌结果用于后续分析。1 3 5实时荧光定量 PC以细菌通用引物 对 27F/1492 扩增菌株 Clos-tridium butyricum D9(本实验室保存)的 16S rNA 序列全长，将该 PC 扩增产物连接 pMD19-T 载体，转化大肠杆菌 JM109，培养后挑选出阳性克隆子。提取阳性克隆子质粒测定其吸光值并换算为质粒拷贝数(copies)，10 倍梯度稀释为 DNA 模板。采用细菌通用引物对 515F/806 和梭菌引物对 SJ-F/SJ-分别对上 述 模 板 进 行 实 时 荧 光 定 量 PC 扩 增，以lg(copies/L)为横坐标，对应的 Ct 值为纵坐标，分别绘制细菌和梭菌的荧光定量标准曲线，其中细菌和梭菌荧光定量 PC 扩增体系及条件分别参照文献 13和 14 的方法。以每个酒醅样品的 DNA 为模板利用上述引物和方法并测定其 Ct 值，分别根据细菌和梭菌荧光定量标准曲线计算每个样品中的细菌及梭菌的含量。1 3 6数据处理使用 SPSS 26 进行理化因子显著性分析，选择Duncan s 多重比较方法，不同小写字母表示处理之间存在显著差异(P 0 05)。Shannon 指数的计算如公式(1)所示。H=si=1(Pi ln Pi)(1)式中:H，Shannon 多样性指数;s，微生物属的总数;Pi，第 i 个属的数量占总数的比例。柱状图、标准曲线均使用 Origin 9 0 绘制。采用HemI 软件绘制热图，Canoco 5 软件绘制主成分分析(principal component analysis，PCA)图和冗余分析(redundancy analysis，DA)图。Shannon 指数箱线图和相关性热图使用 TUtool 云平台完成。食品与发酵工业FOOD AND FEMENTATION INDUSTIES962023 Vol.49 No.5(Total 473)2结果与分析2 1理化性质分析如表 1 所示，根据发酵过程中酒醅的温度变化可以将整个发酵周期分为 3 个阶段，即升温期(3 7 d)、稳定期(7 25 d)和降温期(25 60 d)，这与HU 等15 的研究结果一致。这主要是由于发酵前期酒醅营养充足，有利于微生物的生长繁殖，从而产生大量发酵热16。随着发酵时间的延长，酒醅中的营养成分逐渐减少，微生物生长和代谢减缓，酒醅温度在 34 保持 18 d 左右。在发酵后期，酒醅中可利用养分基本耗尽，加上乙醇和酸的增加，部分不适应环境的微生物逐渐死亡或者休眠17，发酵作用和微生物活动均减小，酒醅温度下降。在发酵前7 d，酒醅中除含水量下降外，其余理化指标均上升;在发酵 7 25 d，乙醇、总酸含量和含水量上升，还原糖和总酯含量开始下降铵态氮略微上升后下降;在发酵 25 60 d，含水量和总酸均在 45 d 时达到最大值，分别为(63.76 3.43)%和(24.63 2)mg/g，然后下降，铵态氮和总酯在 45 d 时略微上升而后下降;还原糖和 pH 则在 45 d 时下降至最低值(6.08 5.82)mg/g 和 3.36;乙醇在 45 d 时略微下降，然后在60 d 时升至最高值(126.42 10.37)mg/g。酒醅 pH 在发酵前 45 d 均呈现下降趋势，发酵结束时略有回升。在发酵过程中，水分是保障酒醅中微生物生长代谢的重要因素，并为其代谢产物提供有效溶剂18。还原糖含量能间接反应白酒的发酵进程19，发酵开始时，原料中的淀粉被分解为还原糖，但此时微生物生长代谢尚不活跃，对还原糖的利用率还不高，导致其含量积累增加，随着发酵进行不断减少直到发酵结束。乙醇是白酒发酵的重要产物，平均生成速率为 0 57 mg/(gd)，其质量分数在 3 7 d 时增加最快，为2 11 mg/(gd)，在发酵至45 d 时略微减少，这主要是由于在发酵后期酸和醇等合成脂类的前体物质大量积累，进而发生酯化反应生成酯类，加上酒醅含水量增加，稀释了乙醇的浓度，导致乙醇减少。另外在封闭的发酵环境中，微生物不断代谢积累有机酸，酒醅中总酸快速上升，pH 开始下降。表 1不同酒醅样品中理化性质变化Table 1Changes of physicochemical properties in different Jiupei samples发酵时间/d温度/含水量/%还原糖/(mgg1)乙醇/(mgg1)铵态氮/(mgkg1)pH总酸/(mgg1)总酯/(mgg1)327 163a58 24 266a3363 21 16ab9306 9 03a12 47 309a3 61 0 06a1143 1 65a45 43 1049ab733 5 108b57 21 02a3673 19 59a10149 9 18ab12 56 27a3 57 0 01ab1184 0 84a52 64 2021b1534 163b58 73 137ab2119 20 76ab10957 5 62abc12 86 23a35 0 03b1677 2 86b30 88 933ab2534 082b61 78 086ab856 10 8ab12484 1177c7 57 18a3 42 0 05c2123 0 16c19 81 92a4529 141a63 76 343b608 5 82b1206 7 81bc88 161a3 36 0 02c2463 2c22 46 79a602733 047a61 35 254ab654 5 83ab12642 1037c8 17 116a3 41 0 02c2314 0 95c20 27 708a注:组间不同字母表示差异显著(P 005)2 2梭菌群落 多样性分析Shannon 指数可以用来表示酒醅中物种多样性，指数越大，多样性越高。如图 1 所示，梭菌 Shannon指数的范围在 0 01 0 27，呈现先减少(3 7 d)，后稳定(7 25 d)，增加(25 45 d)，再下降(45 60 d)的趋势，梭菌群落的多样性在不断发生变化，且在第7 天达到最低，第 45 天最高，这可能与酒醅中氧含量有关。在发酵前7 d，酒醅中氧含量较多，有利于芽孢杆菌等好氧菌的生长繁殖;在发酵中后期(25 45 d)，酒醅溶氧量被消耗，更适合于梭菌等厌氧菌的生长，同时好氧菌生物量减少，进而导致梭菌群落多样性显著增加。2 3酒醅中梭菌群落-多样性解析在属水平上，共注释到 39 个梭菌纲中的属。基于上述 39 个属进行 PCA，探究不同发酵节点之间梭图 1不同酒醅样品中梭菌群落 Shannon 指数变化Fig 1Shannon index of clostridial community indifferent Jiupei samples注:*代表具有显著变化(P 005)菌的群落差异，结果如图 2 所示。18 个样品整体上可分为 4 类，即 3、45、60 d 的样品分别聚为一类，7、15、25 d 样品聚在一起，每一类样品之间具有类似的梭菌群落。由此可知，酒醅梭菌群落变化与发酵时间研究报告2023 年第49 卷第5 期(总第473 期)97有较好的关联性，其中发酵 3、45、60 d 的酒醅分别具有独特的梭菌群落，而 7、15、25 d 酒醅具有相似的梭菌群落。尽管不同发酵时间的酒醅梭菌群落能区分开，但是 B3 平行样品之间具有较大的离散度，这可能是由于:白酒酿造实际生产中的物料混合(大曲、粮食、酒醅等)多为人工拌和，很难保证不同平行窖池中物料拌匀程度完全一致，这也导致不同平行样品之间接入的微生物群落自身存在差异;此外，白酒酿造为固态发酵，平行样品间的酒醅微环境(如酒醅营养成分含量、含水量和溶氧等)可能也会存在差异，使得平行窖池之间酒醅微生物发酵进度不一，进而导致部分平行样品不能很好地聚类20 21。图 2不同酒醅样品中梭菌组成的 PCA 图Fig 2PCA plot based on the clostridia communitystructures in different Jiupei samples2 3 1酒醅梭菌群落组成根据分类学注释结果(图 3-a)，酒醅中共检测到6 个优势菌纲(1%):Bacilli(芽孢杆菌纲，平均相对丰度 78 12%)、Gammaproteobacteria(-变形菌纲，5 96%)、Bacteroidia(拟杆菌纲，5 33%)、Alphapro-teobacteria(-变形菌纲，2 97%)、Clostridia(梭菌纲，1 71%)和 Actinobacteria(放线菌纲，1 11%)。其中Clostridia 由 1 个目、13 个科和 39 个属组成。在目水平上，仅鉴定到单一的 Clostridiales(梭菌目)。在 13 个科中，平均相对丰度 0 1%的科有uminococcaceae(瘤胃菌科，0 71%)、Clostridiaceae(梭菌科，0 35%)和 Lachnospiraceae(毛螺菌科，0.35%)。表明在酒醅发酵过程中梭菌变化主要集中在上述 3 个科，与钱玮等10 的结果一致。在属水平上，共检测到 39 个可鉴定的属，平均相对含量 0.1%的仅有 Caproiciproducens(己酸菌属 0.4%)和Clostridium_sensu_stricto_1(狭义梭菌属 1，0.1%)。与窖泥中梭菌相比10，酒醅中的梭菌群落多样性较低，所检测到的梭菌数量较少，如常见的 Heliobacteri-aceae(太阳杆菌科)在酒醅中也并未检测到。导致这一结果的原因可能有:酒醅中乳杆菌的相对丰度较高，其主要终产物乳酸会抑制没有耐酸特性的梭菌生长，加上乳酸链球菌素和乳酸菌素等细菌素会抑制革兰氏阳性菌的生长，导致酒醅中梭菌较少14，22;本研究使用细菌通用引物进行测序分析，可能会因其扩增偏好性导致部分梭菌或含量低梭菌的无法检测到23。2 3 2酒醅梭菌群落组成的演替规律梭菌的平均相对丰度为 0 12%5 58%，在发酵过程中呈先减少(3 25 d)后增加(25 45 d)再减少(45 60 d)的趋势。在发酵第 7 天含量最低(0.12%)，其次为 15 d 和 25 d，在发酵 45 d 时占比最高(5 58%)(图 3-b)。在科水平上，uminooccace-ae(0 71%)和 Clostridiaceae(0 35%)和 Lachnospir-aceae(0 35%)为主要优势科，其总平均相对丰度占梭菌纲的 73 18%。酒醅梭菌群落的变化主要由上述 3 个科驱动，其中 uminococcaceae 在前 25 d 呈现减少的趋势，在发酵 25 45 d 大幅增加，随后又减少。Clostridiaceae 呈现和 uminococcaceae 相反的规律。Lachnospiraceae 则呈现不稳定的变化过程。a 细菌纲水平;b 梭菌纲水平图 3不同酒醅样品中细菌纲和梭菌纲水平相对丰度Fig 3elative abundance of bacteria and clostridiaat class level in different Jiupei samples在属水平上，挑选相对丰度前 20 个梭菌纲中的属进行热图分析(图 4)。聚类分析结果表明，整个发酵过程可分为 3 类，B3 为类，B7、B15 和 B25 为类，B45 和 B60 为类。说明在发酵过程中，梭菌会随着发酵时间发生整体性的变化，也与图 2 的结果存食品与发酵工业FOOD AND FEMENTATION INDUSTIES982023 Vol.49 No.5(Total 473)在一定的一致性。在发酵 7 25 d，所有梭菌的相对丰度均不同程度地减少，发酵至 45 d 时，除 Saccharo-fermentans(产乙酸糖发酵菌)和 Acetanaerobacterium(延长厌氧醋菌)外，其余 18 个属的相对含量均增加，其中，Caproiciproducens、Clostridium_sensu_stricto_1、Sedimentibacter、Fonticella、Anaerofilum、uminiclos-tridium_1(瘤胃梭菌属 1)和 Butyrivibrio_2(丁酸弧菌属2)的相对丰度在发酵至45 d 时增加至最高。在这7 个属中，Caproiciproducens 含量最多，达到 2 29%。GAO 等24 研究发现 Caproiciproducens 为浓香型白酒窖泥中的主要产己酸细菌，因此推测该属中的某些物种在酒醅中合成己酸起着重要的作用。图 4不同酒醅样品中梭菌属水平组成Fig 4Clostridial composition in different Jiupeisamples at the genus level2 4梭菌荧光定量分析实时荧光定量 PC 结果进一步研究了发酵过程酒醅中梭菌纲微生物的基因拷贝数以及在总细菌中的占比变化(图 5)。以浓度对数为横坐标，Ct 值为纵坐标，拟合绘制总细菌拷贝数标准曲线为:y=3 033 3x+39 561，2=0 997 3，梭菌纲拷贝数标准曲线为:y=3 093 2x+37 519，2=0 998 7，标准曲线均绘制良好。按上述标准曲线计算各样本中总细菌及梭菌微生物的 16S rNA 基因拷贝数(图 5-a)，梭菌纲的基因拷贝数呈现先减少后增加的动态变化，第 25 天时最低(基因拷贝数为 5 21 lg 个/g 窖泥)，第 45 天升至最高(7 13 lg 个/g 窖泥)，发酵结束时(6 01 lg 个/g 窖泥)较发酵开始时有所减少。这也与 16S rNA 所得梭菌相对含量的变化规律一致。通过定量发现酒醅中梭菌和总细菌的基因拷贝数比值为 0 17%4 28%(图 5-b)，与 16S rNA 测序结果类似，梭菌的占比同样在发酵45 d 时最高。a 荧光定量结果;b 基因拷贝数比值图 5不同酒醅样品中总细菌及梭菌纲细菌实时荧光定量结果Fig 5eal-time PC results of total contents of bacteriaand clostridia in different Jiupei samples2 5理化因子与梭菌之间的相关性DA 展示了酒醅中相对丰度前 10 的梭菌(平均相对丰度 004%)与理化因子之间的关系(图 6-a)，2 个主成分的总解释度为 64 68%，主要集中 Aixs1(59 76%)(图 6-a)。理化因子与梭菌之间的相关性贡献率由大到小依次为温度(45.5%)含水量(17.8%)铵态氮(11.0%)还原糖(10.3%)乙醇(8.0%)pH(5.2%)总酯(1.5%)总酸(0.6%)，其中温度与梭菌群落结构呈极显著关系(P 0.01)。不同发酵节点的酒醅样品能有效地被区分开来，再次说明酒醅中梭菌群落的变化受到发酵时间的影响，同时也受酒醅理化性质的影响。18 个酒醅样品分布在不同的象限，例如，发酵3 d 和7 d 的样品分别存在于第 1 和第 2 象限中，与铵态氮、总酯、pH 和还原糖呈现较强的正相关;第 3 象限主要有发酵 15 d 和 25 d 的样品，与温度呈正相关;发酵 45 d和 60 d 的样品主要在第 4 象限，与乙醇、总酸和含水量呈正相关。综上，在不同的发酵节点酒醅中驱动梭菌群落变化的理化因子不一样，温度、含水量和铵态氮是影响较大的环境因子。通过计算相对丰度前 10 的梭菌和理化因子之间的 Pearson 相关性系数进一步分析酒醅理化因子对梭菌的影响。如图 6-b 所示，除 Christensenellaceae_-7研究报告2023 年第49 卷第5 期(总第473 期)99_group 外，其他 9 个梭菌均与总酸和含水量呈正相关。其中，总酸与 Caproiciproducens、Fonticella、Anaer-ofilum 和 Butyrivibrio_2 呈显著正相关(0.01 P 0.05)。含水量与 uminidostridium_1 呈显著正相关(0 001 P 0 01)，与 Caproiciproducens、Fonticella、Anaerofilu 呈极显著正相关(0 000 1 P 0 001)，与Butyrivibrio_2 的相关性更强(0 000 1 P 0 001)。pH 仅与 Christensenellaceae_-7_group 呈正相关，与Caproiciproducens 呈显著负相关(0 01 P 0 05)。此外，温度与所选梭菌均呈负相关，乙醇、总酯、还原糖和铵态氮对梭菌均无显著性影响。综上，通过明晰酒醅梭菌群落与理化因子相关性，有助于通过监控酒醅理化因子(如含水量、pH 值等)间接了解酒醅中梭菌的状态，进而为提高酒醅中梭菌含量提供理论基础。a DA;b Pearson 相关性分析图 6理化因子与梭菌前 10 属相关性分析Fig 6The correlation analysis between the top 10 genera ofClostridia and physicochemical factors in Jiupei samples注:*:0 01 P 0 05;:0001 P 0 01;:0 000 1 P 0 0013结论本研究中，供试酒醅样品梭菌群落由 39 个梭菌纲中的属组成，梭菌群落的相对丰度在发酵至 7 d 时减少，然后保持稳定至 25 d，在 25 45 d 发酵中又增加，然后减少至发酵结束。其群落 Shannon 指数、丰度和基因拷贝数均在发酵 45 d 达到最高。在梭菌群落的演替过程中，Caproiciproducens 和 Clostridium_sen-su_stricto_1 主要的优势梭菌。DA 结果表明，在发酵 3 7 d，铵态氮、总酯、pH 和还原糖是驱动梭菌演替的主要因素;发酵 15 25 d 时温度成为了主要的驱动因素;到发酵 45 60 d 乙醇、总酸和含水量成为新的影响因素。Pearson 相关性分析表明，总酸和含水量具有较大影响，温度与相对丰度前 10 的梭菌呈负相关，pH 仅与 Christensenellaceae_-7_group 呈正相关，与 Caproiciproducens 呈显著负相关。乙醇、总酯、还原糖对梭菌群落均无显著性影响。总的来说，尽管高通量测序能全面地认识酒醅中梭菌的多样性，但由于使用细菌通用引物导致梭菌信息丢失导致梭菌数量较少。在进一步的研究中，应考虑设计使用梭菌特异性引物，并选择更深层次的测序方法检测酒醅中梭菌的具体代谢途径及产物。参考文献 1ZHENG X W，HAN B Z Baijiu(白酒)，Chinese liquor:History，classification and manufactureJ Journal of Ethnic Foods，2016，3(1):19 25 2XU Y Q，ZHAO J，LIU X，et al Flavor mystery of Chinese tradi-tional fermented Baijiu:The great contribution of ester compounds J Food Chemistry，2022，369:130920 3胡晓龙，王康丽，余苗，等 浓香型窖泥原核微生物群落多样性及其在窖池中的空间异质性J 食品与发酵工业，2020，46(11):77 84HU X L，WANG K L，YU M，et al Biodiversity and spatial hetero-geneity of prokaryote community in strong-flavor Baijiu fermentationpit mudsJ Food and Fermentation Industries，2020，46(11):77 84 4DE P Physiology and sporulation in Clostridium J Microbiol-ogy Spectrum，2014，2(4):TBS-0 010 2 012 5ZOU W，YE G B，ZHANG K Z Diversity，function，and applica-tion of Clostridium in Chinese strong flavor Baijiu ecosystem:A re-viewJ Journal of Food Science，2018，83(5):1 193 1 199 6HU X L，DU H，XU Y Identification and quantification of the capro-ic acid-producing bacterium Clostridium kluyveri in the fermentation ofpit mud used for Chinese strong-aroma type liquor production J In-ternational Journal of Food Microbiology，2015，214:116 122 7WANG Q，WANG C D，LI C H，et al Clostridium luticellarii spnov，isolated from a mud cellar used for producing strong aromaticliquors J International Journal of Systematic and Evolutionary Mi-crobiology，2015，65(12):4 730 4 733 8袁华伟，赵猛，蒋君，等 窖泥中产丁酸菌的筛选、鉴定及生长性能研究 J 酿酒科技，2017(5):35 40YUAN H W，ZHAO M，JIANG J，et al Isolation and identificationof a butyric acid-producing strain from pit mud and study on itsgrowth performanceJ Liquor-Making Science Technology，2017(5):35 40 9蒲秀鑫，柴丽娟，徐鹏翔，等 泸型酒窖泥中梭菌的分离及代谢产物分析 J 微生物学报，2019，59(12):2 427 2 436PU X X，CHAI L J，XU P X，et al Isolation and metabolic charac-ters of Clostridium strains from pit mud of Luzhou-flavor baijiuJ食品与发酵工业FOOD AND FEMENTATION INDUSTIES1002023 Vol.49 No.5(Total 473)Acta Microbiologica Sinica，2019，59(12):2 427 2 436 10钱玮，陆震鸣，柴丽娟，等 泸型酒酒醅与窖泥中梭菌群落结构、演替和功能差异 J 生物工程学报，2020，36(6):1 190 1 197QIAN W，LU Z M，CHAI L J，et al Differences of the structure，succession and function of Clostridial communities between Jiupeiand pit mud during Luzhou-flavour Baijiu fermentationJ ChineseJournal of Biotechnology，2020，36(6):1 190 1 197 11李聪聪 浓香型白酒窖泥梭菌群落多样性及高产丁酸梭菌发酵特性研究 D 郑州:郑州轻工业学院，2018LI C C Diversity of Clostridia communities and fermentation char-acteristics of Clostridia with high yield of butyric acid in Chinesestrong-flavor liquor pit mudsD Zhengzhou:Zhengzhou Universityof Light Industry，2018 12沈怡方 白酒生产技术全书M 中国轻工业出版社，1998SHEN Y F Baijiu Production TechnologyM China Light Indus-try Press，1998 13XIAO C，LU Z M，ZHANG X J，et al Bio-heat is a key environ-mental driver shaping the microbial community of medium-tempera-ture Daqu J Applied and Environmental Microbiology，2017，83(23):e01 550 01 517 14HU X L，DU H，EN C，et al Illuminating anaerobic microbialcommunity and cooccurrence patterns across a quality gradient inChinese liquor fermentation pit mudsJ Applied and Environ-mental Microbiology，2016，82(8):2 506 2 515 15HU X L，TIAN J，WANG K L，et al The prokaryotic communi-ty，physicochemical properties and flavors dynamics and their cor-relations in fermented grains for Chinese strong-flavor Baijiu pro-ductionJ Food esearch International(Ottawa，Ont)，2021，148:110626 16陈雪，张永利，闫宗科，等 凤香型酒醅微生物群落演替及其与理化指标的相关性分析J 食品科学，2020，41(22):200 205CHEN X，ZHANG Y L，YAN Z K，et al Correlation analysis ofmicrobial community succession in fermented grains of Xifeng-flavorChinese liquor with its physicochemical indexesJ Food Science，2020，41(22):200 205 17YAN S B，WANG S C，WEI G G，et alInvestigation of the main pa-rameters during the fermentation of Chinese Luzhou-flavour liquor J Journal of the Institute of Brewing