基于油不饱和度的便携式Raman光谱快速无损鉴别山核桃品种

《Food Control》：Rapid and non-destructive discrimination of pecan cultivars by oil unsaturation using portable Raman spectroscopy

　　

利用便携式785 nm拉曼光谱通过油不饱和度快速无损鉴别山核桃品种。便携式785 nm拉曼光谱可非破坏性地按油不饱和度对山核桃品种进行排序；1656、3010和1265 cm^?1处的三个一致谱带标记对排序结果一致；过氧化值(peroxide value, PV)可追踪初级氧化，而拉曼指纹保持稳定；相同的品种梯度在提取油和完整籽仁中均出现；该方法无需试剂，可实现快速的在线筛选。

美国山核桃（Carya illinoinensis）是重要的经济坚果，其油占籽仁质量的65%~75%，营养价值主要来源于油。油富含不饱和脂肪酸，油酸(oleic acid, C18:1)和亚油酸(linoleic acid, C18:2)构成脂肪酸的绝大多数，两者的比例（油酸/亚油酸比）因品种和产地而异，已成为生产与育种中的核心品质指标。亚油酸含两个双键，比单不饱和油酸更易氧化，因此高油酸品种更耐储存，高亚油酸品种易在储存中产生酸败。不饱和度直接决定山核桃的货架期、风味保持能力与市场价值。然而，既有的参考方法（如脂肪酸甲酯气相色谱法、碘值滴定等）需要油提取、化学试剂、专业人员及较长的分析时间，难以满足育种圃、收购现场或加工线上对大量样品进行快速筛选的需求。拉曼光谱可探测分子振动，无需制样、对水不敏感，且已有便携式785 nm仪器，为山核桃品质检测提供了潜在方案。已有研究对山核桃油进行拉曼表征并与其脂肪酸组成关联，但尚未确立基于油不饱和度的品种判别，也未验证能否在完整籽仁上实现非破坏性测量。

研究人员使用便携式785 nm拉曼系统，对四个品种的提取油（n=72）与四个品种的完整籽仁（n=131）进行了系统测量。结果表明，在提取油中四品种形成了可重复的不饱和度梯度（由低到高为Desirable、Schley、Stuart、Moneymaker），三个独立的不饱和度敏感谱带标记（1656、3010、1265 cm^?1）给出的排序完全一致；在完整籽仁上同样恢复出显著梯度（Creek、Cape Fear、Desirable、Sumner），说明无需提取油即可判别。过氧化值测定确认油在储存期间发生了显著初级氧化（PV从3.1升至17.5毫当量/千克），但拉曼指纹和品种排序保持稳定，表明该指纹反映的是内在组成差异而非氧化状态。该研究证明便携式Raman光谱可作为快速、无损、无需试剂的检测工具，为山核桃品质筛选和育种评价提供了实用新途径。

研究采用MarqMetrix All-In-One便携式拉曼分析仪（785 nm激发，450 mW功率，2000 ms积分，20次累加）采集光谱；油样各得一条光谱，籽仁在子叶表面三点平均以抑制表面粗糙度影响。预处理包括非对称最小二乘基线校正、归一化至2800~3025 cm^?1的C-H伸缩区面积、标准正态变量(SNV)缩放与Savitzky-Golay二阶导数。三个不饱和度指数按固定峰面积比计算：1656/1440、3010/CH伸缩区、1265/1440。统计方法包括单因素方差分析(ANOVA)、Welch检验、主成分分析(PCA)、ANOVA-同时成分分析(ASCA)与偏最小二乘判别分析(PLSDA)。油样为四个品种经三种温湿度储存条件、79天内六次冷榨，共72样；籽仁为四个品种共131个光谱，样品来自美国佐治亚和密西西比州多个果园的商业混合批次。

3.1 山核桃油拉曼光谱特征
通过采集四个品种提取油的平均光谱，所有诊断脂质谱带（如1440 cm^?1的CH₂剪式振动、1656 cm^?1的顺式C=C伸缩、3010 cm^?1的=C–H伸缩等）均被清晰分辨，证实了波数归属的正确性和采集质量。

3.2 提取油中的品种不饱和度梯度
在指纹区，四品种平均光谱高度相似，差异仅集中于不饱和度敏感谱带。以1656/1440指数衡量时，品种形成显著且单调的不饱和度梯度（Desirable、Schley、Stuart、Moneymaker），三个标记均极显著（F=75.0、68.2、70.9，p<10^?20），所有两两对比均显著，且三个标记排序完全一致。

3.3 完整籽仁的非破坏性确认
对131个籽仁光谱进行相同分析，同样恢复出显著的不饱和度梯度（Creek、Cape Fear、Desirable、Sumner），三个标记均显著（F约21~27，p<10^?10）且排序一致，表明可在不提取油的情况下直接区分品种。

3.4 提取油中的多变量判别
主成分分析(PCA)中第一主成分承载80%方差，载荷峰值位于1656、1265和1440 cm^?1；偏最小二乘判别分析(PLSDA)的分类准确率为57%（重复交叉验证）和62.5%（留一天法，置换p=0.003），混淆矩阵显示Desirable和Moneymaker分别以83%和72%被正确识别，相邻的Schley和Stuart互相混淆较多。

3.5 ASCA方差划分
ANOVA-同时成分分析(ASCA)将61.9%的光谱方差归于品种因子（p<0.001），储存条件因子仅占1.5%、储存天数因子占6.3%，均不显著，残差为30.3%。品种效应载荷集中于1656和1265 cm^?1谱带，直接支持不饱和度作为判别依据。

3.6 品种指纹对储存氧化的稳健性
过氧化值随储存时间单调上升（Spearman r=0.99），但Raman不饱和度指数不随天数变化，与PV无相关（r=0.11），且各品种间PV无显著差异（F=0.47，p=0.70），说明品种信号来自组成差异而非氧化程度。

3.7 跨基底比较
品种方差占比在油中达62%，籽仁中为10%；ANOVA的F值在油中为70~75，籽仁中为21~27；分类准确率两基底接近（约57%）。提取油是更干净的基底，但完整籽仁测量仍然可行，且具有非破坏性的实用优势。

3.8 负责品种判别的波数
逐波数ANOVA显示最强判别峰位于3010~3020、1255~1265和1656 cm^?1。PLSDA的变量重要性投影(VIP)分数因相关变量共享权重而分散，但仅用三个谱带指数即可达到74%的留一天准确率，高于全谱二阶导数的62%，进一步证明判别由不饱和度特征驱动。

讨论部分指出，三个独立标记一致、ASCA方差划分与多变量分析共同证明品种效应具有明确化学基础；提取油和籽仁结果一致，后者因基质成分复杂导致信号较弱，但油占籽仁约70%，脂质信号仍占主导。过氧化值实验说明在早期氧化水平下总双键含量几乎不变，因此拉曼指纹稳健。局限性包括：样品为商业混合批次，不能完全解耦果园和季节影响；所用为相对不饱和度指数而非绝对脂肪酸百分比；采样时间与测量session存在部分混杂；冷压步骤构成工作流程的一部分。未来方向是将拉曼指数与气相色谱参考值配对，建立定量预测模型和验证分类器。

结论：便携式拉曼光谱能够分辨山核桃油不饱和度的品种水平差异。在提取油中，四个品种形成可重复的不饱和度梯度（Desirable、Schley、Stuart、Moneymaker），三个谱带标记均显著（F=68~75；p<10^?20），ASCA将62%的光谱方差归因于品种（p<0.001），而储存条件和时间不显著。相同的排序在完整籽仁中非破坏性地恢复（F=21~27），表明测量不需要提取油。独立的过氧化值测定确认油在储存期间氧化，平均值在79天内从3.1升至17.5毫当量/千克，但拉曼指纹和品种排序保持稳定，因为这种早期氧化未反映在不饱和度谱带中。因此，品种特征反映每个品种内在的不饱和度，且在测试范围内不受储存氧化的影响。该测量快速、无需试剂、在便携式仪器上运行，为育种和加工中山核桃油质量的无损筛选提供了实用基础。研究人员将此作为可行性研究，其中拉曼指数的不饱和度意义建立在与参考化学的既定校准基础上，排序被视为特定于当前材料。将拉曼指数与气相色谱参考测量配对将能实现下一步，即开发预测模型，既包括从光谱直接估算油不饱和度和脂肪酸组成的定量校准，也包括识别未知样品品种的验证分类器。