昆虫学报 ›› 2023, Vol. 66 ›› Issue (7): 946-958.doi: 10.16380/j.kcxb.2023.07.009

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

靛灰蝶和彩斑尾蚬蝶线粒体基因组及基于线粒体基因的凤蝶总科系统发育关系

闫振天1, 罗斯特2, 范真槐1, 陈斌1,*   

  1. (1. 重庆师范大学昆虫与分子生物学研究所, 媒介昆虫重庆市重点实验室, 重庆 401331; 2. 厦门大学生命科学学院, 厦门 361102)
  • 出版日期:2023-07-20 发布日期:2023-08-17

Mitochondrial genomes of Caerulea coeligena (Lepidoptera: Lycaenidae) and Dodona maculosa (Lepidoptera: Riodinidae) and phylogenetic relationships in Papilionoidea based on mitochondrial genes

YAN Zhen-Tian1, LUO Si-Te2, FAN Zhen-Huai1, CHEN Bin1,*   

  1.  (1. Chongqing Key Laboratory of Vector Insects, Institute of Entomology and Molecular Biology, Chongqing Normal University, Chongqing 401331, China; 2. School of Life Sciences, Xiamen University, Xiamen 361102, China)
  • Online:2023-07-20 Published:2023-08-17

摘要: 【目的】对靛灰蝶Caerulea coeligena和彩斑尾蚬蝶Dodona maculosa线粒体全基因组进行测序和分析,并基于已知线粒体基因组数据探析凤蝶总科(Papilionoidea)的系统发育关系。【方法】利用Illumina NovaSeq 6000测序平台测定靛灰蝶和彩斑尾蚬蝶线粒体全基因组序列,对其进行拼装和注释,并对其tRNA基因二级结构进行预测;基于NCBI数据库中已报道的凤蝶总科6科48种蝴蝶线粒体基因组13个蛋白编码基因(protein-coding genes, PCGs)核苷酸序列,对凤蝶总科50个种的13个PCGs的选择压力进行分析;利用最大似然法(maximum likelihood, ML)和贝叶斯法(Bayesian inference, BI)分析这50个种及其所属科和亚科的系统发育关系。【结果】靛灰蝶和彩斑尾蚬蝶线粒体全基因组长分别是15 164和15 486 bp (GenBank登录号分别为MZ489120和MZ489121),其基因组大小、结构和组成、基因排列、核苷酸组成、密码子使用和tRNA结构预测等均与已知近缘种蝴蝶相似;凤蝶总科线粒体基因组进化经历了高水平的纯化选择,13个PCGs非同义替换率和同义替换率比值(Ka/Ks)明显小于1。凤蝶总科重建的系统发育关系为:(凤蝶科(Papilionidae)+弄蝶科(Hesperiidae)+粉蝶科(Pieridae)+(蛱蝶科(Nymphalidae)+(灰蝶科(Lycaenidae)+蚬蝶科(Riodinidae))))。【结论】基于凤蝶总科已报道的48种和本研究新测序两种线粒体基因组重建的系统发育树表明蚬蝶科与灰蝶科为单系且互为姐妹群关系,且与蛱蝶科具有更近的亲缘关系,凤蝶科可能是凤蝶总科中最早进化的类群。本研究为进一步研究蝴蝶的系统进化和完善其分类系统提供了新的分子数据支撑。

关键词: 靛灰蝶, 彩斑尾蚬蝶, 凤蝶总科, 线粒体基因组, 系统发育

Abstract: 【Aim】 To sequence and analyze the complete mitochondrial genomes (mtgenomes) of Caerulea coeligena and Dodona maculosa, and to explore the phylogenetic relationships in Papilionoidea based on the known mitochondrial genome data. 【Methods】 The complete mtgenomes of C. coeligena and D. maculosa were sequenced using Illumina NovaSeq 6000 sequencing platform, assembled and annotated, and the secondary structure of tRNA genes was predicted. Based on the nucleotide sequences of 13 protein-coding genes (PCGs) from the mtgenomes of 48 reported butterfly species belonging to six families in Papilionoidea obtained from the NCBI database, the selective pressures on these 13 PCGs in 50 species were analyzed. The phylogenetic relationships among these 50 species, as well as their respective families and subfamilies, were analyzed using maximum likelihood (ML) and Bayesian inference (BI) methods. 【Results】 The mtgenomes of C. coeligena and D. maculosa (GenBank accession numbers are MZ489120 and MZ 489121, respectively) are 15 164 and 15 486 bp in length, respectively. Their genome size, structure and composition, gene arrangement, nucleotide composition, codon usage, and predicted tRNA structure of mitogenomes of C. coeligena and D. maculosa were similar to those of the known closely related butterfly species. The mtgenomes in Papilionoidea exhibited a high level of purifying selection, with the ratios of nonsynonymous substitution rate to synonymous substitution rate (Ka/Ks) for the 13 PCGs obviously lower than 1. The phylogenetic relationships in Papilionoidea were (Papilionidae+Hesperiidae+Pieridae+(Nymphalidae+ (Lycaenidae+Riodinidae))).【Conclusion】 The reconstructed phylogenetic tree based on the mtgenomes of 48 reported butterfly species in Papilionoidea and two newly sequenced species in this study suggests that Riodinidae and Lycaenidae be monophyletic, and sister each other. The sister group would be closely related to Nymphalidae in phylogenetics. Moreover, Papilionidae may represent one of the earliest lineages in Papilionoidea. This study provides new molecular data to support further research on the systematic evolution and refinement of the classification system for butterflies.

Key words: Caerulea coeligena, Dodona maculosa, Papilionoidea, mitochondrial genome, phylogeny