Atteva charopis的比较线粒体基因组研究及系统发育分析

doi:10.13560/j.cnki.biotech.bull.1985.2023-0171

摘要/Abstract

摘要：

线粒体基因组有助于对巢蛾总科的分子进化和系统发育的理解。本研究测序、分析了Atteva charopis Turner, 1903的线粒体全基因组，并基于线粒体基因组的核苷酸和氨基酸数据重建了巢蛾总科的系统发育关系。在系统发育分析中，内群取样包含巢蛾总科的10种昆虫，外群选择2个谷蛾总科Tineoidea（Phyllocnistis citrella和Caloptilia theivora）昆虫；分别利用最大似然法和贝叶斯法重建巢蛾总科的系统发育关系。Atteva charopis的线粒体基因组全长为15 163 bp（GenBank序列号：OQ130005），包含37个基因（13个蛋白质编码基因、2个核糖体RNA基因和22个转运RNA基因）和一段非编码控制区。与假定昆虫的trn基因（trnI-trnQ-trnM）相比，该线粒体基因出现了trn基因重排（trnI-trnM-trnQ）。两种不同的系统发育分析方法构建的系统发育关系树显示，潜叶蛾科、Praydidae和斑巢蛾科为单系群，而菜蛾科为非单系群，斑巢蛾科和Praydidae亲缘关系较近，Atteva charopis与Atteva aurea形成姐妹群。本研究结果可为深入研究巢蛾总科的系统发育提供线粒体基因组数据。

关键词: 鳞翅目, 巢蛾总科, 比较线粒体基因组, 系统发育

Abstract:

Mitochondrial genomes（mitogenomes）advance our understanding of molecular evolution and phylogenetic relationships in Yponomeutoidea. In this study, we sequenced and analyzed the mitochondrial genome of Atteva charopis Turner, 1903. We reconstructed the phylogenetic relationships of Yponomeutoidea based on the nucleotide and amino acid data of the mitochondrial genome. In phylogenetic analyses, we selected 10 exemplars of Yponomeutoidea as ingroups，and two species of Tineoidea（Phyllocnistis citrella and Caloptilia theivora）as outgroups. The phylogenetic analyses were conducted using Maximum likelihood and Bayesian inference methods. The mitochondrial genome of Atteva charopis was a circular molecule of 15 163 bp in length（GenBank accession No. OQ130005）, containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a non-coding control region. Compared to the putative insect trn gene（trnI-trnQ-trnM）, there has been a trn gene rearrangement（trnI-trnM-trnQ）in this mitochondrial genome. Both phylogenetic inference methods produced a similar tree topological structure:Family Lyonetiidae, Praydidae and Attevidae were the monophyletic groups. However, the family Plutellidae was the non-monophyletic group. The family Attevidae was closely related to Praydidae. In addition，A. charopis was the sister to A. aurea. Our results provide a mitochondrial genome data for the intensive study of the phylogeny of Yponomeutoidea.

Key words: Lepidoptera, Yponomeutoidea, comparative mitochondrial genome, phylogeny

林兴雨, 宋南. Atteva charopis的比较线粒体基因组研究及系统发育分析[J]. 生物技术通报, 2023, 39(12): 300-310.

LIN Xing-yu, SONG Nan. Comparative Mitochondrial Genome and Phylogenetic Analysis of Atteva charopis Turner, 1903[J]. Biotechnology Bulletin, 2023, 39(12): 300-310.

图/表 11

参考文献 29

[1]	Sohn JC, Regier JC, Mitter C, et al. A molecular phylogeny for yponomeutoidea(insecta, Lepidoptera, ditrysia)and its implications for classification, biogeography and the evolution of host plant use[J]. PLoS One, 2013, 8(1): e55066. doi: 10.1371/journal.pone.0055066 URL
[2]	Slipinski SA, Leschen RAB, Lawrence JF. Order Coleoptera Linnaeus, 1758. In: Zhang, Z.-Q.(Ed.)Animal biodiversity: an outline of higher-level classification and survey of taxonomic richness[J]. Zootaxa, 2011, 3148(1): 203-208.
[3]	Powell JA, Mitter C, Farrell B, et al. Evolution of larval food preferences in Lepidoptera[J]. Kristensen NP, editor, 2013, 4: 403-422.
[4]	Grimaldi D, Engel MS. Evolution of the Insects[M]. Cambridge University Press, 2005.
[5]	Yang G, Zhang YN, Gurr GM, et al. Electroantennogram and behavioral responses of Cotesia plutellae to plant volatiles[J]. Insect Sci, 2016, 23(2): 245-252. doi: 10.1111/ins.2016.23.issue-2 URL
[6]	Leather S R. Insects on bird cherry. I. The bird cherry ermine moth, Yponomeuta evonymellus(L.)(Lepidoptera: Yponomeutidae)[J]. Entomologist's Gazette, 1986, 37(4): 209-213.
[7]	Hoebeke ER. Yponomeuta cagnagella(Lepidoptera: Yponomeutidae)₁: a Palearctic ermine moth in the United States, with notes on its recognition, seasonal history, and habits[J]. Ann Entomol Soc Am, 1987, 80(4): 462-467. doi: 10.1093/aesa/80.4.462 URL
[8]	Cameron SL. How to sequence and annotate insect mitochondrial genomes for systematic and comparative genomics research[J]. Syst Entomol, 2014, 39(3): 400-411. doi: 10.1111/syen.2014.39.issue-3 URL
[9]	Hua JM, Li M, Dong PZ, et al. Comparative and phylogenomic studies on the mitochondrial genomes of pentatomomorpha(Insecta: Hemiptera: Heteroptera)[J]. BMC Genomics, 2008, 9: 610. doi: 10.1186/1471-2164-9-610
[10]	Li H, Liu HY, Song F, et al. Comparative mitogenomic analysis of damsel bugs representing three tribes in the family Nabidae(In-secta: Hemiptera)[J]. PLoS One, 2012, 7(9): e45925. doi: 10.1371/journal.pone.0045925 URL
[11]	Ma C, Yang PC, Jiang F, et al. Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust[J]. Mol Ecol, 2012, 21(17): 4344-4358. doi: 10.1111/j.1365-294X.2012.05684.x pmid: 22738353
[12]	林兴雨, 翟卿, 宋南, 等. 锯谷盗线粒体基因组及扁甲总科系统发育分析[J]. 河南农业大学学报, 2023, 57(1): 109-117.
	Lin XY, Zhai Q, Song N, et al. The mitochondrial genome of Oryzaephilus surinamensis and a phylogenetic analysis of cucujoidea[J]. J Henan Agric Univ, 2023, 57(1): 109-117.
[13]	Patel RK, Jain M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data[J]. PLoS One, 2012, 7(2): e30619. doi: 10.1371/journal.pone.0030619 URL
[14]	Jin JJ, Yu WB, Yang JB, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes[J]. Genome Biol, 2020, 21(1): 241. doi: 10.1186/s13059-020-02154-5
[15]	Bernt M, Donath A, Jühling F, et al. MITOS: improved de novo metazoan mitochondrial genome annotation[J]. Mol Phylogenet Evol, 2013, 69(2): 313-319. doi: 10.1016/j.ympev.2012.08.023 URL
[16]	Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11[J]. Mol Biol Evol, 2021, 38(7): 3022-3027. doi: 10.1093/molbev/msab120 pmid: 33892491
[17]	Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes[J]. J Mol Evol, 1995, 41(3): 353-358. doi: 10.1007/BF00186547 pmid: 7563121
[18]	Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 2009, 25(11): 1451-1452. doi: 10.1093/bioinformatics/btp187 pmid: 19346325
[19]	Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability[J]. Mol Biol Evol, 2013, 30(4): 772-780. doi: 10.1093/molbev/mst010 pmid: 23329690
[20]	Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses[J]. Bioinformatics, 2009, 25(15): 1972-1973. doi: 10.1093/bioinformatics/btp348 pmid: 19505945
[21]	Kück P, Longo GC. FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies[J]. Front Zool, 2014, 11(1): 81. doi: 10.1186/s12983-014-0081-x pmid: 25426157
[22]	Nguyen LT, Schmidt HA, von Haeseler A, et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies[J]. Mol Biol Evol, 2015, 32(1): 268-274. doi: 10.1093/molbev/msu300 URL
[23]	Ronquist F, Teslenko M, van der Mark P, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space[J]. Syst Biol, 2012, 61(3): 539-542. doi: 10.1093/sysbio/sys029 pmid: 22357727
[24]	Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW(OGDRAW)version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes[J]. Nucleic Acids Res, 2019, 47(W1): W59-W64. doi: 10.1093/nar/gkz238 URL
[25]	Boore JL. Animal mitochondrial genomes[J]. Nucleic Acids Res, 1999, 27(8): 1767-1780. doi: 10.1093/nar/27.8.1767 pmid: 10101183
[26]	Huang XY, Chen B, Wei ZH, et al. First report of complete mitochondrial genome in the tribes coomaniellini and dicercini(Coleoptera: Buprestidae)and phylogenetic implications[J]. Genes, 2022, 13(6): 1074. doi: 10.3390/genes13061074 URL
[27]	Xiao LF, Zhang SD, Long CP, et al. Complete mitogenome of a leaf-mining buprestid beetle, Trachys auricollis, and its phylogenetic implications[J]. Genes, 2019, 10(12): 992. doi: 10.3390/genes10120992 URL
[28]	Shen X, Li X, Sha ZL, et al. Complete mitochondrial genome of the Japanese snapping shrimp Alpheus japonicus(Crustacea: Decapoda: Caridea): gene rearrangement and phylogeny within Caridea[J]. Sci China Life Sci, 2012, 55(7): 591-598. doi: 10.1007/s11427-012-4348-1 pmid: 22864833
[29]	Jeong JS, Park JS, Sohn JC, et al. The first complete mitochondrial genome in the family Attevidae(Atteva aurea)of the order Lepidoptera[J]. Biodivers Data J, 2022, 10: e89982. doi: 10.3897/BDJ.10.e89982 URL

科 Family	种名 Species	基因组数据库编号 GenBank No.	长度 Length/bp
潜叶蛾科 Lyonetiidae	Leucoptera malifoliella Lyonetia clerkella	JN790955 MF045483	15 646 15 259
Praydidae	Prays oleae	KM874804	16 499
斑巢蛾科 Attevidae	Atteva aurea Atteva charopis	ON480203 OQ130005	16 391 15 163
菜蛾科 Plutellidae	Acrolepiopsis assectella Plutella armoraciae	MW662615 MW662613	15 369 15 569
	Plutella porrectella	MW662614	16 196
	Plutella australiana	MG787473	15 962
	Plutella xylostella	KM023645	16 014
Scythropiidae	Scythropia crataegella	ON755208	15 350
外群 Outgroup	Caloptilia theivora Phyllocnistis citrella	MK541932 MN792920	15 297 15 416

科 Family	种名 Species	基因组数据库编号 GenBank No.	长度 Length/bp
潜叶蛾科 Lyonetiidae	Leucoptera malifoliella Lyonetia clerkella	JN790955 MF045483	15 646 15 259
Praydidae	Prays oleae	KM874804	16 499
斑巢蛾科 Attevidae	Atteva aurea Atteva charopis	ON480203 OQ130005	16 391 15 163
菜蛾科 Plutellidae	Acrolepiopsis assectella Plutella armoraciae	MW662615 MW662613	15 369 15 569
	Plutella porrectella	MW662614	16 196
	Plutella australiana	MG787473	15 962
	Plutella xylostella	KM023645	16 014
Scythropiidae	Scythropia crataegella	ON755208	15 350
外群 Outgroup	Caloptilia theivora Phyllocnistis citrella	MK541932 MN792920	15 297 15 416

基因 Gene	基因长度 Gene length/bp	起始位置 Start position/bp	终止位置 Stop position/bp	起始密码子 Start codon	终止密码子 Stop codon	编码链 Coding strand
trnI trnM	64 69	399 464	462 532			H H
trnQ	69	604	536			L
nad2	1 026	645	1 670	ATT	TAA	H
trnW	66	1 669	1 734			H
trnC	64	1 790	1 727			L
trnY	66	1 855	1 790			L
cox1	1 536	1 859	3 394	CGA	TAA	H
trnL2	67	3 390	3 456			H
cox2	682	3 457	4 138	ATG	T	H
trnK	71	4 139	4 209			H
trnD	67	4 209	4 275			H
atp8	159	4 276	4 434	ATT	TAA	H
atp6	672	4 428	5 099	ATG	TAA	H
cox3	789	5 099	5 887	ATG	TAA	H
trnG	67	5 890	5 956			H
nad3	354	5 957	6 310	ATA	TAA	H
trnA	66	6 324	6 389			H
trnR	62	6 391	6 452			H
trnN	66	6 453	6 518			H
trnS1	60	6 518	6 577			H
trnE	65	6 579	6 643			H
trnF	67	6 708	6 642			L
nad5	1 734	8 440	6 708	ATT	TA	L
trnH	64	8 504	8 441			L
nad4	1 342	9 846	8 505	ATG	T	L
nad4l	285	10 131	9 847	ATG	TAA	L
trnT	65	10 134	10 198			H
trnP	65	10 263	10 199			L
nad6	502	10 265	10 766	ATT	T	H
cob	1 152	10 799	11 950	ATG	TAA	H
trnS2	66	11 965	12 030			H
nad1	957	12 985	12 029	ATG	TAA	L
trnL1	68	13 054	12 987			L
rrnL	1 330	14 359	13 030			L
trnV	65	14 446	14 382			L
rrnS Control region	697 426	15 143 15 136	14 447 398			L Non-coding sequence

基因 Gene	基因长度 Gene length/bp	起始位置 Start position/bp	终止位置 Stop position/bp	起始密码子 Start codon	终止密码子 Stop codon	编码链 Coding strand
trnI trnM	64 69	399 464	462 532			H H
trnQ	69	604	536			L
nad2	1 026	645	1 670	ATT	TAA	H
trnW	66	1 669	1 734			H
trnC	64	1 790	1 727			L
trnY	66	1 855	1 790			L
cox1	1 536	1 859	3 394	CGA	TAA	H
trnL2	67	3 390	3 456			H
cox2	682	3 457	4 138	ATG	T	H
trnK	71	4 139	4 209			H
trnD	67	4 209	4 275			H
atp8	159	4 276	4 434	ATT	TAA	H
atp6	672	4 428	5 099	ATG	TAA	H
cox3	789	5 099	5 887	ATG	TAA	H
trnG	67	5 890	5 956			H
nad3	354	5 957	6 310	ATA	TAA	H
trnA	66	6 324	6 389			H
trnR	62	6 391	6 452			H
trnN	66	6 453	6 518			H
trnS1	60	6 518	6 577			H
trnE	65	6 579	6 643			H
trnF	67	6 708	6 642			L
nad5	1 734	8 440	6 708	ATT	TA	L
trnH	64	8 504	8 441			L
nad4	1 342	9 846	8 505	ATG	T	L
nad4l	285	10 131	9 847	ATG	TAA	L
trnT	65	10 134	10 198			H
trnP	65	10 263	10 199			L
nad6	502	10 265	10 766	ATT	T	H
cob	1 152	10 799	11 950	ATG	TAA	H
trnS2	66	11 965	12 030			H
nad1	957	12 985	12 029	ATG	TAA	L
trnL1	68	13 054	12 987			L
rrnL	1 330	14 359	13 030			L
trnV	65	14 446	14 382			L
rrnS Control region	697 426	15 143 15 136	14 447 398			L Non-coding sequence

特征 Feature	长度 Size/bp	AT/%	AT偏倚 AT-skew	GC偏倚 GC-skew
PCGs-H	6 870	78.30	-0.145	-0.113
PCGs-L	4 314	81.70	-0.171	0.316
rRNAs-L	2 027	85.10	0.040	0.344
tRNAs-H	921	80.90	0.031	0.023
tRNAs-L	528	83.50	0.016	0.425
全基因组 Whole genome	15 163	81.10	-0.024	-0.196