不同培曲工艺对特香型大曲微生物群落结构的影响

1、PAGE PAGE - 14 -不同培曲工艺对特香型大曲微生物群落结构的影响我国白酒酿造工艺历史悠久，并随培曲技术的发明与改进，推动了白酒生产的规模化、规范化发展。不同白酒大曲的原料及其配比不同，但多以小麦、大麦、豌豆等为原料压制成曲坯，入曲房自然发酵而成。大曲含有多种微生物和各种酶系，构成酿酒的内在动力；同时大曲自身所含的微生物代谢产物及原料分解产物，直接或间接构成白酒的风味物质，使白酒具有各种不同的独特风味1-3，故有“曲为酒之骨”之美誉。堆曲培养(培曲工艺)是决定大曲质量优劣的关键环节，堆曲方式与培曲条件的改变均对大曲质量造成不同程度的影响。李伟安等4采用传统与架式2种大曲培养方式，比较

2、两者大曲理化指标与微生物群落结构差异，表明架式大曲更有利于细菌、乳酸菌和霉菌的生长。陈可丹等5分析不同顶温的特香型大曲理化指标及微生物群落演替差异，表明中顶温曲的酸性蛋白酶、糖化力、酯化力均高于高顶温曲。邢钢等6跟踪测定了3种不同温度大曲制作过程中理化指标变化，结果显示大曲糖化力、酯化力、液化力、发酵力4项指标相差较大，其中中温大曲最优，低温大曲居中，高温大曲最小。明红梅等7通过提高制曲品温来探究制曲温度对酱香型大曲质量的影响，结果表明超高温大曲的感官质量优于普通高温大曲，但糖化力与酯化力明显较普通高温大曲更低。特香型白酒具有“整粒大米为原料、大曲面麸加酒糟、红褚条石垒酒窖”的工艺特征，以及“

3、浓头酱尾清中间、三香俱备犹不靠”的风味特点8。目前对特香型大曲的研究多集中在大曲理化指标及采用可培养手段对大曲进行微生物群落分析9，大曲制作过程中培曲方法对其微生物群落结构影响的报道较少见。本研究以特香型大曲为研究对象，研究不同培曲方法对大曲微生物群落结构的影响，以期解决前期发酵上霉过快、中后期品温下降迅速的问题，为进一步提高特香型大曲质量提供理论与技术参考。1材料与方法1.1试验材料1.1.1大曲样品采集在江西四特酒有限责任公司大曲车间，取4间曲房，其中2间为常规培曲工艺(traditionalculturetechniquesofDaqu，TT)，2间为新培曲工艺(newcultureTe

4、chniquesofDaqu，NT)，4间曲房分别于0(曲坯入房)、5、10、15、20、25d(曲坯出房)5点取样，粉碎后混匀、过筛，置-80冰箱贮存。1.1.2试验试剂碘、干酪素、可溶性淀粉、硫代硫酸钠、乳酸钠、柠檬酸、乳酸、三氯乙酸、羧甲基纤维素钠、L-酪氨酸、3,5-二硝基水杨、己酸，均为分析纯，西陇科学股份有限公司；-淀粉酶试剂盒，南京建成生物科技有限公司；Omega提取试剂盒Mag-BindSoilDNAKit。1.1.3仪器与设备FA2022B电子天平，上海佑科仪器仪表有限公司；PHS-3G雷磁电子pH计，上海精密科学仪器有限公司；TCL20M高速冷冻离心机，常州百密思仪器有限公

5、司；HH-3数显恒温水浴锅，国华电器有限公司；ZFD-5040电热恒温鼓风干燥箱，上海跃进医疗器械厂；T1可见光分光光度计，上海屹谱仪器有限公司；K9860全自动凯氏定氮仪、SH220N石墨消解仪，海能仪器有限公司。1.2试验方法1.2.1培曲工艺大曲原料配方、曲坯制作方法参考文献10。TT与NT的工艺区别在于：推迟曲坯覆盖稻草及翻曲时间，详见表1。表1特香型大曲不同培曲工艺盖稻草与翻曲时间单位：dTable1Straw-coveringtimesandmanualDaqu-turningindifferentTexiangDaquculturingtechniques1.2.2大曲理化指标及

6、酶活力测定大曲温度：每天上午和下午各1次，将温度计插入大曲约2cm深度，随机测定3个曲坯温度。分别计算测定温度平均值即为大曲温度。水分、酸度、酸性蛋白酶活力、纤维素酶活力、糖化力、液化力、酯化力的测定方法参考文献11。1.2.3大曲微生物量碳、氮测定大曲中微生物量碳、氮含量的测定采用熏蒸浸提法12-13。取新鲜的大曲研磨后，置于真空干燥器内用去乙醇氯仿黑暗熏蒸24h，用0.5mol/L硫酸钾溶液浸提，取浸提液，抽滤(滤膜0.5m)，获得滤液，取滤液分别用重铬酸钾氧化法与凯氏法测定碳、氮含量；另将不熏蒸大曲浸提、抽滤，测定滤液碳、氮含量，作为对照。大曲微生物量碳(microbialbiomass

7、carbon，MBC)与大曲微生物量氮(microbialbiomassnitrogen，MBN)测定如公式(1)、公式(2)所示：(1)(2)式中：EC、EN分别为熏蒸与未熏蒸大曲浸提液有机碳、氮的差值；0.38、0.45为校正系数；详细测定、计算方法参考文献14-16。1.2.4大曲总DNA提取及高通量测序样品的前处理方法参考文献17。采用Omega的提取试剂盒Mag-BindSoilDNAKit提取宏基因组DNA。用正向引物338F(5-ACTCCTACGGGAGGCAGCA-3)、反向引物806R(5-GGACTACHVGGGTWTCTAAT-3)扩增细菌核糖体16SrRNA的V3V4

8、区。用正向引物ITS5F(5-GGAAGTAAAAGTCGTAACAAGG-3)、反向引物ITS1R(5-GCTGCGTTCTTCATCGATGC-3)扩增真菌核糖体ITS1。细菌、真菌扩增产物的IlluminaMiSeq测序由派森诺科技有限公司(中国上海)完成。本文细菌、真菌群落测序数据已提交至NCBISRA数据库(http:/./traces/sra)，登录号为PRJNA716114、PRJNA716120。1.2.5数据处理本文所有实验设置3组重复，数据以平均值标准差表示。采用MicrosoftOfficeExcel2022与R语言进行数据计算和分析；单因素ANOVA分析理化指标间显著性

9、，Adonis法分析微生物群落组间显著性，P0.05表示在0.05水平上具有显著性差异；GraphPadPrism绘制折线图、R语言绘制丰度图、Canoco5绘制冗余分析(redundancyanalysis，RDA)图。2结果与分析2.1大曲发酵的理化指标及酶活性变化TT和NT培曲过程中温度变化如图1所示，TT与NT工艺大曲品温均呈现先上升后下降趋势。TT图1不同培菌工艺的曲房与大曲温度变化Fig.1DeviationofDaqucoretemperaturefromroomtemperaturewithdifferentculturetechniques注：当日日均温度数据源于天气网大曲在

10、发酵49d时，品温较NT的上升快，高出26；在第3次翻曲(10d)后，从12d开始TT品温较NT下降迅速，从12d至发酵结束时较NT的品温降低312。表明NT工艺可使曲坯发酵前期减缓升温速度、后期放缓降温速度。大曲理化指标及酶活力变化规律如图2所示。NT工艺大曲的纤维素酶活力、糖化力、液化力、酸性蛋白酶活力均显著高于TT(P0.05)，酸度低于TT，而两者的含水量与酯化力差异不显著。纤维素酶活力和糖化力从发酵5d开始NT高于TT，而液化力和酸性蛋白酶活力从发酵开始后NT一直高于TT，至25d结束分别高出63.46%、78.7%。综上结果表明NT工艺曲的酶活力优于TT。a-含水量；b-酸度；c-

11、纤维素酶活力；d-糖化力；e-液化力；f-酸性蛋白酶活力；g-酯化力图2不同培曲工艺大曲的理化指标及酶活力变化Fig.2EffectsofdifferentculturetechniquesonDaquphysicochemicalindicesandenzymaticactivities2.2大曲发酵的MBC、MBN变化大曲发酵过程中MBC、MBN动态变化如图3所示。TT与NT的MBC、MBN差异明显，但均在15d时达到峰值。NT的MBC在10d前低于TT，15d后高于TT，表明在010d期间，TT曲坯微生物生物量大于NT，并在15d时达到最大值，15d以后后者生物量大于前者。MBN是微生物

12、对氮素矿化和固持作用的综合反映18，由于MBN大小受微生物生物量、环境氮素等多重因素影响，NT的MBN仅在发酵结束时高于TT，提高了59%。2.3不同培曲工艺大曲微生物多样性分析为比较不同培曲工艺对大曲微生物多样性的影响，基于IlluminaHiSeq测序平台，对TT、NT曲坯进行多样性评价，结果见表2。每个样品文库的覆盖率(coverage)高达99.8%以上，说明此次测序结果可靠。比较各样品的指数发现，细菌各项多样性指数均表现上升趋势，而真菌多样性指数则呈现先下降后升高，TT曲坯细菌的Chao1、ACE、Simpson和Shannon指数高于同期的NT，NT真菌各项指数高于同期的TT。综上

13、结果表明，TT的细菌丰度、多样性高于NT，而真菌则相反，NT高于TT。a-大曲MBC；b-大曲MBN图3不同培曲工艺对大曲MBC、MBN的影响Fig.3EffectsofdifferentculturetechniquesonmicrobialcarbonandnitrogenbiomassinDaqu表2不同培曲工艺的大曲微生物群落多样性和丰度Table2EffectofdifferentculturetechniquesonDaqumicrobialabundanceandcommunitydiversity注：OTU为分类操作单元(operationaltaxonomicunits)2.

14、4不同培曲工艺对大曲微生物群落结构的影响选取TT、NT细菌、真菌OTU总数前10的序列，绘制物种(属水平)相对丰度柱状图，如图4-a、图4-b所示。NT和TT大曲的优势细菌属无显著性差异。魏斯氏菌属(Weissella)、乳杆菌属(Lactobacilus)、醋酸杆菌属(Acetobacter)、芽孢杆菌属(Bacillus)和片球菌属(Pediococcus)是整个发酵过程优势菌(图4-a)。其中，Weissella在第0天占据主导地位(相对丰度83.46%)，第20天时TT降低至23.8%，NT降低至13.34%；Lactobacilus在0d时相对丰度为2.2%，随后逐渐成为主要优势菌属

15、，尤其在NT第10天相对丰度达到66.9%，20d降至53.62%，TT中分别为19.28%和16.76%；Acetobacter在NT中由第10天相对丰度20.5%降到20d的4.2%，而在TT中由10d相对丰度0.047%升高到20d的4.86%。TT与NT发酵过程中优势真菌属群落间差异显著。其中嗜热子囊菌(Thermoascus)、曲霉(Aspergillus)、根毛霉(Rhizomucor)、根霉(Rhizopus)为NT、TT大曲的优势真菌(图4-b)。丝孢毕赤酵母(Hyphopichia)和假丝酵母(Candida)在第0天时占据主导地位，相对丰度分别为68.2%和15%，随后不断

16、降低，20d时在NT中的相对丰度分别为2.1%、5.5%，而在TT大曲中仅为0.3%、0.2%。值得注意的是，Thermoascus是TT大曲发酵中后期的主要优势真菌，10d时相对丰度达到95.72%，20d时降至26.7%；在NT10d时仅为0.3%，20d时提升至22.7%，是NT大曲后期的主要优势真菌。Aspergillus、Rhizomucor、Rhizopusa-细菌属水平；b-真菌属水平；c-细菌群落PCoA；d-真菌群落PCoA图4不同培曲工艺大曲微生物群落结构及PCoAFig.4MicrobialcommunitystructureandPCoAofDaquproducedwi

17、thdifferentculturetechniques在第0天相对丰度分别为0.08%、0.01%、0.1%，并随着发酵进行逐渐成为优势菌。Aspergillus在发酵第20天时相对丰度分别为TT68.1%、NT39.6%，是2种大曲发酵后期的主要优势真菌。Rhizomucor、Rhizopus则主要在NT大曲中占据主导地位，第20天时相对丰度为50.6%与13.8%(NT)远高于TT大曲第20天时相对丰度3.6%与0.8%(TT)。基于WeightedUniFrac距离主坐标分析(principalcoordinatesanalysis,PCoA)(图4-c、图4-d)，细菌群落结构的第一

18、主成分得分为64.97%，第二主成分得分为26.97%，前2个主成分积累得分91.94%(解释度大于50%)，故可基于图4-c分析样品之间细菌群落差异。0d时细菌群落结构与TT、NT的10d、20d差异显著，而TT与NT细菌群落间无明显差异，但TT10d与TT20d、NT10d与NT20d的细菌群落结构更相近。对于真菌群落结构，第一主成分得分48.63%，第二主成分得分为35.22%，前2个主成分积累得分83.85%，解释度大于50%。样品0d、TT与NT三者之间差异显著，其中TT10d与TT20d倾向聚为一簇，而NT10d与NT20d间并没有明显的成簇。2.5大曲微生物群落与环境因子、大曲理

19、化指标间的RDA将大曲样品中主要细菌属、真菌属与大曲的糖化力、液化力、酯化力、纤维素酶活力、酸性蛋白酶活力5个酶活力及酸度、含水量、品温、曲房温度4个理化指标进行RDA，如图5所示。如图5-a所示，第一与第二主轴分别解释了主要优势细菌相对丰度变化的86.11%和8.28%，总解释度为94.39%，图5-b中，第一与第二主轴分别解释了主要真核生物群落相对丰度变化的46.16%与28.9%，总解释度为75.06%，说明微生物在属水平上的分布与理化指标及酶活力之间关系密切。大曲酸度、品温、曲房温度的变化与Lactobacilus、Acetobacter、Bacillus、Thermoascus、As

20、pergillus、Rhizomucor和Rhizopus呈正相关，与Weissella、Hyphopichia和Candida呈负相关；而含水量与上述菌群的相关性和酸度、品温等的相反(图5)。酸性蛋白酶活力和液化力主要与细菌属Lactobacilus、Acetobacter呈正相关；在真菌属中Aspergillus、Rhizopus与酸性蛋白酶活力呈正相关，液化力主要与Candida、Rhizopus和Hyphopichia呈正相关。酯化力、糖化力主要与细菌属Weissella及真菌属Hyphopichia、Candida呈正相关；纤维素酶活力主要与细菌属Weissella呈负相关、与真菌属

21、Aspergillus、Rhizopus呈正相关。a-细菌群落；b-真菌群落图5大曲微生物群落与酶活力、理化指标冗余分析Fig.5RedundancyanalysesbetweenmicrobialcommunityandenzymeactivityandphysicochemicalindicesinDaqu注：优势菌群为蓝色线、酶活力为红色线、理化指标为绿色线3讨论为解决特香型大曲前期发酵上霉过快、中后期品温下降迅速的问题，本文采取NT，即延迟覆盖稻草和推迟翻曲时间。与TT相比，发酵前期NT工艺推迟1d覆盖稻草、仅翻曲1次，减缓了品温上升速度；发酵后期NT再经历第2、3次翻曲，减缓了品温下

22、降速度，高温期得到延长，而曲坯后期坚挺的品温有利于曲香生成19。微生物量碳氮常作为土壤微生物生物量及土壤肥力的检测指标20，本文首次尝试将该检测方法用于大曲的微生物生物量分析。由于NT品温前期低于TT，微生物量碳、氮均低于TT，但在15d后发生逆转，微生物生物量高于TT。即NT较TT前期微生物生长减慢，后期生长优于TT。培曲工艺调整后，微生物多样性及群落结构出现明显变化。NT的真菌丰度、多样性高于TT，而细菌则相反，TT高于NT。0d时因曲坯制作过程中添加丟糟，Weissella、Hyphopichia和Candida占据主导地位，随后NT在10d时Lactobacilus、Aspergillus、Hyphopichia和Rhizopus相对丰度高，而TT则是Weissella、Lactobacilus、Thermoascus相对丰度高，特别是Thermoascus高达95.72%，推测与品温关系密切，此时TT正处于高温发酵阶段。在20d时，NT中的真菌是Thermoascus、Rhizomucor、Rhizopus相对丰度偏高，TT是Thermoascus、Aspergillus偏