Study on the Correlation between Microbial Community Succession and Physicochemical Factors during the Storage of Maotai-flavor Liquor in Daqu at High Temperature
-
摘要: 目的：本研究分析了酱香型白酒高温大曲在储存过程中（0~150 d）的理化特性和微生物的变化过程，以期了解大曲的老化成熟规律。方法：该研究基于传统可培养方法和高通量测序技术分析大曲在储存过程中微生物群落结构演替规律，结合相关性分析理化因子和微生物之间的关系。结果：储存过程中的可培养微生物的数量呈现先增加后下降的趋势；温度和糖化力指标呈先上升后下降的趋势，水分、酸度随着存储时间的增加呈现下降的趋势。储存过程中的主要优势细菌属为Kroppenstedtia、Lactobacillus、Staphylococcus、Enterobacter；主要优势真菌属为Thermomyces、Trichomonascus、Byssochlamys。微生物群落与理化指标间存在复杂的作用关系，Bacillus、Kroppenstedtia、Aspergillus、Trichomonascus与温度呈正相关关系；Byssochlamys、Dipodascus与温度、水分、酸度、液化力呈正相关关系；Thermomyces与温度、水分、酸度呈负相关关系。结论：本研究揭示了酱香型高温大曲在储存过程中微生物群落组成情况及理化因子间的相互关系，为进一步阐明酱香型白酒高温大曲后熟机理提供了理论参考。Abstract: Objective: Through analyzing the physicochemical properties and microbial changes of high temperature Daqu of Maotai-flavored during the storage (0~150 days), the paper aimed at comprehending the aging and maturity law of Daqu. Methods: Based on the traditional cultivable method and high-throughput sequencing technology, the research analyzed the succession law of microbial community structure during the storage of Daqu, and the relationships between physicochemical factors and microorganisms by combining the correlation. Results: The number of cultivable microorganisms showed a trend of first increased and then declined with the increased storage time during the storage. The temperature and saccharification index showed a trend of first increased and then declined, moisture and acidity showed a declined trend with increasing storage time. Meanwhile, the main leading bacteria genus during the storage were Kroppenstedtia, Lactobacillus, Staphylococcus and Enterobacter respectively, the main leading fungus genus were Thermomyces, Trichomonascus and Byssochlamys respectively. There was a complicated relationship between microbial community and physicochemical indicators, of which Bacillus, Kroppenstedtia, Aspergillus and Trichomonascus were positively correlated with temperature, Byssochlamys and Dipodascus were positively correlated with temperature, moisture, acidity and liquefaction power. Thermocycles was negatively correlated with temperature, moisture, and acidity. Conclusion: The research revealed the microbial community composition status and the relationships between physicochemical factors of high-temperature Daqu of Maotai-flavored during the storage, which provided a theoretical reference for further introducing the post ripeness mechanism of high-temperature Daqu of Maotai-flavored.
-
“酒乃曲之骨”，酱香型白酒发酵过程中的微生物主要源于高温大曲[1]，大曲是一种富含多酶多菌的微生物发酵剂，具有“糖化、发酵、生香”等功能，是白酒发酵过程的重要物质保障，其品质的好坏决定了酱香型白酒的产量和质量[2-5]。
酱香型高温大曲采用生料制曲、添加母曲自然接种、经40 d发酵逐渐形成了以细菌、霉菌和酵母菌为主的群落结构[6]。酱香型高温大曲制作完成后，需进行拆曲入库储存，使大曲进入“后熟”阶段。合理的储存期对酱香型白酒酿造有至关重要的作用。在大曲储存的理化特性研究方面，张丽等[7]研究结果表明高温大曲储存5~6个月时，液化力、酯化力等理化特性趋于稳定；刘建文等[8]发现储存3~4个月的大曲，其功能酶活力均达到最好状态。在微生物多样性研究方面，粱丽文等[9]研究结果表明高温大曲存储过程中主要优势细菌属为Bacillus，真菌属为Pichia；梁晨等[10]研究发现大曲贮存4~6个月，更有利于白酒的酿造，且Thermoactinomyces、Bacillus等是大曲储存过程中的优势细菌属，未对真菌多样性进行相关研究。陈晓茹等[11]研究表明强化曲在贮存过程中的优势菌属是Bacillus、常规曲的优势菌属为Weissella。目前，对酱香型白酒酿造过程中的微生态结构分析方法主要是传统可培养方法，但该方法耗时长，准确性较低，且具有一定的局限性，远远不能满足对微生物群落结构和多样性的认识，但是传统的可培养方法在微生物检测上的作用是无法替代的。然而随着分子生物学技术的快速发展，能快速精确地揭示微生物种类的多样性，鉴定样本中微生物群落组成的能力，可检测到未发现的微生物，更加全面地了解系统中的微生物，客观地认识微生物的群落结构及其演替变化。对于高温大曲“后熟”质量的评价主要依靠感官和理化指标，且在储存过程其影响因素有待于进一步阐明。在实际生产中，大曲是否成熟主要依靠生产经验判断，存在着盲目性和不确定性。
因此，本文以“后熟”过程的酱香型白酒高温大曲为研究对象，同时采用传统可培养方法和高通量测序技术解析细菌和真菌群落组成，并结合高温大曲中各理化因子动态变化规律，分析微生物群落和理化指标间的相互作用关系，以期为阐明酱香型白酒高温大曲“后熟”过程各因素变化机理提供理论参考。
1. 材料与方法
1.1 材料与仪器
高温大曲样品采集自GZXJ制曲车间不同储存时间成品曲。孟加拉红培养基、营养琼脂 北京奥博星生物技术有限责任公司；WL培养基 上海博微生物科技有限公司；DNA提取试剂盒 美国OMEGA BioTek公司；聚合酶链式反应引物 上海翌圣生物科技股份有限公司。
SW-CJ-2F无菌操作台 苏州净化设备有限公司；ZQPW-70全温振荡培养箱 天津市莱玻特瑞仪器设备有限公司；KG-AP32L自动蒸汽灭菌器 日本ALP公司；GL-88B漩涡混合器 海门市其林贝尔仪器制造有限公司；TND03-H-H混匀型干式恒温器 深圳托能达科技有限公司；ETC811 PCR仪 北京东胜创新生物科技有限公司；DYCZ-21电泳仪 北京市六一仪器厂；FR-1000凝胶成像系统 上海复日科技有限公司。
1.2 实验方法
1.2.1 样品采集
选取GZXJ公司制曲车间某一储存曲库，在储存过程中进行定点跟踪取样，每隔1月取1次，共计6个样品。样品磨碎后置于无菌密封袋中，于−20 ℃下低温保藏，按取样时间依次标记为GX-1（0 d）、GX-2（30 d）、GX-3（60 d）、GX-4（90 d）、GX-5（120 d）、GX-6（150 d）。
1.2.2 理化指标的测定
温度、水分、酸度、糖化力、液化力的测定具体参照QB/T 4257-2011《酿酒大曲通用分析方法》[12]。
1.2.3 可培养微生物的分离保藏、计数及鉴定
a. 分离及保藏：称取10 g样品于90 mL无菌水中，于25 ℃、180 r/min条件下振荡30 min，吸取1 mL悬浮液稀释备用，分别采用营养琼脂、WL培养基和孟加拉红培养基分离筛选细菌、酵母菌和霉菌，观察菌落在培养基上的形态。细菌和酵母菌采用−80 ℃甘油保藏，丝状真菌采用4 ℃斜面保藏，每株菌株平行保藏3管。
b. 计数：采用稀释涂布的方法，分别采用营养琼脂培养基于37 ℃条件下培养细菌；WL培养基和孟加拉红培养基于28 ℃条件下培养酵母菌和丝状真菌培养结束后，用计数器对肉眼观察的单菌落计数，记录稀释倍数和相应的菌落数量。
c. 鉴定：根据菌落形态特征进行初步分类后，每种类别挑选代表性菌株1~3株[13]。分别采用Ezup柱式基因组DNA试剂盒对菌株DNA进行提取，具体操作流程见试剂盒说明书。参照吴成等[14]方法，细菌采用通用引物27F（AGTTTGATCMTGGCTCAG）和1492R（GGTTACCTTGTTACGACTT）对 16S rRNA 基 因 扩 增 测 序；酵母菌采用通用引物NL1（GCATATCAATAAGCGGAGGAAAAG）和NL4（GGTCCGTGTTTCAAGACGG）对26S rRNA 基因D1/D2 片段扩增测序；丝状真菌采用通用引物ITS1（TCCGTAGGTGAACCTGCGG）和ITS4（TCCTCCGCTTATTGATATGC）对内转录间隔区1和2片段扩增测序。
PCR反应体系为：模板DNA 0.5 μL，10×Buffer（with Mg2+）2.5 μL，dNTP2.5 μL，Taq酶0.2 μL，正向及反向引物各0.5 μL，使用ddH2O补齐体系至25 μL。PCR扩增条件为：94 ℃预变性4 min；94 ℃变性45 s；55 ℃退火45 s；72 ℃延伸1 min；30次循环；72 ℃终延伸10 min 降至4 ℃。扩增完成后，回收PCR产物并纯化。
1.2.4 大曲总DNA提取及扩增
总DNA提取：按照DNA试剂盒E.Z.N.A™ Mag-Bind Soil DNA Kit提取试剂盒方法提取DNA。细菌16S V3~V4区域片段，选择341F（5'-CCTACGGGNGGCWGCAG-3'）和805R（5'-GACTACHVGGGTATCTAATCC-3'）为扩增引物；真菌采用ITS3（GCATCGATGAAGAACGCAGC）和ITS4（TCCTCCGCTTATTGATATGC）为扩增引物。
PCR扩增条件为：95 ℃预变性3 min；94 ℃变性20 s；55 ℃退火20 s；72 ℃延伸30 s；25次循环；72 ℃终延伸5 min降至4 ℃。PCR反应体系为：2×Hieff® Robust PCR Master Mix 15 μL，正向及反向引物各1 μL，DNA模板20~30 ng，使用ddH2O补齐体系至30 μL。
1.2.5 高通量测序
采用Illumina Miseq测序平台，分别对细菌16S rRNA基因 V3~V4区、真菌ITS3~ITS4区进行高通量测序分析（由上海生物工程有限公司完成）。
1.3 数据处理
使用UNITE为参考数据库进行物种注释，基于Silva库对97%相似水平的OTU代表序列进行分类学注释分析，细菌和真菌分别采用RDP和UNITE数据库进行比对。采用Microsoft Office Excel 2016进行数据统计和处理；理化因子动态图、可培养微生物动态图和物种相对丰度柱状图采用Origin 2017软件进行处理绘制；优势微生物和理化因子之间相关性热图先利用SPSS 22.0软件计算Spearman相关性系数，再使用OmicShare Tools平台绘制组间动态相关性热图；动态丰度Venn图使用micShare Tools平台绘制；冗余分析（Redundancy analysis，RDA）采用Canoco5.0软件。
2. 结果与分析
2.1 储存过程中理化指标的动态变化
大曲中的水分、酸度、液化力、糖化力等是评价大曲品质好坏常用的理化指标。由图1可知，随着储存时间的变化，储存仓温度先上升后下降，曲坯的水分、酸度逐渐下降，糖化力呈现先上升后下降的趋势，液化力在（0.02±0.01）~（0.04±0.02）g/g·h之间波动。微生物可以利用大曲中淀粉、脂肪以及蛋白质降解成小分子进而代谢产酸，随着储存过程中曲坯水分逐渐挥发，微生物的生长繁殖受到影响，导致酸度下降。有研究表明，大曲的酸度主要来源于生酸微生物降解蛋白质、脂肪、淀粉等物质和有机酸代谢产酸，是微生物综合作用的结果[15~17]。大曲的糖化力受工艺控制参数的影响，并且与产品风格密切相关[18]。由图1可知，随着储存时间的变化，糖化力先升高后下降，在储存120 d（GX-5）时达最高，为248.97±0.02 mg/g·h。大曲的液化力反映大曲中液化型淀粉酶将淀粉水解为糊精的能力，主要是由霉菌和枯草芽孢杆菌产生[19]。
2.2 可培养微生物的动态变化
采用稀释涂布平板法对储存0~150 d的可培养微生物进行计数、分离，其结果分别表1、图2所示。
表 1 可培养微生物菌落总数变化Table 1. Changes in the total number of culturable microorganisms样品 细菌数
[（×104 CFU）·g−1]霉菌数
[（×104 CFU）·g−1]酵母菌数
[（×104 CFU）·g−1]GX-1 1.40±0.11f 1.40±0.26b 12.10±0.31a GX-2 5.50±0.30e 5.53±0.53a 3.14±0.19d GX-3 81.70±1.02a 0.19±0.04d 3.57±0.15c GX-4 59.20±0.58b 0.75±0.09c 4.32±0.23b GX-5 58.80±0.77c 0.53±0.03d 0.48±0.03e GX-6 47.80±0.42d 0.20±0.04d 0.21±0.04e 注：同列不同小写字母表示不同储存时间微生物数量具有显著性差异（P<0.05）。 ![]()
图 2 可培养方法分离筛选微生物形态特征Figure 2. Microbial characteristics revealed by culture-dependent approaches由表1可知，大曲初入储存仓时，细菌和霉菌均在104 CFU/g数量级，酵母菌在105 CFU/g数量级；大曲储存到60 d（GX-3）时，可培养微生物数量最高，90 d（GX-4）时可培养微生物的数量缓慢下降。分析可能原因是在大曲初入曲房时，储存仓温度较低，随着储存时间增加，曲房温度升高，此阶段大多数的芽孢细菌在此时形成芽孢，有利于微生物生长繁殖。后期随着储存时间的变化，曲坯的水分降低，储存仓温度下降，微生物数量缓慢下降。大曲在储存过程中受成品曲仓、曲坯内部温度、水分、酸度等参数的影响，同时也受曲坯内部微生物生态结构变化的调控，储存过程中存在着一个双向调控机制[20-21]。
本研究采用稀释涂布平板法分离筛选到15株细菌、6株酵母菌和12株霉菌，共保藏菌株99株。选择代表性菌株进行分子鉴定为B. licheniformis 、B. amyloliquifaciens、S. gallinarum三种细菌（a~c）；T. ciferrii、M. farinosa、S. cerevisiae三种酵母菌（d~f）；L. ramosa、A. sydowii、A. chevalieri（g~i）三种霉菌，其形态特征如图2所示。相关研究表明，高温大曲中分离出的优势菌种也具有一定的相似性，王晓丹等[22]研究表明茅台地区酱香型高温大曲中的主要优势细菌属为Bacillus；张立强等[23]从泸州老窖大曲中分离出高产纤溶酶的优势菌B. licheniformis：李子键等[24]从中高温大曲中分离出耐高温的微小根毛霉、布氏横梗霉，并且分别具有产糖化酶、蛋白酶活性功能。
2.3 大曲在不同储存时间免培养微生物的变化情况
2.3.1 Alpha多样性指数分析
样品测序完成后，得到的原始序列经去除嵌合体、过滤等处理后共得到498207条有效序列，其中细菌202228条、真菌295979条。细菌和真菌测序覆盖率都为0.99以上，表明测序深度足够，样品中几乎所有样本都被检出，测序结果可真实反映出酱香型高温大曲在储存过程中细菌和真菌群落多样性组成情况。对于细菌而言，Shannon指数呈先升后降，Simpson指数先降低后趋于平稳，表明大曲群落多样性在储存60 d（GX-3）、90 d（GX-4）、120 d（GX-5）时最高；对于真菌而言，Shannon指数呈缓慢下降趋势，Simpson指数呈逐渐升高趋势，表明随着大曲储存时间延长，真菌群落多样性逐渐降低（表2）。
表 2 不同储存时间真菌和细菌Alpha多样性Table 2. Alpha diversity of fungi and bacteria at different storage times测序片段 样品 序列数目 OTUs Shannon Chao Ace Simpson Coverage 16S GX-1 22416 50 1.19 65.00 81.72 0.40 0.9993 GX-2 36379 52 1.97 56.00 57.87 0.20 0.9998 GX-3 34416 63 2.41 105.75 135.57 0.13 0.9994 GX-4 31846 94 2.42 119.20 166.07 0.17 0.9991 GX-5 36316 80 2.34 93.33 93.12 0.16 0.9996 GX-6 40855 78 2.34 117.00 147.45 0.14 0.9993 ITS GX-1 46639 21 1.18 31.00 25.64 0.36 0.9999 GX-2 53990 26 1.11 35.33 56.37 0.42 0.9999 GX-3 42845 23 0.99 30.00 44.15 0.42 0.9998 GX-4 51947 44 0.93 122.00 58.68 0.57 0.9998 GX-5 54880 30 0.87 38.25 59.97 0.48 0.9998 GX-6 45678 71 0.48 76.14 74.46 0.76 0.9998 2.3.2 OTU分布的Venn图
Venn 图可以用来统计样本中共有和独有的OTU的数目，直观地展现样品OTU数目组成的相似性及重叠情况[25]。如图3（a）所示，样品GX-1（0 d）到GX-6（150 d）独有的OTU数目依次为7、8、6、23、8和14，其中GX-4（90 d）样本中独有的物种数所占比例高，那么样本之间的相似性差，该样本的独特性高；6个样品共有的 OTU为11个，其所代表的细菌物种长期存在于储藏期内。图3（b）所示，样品GX-1（0 d）到GX-6（150 d）独有的OTU数目依次为3、1、2、11、5和47，共有的OTU为7个，其中GX-6样本独有的物种数所占比例高，独特性高，GX-1~GX-5（0~120 d）与共有样本的数量接近，说明共有物种的占比高，因此样本的相似性就高。
![]()
图 3 不同样品的细菌（a）、真菌（b） OTU Venn图Figure 3. Venn profile of OTU of bacteria (a) and fungi (b) in different samples2.3.3 不同储存时间细菌、真菌群落动态变化
样品经高通量测序后，细菌共检出12个门、21个纲、43个目、80个科、119个属和179个种，其中Firmicutes为绝对优势细菌门，储存过程中平均占比71.69%；Kroppenstedtia（平均占比30.45%）、Lactobacillus（平均占比11.53%）、Staphylococcus（平均占比10.31%）为主要优势菌属。真菌共检5个门、17个纲、37个目、89个科、142个属和203个种，其中Ascomycota为绝对优势真菌门，平均占比98.11%，Thermomyces（43.87%）、Trichomonascus（29.77%）为主要优势真菌属。图4、图5分别为平均相对丰度至少在一个大曲样品中大于1%的菌群结构。
![]()
图 4 大曲在门水平上细菌（a）、真菌（b）组成情况Figure 4. Composition of bacteria (a) and fungi (b) in high-temperature Daqu at the phylum![]()
图 5 大曲在属水平上细菌（a）、真菌（b）组成情况Figure 5. Composition of bacteria (a) and fungi (b) in high-temperature Daqu at the genus由图5可知，Lactobacillus、Lactococcus、Saccharopolyspora、Virgibacillus细菌属和Thermomyces、Trichomonascus、Aspergillus真菌属物种多样性呈现先增后降，在储存到120 d（GX-5）以后，优势菌群的相对丰度含量逐渐下降。分析造成这种现象可能的原因是储存时间过长的曲块，曲虫增多，使大曲中的淀粉等物质流失，破坏大曲的微生物群落结构，因而导致微生物多样性降低[26]。
酱香制曲过程中的细菌主要来源于空气、场地、原辅料及人体，在酱香型白酒酿造过程中的主要作用是产香，在发酵体系中能产生大量酶类，是形成酱香型白酒风味前体物质的主要菌群[27-28]。相关研究表明，Bacillus是酱香型白酒发酵过程中主要的功能细菌种群，具有代谢水解酶、乙偶姻等风味物质功能，对酱香型白酒的风味形成有重要贡献[29-30]；Lactobacillus属于乳酸菌菌系，大多是兼性厌氧或厌氧菌，在发酵过程具有重要的生物学调控作用[31]。真菌方面，酿酒酵母是产生酒精的主要功能菌种，曲霉菌和根霉菌具有很强的酶活力，且绝大多数的真菌属于嗜热真菌属，可以在高温下进行生长代谢，产酶丰富，具有分解酒醅中大分子物质、产生白酒中风味化合物及其前提物质的作用[32]。
2.3.4 理化因子与微生物群落相关性分析
理化因子的变化对微生物的生长繁殖及其行为具有较强的影响，进而调控微生物群落的结构及其演替[21]。如图6（a）所示，RDA分析结果表明大曲在储存0 d（GX-1）时，大曲的后熟过程与温度、水分、液化力和酸度呈正相关关系，与糖化力呈负相关关系，表明该时期温度和水分适宜，比较适合微生物生长；储存30~60 d（GX-2~GX-3）时，大曲中的微生物群落结构不受水分、酸度等因素的影响，表明该时间段理化指标对微生物群落结构影响不大；储存90~150 d（GX-4~GX-6）时，温度、水分、酸度、液化力与储存过程呈负相关，表明了大曲在存储后期，曲仓温度降低，曲坯水分含量低，不适宜微生物代谢活动，因此，微生物多样性降低。
![]()
图 6 微生物群落结构与大曲理化因子相关性分析注：a：冗余分析，GX-1~GX-6分别代表大曲储存0~150 d，理化因子箭头间的夹角代表正、负相关性；b：优势微生物群落与环境因子的相关性Heatmap图，热图中的X轴和Y轴分别为理化因子和物种（属水平），通过计算获得R值和P值来考察两个（组）变量之间的相关程度。R值在图中以不同颜色展示，热图右侧色卡是不同R值的颜色分区。*：0.01<P≤0.05，**：0.001<P≤0.01。Figure 6. Correlation analysis between microble community structure and physical and chemical factors of Daqu为揭示理化因子对优势微生物群落的影响，选取相对丰度分别排名前10的优势细菌属、优势真菌属（相对丰度含量至少在一个大曲样品中大于1%的菌群结构）与储存过程中理化因子进行相关性分析。结果如图6（b）所示，对细菌而言，Bacillus与水分、酸度、糖化力、液化力呈负相关关系，与温度呈正相关关系，表明了该细菌属的生长代谢等活动受温度的影响，在大曲储存后期，曲仓温度降低，其相对丰度含量降低，高酸的环境也会抑制其生长。Staphylococcus与温度、水分呈正相关关系、与酸度、糖化力、液化力呈负相关关系；Kroppenstedtia与温度、水分、酸度、液化力呈正相关关系，与糖化力呈负相关关系。
对真菌而言，Aspergillus、Trichomonascus等菌属与水分、酸度、糖化力、液化力呈负相关关系，与温度呈正相关关系；Byssochlamys、Dipodascus与温度、水分、酸度、液化力呈正相关关系，与糖化力呈负相关关系；Thermomyces与温度、水分、酸度呈负相关关系；Rhizomucor与水分、温度、酸度、液化力呈正相关关系。研究表明Thermomyces属于嗜热真菌属，具有促进蛋白质降解的作用，使得微生物能充分利用酿酒原料进行生长代谢[33]。Rhizomucor具有糖化发酵能力，能形成大量乳酸[34]。Aspergillus、Thermomyces、Bacillus、Issatchenkia、Lactococcus等是大曲储存过程中的主导菌群，这些菌群的存在有利于提升大曲的发酵性能[35]。
3. 结论
本文同时采用了传统可培养方法和高通量测序技术对酱香型白酒高温大曲在储存（后熟）过程中微生物群落组成情况及理化因子间的相互关系进行了解析。可培养结果表明，随着储存时间的增加，菌落总数呈现先上升后下降趋势；高通量测序结果表明，在储存过程中Kroppenstedtia、Lactobacillus、Staphylococcus、Thermomyces、Trichomonascus等为主要优势菌属；RDA分析结果表明大曲在储存0 d时，大曲的后熟过程与温度、水分、液化力和酸度呈正相关关系，与糖化力呈负相关关系，表明该时期温度和水分适宜，比较适合微生物生长；储存30~60 d时，大曲中的微生物群落结构不受水分、酸度等因素的影响，表明该时间段理化指标对微生物群落结构影响不大；储存90~150 d时，温度、水分、酸度、液化力与储存过程呈负相关，表明了大曲在存储后期，曲仓温度降低，曲坯水分含量低，不适宜微生物代谢活动，因此，微生物多样性降低。目前，对于酱香型高温大曲微生物多样性相关研究层出不穷，在优势微生物研究方面结论基本一致，说明了酱香型高温大曲的微生物生态系统基本保持一致。储存过程中微生物群落与理化因子存在复杂的作用关系，大曲在储存过程中微生物群落的多样性受温度、水分和酸度等理化因子的影响，表明了储存过程中各微生物、生化指标间相互关联，性质上相互影响，使酱香型大曲在最佳的贮存期进行酿造，从而达到良好的产酒、产香风味。通过相关性计算可知微生物与环境因子之间的关系，目前主要是基于数学模型进行统计分析，但其内在的相互作用机制还有待进一步研究。理化因子在储存过程中对微生物的生长繁殖及其行为具有较强的影响，进而调控微生物群落的结构及其演替，但具体的影响过程还需进一步的实验验证。酱香型白酒高温大曲“后熟”时间的长短关系着大曲的品质及产酒的质量，但是目前没有统一的标准来界定具体的储存时间，因此，该研究为进一步阐明酱香型白酒高温大曲后熟机理提供了理论参考。
-
![]()
图 2 可培养方法分离筛选微生物形态特征
Figure 2. Microbial characteristics revealed by culture-dependent approaches
![]()
图 3 不同样品的细菌（a）、真菌（b） OTU Venn图
Figure 3. Venn profile of OTU of bacteria (a) and fungi (b) in different samples
![]()
图 4 大曲在门水平上细菌（a）、真菌（b）组成情况
Figure 4. Composition of bacteria (a) and fungi (b) in high-temperature Daqu at the phylum
![]()
图 5 大曲在属水平上细菌（a）、真菌（b）组成情况
Figure 5. Composition of bacteria (a) and fungi (b) in high-temperature Daqu at the genus
![]()
图 6 微生物群落结构与大曲理化因子相关性分析
注：a：冗余分析，GX-1~GX-6分别代表大曲储存0~150 d，理化因子箭头间的夹角代表正、负相关性；b：优势微生物群落与环境因子的相关性Heatmap图，热图中的X轴和Y轴分别为理化因子和物种（属水平），通过计算获得R值和P值来考察两个（组）变量之间的相关程度。R值在图中以不同颜色展示，热图右侧色卡是不同R值的颜色分区。*：0.01<P≤0.05，**：0.001<P≤0.01。
Figure 6. Correlation analysis between microble community structure and physical and chemical factors of Daqu
表 1 可培养微生物菌落总数变化
Table 1 Changes in the total number of culturable microorganisms
样品 细菌数
[（×104 CFU）·g−1]霉菌数
[（×104 CFU）·g−1]酵母菌数
[（×104 CFU）·g−1]GX-1 1.40±0.11f 1.40±0.26b 12.10±0.31a GX-2 5.50±0.30e 5.53±0.53a 3.14±0.19d GX-3 81.70±1.02a 0.19±0.04d 3.57±0.15c GX-4 59.20±0.58b 0.75±0.09c 4.32±0.23b GX-5 58.80±0.77c 0.53±0.03d 0.48±0.03e GX-6 47.80±0.42d 0.20±0.04d 0.21±0.04e 注：同列不同小写字母表示不同储存时间微生物数量具有显著性差异（P<0.05）。 
下载: 导出CSV
表 2 不同储存时间真菌和细菌Alpha多样性
Table 2 Alpha diversity of fungi and bacteria at different storage times
测序片段 样品 序列数目 OTUs Shannon Chao Ace Simpson Coverage 16S GX-1 22416 50 1.19 65.00 81.72 0.40 0.9993 GX-2 36379 52 1.97 56.00 57.87 0.20 0.9998 GX-3 34416 63 2.41 105.75 135.57 0.13 0.9994 GX-4 31846 94 2.42 119.20 166.07 0.17 0.9991 GX-5 36316 80 2.34 93.33 93.12 0.16 0.9996 GX-6 40855 78 2.34 117.00 147.45 0.14 0.9993 ITS GX-1 46639 21 1.18 31.00 25.64 0.36 0.9999 GX-2 53990 26 1.11 35.33 56.37 0.42 0.9999 GX-3 42845 23 0.99 30.00 44.15 0.42 0.9998 GX-4 51947 44 0.93 122.00 58.68 0.57 0.9998 GX-5 54880 30 0.87 38.25 59.97 0.48 0.9998 GX-6 45678 71 0.48 76.14 74.46 0.76 0.9998
下载: 导出CSV
-
[1] 李静心, 王艳丽, 何宏魁, 等. 基于高通量测序技术解析高温大曲和中高温大曲的真菌群落结构[J]. 食品与发酵工业,2018,44(12):52−59. [LI J X, WANG Y L, HE H K, et al. Analysis of fungal community structure of high temperature Daqu and medium high temperature Daqu based on high-throughput sequencing technology[J]. Food and Fermentation Industry,2018,44(12):52−59. doi: 10.13995/j.cnki.11-1802/ts.018043 LI J X, WANG Y L, HE H K, et al. Analysis of fungal community structure of high temperature Daqu and medium high temperature Daqu based on high-throughput sequencing technology [J]. Food and Fermentation Industry, 2018, 44 (12): 52-59. doi: 10.13995/j.cnki.11-1802/ts.018043
[2] WU Q, CHEN L, XU Y. Yeast community associated with the solid state fermentation of traditional Chinese Maotai-flavor liquor[J]. International Journal of Food Microbiology,2013,166(2):323−330. doi: 10.1016/j.ijfoodmicro.2013.07.003
[3] ZHANG X, ZHAO J, DU X. Barcoded pyrosequencing analysis of the bacterial community of Daqu for light-flavour Chinese liquor[J]. Letters in Applied Microbiology,2014,58(6):549−55. doi: 10.1111/lam.12225
[4] LIU P, ZHANG L, DU X, et al. Dynamic analysis of physicochemical and biochemical indices and microbial communities of light-flavor Daqu during storage[J]. Journal of the American Society of Brewing Chemists,2019,77(4):1−8.
[5] 刘鹏烩. 清香大曲贮存期的变化和大曲糖化酶功能菌及快速检测的相关研究[D]. 太原: 山西师范大学, 2018. LIU P H. Dynamic changes of light-flavor Daqu during storage period and researches on the glucohydrolace functional bacteria and rapid detection[D]. Taiyuan: Shanxi Normal University, 2018.
[6] 王晓丹, 雷安亮, 班世栋等. 酱香型大曲细菌的多样性[J]. 食品与发酵工业,2017,43(5):70−75. [WANG X D, LEI A L, BAN S D, et al. Bacterial diversity of Maotai flavor Daqu[J]. Food and Fermentation Industry,2017,43(5):70−75. doi: 10.13995/j.cnki.11-1802/ts.201705012 WANG X D, LEI A L, BAN S D, et al. Bacterial diversity of Maotai flavor Daqu[J]. Food and Fermentation Industry, 2017, 43 (5): 70-75. doi: 10.13995/j.cnki.11-1802/ts.201705012
[7] 张丽, 卢红梅, 陈莉, 等. 酱香型大曲贮存过程中发酵性能变化的研究[J]. 现代食品科技,2017,33(9):231−237. [ZHANG L, LU H M, CHEN L, et al. Study on changes in the biochemical properties of Maotai-flavor Daqu during storage[J]. Modern Food Science and Technology,2017,33(9):231−237. doi: 10.13982/j.mfst.1673-9078.2017.9.034 ZHANG L, LU H M, CHEN L, et al. Study on changes in the biochemical properties of Maotai-flavor Daqu during storage[J]. Modern Food Science and Technology, 2017, 33(9): 231-237. doi: 10.13982/j.mfst.1673-9078.2017.9.034
[8] 刘建文, 吴生文, 朱力红, 等. 特香型大曲中纤维素酶、半纤维素酶、蛋白酶、酯化酶活力变化规律的研究[J]. 酿酒,2013,40(5):40−42. [LIU J W, WU S W, ZHU L H, et al. Study on the change law of cellulase, hemicellulase, protease and esterase activities in special flavor Daqu[J]. Liquor Making,2013,40(5):40−42. doi: 10.3969/j.issn.1002-8110.2013.05.012 LIU J W, WU S W, ZHU L H, et al. Study on the change law of cellulase, hemicellulase, protease and esterase activities in special flavor Daqu[J]. Liquor Making, 2013, 40(5): 40-42. doi: 10.3969/j.issn.1002-8110.2013.05.012
[9] 梁丽文, 缪礼鸿, 杨团元, 等. 浓酱兼香型白酒不同储存期的高温大曲微生物群落结构与发酵特征分析[J]. 酿酒科技,2016(7):37−41. [LIANG L W, MIAO L H, YANG T Y, et al. Analysis of microbial community structure and fermentation characteristics of high temperature Daqu in different storage periods of thick sauce and flavor Baijiu[J]. Liquor-Making Science & Technology,2016(7):37−41. doi: 10.13746/j.njkj.2016084 LIANG L W, MIAO L H, YANG T Y, et al. Analysis of microbial community structure and fermentation characteristics of high temperature Daqu in different storage periods of thick sauce and flavor Baijiu [J]. Liquor-Making Science & Technology, 2016(7): 37-41 doi: 10.13746/j.njkj.2016084
[10] 梁晨, 杜海, 徐岩. 大曲贮存过程中原核微生物群落结构及风味成分演替规律[J]. 微生物学通报,2017,44(2):384−393. [LIANG C, DU H, XU Y. The succession of procaryotic microbial community and the flavor components in the storage process of Daqu[J]. Journal of Microbiology,2017,44(2):384−393. doi: 10.13344/j.microbiol.china.160163 LIANG C, DU H, XU Y. The succession of procaryotic microbial community and the flavor components in the storage process of Daqu[J]. Journal of Microbiology, 2017, 44(2): 384-393. doi: 10.13344/j.microbiol.china.160163
[11] 陈晓茹, 黄钧, 周荣清, 等. 强化曲贮存过程中群落与代谢组分的变化规律[J]. 食品与发酵工业,2022,48(3):63−69. [CHENG X R, HUANG J, ZHOU Q R, et al. Changes of community and metabolic components during the storage of fortified Daqu[J]. Food and Fermentation Industry,2022,48(3):63−69. doi: 10.13995/j.cnki.11-1802/ts.028640 CHENG X R, HUANG J, ZHOU Q R, et al. Changes of community and metabolic components during the storage of fortified Daqu[J]. Food and Fermentation Industry, 2022, 48(3): 63-69. doi: 10.13995/j.cnki.11-1802/ts.028640
[12] QB/T 4257-2011, 酿酒大曲通用分析方法[S]. 中华人民共和国轻工行业标准, 2011. QB/T 4257-2011, General analysis method of liquor making Daqu[S]. Light Industry Standards of the People's Republic of China, 2011.
[13] 赵金松, 张良, 孙啸涛, 等. 酱香型大曲微生物群落特征与理化指标的相关性研究[J]. 基因组学与应用生物学,2019,38(1):199−204. [ZHAO J S, ZHANG L, SUN X T, et al. Study on the correlation between microbial community characteristics and physical and chemical indicators of Maotai flavor Daqu[J]. Genomics and Applied Biology,2019,38(1):199−204. doi: 10.13417/j.gab.038.000199 ZHAO J S, ZHANG L, SUN X T, et al. Study on the correlation between microbial community characteristics and physical and chemical indicators of Maotai flavor Daqu [J]. Genomics and Applied Biology, 2019, 38(1): 199-204. doi: 10.13417/j.gab.038.000199
[14] 吴成, 程平言, 谢丹, 等. 酱香型白酒四轮次堆积发酵理化因子、风味物质与微生物群落相关性研究[J/OL]. 食品科学: 1−12 [2022-9-13]. http://kns.cnki.net/kcms/detail/11.2206.TS.20220622.0835.002.html WU C, CHENG P Y, XIE D, et al. Relationship analysis among physicochemical parameters, flavor compounds and microbial community during fourth round heap fermentation of Jiang-flavored Baijiu[J/OL]. Food Science: 1−12 [2022-9-13]. http://kns.cnki.net/kcms/detail/11.2206.TS.20220622.0835.002.html
[15] WANG C, WU C, QIU S. Yeast diversity investigation of Vitis davidii Föex during spontaneous fermentations using culture-dependent and high-throughput sequencing approaches[J]. Food Research International,2019,126:108582. doi: 10.1016/j.foodres.2019.108582
[16] 王颖, 邱勇, 王隆, 等. 不同产区酱香型高温大曲黑、白、黄曲的理化、挥发性成分差异性分析[J]. 中国调味品,2022,47(6):155−159. [WANG Y, QIU Y, WANG L, et al. Analysis on the difference of physicochemical indexes and volatile components of black, white and yellow koji at hight temperature with Maotai flavor in different production areas[J]. Chinese Condiments,2022,47(6):155−159. doi: 10.3969/j.issn.1000-9973.2022.06.029 WANG Y, QIU Y, WANG L, et al. Analysis on the difference of physicochemical indexes and volatile components of black, white and yellow koji at hight temperature with Maotai Flavor in different production areas [J]. Chinese Condiments, 2022, 47 (6): 155-159. doi: 10.3969/j.issn.1000-9973.2022.06.029
[17] ZHANG W X, WU Z Y, ZHANG Q S, et al. Combination of newly developed high quality Fuqu with traditional Daqu for Luzhou-flavor liquor brewing[J]. World Journal of Microbiology & Biotechnology,2009,25(10):1721−1726.
[18] 章肇敏, 吴生文, 李科发. 特香型白酒酿造过程中物质变化规律研究[J]. 食品工业科技,2012,33(9):108−110, 113. [ZHANG Z M, WU S W, LI K F. Study on the change rules of chemical composition of site-flavor fermented grains in fermentation process[J]. Food Industry Science and Technology,2012,33(9):108−110, 113. doi: 10.13386/j.issn1002-0306.2012.09.082 ZHANG Z M, WU S W, LI K F. Study on the change rules of chemical composition of site-flavor fermented grains in fermentation process[J]. Food Industry Science and Technology, 2012, 33(9): 108-110, 113. doi: 10.13386/j.issn1002-0306.2012.09.082
[19] 范文来, 徐岩. 大曲酶系研究的回顾与展望[J]. 酿酒,2000(3):35−40. [FAN W L, XU Y. Review and prospect of Daqu enzyme system research[J]. Liquor Making,2000(3):35−40. FAN W L, XU Y. Review and prospect of Daqu enzyme system research[J]. Liquor Making, 2000(3): 35-40.
[20] 李洁, 武轩. 传统凤型大曲后熟期成熟度评价及最佳贮藏期研究[J]. 食品与发酵工业,2020,46(24):97−102. [LI J, WU X. Maturity evaluation and optimum storage time of traditional Feng-type Daqu in post-maturity period[J]. Food and Fermentation Industries,2020,46(24):97−102. doi: 10.13995/j.cnki.11-1802/ts.024568 LI J, WU X. Maturity evaluation and optimum storage time of traditional Feng-type Daqu in post- maturity period[J]. Food and Fermentation Industries, 2020, 46(24): 97-102. doi: 10.13995/j.cnki.11-1802/ts.024568
[21] 左乾程, 黄永光, 郭敏, 等. 酱香型白酒机械化制曲发酵细菌群落的演替[J]. 食品科学,2021,42(18):150−156. [ZUO Q C, HUANG Y G, GUO M, et al. The succession of bacterial community in Maotai flavor Baijiu mechanized koji fermentation[J]. Food Science,2021,42(18):150−156. doi: 10.7506/spkx1002-6630-20200705-059 ZUO Q C, HUANG Y G, GUO M, et al. The succession of bacterial community in Maotai flavor Baijiu mechanized koji fermentation [J]. Food Science, 2021, 42 (18): 150-156. doi: 10.7506/spkx1002-6630-20200705-059
[22] 王晓丹, 班世栋, 周鸿翔, 等. 贵州省遵义地区3个酱香型大曲细菌群落的比较分析[J]. 食品科学,2016,37(7):110−116. [WANG X D, BAN S D, ZHOU H X, et al. Comparative analysis of bacterial communities of three Maotai flavor Daqu in Zunyi area, Guizhou Province[J]. Food Science,2016,37(7):110−116. WANG X D, BAN S D, ZHOU H X, et al. Comparative analysis of bacterial communities of three Maotai flavor Daqu in Zunyi area, Guizhou Province[J]. Food Science, 2016, 37(7): 110-116.
[23] 张立强, 冉茂芳, 王松涛, 等. 泸州老窖大曲中高产纤溶酶菌株的分离及鉴定[J]. 酿酒科技,2021(3):17−24. [ZHANG L Q, RAN M F, WANG S T, et al. Isolation and identification of high fibrinolytic enzyme producing strains from Luzhou Laojiao Daqu[J]. Liquor-Making Science & Technology,2021(3):17−24. doi: 10.13746/j.njkj.2021033 ZHANG L Q, RAN M F, WANG S T, et al. Isolation and identification of high fibrinolytic enzyme producing strains from Luzhou Laojiao Daqu[J]. Liquor-Making Science & Technology, 2021(3): 17-24. doi: 10.13746/j.njkj.2021033
[24] 李子健, 黄丹, 张曼, 等. 中高温大曲中两株耐热霉菌的分离鉴定及产酶特性比较[J]. 食品与发酵科技,2020,56(3):1−7. [LI Z J, HUANG D, ZHANG M, et al. Isolation and Identification of two thermostable fungi from medium and high temperature Daqu and comparison of their enzymatic characteristics[J]. Food and Fermentation Technology,2020,56(3):1−7. LI Z J, HUANG D, ZHANG M, et al. Isolation and Identification of two thermostable fungi from medium and high temperature Daqu and comparison of their enzymatic characteristics[J]. Food and Fermentation Technology, 2020, 56(3): 1-7.
[25] 施思, 彭智辅, 乔宗伟, 等. 浓香型大曲贮藏过程中糖化力发酵力变化及真菌多样性分析[J]. 食品与发酵工业,2017,43(5):76−79. [SHI S, PENG Z F, QIAO Z W, et al. Changes of saccharifying and fermenting power of Luzhou flavor Daqu during storage and analysis of fungal diversity[J]. Food and Fermentation Industries,2017,43(5):76−79. doi: 10.13995/j.cnki.11-1802/ts.201705013 SHI S, PENG Z F, QIAO Z W, et al. Changes of saccharifying and fermenting power of Luzhou flavor Daqu during storage and analysis of fungal diversity[J]. Food and Fermentation Industries, 2017, 43(5): 76-79. doi: 10.13995/j.cnki.11-1802/ts.201705013
[26] 董大伟, 王晓慧, 苏葛, 等. 中高温大曲存储后的感官质量研究[J]. 酿酒科技,2018(12):39−44. [DONG D W, WANG X H, SU G, et al. Research on sensory quality of medium high-temperature Daqu after storage[J]. Liquor-Making Science & Technology,2018(12):39−44. doi: 10.13746/j.njkj.2018131 DONG D W, WANG X H, SU G, et al. Research on sensory quality of medium high-temperature Daqu after storage [J]. Liquor-Making Science & Technology, 2018(12): 39-44. doi: 10.13746/j.njkj.2018131
[27] NURSTEN H E. The Maillard reaction[electronic resource]: Chemistry, biochemistry, and implications/. 2005.
[28] MARCHAND S, REVELl G D, VERCAUTEREN J, et al. Possible mechanism for involvement of cysteine in aroma production in wine[J]. Journal of Agricultural & Food Chemistry,2002,50(21):6160−6164.
[29] CHENG L, XU Y, WU Q, et al. Improving flavor metabolism of Saccharomyces cerevisiae by mixed culture with Bacillus licheniformis for Chinese Maotai-flavor liquor making[J]. J Ind Microbiol Biotechnol,2015,42(12):1601−1608. doi: 10.1007/s10295-015-1647-0
[30] 郭敏. 基于高通量测序对酱香大曲制曲微生态多样性的研究[D]. 贵阳: 贵州大学, 2018. GUO M. Study on microecological diversity of Maotai flavor Daqu based on high-throughput sequencing[D]. Guiyang: Guizhou University, 2018.
[31] 童炳丽, 刘济明, 陈敬忠, 等. 米槁根际土壤真菌多样性及其与果实药用活性成分含量的相关性分析[J]. 菌物学报,2019,38(7):1058−1070. [TONG B L, LIU J M, CHEN J Z, et al. Fungal diversity in rhizosphere soil of Migao and its correlation with the content of medicinal active components in fruit[J]. Acta mycologica Sinica,2019,38(7):1058−1070. doi: 10.13346/j.mycosystema.190062 TONG B L, LIU J M, CHEN J Z, et al. Fungal diversity in rhizosphere soil of Migao and its correlation with the content of medicinal active components in fruit [J]. Acta mycologica Sinica, 2019, 38(7): 1058-1070. doi: 10.13346/j.mycosystema.190062
[32] 朱治宇, 黄永光. 茅台镇酱香白酒不同轮次主酿区可培养霉菌种群结构多样性[J]. 食品科学,2020,41(22):184−192. [ZHU Z Y, HUANG Y G. Structure and diversity of culturable mold communities in different fermentation rounds in main Maotai-flavored liquor-brewing areas of Maotai Town[J]. Food Science,2020,41(22):184−192. doi: 10.7506/spkx1002-6630-20190918-226 ZHU Z Y, HUANG Y G. Structure and diversity of culturable mold communities in different fermentation rounds in main Maotai-flavored liquor-brewing areas of Maotai town[J]. Food Science, 2020, 41(22): 184-192. doi: 10.7506/spkx1002-6630-20190918-226
[33] 唐小飞. 西藏热泉嗜热真菌分离鉴定、产酶及生物活性初步研究[D]. 汉中: 陕西理工大学, 2019. TANG X F. Isolation, identification, enzyme production and biological activity of thermophilic fungi from hot springs in Tibet[D]. Hanzhong: Shaanxi University of Technology, 2019.
[34] 张芸曌. 中高温大曲主发酵期微生物群落与环境因子及理化性质的关联性研究[D]. 自贡: 四川轻化工大学, 2019. ZHANG Y Z. Study on the correlation between microbial community and environmental factors and physicochemical properties during the main fermentation period of medium-high temperature Daqu[D]. Zigong: Sichuan University of Science &Engineering, 2019.
[35] 杜向军, 明红梅, 马浩, 等. 机制与人工大曲曲块部位间微生物群落结构差异及关键影响因素分析[J/OL]. 食品与发酵工业: 1−10 [2022-7-22]. DOI: 10.13995/j.cnki.11-1802/ts.031775. DU X J, MING H M, MA H, et al. Analysis of microbial community structure differences and key influencing factors among different parts of mechanical Daqu and artificial Daqu[J/OL]. Food and Fermentation Industries: 1−10 [2022-7-22]. DOI: 10.13995/j.cnki.11-1802/ts.031775.
-
期刊类型引用(1)
1. 赵康妤，杨萍，马俊坤，舒文静，杨凤，谢伊莎，刘庆庆. 超声协同热处理对菜籽蛋白结构及溶解性的影响. 食品科学. 2025(01): 100-107 . 
百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 254
- HTML全文浏览量: 32
- PDF下载量: 28
- 被引次数: 1
