昆虫学报 ›› 2017, Vol. 60 ›› Issue (1): 84-94.doi: 10.16380/j.kcxb.2017.01.010

• 研究论文 • 上一篇    下一篇

烟蚜线粒体基因组全序列及系统发育分析(英文)

杨洪1,2, 汪佳佳2, 戴仁怀2,*, 喻会平3, 杨茂发1,2   

  1. (1. 贵州大学烟草学院, 贵阳 550025; 2. 贵州大学昆虫研究所, 贵阳 550025; 3. 贵州省烟草公司毕节市公司, 贵州毕节 551700)
  • 出版日期:2017-01-20 发布日期:2017-01-20

Sequence and phylogenetic analysis of the complete mitogenome of Myzus persicae (Hemiptera: Aphididae)(In English

YANG Hong1,2, WANG Jia-Jia2, DAI Ren-Huai2,*, YU Hui-Ping3, YANG Mao-Fa1,2   

  1. (1. College of Tobacco Science, Guizhou University, Guiyang 551700, China; 2. Entomology Institute, Guizhou University, Guiyang 550025, China; 3. Bijie Branch of Guizhou Tobacco Corporation, Bijie, Guizhou 551700, China)
  • Online:2017-01-20 Published:2017-01-20

摘要: 【目的】蚜虫是杂食性农业害虫。本研究旨在通过线粒体基因组分析更好地了解蚜科昆虫的系统发育关系。【方法】结合第二代测序和PCR扩增技术获得了烟蚜Myzus persicae线粒体基因组全序列,与蚜科其他昆虫进行了对比分析;以贝叶斯法和最大似然法基于13个蛋白编码基因对蚜科进行了系统发育分析。【结果】烟蚜线粒体基因组(GenBank登录号: KU_236024)序列全长17 832 bp,A+T含量84.1%,AT偏斜为0.094,GC偏斜为-0.296。包含13个蛋白编码基因(cox1-3, nad1-6, nad4L, atp6, atp8和cytb),22个tRNA,2个rRNA基因(rrnLrrnS)和2个长的非编码区,其基因排列顺序与已知的蚜科昆虫相似,除了nad4以单独的T结尾,所有的蛋白编码基因均以ATN作为起始密码子,TAA作为终止密码子。在烟蚜线粒体基因组中,tRNAGlu 和 tRNAPhe中间有一段307 bp的非编码区,该编码区包含2个重复单元,烟蚜的控制区长2 531 bp,是所有测序蚜虫线粒体基因组中最长的。rrnL的二级结构包含6个结构域,44个茎环结构;rrnS的二级结构有3个结构域,24个茎环结构。基于烟粉虱和其他20种昆虫的13个蛋白编码基因重建的BI和ML系统发育树,与传统形态学分类结果一致。【结论】蚜亚科和长管蚜亚科的单系性得到了很好的支持;在长管蚜亚科的分支中,M. persicaeD. noxia聚成一支,并且C. salicicola位于进化枝的底部。本研究结果为蚜科系统发生关系重建积累了有价值的数据资料。

关键词:  , 蚜虫; 烟蚜; 线粒体基因组; RNA; 二级结构; 系统发育

Abstract: 【Aim】 Aphididae insects are extremely polyphagous agricultural pests. The present study aims to better understand the phylogenetic relationship within Aphididae by mitogenome analysis. 【Methods】 The complete mitogenome of Myzus persicae (Hemiptera: Aphididae) was sequenced by next-generation sequencing (NGS) and PCR amplification, and sequence analysis in comparison with other Aphididae insects was made. The Bayesian and Maximum Likelihood phylogenetic analyses based on the 13 PCG sequences were performed. 【Results】 The mitogenome of M. persicae (GenBank accession no. KU_236024) is 17 832 bp in size and harbors an A+T content of 84.1%, an AT-skew of 0.094 and a GC-skew of -0.296. It contains 13 protein-coding genes (PCGs) (cox1-3, nad1-6, nad4L, atp6, atp8 and cytb), 22 tRNA genes, two rRNA genes (rrnL and rrnS) and two long non-coding regions. The gene order and orientation of the mitogenome of M.persicae are similar to those of other aphids. All protein-coding genes (PCGs) start with an ATN initiation codon and terminate with a TAA stop codon, except for nad4, which terminates with a single T nucleotide as an incomplete stop codon. In the mitogenome of M.persicae, there is a non-coding region (307 bp) between tRNAGlu and tRNAPhe, which includes variable numbers of tandem repeats in a lineage-specific manner. The mitogenome of M.persicae has the control region of 2 531 bp which is the longest in all the sequenced aphid mitogenomes. In the mitogenome of M. persicae, six domains with 44 helices are present in the secondary structure of rrnL and three domains with 24 helices are present in the secondary structure of rrnS. Phylogenetic analysis based on 13 PCGs from M. persicae and 20 other insect species showed that every superfamily in the tree forms a monophyletic clade. 【Conclusion】 Within the Aphidinae subfamily, the monophyly of Aphidini and Macrosiphini had statistically high values. M.p ersicae was clustered with D. noxia, and C.s alicicola was positioned most basally within the clade of Macrosiphini. The results of this study provide some valuable molecular data for reconstructing the phylogenetic relationship in Aphididae.

Key words: Aphid, Myzus persicae, mitogenome, RNA, secondary structure, phylogeny