枸杞棉蚜抗、感吡虫啉品系转录组比较分析  

Comparative Transcriptome Analysis of Resistant and Sensitive Strains to Imidacloprid of Aphis gossypii(Hemiptera:Aphisidae)

作　　者：王新霞 张慧斌 刘雲祥 刘思雨 来有鹏[1] 李秋荣[1] WANG Xinxia;ZHANG Huibin;LIU Yunxiang;LIU Siyu;LAI Youpeng;LI Qiurong(Academy of Agriculture and Forestry Sciences,Qinghai University,State Scientific Observing and Experimental Station of Crop Pest in Xining,Ministry of Agriculture and Rural Affairs,Provincial Key Laboratory of Agricultural Integrated Pest Management in Qinghai,Xining 810016,China)

机构地区：[1]青海大学农林科学院,农业农村部西宁作物有害生物科学观测实验站,青海省农业有害生物综合治理重点实验室,青海西宁810016

出　　处：《华北农学报》2025年第1期177-188,共12页Acta Agriculturae Boreali-Sinica

基　　金：青海省“昆仑英才·乡村振兴人才”项目;青海省自然科学基金项目(2017-ZJ-934Q)。

摘　　要：为挖掘枸杞棉蚜对杀虫剂吡虫啉产生抗药性的候选基因,以枸杞棉蚜对吡虫啉的室内抗药性品系和相对敏感性品系为材料,利用Illumina HiSeq 2500高通量测序技术获得2个品系的转录组数据,通过NCBI数据库进行基因注释,在转录水平上对2个品系的差异基因GO功能及KEGG代谢途径等生物信息学进行分析,采用qRT-PCR技术检测其中8个候选差异基因(CYP6a2、CYP6a13、CYP6k1、CYP6j1、CYP4c1、AChE2、CarE和ALP3)的相对表达情况,并分析了抗药性相关基因的进化关系。经测序及序列拼接,共获得70 101条Unigene(单基因簇),平均长度为654.37 bp,在NR、GO和KEGG数据库中分别注释到29 131,27 861,2 993条Unigene。通过对两品系的差异基因进行NR注释分析,共发现289个差异基因,其中,22个可能与抗药性相关,包括9个解毒酶系基因(CYP6a13、CYP6k1、CYP6j1、CYP4c1和ALP3各1个,CYP6a2和CarE各2个),8个表皮蛋白基因(CP)及其前体(Cuticular protein precursor, CPP),1个靶标酶系基因(Alkaline phosphatase,AChE2),2个转录因子基因(WRKY1、LZTFL1,Leucine zipper transcription factor-like protein 1 gene),1个胰脂肪酶相关蛋白2基因(PLRP2)和1个多抗相关蛋白基因(Multidrug resistance-associated protein, MRP)。定量结果表明,差异基因CYP6a2、CYP6a13、CYP6k1、CYP6j1、CarE和ALP3在抗药性品系的表达量均显著高于敏感性品系。对抗药性相关重要基因CYP、GST、ALP和CP构建系统进化树,通过分析基因的遗传关系,发现枸杞棉蚜与大豆蚜、豌豆蚜的亲缘关系相对较近。获得了枸杞棉蚜对吡虫啉的抗药性与敏感性品系转录组的整体表达模式,筛选到抗药性相关差异表达基因,发现2个品系中87.50%的候选差异基因在转录水平和mRNA水平的表达量变化一致。To explore candidate genes for resistance of Aphis gossypii to the insecticide imidacloprid,the transcriptome data of the two strains were obtained and compared with imidacloprid indoor resistant and relatively sensitive lines of Aphis gossypii by using Illumina HiSeq 2500 high-throughput sequencing technology.The gene annotation was performed using the NCBI database,and bioinformatics analysis of the differential genes at the transcriptional level included GO function,KEGG metabolic pathway,and other analyses.The relative expression of eight candidate differentially expressed genes(CYP6a2,CYP6a13,CYP6k1,CYP6j1,CYP4c1,AChE2,CarE and ALP3)was detected using qRT-PCR technology,and the evolutionary relationships of resistance related genes were analyzed.After sequencing and sequence splicing,a total of 70101 Unigenes were obtained,with an average length of 654.37 bp.29131,27861 and 2993 Unigenes were annotated in NR,GO and KEGG databases,respectively.According to the NR annotation analysis of the differential genes of insecticide resistance and sensitivity strains,a total of 22 differential genes that may be related to insecticide resistance were found,including 9 detoxification enzyme genes(CYP6a13,CYP6k1,CYP6j1,CYP4c1 and ALP3 each,two CYP6a2,two CarE),8 cuticle protein genes(CP)and their precursors(CPP),target enzyme genes(AChE2),2 transcription factors(WRKY1,leucine zipper transcription factor-like protein 1 gene,LZTFL1),one pancreatic lipase-related protein 2 gene(PLRP2)and one multidrug resistance-associated protein gene(MRP).The results of qRT-PCR indicated that the expression levels of CYP6a2,CYP6a13,CYP6k1,CYP6j1,CarE and ALP3 genes in insecticide resistance strains were significantly higher than those in sensitive strains.The phylogenetic tree analysis of CYP,ALP,GST and CP genes obtained by NR annotation indicated that the genetic relationship between the Aphis gossypii and Aphis glycines and Acyrthosiphon pisum was relatively close.We found that 87.50%of the candidate differentially expressed genes in th

关 键 词：枸杞棉蚜 吡虫啉 抗药性 敏感性 转录组 差异表达基因 

分 类 号：Q963[生物学—昆虫学] S433.3[农业科学—农业昆虫与害虫防治]