不同品种鸡蛋贮期S-卵白蛋白含量分析及其可见-近红外光谱无损检测模型研究.docx

不同品种鸡蛋贮期S-卵白蛋白含量分析及其可见/近红外光谱无损检测模型研究 付丹丹 王巧华 高升 马美湖摘要利用可见/近红外光纤光谱采集罗曼粉壳和海蓝褐壳两个品种的鸡蛋在349-1000 nm的透射光谱， 对270枚鸡蛋的天然卵白蛋白的S-型空间构象异构体（S-卵白蛋白， S-ovalbumin， S-ova）含量进行了定量分析， 实现了不同品种鸡蛋中S-卵白蛋白含量的快速无损检测通过比较贮期不同品种鸡蛋的平均光谱发现， 两个品种鸡蛋的光谱吸收峰位置相同， 仅可见光范围内的光谱吸收能量值有所不同通过标准正态变量校正（SNV）对原始光谱进行预处理， 并利用无信息变量消除算法（UVE）从500～950 nm的全光谱中提取了67个特征波长， 建立的偏最小二乘（PLS）回归模型可以很好地预测不同品种的S-卵白蛋白含量为了更进一步消除特征波长之间的多重共线性， 利用逐步回归（Stepwise regression）算法对特征波长进行二次筛选， 最终筛选出了16个特征波长， 建立多元回归模型， 其校正集的决定系数（R2）为0.9511， 均方根误差（RMSE）为0.0478， 预测集的R2为0.8380， RMSE为0.1116， 预测集相对分析误差（RPD）为2.2620。

此模型对预测集中50个罗曼粉壳鸡蛋和40个海蓝褐壳鸡蛋样本的R2分别为0.8119和0.9116， RMSE分别为0.1298和0.0834， 模型适用性更佳本研究结果表明， 可见/近红外光谱能够对不同品种的S-卵白蛋白含量进行无损检测， 建立的通用预测模型为开发便携式蛋白含量无损检测装置奠定了基础关键词可见/近红外光谱; 鸡蛋; S-卵白蛋白; 相关性; 通用模型1引 言鸡蛋作为人们日常生活中重要食物之一， 其品质的评价和预测一直是食品加工等领域的研究热点卵白蛋白（Ovalbumin， Ova）是蛋清蛋白的主要成分， 占蛋清总蛋白的54%～63%， 卵白蛋白在贮存期间不可逆地转化为一种热稳定形式S-型空间构象异构体S-卵白蛋白（S-ovalbumin， S-ova）， 贮存期间蛋清中S-卵白蛋白相对全部卵白蛋白的质量分数由新产蛋的5%上升至冷藏6个月后的81%[1]， 且其形成不依赖于鸡龄、营养状况和蛋重等因素， 显示出高的重复性与低的自然变异性， 转化速率也仅受pH值和温度的影响[2]研究发现， S-卵白蛋白与贮存时间及鸡蛋鲜度评价常用指标之间具有高度相关性[3]因此与鸡蛋品质密切相关的S-卵白蛋白被考虑作为评价商业鸡蛋鲜度的参考指标。

目前， 检测鸡蛋S-卵白蛋白含量主要有3种方法一是将蛋清稀释后水浴加热， 测定吸光度计算其含量[2]; 二是采用电化学传感器， 如分子印迹技术对蛋清中S-卵白蛋白进行特定的分子识别， 继而测定其含量[4]这两种方法均需把蛋打破， 经检测才能准确判断， 因此检测周期长， 程序复杂， 且是破坏性检测， 不能满足快速检测的要求第三种方法利用光谱无损检测技术， 如本研究组前期利采用高光谱成像技术， 利用已建立的模型对S-卵白蛋白含量进行检测[5]， 但该研究只对海蓝褐壳鸡蛋进行了研究， 且预测模型的适用性还有待进一步的考察可见/近红外光谱技术是一种成本低、速度快的无损检测技术[6～8]， 已广泛应用于鸡蛋内部品质的检测[9]和静态检测中段宇飞等[10]建立了可见/近红外光谱数据与鸡蛋新鲜度之间的非线性回归模型; Dong等[11]利用可见/近红外光谱在鸡蛋新鲜度预测的过程中比较了蛋白和整蛋的酸碱度; Syduzzaman等[12]使用600～900 nm波长范围的光谱建立了蛋黄含量预测模型; Dong等[13]指出， 可见/近红外透射光谱是检测非受精蛋的有效手段目前对鸡蛋蛋白含量的无损检测主要是通过氮元素的含量来测定整个蛋清中蛋白质的含量， 吴建虎等[14]基于可见/近红外反射光谱技术， 利用多元散射校正预处理方法和多元回归模型对新鲜蛋中的蛋白质含量进行了较好的预测。

李海峰等[15]指出， 基于400～1000 nm波长范围的光谱能够对新鲜鸡蛋的蛋白质含量进行较好的预测但是目前还没有利用可见/近红外光谱对贮期的蛋白含量， 特别是对某一种特定的主要蛋白含量进行无损检测的报道 此外， 对蛋白含量的无损检测都只是针对某一个品种， 其建立的模型适用性有限本研究采用可见/近红外光谱分析技术， 应用无信息变量消除算法（Uninformative variables elimination， UVE）及逐步回归（Stepwise regression）算法筛选出了S-卵白蛋白含量的特征波长， 并应用偏最小二乘回归分析（Partial least squares， PLS）和多元线性回归分析方法， 构建同时适用于两个品种鸡蛋S-卵白蛋白含量的无损预测模型2实验与方法2.1仪器与设备USB2000+可见/近红外光纤光谱仪（美国Ocean Optics公司）; 电子天平、JA2002数显游标卡尺（上海浦春计量仪器有限公司）; DU700紫外/可见光分光光度计（美国贝克曼库尔特有限公司）图1为测定鸡蛋的可见/近红外光谱检测系统示意图2.2实验样本和实验方法实验样品为某养鸡场的当日产的两个品种新鲜鸡蛋， 共270枚， 所有鸡蛋均无裂纹且干净。

其中海蓝褐壳鸡蛋和罗曼粉鸡蛋各135枚， 均贮藏在恒温恒湿培养箱中， 贮藏条件为22℃， 相对湿度为65%每隔5天随机取样， 即贮藏第1、6、11、16、21、26、31、36和41天时， 随机挑选出30枚鸡蛋， 其中海蓝褐壳鸡蛋和罗曼粉壳鸡蛋各15枚， 依次编号2.2.1可见/近红外光谱的采集实验采集每一枚鸡蛋的可见/近红外光谱采集参数的设定分别为：为了避免采集到的光谱数据不失真， 积分时间设置为60 ms; 为了在保证在尽量短的时间内获得更可靠的光譜数据， 平均扫描次数为3; 为了在避免较大的光谱噪音的同时不丢失光谱数据的细节特征， 平滑宽度设置为5; 光谱波段范围为原始范围， 即349～1000 nm， 光谱间隔为1 nm每次实验时， 将随机挑选出来的鸡蛋逐个放置， 以鸡蛋的长轴与地面平行的方式每次放置一枚， 以保证在采集透射光谱的时候无漏光2.2.2鸡蛋新鲜度指标的测定采集完光谱数据后， 用电子天平对每个鸡蛋称重后， 将鸡蛋打破， 放置在玻璃平板上， 然后用游标卡尺测量3个不同位置的浓蛋白（新鲜鸡蛋打破以后， 将其平摊放置， 有明显的厚度那部分蛋清）高度、蛋黄直径和蛋黄高度， 分别取其平均值作为最终每枚鸡蛋浓蛋白高度值、蛋黄直径值和蛋黄高度值。

哈夫值（Haugh unit， HU）[16]是美国农业部蛋品标准规定的检验和表示鸡蛋新鲜度的指标， 利用公式（1）求出 蛋黄指数（Yolk index， YI）是蛋黄高度与蛋黄直径的比值， 也可用于表示鸡蛋的新鲜度， 可利用公式（2）求出HU=100lg（h+7.57-1.7m0.37）（1）其中， m为蛋重（g）， h为浓蛋白高度（mm）YI=H/D（2）其中，H为蛋黄高度（mm）; D为蛋黄直径（mm）2.2.3S-卵白蛋白含量的测定将采集完光谱数据的鸡蛋打破， 参照文献[2]中的方法， 分别测定两个品种鸡蛋的S-卵白蛋白的含量 首先分离蛋清的蛋黄， 每枚鸡蛋分别取（5 0.05）g蛋清液， 各加入25 mL磷酸盐缓冲液（6.85 g NaH2PO42H2O与21.55 g Na2HPO4溶解于去离子水中， 定容至500 mL每次实验时稀释20倍， 调至pH 7.5）; 经磁力搅拌后， 每个样本各取两份5 mL混合液， 其中一份于75℃水浴加热30 min后冷却， 另一份未加热的样本作为对照分别向未加热及加热的两份混合溶液中加入沉淀剂并离心， 上清液加入双缩脲溶液混合液， 静置。

测定每个样品在540 nm处的吸光度（OD）， 以2 mL去离子水与4 mL双缩脲溶液混合液做对照按照式（3）计算S-卵白蛋白含量（S-OC）[2]：S-OC（%）=ODTOD0100% （3）其中， OD0为未加热时样品的吸光度， ODT为加热后样品的吸光度2.3样本划分与光谱预处理两个品种鸡蛋样本均为135个， 将其合并后进行建模分析， 共270个样本， 按照2∶1的原则利用基于联合X-Y距离的样本集划分（Sample set partitioning based on joint X-Y distance， SPXY）方法将所有的样本划分为校正集和预测集最终校正集样本有180个， 其中，罗曼粉壳鸡蛋85个， 海兰褐壳鸡蛋95个; 预测集样本有90个， 其中， 罗曼粉壳鸡蛋50个， 海兰褐壳鸡蛋40个各样本集的S-卵白蛋白含量描述性统计分析如表1所示由于光谱数据会受到仪器噪声及随机误差等的影响， 进而影响建模的准确性和精度， 因此需要采取合适的光谱预处理方法提高预测模型的预测能力和稳健性2.4特征波长的筛选UVE是基于PLS回归系数B的算法， 用于消除不提供信息的变量[17，18]。

依据PLS的原理， 在利用光谱矩阵X和S-卵白蛋白含量矩阵进行建模时， 两者之间存在关系Y=XB+e， 其中B是系数向量， e是误差向量， UVE是把相同与自变量数的随机变量矩阵加入到光谱矩阵中， 继而通过交叉验证建立PLS模型， 得到相应的系数矩阵、分析系数向量b的平均值和标准偏差的商t的可靠性t（i） =mean（bi）/sd（bi）（4）式中， mean（b）是系数向量b的平均值， sd（b）是系数向量b的标准偏差， i是光谱矩阵中的第i列向量依据t的绝对值大小确定是否把第i列变量用到随后的PLS模型中逐步回归法筛选特征波长是根据对回归贡献最大的解释变量相应的回归方程逐步地引入其它的解释变量[19]， 经过逐步回归筛选出来的特征波长重要且没有严重的多重共线性2.5模型的建立与评价PLS在普通的多元回归模型基础上进行了优化[20]， PLS在建模时， 先分解自变量和因变量， 同时从两类变量中提取成分， 然后按照相关性大小将提取的成分排列， 最后选择能够更好解释整个系统的新主成分参与模型的建立[21]由于在利用全波段建模时， 自变量数过多， 故选择建立PLS模型， 比较原始光谱和各种光谱预处理方法的建模效果。

对于建立的模型， 通过比较决定系数（Coefficient of determination， R2）及均方根误差（Root mean square error， RMSE）评价模型的预测性能， 其中， Rc2为校正集决定系数（Coefficient of determination of calibration set， Rc2）， RMSEC为校正集均方根误差（Root mean square error of calibraiton set， RMSEC）， Rp2为预测集决定系数（Coefficient of determination of prediction set， Rp2）， RMSEP为校正集均方根误差（Root mean square error of prediction set， RMSEP）， 同时，使用相对分析误差（Residual predictive deviation， RPD）表征模型的预测能力[22]各评价指标的相关计算公式如下：R2=ni=1（yi-）2/ni=1（yi-）2（5）RMSEC（RMSEP）=1nni=1（yi-）2（6）RPD=SDRMSEP（7）其中， yi和i分別为校正集或预测集中的第i个样本的S-卵白蛋白含量实测值或预测值， 为对应校正集或预测集中总样本的S-卵白蛋白含量平均值， SD为预测集样本的S-卵白蛋白。