昆虫学报 ›› 2019, Vol. 62 ›› Issue (1): 101-116.doi: 10.16380/j.kcxb.2019.01.011

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

林氏按蚊线粒体全基因组序列的测定及基于线粒体基因组的按蚊属系统发育分析(英文)

毛启萌, 李廷景, 付文博, 闫振天, 陈斌   

  1.  (重庆师范大学昆虫与分子生物学研究所, 重庆市媒介昆虫重点实验室, 重庆 401331)
  • 出版日期:2019-01-20 发布日期:2019-01-24

Sequencing of the complete mitochondrial genome of Anopheles lindesayi and a phylogenetic analysis of the genus Anopheles (Diptera: Culicidae) based on mitochondrial genomes (In English)

MAO Qi-Meng, LI Ting-Jing, FU Wen-Bo, YAN Zhen-Tian, CHEN Bin*   

  1. (Chongqing Key Laboratory of Vector Insects, Institute of Entomology and Molecular Biology, Chongqing Normal University, Chongqing 401331, China)
  • Online:2019-01-20 Published:2019-01-24

摘要:

【目的】对林氏按蚊Anopheles lindesayi完整的线粒体基因组进行测序及分析,依据已知的线粒体基因组构建并讨论按蚊属蚊虫的分子系统发育关系。【方法】对林氏按蚊线粒体基因组进行测序、注释,并对其基本特征和基本组成进行分析。基于串联的13个蛋白质编码基因的核苷酸序列和氨基酸序列,用ML法和贝叶斯法构建林氏按蚊和按蚊属其他32种蚊虫的系统发育树,据此探讨按蚊属蚊虫的系统发育关系和系统分类。【结果】林氏按蚊线粒体基因组全长为15 366 bp,包含13个蛋白质编码基因,22个tRNA基因,2个rRNA基因和一段控制区。林氏按蚊线粒体基因组呈现明显的AT偏斜和GC偏斜,AT偏斜为正,GC偏斜为负。除了COX1使用TCG和ND5使用GTG作为起始密码子以外,其他蛋白质编码基因的起始密码子均遵循ATN原则;终止密码子为TAA或者T。除了tRNASer(AGN)以外,其他的tRNA基因均呈现典型的三叶草二级结构。控制区AT含量最高,为94.54%。滑窗分析显示蛋白质编码基因是用于构建亚属或属水平系统发育关系的最佳分子标记。系统发育树强烈支持塞蚊亚属Cellia、按蚊亚属Anopheles、徕蚊亚属Nyssorhynchus和柯特蚊亚属Kerteszia均为单系群。小五斑按蚊An. atroparvus和四斑按蚊An. quadrimaculatus A这两个种聚到一起,从传统的形态分类上讲,它们和林氏按蚊均属于按蚊亚属按蚊系蚊虫。但本研究构建的4个系统发育树均显示,(小五斑按蚊An. atroparvus+四斑按蚊An. quadrimaculatus A)和林氏按蚊被属于迈蚊系的中华按蚊分开,这为两个系的分类提供了新的论点。【结论】本研究获得了林氏按蚊的完整的线粒体基因组,探析了按蚊属的线粒体基因组特征和系统发育关系,为进一步研究蚊科线粒体基因组和系统发育关系提供了依据。

关键词:  蚊科, 按蚊属, 林氏按蚊, 线粒体基因组, 核苷酸组成, 系统发育

Abstract:

【Aim】 To sequence and analyze the complete mitochondrial genome (mtgenome) of Anopheles lindesayi, and to construct and discuss the molecular phylogenetic relationships of the genus Anopheles based on the known mtgenomes. 【Methods】 The complete mtgenome of An. lindesayi was sequenced and annotated, and its general features and base composition were analyzed. The phylogenetic relationships between An. lindesayi and other 32 species in the genus Anopheles were reconstructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the concatenated nucleotide sequences and amino acid sequences of 13 protein-coding genes (PCGs), and the systematics of the genus Anopheles was discussed based on the phylogenetic analysis. 【Results】 The complete mtgenome sequence of An. lindesayi is 15 366 bp in length, which contains 13 PCGs, 22 tRNA genes, two rRNA genes, and one control region. The mtgenome has a clear bias in nucleotide composition with a positive AT-skew and a negative GC-skew. The initiation codons of PCGs comply with the ATN rule except that COX1 uses TCG and ND5 uses GTG as the start codon, and the termination codon is TAA or incomplete T. The tRNAs have the typical clover-leaf structure, but  tRNASer(AGN) has a large loop instead of the conserved stem-and-loop structure. The control region has the highest AT content 94.54%). The sliding window analysis suggested that the PCGs are the most suitable markers to elucidate the phylogenetic relationships at the subgenus and genus levels. The phylogenetic analysis strongly supported each of the subgenera Cellia, Anopheles, Nyssorhynchus and Kerteszia to be a monophyletic group. An. lindesayi and An. atroparvus+An. quadrimaculatus A, the latter two traditionally classified as Anopheles Series, are separated by An. sinensis in Myzorhynchus Series on all the four phylogenetic trees, raising a new debate on the taxonomy of the two Series. 【Conclusion】This study obtained the complete mtgenome sequence of An. lindesayi and explored the mtgenome characteristics and phylogeny of the genus Anopheles, establishing the information frame for further studying the mtgenomes and phylogeny of the family Culicidae.

Key words:  Culicidae, Anopheles, Anopheles lindesayi, mitochondrial genome, nucleotide composition, phylogeny