麸皮是谷物加工的主要副产物,我国每年产出的小麦麸皮在2 000万t以上
[1]。麸皮是麦粒的外层,约占整个麦粒的15%,小麦麸皮由13.18%蛋白质、14.25%淀粉、3.4%脂肪、3.8%矿物质和55.60%非淀粉碳水化合物组成
[2-3]。木聚糖是小麦麸皮中碳水化合物的主要成分,约占小麦麸皮的20%,是制备小麦麸皮还原性糖的良好资源。目前麸皮大部分用作动物饲料,不仅经济利用效益低,而且食用和商用价值均较低
[4-5]。近年来,研究发现麦麸不仅具有对人体必不可少的生理作用,包括降血糖、抗氧化、抗炎
[6]、抗肿瘤和免疫调节作用
[7-8],同时也可以提高畜禽的生长性能、改善动物产品品质和提高生产效益
[9]。还原糖是指能够被还原的糖,包括葡萄糖和果糖。发酵麸皮中还原糖含量不仅反映了乙醇的生产和发酵状态,而且也反映了发酵麸皮中的微生物代谢活性
[10],从而间接的反应发酵过程。麦麸蛋白质含量丰富,且氨基酸构成比例均衡,包含人体所有必需氨基酸,是一种优质蛋白质资源
[11]。小麦麸皮蛋白可以制备成抗氧化活性较高的小分子小麦麸皮多肽,具有多种生物功效
[12-13]。相关研究证实,畜禽对饲粮的偏爱程度与其营养物质种类和含量相关,故对麦麸中营养物质的快速检测尤为重要
[14-15]。然而,在传统的饲料生产过程中,发酵料的测定依赖于工人或常规的物理和化学方法,还原糖在试剂中溶解度较小,所测结果偏低;蛋白质的测定需要浓硫酸等危险试剂,这些检测方法存在着耗时、及操作复杂等问题,难以满足通过快速简便准确测定发酵过程中的活性成分含量的需求。
近年来,光谱学技术因其简单、快速、无损和低成本的优势
[16],已被用于以快速和无损的方式评价食品的质量,如核桃
[17]、苹果
[18]、桃果
[19]。特别是近红外(Near infraredred, NIR)光谱定量分析技术的研究最为广泛。NIR通过主成分分析、人工神经网和偏最小二乘法(Partial least-square method, PLS)等化学计量学方法,建立物质光谱与待测成分含量的线性或非线性预测模型,从而实现对物质中待测成分的快速定量检测
[20]。据报道,近红外光谱可用于水稻水分和蛋白质含量的高精度测量
[20-24]。祁梅等
[25]利用NIR技术建立麦冬5种代表性成分的PLS定量分析模型,结果表明预测结果准确,可用于麦冬5种成分含量的快速测定。Xie等
[26]建立基于全波段和PLS的近红外模型预测油炸和烘烤马铃薯片中丙烯酰胺含量。魏尊苗等
[27]利用近红外光谱技术建立了关于粗脂肪、蛋白质、淀粉和总糖的油莎豆的品质检测。但目前还没有针对发酵麦麸的营养指标建立快速预测方法相关研究。
因此,本研究首先借助近红外光谱技术,建立发酵麦麸中还原力和可溶性性蛋白质含量的快速检测方法,并进一步通过光谱预处理及特征波段筛选等方法优化模型,增强检测准确性。本研究为实现高效便捷的发酵麦麸特征营养成分含量的检测提供参考和技术支持。
1 材料与方法
1.1 材料及仪器
麸皮、玉米和豆粕采购自当地市场(呼和浩特);商业复合菌(枯草芽孢杆菌1∶地衣芽孢杆菌1∶植物乳杆菌1∶产阮假丝酵母3)购买于市场、纤维素酶(酶活10 000 U)购自下盛夏(北京)生物技术发展有限公司;DNS试剂(上海科丰化学试剂有限公司),均为分析纯。
便携式近红外光谱仪(MicroNIR Onsite) VIAVI Solutions公司;智力光照培养箱 宁波东南仪器有限公司;高速粉碎机 浙江屹立工贸有限公司;电热鼓风干燥箱 上海一恒科学仪器有限公司;Epoch2 酶标仪(微孔板分光光度计)美国伯腾仪器有限公司。
1.2 方法
1.2.1 样品的制备与处理
发酵麸皮的制备参考杜涓等
[28];以小麦麸皮(80.46%)为原料,以豆粕粉(9.32%)、玉米粉(10.22%)为辅料,料水比1∶1.16;加入1 000 U/g 0.1%的复合菌,2 g纤维素酶,于28 ℃下发酵4 d,每12 h取样,为保证样品的均匀性,把收集到的麸皮样本烘干后用高速粉碎机粉碎后过筛,共得到发酵麸皮样本310份。
1.2.2 还原糖及可溶性蛋白的测定
发酵产物还原糖含量采用DNS显色法测定
[28],建立标准曲线方程
y=0.838 4
x+0.026 6。以1∶20的料水比称取发酵麸皮粉末,80 ℃热水提取40 min,冷却后离心,取上清稀释相应倍数。取稀释样液1.0 mL,DNS试剂2 mL,摇匀,沸水浴加热5 min,立即取出,用流水迅速冷却试管。然后加12 mL蒸馏水定容并摇匀后,在540 nm波长处测定吸光度。
采用BCA蛋白浓度测定试剂盒测定发酵产物的可溶性蛋白含量。
1.2.3 近红外光谱测定条件
本试验使用VIAVI Solutions公司便携式近红外光谱仪(MicroNIR Onsite)对发酵麸皮样品进行近红外光谱扫描。为保证仪器的稳定性,开机后预热30 min;在采集过程中,使样品粉末自然填满样品容器,保证样品表面平整。仪器先以白板作为参比,在采集样品光谱过程中,保持垂直照射并避光检测,以漫反射方式设置光谱扫描波长范围为908~1 676 nm,测定单个样品的近红外光谱数据3~5次,积分时间为50 ms,对光谱值进行均值化处理,用于建模分析。通过内部交叉验证的方法,对异常样本剔除并采取逐一回收的方式,剔除时观察模型评价指标,判断是否剔除,以保证最大程度的减小异常样本对模型预测能力的影响,提高构建发酵麸皮还原糖和蛋白含量预测模型的精准度。
1.2.4 校正集和验证集样品的选择
样品按4∶1的比例随机分为校正集和验证集,校正集用于建立测定样品中还原糖和可溶性蛋白定量校正模型,验证集用于验证所建立校正模型的预测能力。
1.2.5 近红外光谱数据建模预处理方法
近红外光谱易受到样品的杂散光以及噪音信号等影响,为保证谱图基线平稳和所建立校正模型的准确性,需对原始光谱进行预处理
[29]。本研究采用预处理方法有无光谱处理(None)、去趋势(Detrend)、标准正规变换(SNV)、二阶导数(SD)、一阶导数(FD)及4种方法任意组合。
1.2.6 特征光谱波段的选择
偏最小二乘法(Partial least squares,PLS)可以处理全光谱信息,但这些信息中包含了过多的冗余信息,通过特征波长的筛选可有效地降低建模的复杂性,提高了模型的预测精度
[30]。因此,在建立模型时,需要找到能充分反映出样品中各指标成分的最佳建模波段来改善模型性能,提高计算速度。当RMSEC、SEC最小,
R2、RPD值最大时,对应的波段为最优特征波段。
1.2.7 主成分因子数的选择
通过PLS建立NIR分析模型时,主因子数也对模型预测准确度有较大影响。主成分太小,则导致建模信息不全,模型预测能力相对较差;数值太大,所建模型包含太多的噪音,会出现过拟合现象,当验证集预测标准偏差最小时,对应的主成分最佳
[31-32]。
1.2.8 模型的质量评价
首先用内部交互验证法检验模型的稳健性,通过校正集的决定系数(Correction determination coefficients,
Rc
2)和校正标准差(Standard error of the calibration,SEC)来判断模型的质量,一般校正决定系数
Rc
2越大,近红外分析结果与化学分析结果越吻合,可信度越高。用验证集的决定系数(Determination coefficient,
Rp
2)和预测均方根误差(Root mean square error of prediction,SEP)来判断模型的精确性,当
Rp
2越大,SEP越小,说明定标模型的准确度越高,则模型预测性能越好。最后根据相对分析误差(Relative percent deviation,RPD)对预测精度进行进一步的评价。当RPD值<2.5时,说明比较困难用该模型进行定量分析,故无法用近红外光谱技术进行分析;若2.5≤RPD<3,说明模型对该成分进行定量分析是可行,但其精度有待进一步提高;若RPD值≥3,说明建立的定标模型预测效果良好,可用于实际检测
[33-34]。
1.3 数据处理
使用Origin2021作图;采用Excel 2019进行数据汇总处理;以SAS 9.2 各样本中发酵麸皮还原糖和蛋白质进行统计分析,使用近红外光谱仪自带软件 MicroNIR TM Pro v3.1对近红外光谱图进行处理和模型建立,用IBM SPSS Statistics 25分析显著性差异。
2 结果与讨论
2.1 麸皮样品化学值
由
图1可知,随发酵时间的延长,麸皮中的还原糖呈先上升后下降再上升的趋势,0~48 h还原糖含量由12.82 mg/g上升至40.61 mg/g(
P<0.001),在60~72 h显著下降随后缓慢上升;可溶性蛋白呈先上升后下降的趋势,其中发酵时间为0~48 h时可溶性蛋白由11.18 mg/g上升至40.45 mg/(
P<0.001);发酵到60 h之后,可溶性蛋白的含量均随着发酵进行显著下降(
P<0.05)还原糖是发酵体系中微生物可以直接利用的碳源,其在发酵过程中的含量变化趋势直接反映微生物的生长情况。本试验中说明发酵12 h以上时,发酵麦麸中的微生物大量繁殖,并分泌糖化酶,酶分解淀粉使得发酵体系中的还原糖含量急剧上升。发酵48 h之后,一方面糖化酶活力下降,另一方面还原糖被微生物用于乙醇发酵。在84 h蛋白被分解利用还原糖含量升高。麸皮在发酵过程中,蛋白质被蛋白酶分解成多肽和一些小分子物质,增加了蛋白质的溶解性,所以可溶性蛋白质的含量增加;而在发酵后期可溶性蛋白更多地被降解成氨基酸和小肽物质,从而造成后期可溶性蛋白含量呈下降趋势。
建立模型的关键因素是建模所用的定标集化学值是否覆盖广泛,一个覆盖范围广的定标集能够提高检测模型的精度和稳定性
[35]。由
表1可知,试验所选取的发酵麸皮样品中还原糖含量范围为12.82~40.61 mg/g,可溶性蛋白含量范围为11.18~40.45 mg/g,覆盖广泛,能够较好的满足模型构建的需求,基于此建立的近红外光谱模型更具有普适性。因此,选取的发酵麦麸样品可用于建立其还原糖和可溶性蛋白的近红外光谱预测模型。
2.2 麸皮样品的近红外原始光谱
239个发酵麸皮还原糖样品的原始光谱反射率曲线如
图2(a)所示,179个可溶性蛋白样品原始光谱反射率曲线如
图2(b)所示,由图中可以看出近红外光谱图走向基本一致,但响应信号存在一定差异,说明样品中的各成分含量存在差异。在1 200 nm和1 450 nm附近拥有较为明显的吸收峰,这与发酵麸皮中大量含氢基团振动有关。
2.3 光谱预处理方法的选择
由
图2可以看出,原始光谱存在明显的基线漂移情况,为确保模型的预测性能,需要对其进行预处理。光谱数据预处理方法包括无处理(None)、去趋势化(Detrend)、标准正态变换(SNV)、二阶导数(SD)、一阶导数(FD)5种单一预处理方法,或上述2~3种预处理方法随机组合的方法。其中,SEC越小,
Rc
2越接近1,RPD越大,表明所建立模型的预测效果更好、适用性更强(
图3)。
表2和
表3为最终筛选最佳预处理方法建立定量模型。其中还原糖和可溶性蛋白模型选择一阶导数(FD),二阶导数(SD)和标准正态(SNN)变换获得的PLS模型最优。其模型的决定系数
Rc2分别为0.854 8和0.833 2,验证均方差SEC分别为1.947和3.29,且2种模型的RPD均>2,说明模型具备较高的可靠性,其中还原糖模型RPD为2.624>2.500能够用于还原糖含量的实际检测,而可溶性蛋白模型RPD为2.448 6<2.500,说明该模型预测结果可达到粗估的效果,但预测精度需进一步提高。
2.4 最适波段的选择
变量数过多会增加预测模型复杂度,过少则易遗漏细节信息,因此采用特征波长提取能够提高其准确度与精度
[36]。建模区间的筛选是将还原糖光谱区间划分为6个子区间(908.1~1 676.2 nm、908.1~945.27 nm、945.27~1 075.35 nm、1 075.35~1 341.7 nm、1 341.7~1 577.09 nm、1 577.09~1 676.2 nm);建模区间的筛选是将可溶性蛋白光谱区间划分为5个子区间(908.1~1 676.2 nm、908.1~1 261.18 nm、1 261.18~1 397.45 nm、1 397.45~1 546.12 nm、1 546.12~1 676.2 nm),通过比较不同波段所建立模型的校正样品预测标准偏差(Standard error of calibration, SEC)、决定系数(Determination coefficient,
Rc
2)和相对分析误差(Relative percent deviation,RPD)筛选最优波段参与建模。由
表4可知,还原糖及可溶性蛋白最优波段为908.1~1 676.2,且RPD>2.5,模型具备较高可靠性,能够用于模型分析。
2.5 还原糖和蛋白SEP随主成分数的变化图
在建模过程中,模型的主因子数是模型的关键参数之一。主因子数过低,模型的拟合度不足、性能不佳;主因子数过高,模型则容易出现过拟合。因此,合适的主因子数对定量分析模型性能至关重要。一般情况下,以SEP为指标,当验证集预测标准偏差最小时,对应的主成分数最佳
[37]。不同主因子下,模型的SEP值如
图4所示。由
图4可以看出,可溶性蛋白模型主因子数为10时模型的SEP值降至最低,随后随着主因子数增加模型的SEP值平缓上升;而还原糖模型主因子数为7时模型的SEP值降至最低,随后随着主因子数增加模型的SEP值逐步上升。因此,还原糖和可溶性蛋白模型的最佳主因子数为7和10。
2.6 定标模型的建立
以筛选的最优光谱预处理方法,最适波段和最佳主因子数作为参数评价指标,对模型进行智能优化,然后采用PLS法建立定量分析模型。各模型实测值与预测值的相关性见
图5,模型的建模结果如
表5所示。发酵麸皮的还原糖和可溶性蛋白定量模型的决定系数
Rc
2分别为0.904 8和0.938 2。
2.7 定标模型的外部验证及评价
还原糖和可溶性蛋白模型外部验证集预测值和实测值相关关系图如
图6所示。2种模型决定系数
R2分别为0.954 5和0.994 4均>0.80,验证集样品的预测值与实测值均围绕在
Y=
x的参考线附近,二者十分接近模型性能理想。
将验证集(
n=20)样品NIR光谱输入定量模型,得出预测值,同时借助SPSS统计学软件将对应预测值与实测值的比值与数值1作显著性分析,结果如
表7和
表8所示,二者不存在显著差异(
P分别为0.959 5和0.901 9)。说明所建立的发酵麸皮2个活性成分定量模型稳定性和准确性较好,且模型的预测准确度良好。
3 讨论
小麦麸皮为小麦加工后的副产品,具有悠久的历史。麦麸中有效成分丰富,包括膳食纤维、麦麸蛋白、酚类(阿魏酸、黄酮类和木酚素等)、植酸和γ-氨基丁酸等,具有抗氧化、降糖、降脂、抗癌、促进肠道蠕动等作用
[38]。但成分测定方法前处理过程繁琐,且溶剂会造成环境污染。因此,开发快速、可靠、无损检测麸皮含量的方法尤为重要。本研究中所采集的麸皮样品来自不同发酵时间段,时间跨度较长,并且样品的活性成分含量变化幅度大,覆盖面广泛,有利NIRS与化学分析方法联用,一定程度上能够提高模型预测的准确性和可靠程度。为了进一步优化模型,人们采用NIRS与多元线性回归法(MLR)、主成分回归法(PCR)、改进偏最小二乘回归法(MPLS)和偏最小二乘回归法(PLS)相结合使预测的结果更为可靠
[39-40]。本研究应用NIRS与偏最小二乘法(PLSR)建立麸皮还原糖和可溶性蛋白含量的近红外模型,通过选择不同波段,不同主成分因子,采用不同预处理方式及组合对模型进行优化,发现全谱图FD+SD+SNV处理所得的模型结果最优,通过交叉验证评价模型预测效果,最终得到准确度较高的还原糖和可溶性蛋白含量的预测模型。用优化后的模型的预测值对20个未参与建模的检验样本的化学值进行配对
t检验,
P=0.959 5、
P=0.901 9>0.05(置信区间为95%),说明2种方法得到的检测结果无显著差异,这与陆登俊
[41]的研究结果一致。而田丽梅等
[42]利用NIRS技术对燕麦粗蛋白质进行预测发现,在950~1 650 nm全波段对光谱进行SNV+MSC+Detrend+FD预处理,模型的
R2和RPD分别可以达到0.88和3.00,说明对光谱进行预处理使十分必要的。苏鹏飞等
[43-44]利用近红外光谱技术结合常规分析方法,建立了高粱水分和淀粉含量的快检模型,可用于实际生产分析检测,提高样品分析效率。这也从侧面佐证了通过近红外光谱技术可以快速准确的检测样品成分含量。利用近红外方法分析预测依赖于传统化学分析方法,化学值测定的准确性将直接决定模型检验的可靠性与适用性。因此,通过改进传统的还原糖和可溶性蛋白含量测试方法,提高原始数据准确度的同时,在一定程度上也能够提高近红外技术预测的准确度,增强模型预测能力。本研究中,用于建模的样品还原糖和可溶性蛋白含量范围为12.82~40.61 mg/g和11.18~40.45 mg/g,可覆盖大部分麸皮还原糖和可溶性蛋白含量范围。通过适当增加样品数量,或者对样品还原糖和可溶性蛋白含量进行分段建模,可建立适用范围更广、精确度更高的模型。本研究对比了不同波段及全波段的影响发现,在908~1 670 nm时,模型的预测效果最优,但是否可以使用部分相关性显著的波长段或者波长段组合通过不同的数据处理方式建立相关性更高的近红外光谱模型,还有待进一步研究。除此之外,主因子数是模型的关键参数之一,唐保山等
[45]利用近红外光谱技术对五倍子主成分的含量进行快速检测时,以相关系数(
Rc
2、
Rp
2)和误差(SEC、SEP和RPD)评估模型预测性能,确定了五倍子中水分、单宁酸和没食子酸含量的近红外模型最佳因子数为9、10和5。
4 结论
本研究以FD+SD+SNV光谱预处理,基于PLSR优选出全波段光谱区域主成分变量数为7和10,成功构建近红外法测定麸皮还原糖和可溶性蛋白含量。还原糖PLSR定量模型Rc2、SEC、SEP、Rp2和RPD分别为0.904 8、1.576 1、2.408 9、 0.954 9和2.791 4;可溶性蛋白PLSR定量模型Rc2、SEC、SEP、Rp2和RPD分别为0.938 2、2.003、1.547 3、0.994 4和2.970 4,模型预测准确性高。本研究表明近红外光谱法结合合适的化学计量学技术可快速准确预测麸皮还原糖和可溶性蛋白含量,此外,该模型可以实时检测发酵过程中麦麸有效成分的动态变化,有望为麸皮质量评价提供新的技术支撑。
内蒙古自治区科技重大专项(2021ZD0023-3)
内蒙古自治区科技重大专项(2021ZD0024-4)
国家乳业技术创新中心创新能力建设重点项目(2022-科研攻关-2)
内蒙古自治区科技计划项目(2022YFSJ0029)
内蒙古自治区科技计划项目(2022YFHH0072)