Abstract：Agrilus zanthoxylumi is a significant and destructive branch-borers pest in the Zanthoxylum bungeanum planting areas of northern China. This insect mainly recognizes volatile odor molecules through sensitivity and specificity via antenna, which plays an important role in its survival and reproduction. To select Agrilus zanthoxylumi adult male and female antenna of differentially expressed genes, preliminary parsing and Agrilus zanthoxylumi candidate differences related to the behavior and the biological prevention and control in terms of seeking mate genes, this study selected adult female and male antenna transcriptome sequencing as the research object, using bioinformatics software to analyze the transcriptome sequencing data obtained and screening of differentially expressed genes between the 2 samples of male and female, further analysis of differentially expressed genes. The results showed that 36 209 unigene were obtained from the antenna transcriptome data of the male and female adults of Agrilus zanthoxylumi. A total of 726 genes were detected with significant differences in expression, among which 459 genes were up-expressed and 267 genes were down-expressed. Through GO (Gene Ontology) enrichment analysis, the differentially expressed genes were found in binding function, catalytic activity, cell process, metabolic process, cell components and organelles which were significantly enriched. A total of 472 differentially expressed genes were annotated in the KEGG database, and Pathway enrichment analysis showed that they were involved in 165 metabolic pathways, with 394 up-expressed and 78 down-expressed genes. This study obtained Agrilus zanthoxylumi adult antenna of the genetic basis of enrichment of the active ingredient based on the transcriptome sequencing technologies, discussed the adult male and female antenna differentially expressed genes, revealed the gene function, understood the relationship between genes. It doesn't only further improve the mechanism of tentacles feelings outside odor material plays an important role, but also provide the basic data for pest prevention and control technology.
巩雪芳, 谢寿安, 杨平, 胡洪焰, 郭莉, 车显荣, 陈迪, 王延来, 吕淑杰. 花椒窄吉丁雌、雄成虫触角差异表达基因的初步筛选[J]. 农业生物技术学报, 2019, 27(7): 1233-1245.
GONG Xue-Fang, XIE Shou-An, YANG Ping, HU Hong-Yan, GUO Li, CHE Xian-Rong, CHEN Di, WANG Yan-Lai, LV Shu-Jie. Preliminary Screening of Differentially Expressed Genes of Female and Male Antenna of Agrilus zanthoxylumi. 农业生物技术学报, 2019, 27(7): 1233-1245.
[1] 安涛, 张洪志, 韩艳华, 等. 2017. 烟蚜茧蜂滞育关联基因的转录组学分析[J]. 中国生物防治学报, 33(5): 604-611. (An T, Zhang H Z, Han Y H, et al.2017. Transcriptome analysis of diapause related genes in aphidiophora tabaci[J]. Chinese Journal of Biological Control, 33(5): 604-611.) [2] 车显荣. 2018. 花椒园昆虫种群群落结构及田间诱捕效果研究[D]. 硕士学位论文, 西北农林科技大学, 导师: 谢寿安, pp.13-28. (Che X R.Study on insect population and community structure and field trapping effect in Zanthoxylum bungeanum[D]. Thesis for M.S., Northwest A&F university, Supervisor: Xie S A, pp. 13-28.) [3] 仇占南, 付艳, 杜传慧, 等. 2018. 苹果柱型基因候选基因的筛选[J]. 农业生物技术学报, 26(1): 53-63. (Qiu Z N, Fu Y, Du C H, et al.2018. Selection of candidate genes for apple stem type genes[J]. Journal of Agricultural Biotechnology, 26(1): 53-63) [4] 崔元, 李晓轩, 孙继国, 等. 2018. 基于RNA-Seq筛选新型鸭细小病毒感染鸭胚成纤维细胞的差异表达基因[J]. 现代畜牧兽医, 354(05): 6-11. (Cui Y, Li X X, Sun J G, et al.2018. Screening the differentially expressed genes of duck embryo fibroblasts infected by novel parvo based on rna-seq[J]. Modern Animal Husbandry Veterinarian, 354(05): 6-11. [5] 樊莹, 董霞. 2015. 抑制消减杂交技术在昆虫差异表达基因研究中的应用[J]. 贵州农业科学, 43(8): 122-125. (FanY, Dong X.2015. Application of subtraction hybridization in insect differential expression gene research[J]. Guizhou Agricultural Science, 43(8): 122-125.) [6] 高焕婷, 张国龙. 2007. 花椒窄吉丁虫的生物学特性及其综合防治措施[J]. 陕西农业科学, (2): 183-184. (Gao H T, Zhang G L. 2007. Biological characteristics and comprehensive control measures of A#. zanthoxylumi#[J]. Shaanxi Agricultural Science, (2): 183-184.) [7] 高莺, 李继华, 韩立赤, 等. 2009. 张应力诱导大鼠骨髓间充质干细胞骨向分化及其差异基因表达分析[J]. 华西口腔医学杂志, 27(2): 213-216. ( Gao Y, Li J H, Han L C, et al.2009. Analysis of bone differentiation and differential gene expression of bone marrow mesenchymal stem cells induced by zhang stress in rats[J]. Journal of Western China Dental Medicine, 27(2): 213-216.) [8] 谷少华. 2010. 苜蓿盲蝽嗅觉相关基因的发掘及功能分析[D]. 硕士学位论文, 中国农业科学院, 导师: 张永军. pp. 1-4. (Gu S H.2010. Discovery and functional analysis of olfactory related genes of Adelphocoris suturalis[D]. Thesis for M.S., Chinese Academy of Agricultural Sciences, Supervisor: Zhang Y J, pp. 1-4.) [9] 胡烨, 李影, 钱爱东. 2010. 几种基因差异表达分析技术在原核生物差异基因筛选中的应用[J]﹒动物医学进展, 31(2): 103-106. (Hu Y, LiY, Qian A D.2010. Application of several gene expression differential analysis techniques in screening prokaryotic differential genes[J]. Animal Medical Advances, 31(2): 103-106.) [10] 姜峰, 王磊, 严会超, 等. 2006. 雌、雄家蚕血液中的保护酶差异研究[J]. 广东农业科学, 8(4): 74-76. ( Jiang F, Wang L, Yan H C, et al.2006. Differential studies on protective enzymes in the blood of male and female silkworms[J]. Guangdong Agricultural Science, 8(4): 74-76.) [11] 姜洪真, 马伯军, 钱前, 等. 2018. 全基因组关联分析(GWAS)在作物农艺性状研究中的应用[J]. 农业生物技术学报, 26(7): 1244-1257. (Jiang H Z, Ma B J, Qian J, et al.2018. Application of genome-wide association analysis (GWAS) in agronomic traits of crops[J]. Journal of Agricultural Biotechnology, 26(7): 1244-1257.) [12] 李佼旬, 邱林, 王选琴, 等. 2018. 基因芯片筛选非综合征型唇腭裂差异表达基因的初步研究[J]. 重庆医科大学学报, 43(8): 1102-1108. (Li J X, Qiu L, Wang X Q, et al.2018. A preliminary study on gene chip selection of differential expression genes for non-syndromic cleft lip and palate[J]. Journal of Chongqing Medical University, 43(8): 1102-1108.) [13] 李欢, 杨新杰, 王维戚, 等. 2018. 应用RNA-Seq技术筛选唾液腺腺样囊性癌嗜神经侵袭相关差异表达基因[J]. 中国口腔颌面外科杂志, 16(2): 114-119. (Li H, Yang X J, Wang W Q, et al.2018. Application of rna-seq technology to screen the differentially expressed genes related to neurotrophic invasiveness in adenoid cystic carcinoma of salivary gland[J]. Chinese Journal of Oral and Maxillofacial Surgery, 16(2): 114-119.) [14] 刘海涛, 李兵, 赵国栋, 等. 2010. 家蚕幼虫不同龄期体内主要解毒酶及其基因表达的性别差异[J]. 昆虫学报, 53(5): 479-486. (Liu H T, Li B, Zhao G D,et al.2010. Sexual differences in main detoxification enzymes and their gene expression in different instars of Bombyx morilarvae[J]. Acta Entomlogica Sinica, 53(5): 479-486.) [15] 刘芳, 宗超, 余林生. 2015. 意大利蜜蜂哺育蜂与采集蜂行为转变相关基因的表达差异研究[J]. 应用昆虫学报, 52(2): 300-307. (Liu F, Zong C, Yu L S.2015. Study on differences in the expression of genes related to changes in bee feeding and collecting behavior in Apis mellifera ligustica Spin[J]. Journal of Applied Insect, 52(2): 300-307.) [16] 刘斌. 2012. 柑橘全抓螨不同品系转录组及基因表达谱差异分析[D]. 硕士学位论文, 西南大学, 导师: 冉春. pp. 8-14. (Liu B.2012. Analysis on the differences of transcriptome and gene expression profile between different lines of Panonchus citri Mc Gregor[D]. Thesis for M.S., Southwest university, Supervisor: Ran C, pp. 8-14.) [17] 陆雪莹. 2016. 准噶尔小胸鳖甲短时低温胁迫响应的转录组分析[J]. 昆虫学报, 59(6): 581-591. (Transcriptome analysis of short-time response to low temperature stress in Microdera punctipennis dzungarica Kasz[J]. Journal of Insect, 59(6): 581-591.) [18] 龙桂友, 饶力群. 2004. 基因差异表达分析技术进展[J]. 生命科学研究, (s2): 48-52. (Long G Y, Rao L Q. 2004. Advances in gene expression analysis techniques[J]. Life science research, (s2): 48-52.) [19] 莫之婧, 马义丽, 李康智, 等. 2016. 以LPAR3为中心的2型糖尿病相关基因功能富集及蛋白互作分析[J]. 基因组学与应用生物学, (11): 34-40. (Mo Z J, Ma Y L, Li K Z, et al. 2016. Functional enrichment and protein interaction analysis of related genes in type 2 diabetes mellitus centered on LPAR3[J]. Genomics & Applied Biology, (11): 34-40) [20] 权丽君. 2016. 花椒窄吉丁的发生及其防治措施[J]. 陕西农业科学, 62(12): 116-116. (Quan L J.2016. Occurrence and control measures of A. zanthoxylumi[J]. Shaanxi Agricultural Science, 62(12): 116-116) [21] 任守杰, 张志娥, 陈晓玲, 等. 2006. 基于醇溶蛋白的20份小麦种质遗传完整性分析[J]. 植物遗传资源学报, 7(1): 44-48. (Ren S J, Zhang Z E, Chen X L, et al.2006. Genetic integrity analysis of 20 wheat germplasm based on alcohol soluble protein[J]. Journal of Plant Genetic Resources, 7(1): 44-48.) [22] 孙立军, 陈忠平, 黄强. 2002. 基因差异表达的高通量分析及其在肿瘤研究中的应用[J]. 癌症, 21(5): 571-575. (Sun L J, Chen Z P, Huang Q.2002. High-throughput analysis of gene differential expression and its application in tumor research[J]. Cancer, 21(5): 571-575.) [23] 宋志强, 郝飞. 2001. 差异表达基因分析[J]. 医学分子生物学杂志, 23(4): 234-237. (Song Z Q, Hao F.2001. Analysis of differentially expressed genes[J]. Journal of Medical Molecular Biology, 23(4): 234-237.) [24] 王桂荣, 郭予元, 吴孔明. 2002a. 棉铃虫普通气味结合蛋白Ⅱ基因的表达及鉴定[J]. 昆虫学报, 45(3): 285-289. (Wang G R, Guo Y Y, Wu K M.2002. Helicoverpa armigera general Ⅱ odorant binding protein gene expression and identification[J]. Journal of Insects, 45(3): 285-289.) [25] 王桂荣, 郭予元, 吴孔明. 2002b. 昆虫触角气味结合蛋白的研究进展[J]. 昆虫学报, 45(1): 131-137. (Wang G R, Guo Y Y, Wu K M.2002. Research progress of insect antenna odor binding protein[J]. Journal of Entomology, 45(1): 131-137.) [26] 吴运梅, 张丽红, 王国秀. 2012. 中华卵索线虫雌、雄寄生后期幼虫基因差异表达分析[J]. 动物学研究, 33(5): 487-492. (Wu Y M, Zhang L H, Wang G X.2012. Analysis of the gene differential expression of the later larvae of Ovomermis sinensis[J]. Zoological Studies, 33(5): 487-492.) [27] 肖海军, 邹超, 陈俊晖, 等. 2014. 昆虫滞育关联差异基因的筛选方法[J]. 江西农业大学学报, 36( 6): 1240-1246. (Xiao H J, Zou C, Chen J H, et al.2014. Screening methods of diapause related genes in insects[J]. Journal of Jiangxi Agricultural University, 36(6): 1240-1246.) [28] 杨婧. 2013. 短额负蝗三种虫态的比较转录组及线粒体转录组作图研究[D]. 硕士学位论文, 陕西师范大学, 导师: 黄原, pp.13-28. (Yang J.2013. A comparative study on the three insect states of the Atractomorpha sinensis Bolivar: Transcriptome and mitochondrial transcriptome mapping[D]. Thesis for M.S., Shaanxi Normal University, Superviser: Huang Y, pp.13-28.) [29] 赵林川, 姚祥, 沈卫德, 等. 1999. 家蚕过氧化氢代谢的性别差异[J]. 蚕业科学, 25(2): 127-128. ( Zhao L C, Yao X, Shen W D, et al.1999. Gender differences in hydrogen peroxide metabolism of silkworm[J]. Sericulture Science, 25(2): 127-128.) [30] 张春兰, 秦孜娟, 王桂芝, 等. 2012. 转录组与RNA-Seq技术[J]. 生物技术通报, 12: 51-56. (Zhang C L, Qin Z J, Wang G Z, et al.2012. Transcriptome and rna-seq technology[J]. Bulletin of Biotechnology, 12: 51-56.) [31] 周作勇, 邵骜骏, 徐婷婷, 等. 2016. 基于高通量测序的猪鞭虫雌、雄虫差异基因表达分析[J]. 畜牧兽医学报, 47(9): 1914-1923. (Zhou Z Y, Shao A J, Xu T T, et al.2016. Differential gene expression analysis of female and male Trichuris suis based on high-throughput sequencing[J]. Journal of Animal Husbandry and Veterinary Medicine, 47(9): 1914-1923.) [32] 朱海生, 刘建汀, 王彬, 等. 2018. 丝瓜过氧化物酶基因的克隆及其在果实褐变中的表达分析[J]. 农业生物技术学报, 26(1): 43-52. (Zhu H S, Liu J T, Wang B, et al.2018. Cloning and expression of peroxidase gene in fruit Browning of towel gourd[J]. Journal of Agricultural Biotechnology, 26(1): 43-52.) [33] 宗超, 刘芳, 余林生. 2014. 中华蜜蜂哺育蜂与采集蜂行为转变相关基因的表达差异研究[J]. 应用昆虫学报, 51(2): 440-447. (Zong C, Liu F, Yu L S.2014. Differences in the expression of genes related to changes in bee feeding and collecting behavior of Apis cerana[J]. Journal of Applied Insect, 51(2): 440-447.) [34] Burkett-Cadena N D, Mullen G R, 2007. Field comparison of Bermuda-hay infusion to infusions of emergent aquatic vegetation for collecting female mosquitoes[J]. Journal of the American Mosquito Control Association, 23(2): 117-123. [35] Cilek J E, Ikediobi C O, Hallmon C F, et al.2011. Semi-field evaluation of several novel alkenol analogs of 1-octen-3-ol as attractants to adult Aedes albopictus and Culex quinquefasciatus[J]. Journal of the American Mosquito Control Association, 27(3): 256-262. [36] Etebari K, Palfreyman R W, Schlipalius D, et al.2011. Deep sequencing-based transcriptome analysis of Plutella xylostella larvae parasitized by Diadegma semiclausum[J]. BMC Genemics, 12(1): 446-463. [37] Howard B E, Heber S.2010. Towards reliable isoform quantification using RNA-SEQ data[J]. BMC Bioinformatics, 11(Suppl3): 6. [38] Hao Y J, Li W S, He Z B, et al.2012. Differential gene expression between summer and winter diapause pupae of the onion maggot Delia antiqua, detected by suppressive subtractive hybridization[J]. Journal of Insect Physiology, 58(11): 1444-1449. [39] Hull J J, Geib S M, Fabrick J A, et al.2013. Sequencing and de novo assembly of the Western tarnished plant bug Lygushesperus Transcriptome[J]. PLoS ONE, 8(1): e55105. [40] Haas B, Zody M.2010. Advancing RNA-Seq analysis[J]. Nature Biotechnology, 28(5): 421-423. [41] Westfall P H, Young S S.1989. P Value adjustments for multiple tests in multivariate binomial models[J]. Journal of the American Statistical Association, 84: 780-786. [42] Hoddle M S, Millar J G, Hoddle C D, et al.2011. Field optimization of the sex pheromone of Stenoma catenifer (Lepidoptera: Elachistidae): Evaluation of lure types, trap height, male flight distances, and number of traps needed per avocado orchard for detection[J]. Bulletin of Entomological Research, 101(2): 145-152. [43] Leal W S.2005.Pheromone reception[J]. Topics in Current Chemistry, 240: 341-360. [44] Lin Z, He D, Zhang X, et al.2005. Linkage map construction and map-ping QTL for cotton fibre quality using SRAP, SSR and RAPD[J]. Plant Breeding, 124: 180-187. [45] Martinez-Barnetche J, Gomez-Barreto R, Ovilla-Muoz M, et al.2012. Transcriptome of the adult female malaria mosquito vector Anopheles albimanus[J]. BMC Genomics, 13(1): 207. [46] Mortazavi A, Williams B A, McCue K, 2008. Mapping and quantifying mammalian Transcriptomes by rna-seq[J]. Nat Methods, 5(7): 621-628. [47] Maxwell C, Antoshechkin I, Kurhanewicz N, et al.2012. Nutritional control of mRNA isoform expression during developmental arrest and recovery in C.elegans[J]. Genome Research, 22(10): 1920-1929. [48] Nookaew I, Papini M, Pornputtpong N, et al.2012. A comprehensive comparison of RNA-Seq-based transcriptome analysis fromreads to differential gene expression and cross-comparison with microarrays: A case study in Saccharomyces cerevisiae[J]. Nucleic Acids Research, 40(20): 10084-10097. [49] Poelchau M, Reynolds J, Denlinger D L, et al.2011. A de novo transcriptome of the Asian tiger mosquito, Aedesalbopictus, to identify candidate transcripts for diapause preparation[J]. BMC Genomics, 12(1): 619. [50] Zhang T, Feng J, Cai C, et al.2011. Synthesis and field test of three candidates for soybean pod borer's sex pheromone[J]. Natural Product Communications, 6(9): 1323-1326.
