窖泥中挥发性物质和微生物群落的空间分布规律及其关系

张朝正，张天爽，董思文，孙 伟，赵 华,

（1.天津科技大学生物工程学院，工业发酵微生物教育部重点实验室，天津市工业微生物重点实验室，天津 300457；2.安徽文王酿酒股份有限公司，安徽临泉 236400）

窖池是浓香型白酒的发酵装置，其中窖泥的质量决定了浓香型白酒的品质[1−2]。在浓香型白酒发酵的过程中，窖泥中栖息的微生物在糟醅与窖泥接触的界面长期接触与反应，进行物质交换，进而产生各种香气成分[3−4]。研究窖泥不同空间位置微生物分布与香气成分的构成对深入理解窖泥在浓香型白酒生产中的作用具有重要意义。

窖泥中的挥发性物质的含量是表征窖泥优劣的一项重要指标，是与糟醅相互交换的物质基础。顶空固相微萃取（head-space solid-phase microextraction,HS-SPME）偶联气质联用法（gas-chromatography mass spctromrtry, GC-MS）具有操作简便、结果易于处理等优点，近年来在窖泥和酒糟中挥发性物质的研究上应用广泛[5−6]。对糟醅中香气成分的空间分布研究发现，同层糟醅中距离窖泥越近的部分酯类和酸类含量越高，且呈逐步递减趋势，而醇类和醛酮类香气成分的含量并未因空间位置的不同而呈现出显著的差异[5]，表明糟醅酯类和酸类的形成和窖泥相关，而醇类和醛酮类的形成和窖泥关系不大。另一项研究表明窖泥中有机酸含量随窖池深度的加深而增加[6]，结果和糟醅中的有机酸含量变化一致。对五粮液窖池中的窖泥和对应位置的出窖糟醅的香气成分进行定性定量分析，发现窖泥和对应的窖边糟中的优势香气成分高度相似，进一步研究发现窖边糟中香气成分含量显著高于中心糟，表明窖泥和窖边糟在白酒发酵过程中一直在相互作用，发生物质交换。

窖泥微生物中的优势菌群是维持窖泥微生物网络结构稳定、维持微生态平衡、进行微生物代谢和种间物质交换、产生重要香气物质及其前体的主要微生物，并随着空间深度的增加，不断反应和富集，形成与窖泥的空间位置密切相关的规律，进而影响着酒体中的呈香物质组成。窖泥中微生物群落的组成受空间分布的影响，与窖泥所在的空间位置密切相关。利用16S rRNA基因克隆文库技术分析了窖池上、中、底三个不同部位的窖泥中原核微生物，发现上层窖泥中微生物群落组成较为单一，中层和底层较为丰富[7]；高通量测序技术解析窖泥中微生物群落的多样性，同样得出在上部和下部的窖泥之间微生物群落组成存在差异的结论[8]。这些研究均表明，窖池中不同部位的窖泥中微生物组成的差异和窖泥中挥发性成分含量，尤其是与酸类和酯类的含量特别相关。

本研究以浓香型白酒窖池中不同空间位置的窖泥为研究对象，采取固相微萃取偶联气质联用法，分析不同空间位置的窖泥中挥发性物质的种类及分布规律，利用高通量测序技术分析不同空间位置窖泥中的微生物种类和分布规律，探索微生物分布规律和挥发性风味物质的分布规律之间的关系，以期为研究不同空间位置的窖泥对酒质的影响提供数据支持，为浓香型白酒的实际生产提供指导。

1 材料与方法

1.1 材料与仪器

窖泥 取自安徽文王酒业老车间窖池3-122，分别取窖池的上部、中部、下部、池底的窖泥，每个部位按照前后左右取四个长约5 cm，直径1.5 cm的圆柱形样品，将四个样品在低温下均匀混合作为一个试样。

手动SPME进样器、50/30 μm DVB/CAR/PDMS Stableflex固相微萃取头 默克Supelco公司；Agilent 7890-5975C气质联用仪 美国安捷伦科技公司。

1.2 实验方法

1.2.1 固相微萃取操作条件 先将萃取头置于气相质谱的进样口处250 ℃老化30 min[9]，取窖泥样品5 g置于20 mL顶空采样瓶中，45 ℃恒温水浴中加热平衡30 min，然后顶空吸附萃取60 min，吸附结束后，将萃取头插入GC-MS进样口250 ℃解吸3 min。每个样品重复取样三次分析。

1.2.2 气质联用分析条件 色谱柱：CP-Wax 57 CB（60 m×0.25 mm×0.4 μm）；载气：He；恒压：33 cm/s；进样口温度：250 ℃；升温程序：初始温度50 ℃，保持1 min，然后以3 ℃/min 升温速率升至180 ℃，保持1 min；再以8 ℃/min上升至225 ℃，保持3 min。

质谱电离方式：EI；电离能量：70 eV；传输线温度：280 ℃；离子源温度：220 ℃；扫描方式：全扫描（full）；质量扫描范围：m/z 33~450 amu；发射电流：100 μA。

1.2.3 窖泥高通量测序分析 将混合好的窖泥样品取样5 g，样品送苏州GENEWIZ生物科技有限公司做高通量测序，分析样品中含有的细菌的种类和数量，所得数据通过软件进行分析。

使用金唯智生物科技公司提供的Magen Hipure Soil DNA Kit试剂盒从窖泥样本中提取DNA，使用Qubit dsDNA HS Assay Kit试剂盒检测DNA浓度。PCR扩增及文库构建，以20~30 ng DNA为模板，使用金唯智生物科技公司设计的一系列PCR引物扩增原核生物16S rDNA上包括V3及V4的2个高度可变区。采用包含"CCTACGGRRBGCAS CAGKVRVGAAT”序列的上游引物和包含"GGA CTACNVGGGTWTCTAATCC"序列的下游引物扩增V3和V4区。另外，通过PCR向16S rDNA的PCR产物末端加上带有Index的接头，以便进行NGS测序。

PCR反应参数：预变性参数94 ℃ 3 min；变性参数94 ℃ 5 s，退火参数57 ℃ 90 s，延伸72 ℃ 10 s，共进行24个循环；终延伸参数72 ℃ 5 min。PCR扩增结果鉴定：PCR产物用1.5%琼脂糖凝胶电泳检测（表1）。

表1 PCR扩增条件Table 1 PCR amplification conditions

上机测序，通过Qubit3.0 Fluorometer（Invitrogen, Carlsbad, CA）检测文库浓度。将文库定量到10 nmol/L，按Illumina MiSeq（Illumina, San Diego,CA, USA）仪器说明进行PE250/FE300双端测序，由MiSeq自带的MiAeq Control Software（MCS）读取序列信息。

1.3 数据处理

本文所用数据平行及重复次数均为3次，绘图软件使用R语言和GraphPad Prism 8。

挥发性物质的定性与定量分析利用随机携带工作站，在NIST08和AMDIS标准谱库检索分析各组分质谱数据，确定分离得到的每个峰对应的风味物质，利用面积归一化法计算各组分的相对含量，利用面积加和比较法计算不同空间位置窖泥中的总酸、总酯、总醇、总醛酮及总挥发性物质的相对含量。

使用Bcl2fastq（v2.17.1.14）软件对测序结果原始图像数据进行基本调用，使用CutAdapt（v1.9.1）、Vsearch（1.9.6）和Qiime（1.9.1）软件对测序数据进行质量优化，去除其中的嵌合体序列，得到有效数据。对序列进行归类操作，使用Vsearch（1.9.6）按照97%相似性对unique序列进行OUT聚类，在聚类过程中进一步去除嵌合体序列，将优化后的序列与OUT代表序列进行比对，与OUT代表序列相似性在97%以上的序列为同一OUT[10]，统计生成OUT丰富表。根据OUT聚类分析结果，统计不同样本共有和特有的OUT，使用R语言绘制韦恩图。使用RDP classifier贝叶斯算法对97%相似水平的OUT代表序列进行分类学分析[11]，得到每个样本在各个水平的群落组成。使用Qiime（1.9.1）软件对有效序列进行随机抽样进行OUT分析，分别计算各α-多样性指数，估计环境群落的物种丰度和多样性。

2 结果与讨论

2.1 窖池中不同空间位置窖泥中挥发性物质的分析

对不同空间位置处的窖泥样品进行分析，总离子流图如图1所示，通过检索谱库NIST08和AMDIS，分析确定窖泥中挥发性物质成分。不同空间位置窖泥中各类挥发性物质成分及相对含量见表2。

表2 挥发性物质成分表及相对含量（%）Table 2 Volatile substance composition table and relative content (%)

图1 窖泥中挥发性物质的总离子流图Fig.1 Total ion chromatogram of volatile components in pit mud

由图2和表3得到的分析结果可知，窖泥中含有丰富的挥发性风味物质，但是不同空间位置的窖泥中这些挥发性物质的分布差异明显。初步鉴定出池壁上部窖泥样品中主要挥发性物质组分有25种，其中醛酮类4种（49.48%）、烃类10种（30.47%）、酚类2种（5.21%）、酯类5种（7.68%）、酸类1种（4.08%）、醇类3种（3.10%），其中最主要的挥发性物质是苯甲醛（40.03%）和5-（1,5-二甲基-4-己烯基）-2-甲基-，[S-（R*，S*）]-1,3-环己二烯（15.54%）。池壁中部窖泥样品中的主要挥发性物质组分有15种，其中酯类3种（67.13%）、醇类3种（19.08%）、酸类2种（5.88%）、酚类2种（4.26%）、醛酮类4种（3.46%）、烃类1种（0.17%），其中最主要的挥发性风味物质是己酸己酯（60.98%），其次是1-己醇（11.43%）。池壁下部窖泥样品中的主要挥发性物质组分有25种其中酯类14种（66.56%）、酸类6种（31.80%）、醇类3种（1.02%）、酚类2种（0.63%），其中最主要的挥发性物质是己酸乙酯（19.94%）、己酸（17.79%）、辛酸乙酯（15.52%）、己酸己酯（14.58%）；池底窖泥样品中的主要挥发性物质组分有31种，其中酯类18种（76.35%），酸类5种（19.37%）、醇类4种（3.42%）、酚类2种（0.57%）、醛酮类1种（0.26%）、烃类1种（0.04%），其中最主要的挥发性物质是己酸己酯（28.62%）、辛酸乙酯（14.98%）、己酸乙酯（13.02%）、己酸（12.26%）。

图2 不同空间位置窖泥中各类挥发性物质的相对含量Fig.2 The relative percentages and types of various volatile substances in pit mud in different spatial locations

2.1.1 酯类物质 由图2-a可知，窖泥池壁上部含有的酯类相对含量较少，中下底部远远高于池壁上部，且图3中池壁下部和底部酯类的种类数也较多，池底最多为18种。池壁中部虽然含有较多的酯类物质，但种类数只有3种。窖泥中所含酯类的峰面积总量随窖池深度的增加呈现大幅递增趋势，从池壁上部峰面积1.66×108增加到池底部276.63×108。表2中，己酸己酯、辛酸乙酯、苯丙酸乙酯为四组窖泥样品中共有成分，其中，己酸己酯由池壁中部位置处开始大量生成，成为池壁中部的主要酯类，占比60.98%，随着窖池深度的增加，酯类物质种类逐渐丰富，己酸己酯在酯类中的占比降低，池壁下部14.58%、池底部28.62%，但生成己酸己酯的量在增加。原因可能是在池壁上部，未能构成良好的厌氧环境，厌氧菌丰度不够未能成为优势菌属，不能将营养物质完全分解转化，产生的酸类和醇类等酯类的前体物质较少，故生成的酯类物质较少。

图3 不同空间位置窖泥中各类挥发性物质的种类数Fig.3 The relative types of various volatile substances in pit mud in different spatial locations

2.1.2 酸类物质 由图2-b和图3可知，窖泥池壁下部所含酸类物质的种类最多为6种，相对含量与其他位置相比最高，池底部次之，池壁上部和池壁中部所含酸类物质种类数和相对含量都较少。窖泥中所含酸类物质的峰面积总量随窖池深度的增加而增加，下底部增长较大，由池壁上部0.88×108增加到中部4.48×108再到池壁下部65.67×108和池底70.14×108。表2显示己酸均在四组窖泥样品中检测到，这与己酸己酯的分布一致，且己酸的含量呈现大幅增加趋势，这可能与产生己酸的菌群为厌氧菌群有关。

2.1.3 醇类物质 在图3中，窖泥各空间位置所含醇类物质种类数无明显差异，其中，池壁上、中、下部均含有3种，池底含有4种醇类物质。图2-c中池壁中部所含醇类物质的相对含量远远高于其它位置处，为19.08%。由表2可知，1-辛醇和1-己醇在四组窖泥样品中均可被检测到， 池底的醇类物质含量较池壁下部多，池壁中部最多。可能是因为在池壁中部，因营养丰富、环境适宜，大量微生物生长繁殖，丰度提高，对营养物质的分解较为充分，醛酮类和烃类物质大幅减少，醇类大量生成。

2.1.4 酚类物质 在图2-d中，窖池不同深度窖泥所含酚类物质均为2种，为苯酚和4-甲基苯酚。池壁上、中部酚类物质的相对含量大于池壁下部和池底。池壁上部的4-甲基苯酚含量高于其余部位，4-甲基苯酚可能与厌氧程度相关[12]。

2.1.5 醛酮类物质 图2-e和图3中，池壁上部所含醛酮类物质相对含量远远高于其它位置处窖泥，达到了49.48%，池壁中部含有四种醛酮类物质，相对含量为3.46%；池壁底部仅有一种，相对含量为0.26%；池壁下部则无醛酮类物质。浓香型白酒生产原料主要为粮食，其细胞壁中含有木质素。窖泥中的微生物将木质素降解，氧化生成苯酚类物质和芳香醛类物质[13]。故而丁香醛等物质多在池壁上部和中部被检测到。

2.1.6 烃类物质 由图2-f和图3可知，烃类物质较多的存在于窖泥池壁上部，10种相对含量为30.47%；池壁的中部和底部都只含有一种，相对含量各自为0.17%和0.04%；池壁下部未检测出烃类物质。酿酒生产过程中，酿酒的原料主要与池壁上部接触，木质素含碳量较高，被降解形成烃类物质，未能构成良好的厌氧环境，厌氧菌丰度不够未能成为优势菌属，未能将烃类进行进一步的转化，导致烃类物质较多地在上部被检测到。

窖池不同空间位置处窖泥含有的挥发性物质主要有酯类化合物、醇类化合物、酸类化合物、酚类化合物、醛酮类化合物、烯烃类化合物，这些物质与白酒中呈香化合物基本一致[14]。其中，酯类、酸类和醇类的种类数较多，且含量呈现出随着窖池深度的增加逐步增加的趋势。这可能是因为窖池中糟醅的水分由上自下逐渐浸落，部分微生物逐渐在窖池的下底部富集，因此窖池下底部的原核微生物群落与中上层存在明显差异的同时也具有一定的相似性[15]，也进一步证明了白酒中风味物质的产生位置主要在窖池的中部以下。不同空间位置的窖泥中优势菌群构成存在差异，影响着代谢反应循环，进而影响着酒体中的呈香物质组成，例如己酸菌发酵产生己酸、丁酸和少量的乙酸[16]；有机酸代谢过程中会产生氢，存在着微生物代谢控制的底物抑制现象，窖泥中的甲烷菌和硝酸盐还原菌均具有解除产酸菌的氢抑制现象[17]，促使有机酸的产生。

综上所述，窖泥中含有的烯烃类物质和醛酮类物质较多分布在池壁的上部，在池壁中下部和池底窖泥中存在较多的挥发性物质主要是酯类、醇类和酸类。不同位置主要挥发性成分存在差异是由于窖池中不同部位的厌氧程度决定的[18]。窖池的中部及以下厌氧性明显强于窖池上部，利于窖泥中厌氧菌的生长和代谢，对有机质的分解代谢较彻底[19]，利用营养物质代谢生酸并酯化，因此烯烃类物质较少，而醇、酸类物质较多，进一步反应生成酯类物质也较多[20]。在池壁下部和池底，由于厌氧环境良好，厌氧菌属和耐高酸等恶劣环境的菌属如梭菌属等成为优势菌属，能将纤维素、半纤维素等降解为含碳物质[21]，还能为酒香中重要的风味物质己酸乙酯的合成提供前体物质，对窖池中微生物群落的结构组成拥有调控能力，故而在这一空间位置能充分生成酯类和酸类等呈香性物质。说明长期低氧环境更有利于窖泥中的微生态系统的完善[22]和有机质的分解。

2.2 窖池中不同空间位置窖泥中的微生物群落分布分析

2.2.1 窖泥中微生物群落的α-多样性 群落生态学中，α-多样性通过一系列统计学指数的分析来估计环境群落的物种丰度和多样性[23]。计算菌群丰度的指数有ACE和Chao指数，ACE是用来估计群落中含有OUT数目的指数，Chao是用Chao1算法估计样本中所含OUT数目的指数[24]，ACE和Chao指数越大，说明群落中含有的OUT数目越多，群落的丰富度越大。计算菌群多样性的指数有Shannon和Simpson指数，Shannon用来估算样本中微生物多样性[25]，值越大说明群落多样性越高，Simpson在生态学中反映优势种在群落中的地位和作用。通过对原始序列进行质量控制和优化，共得到了65688条有效序列，平均长度为449.72 bp。每个窖泥样品中的序列数量从126576~163696个不等。且池壁中部窖泥中的序列数量最多，池底窖泥中的序列数量最少。

Observed OUT表示含有的物种数目，由表3可知，池底部和池壁中部的窖泥物种丰富度较高，池壁下部稍少于池壁中部，池壁上部的窖泥中物种丰富度最低。Shannon指数从2.130增加到4.970，Simpson指数从0.557增加到0.932，均随着窖泥深度的增加而增加，表明不同空间位置窖泥的物种多样性随窖泥深度的增加呈现增加的趋势，池壁中下部和底部的Simpson指数增长较为平缓，表明中下底部窖泥物种多样性增加的同时趋于稳定，且池壁中下部和池底的窖泥中微生物多样性显著高于池壁上部。ACE指数由70.026增加到112.271，随后有小幅度降低，又由99.510增加到114.523。Chao1指数由67.51增加到112.00，小幅降低至98.50，再增加到114.50。ACE指数和Chao1指数表明窖泥微生物物种丰富度呈现先上升后下降再上升的趋势，池底的物种丰富度最高。Chao1、Shannon、Simpson、ACE指数均表明，在整个窖池中，池底的窖泥中微生物物种多样性和物种丰富度最高。

表3 不同空间位置窖泥中微生物群落的α-多样性Table 3 α-Diversity of microbial communities in pit mud in different spatial locations

2.2.2 不同空间位置窖泥中微生物群落门、纲水平的变化 根据图4 OTU注释，共鉴定出10个门，其中细菌门9个，古细菌门1个。共发现8个优势菌群（平均丰度>0.1%），占每个样品细菌总数的98.86%~99.93%，其中厚壁菌门（Firmicutes）为绝对优势菌群。窖池上部优势菌群由高到低的顺序为：厚壁菌门（Firmicutes）>广古菌门（Euryarchaeota）> 拟杆菌门（Bacteroidetes）>变形菌门（Proteobacteria）；池壁中部和下部优势菌群由高到低的顺序较为相似：厚壁菌门（Firmicutes）>拟杆菌门（Bacteroidetes）>广古菌门（Euryarchaeota）>互养菌门（Synergistetes）>软壁菌门（Tenericutes）>放线菌门（Actinobacteria）；窖池底部由高到低分别为：厚壁菌门（Firmicutes）>广古菌门（Euryarchaeota）>放线菌门（Actinobacteria）>拟杆菌门（Bacteroidetes）>暗黑菌门[26]（Atribacteria）>互养菌门（Synergistetes）>软壁菌门（Tenericutes）。随着窖池深度的增加，窖泥中的优势菌群逐渐由4个类群向多个菌群转变，最后又回到单一的绝对优势菌群，厚壁菌门数量先减少后增加，其他优势菌群均有不同程度的先增加后减少最后在池壁底处最少。

图4 窖池不同空间位置微生物群落门水平（a）和纲水平（b）的分布Fig.4 Levels of phylum microflora（a）and class microflora（b）in pit mud samples at different spatial locations

在纲水平上，从所有窖泥样品中检测到15类原核微生物。共发现6个优势菌群（平均丰度>1%），从高到低的顺序为：梭菌纲（Clostridia）>杆菌纲（Bacilli）>拟杆菌纲（Bacteroidia）>甲烷 微菌纲（Methanomicrobia）>互养菌纲（Synergistia）>甲烷杆菌纲（Methanobacteria）。池壁中部和池底窖泥中都有丰富的细菌分布，池壁下部和池壁上部次之。梭菌纲的含量随着窖池深度的增加而逐渐增加，在池壁下部和池底窖泥中达到最高；杆菌纲的含量则随着窖池深度的增加而逐渐减少，在池壁下部和池底的窖泥中接近于0；拟杆菌纲和互养菌纲的含量随着窖池深度的增加逐渐增加，并在池底部时降到最低；随着窖池深度的增加，甲烷微菌纲和甲烷杆菌纲的数量呈现先降低再增加再降低的规律。窖泥微生物的群落组成因处于窖池的不同空间位置而具有明显的空间差异性。

2.2.3 不同空间位置窖泥中微生物群落优势菌属水平的变化 在属水平上，从4个不同空间位置的窖泥样品中检测到64个可识别属和未分类属。对每份窖泥样品中1%以上核心属进行统计分析，共有21个核心属，含量分别为98.43%（池壁上部）、85.01%（池壁中部）、85.65%（池壁下部）、77.73%（池底），表明池壁底是窖泥中微生物最复杂的位置。

本研究选取每个样品的前30个属进行热谱分析，从较小的分类水平上揭示窖泥中微生物群落组成的差异。结果如图5所示，通过热图的颜色差异可直观显示物种的相似性和差异性。前30个属占比分别占每个样品的88.79%~99.2%。根据图6可知，池壁上和池壁中、池壁下和池壁底明显地被分为两个聚类，下底部窖泥中微生物结构与中上部存在部分相似性又存在较大的差异性。在窖池上部，芽孢八叠球菌属（Sporosarcina）的含量最高，随着窖池空间位置的增加逐渐降低，在窖池的下底部含量几近为零，符合多分布于泥土表面的习性。Sedimentibacter、Unclassified_Unclassified、f_Syntrophomonadaceae_Unclassified、甲烷囊菌属（Methanoculleus）、Caproiciproducens等属的变化规律大体相同，都随着窖池空间位置的增加逐渐增加；部分菌如互营单胞菌属（Syntrophomonas）、Petrimonas、甲烷八叠球菌属（Methanosarcina）、Proteiniphilum等则呈现从池壁上至池壁下先增加后在池壁底部降低的规律。说明窖池中微生物种类的空间分布与窖池的空间位置有关。

图5 每个样品中TOP30微生物属水平分布热图Fig.5 Heat map of the horizontal distribution of TOP30 microorganisms in each sample

图6 不同空间位置窖泥基于OUT水平的韦恩图Fig.6 Venn diagram of OTU level in pit mud samples at different spatial locations

梭菌属（Clostridium）、互营单胞菌属（Syntrophomonas）、Sedimentibacter、泰氏菌属（Tissierella）、Garciella等梭菌科属是生成丁酸、己酸等酯类前体物质的重要菌群，其他优势菌属Petrimonas、Proteiniphilum等微生物的代谢产生物也构成白酒中的重要风味物质或风味物质前体。本研究中检出的优势古菌属中，甲烷菌属（Methanobacterium）、甲烷短杆菌属（Methanobrevibacter）和甲烷囊菌属（Methanoculleus）属氢营养型甲烷菌[27]，而氢营养型产甲烷古菌可以H2为代谢底物促进己酸的合成[15]，甲烷八叠球菌属（Methanosarcina）可同时利用乙酸和H2，这些甲烷菌可以通过种间作用来实现己酸乙酯产量的提升。这些优势菌群都是维持窖泥微生物网络结构稳定、维持微生态平衡、进行微生物代谢和种间物质交换、产生重要香气物质及其前体的主要微生物，并随着空间深度的增加，不断反应和富集，形成与窖泥的空间位置密切相关的规律。

韦恩图中不同颜色的圈表示不同的样本分组，图中数字代表每个样本分组特有或共有的OUT数目，最中间数字表明所有样本分组共有的OTU数目。根据OUT聚类分析结果，由图6可知，共有属部分数量为51种，占43.59%，池壁上含有特种属，数量和占比分别为2和1.71%，池壁中下部和池底窖泥中微生物群落较池壁上部具有更大的相似性。

2.2.4 多样本比较分析β-多样性值为两个样本间的相异系数，反映不同样本间的多样性的差异。

样本点之间距离表明了样本中功能分类分布的相似性，距离越近，相似度越高。图7-a中，在主成分1（PC1）的条件下，样本中池壁上部和池壁中部极为接近，相似度较高，与池壁下部和池底部存在差异。在主成分2（PC2）的条件下，样本中的池壁下部和池底距离极为接近，与池壁中部较为相近，相似度较高，与池壁上部距离较远，存在较大差异性。非度量多维尺度法是根据样本中包含的物种信息，以点的形式反应在多维空间上，通过点与点之间的距离来反映不同样本间的差异程度。图7-b中，每一个点代表一个样本，表明窖泥池壁上部的菌群组成与窖泥池壁中下底部的菌群组成存在明显差异，而窖泥池壁中部和池壁下部菌群结构组成类似。β-多样性分析表明，窖泥中的群落结构随着窖泥空间位置的加深与环境因子相互作用而逐渐趋于稳定，空间位置在一定程度上影响了窖泥中原核微生物群落的α-和β-多样性，这主要是由于不同位置窖泥的理化性质不同所致，这是影响微生物群落多样性变化的根本因素[28]。

图7 PCA图（a）和非度量多维标度分析（b）Fig.7 PCA diagram(a) and non-metric multidimensional scaling analysis(b)

2.3 窖池中不同空间分布窖泥中挥发性物质分布和微生物群落分布之间的关系

本研究发现窖池下部和池底窖泥中的原核菌群的丰度和多样性[29]要高于窖池上部和中部的窖泥，这与挥发性物质的种类分布规律相吻合，与窖池下部和底部长期良好的厌氧环境有关，与高酸环境也有一定关联。生孢产氢菌（Hydrogenispora）在窖池的各个空间位置均有分布，主要产物是乙酸、乙醇和氢等[30]。在池壁上部，芽孢八叠球菌属（Sporosarcina）为上层的优势菌属，能发酵产生二氧化碳和乙酸等，但因厌氧环境较差，厌氧菌丰度不够，不能将营养物质完全分解，导致烯烃类和醛类物质较多；池壁中层窖泥拥有丰富的营养物质，环境因子较为适宜，原核微生物群落的多样性也反映着丰富的风味物质组成，优势菌群如柯生菌属（Fastidiosipila）等梭菌属菌群丰度增加，白酒呈香性物质中重要的己酸乙酯及其前体物质等含量也开始增加；池壁下部和底部窖泥因长期处于良好的厌氧环境中，适宜在极端生长环境中的原核微生物充分代谢反应，甲烷细菌属（Methanobacterium）、甲烷短杆菌属（Methanobrevibacter）和甲烷囊菌属（Methanoculleus）属等甲烷菌，在促进己酸的合成的同时，还可以通过种间作用与己酸菌共生共养[31]，实现己酸和己酸乙酯产量的提升，在此空间位置丰度较高的微生物物种，对窖池中的微生物群落结构组成拥有调控能力，故而在这一空间位置能充分生成酯类和酸类等呈香性物质。随着窖泥空间位置的加深，越靠近下层厌氧环境越好，厌氧菌在下底部大量沉积繁殖，在发酵过程中产生的醇、醛、酸、酮、酯等代谢产物不仅为酒体的呈香物质，部分还能作为生长因子或代谢前体物质促进微生物的生长繁殖，为多样的风味物质组成提供了物质基础。

3 结论

随窖池深度的增加，挥发性物质总量也在逐渐增加,峰面积总量由池壁上部21.65×108增加到池底部362.31×108。其中酯类和酸类种类较多，占比达到73.01%~98.36%。窖泥样品中池壁下部和池底部窖泥中的微生物菌群丰度和多样性高于池壁上部和中部，这与挥发性物质的种类和含量分布规律相吻合。微生物门水平上优势菌群占细菌总数的98.86%~99.93%，厚壁菌门为绝对优势菌群。梭菌属和甲烷菌属等优势菌属是进行微生物代谢和种间物质交换、产生重要香气物质及其前体的主要微生物，并随着空间深度的增加，不断反应和富集，形成与窖泥的空间位置密切相关的规律。空间差异性是导致窖泥不同空间位置原核微生物菌群和挥发性物质成分存在差异性的重要原因。

窖泥生态系统中蕴藏着一个与窖泥密切相关的复杂原核生物群落[32−33]，微生物群落与风味物质之间并非简单的一一对应的关系。本项研究联系了窖泥空间特征对原核生物群落的影响及挥发性物质的空间分布规律，研究结果对于不同空间位置窖泥中微生物群落的分布规律及其间的原核生物群之间的相互作用起到一定的参考，对酒质的影响提供数据支持。