作者简介: 陈月飞,女, 1989年生,中国农业大学理学院硕士研究生 e-mail: 949454851@qq.com
利用衰减全反射傅里叶变换红外光谱(ATR-FTIR)结合偏最小二乘法(PLS)对苏云金杆菌制剂中违禁添加的高效氯氟氰菊酯含量进行了快速定量检测。 在3种不同来源的苏云金杆菌制剂中加入不同质量的95%(w/w)高效氯氟氰菊酯原药, 配制成浓度为0.00%~5.00%的混合样品153个, 乙腈作为萃取剂, 选用较少萃取剂用量提高隐性成分在萃取液中浓度以增强红外吸收。 选择3种预处理方法: 平滑法(smooth)、 标准正态变换法(SNV)和多元散射校正法(MSC)以及6种变量选择算法: 无信息变量消除法(UVE)、 区间偏最小二乘法(iPLS)、 移动窗口偏最小二乘法(MWPLS)、 竞争性自适应重加权采样(CARS)、 自助软收缩法(BOSS)、 区间组合优化(ICO), 考察了不同预处理方法及变量选择方法对模型结果的影响。 其中, MSC预处理结合BOSS算法获得最优的模型效果: 交叉验证均方根误差(RMSECV)为1.18×10-3, 交叉验证集决定系数(
Rapid quantitative determination of the illegal addition of lambda-cyhalothrin in Bacillus thuringiensis preparations was carried out using attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) coupled with partial least squares (PLS) method. Three different sources of Bacillus thuringiensis preparations with different masses of 95% (w/w) lambda-cyhalothrin prodrug were added to prepare 153 mixed samples, concentrations ranging from 0.00% to 5.00%, acetonitrile was used as the extractant, and less extractant was used to increase the concentration of lambda-cyhalothrin in the extract to enhance the infrared absorption. Three pretreatment methods (smoothing, standard normal variation (SNV) and multiple scattering correction (MSC)) and six variable selection algorithms(uninformative variable elimination (UVE), interval partial least squares (iPLS), moving window partial least squares (MWPLS), competitive adaptive reweighted sampling (CARS), the bootstrapping soft shrinkage (BOSS) and interval combination optimization (ICO)) , were used to examine the effects of different pretreatment methods and variable selection methods on the model results. Among them, the MSC pretreatment method combined with the BOSS algorithm obtained the optimal model results, with this method, RMSECV=1.18×10-3,
高效氯氟氰菊酯(lambda-cyhalothrin)是新一代Ⅱ 型广谱高效拟除虫菊酯杀虫剂[1], 因活性高、 杀虫谱广、 速效性强被广泛使用。 然而, 其对水生生物具有高毒性[2], 在非靶生物体内与多种毒性缺陷(肝毒性、 肾毒性、 神经毒性和生殖毒性等)密切相关[3]。 苏云金杆菌(Bacillus thuringiensis, Bt)是一种昆虫病原细菌, 是我国目前应用范围较广的微生物杀虫剂, 对非靶标生物安全, 然而Bt制剂存在见效慢、 防效不高、 稳定性较差等问题[4]。 农业农村部近三年农药监督抽查结果显示, 杀虫剂中违规添加隐性成分情况严重, 其中苏云金杆菌制剂的隐性添加成分中, 有多个样品检出含有高效氯氟氰菊酯, 检出含量均在2.00%~5.00%之间[5, 6, 7]。 这些含有高效氯氟氰菊酯的Bt农药制剂的使用对环境和非靶标生物带来诸多不良影响, 迫切需要开发一种快速检测苏云金杆菌中违禁添加高效氯氟氰菊酯的方法。
常规的农药分析方法主要为色谱法和色谱-质谱联用法[8, 9, 10], 这类方法成本高、 检测周期长、 前处理要求高等问题, 无法适用于快速大量筛查和在线检测的需求。 红外光谱技术因成本更低、 绿色环保、 样品前处理简单、 分析速度快, 近年来在食品[11]、 医学[12]、 农业[13]等方面广泛应用。 在农药隐性添加成分检测领域, 同样有红外光谱技术的报道: 高鑫[14]等利用偏最小二乘法(partial least squares, PLS)建立了毒死蜱乳油中敌敌畏的定量模型并优化, Li[15]等采用中红外光谱(middle infrared, MIR)结合区间组合优化(interval combination optimization, ICO)和区间偏最小二乘(interval partial least-squares, iPLS)两种变量选择方法建立了啶虫脒乳油中氟虫腈含量的定量模型, 杜夏瑜[16]等利用红外衰减全反射技术(attenuated total reflection-Fourier transform-infrared, ATR-FT-IR)测定啶虫脒固体制剂中违禁添加氟虫腈含量, 建立固体样品及萃取液的定量模型, 韦芸[17]等通过不同预处理结合变量选择算法考察了苏云金杆菌粉剂中溴虫腈的含量, 建立最优定量模型。 以上研究结果均表明, 红外光谱结合化学计量学在不同农药剂型的违禁添加成分检测中有巨大的应用潜力。
目前, 国内外利用红外光谱定量测定农药制剂中拟除虫菊酯类违禁添加的报道较少。 本研究选取近年来违禁事件频繁涉及的生物农药苏云金杆菌制剂, 采用ATR-FT-IR结合化学计量学方法对苏云金杆菌可湿性粉剂中违禁添加高效氯氟氰菊酯含量进行定量分析。 采用Kennard-Stone(KS)法划分数据集, 以保证校正集和验证集样品按照空间距离分布均匀; 选用3种预处理方法: 平滑去噪算法(smoothing)、 标准正态变换法(standard normal variation, SNV)、 多元散射校正法(multivariate scattering correction, MSC), 6种变量选择算法: 无信息变量消除法(elimination of uninformative variables, UVE)、 区间偏最小二乘法(iPLS)、 移动窗口偏最小二乘法(moving window partial least squares, MWPLS)、 竞争性自适应重加权采样(competitive adaptive reweighted sampling, CARS), 自助软收缩法(bootstrapping soft shrinkage, BOSS), 区间组合优化(ICO), 考察了不同预处理方法和变量选择算法对PLS模型结果的影响。 通过对比模型的预测集决定系数
Nicolet Summit PRO型傅里叶变换红外光谱仪(美国Thermo Fisher公司, 配有Omni iD5-ATR金刚石附件); SCILOGEX MX-S可调式混匀仪(上海科雅生物科技有限公司); TDL-60B型台式高速离心机(上海安亭科学仪器厂)。
苏云金杆菌可湿性粉剂(16 000 IU· mg-1, 康欣生物科技有限公司; 农药登记证号: PD86109-33); 苏云金杆菌可湿性粉剂(32 000 IU· mg-1, 武汉科诺生物科技股份有限公司; 农药登记证号: PD20084969); 苏云金杆菌可湿性粉剂(16 000 IU· mg-1, 乳山韩威生物科技有限公司; 农药登记证号: PD20096222); 高效氯氟氰菊酯原药(95%, w/w, 山东省联合农药工业有限公司); 乙腈(分析纯, 西陇化工股份有限公司); 无水乙醇(分析纯, 昆山金城试剂有限公司)。 其他试剂均为国产分析纯试剂。
1.2.1 样品制备
在3种不同来源的1 g苏云金杆菌可湿性粉剂中分别加入不同质量的95%(w/w)高效氯氟氰菊酯原药, 研细后, 配制成3个浓度系列(0.00%~5.00%, 质量浓度, w/w)高效氯氟氰菊酯固体样品, 每个系列51个样品, 共153个, 样品质量百分数均值为2.50%。 分别用4种溶剂(乙腈、 二甲亚砜、 N, N-二甲基酰胺和四氯化碳)溶解适量高效氯氟氰菊酯原药, 其中乙腈对高效氯氟氰菊酯溶解效果较好, 且溶解后红外光谱图基线较平稳, 故选用其作为混合固体样品的萃取溶剂, 对上述153个样品进行萃取。 取不同高效氯氟氰菊酯添加浓度的样品各1.00 g, 加入0.8 mL乙腈, 置于混匀仪上涡旋振荡2 min, 再以4 000 r· min-1离心5 min, 上清液保存备用。
1.2.2 光谱采集
使用傅里叶变换红外光谱仪(ATR 附件)采集153个样品上清液的中红外光谱图, 每个样品采集 3 条光谱。 仪器参数: 扫描范围4 000~600 cm-1, 分辨率4 cm-1, 扫描次数32次。
1.2.3 数据处理
数据处理软件选用9.7版OMNIC、 10.3版The Unscrambler X、 R2019b版Matlab、 2019版Excel和2018版Origin。
以高效氯氟氰菊酯在苏云金杆菌制剂中的含量(质量百分数)作为分析指标, 对样品上清液的光谱数据建立定量分析模型。 所用光谱为每个样品重复采集的3条光谱的直接平均值, 采用3种预处理方法(平滑去噪、 SNV、 MSC)结合6种变量选择方法(UVE、 IPLS、 MWPLS、 CARS、 BOSS、 ICO)优化定量模型, 其中UVE 是基于波长稳定性的波长选择算法, IPLS、 MWPLS是基于波长区间的波长选择算法, CARS、 BOSS可在全光谱中选择出最佳组合波长, ICO 也是基于model population analysis(MPA)框架开发的区间优化选择方法。 光谱数据导入和矩阵整理在Unscrambler X中进行, 光谱平均、 预处理、 变量选择及PLS模型建立在Matlab中进行。
1.2.4 模型评价指标
由于通过波长选择选取部分重要变量进行建模, 可避免模型出现过拟合风险, 提高模型预测性能, 因此, 除采用
图1为高效氯氟氰菊酯和苏云金杆菌制剂的红外光谱图。 高效氯氟氰菊酯原药在1 720 cm-1处有较强的尖锐吸收峰, 为羰基C═O伸缩振动峰; 1 584、 1 483和1 445 cm-1处有尖锐的吸收峰, 为苯环骨架C— C伸缩振动峰; 1 408~1 378 cm-1处峰为甲基C— H弯曲振动; 1 293~1 270 cm-1处为芳环上C— H面内弯曲振动, 1 232~1 049 cm-1处峰为C— O— C伸缩振动; C— F伸缩振动也在1 400~1 200 cm-1范围内; 900~650 cm-1为芳环C— H面外弯曲振动。 苏云金杆菌可湿性粉剂在1 160~850 cm-1处有一个宽的强吸收峰。
![]() | 图1 (a)高效氯氟氰菊酯原药红外光谱; (b)苏云金杆菌制剂红外光谱Fig.1 (a) ATR-FTIR spectra of 95% lambda-cyhalothrin; (b) ATR-FTIR spectra of Bacillus thuringiensis preparation |
图2为配制样品的上清液的红外光谱图, 可见随高效氯氟氰菊酯浓度的变化, 1 720、 1 350~1 050, 959、 900~780和700~600 cm-1几处的吸收峰值出现了明显差异。
图3为苏云金杆菌在乙腈中的谱图, 与纯乙腈的谱图差异很小, 添加了5%(w/w)高效氯氟氰菊酯的苏云金杆菌在乙腈中的谱图在1 720、 1 350~1 050、 959、 900~780和700~600 cm-1几处明显有吸收峰, 与高效氯氟氰菊酯的强吸收峰相符, 因此, 采用乙腈作为萃取液是可行的。
对153个样品对应的153条平均光谱, 以苏云金杆菌可湿性粉剂中高效氯氟氰菊酯含量为化学值, 采用KS算法划分样品校正集与验证集(校正集∶ 验证集=115∶ 38), 采用PLS法建立苏云金杆菌制剂中高效氯氟氰菊酯含量的定量模型。 具有大量波长变量和相对较少样本的光谱数据情况下, 建模具有过度拟合的高风险, 这会导致多变量校准模型获得不良或低效的预测结果; 变量选择可以用于消除无关或无信息变量, 简化校准建模, 提高模型的预测性能。 因此, 采用3种预处理方法(平滑去噪、 SNV、 MSC)结合6种变量选择方法(UVE、 IPLS、 MWPLS、 CARS、 BOSS、 ICO)对模型进行优化。 由表1可见, MSC预处理方法与不同变量选择方法结合都能获得较好的预测结果, MSC主要用于消除颗粒分布不均匀及颗粒大小产生的散射影响, 适用于固体漫反射和浆状物反射光谱。 实验采用乙腈萃取可湿性粉剂中高效氯氟氰菊酯, 离心后吸取上清液, 肉眼观察无法保证上清液中完全不含不均匀微粒, 通过MSC方法可消除颗粒导致的散射影响。 BOSS变量选择方法与ICO变量选择方法与不同预处理结合均有较好的结果, 其中, 最优模型为MSC预处理结合BOSS变量方法选择的变量建立的模型, 其参数为: RMSECV=1.18× 10-3,
![]() | 表1 高效氯氟氰菊酯模型参数及优化结果 Table 1 Parameters and optimization results of PLS model for lambda-cyhalothrin |
![]() | 表2 MSC-BOSS模型高效氯氟氰菊酯样品校正集与验证集分布 Table 2 Distribution of calibration and validation sets of lambda-cyhalothrin samples for MSC-BOSS modeling |
![]() | 图4 MSC结合BOSS的定量模型效果图 (a): MSC结合BOSS的校正集数据模型; (b): MSC结合BOSS的验证集数据模型Fig.4 Quantitative modeling results of MSC combined with BOSS===(a): Calibration set data; (b): Validation set data |
MSC结合6种变量选择算法选出的变量范围如图5所示。
除MWPLS在约1 600 cm-1附近选取变量外, 其余各变量选择方法均在600~800 cm-1之间选取了变量。 与高效氯氟氰菊酯结构中芳环C— H振动有关, 其中, 与各预处理方法结合模型结果均较好的BOSS与ICO变量选择方法选取的波数变量区间均包含2 000、 1 720、 1 230、 1 000和630 cm-1附近处, 除2 000 cm-1处外, 其他位置随高效氯氟氰菊酯含量增高而红外吸收增强, 体现出高效氯氟氰菊酯红外光谱中C═O伸缩振动、 C— O— C伸缩振动、 C— F伸缩振动和芳环C— H面外弯曲振动, 其中BOSS选择方法仅选取30个变量数即能达到RMSEP=1.01× 10-3,
以KS分集法选出的38个样品作为验证集,
![]() | 表3 高效氯氟氰菊酯PLS模型验证结果 Table 3 PLS model validation results for lambda-cyhalothrin |
采用衰减全反射傅里叶变换红外光谱结合化学计量学方法测定了苏云金杆菌可湿性粉剂中高效氯氟氰菊酯的含量。 选用3种光谱预处理方法(平滑、 SNV、 MSC)和6种变量选择算法(UVE、 iPLS、 MWPLS、 CARS、 BOSS、 ICO)考察预处理和变量选择对模型效果的影响。 其中, MSC预处理结合BOSS算法获得最优的模型效果, 且该模型仅选取30个变量数即获得稳定的预测效果, 样品平均相对误差为4.44%, 高效氯氟氰菊酯含量在2.00%~5.00%的平均相对误差仅2.64%, 对近年来农药监督抽查苏云金杆菌制剂中高效氯氟氰菊酯违规添加含量的检测适用; 低含量样品预测效果也较好, 高效氯氟氰菊酯含量在0.10%~1.90%的样品平均相对误差为7.08%, 该方法也可用于较低含量高效氯氟氰菊酯违规添加的现场快速筛查。 该方法为苏云金杆菌可湿性粉剂中高效氯氟氰菊酯含量的快速测定提供了理论基础, 但实际工作中农药制剂中隐性添加成分可能较为复杂, 仍需结合实际工作进一步探索。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|