中国白酒是世界上最古老的蒸馏酒之一,大曲在白酒酿造过程中影响出酒率,酒质和酒的香型[1]。大曲质量评价主要根据感官、理化和微生物指标[2-3]。优质大曲取决于在生产中对关键因素的有效控制[4]。小麦具有丰富的淀粉和蛋白质[5],可为微生物的生长、繁殖及代谢提供基础,是大曲的主要原料[6]。近年来,小麦品质对大曲和酿酒的影响日益受到重视。如花小麦所制曲块在多个质量指标上优于白小麦[7],有机原粮大曲较普通大曲更有利于酿酒产量与质量的提高[8],而添加糯小麦可提升粳高粱和普通小麦的酿酒品质[9]。
由于对原粮品质与大曲质量间关系了解较少,对专用型品种选育的目标性状尚不明确。因此,本研究选用不同品质特征小麦品种,对其品质与大曲质量指标间关系进行分析,为原粮选择和专用品种选育提供参考。
供试小麦品种为昌麦34、川麦39、蜀麦691、蜀麦482和绵麦51,均由四川农业大学小麦研究所提供(2018—2019)。
甲醛、可溶性淀粉、NaOH、己酸、无水碳酸钠、酚酞指示剂、邻苯二甲酸氢钾、无水乙醇、碘、KI、甲醛、重铬酸钾、浓H2SO4和葡萄糖,西陇化工有限公司。
JYDB 100x40小麦硬度指数测定仪,布勒粮食检测仪器无锡有限公司;UPH-Ⅲ-10T型超纯水机,四川优普超纯科技有限公司;Allegra X-30R Centrifuge 型离心机,贝克曼库尔特商贸(中国)有限公司;CP224S电子天平,德国赛多利斯公司;UV-2000型紫外可见分光光度计,上海尤尼柯有限公司;Infratec TM 1241近红外谷物籽粒分析仪,福斯分析仪器(丹麦)有限公司;FW80高速万能粉碎机,天津市工业电器厂。
1.3.1 小麦籽粒原粮品质测定方法
利用近红外谷物籽粒分析仪[10]测定小麦籽粒的品质性状,包括容重、脂肪、湿面筋含量等9个品质性状。室内测定千粒重,进行发芽试验,小麦硬度指数测定仪测定籽粒硬度(hardness index,HI),粉质率采用目测法(GB1351—2008)测定。
1.3.2 大曲发酵的方法及取样要点
将小麦除杂,润料后粉碎,加水拌料,用模具压制成型,然后进行入房培养。在0、3、6、9、15 d取样,同一时期的曲随机取3块,取3个批次的1/4粉碎混合,过40目筛于-20 ℃保存待测。具体发酵工艺参数及要点参考文献[11-12]。
1.3.3 大曲分析测定方法
水分、酸度、还原糖、氨基酸态氮、淀粉消耗、糖化力、发酵力、液化力和酯化力的测定参照行业标准《酿酒大曲通用分析方法》[13];失重变化为取样时曲块质量和0 d的差值所占的百分比;乙醇的制备为底物为8%的蔗糖溶液下,加入各发酵时期的1 g干曲粉,于30 ℃发酵72 h,其测定方法参照文献[14]进行;蛋白质质量增加量利用凯氏定氮法测定[15]。
1.3.4 数据分析
数据经过Excel初步处理后,软件Origin 6.0做图,采用SPSS 19.0进行单因素方差分析(ANOVA)、Turkey差异显著性检验、简单相关分析、回归分析和主成分分析(principal component analysis,PCA)。利用主成分因子计算主成分得分并和相应的权重进行线性加权求和[16],计算出各品种的综合评价得分。
供试品种在12个品质性状上存在显著或极显著差异(表1)。其中,昌麦34千粒重最大,蜀麦482容重最大,绵麦51在千粒重和容重上均表现为最小值。川麦39蛋白质含量最高,吸水率、湿面筋含量和沉降值也较高。昌麦34和蜀麦691的淀粉含量较高,而川麦39最低。方差分析表明(表2),蛋白质质量分数增加量和糖化力在品种间差异达到显著和极显著,水分、酸度等多个理化指标在不同发酵阶段差异显著和极显著,并且品种与发酵时间的互作对大曲品质的影响达到显著和极显著。
表1 小麦品种品质性状表现
Table 1 Qualities of different wheat varieties
材料昌麦34川麦39蜀麦691蜀麦482绵麦51籽粒颜色白白白红红发芽率/%75.66±4.04b81.33±3.05a82.33±2.52a78.33±4.04b83.33±2.08a千粒重/g56.65±0.87a44.34±0.83b41.86±0.43c44.77±0.53b40.77±0.51ac硬度65.93±0.58a69.53±0.40a58.37±0.40b51.47±0.12c44.13±0.76d粉质率/%12.67±2.08c11.33±1.53c85.67±2.08a86±1.73a74.66±3.05b容重/(g·L-1)808.4±0.60b821.0±0.63a820.6±0.78a823.1±0.81a778.7±0.39c脂肪∗/%0.63±0.02b1.22±0.01a0.4±0.03c0.95±0.02a0.65±0.03b纤维∗/%2.15±0.04c1.29±0.03d2.75±0.04b2.42±0.03b3.38±0.03a吸水率/%55.21±0.04b60.11±0.03a48.37±0.05c55.22±0.06b53.55±0.06b湿面筋∗/%31.17±0.10a33.64±0.05a20.45±0.03c27.18±0.14b26.34±0.09b沉降值/mL25.59±0.53b31.96±0.94a18.20±0.31c25.99±0.76b9.84±0.17d淀粉∗/%67.33±1.15a43.33±1.53d67.33±2.08a51.33±1.53c58.33±1.53b蛋白质∗/%13.14±0.02b14.88±0.15a10.79±0.03c11.87±0.01c8.55±0.03d
注:*,质量分数;表中不同小写字母在同行水平上表示差异显著(P<0.05)
大曲发酵过程中水分均缓慢下降(图1-a),发酵前9 d下降较快,之后趋于稳定。所有材料曲块重量降低了约25%(图1-b),蜀麦691的失重变化较平缓。酸度和氨基酸态氮与曲香呈极显著正相关[17],曲块酸度先增加后减小(图1-c)。氨基酸态氮含量呈升高趋势(图1-d),其可衡量蛋白质利用率和氨基酸转化率[3, 18-19],原料高蛋白质含量有利于微生物的生长从而提高曲块质量[20-21]。还原糖的变化差异较大(图1-e),蜀麦691持续增加后略有下降。材料间乙醇产量差异较大,受基因型影响[22],昌麦34和绵麦51先增加后下降并再次上升(图1-f),蜀麦691持续上升并稳定。淀粉消耗量持续下降并稳定(图1-g),在第15天蜀麦691消耗最多,川麦39最低。淀粉消耗能保证微生物和生化酶的富集[15, 23]。绵麦51曲块蛋白质含量增加较快,而川麦39和昌麦34变化幅度较小(图1-h)。高蛋白含量小麦通常硬度较大[21, 24],制曲中较难压制成型。
表2 不同小麦制曲发酵的品质性状方差分析(F值)
Table 2 Variance analysis of the traits among the different raws of Daqu
品质性状变异因子品种发酵时间品种×发酵时间水分0.07253.14∗∗105.29∗∗酸度0.404.83∗∗1.39∗氨基酸态氮1.787.55∗∗13.94∗∗还原糖1.412.093.03乙醇0.679.52∗∗5.17∗失重变化0.14100.17∗∗95.62∗∗蛋白增加含量15.41∗∗0.674.93∗淀粉消耗0.819.89∗∗16.56∗∗液化力1.069.94∗∗13.09∗∗糖化力3.33∗2.486.76∗∗发酵力1.694.20∗6.63∗∗酯化力0.075.14∗∗10.63∗∗
注:*表示差异显著(P<0.05);**表示差异极显著(P<0.01)(下同)
a-水分;b-失重变化;c-酸度;d-氨基酸态氮;e-还原糖;f-乙醇;g-淀粉消耗;h-蛋白质质量分数增加量
图1 大曲理化指标的变化
Fig.1 Relationship between fermentation time and the physical and chemical indicators of Daqu
供试材料在发酵中糖化力均先减小后增加并趋于稳定的变化(图2-a),蜀麦691呈波浪式变化,绵麦51变化较平缓。液化力呈现出增加的趋势(图2-b),绵麦51最高(4.20 U/g),昌麦34最低(1.64 U/g)。发酵力均呈先增加后下降的趋势(图2-c),前3 d均急剧增加,到第9天有所下降,但昌麦34表现持续增加。酯化力变化均先增加后降低(图2-d),其中绵麦51最大,蜀麦482最小,其余材料表现较为相近。
a-糖化力;b-液化力;c-发酵力;d-酯化力
图2 发酵时间与大曲质量指标的变化关系
Fig.2 Relationship between fermentation time and dynamic contents of the Daqu
简单相关分析表明(表3),小麦品质性状与大曲质量指标间存在复杂的相关关系,有12对性状间显著或极显著正相关,12对指标间存在显著或极显著负相关。籽粒容重、脂肪含量、吸水率、湿面筋含量、沉降值与大曲中蛋白增加含量极显著正相关;籽粒蛋白质含量与液化力、酯化力呈极显著负相关。以小麦品质为自变量(X)(表3),以大曲质量指标(YD)为因变量得到多元逐步回归方程。其中,糖化力回归方程为YD1=98.06+72.01X5 +4.23X9+3.53X10(R2=0.58),表明脂肪含量、沉降值和淀粉含量是影响糖化力的主要因素。液化力回归方程为YD2=32.15-0.34X4-0.07X10(R2=0.44),表明容重和淀粉含量是影响液化力的主要因素。发酵力回归方程为YD3=6.208+0.03X1-0.01X4-0.007X10(R2=0.44),表明千粒重、容重和淀粉含量是影响发酵力的主要因素。酯化力回归方程为YD4=201.38-0.21X4-0.58X10(R2=0.49),表明容重和淀粉含量是影响酯化力的主要因素。
表3 大曲质量指标与小麦品质的相关系数
Table 3 Correlation between Daqu traits and wheat qualities
品质指标X1X2X3X4X5X6X7X8X9X10X11D10.05-0.030.02-0.03-0.100.07-0.08-0.05-0.06 0.73∗∗-0.35D20.150.26-0.260.080.12-0.230.170.200.19-0.370.38D30.21-0.12-0.07-0.100.21-0.010.260.270.09-0.41∗0.37D40.360.37-0.240.29-0.06-0.26-0.0020.090.28-0.48∗0.42∗D5-0.16-0.050.003-0.21-0.110.12-0.08-0.07-0.21-0.62∗∗0.09D60.240.60∗∗-0.48∗0.62∗∗0.73∗∗-0.82∗∗0.67∗∗0.61∗∗0.85∗∗-0.42∗1D70.08-0.240.06-0.28-0.090.24-0.030.003-0.211-0.42∗D8-0.14-0.090.110.020.040.04-0.01-0.06-0.03-0.77∗∗0.34YD10.260.23-0.280.270.51∗∗-0.43∗0.51∗∗0.48∗0.50∗0.42∗0.31YD2-0.27-0.180.03-0.270.120.100.110.05-0.18-0.54∗∗0.20YD30.41∗0.08-0.31-0.170.05-0.060.230.330.08-0.40∗0.25YD4-0.150.14-0.13-0.17-0.170.02-0.12-0.08-0.18-0.63∗∗0.002
注:D1,水分;D2,酸度;D3,氨基酸态氮;D4,还原糖;D5,乙醇;D6,蛋白增加含量;D7,淀粉消耗;D8,失重变化;YD1,糖化力;YD2,液化力;YD3,发酵力;YD4,酯化力;X1,千粒重; X2,硬度;X3,粉质率;X4,容重;X5,脂肪;X6,纤维;X7,吸水率;X8,湿面筋含量;X9,沉降值;X10,淀粉含量;X11,蛋白含量(下同)
大曲质量指标间相关分析表明(表4),有24对性状间呈显著或极显著正相关,18对指标间呈显著或极显著负相关。酸度、液化力、发酵力与氨基酸态氮含量呈极显著正相关,液化力、乙醇与酯化力呈极显著正相关,淀粉消耗与糖化力、液化力和酯化力呈极显著负相关,这与前人结论一致[25]。适量的酸度可抑制有害杂菌生长繁殖,也为有益微生物提供营养和参与酯化反应生成香味物质。而水解淀粉生成的还原糖既可为微生物生长提供能量,又可生成酒精,反映了糖化与发酵速度的协调关系。还原糖和氨基酸态化合物既是微生物的碳源、氮源和能源物质,也是美拉德反应的主要反应物,可增加复合曲香。大曲糖化力和液化力影响出酒率并受淀粉消耗影响[26]。对所有指标标准化后进行主成分分析(表5),其中前3个主成分的累积贡献率为74.57%。第1主成分的贡献率为50.97%,其中主要的是水分、氨基酸态氮、乙醇、淀粉消耗、失重变化、液化力和酯化力;第2主成分指标中主要是还原糖、酸度和发酵力,贡献率为15.09%;第3主成分指标中主要是蛋白增加含量、乙醇产量、糖化力,贡献率为8.51%。
表4 大曲品质性状的相关系数
Table 4 Correlation of quality traits in Daqu
品质指标D1D2D3D4D5D6D7D8YD1YD2YD3YD4D11D2-0.281D3-0.71∗∗0.46∗1D4-0.40∗0.48∗0.381D5-0.59∗∗0.55∗∗0.360.221D6-0.350.380.37∗0.42∗0.091D70.73∗∗-0.37-0.41∗0.48∗-0.62∗∗-0.42∗1D8-0.91∗∗0.370.75∗∗0.42∗0.66∗∗0.35-0.77∗∗1YD10.47∗-0.16∗0.22∗-0.23-0.53∗∗0.310.42∗0.51∗∗1YD2-0.77∗∗0.330.61∗∗0.080.68∗∗0.20-0.54∗∗0.80∗∗-0.371YD3-0.53∗∗-0.45∗0.59∗∗0.47∗0.47∗0.25-0.40∗0.53∗∗-0.150.391YD4-0.67∗∗0.240.250.200.66∗∗0.002∗-0.63∗∗0.64∗∗-0.56∗∗0.68∗∗0.371
表5 品质评价因子成分矩阵和特征值
Table 5 Component matrix and eigenvalues of quality evaluation factors
品质指标主成分F1F2F3D1-0.89-0.260.17D20.180.78-0.09D30.720.380.17D40.090.87-0.04D50.560.34-0.54D60.380.480.62D7-0.64-0.430.25D80.880.32-0.20YD1-0.28-0.140.83YD20.900.02-0.21YD30.440.58-0.03YD40.73-0.460.03特征值6.121.811.02方差贡献率/%50.9715.098.51累积方差贡献率/%37.1558.5574.57
对各品质指标进行标准化处理,转化成均值为0、标准差为1 的无量纲数据。根据标准化后的各指标与因子载荷矩阵计算各主成分得分,公式如下:
PCA1=-0.37D1+0.23D2+0.30D3+0.21D4+0.32D5+0.15D6-0.33D7+0.38D8-0.22YD1+0.33YD2+0.26Y D3+0.29YD4;
PCA2=0.04D1+0.29D2+0.22D3+0.36D4-0.20D5+0.55D6-0.002D7-0.04D8+0.45YD1-0.18YD2+0.21YD3-0.34YD4;
PCA3=0.25D1+0.43D2-0.26D3+0.50D4+0.19D5-0.25D6-0.04D7-0.19D8-0.38YD1-0.38YD2+0.13YD3+0.03YD4;
以3个主成分对应的方差相对贡献率为权重,计算各品种的综合评价得分,公式为 F综=0.68PCA1+0.20 PCA2+0.11 PCA3。
根据综合得分模型计算出各主成分的得分,第1主成分得分川麦39和绵麦51较高,第2主成分得分蜀麦482和昌麦34较高,第3主成分得分蜀麦691和昌麦34较高。而综合得分最好的是昌麦34和蜀麦691,表明其制曲的综合品质相对较好。但昌麦34的粉质率较低,硬度较高,不利于曲块的压制。
本研究利用5个品质特征差异较大的小麦品种制曲,对小麦的11个品质性状和大曲发酵中5个阶段的12个质量指标进行了测定,并采用方差分析、多元回归分析和主成分分析对性状间关系进行了深入研究。结果表明使用不同品质特征的小麦品种,对大曲质量指标有显著的影响,因此生产中有必要选育和使用制曲专用型小麦品种。大曲主要质量指标受小麦淀粉含量、千粒重、容重、沉降值和脂肪含量等的影响,因而在品种选择时要关注多个品质性状的相互组合。尽管多元回归中未检测到小麦蛋白质含量与大曲主要质量指标间的显著关系,这可能是因为材料间蛋白含量区间较小,材料数有限。通过主成分分析,利用性状各主成分值和权重所得主成分综合得分的方法来评价小麦品种制曲优劣,能更好地筛选小麦制曲原粮,克服了传统方法根据少数指标选择的缺点,可应用于生产实践。
[1] 张立新.小议制酒大曲原料[J].中国酿造,2005(3):36-39.
[2] 邢 钢,敖宗华,王松涛,等.不同温度大曲制曲过程理化指标变化分析研究[J].酿酒科技,2014(6):20-24.
[3] XIAO C,LU Z M,ZHANG X J,et al. Bio-heat is a key environmental driver shaping the microbial community of medium-temperature Daqu[J].Applied and Environmental Microbiology,2017,83(23):1-13.
[4] 沈毅,陈波,张亚东,等.影响高温大曲质量的关键控制点[J].酿酒科技,2019(8):1-6.
[5] 石冬梅,李渐鹏,邱斌,等.小麦精深加工产业发展概况[J].中国食物与营养,2016,22(15): 41-43.
[6] 安明哲,李杨华,宋延富,等.酿酒小麦流变学特性的研究[J].酿酒科技,2019(4):137-141.
[7] 余有贵,罗俊,熊翔,等.小麦品种对偏高温大曲发酵动态影响的研究[J].食品科学,2007,28(10):405-409.
[8] 涂荣坤,蔡小波,沈才萍,等.有机大曲与普通大曲的差异性研究[J].酿酒科技,2014(1):1-5.
[9] 赵国君,李斌,徐智斌,等.添加糯小麦对普通小麦和粳高粱小曲酒酿造特性的影响[J].麦类作物学报,2013,33(5):942-945.
[10] 姬玉梅,孙华.小麦主要品质指标的近红外光谱分析模型建立[J].食品研究与开发,2017,38(21): 142-145.
[11] 马加军,吕月明.包包曲生产工艺条件控制[J].酿酒科技,2011(11):92-95.
[12] 王贤,丁琳.中高温大曲制曲过程关键控制点探析[J].酿酒科技,2015(3):87-90.
[13] QB/T 4257—2011 酿酒大曲通用分析方法[S].北京:中国标准出版社,2011.
[14] 岳元媛,向文良,张文学,等.白酒发酵糟醅中残余酒精含量检测方法的探讨[J].中国酿造,2005(7):55-58.
[15] 任道群,唐玉明,姚万春,等.泸州老窖大曲化学特性差异研究[J].酿酒科技,2005(5):51-55.
[16] 吴松霞,郜海燕,韩延超,等.基于主成分分析的不同品种茭白品质评价[J].中国食品学报,2019,19(7):1-11.
[17] 曹振华,沈才洪,秦辉,等.泸型大曲中曲香与微生物和生理生化指标间的关系[J].酿酒科技,2016(3):42-45.
[18] 刘延波,刘润雨,赵志军,等.赊店老酒大曲质量差异研究[J].食品工业科技,2019,40(23):77-80;87.
[19] WANG Z M,LU Z M,YU Y J,et al. Batch-to-batch uniformity of bacterial community succession and flavor formation in the fermentation of Zhenjiang aromatic vinegar[J].Food Microbiology,2015,50:64-69.
[20] 赵国君,徐智斌,冯波,等.糯小麦的酿酒特性研究[J].中国农业科学,2013,46(6):1 127-1 135.
[21] 冯雨.小麦软硬度对高温大曲的影响[J].酿酒,2019,2(46):88-90.
[22] 丁丽,张志鹏,彭远松,等.小麦籽粒品质性状与葡萄糖和乙醇产量的关系分析[J].食品与发酵工业,2019,45(23):73-78.
[23] 李斌,徐智斌,冯波,等.糯小麦与普通小麦的固态发酵特性比较[J].麦类作物学报,2012,32(5):949-954.
[24] 石玉,张永丽,于振文.小麦籽粒蛋白质组分含量及其与加工品质的关系[J].作物学报,2009,35(7):1 306-1 312.
[25] 炊伟强,敖宗华,张春林,等.泸州老窖大曲感官特征与微生物、理化指标和生化性能的关联研究[J].食品与生物技术学报,2011,30(5):761-766.
[26] 陈泽军,周瑞平,尹礼国,等.酿酒专用小麦在多粮浓香型酒厂的应用[J].酿酒科技,2008(5):65-68.