收稿日期: 2020-01-17
网络出版日期: 2020-09-09
基金资助
国家科技重大专项(2017ZX10303404)
Analysis of deltamethrin resistance-related genes based on the transcriptome of Culex pipiens quinquefasciatus
Received date: 2020-01-17
Online published: 2020-09-09
Supported by
National Science and Technology Major Project(2017ZX10303404)
目的 通过对致倦库蚊溴氰菊酯敏感株和抗性株进行转录组分析,筛选可能与抗性相关的差异表达基因。方法 对致倦库蚊实验室敏感株（SS株）、海南株（HN株）和溴氰菊酯药筛抗性株（RR株）3个不同抗性水平的株系进行生物测定,分别构建3个株系的cDNA文库,提取RNA,进行转录组测序。基于FPKM值进行差异性分析,对差异基因进行GO和KEGG功能注释与富集分析。筛选出致倦库蚊不同株系中表达量有显著性差异的基因,并采用实时定量PCR验证。结果 生物测定结果显示,与致倦库蚊SS株相比,HN株和RR株的抗性倍数分别达到了4 828倍和197 241倍,为高等抗性水平。RNA-seq测序结果显示,致倦库蚊SS株、HN株和RR株转录组共注释基因转录本197 079个。GO富集分析结果显示,在结合、催化活性、细胞过程、代谢过程、单一生物过程、细胞、细胞连接、细胞器中基因富集数量较多。KEGG富集分析结果显示,在富集度最高的前20个条目中,RR株主要富集于剪接体、核糖体等通路,HN株则主要富集于剪切体、RNA转运、hedgehog信号通路、细胞周期、核苷酸切除修复、DNA复制等方向。基于FPKM值的差异性分析结果显示,致倦库蚊RR株和HN株较SS株均上调的有5 357个,均下调的有10 978个。与溴氰菊酯抗性相关基因的靶向性筛选结果显示,与SS株相比,RR株和HN株注释的细胞色素P450基因、谷胱甘肽转移酶基因、气味结合蛋白基因、气味受体基因、视觉基因、核糖体蛋白基因和碳酸酐酶基因转录本中,均上调的分别为100、17、13、4、16、191和1个,均下调的分别为253、48、43、4、15、504和9个。从上调基因中挑选的12个基因（2个细胞色素P450基因、1个谷胱甘肽转移酶基因、1个视觉基因、1个核糖体蛋白基因、1个碳酸酐酶基因、1个核酸内切酶基因、2个气味结合蛋白基因、3个气味受体基因）进行实时定量PCR验证结果显示,在RR株和HN株中,这12个基因较SS株均显著上调,其中RR株较SS株上调1.62~80.26倍,HN株较SS株上调2.01~189.48倍。相关性分析结果显示,12个目标基因的RNA-seq结果与实时定量PCR结果在HN株和RR株中均呈正相关（r = 0.802、0.710）。证明这12个目标基因的实时定量PCR结果与转录组测序结果趋势一致。结论 基于高通量测序结果和定量PCR验证分析,致倦库蚊中的某些细胞色素P450基因、谷胱甘肽S转移酶基因、视觉基因、嗅觉蛋白基因可能在其对溴氰菊酯的抗性中起一定作用。
沈瑞鑫 , 王意婷 , 李春晓 , 吴家红 , 赵彤言 , 陈艳 . 基于致倦库蚊转录组的溴氰菊酯抗性相关基因分析[J]. 中国寄生虫学与寄生虫病杂志, 2020 , 38(4) : 453 -463 . DOI: 10.12140/j.issn.1000-7423.2020.04.010
Objective To screen for differentially-expressed genes related to deltamethrin resistance by transcriptome analysis of deltamethrin-susceptible and -resistant strains of Culex pipiens quinquefasciatus.Methods Bioassays were performed on deltamethrin-susceptible (SS) strain, Hainan (HN) strain, and deltamethrin-resistant (RR) strain of Cx. pipiens quinquefasciatus. cDNA libraries of each strain were constructed for transcript sequencing. The differentially-expressed genes obtained by differential analysis base on the FPKM value (number of fragments per kilobase of transcript sequence per millions base pairs sequenced) were subjected to GO enrichment analysis and KEGG (kyoto encyclopedia of genes and genomes) functional annotation and enrichment analysis. Genes with significant differences in expression level were screened out and verified by real-time quantitative PCR.Results The resistance levels of the HN and RR strains were 4 828 and 197 241 folds of the SS strain, respectively, both showing high resistance to deltamethrin. A total of 197 079 transcripts were annotated among the transcriptomes of the three Cx. pipiens quinquefasciatus strains. The results of GO enrichment analysis showed enrichment of genes in terms of binding, catalytic activity, cellular process, metabolic process, single organism process, cellular process, metabolic process, cell, cell junction, and organelle. The results of KEGG enrichment analysis showed that among the top 20 items, the RR strain showed enrichment of genes in terms of spliceosome and ribosome pathways, while the HN strain showed significant enrichment of genes in spliceosome, RNA transport, hedgehog signaling pathway, cell cycle, nucleotide excision repair, and DNA replication.Results of differential analysis based on FPKM value showed that there were 5 357 up-regulated and 10 978 down-regulated genes in both the RR and HN strains, compared to the SS strain. The targeting screening for genes related the genes with differential expression in RR and HN strains including cytochrome P450 (P450), glutathione transferase (GST), odorant binding protein (OBP), odorant receptor (OR), opsin gene, ribosomal protein gene, and carbonic anhydrase, in comparison with the SS strain. The annotated transcriptomes of these differential genes in RR and HN strains were found having 100, 17, 13, 4, 16, 191 and 1 up-regulated transcripts, and 253, 48, 43, 4, 15, 504 and 9 down-regulated transcripts, respectively, compared to the SS strain. To verify, twelve up-regulated genes were selected for real-time quantitative PCR, including 2 P450 genes, 1 GST gene, 1 opsin gene, 1 ribosomal protein gene, 1 carbonic anhydrase (CA) gene, 1 nucleic acid endonuclease gene, 2 OBP genes and 3 OR genes. The results showed that in RR and HN strain, the 12 genes were all significantly up-regulated compared to the SS strain, among them the expressions of those in the RR strain were up-regulated by 1.62-80.26 folds, and those in the HN strain up-regulated by 2.01-189.48 folds. The correlation analysis on the 12 targeting genes showed a positive correlation between the results of RNA-seq and the results of real-time quantitative PCR for NH and RR strain, with a correlation coefficient r of 0.802 and 0.710, respectively. It demonstrates that for the 12 genes, the results of real-time quantitative PCR were consistent with the transcriptome sequencing data.Conclusion Based on the high-throughput sequencing results and quantitative PCR verification analysis, some genes including cytochrome P450 genes, glutathione S-transferase genes, visual genes, and olfactory protein genes in Cx. pipiens quinquefasciatus may play certain roles in the deltamethrin resistance.
| [1] | Nasci RS, Miller BR. Culicine mosquitoes and the agents they transmit[J]. Biol Dis Vectors, 1996,18(2):85-97. |
| [2] | Arensburger P, Megy K, Waterhouse RM, et al. Sequencing of Culex pipiens quinquefasciatus establishes a platform for mosquito comparative genomics[J]. Science, 2010,330(6000):86-88. |
| [3] | Liu NN, Liu HQ, Zhu F, et al. Differential expression of genes in pyrethroid resistant and susceptible mosquitoes, Culex quinquefasciatus (S.)[J]. Gene, 2007,394(1/2):61-68. |
| [4] | Ding YR, Yan ZT, Si FL, et al. Mitochondrial genes associated with pyrethroid resistance revealed by mitochondrial genome and transcriptome analyses in the malaria vector Anopheles sinensis (Diptera ∶ Culicidae)[J]. Pest Manag Sci, 2020,76(2):769-778. |
| [5] | Liu BQ, Chen B, Qiao L. The progress of quantitative trait loci for mosquitoes pyrethroid resistance[J]. J Chongqing Norm Univ Nat Sci, 2016,33(6):26-31. (in Chinese) |
| [5] | ( 刘柏琦, 陈斌, 乔梁. 蚊虫抗拟除虫菊酯数量性状位点的研究进展[J]. 重庆师范大学学报(自然科学版), 2016,33(6):26-31.) |
| [6] | Cai YS, Deng ML, Yuan MT, et al. Resistance of Culex quinquefasciatus to normal used insceticides in Guang’an City[J]. Chin J Hyg Insect Equip, 2018,24(4):410-411. (in Chinese) |
| [6] | ( 蔡运山, 邓茂玲, 袁茂涛, 等. 广安市致倦库蚊对常用杀虫剂的抗性监测[J]. 中华卫生杀虫药械, 2018,24(4):410-411.) |
| [7] | Yang XY, Liu C, Chen XW, et al. Resistance to 3 kind of pyrethroid insecticides in Culex pipiens quinquefasciatus from Haikou[J]. China Trop Med, 2019,19(11):1092-1094. (in Chinese) |
| [7] | ( 杨新艳, 刘超, 陈学文, 等. 海口市致倦库蚊成蚊对3种菊酯类杀虫剂的抗性调查[J]. 中国热带医学, 2019,19(11):1092-1094.) |
| [8] | Liang QG, Wang ZY, Yang X, et al. Investigation on resistance of Culex pipiens quinquefasciatus to four commonly used insecticides in Guiyang[J]. Chin J Endem, 2019,8(6):476-480. (in Chinese) |
| [8] | ( 梁秋果, 王政艳, 杨茜, 等. 贵阳市致倦库蚊成蚊对4种常用杀虫剂抗药性调查[J]. 中华地方病学杂志, 2019,8(6):476-480.) |
| [9] | Liu Y, Zhang SH, Qin YM, et al. Investigation on current status of resistance of Culex pipiens quinquefasciatus to commonly used insecticides in Shenzhen City, Guangdong Province[J]. Chin J Vect Biol Contr, 2020,31(3):362-365. (in Chinese) |
| [9] | ( 刘阳, 张韶华, 秦彦珉, 等. 广东省深圳市致倦库蚊对常用杀虫剂抗药性现状调查[J]. 中国媒介生物学及控制杂志, 2020,31(3):362-365.) |
| [10] | Phillips RS. Current status of malaria and potential for control[J]. Clin Microbiol Rev, 2001,14(1):208-226. |
| [11] | Sawicki RM, Denholm I. Adaptation of insects to insecticides[M] //Evered D, Collins GM. Ciba Foundation Symposium--Origins and Development of Adaptation. Chichester, UK: John Wiley & Sons, Ltd., 2008: 152-166. |
| [12] | Brattsten LB, Holyoke CW, Leeper JR, et al. Insecticide resistance: challenge to pest management and basic research[J]. Science, 1986,231(4743):1255-1260. |
| [13] | Xu X, Qian K. Research progress on the resistance of mosquitoes to pyrethroid insecticides[J]. Cap J Pub Health, 2018,12(1):9-12. (in Chinese) |
| [13] | ( 徐鑫, 钱坤. 蚊虫对拟除虫菊酯类杀虫剂抗性研究进展[J]. 首都公共卫生, 2018,12(1):9-12.) |
| [14] | Wang CB, Lu WH, Lin Y, et al. Development and application of transcriptome sequencing[J]. Eucalypt Sci Technol, 2018,35(4):20-26. (in Chinese) |
| [14] | ( 王楚彪, 卢万鸿, 林彦, 等. 转录组测序的发展和应用[J]. 桉树科技, 2018,35(4):20-26.) |
| [15] | Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics[J]. Nat Rev Genet, 2009,10(1):57-63. |
| [16] | Br?utigam A, Gowik U. What can next generation sequencing do for you? Next generation sequencing as a valuable tool in plant research[J]. Plant Biol (Stuttg), 2010,12(6):831-841. |
| [17] | Bonizzoni M, Ochomo E, Dunn WA, et al. RNA-seq analyses of changes in the Anopheles gambiae transcriptome associated with resistance to pyrethroids in Kenya: identification of candidate-resistance genes and candidate-resistance SNPs[J]. Parasit Vectors, 2015,8:474. |
| [18] | Lv Y, Wang WJ, Hong SC, et al. Comparative transcriptome analyses of deltamethrin-susceptible and-resistant Culex pipiens pallens by RNA-seq[J]. Mol Genet Genomics, 2016,291(1):309-321. |
| [19] | Li CX, Guo XX, Zhang YM, et al. Identification of genes involved in pyrethroid-, propoxur-, and dichlorvos-insecticides resistance in the mosquitoes, Culex pipiens complex (Diptera ∶ Culicidae)[J]. Acta Trop, 2016,157:84-95. |
| [20] | Bonizzoni M, Afrane Y, Dunn WA, et al. Comparative transcriptome analyses of deltamethrin-resistant and-susceptible Anopheles gambiae mosquitoes from Kenya by RNA-seq[J]. PLoS One, 2012,7(9):e44607. |
| [21] | WHO. Report of the WHO Informal Consultationon the evaluation and testing of insecticides[R]. Geneva: WHO, 1996: 7-11. |
| [22] | Mortazavi A, Williams BA, McCue K, et al. Mapping and quantifying mammalian transcriptomes by RNA-seq[J]. Nat Methods, 2008,5(7):621-628. |
| [23] | Garber M, Grabherr MG, Guttman M, et al. Computational methods for transcriptome annotation and quantification using RNA-seq[J]. Nat Methods, 2011,8(6):469-477. |
| [24] | Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol, 2014,15(12):550. |
| [25] | Young MD, Wakefield MJ, Smyth GK, et al. Gene ontology analysis for RNA-seq: accounting for selection bias[J]. Genome Biol, 2010,11(2):R14. |
| [26] | Sun HB, Sun LC, He J, et al. Cloning and characterization of ribosomal protein S29, a deltamethrin resistance associated gene from Culex pipiens pallens[J]. Parasitol Res, 2011,109(6):1689-1697. |
| [27] | Yu J, Hu SL, Ma K, et al. Ribosomal protein S29 regulates metabolic insecticide resistance through binding and degradation of CYP6N3[J]. PLoS One, 2014,9(4):e94611. |
| [28] | Liu HQ, Liu YB, Song YP, et al. Cloning and bioinformatics analysis of full length cDNA of permethrin-resistance associated opsin gene PR-OP of Culex pipens pallens[J]. J Northwest A F Univ Nat Sci Ed, 2010,38(9):109-117. (in Chinese) |
| [28] | ( 刘虎岐, 刘应保, 宋云鹏, 等. 淡色库蚊氯菊酯抗性相关基因PR-OP全长cDNA的克隆及生物信息学分析[J]. 西北农林科技大学学报(自然科学版), 2010,38(9):109-117.) |
| [29] | Zhou D, Duan BY, Xu Y, et al. NYD-OP7/PLC regulatory signaling pathway regulates deltamethrin resistance in Culex pipiens pallens (Diptera ∶ Culicidae)[J]. Parasit Vectors, 2018,11(1):419. |
| [30] | Hu XB, Sun Y, Wang WJ, et al. Cloning and characterization of NYD-OP7, a novel deltamethrin resistance associated gene from Culex pipiens pallens[J]. Pestic Biochem Physiol, 2007,88(1):82-91. |
| [31] | Tang AH, Tu CP. Biochemical characterization of Drosophila glutathione S-transferases D1 and D21[J]. J Biol Chem, 1994,269(45):27876-27884. |
| [32] | Gong YH, Li T, Zhang L, et al. Permethrin induction of multiple cytochrome P450 genes in insecticide resistant mosquitoes, Culex quinquefasciatus[J]. Int J Biol Sci, 2013,9(9):863-871. |
| [33] | Tene BF, Poupardin R, Costantini C, et al. Resistance to DDT in an urban setting: common mechanisms implicated in both M and S forms of Anopheles gambiae in the City of Yaoundé Cameroon[J]. PLoS One, 2013,8(4):e61408. |
| [34] | Menze BD, Riveron JM, Ibrahim SS, et al. Multiple insecticide resistance in the malaria vector Anopheles funestus from northern Cameroon is mediated by metabolic resistance alongside potential target site insensitivity mutations[J]. PLoS One, 2016,11(10):e0163261 |
| [35] | Djègbè I, Agossa FR, Jones CM, et al. Molecular characterization of DDT resistance in Anopheles gambiae from Benin[J]. Parasit Vectors, 2014,7(1):409. |
| [36] | Paramasivan R, Sivaperumal R, Dhananjeyan KJ, et al. Prediction of 3-dimensional structure of salivary odorant-binding protein-2 of the mosquito Culex pipiens quinquefasciatus, the vector of human lymphatic filariasis[J]. In Silico Biol, 2007,7(1):1-6. |
| [37] | Xu PX, Choo YM, de la Rosa A, et al. Mosquito odorant receptor for DEET and methyl jasmonate[J]. Proc Natl Acad Sci USA, 2014,111(46):16592-16597. |
| [38] | Hallem EA, Dahanukar A, Carlson JR. Insect odor and taste receptors[J]. Annu Rev Entomol, 2006,51:113-135. |
| [39] | Ingham VA, Anthousi A, Douris V, et al. A sensory appendage protein protects malaria vectors from pyrethroids[J]. Nature, 2020,577(7790):376-380. |
| [40] | Yu TY, Garcia VE, Symington LS. CDK and Mec1/Tel1-catalyzed phosphorylation of Sae2 regulate different responses to DNA damage[J]. Nucleic Acids Res, 2019,47(21):11238-11249. |
/
| 〈 |
|
〉 |