虞龙涛^,1, 杨何妍^,1, 苏宇晨², 颜伟玉¹, 吴小波^,11江西农业大学蜜蜂研究所/江西省蜜蜂生物学与饲养重点实验室,南昌 330045
²江西农业大学动物科学技术学院,南昌 330045

The Effect of Flumethrin on Metabolism of Worker Larvae of Apis mellifera with LC-MS Technique

YU LongTao^,1, YANG HeYan^,1, SU YuChen², YAN WeiYu¹, WU XiaoBo^,11Honeybee Research Institute, Jiangxi Agricultural University/Jiangxi Province Key Laboratory of Honeybee Biology and Beekeeping, Nanchang 330045
²College of Animal Science and Technology, Jiangxi Agricultural University, Nanchang 330045

通讯作者: 吴小波,E-mail：wuxiaobo21@163.com

虞龙涛和杨何妍为同等贡献作者。
责任编辑: 岳梅
收稿日期:2020-09-20接受日期:2020-11-9网络出版日期:2021-06-16

基金资助:

国家自然科学基金.31760714 江西省主要学科学术和技术带头人项目.20194BCJ22007

Received:2020-09-20Accepted:2020-11-9Online:2021-06-16
作者简介 About authors
虞龙涛,E-mail：1208789722@qq.com。

杨何妍,E-mail：1035361586@qq.com。






摘要
【目的】氟氯苯氰菊酯（flumethrin）属于第二代拟除虫菊酯类杀虫、杀螨剂,应用于蜂群除螨。由于杀螨剂具有一定的毒性,在杀死螨虫的同时也会威胁到蜜蜂的健康。本研究采用液相色谱-质谱联用技术（LC-MS）检测氟氯苯氰菊酯处理后的西方蜜蜂工蜂幼虫代谢物,筛选差异代谢物,并分析涉及到的代谢通路,研究氟氯苯氰菊酯对蜜蜂的毒理作用,为蜜蜂饲养提供参考。【方法】控制蜂王在空的工蜂巢脾上产卵12 h,并将产卵区域分为4组,从第5天开始分别给各组小幼虫饲喂含氟氯苯氰菊酯的糖水（0、0.5、5、50 mg·kg^-1）,剂量从第5天至第8天逐日递增（1.5、2、2.5、3 μL）,第9天采集幼虫的淋巴液。运用代谢组学中的LC-MS技术分析蜜蜂幼虫的代谢物,结合主成分分析（PCA）、偏最小二乘法-判别分析（PLS-DA）,筛选差异显著的代谢物,并对氟氯苯氰菊酯处理组和对照组间共同差异代谢物进行代谢通路分析。【结果】0.5 mg·kg^-1组与对照组相比,共有差异代谢物190种,共鉴定87种;5 mg·kg^-1组与对照组相比,共有差异代谢物275种,共鉴定97种;50 mg·kg^-1组与对照组相比,共有差异代谢物275种,共鉴定131种。从鉴定的差异代谢物中筛选到29种氟氯苯氰菊酯处理组共同差异代谢物,其中16种代谢物上调,12种代谢物下调,1种代谢物在0.5和50 mg·kg^-1组中下调,在5 mg·kg^-1组中上调。这些差异代谢物包括核糖、嘌呤及其衍生物、脂肪酸及其偶联物等。经代谢通路富集分析,发现差异显著（P＜0.05）的代谢通路包括氨基糖和核苷酸糖代谢、药物代谢-其他酶类、α-亚麻酸代谢等途径。【结论】LC-MS技术能够有效地分析氟氯苯氰菊酯影响下蜜蜂幼虫代谢物的变化,氟氯苯氰菊酯会引起蜜蜂幼虫体内尿苷二磷酸-N-乙酰葡糖胺、硫唑嘌呤、愈伤酸、9-羟基壬酸、氢过氧化亚油酸等物质含量异常,这些差异代谢物的变化证实了氟氯苯氰菊酯引起蜜蜂体内各种物质代谢紊乱,对这些代谢过程分析可进一步阐释蜜蜂代谢有毒化合物的机制,并为杀螨剂及其他有毒化合物对蜜蜂的胁迫提供理论依据。
关键词： 西方蜜蜂;LC-MS;代谢组学;氟氯苯氰菊酯;工蜂幼虫

Abstract
【Objective】Flumethrin belongs to the second generation pyrethroid insecticides and acaricides, which is used for the control of honeybee mites. Because of the toxicity of acaricides, it can not only kill mites, but also threat the health of honeybees. The different metabolites of Apis mellifera worker larvae treated with different concentrations of flumethrin were tested using liquid chromatography-mass spectrometry (LC-MS) technology and the metabolic pathways involved of different metabolites were analyzed, so as to explore the toxicological effect of flumethrin on honeybees and provide references for scientific using in beekeeping.【Method】The queen was controlled to lay eggs on an empty worker frame for 12 h, and the spawning area was divided into four groups. From the 5th day, the small larvae of each group were fed with sugar water containing different concentrations of flumethrin (0, 0.5, 5, 50 mg·kg^-1), the dose was increased daily from day 5 to day 8 (1.5, 2, 2.5, 3 μL), and lymph fluid from larvae was collected on day 9. The metabolites of A. mellifera larvae were analyzed by LC-MS and the metabolites with significant difference were screened by principal component analysis (PCA) and partial least squares-discriminant analysis (PLS-DA), while the metabolic pathways of common differential metabolites in flumethrin treatment groups were analyzed.【Result】Compared with the control group, a total of 190 different metabolites were found and 87 types were identified in 0.5 mg·kg^-1 group, and a total of 275 different metabolites were identified and 97 types were identified in 5 mg·kg^-1 group, while there were a total of 275 different metabolites and a total of 131 species were identified in 50 mg·kg^-1 group. Meanwhile, 29 common differential metabolites in treatment groups were screened, of which 16 metabolites were up-regulated, 12 metabolites were down-regulated, and 1 metabolite was down-regulated in 0.5 and 50 mg·kg^-1 groups while it was up-regulated in 5 mg·kg^-1 group. These differential metabolites include ribose, purine and its derivatives, fatty acids with their conjugates. After enrichment analysis of metabolic pathways, significant differences (P＜0.05) were found in the metabolic pathways, which include amino sugar and nucleotide sugar metabolism, drug metabolism-other enzymes, α-linolenic acid metabolism and other pathways.【Conclusion】The LC-MS technology can effectively analyze the changes of metabolites in the honeybee larvae treated with flumethrin, and flumethrin can cause the contents of UDP-N-acetylglucosamine, azathioprine, traumatic acid, 9-oxononanoic acid and 13(s)-HPODE abnormal in honeybee larvae. The changes of these different metabolites confirm that flumethrin causes various substance metabolism disorders in honeybees. The analysis of these metabolic processes can further explain the mechanism of honeybees metabolizing toxic compounds, and provide a theoretical basis for the stress of acaricides and other toxic compounds on honeybees.
Keywords：Apis mellifera;LC-MS;metabolomics;flumethrin;worker larva


PDF (1581KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
虞龙涛, 杨何妍, 苏宇晨, 颜伟玉, 吴小波. 基于LC-MS技术研究氟氯苯氰菊酯对西方蜜蜂工蜂幼虫代谢的影响[J]. 中国农业科学, 2021, 54(12): 2689-2698 doi:10.3864/j.issn.0578-1752.2021.12.018
YU LongTao, YANG HeYan, SU YuChen, YAN WeiYu, WU XiaoBo. The Effect of Flumethrin on Metabolism of Worker Larvae of Apis mellifera with LC-MS Technique[J]. Scientia Acricultura Sinica, 2021, 54(12): 2689-2698 doi:10.3864/j.issn.0578-1752.2021.12.018


开放科学（资源服务）标识码（OSID）：

0 引言

【研究意义】蜜蜂是一种非常重要的授粉昆虫,许多农作物都依赖蜜蜂授粉,蜜蜂给全世界农业创造了巨大的经济价值^[1]。然而蜜蜂也受到许多危害,如病毒、寄生虫、螨虫等,其中蜂螨作为主要危害,给全球蜜蜂产业造成严重的损失,目前治理蜂螨最为简便的方法是使用杀螨剂。氟氯苯氰菊酯（flumethrin）属于第二代拟除虫菊酯类杀虫、杀螨剂,使用时将其挂入蜂箱巢脾中,通过蜜蜂活动时接触药条并互相传递,使药发挥作用^[2],由于杀螨剂具有一定的毒性,在杀死螨虫的同时也会威胁到蜜蜂的健康。目前,许多研究表明杀螨剂会给蜜蜂生理和行为带来不可逆转的影响^[3,4,5,6],但是关于氟氯苯氰菊酯影响蜜蜂的健康状况机制研究则相对滞后。因此,测定氟氯苯氰菊酯影响下工蜂幼虫代谢物的变化,并探讨涉及的代谢通路,可为进一步阐释杀螨剂影响蜜蜂发育的机制提供参考。【前人研究进展】氟氯苯氰菊酯作为一种毒性药物,直接暴露于氟氯苯氰菊酯下的蜜蜂健康状况会受到其明显的影响。有研究表明,亚致死剂量的氟氯苯氰菊酯会引起西方蜜蜂（Apis mellifera）机体紊乱^[3,4,5],QI等用1、10、100 μg·L^-1氟氯苯氰菊酯处理西方蜜蜂后发现,氟氯苯氰菊酯在幼虫变态和成虫羽化过程中都有明显的致死作用。在幼虫期连续暴露后,幼虫的抗氧化酶（SOD和CAT）、脂质过氧化（MDA、LPO和POD）和解毒酶（GSH、GST和GR）发生了显著变化^[3]。而且经0.01、0.1和1.0 mg·L^-1氟氯苯氰菊酯慢性染毒14 d后,蜜蜂中肠出现明显的抗氧化反应、解毒作用、免疫反应和细胞凋亡^[4]。此外,牛新月^[5]研究发现,给西方蜜蜂幼虫饲喂含氟氯苯氰菊酯的糖水,蜜蜂幼虫的发育受到影响,并且刚出房的蜜蜂出现畸形,蜜蜂解毒关键酶的含量以及活性也有明显差异。氟氯苯氰菊酯还会影响中华蜜蜂（Apis cernana）的寿命以及死亡率,TAN等^[6]研究发现,摄入亚致死剂量含氟氯苯氰菊酯的糖水会显著影响中华蜜蜂的寿命并影响到蜜蜂的采集行为。江武军等^[7]采用10、100 mg·kg^-1的氟氯苯氰菊酯糖水饲喂中华蜜蜂,发现氟氯苯氰菊酯处理组蜜蜂中肠组织AcCYP9E2表达量显著高于对照组,CYP9E2参与了外源物质代谢解毒过程,表明氟氯苯氰菊酯会影响蜜蜂解毒功能。【本研究切入点】多数研究表明,氟氯苯氰菊酯对蜜蜂幼虫以及成虫带来不可逆转的影响,但蜜蜂对氟氯苯氰菊酯代谢的研究相对滞后,有关氟氯苯氰菊酯对蜜蜂血淋巴代谢物影响的研究未见报道,蜜蜂代谢氟氯苯氰菊酯的机理尚不明确。【拟解决的关键问题】采用液相色谱-质谱联用技术（liquid chromatography-mass spectrometry,LC-MS）检测氟氯苯氰菊酯处理后蜜蜂幼虫代谢物,并筛选差异代谢物,分析涉及到的代谢通路,结合代谢通道数据库（Kyoto encyclopedia of genes and genomes,KEGG）,探究氟氯苯氰菊酯影响下蜜蜂代谢物的变化规律,为进一步阐明蜜蜂代谢氟氯苯氰菊酯的机制提供理论依据。

1 材料与方法

1.1 供试蜜蜂

利用群势相近的3群西方蜜蜂,饲养于江西农业大学蜜蜂研究所内,试验开展时间为2019年。

1.2 药物配置

根据氟氯苯氰菊酯在蜂产品中的残留情况^[8,9,10,11],结合氟氯苯氰菊酯对蜜蜂的亚致死剂量浓度^[6,12],氟氯苯氰菊酯处理浓度为0.5、5、50 mg·kg^-1。用50%糖水稀释氟氯苯氰菊酯至0.5、5、50 mg·kg^-1,以饲喂50%糖水作为对照组。

1.3 限王产卵及淋巴液的采集

将蜂王限制在巢脾上产卵12 h,随后将蜂王产卵区域分为4组,并将产卵的巢脾放在继箱上孵化,第5天开始,在巢脾不同区域每天用移液枪分别给小幼虫饲喂含有氟氯苯氰菊酯的糖水,饲喂量从第5天至第8天逐日递增（1.5、2、2.5、3 μL）,第9天将幼虫从脾中挑出,置于经高温灭菌的培养皿内,用高压双蒸水冲洗3遍直到其体表被冲洗干净,之后用滤纸吸干幼虫体表水分。将幼虫置于PE手套上,用剪刀在幼虫的头部剪一个小口,用移液枪吸取流出的血淋巴,转移至1.5 mL的EP管。每个EP管收集的蜜蜂幼虫血淋巴体积为200 μL。然后14 000×g、4℃离心15 min。离心后,取120—150 μL上清于新的1.5 mL EP管中为1个样品,并快速置于液氮中冻存,每组6个重复。

1.4 代谢物的提取

将样本在4℃条件下缓慢解冻后,取100 μL置于96孔板中,加入300 μL提取液（甲醇﹕乙腈=2﹕1,-20℃预冷）、10 μL内标1、10 μL内标2,涡旋混匀1 min,-20℃静置2 h后,4℃,4 000×g离心20 min。离心后取300 μL上清,置于冷冻真空浓缩仪抽干后,加入150 μL复溶液（甲醇﹕水=1﹕1）进行复溶,涡旋振荡1 min,4℃,7 000×g离心30 min,取上清置于上样瓶中。每个样本的上清液各取10 μL混合成QC质控样本,用于评估LC-MS分析过程的重复性和稳定性,通过QC样本检测的重复性对数据质量进行评估,内容包括QC样本的基峰离子流图（base peak chromatogram,BPC）、主成分分析（principle component analysis,PCA）、提取到的峰数量和峰相应强度差异。

1.5 色谱条件

本研究期望在正、负离子模式下寻找3种剂量氟氯苯氰菊酯处理组与对照组共同的差异代谢物,但负离子模式下未检测出共同的差异代谢物,因此只分析正离子模式下的结果。所使用的色谱柱为BEH C18色谱柱（1.7 μm,2.1×100 mm,Waters,USA）。流动相为含0.1%甲酸的水溶液（A液）和含0.1%甲酸的100%甲醇（B液）。采用以下梯度进行洗脱：0—1 min,2% B液;1—9 min,2%—98% B液;9—12 min,98% B液;12—12.1 min,98% B液—2% B液;12.1—15 min,2% B液。流速为0.35 mL·min^-1,柱温45℃,进样量为5 μL。

1.6 质谱条件

利用Q Exactive质谱仪（Thermo Fisher Scientific,USA）进行一级、二级质谱数据采集。质谱扫描质核比范围为70—1 050,一级分辨率为70 000,AGC为3e⁶,最大注入时间为100 ms。按照母离子强度,选择Top3进行碎裂,采集二级信息,二级分辨率为17 500,AGC为1e⁵,最大注入时间为50 ms,碎裂能量设置为20、40、60 eV。

离子源参数设置：鞘气流速为40 L·min^-1,辅助气流速为10 arb,喷雾电压正离子模式为3 800 V,负离子模式为3 200 V,离子传输管温度为320℃,辅助气加热温度为350℃。

1.7 数据分析

采用LC-MS/MS技术进行非靶向代谢组学分析,使用Compound Discoverer 3.0（Thermo Fisher Scientific,USA）软件进行LC-MS/MS数据处理,完成峰提取、峰对齐和化合物鉴定。利用代谢组信息分析流程进行数据预处理、统计分析及代谢物分类注释和功能注释。通过主成分分析（PCA）对检测到的代谢物的原始数据进行降维,以此分析该数据集中所观测变量的分组、趋势（样本组内和组间相似性和差异性）及离群值（是否存在异常样本）。使用偏最小二乘法-判别分析（partial least squares-discriminant analysis,PLS-DA）分析模型前两个主成分的VIP值,结合单变量分析获得的差异变化倍数（fold change,FC）和T检验（Student’s t test）的结果来筛选差异代谢物。

差异代谢物筛选条件：PLS-DA模型前两个主成分的VIP≥1;FC≥1.2或≤0.83;q-value＜0.05。筛选得到的所有代谢物的VIP值均＞1,表明该变量对样本类别的区分有显著作用;FC为该差异代谢物在两组中的平均含量的比值,本研究中的FC为处理组平均含量/对照组平均含量。对数据进行T检验得到P-value,P-value进行错误发现率（false discovery rate,FDR）校正得到q-value。

筛选出氟氯苯氰菊酯处理组与对照组之间差异代谢物,并根据差异代谢物中性质量数以及化学式等信息筛选处理组间的共同差异代谢物。

基于HMDB人类代谢组学数据库（The Human Metabolome Database）和KEGG京都基因与基因组百科全书,对差异代谢物进行分类注释。根据差异代谢物的富集情况确定差异代谢物代谢通路。

2 结果

2.1 样本总离子流图

如图1,BPC（基峰离子流图）是将每个时间点质谱图中最强的离子连续描绘得到的图谱,将处理组与对照组总离子流图进行重叠,谱图重叠良好,保留时间和峰响应强度均波动小,表明仪器在整个样本检测分析过程中状态良好,信号稳定,数据可信,可进行下一步分析。

图1

新窗口打开|下载原图ZIP|生成PPT
图1对照组与氟氯苯氰菊酯处理组的BPC（基峰离子流图）重叠图

Fig. 1Overlapping BPC (base peak chromatogram) plot of the control group and flumflethrin-treated groups

2.2 主成分分析（PCA）

不同剂量氟氯苯氰菊酯处理组与对照组PCA分析得分图如图2所示,图中每个点代表一个样本,点与点的距离代表样本之间的相似性,椭圆为95%的置信区间。由图2可知,对照组与低剂量组样本分布在PC2的上下两侧,对照组与中剂量组样本主要分布在PC1左右两侧,对照组与高剂量组样本主要分布在PC2的上下两侧,说明对照组与氟氯苯氰菊酯处理组存在组间差异。

图2

新窗口打开|下载原图ZIP|生成PPT
图2对照组与氟氯苯氰菊酯处理组间蜜蜂代谢物 PCA得分图

A：0.5 mg·kg^-1组与对照组之间蜜蜂代谢物PCA得分图PCA scores of honeybee metabolites between the 0.5 mg·kg^-1 group and control group;B：5 mg·kg^-1组与对照组之间蜜蜂代谢物PCA得分图PCA scores of honeybee metabolites between 5 mg·kg^-1 group and control group;C：50 mg·kg^-1组与对照组之间蜜蜂代谢物PCA得分图PCA scores of honeybee metabolites between the 50 mg·kg^-1 group and control group。横轴为第一主成分,纵轴为第二主成分。图3同The horizontal-axis shows the first principal component, and the vertical-axis shows the second principal component. The same as Fig. 3
Fig. 2PCA scores plots of honeybee metabolites between the control group and flumflethrin-treated groups

2.3 偏最小二乘判别分析（PLS-DA）

处理组与对照组之间的PLS-DA 得分图见图3。从第一主成分看,处理组分布在负轴位置,对照组主要分布正轴位置;从第二主成分看,处理组主要分布在负轴位置上,对照组主要分布于正负轴两侧。表明处理组与对照组的代谢物分离较好,处理组的每个样品的代谢物均与对照组有明显差异。

图3

新窗口打开|下载原图ZIP|生成PPT
图3对照组与氟氯苯氰菊酯处理组间蜜蜂代谢物PLS-DA得分图

Fig. 3PLS-DA scores plots of honeybee metabolites between the control group and flumflethrin-treated groups

2.4 PLS-DA分析模型响应排序检验

PLS-DA模型的响应排序检验图验证结果如图4所示。结果表明,对照组与试验组代谢物R²和Q²均＞0.5且其比值接近1,说明 PLS-DA 模型能够较好地预测处理组与试验组代谢物的差异。Q²回归线始终在R²的下方,且Q²回归线的斜率较大,其与纵坐标的交点均在负半轴,说明模型质量较好,没有出现过于拟合的状况,结果稳定可靠。

图4

新窗口打开|下载原图ZIP|生成PPT
图4对照组与氟氯苯氰菊酯处理组间PLS-DA分析模型的响应排序检验图

R²：模型的可解释率,可评估模型的拟合能力Explainable rate of the model, which can be used to assess the fit ability of the model;Q²：模型预测率,可评估模型的预测能力Prediction rate of the model, which can be used to assess the predictive power of the model
Fig. 4Plots of response ranking test for PLS-DA analysis model between the control group and flumethrin-treated groups

2.5 差异代谢物鉴定

利用代谢物的分子质量和二级质谱信息在KEGG、HMDB数据库搜索比对,0.5 mg·kg^-1组与对照组相比,共有差异代谢物190种,共鉴定87种;5 mg·kg^-1组与对照组相比,共有差异代谢物275种,共鉴定97种;50 mg·kg^-1 组与对照组相比,共有差异代谢物275种,共鉴定131种（表1）。从鉴定的差异代谢物中筛选到29种共同的差异代谢物（表2）,其中16种代谢物上调,12种代谢物下调,而1种代谢物在50 mg·kg^-1组和0.5 mg·kg^-1组中均下调,在5 mg·kg^-1 组中表达上调。差异代谢物包括甲苯、核糖、氨基酸及其类似物、脂肪酸及其偶联物、嘌呤及其衍生物等。

Table 1
表1
表1差异代谢物及鉴定结果
Table 1Differential metabolites and identification results

组别 Group	差异代谢物总数 Total number of differential metabolites	差异代谢物鉴定数目 Identified number of differential metabolites	Up	Down
0.5 mg·kg-1-Control	190	87	96	94
5 mg·kg-1-Control	275	97	209	66
50 mg·kg-1-Control	275	131	147	128

0.5 mg·kg^-1-Control：0.5 mg·kg^-1组与对照组之间的差异代谢物The differential metabolites between the 0.5 mg·kg^-1 and control group;5 mg·kg^-1-Control：5 mg·kg^-1组与对照组之间的差异代谢物The differential metabolites between the 5 mg·kg^-1 and control group;50 mg·kg^-1-Control：50 mg·kg^-1组与对照组之间的差异代谢物The differential metabolites between the 50 mg·kg^-1 and control group。表2同The same as Table 2

新窗口打开|下载CSV

Table 2
表2
表2对照组与氟氯苯氰菊酯处理组之间差异代谢物
Table 2Differential metabolites between the control group and flumethrin-treated groups

化合物名称 Compound name	50 mg·kg-1-Control				5 mg·kg-1-Control				0.5 mg·kg-1-Control
	VIP	FC	P	趋势 Trend	VIP	FC	P	趋势 Trend	VIP	FC	P	趋势 Trend
尿苷二磷酸-N-乙酰葡糖胺 UDP-N-acetylglucosamine	2.1451	3.7838	0.0175	↑	1.0349	2.0512	0.0494	↑	1.4664	2.4830	0.0277	↑
硫唑嘌呤Azathioprine	3.5178	22.9131	0.0038	↑	2.2240	6.6142	0.0216	↑	2.8184	12.1241	0.0135	↑
反玉米素Trans-zeatin	3.4573	27.3706	0.0012	↑	2.1436	10.783	0.0109	↑	3.5919	29.5094	0	↑
焦炭酸二乙酯Diethylpyrocarbonate	1.8239	2.6369	0.0126	↑	1.3799	2.1616	0.0312	↑	1.8167	2.6558	0.0272	↑
阿斯巴甜 Methyl alpha-aspartylphenylalaninate	2.4076	6.6251	0.0210	↑	1.5393	3.6759	0.0056	↑	2.5201	7.6114	0.0272	↑
草完隆Noruron	2.5040	5.7411	0.0237	↑	1.8949	5.3422	0.0312	↑	2.4082	5.7124	0.0272	↑
四聚乙醛Ethanal tetramer	3.1358	10.8657	0.0022	↑	1.6834	4.3428	0.0168	↑	1.9323	4.0922	0.0240	↑
扎莫特罗Xamoterol	1.7021	2.5619	0.0303	↑	1.6402	2.9414	0.0271	↑	2.3289	4.8191	0.0272	↑
二甲弗林Dimetofrine	2.0830	5.3298	0.0337	↑	1.9187	6.1752	0.0216	↑	2.3776	6.4440	0.0434	↑
细辛脑Asarone	3.5178	22.9131	0.0038	↑	2.2168	7.7327	0.0301	↑	2.4104	5.5813	0.0399	↑
异丙酚Propofol	1.9412	3.3387	0.0169	↑	1.7439	3.4319	0.0121	↑	2.2298	4.1785	0.0135	↑
芹黄素Apigetrin	2.7393	7.0597	0	↑	2.1954	5.8356	0.008	↑	2.2192	5.1848	0.0070	↑
4,4'-二硝基联苄4,4'-dinitrobibenzyl	3.5924	36.459	0.0024	↑	1.7819	5.2543	0.0073	↑	3.1705	15.6981	0.0019	↑
早熟素Precocene ii	2.8581	11.2521	0.0012	↑	1.0312	2.0035	0.0439	↑	2.1354	3.8515	0.0428	↑
贝那普利拉Benazeprilat	1.3429	2.0242	0.0438	↑	1.8099	6.4223	0.0188	↑	2.6353	10.2252	0.0135	↑
氟扎可特Azacortid	2.6960	12.7132	0.0305	↑	2.4081	15.5769	0.0188	↑	3.5251	35.6758	0.0048	↑
甲苯Toluene	1.5881	0.4467	0.0226	↓	1.4841	0.2876	0.0301	↓	1.9802	0.2721	0.0235	↓
愈伤酸Traumatic acid	2.4552	0.1909	0.0021	↓	1.5634	0.2683	0.0121	↓	2.0103	0.2575	0.0079	↓
9-羟基壬酸9-oxononanoic acid	2.9634	0.0980	0.0006	↓	1.3122	0.2638	0.0139	↓	2.7705	0.1129	0.0015	↓
普瑞巴林Pregabalin	1.6876	0.4322	0.0028	↓	1.2720	0.4581	0.0178	↓	1.0768	0.5848	0.0389	↓
阿伐那非Avanafil	1.9936	0.3237	0.0253	↓	1.1258	0.4234	0.0347	↓	1.7236	0.3553	0.0235	↓
1-壬酸1-nonanoic acid	2.6990	0.1060	0.0104	↓	2.0052	0.1232	0.0238	↓	2.5335	0.1151	0.0260	↓
匹伐加宾Pivagabine	2.8956	0.1266	0.0057	↓	2.0864	0.1468	0.0560	↓	2.1990	0.1935	0.0135	↓
美噻吨Metixene	2.6297	0.1776	0.0012	↓	1.1710	0.3815	0.0271	↓	1.4653	0.3801	0.0277	↓
樟脑(+/-)-camphor	2.2568	0.2106	0.0169	↓	1.4478	0.2932	0.0247	↓	2.2074	0.2177	0.0216	↓
扎莫特罗Xamoterol [usan:ban:inn]	1.8083	0.2538	0.0229	↓	2.1189	0.1402	0.0056	↓	1.8550	0.2617	0.0235	↓
二甲福林Dimetofrine	2.5825	0.1265	0.0139	↓	2.1321	0.1535	0.0238	↓	2.6898	0.1118	0.0199	↓
苯乙酮Acetophenone	1.7796	0.3617	0.0213	↓	1.8288	0.2989	0.0139	↓	1.5460	0.4091	0.0387	↓
氢过氧化亚油酸13(s)-HPODE	1.1855	0.5253	0.0414	↓	1.6963	4.7824	0.0331	↑	1.5777	0.4183	0.0235	↓

VIP：变量影响重要性因子,该值越大,表明该成分对模型的贡献度越大Variable importance in the projection, the larger the value, the greater the contribution of the component to the model;FC：变化倍数,该值反映了处理组与对照组的差异Fold change, which reflects the difference between the treated groups and control group。趋势：“↑”,处理组与对照组相比表达上调; “↓”,处理组与对照组相比表达下调Trend: “↑”, metabolites in the group treated with flumethrin are up-regulated expression compared to the control group; “↓”, metabolites in the group treated with flumethrin are down-regulated expression compared to the control group

新窗口打开|下载CSV

将共同差异代谢物进行富集分析,最终确定到氨基糖与核糖代谢、药物代谢-其他酶类代谢、α-亚麻酸代谢、其他代谢途径4条代谢通路（表3）,其中尿苷二磷酸-N-乙酰葡糖胺、硫唑嘌呤、反玉米素表达上调,其对应的氨基糖与核糖代谢和药物代谢-其他酶类代谢等代谢通路加强;甲苯、愈伤酸、9-羟基壬酸表达下调,α-亚麻酸代谢通路减弱;氢过氧化亚油酸在50 mg·kg^-1组和0.5 mg·kg^-1组表达下调,在5 mg·kg^-1 组中表达上调,因此氢过氧化亚油酸在50 mg·kg^-1 组和0.5 mg·kg^-1组中代谢通路减弱,在5 mg·kg^-1 组中代谢通路加强。

Table 3
表3
表3差异代谢物聚集的代谢通路
Table 3Metabolic pathways of differential metabolite aggregation

化合物名称Compound name	Pathway ID	通路Pathway
尿苷二磷酸-N-乙酰葡糖胺 UDP-N-acetylglucosamine	Map00520	氨基糖与核糖代谢 Amino sugar and nucleotide sugar metabolism
硫唑嘌呤 Azathioprine	Map00983	药物代谢-其他酶类代谢 Drug metabolism-other enzymes
反玉米素 Trans-zeatin	Map01100	其他代谢途径 Metabolic pathways
甲苯 Toluene	Map01100	其他代谢途径 Metabolic pathways
愈伤酸 Traumatic acid	Map00592	ɑ-亚麻酸代谢 alpha-Linolenic acid metabolism
9-羟基壬酸 9-oxononanoic acid	Map00592	ɑ-亚麻酸代谢 alpha-Linolenic acid metabolism
氢过氧化亚油酸 13(s)-HPODE	Map01100	其他代谢途径 Metabolic pathways

新窗口打开|下载CSV

3 讨论

氟氯苯氰菊酯作为第二代拟除虫菊酯类杀虫、杀螨剂,已被广泛应用到养蜂生产中,用于杀死蜂群中的螨虫。然而,氟氯苯氰菊酯在杀螨虫的过程中,对非靶标生物——蜜蜂也产生一定的危害,甚至污染蜂产品以及蜜蜂的食物,间接地影响着蜜蜂的健康^[3,4,5,6]。本研究利用代谢组学的方法,探究了亚致死剂量的氟氯苯氰菊酯对西方蜜蜂工蜂幼虫代谢物的影响。研究发现LC-MS技术能够分析蜜蜂幼虫的代谢物,谱图重叠良好;主成分分析和偏最小二乘判别显示各组平行样聚在一起,与对照组分开,说明氟氯苯氰菊酯影响着蜜蜂幼虫的代谢。图2-C中对照组样品中有一个未落在95%的置信区间内,可能是由于样品衍生化过程中的试验误差。PLS-DA分析模型显示结果稳定可靠。经过对代谢物进行检测发现,随着氟氯苯氰菊酯浓度的增加,差异代谢物也明显增加,鉴定到的差异代谢物也逐渐增加,说明氟氯苯氰菊酯浓度越高,对蜜蜂幼虫代谢影响越大。根据氟氯苯氰菊酯处理组与对照组之间差异代谢物的分子质量和二级质谱信息,筛选到29种氟氯苯氰菊酯处理组与对照组间共同的差异代谢物,其中,尿苷二磷酸-N-乙酰葡糖胺、硫唑嘌呤、愈伤酸、9-羟基壬酸、氢过氧化亚油酸等代谢物变化较为明显,可能是由于氟氯苯氰菊酯对蜜蜂幼虫体内这些物质富集的代谢通路影响较大,表现为差异代谢物含量异常。

尿苷二磷酸-N-乙酰葡糖胺是蜜蜂脂多糖合成的重要前体。蜜蜂翅翼由双层体壁组成,每层体壁含有多层几丁质^[13]。而在昆虫中,几丁质的合成对于其生长发育至关重要。蜜蜂幼虫生长发育到一定时期就要脱去旧表皮的束缚,躯体才会显著长大^[14]。几丁质合成酶利用尿苷二磷酸-N-乙酰葡糖胺合成几丁质^[15],经氟氯苯氰菊酯处理,造成尿苷二磷酸-N-乙酰葡糖胺合成系统紊乱,几丁质合成异常,蜜蜂蜕皮羽化过程失调,因而蜜蜂出现残翅现象,这从另一方面解释了牛新月^[5]的研究结果中蜜蜂出现的残翅现象。

硫唑嘌呤与还原性谷胱甘肽反应最终转变为硫唑嘌呤活性代谢物6-硫鸟嘌呤核苷酸,其整合到细胞,影响DNA的复制及RNA的表达而发挥免疫抑制作用^[16]。硫唑嘌呤在50 mg·kg^-1组中的FC值高达22.9131,而0.5 mg·kg^-1组仅为12.1241,表明随着氟氯苯氰菊酯剂量增加,硫唑嘌呤相对含量显著增加,具有明显的剂量依赖性。HERNANDEZ等^[17]研究表明,谷胱甘肽S-转移酶（GST）是昆虫降解氟氯苯氰菊酯等杀螨剂的一种机制,亚致死剂量的氟氯苯氰菊酯会使得GST基因表达上调。除此之外,蜜蜂代谢氟氯苯氰菊酯引起硫唑嘌呤表达上调,硫唑嘌呤代谢需要大量的还原性谷胱甘肽参与反应,引起还原性谷胱甘肽表达上调,这从另一方面解释了蜜蜂的解毒机制。

愈伤酸、9-羟基壬酸为卵磷脂的代谢产物,有研究表明,卵磷脂有分解油脂、清除过氧化物的作用^[18]。在蜜蜂体内,卵磷脂代谢为愈伤酸和9-羟基壬酸受到多种酶的调节,其中一个酶是过氧化物裂解酶（hydroperoxide lyase,HPL）,HPL属于细胞色素P450（CYP450）类蛋白质家族中的一种^[19]。CYP450在外来化合物,如农药、环境污染物和致癌物质的新陈代谢中起着至关重要的作用^[20]。此外,CYP450在维生素代谢、不饱和脂肪酸的氧化和胆固醇的生物合成中起主要作用^[21]。CYP450酶可以通过受体依赖机制被各种外源物质和内源底物转录激活^[20]。药物在蜜蜂体内通过两个步骤：经CYP450酶的氧化还原反应及水解,再经GST的乙酰化、甲基化反应或毒物及其代谢物与内源性物质结合,加速毒物代谢或外排^[22]。CYP450酶系具有降低过氧化水平的作用,而9-羟基壬酸是脂肪酸过氧化的主要产物之一^[23]。9-羟基壬酸是脂肪酸代谢产生的主要醛类物质,口服9-羟基壬酸可显著提高脂质过氧化水平,同时显著降低大鼠肝脏中脂肪的生成^[24]。50 mg·kg^-1组9-羟基壬酸FC值仅为0.098,表明对照组9-羟基壬酸的相对含量是50 mg·kg^-1组的10倍。并且随着剂量的增加,FC值明显下降,表明9-羟基壬酸和愈伤酸具有明显的剂量依赖性。而氟氯苯氰菊酯处理会使9-羟基壬酸和愈伤酸含量下降,这可能是GST和CYP450酶共同作用的结果。

亚油酸是蜜蜂体内极为重要的不饱和脂肪酸之一,当蜜蜂缺乏亚油酸时,会出现蜕皮困难、发育不良、生长缓慢等症状^[25]。亚油酸在脂肪氧合酶和CYP450的作用下容易被氧化^[26],在脂肪氧合酶作用下,亚油酸主要形成氢过氧化十八碳二烯酸,而在CYP450作用下,亚油酸则经环氧化形成环氧化物^[27]。GST是蜜蜂体内极为重要的解毒酶之一,其可以催化谷胱甘肽和有毒化合物的轭合,增加其水溶性,降低产物毒性,同时还能缓解氟氯苯氰菊酯等杀螨剂导致的氧化应激并减少脂质过氧化物的产生^[28]。有研究表明,谷胱甘肽能抑制亚油酸的氧化^[29]。除此之外,谷胱甘肽过氧化物酶可以分解机体产生的异常过氧化物,可催化还原型谷胱甘肽转化为氧化型谷胱甘肽,将毒性过氧化物转化为羟基化合物,以保护细胞膜免受氧化应激损伤^[30]。而氢过氧化亚油酸是亚油酸过氧化的产物之一,氢过氧化亚油酸相对含量表明了氟氯苯氰菊酯对蜜蜂的氧化应激水平。氢过氧化亚油酸在0.5、50 mg·kg^-1组的FC值分别为0.4183、0.5253,但是5 mg·kg^-1组FC值高达4.7824,表明氢过氧化亚油酸在低、高剂量组相对含量远低于对照组,中剂量组相对含量远高于对照组,可能是谷胱甘肽在发挥作用,阻止了亚油酸氧化,从而造成3组间氢过氧化亚油酸含量存在差异。氢过氧化亚油酸在不同剂量氟氯苯氰菊酯处理组中的相对含量与氧化型谷胱甘肽相对含量的趋势相似,这一结果表明谷胱甘肽在缓解脂质氧化应激中的作用,进一步表明在蜜蜂体内氟氯苯氰菊酯引起的氧化应激可能主要是由还原性谷胱甘肽参与消除的。

4 结论

运用液相色谱-质谱联用技术对氟氯苯氰菊酯处理后的蜜蜂工蜂幼虫的差异代谢物进行分析,筛选到29种对照组与处理组间共同的差异代谢物,其中有7种差异代谢物显著聚集于氨基糖和核苷酸糖代谢、药物代谢-其他酶类、α-亚麻酸代谢等途径,对这些结果进一步分析发现氟氯苯氰菊酯对蜜蜂幼虫的健康存在一定的威胁,研究结果可为进一步阐释蜜蜂的解毒机制提供参考。

参考文献 View Option 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1] 刘朋飞, 吴杰, 李海燕, 林素文. 中国农业蜜蜂授粉的经济价值评估
中国农业科学, 2011,44(24):5117-5123.
[本文引用: 1]

LIU P F, WU J, LI H Y, LIN S W. Economic values of bee pollination to China’s agriculture
Scientia Agricultura Sinica, 2011,44(24):5117-5123. (in Chinese)
[本文引用: 1]

[2] 赵静, 李熠, 薛晓峰. 治螨新药——氟氯苯氰菊酯的性能和应用
中国蜂业, 2004,55(4):17-18.
[本文引用: 1]

ZHAO J, LI Y, XUE X F. New drug for mite treatment—— Performance and application of flumethrin
Apiculture of China, 2004,55(4):17-18. (in Chinese)
[本文引用: 1]

[3] QI S Z, ZHU L Z, WANG D H, WANG C, CHEN X F, XUE X F, WU L M. Flumethrin at honey-relevant levels induces physiological stresses to honey bee larvae (Apis mellifera L.) in vitro
Ecotoxicology and Environmental Safety, 2020,190:110101.
DOI:10.1016/j.ecoenv.2019.110101URL [本文引用: 4]

[4] QI S Z, NIU X Y, WANG D H, WANG C, ZHU L Z, XUE X F, ZHANG Z Y, WU L M. Flumethrin at sublethal concentrations induces stresses in adult honey bees (Apis mellifera L.)
The Science of the Total Environment, 2020,700:134500.
DOI:10.1016/j.scitotenv.2019.134500URL [本文引用: 4]

[5] 牛新月. 氟氯苯氰菊酯对意大利蜜蜂(Apis mellifera ligustica L.)的毒性作用研究
[D]. 新乡: 河南科技学院, 2019.
[本文引用: 5]

NIU X Y. Toxic effects of flumethrin on Apis mellifera ligustica L. (Hymenoptera: Apidae)
[D]. Xinxiang: Henan Institute of Science and Technology, 2019. (in Chinese)
[本文引用: 5]

[6] TAN K, YANG S, WANG Z W, MENZEL R. Effect of flumethrin on survival and olfactory learning in honeybees
PLoS ONE, 2013,8(6):e66295.
DOI:10.1371/journal.pone.0066295URL [本文引用: 4]

[7] 江武军, 何旭江, 王子龙, 颜伟玉, 曾志将, 吴小波. 中华蜜蜂细胞色素CYP9E2基因克隆及其表达分析
昆虫学报, 2016,59(10):1050-1057.
[本文引用: 1]

JIANG W J, HE X J, WANG Z L, YAN W Y, ZENG Z J, WU X B. Cloning and expression analysis of cytochrome CYP9E2 gene in the Chinese honeybee, Apis cerana cerana
Acta Entomologica Sinica, 2016,59(10):1050-1057. (in Chinese)
[本文引用: 1]

[8] YU L S, LIU F, WU H, TAN H R, RUAN X C, CHEN Y, CHAO Z. Flumethrin residue levels in honey from apiaries of China by high-performance liquid chromatography
Journal of Food Protection, 2015,78(1):151-156.
DOI:10.4315/0362-028X.JFP-14-297URL [本文引用: 1]

[9] 韩爱华. 氯氰菊酯和氟氯苯氰菊酯免疫分析化学研究
[D]. 扬州: 扬州大学, 2007.
[本文引用: 1]

HAN A H. Studies on immunochemistry for analysis of cypermethrin and flumethrin
[D]. Yangzhou: Yangzhou University, 2007. (in Chinese)
[本文引用: 1]

[10] ROSENKRANZ P, AUMERIER P, ZIEGELMANN B. Biology and control of Varroa destructor
Journal of Invertebrate Pathology, 2010,103:S96-S119.
DOI:10.1016/j.jip.2009.07.016URL [本文引用: 1]

[11] BOGDANOV S. Current status of analytical methods for the detection of residues in bee products
Apiacta, 2003,38:190-197.
[本文引用: 1]

[12] JOHNSON R M, ELLIS M D, MULLIN C A, FRAZIER M. Pesticides and honey bee toxicity—USA
Apidologie, 2010,41(3):312-331.
DOI:10.1051/apido/2010018URL [本文引用: 1]

[13] 李帅伟, 杨庆生, 苗博钧. 蜜蜂翅翼几丁质结构的粘弹性力学性能研究//中国力学大会论文集, 2019: 3343-3349.
[本文引用: 1]

LI S W, YANG Q S, MIAO B J. Study on viscoelastic properties of chitin structure of bee wings//The Chinese Congress of Theoretical and Applied Mechanics
2019: 3343-3349. (in Chinese)
[本文引用: 1]

[14] 曾志将. 蜜蜂生物学. 北京: 中国农业出版社, 2007: 2-3.
[本文引用: 1]

ZENG Z J. Bee Biology. Beijing: China Agriculture Press, 2007: 2-3. (in Chinese)
[本文引用: 1]

[15] GLASER L, BROWN D H. The synthesis of chitin in cell-free extracts of Neurospora crassa
Journal of Biological Chemistry, 1957,228(2):729-742.
DOI:10.1016/S0021-9258(18)70655-8URL [本文引用: 1]

[16] KARNER S, SHI S J, FISHER C, SCHAEFFELER E, NEURATH M F, HERRLINGER K R, HOFMANN U, SCHWAB M. Determination of 6-thioguanosine diphosphate and triphosphate and nucleoside diphosphate kinase activity in erythrocytes: novel targets for thiopurine therapy?
Therapeutic Drug Monitoring, 2010,32(2):119-128.
DOI:10.1097/FTD.0b013e3181d12f19URL [本文引用: 1]

[17] HERNANDEZ E P, KUSAKISAKO K, TALACTAC M R, GALAY R L, HATTA T, FUJISAKI K, TSUJI N, TANAKA T. Glutathione S-transferases play a role in the detoxification of flumethrin and chlorpyrifos in Haemaphysalis longicornis
Parasites and Vectors, 2018,11(1):460.
DOI:10.1186/s13071-018-3044-9URL [本文引用: 1]

[18] WOJCICKI J, PAWLIK A, SAMOCHOWIEC L, KALDONSKA M, MYSLIWIEC Z. Clinical evaluation of lecithin as a lipid-lowering agent
Phytotherapy Research, 1995,9(8):597-599.
DOI:10.1002/(ISSN)1099-1573URL [本文引用: 1]

[19] HUSSON F, BOMPAS D, KERMASHA S, BELIN J M. Biogeneration of 1-octen-3-ol by lipoxygenase and hydroperoxide lyase activities of Agaricus bisporus
Process Biochemistry, 2001,37(2):177-182.
DOI:10.1016/S0032-9592(01)00201-1URL [本文引用: 1]

[20] MANIKANDAN P, NAGINI S. Cytochrome P450 structure, function and clinical significance: A review
Current Drug Targets, 2018,19(1):38-54.
[本文引用: 2]

[21] OMURA T. Forty years of cytochrome P450
Biochemical and Biophysical Research Communications, 1999,266(3):690-698.
DOI:10.1006/bbrc.1999.1887URL [本文引用: 1]

[22] 郭晓波. 甘草解毒作用及机制的研究
[D]. 西安: 陕西师范大学, 2016.
[本文引用: 1]

GUO X B. Research on detoxification effect and mechanism of glycyrrhiza
[D]. Xi’an: Shaanxi Normal University, 2016. (in Chinese)
[本文引用: 1]

[23] REN R, HASHIMOTO T, MIZUNO M, TAKIGAWA H, YOSHIDA M, AZUMA T, KANAZAWA K. A lipid peroxidation product 9- oxononanoic acid induces phospholipase A₂ activity and thromboxane A₂ production in human blood
Journal of Clinical Biochemistry and Nutrition, 2013,52(3):228-233.
DOI:10.3164/jcbn.12-110URL [本文引用: 1]

[24] KANAZAWA K, NATAKE M. Identifications of 9-oxononanoic acid and hexanal in liver of rat orally administered with secondary autoxidation products of linoleic acid
Agricultural and Biological Chemistry, 1986,50(1):115-120.
[本文引用: 1]

[25] 王帅. 意大利蜜蜂饲粮中适宜赖氨酸添加水平的研究
[D]. 泰安: 山东农业大学, 2018.
[本文引用: 1]

WANG S. Studies of optimal lysine levels in diet of honeybee (Apis mellifera L.)
[D]. Taian: Shandong Agricultural University, 2018. (in Chinese)
[本文引用: 1]

[26] 姜春姣, 江芸, 耿志明, 张牧焓, 孙冲, 卞欢, 王道营, 徐为民. 亚油酸氧化产物——羟基十八碳二烯酸的研究进展
食品科学, 2018,39(7):278-284.
[本文引用: 1]

JIANG C J, JIANG Y, GENG Z M, ZHANG M H, SUN C, BIAN H, WANG D Y, XU W M. Progress in research on hydroxyoctadecaenoic acids as oxidation products of linoleic acid
Food Science, 2018,39(7):278-284. (in Chinese)
[本文引用: 1]

[27] NIKI E, YOSHIDA Y, SAITO Y, NOGUCHI N. Lipid peroxidation: Mechanisms, inhibition, and biological effects
Biochemical and Biophysical Research Communications, 2005,338(1):668-676.
DOI:10.1016/j.bbrc.2005.08.072URL [本文引用: 1]

[28] 徐莉, 王建华, 梅宇, 李冬植. 解毒酶和转运蛋白介导的害虫抗药性分子机制研究进展
农药学学报, 2020,22(1):1-10.
[本文引用: 1]

XU L, WANG J H, MEI Y, LI D Z. Research progress on the molecular mechanisms of insecticides resistance mediated by detoxification enzymes and transporters
Chinese Journal of Pesticide Science, 2020,22(1):1-10. (in Chinese)
[本文引用: 1]

[29] 张娟. 谷胱甘肽对乳酸菌胁迫抗性的调控机制研究
[D]. 无锡: 江南大学, 2008.
[本文引用: 1]

ZHANG J. Regulation mechanism of glutathione on stress resistance of lactic acid bacteria
[D]. Wuxi: Jiangnan University, 2008. (in Chinese)
[本文引用: 1]

[30] 李春燕, 孙传政, 宋鑫. 肿瘤细胞死亡的一种新形式——铁死亡
中国生物化学与分子生物学报, 2019,35(11):1208-1214.
[本文引用: 1]

LI C Y, SUN C Z, SONG X. A new form of tumor cell death: Ferroptosis
Chinese Journal of Biochemistry and Molecular Biology, 2019,35(11):1208-1214. (in Chinese)
[本文引用: 1]