删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

新疆马铃薯甲虫对噻虫嗪的抗性监测及其细胞色素P450基因表达分析

本站小编 Free考研考试/2021-12-26

石鑫,1, 李莎1, 王志敏1, 付开赟2, 付文君3, 姜卫华,11南京农业大学植物保护学院,南京 210095
2新疆农业科学院植物保护研究所/农业农村部西北荒漠绿洲作物有害生物综合治理重点实验室,乌鲁木齐 830091
3新疆伊犁州农业技术推广总站,新疆伊宁 835000

Resistance Monitoring to Thiamethoxam and Expression Analysis of Cytochrome P450 Genes in Leptinotarsa decemlineata from Xinjiang

SHI Xin,1, LI Sha1, WANG ZhiMin1, FU KaiYun2, FU WenJun3, JIANG WeiHua,11College of Plant Protection, Nanjing Agricultural University, Nanjing 210095
2Research Institute of Plant Protection, Xinjiang Academy of Agricultural Sciences/Key Laboratory of Intergraded Management of Harmful Crop Vermin of China North-Western Oasis, Ministry of Agriculture and Rural Affairs, Urumqi 830091
3Agricultural Technology Extension Master Station of Yili Prefecture, Yining 835000, Xinjiang

通讯作者: 姜卫华,E-mail: jwh@njau.edu.cn

责任编辑: 岳梅
收稿日期:2020-10-30接受日期:2020-11-21
基金资助:国家重点研发计划.2018YFD0200802


Received:2020-10-30Accepted:2020-11-21
作者简介 About authors
石鑫,E-mail: 907094441@qq.com








摘要
【目的】以世界性害虫马铃薯甲虫(Leptinotarsa decemlineata)为研究对象,筛选其参与噻虫嗪解毒的细胞色素P450主效基因。【方法】于2018、2019年利用点滴法监测新疆察布查尔县、伊宁县、塔城市、乌鲁木齐市和吉木萨尔县11个马铃薯甲虫田间种群对噻虫嗪的抗性水平,并测定分析噻虫嗪LD50处理72 h后成虫中3种主要解毒酶细胞色素P450酶(P450)、谷胱甘肽S-转移酶(GST)和酯酶(EST)的活性变化。利用Illumina HiSeqTM 2500 高通量测序技术进行转录组测序,获得抗性和敏感种群之间的差异表达基因(DEG),运用实时荧光定量PCR(qRT-PCR)技术验证3个P450基因并分析6个P450基因(CYP4C1、CYP9e2CYP305a1CYP9Z25、CYP9YaCYP9Yc)在不同种群及噻虫嗪处理后的表达情况。【结果】抗性监测发现,2018年的YN1、URMQ1、TC和QPQL1以及2019年的YN2和JMS种群成虫对噻虫嗪抗性倍数(RR)分别达到5.99、8.81、10.86和20.33倍和11.82、14.05倍,为低至中水平抗性。解毒酶活性测定结果表明,2个敏感种群URMQ2、URMQ3以及2个抗性种群YN2、JMS成虫分别经噻虫嗪LD50处理72 h后其P450活性均显著增加,分别为对照的1.76、2.75、1.91和1.66倍。另外,URMQ3种群的GST(CDNB和DCNB为底物)及URMQ2、YN2种群的EST活性显著升高,分别为对照的1.19、2.10倍和1.35、1.91倍。转录组分析表明,通过合并组装分别获得噻虫嗪敏感和抗性种群样品56 872 051和62 249 136个原始序列数据以及55 903 706和61 082 076个过滤后的序列数据,过滤后的序列长度分别为8.39和9.16 G,碱基错误率均为0.03%。抗性种群中差异表达的P450基因13个,其中2个基因显著上调。3个上调的P450基因CYP4C1、CYP9e2CYP305a1的qRT-PCR检测结果与转录组测序结果基本一致,qRT-PCR分析还发现CYP9Ya、CYP9Yc、CYP4C1、CYP305a1CYP9Z25在抗性种群成虫中表达量显著增加,噻虫嗪LD50处理均可引起URMQ2种群4龄幼虫和成虫CYP9Ya、CYP9Yc表达量显著上调。相关性分析发现成虫对噻虫嗪的抗性水平与CYP9Ya表达量呈显著正相关。【结论】CYP9Ya可能在马铃薯甲虫中对噻虫嗪具有重要的解毒作用,其他基因的作用也不能排除。
关键词: 马铃薯甲虫;噻虫嗪;抗药性;细胞色素P450;基因表达

Abstract
【Objective】Colorado potato beetle (CPB), Leptinotarsa decemlineata, is a significant agricultural pest in potato worldwide. The objective of the study is to identify the main cytochrome P450 genes involved in metabolic regulation of neonicotinoid insecticide thiamethoxam in L. decemlineata.【Method】In this study, the resistance level to thiamethoxam was assayed using topical applications in 11 field populations of L. decemlineata from Qapqal (QPQL), Yining (YN), Tacheng (TC), Urumqi (URMQ), Jimusar (JMS) of Xinjiang in 2018 and 2019. The activities of three detoxifying enzymes including cytochrome P450 monooxygenase (P450), glutathioneS-transferase (GST) and esterase (EST), exposed to LD50of thiamethoxam for 72 h were analyzed by in vitro enzyme activity assay. The differentially expressed genes (DEGs) of L. decemlineata adults susceptible and resistant to thiamethoxam were detected by Illumina HiSeqTM 2500 sequencing platform. The expression verification of three P450 genes and the expression change of six P450 genes in different populations and the beetles exposed to LD50 of thiamethoxam for 72 h were performed by real-time quantitative PCR (qRT-PCR).【Result】Resistance monitoring results showed that the adults of YN1, URMQ1, TC and QPQL1 populations in 2018 and YN2 and JMS in 2019 developed low to moderate level resistance to thiamethoxam with resistance ratio (RR) ranging from 5.99 to 20.33 times. The P450 activities of adults from two sensitive populations URMQ2 and URMQ3 and two resistant populations YN2 and JMS after treatment with thiamethoxam LD50for 72 h were significantly increased to 1.76, 2.75, 1.91 and 1.66 times as high as the control, respectively. In addition, the GST (CDNB and DCNB as substrates) of URMQ3 population and the EST activities of URMQ2 and YN2 populations were obviously increased to 1.19, 2.10 and 1.35, 1.91 times of the control, respectively. Based on Illumina RNA sequencing, the raw reads (56 872 051 and 62 249 136), clean reads (55 903 706 and 61 082 076), clean bases (8.39 and 9.16 G) and error of basic group (0.03%) were obtained for the thiamethoxam-susceptible and -resistant samples of adult populations, respectively. Thirteen differentially expressed P450 genes were identified in the thiamethoxam-resistant sample, of which two genes were significantly up-regulated. The expression levels of three P450 genes, namely CYP4C1,CYP9e2and CYP305a1by RNA sequencing were validated by qRT-PCR analysis. Furthermore, CYP9Ya, CYP9Yc, CYP4C1, CYP305a1 and CYP9Z25 were up-regulated in the thiamethoxam-resistant adults, and the expression of CYP9Ya, CYP9Yc increased significantly in the 4th instar larvae and adults under thiamethoxam treatment. It was found by correlation analysis that there was a significant positive correlation between resistance level to thiamethoxam and CYP9Ya expression of adults. 【Conclusion】CYP9Ya may be involved in the production of thiamethoxam resistance in L. decemlineata, and the role of other genes cannot be ruled out.
Keywords:Leptinotarsa decemlineata;thiamethoxam;resistance;cytochrome P450;gene expression


PDF (615KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
石鑫, 李莎, 王志敏, 付开赟, 付文君, 姜卫华. 新疆马铃薯甲虫对噻虫嗪的抗性监测及其细胞色素P450基因表达分析[J]. 中国农业科学, 2021, 54(14): 3004-3016 doi:10.3864/j.issn.0578-1752.2021.14.007
SHI Xin, LI Sha, WANG ZhiMin, FU KaiYun, FU WenJun, JIANG WeiHua. Resistance Monitoring to Thiamethoxam and Expression Analysis of Cytochrome P450 Genes in Leptinotarsa decemlineata from Xinjiang [J]. Scientia Acricultura Sinica, 2021, 54(14): 3004-3016 doi:10.3864/j.issn.0578-1752.2021.14.007


开放科学(资源服务)标识码(OSID):

0 引言

【研究意义】马铃薯甲虫(Leptinotarsa decemlineata)属鞘翅目叶甲科,是世界公认的马铃薯毁灭性检疫害虫。除直接取食危害寄主植物叶片外,成虫还能传播褐斑病、环腐病等病害,造成农作物产量和品质的巨大损失[1,2]。马铃薯甲虫自20世纪90年代初经哈萨克斯坦口岸传入我国新疆,现已扩散至天山以北地带,对当地马铃薯安全生产构成持续的威胁[3]。作为第4大类杀虫剂[4],新烟碱类杀虫剂以其作用机制独特、高效、低毒、对环境安全成为防治马铃薯甲虫的主要杀虫剂之一,其中1995年问世的吡虫啉和1998年市场化的噻虫嗪分别为第1代和第2代新烟碱类药剂的代表品种[5]。随着使用的频繁,国内外抗性调查发现马铃薯甲虫已对吡虫啉等产生了不同水平的抗性。因此,阐明相关抗药性分子机制对于防止抗性加剧及延长该类药剂的使用寿命具有重要意义。【前人研究进展】ZHAO等[6]2000年报道纽约州长岛马铃薯甲虫成虫对吡虫啉产生了100—150倍的极高水平抗性;MOTA-SANCHEZ等[7]于2006年报道该地区种群对吡虫啉的抗性已发展到309倍,另外美国特拉华、爱达荷、伊利诺斯等地区以及加拿大的安大略、爱德华王子岛、魁北克田间种群也对吡虫啉和噻虫嗪产生不同程度的抗性。在我国,王志田等[8]调查发现新疆奇台和乌鲁木齐马铃薯甲虫成虫对吡虫啉的抗性个体达到30%左右;刘萍等[9]监测发现2009年新疆6个马铃薯甲虫成虫种群中有3个对啶虫脒和噻虫嗪产生低水平抗性,2010年监测的6个种群均对噻虫嗪产生了抗性,其中低抗和中抗种群各3个。由此可见,新疆马铃薯甲虫对新烟碱类药剂尤其是噻虫嗪的抗性逐年发展。解毒代谢作用的增强是昆虫抗药性产生的重要机制之一。昆虫体内的解毒酶系主要包括细胞色素P450单加氧酶系(cytochrome P450 monooxygenases,P450)、谷胱甘肽S-转移酶(glutathione S-transferase,GST)和酯酶(esterase,EST),其中由P450介导的昆虫抗药性是主要的抗性机制。研究表明害虫对新烟碱类杀虫剂抗性的产生主要与P450基因过表达有关。KAPLANOGLU等[10]对吡虫啉抗性马铃薯甲虫种群中显著上调的CYP4Q3 进行RNA干扰(RNA interference,RNAi),验证了该基因在抗性产生中的作用;CLEMENTS等[11,12]发现美国威斯康星州吡虫啉抗性马铃薯甲虫成虫CYP9Z26过表达,通过RNAi可使马铃薯甲虫对吡虫啉的抗性水平降低,作者随后在该抗性种群中又发现另一个过表达的P450基因CYP6K1[13]。【本研究切入点】近年来抗性监测已发现新疆马铃薯甲虫对噻虫嗪的敏感性逐年下降,抗性分布不断扩展且水平较高,但有关抗性分子机制尚不明确。【拟解决的关键问题】监测新疆不同地区马铃薯甲虫对噻虫嗪的抗性水平,测定噻虫嗪处理对主要解毒酶活性的影响,利用转录组分析抗性和敏感种群的差异表达基因(differentially expressed gene,DEG),定量分析不同田间种群以及噻虫嗪处理后的6个P450基因的表达水平,筛查与噻虫嗪抗性相关的P450基因,为阐明马铃薯甲虫对噻虫嗪的代谢抗性分子机制打下基础,同时为马铃薯甲虫的抗性治理提供依据。

1 材料与方法

试验于2018—2019年在新疆农业科学院植物保护研究所安宁渠试验基地和南京农业大学植物保护学院昆虫生理生化与分子生物学实验室完成。

1.1 供试昆虫

马铃薯甲虫分别于2018年和2019年6—7月采自新疆察布查尔县(Qapqal,QPQL)、伊宁县(Yining,YN)、塔城市(Tacheng,TC)、乌鲁木齐市(Urumqi,URMQ)和吉木萨尔县(Jimusar,JMS)的田间马铃薯苗上,种群具体采集信息见表1。室内用马铃薯苗饲养,养虫室条件为(26±1)℃,相对湿度50%—60%,光周期为光/暗=16 h/8 h。选取大小一致、生长状态良好的4龄幼虫和成虫作为试虫。2018年抗性测定参照的敏感种群QPQL-S于2015年6月采自察布查尔县扎库齐牛录乡,不接触药剂下饲养3年。

Table 1
表1
表1马铃薯甲虫种群采集信息
Table 1Sampling information of L. decemlineatapopulations
采集时间 Sampling time种群 Population采集地点 Sampling location
2018-06察布查尔QPQL1察布查尔县扎库齐牛录乡Zhakuqiniulu Town, Qapqal County
2018-06伊宁YN1伊宁县英塔木乡Yingtamu Town, Yining County
2018-06塔城TC塔城市阿西尔乡Axier Town, Tacheng City
2018-07乌鲁木齐URMQ1乌鲁木齐市安宁渠镇Anningqu Town, Urumqi City
2019-06察布查尔QPQL2察布查尔县扎库齐牛录乡Zhakuqiniulu Town, Qapqal County
2019-06伊宁YN2伊宁县英塔木乡Yingtamu Town, Yining County
2019-06伊宁YN3伊宁县胡地于孜乡Hudiyuzi Town, Yining County
2019-07乌鲁木齐URMQ2乌鲁木齐市萨尔达坂乡Saerdaban Town, Urumqi City
2019-07乌鲁木齐URMQ3乌鲁木齐市头屯河区Toutunhe District, Urumqi City
2019-07乌鲁木齐URMQ4乌鲁木齐市板房沟乡Banfanggou Town, Urumqi City
2019-07吉木萨尔JMS吉木萨尔县泉子街镇Quanzijie Town, Jimusar County

新窗口打开|下载CSV

1.2 供试药剂

97%吡虫啉原药(江苏红太阳有限公司);95%噻虫嗪原药(盐城双宁农化有限公司)。

1.3 毒力测定

采用微量点滴法。使用丙酮溶解原药并稀释制备1、0.5、0.25、0.125、0.0625、0.03125、0.015625、0.007813 mg·mL-1系列梯度浓度的药液。采用Hamilton微量进样器进行药剂点滴,药剂点滴部位:成虫腹部腹面第二节,4龄幼虫腹部背面的第3、4节;点滴量:4龄幼虫0.22 μL/头,成虫1.1 μL/头。点滴同体积的丙酮作为对照。每个处理试虫10头,共3个重复。处理后置于塑料培养皿(直径9 cm),在室内用新鲜马铃薯叶饲养,72 h后检查死亡虫数。马铃薯甲虫死亡标准参考刘萍等[9]和卢伟平等[14]

1.4 解毒酶活性测定

1.4.1 酶液制备 收集URMQ2、URMQ3、YN2和JMS种群成虫分别经噻虫嗪LD50点滴法处理72 h后的存活个体,对照为丙酮处理,将试虫单头放于匀浆器内,加入1.2 mL 0.1 mol·L-1磷酸缓冲液(pH 7.6,其中含1 mmol·L-1 EDTA,1 mmol·L-1 DTT,1 mmol·L-1 PTU,1 mmol·L-1 PMSF)匀浆用于细胞色素P450单加氧酶活性测定;加入1.2 mL 0.1 mol·L-1磷酸缓冲液(pH 7.6)匀浆用于谷胱甘肽S-转移酶活性测定;加入1.2 mL 0.02 mol·L-1磷酸缓冲液(pH 7.0)匀浆用于酯酶活性测定。匀浆液于4℃,13 000×g下离心15 min后,将上清液转移到新的1.5 mL离心管,重复离心一次后,取上清液作为酶源。每个处理6个生物学重复,3个技术性重复。

1.4.2 P450活性测定 在96孔酶标板中加入100 μL 2 mmol·L-1对硝基苯甲酚(PNA)和90 μL酶液,混合液在30℃温育3 min,然后加入10 μL 9.6 mmol·L-1 NADPH启动反应。利用SPECTRA max®340-PC型酶标仪(美国Molecular Devices公司)在波长405 nm下记录15 min内的吸光值变化,酶促反应温度为30℃。

1.4.3 GST活性测定 在96孔酶标板中分别加入10 μL稀释10倍的酶液(CDNB为底物)或25 μL酶液(DCNB为底物)、1.2 mmol·L-1CDNB或DCNB 100 μL,以及6 mmol·L-1谷胱甘肽(GSH)100 μL。在波长340 nm下记录10 min内的吸光值变化,酶促反应温度为27℃。

1.4.4 EST活性测定 在96孔酶标板中分别加入20 μL稀释10倍的酶液,0.2 mol·L-1 PBS(pH 6.0)与固蓝RR盐和α-乙酸萘酯(α-NA)的乙醇溶液(100 mmol·L-1)以1﹕2﹕0.01比例的混合液205 μL,在波长450 nm下记录10 min内的吸光值变化,酶促反应温度为27℃。

1.4.5 蛋白含量测定 采用Bradford考马斯亮蓝法测定粗酶液的蛋白质含量,以BSA做标准曲线[15]

1.5 转录组测序

URMQ3和JMS种群成虫用于转录组测序分析,每个种群设置3个重复,每个重复3头试虫。

转录组测序数据组装及分析:RNA的提取、cDNA文库的建立及测序工作由北京诺禾致源科技股份有限公司完成,采用Illumina HiSeqTM2500系统进行测序。使用HISAT2v2.0.5构建参考基因组的索引,并使用HISAT2 v2.0.5将配对末端clean reads与参照基因组比对,采用StringTie软件新转录本进行组装[16]

转录组差异表达基因的筛选:首先对原始的readcount进行标准化(normalization),主要是对测序深度的校正。然后通过统计学模型进行假设检验概率(P-value)的计算,最后进行多重假设检验校正(BH),得到FDR值(错误发现率)。以P<0.05、|log2fold change(FC)|≥1作为筛选标准[17],|log2FC|越大,则基因表达差异越明显。筛选出解毒代谢酶、表皮蛋白、ABC转运蛋白、烟碱型乙酰胆碱受体等与抗性相关的差异表达基因。

1.6 马铃薯甲虫种群P450基因表达量分析

1.6.1 总RNA提取及反转录 将马铃薯甲虫4龄幼虫和成虫单头分别按照Trizol试剂说明进行总RNA提取,于-80℃冰箱贮存备用,每个处理6个重复。按照PrimeScriptTM RT Reagent Kit with gDNA Eraser(日本TaKaRa公司)反转录试剂盒说明书合成cDNA。

1.6.2 实时荧光定量PCR 根据GenBank上马铃薯甲虫CYP9Ya、CYP9Yc、CYP9Z25、CYP4C1、CYP9e2、CYP305a1的序列(NCBI登录号分别为XP_023027514、XP_023018115、XP_023018114、XP_023021616、XP_ 023023178、XP_023015477),利用Beacon Designer 7.0软件设计荧光定量PCR引物,委托南京金斯瑞公司合成,以RPL4Ef1α作为内参基因[18],引物序列见表2。采用ABI 7500实时PCR系统(美国Applied Biosystems公司),以SYBR Green为染料,用稀释5倍的cDNA作为模板,测定P450基因的相对表达水平。反应体系:SYBR Premix Ex TaqTM(2×)10 µL,cDNA模板1.0 µL,上、下游引物各0.8 µL,ROX Reference Dye2(50×)0.4 µL,用无核酶水补足至20 µL。PCR程序采用两步法:95℃ 30 s;95℃ 5 s;60℃ 34 s,共40个循环,然后95℃ 15 s;60℃ 60 s;95℃ 15 s用于记录熔解曲线。每组试验重复3次。

Table 2
表2
表2马铃薯甲虫实时荧光定量PCR引物
Table 2The primers ofL. decemlineata for qRT-PCR
基因 Gene上游引物序列 Forward primer sequence (5′-3′)下游引物序列 Reverse primer sequence (5′-3′)
CYP9YaTGGACCAAAGAGTACCCGAAGGAATCGACAATTCAGTCCGAGGTACGA
CYP9YcATTTGGAAGTGGACCACGCAATTGGCCAAATTCCTTCTCGTCCAGCAAG
CYP9Z25ACTTGCACACCTTTGAGAGCATCCTTGTCAACTTTCTCATGCCCACCC
CYP4C1AGCCTGACATTCCACTTCCTGCTCCTCACTAACATCT
CYP9e2CATCATAGCACAAGCATTCTTTTCCTCCACATATTTCATCAG
CYP305a1TTCCGAAGAATACGATGGTGCGAAGATGAAGAGTTACAA
RPL4AAAGAAACGAGCATTGCCCTTCCGTTGTCGCTGACACTGTAGGGTTGA
Ef1αAAGGTTCCTTCAAGTATGCGTGGGGCACAATCAGCTTGCGATGTACCA

新窗口打开|下载CSV

1.7 数据处理与分析

死亡率和校正死亡率的计算使用Abbott公式,采用机率值分析法计算毒力回归曲线、LD50、相关系数及95%置信限,抗性倍数(RR)=田间测试种群的LD50/(相对)敏感种群的LD50。抗药性水平参照以下标准:敏感(RR<3.0)、敏感性降低(RR 3.1—5.0)、低水平抗性(RR 5.1—10.0)、中水平抗性(RR 10.1—40.0)、高水平抗性(RR 40.1—160.0)和极高水平抗性(RR>160.0)[19]。利用2-ΔΔCt方法[20]计算相对表达量,数据用Excel 2003进行分析。利用SPSS 18.0软件(One-way ANOVA)进行数据的差异显著性分析(Tukey检验,P<0.05),以及对抗性水平和基因表达量进行相关性分析。

2 结果

2.1 新疆不同马铃薯甲虫种群对噻虫嗪的抗性

对于4龄幼虫,2018年测试种群中,与敏感种群QPQL-S相比,YN1、TC和URMQ1种群对噻虫嗪均表现敏感(RR分别为1.88、2.14和2.42倍),而QPQL1种群为敏感性降低(RR为4.23倍);2019年测试种群中,URMQ2种群对噻虫嗪的LD50最低(5.52 ng/头),以此作为相对敏感种群,URMQ3、URMQ4、JMS、QPQL2、YN3种群对噻虫嗪均保持敏感(RR<3倍),仅YN2种群敏感性降低(RR为3.71倍)(表3)。

Table 3
表3
表3新疆不同马铃薯甲虫种群4龄幼虫和成虫对噻虫嗪的敏感性
Table 3Susceptibility to thiamethoxam of L. decemlineata 4th instar larvae and adults from Xinjiang
虫态
Stage
年份
Year
种群
Population
斜率Slope
(b±SE)
致死中量(95%置信限)
LD50(95% FL) (ng/larva)
抗性倍数
Resistant ratio (RR)
4龄幼虫
4th instar larva
2018QPQL-S3.456±0.4196.64 (5.12-8.61)1.00
YN12.588±0.09012.40 (10.01-15.43)1.88
TC3.236±0.21114.08 (11.53-17.20)2.14
URMQ13.438±0.25516.04 (13.08-19.67)2.42
QPQL11.255±0.08627.94 (17.02-45.81)4.23
2019URMQ210.911±2.5385.52 (5.21-5.84)1.00
URMQ32.829±0.2705.85 (4.69-7.29)1.05
JMS11.508±3.5715.93 (5.60-6.28)1.07
QPQL22.514±0.2409.64 (7.49-12.42)1.75
URMQ42.483±0.20512.97 (9.62-17.49)2.36
YN32.095±0.29314.99 (11.38-19.76)2.73
YN22.607±0.10520.41 (14.76-28.22)3.71
成虫
Adult
2018QPQL-S1.648±0.0989.42 (6.73-13.39)1.00
YN12.665±0.07556.33 (45.14-70.22)5.99
URMQ15.934±0.32182.76 (74.50-91.94)8.81
TC6.857±0.460102.11 (94.00-110.92)10.86
QPQL11.719±0.085191.11 (139.32-262.28)20.33
2019URMQ211.203±1.60724.69 (22.81-26.72)1.00
URMQ33.707±0.55137.65 (30.57-46.38)1.53
URMQ43.023±0.40658.31 (46.96-72.41)2.36
QPQL21.636±0.09185.25 (61.06-119.00)3.45
YN32.802±0.102111.19 (89.43-138.25)4.50
YN21.809±0.212291.97 (161.07-529.24)11.82
JMS1.384±0.138347.01 (204.51-588.80)14.05

新窗口打开|下载CSV

在成虫中,与QPQL-S相比,2018年测试种群均对噻虫嗪产生抗性,其中YN1和URMQ1种群为低水平抗性,而TC和QPQL1种群已达中等水平抗性,RR分别为5.99、8.81、10.86和20.33倍;2019年所有种群中,URMQ2种群对噻虫嗪的LD50最低,为24.69 ng/头,与之相比,URMQ3和URMQ4种群对噻虫嗪保持敏感(RR<3倍),QPQL2、YN3种群为敏感性降低,RR分别为3.45和4.50倍,YN2和JMS种群则达到中水平抗性,RR分别为11.82和14.05倍(表3)。

从两年的抗性监测结果可以看出新疆马铃薯甲虫成虫对噻虫嗪的敏感性普遍较低,抗性达低至中等水平,而4龄幼虫对噻虫嗪普遍敏感。

2.2 噻虫嗪胁迫下马铃薯甲虫的解毒酶活性

噻虫嗪可引起3种解毒酶的活性发生改变。其中P450活性在4个测试种群URMQ2、URMQ3、YN2、JMS处理组的活性均显著提高,分别为对照的1.76、2.75、1.91和1.66倍;URMQ3种群处理组的GST活性(CDNB和DCNB为底物)活性显著升高,分别为对照的1.19和2.10倍,而其他种群的GST活性无显著变化;URMQ2和YN2种群处理组的EST活性显著增加,分别为对照的1.35和1.91倍,其他种群的EST活性与对照相比差异不显著(图1)。

图1

新窗口打开|下载原图ZIP|生成PPT
图1噻虫嗪LD50处理72 h马铃薯甲虫成虫解毒酶活性的变化

图中数据为平均数±标准误,柱上不同字母表示差异显著(P<0. 05)。下同
Fig. 1Changes of detoxification enzyme activities of L. decemlineata adults exposed to LD50 of thiamethoxam for 72 h

Data in the figure are mean±SE. Different letters on the bars indicate significant difference atP<0. 05 level. The same as below


2.3 马铃薯甲虫敏感和抗性种群的转录组分析及验证

2.3.1 转录组数据组装 对噻虫嗪敏感种群UMRQ3和抗性种群JMS成虫进行转录组分析,分别获得56 872 051和62 249 136个原始序列数据以及55 903 706和61 082 076个过滤后的序列数据;过滤后的序列长度分别为8.39和9.16 G;碱基错误率均为0.03%;Phred数值>20的碱基占总碱基的百分比(Q20)分别为97.77%和97.75%;Phred数值>30的碱基占总碱基的百分比(Q30)分别为93.28%和93.25%;鸟嘌呤+胞嘧啶(guanine cytosine,GC)含量分别为38.95%和38.68%,表明转录组测序质量较高(表4)。

Table 4
表4
表4敏感和抗性马铃薯甲虫种群成虫的转录组测序结果
Table 4Sequencing results of susceptible and resistant adults of L. decemlineata from Xinjiang
样品
Sample
原始序列个数
Number of raw sequence
过滤后的序列个数
Number of filtered sequence
过滤后的序列长度
Clean bases
碱基错误率
Error of basic group (%)
Q20 (%)Q30 (%)GC含量GC content (%)
UMRQ356872051559037068.39 G0.0397.7793.2838.95
JMS62249136610820769.16 G0.0397.7593.2538.68

新窗口打开|下载CSV

2.3.2 抗药性相关的差异表达基因 按照解毒机制进行分类,鉴定到差异表达的P450基因13个,谷胱甘肽S-转移酶基因1个,酯酶基因7个,尿苷二磷酸-葡萄糖基转移酶基因9个,表皮蛋白基因6个,ABC转运蛋白基因12个,烟碱型乙酰胆碱受体基因11个(表5)。

Table 5
表5
表5抗性种群和相对敏感种群的差异表达基因
Table 5DEGs between the resistant population and relatively sensitive population
差异表达基因分类
DEG classification
所有基因
Whole transcriptome
上调基因
Up-regulated gene
下调基因
Down-regulated gene
显著上调基因Significantly up-regulated gene显著下调基因Significantly down-regulated gene
细胞色素P450 Cytochrome P4501331023
谷胱甘肽S-转移酶Glutathione S-transferase 10100
酯酶Esterase73400
尿苷二磷酸-葡萄糖基转移酶UDP-glucuronosyltransferase94513
表皮蛋白Cuticular protein63300
ABC转运蛋白ABC transporter124830
烟碱型乙酰胆碱受体Nicotinic acetylcholine receptor111300
log2FC>1为判断基因上调的标准,P<0.05即为显著。下同
Up-regulated transcripts were classified using a log2FC>1 and significant difference atP<0.05 level. The same as below

新窗口打开|下载CSV

抗性种群JMS中表达量上调的基因见表6,其中,显著上调的P450基因有2个(111509985、111511395)、ABC转运蛋白基因3个(有2个基因为新预测的基因,NCBI未见相关序列,由转录组测序公司编号)、尿苷二磷酸-葡萄糖基转移酶基因1个,显著上调的基因用加粗字体表示。

Table 6
表6
表6抗性种群中表达量上调的基因
Table 6Up-regulated genes in the resistant population
基因名称Gene symbol基因号 Gene IDlog2FC PPvalue
CYP4C1类基因 P450 4C1-like 1115099853.759.82E-17
CYP9e2类基因 P450 9e2-like 1115113952.150.001
可能为CYP305a1类基因 Probable P450 305a1 1115049771.850.449
酯酶E4类基因Esterase E4-like1115137092.410.051
酯酶B1类基因%2C转录本变种X1
Esterase B1-like%2C transcript variant X1
1115020421.180.190
羧酸酯酶6类基因 Carboxylesterase-6-like1115145301.070.553
2-羟乙酰鞘氨1-β类半乳糖糖基转移酶基因
2-hydroxyacylsphingosine 1-beta-galactosyltransferase-like
1115166321.690.444
2-羟乙酰鞘氨1-β类半乳糖糖基转移酶基因
2-hydroxyacylsphingosine 1-beta-galactosyltransferase-like
1115123271.900.457
未鉴定的LOC111503457 Uncharacterized LOC1115034571115034571.520.543
尿苷二磷酸-葡萄糖基转移酶2B30类基因
UDP-glucuronosyltransferase 2B30-like
1115133162.520.023
蛹表皮蛋白36类基因 Pupal cuticle protein 36-like1115018893.280.073
幼虫表皮蛋白A2B类基因 Larval cuticle protein A2B-like1115031501.160.173
内表皮结构糖蛋白SgAbd-8类基因
Endocuticle structural glycoprotein SgAbd-8-like
1115174843.910.173
ABC转运蛋白跨膜区 ABC transporter transmembrane regionnovel.11207.536.85E-06
多药耐药相关蛋白1类基因
Multidrug resistance-associated protein 1-like
1115164597.987.63E-06
ABC转运蛋白 ABC transporternovel.12737.350.003
ATP-结合盒G亚家族类基因1%2C 转录本变种X1
ATP-binding cassette sub-family G member 1-like%2C transcript variant X1
1115109791.790.553
烟碱型乙酰胆碱受体α3亚基 α3 1115126511.000.683

新窗口打开|下载CSV

2.3.3 抗性和敏感种群P450差异表达基因的验证 DEG测序文库结果显示,抗性种群JMS的CYP4C1CYP9e2基因表达量显著上调(P<0.05),分别为相对敏感种群的13.45和4.44倍,CYP305a1在抗性种群中的表达量为相对敏感种群的3.60倍,但差异不显著。qRT-PCR检测结果显示这3个P450基因CYP4C1、CYP9e2CYP305a1在JMS种群中的表达量显著上调(P<0.05),分别为相对敏感种群URMQ3的5.03、5.32和8.01倍(图2)。qRT-PCR检测结果与DEG测序文库显示的结果基本一致。

图2

新窗口打开|下载原图ZIP|生成PPT
图2qRT-PCR技术验证3个P450基因的表达

Fig. 2The relative level of three up-regulated P450 genes by qRT-PCR validation



2.4 新疆不同马铃薯甲虫种群P450基因的表达水平

2018年4龄幼虫中,与QPQL-S相比,CYP9YaCYP9Yc在QPQL1和YN1种群、CYP9Z25在TC种群的表达量均显著增加,表达量分别为QPQL-S的1.90、1.58、1.53、1.44和1.61倍(图3-A);2019年4龄幼虫中,与相对敏感种群URMQ2相比,YN3种群中CYP4C1、CYP9Z25的表达量最高,显著增加为URMQ2的6.74和4.15倍;JMS种群CYP9e2、CYP305a1以及QPQL2的CYP9Ya、CYP9Yc的表达量在所有种群中最高,显著上调,表达量分别是URMQ2的5.35、21.95、3.51和3.11倍。此外,其他种群的其他基因表达量无显著变化(图3-B)。

图3

新窗口打开|下载原图ZIP|生成PPT
图3马铃薯甲虫田间种群4龄幼虫P450基因的相对表达量(A:2018;B:2019)

Fig. 3Relative expression of P450 genes in the field populations of L. decemlineata 4th instar larvae from Xinjiang in 2018 (A) and 2019 (B)



2018年成虫中,CYP9YaCYP9YcCYP9Z25分别在QPQL1、URMQ1种群中显著上调,表达量分别为QPQL-S的1.53、1.50和2.22倍,而CYP9YaCYP9Yc在TC种群中的表达量均显著降低,表达量是QPQL-S的49.45%和52.55%(图4-A);2019年成虫中,与相对敏感种群URMQ2相比,YN2种群CYP4C1CYP9Yc的表达最高且显著增加,表达量分别为URMQ2的5.56、4.03倍,CYP305a1的表达量在JMS中最高,显著高于URMQ2(6.6倍),其次是QPQL2,其他4个种群的表达量较低;CYP9Z25在YN3的表达量最高,显著上调,表达量是URMQ2的7.49倍,其次是URMQ4和YN2种群,JMS、QPQL2和URMQ3种群的表达量较低;CYP9e2的表达量在各种群中无显著变化(图4-B)。

图4

新窗口打开|下载原图ZIP|生成PPT
图4马铃薯甲虫田间种群成虫P450基因的相对表达量(A:2018;B:2019)

Fig. 4Relative expression of P450 genes in field populations of L. decemlineata adults from Xinjiang in 2018 (A) and 2019 (B)



2.5 噻虫嗪胁迫下马铃薯甲虫P450基因的表达分析

将2019年马铃薯甲虫相对敏感种群URMQ2的4龄幼虫和成虫经噻虫嗪LD50处理72 h,利用qRT-PCR测定其P450基因的表达变化,结果如图5所示。4龄幼虫中,以未处理同日龄的试虫为对照,处理组CYP9e2、CYP305a1、CYP9Ya、CYP9Yc的表达量显著上升,分别为对照的5.08、4.25、2.38和2.79倍;而成虫处理组CYP9Ya、CYP9Yc显著上调,表达量分别为对照的1.80和1.68倍,CYP9e2、CYP305a1显著下调。可以看出,经噻虫嗪LD50处理均可诱导4龄幼虫和成虫CYP9Ya、CYP9Yc表达量显著增加,CYP9e2、CYP305a1在4龄幼虫和成虫的表达变化则相反,而CYP4C1、CYP9Z25的表达量均无显著变化。

图5

新窗口打开|下载原图ZIP|生成PPT
图5相对敏感种群4龄幼虫和成虫经噻虫嗪处理后P450基因的相对表达量

L4:4龄幼虫4th instar larva;A:成虫 Adult;TMX:噻虫嗪Thiamethoxam
Fig. 5Relative expression of P450 genes in 4th instar larvae and adults of relatively sensitive population of L. decemlineata exposed to thiamethoxam



2.6 新疆马铃薯甲虫对噻虫嗪的抗性水平与P450基因表达量的相关性分析

利用SPSS分析2018年和2019年新疆马铃薯甲虫不同种群对噻虫嗪的抗性水平与上述P450基因表达量之间的相关性,结果如表7所示,马铃薯甲虫成虫对噻虫嗪的抗性水平与CYP9Ya表达量之间的相关系数(R)为0.810(P=0.027),二者之间为高度显著正相关。

Table 7
表7
表7新疆马铃薯甲虫对噻虫嗪的抗性水平与P450基因表达量的相关性分析
Table 7Correlation analysis between expression of P450 genes and resistance level to thiamethoxam of L. decemlineata from Xinjiang
采集时间
Sampling time
基因
Gene
RR value PP value
4龄幼虫4th instar larva成虫Adult4龄幼虫4th instar larva成虫Adult
2018CYP9Ya0.5490.5480.3380.339
CYP9Yc0.7900.4870.1120.406
CYP9Z250.0490.3090.9380.613
2019CYP9Ya-0.0960.810*0.8380.027
CYP9Yc-0.1010.6460.8290.117
CYP9Z250.2700.1890.5580.685
CYP4C10.6670.6790.1020.094
CYP9e2-0.1980.4260.6700.340
CYP305a10.1440.1790.7580.702

新窗口打开|下载CSV

3 讨论

马铃薯甲虫对杀虫剂容易产生抗性,对于一种新的杀虫剂,往往在使用2—4年就产生明显的抗药性[2]。噻虫嗪是新疆地区防治马铃薯甲虫的主推药剂之一,随着使用周期的累积,其抗性发展不可避免。2009、2010年的调查发现新疆马铃薯甲虫4龄幼虫对噻虫嗪普遍敏感,仅2010年特克斯(TKS)种群产生8.03倍的低水平抗性[14];同时马铃薯甲虫成虫的抗性监测显示2009年昌吉(CJ)种群对噻虫嗪具有5.22倍的低水平抗性,2010年CJ、QPQL、TC和URMQ1种群对噻虫嗪发展了低到中等水平抗性(7.48、5.98、10.15和10.39倍)[9]。本研究抗性监测表明2018年YN1、URMQ1、TC和QPQL1以及2019年YN2、JMS田间种群成虫对噻虫嗪产生了5.99—20.33倍的低至中水平抗性。其中2019年YN2种群对噻虫嗪的抗性水平比2018年YN1种群增加了1倍。对上述采集地用药情况的调查发现,察布查尔、伊宁地区近3年均施用噻虫嗪、吡虫啉等药剂,吉木萨尔县近10年持续使用阿立卡(12.6%噻虫嗪·9.4%高效氯氟氰菊酯)防治马铃薯甲虫。可见长期、频繁的用药导致当地种群对新烟碱类杀虫剂的敏感性明显下降。另外,还可以发现本文的监测结果与刘萍等[9]、卢伟平等[14]的调查结果具有相似的趋势,即成虫对噻虫嗪的敏感性较低,而幼虫的敏感性较高,且随着时间发展田间最高抗性水平有所上升。虽然新疆马铃薯甲虫对噻虫嗪的抗性还处于低至中等水平,但其发展不容忽视,注意合理轮换使用不同作用机制的药剂,尽量避免在高龄期或成虫期用药,同时开展持续的抗性监测及抗性机制研究是马铃薯甲虫抗性治理的重要前提。

细胞色素P450单加氧酶是马铃薯甲虫代谢吡虫啉的主要酶系。ZHAO等[6]研究发现,胡椒基丁醚在纽约州长岛马铃薯甲虫成虫和幼虫中对吡虫啉都有明显增效作用;MOTA-SANCHEZ等[7,21]在其他种群的增效研究也获得同样的结果。在本研究中,噻虫嗪LD50处理72 h都可引起敏感和抗性马铃薯甲虫4个种群的P450活性显著升高,另外URMQ2和YN2种群的EST及URMQ3种群的GST活性也显著增加,这表明P450可能在马铃薯甲虫对噻虫嗪的解毒代谢中起主要作用。类似的研究发现褐飞虱(Nilaparvata lugens)噻虫嗪抗性种群的细胞色素P450活性也显著高于室内敏感品系(2.13倍)[22]

近年来随着分子生物学及生物信息学的快速发展,有关昆虫对新烟碱类杀虫剂抗性的解毒代谢分子机制也取得了一系列进展[10,11,12,13]。为了筛查新疆马铃薯甲虫对噻虫嗪抗性相关的解毒酶基因,本研究就马铃薯甲虫噻虫嗪抗性和敏感成虫进行转录组分析,结果发现与抗性相关的显著差异表达的上调解毒酶基因有P450基因2个、尿苷二磷酸-葡萄糖基转移酶基因1个、ABC转运蛋白基因3个,但未见有显著差异表达的酯酶和谷胱甘肽S-转移酶基因。CLEMENTS等[23]将吡虫啉处理后的马铃薯甲虫与对照进行转录组测序,按照解毒机制进行分类,发现吡虫啉可诱导4个P450基因、1个尿苷二磷酸-葡萄糖基转移酶基因以及1个羧酸酯酶基因上调。可见不同背景的种群其解毒酶基因的表达量变化不同。本文利用qRT-PCR验证转录组数据中3个上调的P450基因CYP4C1、CYP9e2CYP305a1在抗性种群的表达量,结果与DEG数据基本一致,另外结合ZHU等[18]在美国长岛马铃薯甲虫吡虫啉抗性种群中发现的过表达基因CYP9Z25、CYP9YaCYP9Yc,本研究进一步对上述6个P450基因在新疆马铃薯甲虫不同田间种群的表达进行了比较分析,结果表明,4龄幼虫中,对噻虫嗪敏感性降低的QPQL1种群中CYP9YaCYP9Yc的表达量显著增加;成虫中,对噻虫嗪低水平抗性的URMQ1种群中CYP9Z25高表达,CYP9YaCYP9Yc、CYP4C1CYP9Yc以及CYP305a1分别在QPQL1、YN2以及JMS种群中显著上调,这些种群对噻虫嗪达到中水平抗性。综上,可以看出CYP9YaCYP9Yc在敏感性降低的幼虫和抗性成虫种群中表达量均显著提高,除CYP9e2以外其他5个基因在抗性成虫中都显著上调,表明这些基因可能参与了新疆马铃薯甲虫对噻虫嗪抗性的产生。

吡虫啉对昆虫P450基因表达具有诱导作用。ZHU等[18]研究显示对吡虫啉约30倍抗性的美国长岛马铃薯甲虫种群中有21个P450基因可被吡虫啉诱导上调表达;吡虫啉处理蜜蜂可引起6个P450基因的表达量增加[24];大草蛉(Chrysopa pallens)幼虫经吡虫啉诱导4 h其体内P450基因表达量略微上调,16 h后转录水平达最高值,随后逐渐下降[25]。本文qRT-PCR结果表明,经噻虫嗪LD50处理72 h,URMQ2种群4龄幼虫和成虫的CYP9YaCYP9Yc均显著上调表达。ZHANG等研究发现褐飞虱中13个可被吡虫啉诱导的P450基因在抗性种群中却无差异表达,而CYP4FB2、CYP439B1CYP6BD12不能被吡虫啉诱导,但在抗性种群中表达上调[26],这表明可被药剂诱导的基因与抗性相关的过表达基因并非完全一致。因此,综合本研究中诱导试验及不同田间种群P450基因的表达分析结果,在噻虫嗪抗性种群中显著上调同时又能被噻虫嗪诱导的基因即为CYP9YaCYP9Yc。另外相关性分析显示成虫对噻虫嗪的抗性与CYP9Ya表达量呈显著正相关。而真正参与马铃薯甲虫对噻虫嗪抗性形成的基因有待进一步的功能验证。

4 结论

通过抗性监测发现2018年YN1、URMQ1、TC和QPQL1以及2019年YN2和JMS种群成虫对噻虫嗪产生了低至中水平抗性。噻虫嗪处理可引起测试种群P450活性的升高,表明其在马铃薯甲虫对噻虫嗪抗性产生中可能起到重要作用。抗性和敏感种群的转录组分析以及qRT-PCR结果显示包括CYP9YaCYP9Yc的5个P450基因在抗性种群中显著上调,同时噻虫嗪又可诱导这两个基因在4龄幼虫和成虫中的表达。结合相关性分析推测CYP9Ya与马铃薯甲虫对噻虫嗪的抗性形成密切相关,但不排除其他P450基因以及其他解毒酶基因的作用。

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

ALYOKHIN A. Colorado potato beetle management on potatoes: Current challenges and future prospects//Fruit,
Vegetable and Cereal Science and Biotechnology, 2009, 3(1):10-19.

[本文引用: 1]

ALYOKHIN A, BAKER M, MOTA-SANCHEZ D, DIVELY G, GRAFIUS E. Colorado potato beetle resistance to insecticides
American Journal of Potato Research, 2008, 85:395-413.

DOI:10.1007/s12230-008-9052-0URL [本文引用: 2]

郭文超, 吐尔逊, 许建军, 刘建, 何江, 李晶, 马德成, 王俊. 马铃薯甲虫识别及其在新疆的分布、传播和危害
新疆农业科学, 2010, 47(5):906-909.

[本文引用: 1]

GUO W C, TU E R, XU J J, LIU J, HE J, LI J, MA D C, WANG J. Research on the identification of Colorado potato beetle and its distribution, dispersal and damage in Xinjiang
Xinjiang Agricultural Sciences, 2010, 47(5):906-909. (in Chinese)

[本文引用: 1]

须志平. 新烟碱类杀虫剂在作物保护方面的应用
世界农药, 2009, 31(1):18-21.

[本文引用: 1]

XU Z P. Application of neonicotinoids insecticides in crop protection
World Pesticides, 2009, 31(1):18-21. (in Chinese)

[本文引用: 1]

王彦华, 王鸣华. 害虫对噻虫嗪抗药性及其治理
世界农药, 2008, 30(4):42-45.

[本文引用: 1]

WANG Y H, WANG M H. Pest resistance to thiamethoxam and its management
World Pesticides, 2008, 30(4):42-45. (in Chinese)

[本文引用: 1]

ZHAO J Z, BISHOP B A, GRAFIUS E J. Inheritance and synergism of resistance to imidacloprid in the Colorado potato beetle (Coleoptera: Chrysomelidae)
Journal of Economic Entomology, 2000, 93(5):1508-1514.

DOI:10.1603/0022-0493-93.5.1508URL [本文引用: 2]

MOTA-SANCHEZ D, HOLLINGWORTH R M, GRAFIUS E J, MOYER D D. Resistance and cross-resistance to neonicotinoid insecticides and spinosad in the Colorado potato beetle, Leptinotarsa decemlineata (Say) (Coleoptera: Chrysomelidae)
Pest Management Science, 2006, 62(1):30-37.

DOI:10.1002/(ISSN)1526-4998URL [本文引用: 2]

王志田, 姜卫华, 李国清. 新疆北疆马铃薯甲虫成虫抗药性水平监测
农药, 2010, 49(3):206-208.

[本文引用: 1]

WANG Z T, JIANG W H, LI G Q. Insecticide resistance in adult of the Colorado potato beetle among north Xinjiang Uygur autonomous region
Agrochemicals, 2010, 49(3):206-208. (in Chinese)

[本文引用: 1]

刘萍, 姜卫华, 卢伟平, 李国清. 新疆北疆马铃薯甲虫成虫对新烟碱类杀虫剂的敏感性变化
农药学学报, 2011, 13(3):271-275.

[本文引用: 4]

LIU P, JIANG W H, LU W P, LI G Q. Susceptibility of Colorado potato beetle Leptinotarsa decemlineata adults from northern Xinjiang Uygur autonomous region to 4 neonicotinoids
Chinese Journal of Pesticide Science, 2011, 13(3):271-275. (in Chinese)

[本文引用: 4]

KAPLANOGLU E, CHAPMAN P, SCOTT I M, DONLY C. Overexpression of a cytochrome P450 and UDP-glycosyltransferase is associated with imidacloprid resistance in the Colorado potato beetle, Leptinotarsa decemlineata
Scientific Reports, 2017, 7:1762.

DOI:10.1038/s41598-017-01961-4URL [本文引用: 2]

CLEMENTS J, SCHOVILLE S, PETERSON N, LAN Q, GROVES R L. Characterizing molecular mechanisms of imidacloprid resistance in select populations of Leptinotarsa decemlineata in the central sands region of Wisconsin
PLoS ONE, 2016, 11(1):e0147844.

DOI:10.1371/journal.pone.0147844URL [本文引用: 2]

CLEMENTS J, SCHOVILLE S, PETERSON N, HUSETH A S, LAN Q, GROVES R L. RNA interference of three up-regulated transcripts associated with insecticide resistance in an imidacloprid resistant population of Leptinotarsa decemlineata
Pesticide Biochemistry and Physiology, 2017, 135:35-40.

DOI:10.1016/j.pestbp.2016.07.001URL [本文引用: 2]

CLEMENTS J, OLSON J M, SANCHEZ-SEDILLO B, BRADFORD B, GROVES R L. Changes in emergence phenology, fatty acid composition, and xenobiotic-metabolizing enzyme expression is associated with increased insecticide resistance in the Colorado potato beetle
Archives of Insect Biochemistry and Physiology, 2020, 103(3):e21630.

[本文引用: 2]

卢伟平, 刘萍, 姜卫华, 李国清. 4种新烟碱类杀虫剂对马铃薯甲虫的触杀毒力比较
农药, 2011, 50(2):137-140.

[本文引用: 3]

LU W P, LIU P, JIANG W H, LI G Q. Comparison of contact toxicities of 4 neonicotinoids against Leptinotarsa decemlineata
Agrochemicals, 2011, 50(2):137-140. (in Chinese)

[本文引用: 3]

BRADFORD M M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding
Analytical Biochemistry, 1976, 72:248-254.

DOI:10.1016/0003-2697(76)90527-3URL [本文引用: 1]

祝海栋, 李瑞琳, 何小雨, 赵丹, 韩鑫胤, 牛北方. 六倍体小麦基因组注释流程构建与优化
计算机系统应用, 2019, 28(8):222-228.

[本文引用: 1]

ZHU H D, LI R L, HE X Y, ZHAO D, HAN X Y, NIU B F. Construction and optimization of hexaploid wheat genome annotation process
Computer Systems and Applications, 2019, 28(8):222-228. (in Chinese)

[本文引用: 1]

LANGMEAD B, TRAPNELL C, POP M, SALZBERG S L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome
Genome Biology, 2009, 10(3):R25.

DOI:10.1186/gb-2009-10-3-r25URL [本文引用: 1]

ZHU F, MOURAL T W, NELSON D R, PALLI S R. A specialist herbivore pest adaptation to xenobiotics through up-regulation of multiple cytochrome P450s
Scientific Reports, 2016, 6:20421.

DOI:10.1038/srep20421URL [本文引用: 3]

王彦华, 陈进, 沈晋良, 高聪芬, 黄悦, 张久双, 李文红, 周威君. 防治褐飞虱的高毒农药替代药剂的室内筛选及交互抗性研究
中国水稻科学, 2008, 22(5):519-526.

[本文引用: 1]

WANG Y H, CHEN J, SHEN J L, GAO C F, HUANG Y, ZHANG J S, LI W H, ZHOU W J. Laboratory screening and cross-resistance analysis of alternative insecticides for highly-toxic pesticides for controlling brown planthopper, Nilaparvata lugens
Chinese Journal of Rice Science, 2008, 22(5):519-526. (in Chinese)

[本文引用: 1]

LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method
Methods, 2001, 25(4):402-408.

DOI:10.1006/meth.2001.1262URL [本文引用: 1]

MOTA-SANCHEZ D, WHALON M, GRAFIUS E, HOLLINGWOETH R. Resistance of Colorado potato beetle to imidacloprid
Resistant Pest Management, 2000, 11(1):31-34.

[本文引用: 1]

张钰明, 向兴, 王学贵. 噻虫嗪对褐飞虱的毒力及解毒代谢酶活性的影响
华南农业大学学报, 2020, 41(3):80-85.

[本文引用: 1]

ZHANG Y M, XIANG X, WANG X G. Effects of thiamethoxam on toxicity and detoxification metabolic enzyme activity of Nilaparvata lugens
Journal of South China Agricultural University, 2020, 41(3):80-85. (in Chinese)

[本文引用: 1]

CLEMENTS J, SANCHEZ-SEDILLO B, BRADFIELD C A, GROVES R L. Transcriptomic analysis reveals similarities in genetic activation of detoxification mechanisms resulting from imidacloprid and chlorothalonil exposure
PLoS ONE, 2018, 13(10):e0205881.

DOI:10.1371/journal.pone.0205881URL [本文引用: 1]

WU M C, CHANG Y W, LU K H, YANG E C. Gene expression changes in honey bees induced by sublethal imidacloprid exposure during the larval stage
Insect Biochemistry and Molecular Biology, 2017, 88:12-20.

DOI:10.1016/j.ibmb.2017.06.016URL [本文引用: 1]

刘昌燕, 曾凡荣. 大草蛉P450基因cDNA片段的克隆及吡虫啉诱导表达
中国生物防治学报, 2014, 30(3):427-433.

[本文引用: 1]

LIU C Y, ZENG F R. Cloning of Chrysopa pallens P450 cDNA fragment and the expression induced by imidacloprid
Chinese Journal of Biological Control, 2014, 30(3):427-433. (in Chinese)

[本文引用: 1]

ZHANG Y X, YANG Y X, SUN H H, LIU Z W. Metabolic imidacloprid resistance in the brown planthopper, Nilaparvata lugens, relies on multiple P450 enzymes
Insect Biochemistry and Molecular Biology, 2016, 79:50-56.

DOI:10.1016/j.ibmb.2016.10.009URL [本文引用: 1]

相关话题/基因 新疆 序列 数据 种群