中国农业科学 ›› 2016, Vol. 49 ›› Issue (13): 2534-2543.doi: 10.3864/j.issn.0578-1752.2016.13.009

所属专题: 烟粉虱

• 植物保护 • 上一篇    下一篇

烟粉虱溶菌酶基因的鉴定及表达分析

于 洁1,王登杰2,雷仲仁1,王海鸿1

 
  

  1. 1中国农业科学院植物保护研究所植物病虫害生物学国家重点实验室,北京 100193
    2四川省达州市达川区植保植检站,四川达州 635000
  • 收稿日期:2016-02-04 出版日期:2016-07-01 发布日期:2016-07-01
  • 通讯作者: 王海鸿,E-mail:wanghaihong2020@sina.com
  • 作者简介:于洁,Tel:010-62815930;E-mail:candyyudi@163.com。王登杰,E-mail:wdj198974@163.com。于洁和王登杰为同等贡献作者。
  • 基金资助:
    国家科技支撑计划(2012BAD19B06)、国家公益性行业(农业)科研专项(201303019-02)

Identification and Expression Analysis on Lysozyme Gene of Bemisia tabaci

YU Jie1, WANG Deng-jie2, LEI Zhong-ren1, WANG Hai-hong1   

  1. 1State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193
    2Dachuan Plant Protection and Quarantine Station in Dazhou City of Sichuan Province, Dazhou 635000, Sichuan
  • Received:2016-02-04 Online:2016-07-01 Published:2016-07-01

摘要: 【目的】鉴定烟粉虱(Bemisia tabaci)体内的溶菌酶基因,并对其序列特征、进化关系和表达模式进行分析,探索该基因在烟粉虱应对球孢白僵菌侵染过程中的作用,为进一步解析烟粉虱先天免疫过程提供理论依据。【方法】以第2代高通量测序的结果为依据,比对非冗余核苷酸数据库,筛选E值<10-5的基因为烟粉虱溶菌酶基因(lysozyme gene of Bemisia tabaci,BtLyz),利用ORF Finder预测BtLyz的开放阅读框,并得到氨基酸序列;利用Pfam和SMART预测BtLyz的蛋白质结构域;利用SinglP 4.1 Server预测BtLyz序列的信号肽;利用MEGA6.0软件进行BtLyz的氨基酸多序列比对;利用MEGA6.0软件,对BtLyz与其他昆虫的31个同源蛋白序列进行系统进化分析,并利用邻接法构建系统进化树来进一步分析鉴定该基因;采用实时荧光定量PCR(qRT-PCR)分析球孢白僵菌侵染烟粉虱卵、若虫、成虫后,在12、24、36、48、60 h时目标基因的时空表达模式。【结果】鉴定出烟粉虱4个溶菌酶基因,分别命名为BtLyz1、BtLyz2、BtLyz3和BtLyz4,基因序列全长分别为1 819、1 149、829和928 nt,分别编码159、160、148和160个氨基酸,且均含有信号肽。蛋白质结构分析、氨基酸序列比对和系统进化分析发现,BtLyz1和BtLyz4属于i型溶菌酶,BtLyz2和BtLyz3属于c型溶菌酶。BtLyz-1与豌豆长管蚜ApLyz-i1位于同一进化支上,BtLyz4与柑橘木虱DcLyz-i3和豌豆长管蚜ApLyz-i2位于同一进化支上。BtLyz2和BtLyz3与褐飞虱NlLyz-c1和异色瓢虫HaLyz-c3位于同一进化支上。与对照相比,卵期4个溶菌酶基因和若虫期BtLyz4的相对表达量没有显著变化;若虫期BtLyz1的相对表达量先上调再下调,而后又有所上调的波动趋势,24 h时上调表达最明显,为对照组的4.55倍;成虫期BtLyz1和BtLyz4的相对表达量均为先上调再下调,而后又有所上调的波动趋势,60 h时上调表达最明显,为对照组的11.31和4.21倍;若虫期BtLyz2和BtLyz3的相对表达量均为先上调再下调表达,在诱导24 h时上调表达最明显,为对照组的5.09和8.31倍,而后急剧下调,在60 h时下调表达最明显,为对照组的0.19和0.13倍;成虫期BtLyz2和BtLyz3的相对表达量均为先上调再下调的表达趋势,在24 h时上调最明显,为对照组的5.56和8.84倍。【结论】从烟粉虱体内鉴定了4个溶菌酶基因序列,其中2个为c型溶菌酶序列,2个为i型溶菌酶序列;在球孢白僵菌侵染烟粉虱不同虫态、不同时间后,卵期的相对表达量与对照相比没有显著变化;在若虫和成虫期,不同BtLyzs表现出不同的表达趋势。4个BtLyzs可能参与了烟粉虱对真菌侵染的免疫反应,有可能成为烟粉虱生物防治的潜在新靶标。

关键词: 烟粉虱, 溶菌酶, 生物信息学分析, 实时荧光定量PCR

Abstract: 【Objective】The objective of this study is to provide a basis for understanding the important role of lysozyme genes in Bemisia tabaci infected by the fungus Beauveria bassiana , identify lysozyme genes and analyze their sequence features, evolutionary relationships and expression pattern, and to provide a theoretical foundation for charifying innate immunity in B. tabaci. 【Method】Target genes were screened from the results of the second-generation high-throughput sequencing against the redundancy nucleotide databases when the E value <10-5. The open reading frame and amino acid sequence of BtLyzs were found using the ORF Finder software; the protein domains of BtLyzs were predicted using Pfam and SMART software; the signal peptide sequence of BtLyzs was predicted using SignalP 4.1 Server; the amino acid sequence alignment of BtLyzs was performed using the MEGA6.0 software; BtLyzs were characterized by using phylogenetic analysis with homologous genes of 31 other insects, a neighbor-joining of phylogenetic tree was constructed using the MEGA 6.0 software for further analysis and identification of the genes. The expression pattern of screened genes were determined at egg, nymph and adult stages of B. tabaci which were infected by the fungus Be. bassiana after 12, 24, 36, 48, 60 h using the quantitative real-time PCR (qRT-PCR). 【Result】 Four lysozyme genes were identified and designated as BtLyz1, BtLyz2, BtLyz3 and BtLyz4, which were 1 819, 1 149, 829, 928 nt, respectively. They were predicted to encode proteins of 159, 160, 148, 160 amino acids, respectively. Amino acid sequence alignment, phylogenetic analysis and protein tertiary structure prediction showed that BtLyz1 and BtLyz4 belong to i-type lysozymes, BtLyz2 and BtLyz3 belong to c-type lysozymes. BtLyz-1 formed a clade specific with Acyrthosiphon pisum of ApLyz-i1, BtLyz4 formed a clade specific with Diaphorina citri of DcLyz-i3 and Acyrthosiphon pisum of ApLyz-i2; BtLyz2 and BtLyz3 formed a clade specific with Nilaparvata lugens of NlLyz-c1 and Harmonia axyridis of HaLyz-c3. Compared to the control, the relative expression of all four genes in egg and BtLyz4 in nymph were not induced significantly by fungi-infected; the transcription of BtLyz1 nymph stage underwent the fluctuation of up-regulation, down-regulation, and second up-regulation and peaked at 24 h, it was increased 4.55 folds compared to the control; the transcription of BtLyz1 and BtLyz4 adult stage underwent the fluctuation of up-regulation, down-regulation, and second up-regulation and peaked at 60 h, they were increased 11.31 and 4.21 folds compared to the control. The transcription of BtLyz2 and BtLyz3 nymph stage underwent the fluctuation of up-regulation and down-regulation and peaked at 24 h, they were increased 5.09 and 8.31 folds compared to the control, then down-regulated obviously at 60 h, they were reduced 0.19 and 0.13 folds compared to the control. The transcription of BtLyz2 and BtLyz3 adult stage underwent the fluctuation of up-regulation and down-regulation and peaked at 24 h, they were increased 5.56 and 8.84 folds compared to the control. 【Conclusion】 Four BtLyz genes were identified in B. tabaci. Among them, two are c-type and two are i-type lysozymes. Transcriptional levels of four BtLyzs genes in B. tabaci were induced through different developmental stages and different time points in fungi-treated individuals compared to the control, they were not induced significantly in eggs, and showed different expression trends in nymphs and adults. The four BtLyzs genes probably might participated in the innate immune responses to fungus infection, and could be a new potential target for biocontrol of B. tabaci.

Key words: Bemisia tabaci, lysozymes, bioinformatic analyses, quantitative real-time PCR