- Open Access
The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance
BMC Genomics volume 22, Article number: 755 (2021)
Mitochondrial genomes (mitogenomes) have greatly improved our understanding of the backbone phylogeny of Lepidoptera, but few studies on comparative mitogenomics below the family level have been conducted. Here, we generated 13 mitogenomes of eight tortricid species, reannotated 27 previously reported mitogenomes, and systematically performed a comparative analysis of nucleotide composition, gene variation and phylogenetic performance.
The lengths of completely sequenced mitogenomes ranged from 15,440 bp to 15,778 bp, and the gene content and organization were conserved in Tortricidae and typical for Lepidoptera. Analyses of AT-skew and GC-skew, the effective number of codons and the codon bias index all show a base bias in Tortricidae, with little heterogeneity among the major tortricid groups. Variations in the divergence rates among 13 protein-coding genes of the same tortricid subgroup and of the same PCG among tortricid subgroups were detected. The secondary structures of 22 transfer RNA genes and two ribosomal RNA genes were predicted and comparatively illustrated, showing evolutionary heterogeneity among different RNAs or different regions of the same RNA. The phylogenetic uncertainty of Enarmoniini in Tortricidae was confirmed. The synonymy of Bactrini and Olethreutini was confirmed for the first time, with the representative Bactrini consistently nesting in the Olethreutini clade. Nad6 exhibits the highest phylogenetic informativeness from the root to the tip of the resulting tree, and the combination of the third coding positions of 13 protein-coding genes shows extremely high phylogenetic informativeness.
This study presents 13 mitogenomes of eight tortricid species and represents the first detailed comparative mitogenomics study of Tortricidae. The results further our understanding of the evolutionary architectures of tortricid mitogenomes and provide a basis for future studies of population genetics and phylogenetic investigations in this group.
The mitochondrial genome (mitogenome) of insects is a circular double-stranded molecule and generally consists of 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs) and 22 transfer RNA genes (tRNAs) [1, 2]. In addition, several noncoding elements, including the control region regulating the replication and transcription of the mitogenome, are present . In recent years, the number of insect mitogenomes sequenced has dramatically increased, greatly improving our understanding of phylogenetics, species delimitation and identification, and population genetics for multiple insect groups [2, 4,5,6,7,8].
In insects, Lepidoptera is the second largest order after Coleoptera, with more than 157,000 described extant species in 43 superfamilies [9,10,11]. To date, the mitogenomes of approximately 550 lepidopteran species or subspecies have been sequenced and submitted to GenBank (https://www.ncbi.nlm.nih.gov/; last visited in September 2021). Utilizing thisdata, recent efforts of comparative mitogenomics have greatly advanced our understanding of evolution particularly the backbone phylogeny of Lepidoptera [4, 12,13,14,15,16]. However, studies on comparative mitogenomics below the family level are limited in Lepidoptera [2, 17].
The family Tortricidae, the sole member of the superfamily Tortricoidea, consists of approximately 11,365 described extant species, representing one of the most species-rich families in Lepidoptera [9, 18]. Despite the wide acceptance of the three-subfamily classification [19, 20], phylogenetic instabilities exist among prior investigations at the tribe level or below. In contrast to morphology [19,20,21], for instance, a multilocus study  suggests that Bactrini and Endothenini should be synonymized with Olethreutini of Olethreutinae. Using more taxon samplings than Regier et al. , Fagua et al.  found that the nonmonophyly of Olethreutini still exists. In addition, Tortricidae is notable for containing numerous important pest species that cause large losses in crop and forest production [24, 25], and many pest species (e.g., Grapholita spp.) are difficult to distinguish due to their morphological similarity . Thus, in studies addressing species delimitation and the population genetics of these species, molecular markers have been traditionally employed [26,27,28,29,30]. Among these markers, however, only mitochondrial cox1 and cox2 are predominately used [26,27,28,29,30], indicating the necessity of screening potential marker candidates through sequencing more tortricid mitogenomes, especially in the context of increasing challenges with the standard cox1 barcoding marker [31,32,33].
In Tortricidae, the mitogenomes of only 23 species from five tribes (Archipini, Eucosmini, Grapholitini, Olethreutini and Tortricini) have been sequenced (GenBank, August 2020). Moreover, most studies on tortricid mitogenomes focus on the description of a single mitogenome [34, 35]. Fagua et al.  sequenced and annotated six Choristoneura species of Archipini and emphasized the mitogenomic divergences of Choristoneura spp. associated with habit cooling cycles of the Northern Hemisphere in the Pliocene. Furthermore, the mitogenomic phylogeny based on 19 species from four tortricid tribes shows Eucosmini as paraphyletic.
In this study, we performed a thorough comparative mitogenomic analysis using 13 newly generated mitogenomes of eight tortricid species together with 28 previously reported mitogenomes, of which 27 were reannotated herein where necessary to 1) evaluate the nucleotide composition across major groups of Tortricidae; 2) compare the mitochondrial genetic variation at various taxonomic levels; and 3) assess the phylogenetic performances of various data partitions by constructing the preliminary phylogeny of Tortricidae. The evolutionary architectures of tortricid mitogenomes analysed herein will effectively facilitate further studies on the phylogeny and population genetics of Tortricidae and related groups.
Results and discussion
Generation of mitogenome data
Nine complete and four nearly complete mitogenomes were generated and annotated for eight species of Tortricidae. In the nearly complete genomes, we failed to assemble the control regions characterized by highly biased base composition, which was probably due to the disruption of sequencing reactions. All newly generated mitogenomes have been submitted to GenBank with the accession numbers shown in Additional file 1: Table S1.
The 13 mitogenomes each contained 37 typical mitochondrial genes in insects and showed identical gene organization to other reported tortricid mitogenomes which is also typical of Lepidoptera . The lengths of the completely sequenced mitogenomes ranged from 15,440 bp (L. koenigiana) to 15,778 bp (Olethreutes sp.), compared to other reported tortricid mitogenomes, which have ranged from 15,224 bp (Lobesia sp.) to 15,933 bp (A. fimbriana). The annotation information of mitogenomes sequenced herein is summarized in Additional file 2: Table S2.
When alignments were conducted among tortricid mitogenomes, some gene boundaries were ambiguous, as previously reported for other mitogenomes. To eliminate their potential impact on subsequent analyses, we carefully checked and revised the annotations of these genes in GenBank mainly according to the methods of Cameron . The reannotation information of 27 previously reported tortricid mitogenomes is summarized in Additional file 3: Table S3.
The overall nucleotide composition was A (40.5%), G (7.8%), C (11.5%) and T (40.2%), with a highly biased A + T content (80.7%), which is commonly present in insect mitogenomes . Among the six tribes of Tortricidae (Fig. 1a), the A + T contents ranged from 80.2% (Olethreutini) to 80.8% (Archipini and Grapholitini), showing little heterogeneity in nucleotide composition among the tortricid groups analysed in the present study. This is in contrast to some insects from the family to order taxonomic levels [5, 37,38,39,40]. Among the three codon positions within the 13 PCGs, the lowest A + T content was found for the first codon position, followed by the second and third codon positions, in accordance with most groups of insects, such as Zygaenoidea of Lepidoptera [14, 15] and Cimicomorpha of Hemiptera . Overall, the rRNAs showed a higher A + T content than the PCGs and tRNAs. The AT-skew and GC-skew are commonly used for evaluating the nucleotide composition of insect mitogenomes [41, 42]. In Tortricidae, negligible AT-skews and negative GC-skews were recognized (Fig. 1b, Additional file 4: Table S4), and six tortricid tribes consistently showed that the second codon positions of 13 PCGs and rRNAs had the lowest values of AT-skew and GC-skew, respectively, a feature commonly present in Lepidoptera .
Synonymous codons are generally used with different frequencies . For each amino acid, at least one synonymous codon showed a relative synonymous codon usage (RSCU) value greater than one in tortricid mitogenomes (Additional file 5: Table S5). The RSCU value of UUA encoding Leu, for instance, was higher than the sum of five other synonymous codons, whereas CUC was not used in Enarmoniini. The effective number of codons (ENC) and the codon bias index (CBI) were 30.52–35.22 and 0.65–0.75, respectively, and both exhibited codon usage bias among tortricid mitogenomes to some extent (Fig. 1c, d). Moreover, the positive correlation between the ENC and GC3s and the negative correlation between the CBI and GC3s (Fig. 1e, f) indicate that the genomic G + C content is a significant factor in determining codon bias among organisms [44,45,46]. In addition, UUA (Leu), AUU (Ile), UUU (Phe), AUA (Met) and AAU (Asn) represented the five most frequently used codons for all tribes, of which a high A + T content obviously contributed to the overall A + T bias of the whole mitogenome.
Mitochondrial gene variation of Tortricidae
To evaluate the variation patterns of 13 PCGs, genetic distance and nucleotide diversity were calculated. Within species, among species within genera and among genera within tribes, the averaged genetic distances for all PCGs were obviously separated and increased as the taxonomic ranks increased (Fig. 2a). However, the distance gaps were narrower between ranks above the genus level, and nad1, nad2 and atp8 even showed lower values between subfamilies than between tribes within subfamilies. For instance, the intraspecific distances of 13 PCGs in G. delineana ranged from 0 to 0.002. The distances among the three Archips species analysed were also variable across 13 PCGs, with the highest value being 0.1 for nad6. In contrast, the highest distance of 0.108was for cox3 between the two Adoxophyes species. Our results reflect the variation in rates of divergence among both PCGs and tortricid subgroups, as indicated by Fagua et al.  and in other insect groups such as the Miridae of Hemiptera .
Regarding the nucleotide diversity of 13 PCGs, sliding window analysis revealed a variable nucleotide diversity both within and among PCGs (Fig. 2b). The average values of nucleotide diversity for individual genes varied from 0.087 (nad5) to 0.165 (nad6). Similar to nad5, cox2, nad4, cox1, nad4l and nad1 showed relatively low nucleotide diversity, whereas the remaining PCGs, as well as nad6, exhibited relatively high nucleotide diversity. Nucleotide diversity is commonly used for identifying regions with high nucleotide divergence and could provide guidelines for selecting species-specific markers, especially in taxa with high morphological similarity [47, 48]. The genes with higher levels of divergence identified herein would provide potential marker candidates for population genetics and species delimitation in Tortricidae.
To characterize the variation distribution of the 24 RNAs, their secondary structures were predicted and comparatively illustrated using three Archips species as an example. Comparative tRNA analyses (Additional file 6: Fig. S1) showed that the tRNA structures, including the loss of the DHU arm in trnS1 (AGN), were highly conserved across the three congeneric mitogenomes, a feature commonly present in lepidopteran insects as well as other metazoans [49, 50]. The nucleotide substitution rates of the 22 tRNAs were different, with trnW showing no variation and trnD and trnM having higher nucleotide substitutions. In addition, for each tRNA, more variable sites were generally present in the TψC arm, TψC loop and DHU loop, a variation distribution similar to the Macroheterocera of Lepidoptera [14, 15]. The secondary structures of the two rRNAswere generally identical to those proposed for some lepidopterans, especially tortricid species [36, 43, 51, 52]. In rrnS (Fig. 3), three domains (I–III) were recognized in three Archips species (Fig. 4), and the variable sites showed a scatter distribution across three domains. H1047 and H1074 in rrnS are highly variable among insect orders [43, 53]. In our analyses, H1074 was highly conserved among the three Archips species, indicating that its variation level may increase with higher taxonomic ranks. Regarding the secondary structure of rrnL (Fig. 4), five domains (I–II, IV–VI) were detected, with domain III absent, as in other insects [43, 54]. In rrnL, a typical feature is the existence of a microsatellite sequence of (TA) n in the stem region of H2347, and the difference in the repeat number among three Archips species makes it a highly variable region. This feature has also been found in some Grapholita species of Tortricidae .
By adding 13 newly sequenced mitogenomes to 28 existing mitogenomes from GenBank, we performed various phylogenetic analyses to test their phylogenetic implications and to evaluate the phylogenetic performance of different data partitions.
Using four datasets, maximum likelihood (ML) and Bayesian inference (BI) analyses consistently recovered the six tribes and two subfamilies as monophyletic, except for Eucosmini in the ML analysis of the P12 dataset (Fig. 5, Additional file 7–12: Figs. S2–S7). In the subfamily Olethreutinae, four tribes were included in this study, and Enarmoniini was sampled for the first time. (Enarmoniini + (Olethreutini + (Eucosmini + Grapholitini))) was recovered, of which the sister group between Eucosmini and Grapholitini was also supported by Regier et al.  and Fagua et al.  based on multilocus data. The present study consistently showed Grapholita as nonmonophyletic, confirming our recent study  and previous multilocus studies [22, 23]. Our results consistently placed B. venosana, historically belonging to Bactrini, into Olethreutini, reinforcing the synonymy of the two tribes , based on mitogenome evidence for the first time. Six G. delineana individuals from different locations were sampled to compare mitochondrial gene variation at the species level, and the phylogenetic results confirmed their position in Grapholitini. Regarding Tortricinae, two tribes were sampled in the present study. In our recent study , the position of Epiphyas in Archipini was not unstable, being sister to either Choristoneura or to Adoxophyes. The inclusion of three Archips spp. revealed the Epiphyas was consistently sister to the (Choristoneura + Archips) across all datasets, although this relationship remains to be clarified with increased sampling for Archipini. In brief, although the Enarmoniini and Olethreutini were represented for the first time, the present mitogenomic phylogeny only included six of the 19 tribes in Tortricidae . Thus, the relationships among the six tortricid tribes recovered herein definitely need further confirmation with increased sampling.
Despite the high efficiency of using mitogenomes in insect phylogenetic inference, systematic evaluation of the relative contribution of each mitochondrial gene or data partition to the resulting trees has seldom been conducted, especially below the family level. In this study, phylogenetic informativeness (PI) was calculated to assess the contribution of each data partition to the phylogenetic tree from the P123RT dataset (Fig. 6). The PI curves for each data partition were similar in shape, with a steady increase from the root to a peak in the tree showing genus-level relationships and a rapid decrease closer to the tips. Among the 13 PCGs, nad6 exhibited the highest PI from the root to the tip of the resulting tree, whereas nad4, nad5, and cox3 showed lower PI. The high PI of nad6 was also revealed by Nie et al. , who performed a phylogenetic analysis of Galerucinae in Coleoptera. Mitochondrial rrnL has been demonstrated to be highly informative in inferring tribe- and subtribe-level relationships in Satyrinae by Yang and Zhang [57, 58]. A similar result was revealed in the present study, and 16S rDNA had a relatively high PI along with nad6. Overall, the 22 tRNAs showed the lowest PI. Furthermore, the three coding positions of 13 PCGs were compared, and the combination of the third coding positions of 13 PCGs showed extremely high PI, as expected. This result was also demonstrated by our phylogenetic analyses which showed higher support values on most nodes of the P123 dataset than on the P12 dataset.
Samples, DNA extraction and mitogenome sequencing
Adults were collected at various locations in China from 2018 to 2019. Specimens were preserved in 100% ethanol in a − 80 °C environment until they were used for DNA extraction. Species identification was conducted through morphology [59, 60] or/and standard mitochondrial cox1 barcoding . A total of 13 samples representing eight species were selected for sequencing. Among them, L. koenigiana is the first sequenced species for Enarmoniini. Three congeneric species, A. podana, A. betulanus and Archips sp. were sampled as an example to evaluate mitochondrial gene variation at the genus level. Six and one samples of G. delineana and A. fimbriana were included respectively to test the intraspecific variation of mitochondrial genes among population individuals. In addition, two Olethreutini species were selected. Olethreutes sp. represented the first sequenced species of Olethreutes, the nominal genus of Olethreutini, and B. venosana, historically belonging to Bactrini, was included to test the synonym of Bactrini and Olethreutini proposed by Regier et al. .
Total genomic DNA was extracted from head and thorax tissues of a single sample using a DNeasy tissue kit (Qiagen, Germany), following the manufacturer’s instructions. The libraries of 13 samples were individually constructed and sequencing was conducted using an Illumina HiSeq 2500 platform with a 150-bp paired-end strategies. Voucher specimens were deposited in the Biology Laboratory of Zhoukou Normal University, China (Additional file 13: Table S6).
Mitogenome assembly, annotation and analysis
FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) was used for quality control of the raw sequences. Clean paired reads were obtained using AdapterRemoval version 2  and SOAPdenovo version. 2.01 . Geneious R11  was used for mitogenome assembly with default settings. Briefly, the “map to reference” strategy was selected to map all cleaned reads to an “anchor” that represents the standard mitochondrial cox1 barcoding sequence amplified earlier using the general insect primer pair Lco1490 (F) and Hco2198 (R) . After iteration up to 100 times with custom sensitivity, a target sequence with high coverage contigs was generated. Then, MEGA X  was used to check the beginning and end of the contig sequence to circularize a complete mitochondrial genome after deleting the overlapping sequence.
The mitogenome sequence was annotated using the MITOS webserver with invertebrate genetic code . MEGA X was used to reconfirm gene boundaries by aligning the new mitogenome with previously reported tortricid mitogenomes available in GenBank. The 22 tRNAs and their secondary structures were reidentified using tRNAScan-SE server version 1.21 . The secondary structures of the two rRNAs were inferred following the models of three other lepidopterans with minor modifications [36, 43, 51].
Nucleotide composition and averaged genetic distances under the Kimura-2-parameter model of 13 PCGs were calculated using MEGA X. Strand asymmetry was calculated according to the formulas: AT-skew = [A – T]/[A + T] and GC-skew = [G – C]/[G + C] . Sliding window analysis (a sliding window of 300 bp and a step size of 25 bp) exhibiting nucleotide diversity for each of 13 PCGs from all samples was performed using DNASP version 6.0 . Additionally, DNASP version 6.0 was used to calculate the ratios of nonsynonymous substitution (Ka) and synonymous substitution (Ks) for PCGs. In addition, the effective number of codons (ENC) and the codon bias index (CBI) were measured using CodonW version 1.4.2 .
A total of 41 mitogenomes from 30 tortricid species were used in phylogenetic analyses, including 13 newly sequenced mitogenomes and 28 mitogenomes downloaded from GenBank (Table S1). In addition, two species from Cossoidea served as outgroup taxa. Thirteen PCGs were individually aligned with codon-based mode in the TranslatorX online platform . Two rRNAs and 22 tRNAs were independently aligned with the Q-INS-i algorithm as implemented in the MAFFT online platform . MEGA X was used to check all alignments, and PhyloSuite version 1.2.1  was employed to generate four datasets: 1) P12: the first and second codon positions of PCGs; 2) P123: all codon positions of PCGs; 3) P123R: all codon positions of PCGs plus two RNAs; and 4) P123RT: all codon positions of PCGs plus two RNAs and 22 tRNAs. In addition, independent runs for 13 PCGs and two rRNAs were performed to comparatively evaluate their phylogenetic performance.
Maximum likelihood (ML) analyses were conducted using IQ-TREE 2.0.4  under the partitioning schemes and corresponding substitution models (Additional file 14: Table S7) determined by ModelFinder . Branch supports were calculated using 1000 ultrafast bootstrap replicates . Bayesian inference (BI) analyses were performed with MrBayes version 3.2.6  with the partitioned models (Additional file 15: Table S8) determined by PartitionFinder version 2.1.1 . Twelve processors were used to perform three independent runs each with four chains (three heated and one cold) simultaneously for more than 10,000,000 generations sampled every 100 generations. Convergences were considered to be reached when the estimated sample size (ESS) value was above 200 established by Tracer version 1.7  and the potential scale reduction factor (PSRF) approached 1.0 . The first 25% of samples were discarded as burn-in and the remaining trees were used to calculate posterior probabilities in a 50% majority-rule consensus tree.
Phylogenetic informativeness (PI) profiles were used to quantify the relative contributions of 13 PCGs and two rRNAs to the resulting tree. The peak of the PI distribution is suggested to predict maximum phylogenetic informativeness for corresponding gene partitioning . To obtain PI profiles, PhyDesign [81, 82] was used with the aligned sequences and an ultrametric tree as input files. The ultrametric tree was constructed using BEAST version 1.7.5 .
Availability of data and materials
The data that support the findings of this study are openly available in GenBank of NCBI at https://www.ncbi.nlm.nih.gov, with the reference numbers MH013482, MW936632–MW936633, MW924656–MW924665.
- atp6 and atp8 :
ATPase subunits 6 and 8
- cob :
- cox1–cox3 :
Cytochrome c oxidase subunits I-III
- nad1–nad6 and nad4l :
NADH dehydrogenase subunits 1–6 and 4 L
- rrnL and rrnS :
Large and small subunit ribosomal RNA (rRNA)
Relative synonymous codon usage
Curole JP, Kocher TD. Mitogenomics: digging deeper with complete mitochondrial genomes. Trends Ecol Evol. 1999;14(10):394–8. https://doi.org/10.1016/S0169-5347(99)01660-2.
Cameron SL. Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014;59(1):95–117. https://doi.org/10.1146/annurev-ento-011613-162007.
Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27(8):1767–80. https://doi.org/10.1093/nar/27.8.1767.
Timmermans MJTN, Lees DC, Simonsen TJ. Towards a mitogenomic phylogeny of Lepidoptera. Mol Phylogenet Evol. 2014;79:169–78. https://doi.org/10.1016/j.ympev.2014.05.031.
Liu Y, Song F, Jiang P, Wilson JJ, Cai W, Li H. Compositional heterogeneity in true bug mitochondrial phylogenomics. Mol Phylogenet Evol. 2018;118:135–44. https://doi.org/10.1016/j.ympev.2017.09.025.
Du ZY, Hasegawa H, Cooley JR, Simon C, Yoshimura J, Cai W, et al. Mitochondrial genomics reveals shared phylogeographic patterns and demographic history among three periodical cicada species groups. Mol Biol Evol. 2019;36(6):1187–200. https://doi.org/10.1093/molbev/msz051.
Song N, Zhang H, Zhao T. Insights into the phylogeny of Hemiptera from increased mitogenomic taxon sampling. Mol Phylogenet Evol. 2019;137:239–49. https://doi.org/10.1016/j.ympev.2019.05.009.
Nie R, Vogler AP, Yang XK, Lin M. Higher-level phylogeny of longhorn beetles (Coleoptera: Chrysomeloidea) inferred from mitochondrial genomes. Syst Entomol. 2021;46(1):56–70. https://doi.org/10.1111/syen.12447.
van Nieukerken EJ, Kaila L, Kitching IJ, Kristensen NP, Lees DC. Minet, Joël; Mitter JM, et al. Order Lepidoptera Linnaeus, 1758. Zootaxa. 2011;3148:212–21.
Mitter C, Davis DR, Cummings MP. Phylogeny and evolution of Lepidoptera. Annu Rev Entomol. 2017;62(1):265–83. https://doi.org/10.1146/annurev-ento-031616-035125.
Triant DA, Cinel SD, Kawahara AY. Lepidoptera genomes: current knowledge, gaps and future directions. Curr Opin Insect Sci. 2018;25:99–105. https://doi.org/10.1016/j.cois.2017.12.004.
Kim MJ, Kang AR, Jeong HC, Kim KG, Kim I. Reconstructing intraordinal relationships in Lepidoptera using mitochondrial genome data with the description of two newly sequenced lycaenids, Spindasis takanonis and Protantigius superans (Lepidoptera: Lycaenidae). Mol Phylogenet Evol. 2011;61(2):436–45. https://doi.org/10.1016/j.ympev.2011.07.013.
Yang X, Cameron SL, Lees DC, Xue D, Han H. A mitochondrial genome phylogeny of owlet moths (Lepidoptera: Noctuoidea), and examination of the utility of mitochondrial genomes for lepidopteran phylogenetics. Mol Phylogenet Evol. 2015;85:230–7. https://doi.org/10.1016/j.ympev.2015.02.005.
Yang M, Zhang H, Song L, Shi Y, Liu X. The complete mitochondrial genome of Mahanta tanyae compared with other zygaenoid moths (Lepidoptera: Zygaenoidea). J Asia Pac Entomol. 2019;22(2):513–21. https://doi.org/10.1016/j.aspen.2019.03.010.
Yang M, Song L, Shi Y, Li J, Zhang Y, Song N. The first mitochondrial genome of the family Epicopeiidae and higher-level phylogeny of Macroheterocera (Lepidoptera: Ditrysia). Int J Biol Macromol. 2019;136:123–32. https://doi.org/10.1016/j.ijbiomac.2019.06.051.
Kim MJ, Kim J, Cameron SL. How well do multispecies coalescent methods perform with mitochondrial genomic data? A case study of butterflies and moths (Insecta: Lepidoptera). Syst Entomol. 2020;45:857–73. https://doi.org/10.1111/syen.12431.
Fagua G, Condamine FL, Brunet BMT, Clamens AL, Laroche J, Levesque RC, et al. Convergent herbivory on Pinaceae by Choristoneura moths after boreal forest formation. Mol Phylogenet Evol. 2018;123:35–43. https://doi.org/10.1016/j.ympev.2018.01.013.
Gilligan TM, Baixeras J, Brown JW. T@RTS: Online World Catalogue of the Tortricidae (Version 4.0). 2018. http://www.tortricidae.com/catalogueGenusList.asp?gcode=448.
Horak M, Brown RL. 1.2 Taxonomy and phylogeny. In: van der Geest LPS, HH Evenhuis, eds. Tortricid pests, their biology, natural enemies and control. Amsterdam: Elsevier Science Publishers B.V. 1991; 23–48.
Horak M. Olethreutine moths of Australia (Lepidoptera: Tortricidae), Monographs on Australian Lepidoptera; 2006. https://doi.org/10.1071/9780643094086.
Brown JW. World catalogue of insects, Vol. 5: Tortricidae (Lepidoptera). Stenstrup: Apollo Books; 2005. https://doi.org/10.1163/9789004475380.
Regier JC, Brown JW, Mitter C, Baixeras J, Cho S, Cummings MP, et al. A molecular phylogeny for the leaf-roller moths (Lepidoptera: Tortricidae) and its implications for classification and life history evolution. PLoS One. 2012;7(4):e35574. https://doi.org/10.1371/journal.pone.0035574.
Fagua G, Condamine FL, Horak M, Zwick A, Sperling FAH. Diversification shifts in leafroller moths linked to continental colonization and the rise of angiosperms. Cladistics. 2017;33(5):449–66. https://doi.org/10.1111/cla.12185.
van der Geest LPS, Evinhuis HH. Tortricid pests, their biology, natural enemies and control. Amsterdam: Elsevier Science Publishers B.V; 1991.
Dombroskie J, Sperling FAH. Phylogeny of the tribe Archipini (Lepidoptera: Tortricidae: Tortricinae) and evolutionary correlates of novel secondary sexual structures. Zootaxa. 2013;3729(1):001–62. https://doi.org/10.11646/zootaxa.3729.1.1.
Zheng Y, Wu RX, Dorn S, Chen MH. Diversity of tortricid moths in appleorchards: evidence for a cryptic species of Grapholita (Lepidoptera: Tortricidae) from China. Bull Entomol Res. 2017;107(2):268–80. https://doi.org/10.1017/S0007485316000973.
Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT. Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu Rev Ecol Syst. 2006;37(1):545–79. https://doi.org/10.1146/annurev.ecolsys.37.091305.110018.
Zheng Y, Peng X, Liu G, Pan H, Dorn S, Chen M. High genetic diversity and structured populations of the oriental fruit moth in its range of origin. PLoS One. 2013;8(11):e78476. https://doi.org/10.1371/journal.pone.0078476.
Wei S, Cao L, Gong Y, Shi B, Wang S, Zhang F, et al. Population genetic structure and approximate Bayesian computation analyses reveal the southern origin and northward dispersal of the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae) in its native range. Mol Ecol. 2015;24(16):4094–111. https://doi.org/10.1111/mec.13300.
Wang K, Li Y, Zheng Y, Duan X, Zhang M, Peng X, et al. Molecular identification of Grapholita dimorpha Komai (Lepidoptera: Tortricidae) and Grapholita molesta busk (Lepidoptera: Tortricidae) based on mitochondrial COI and COII genes. J Northwest A&F Univ. 2016;44:156–64.
Sheffield NC, Song H, Cameron SL, Whiting MF. A comparative analysis of mitochondrial genomes in Coleoptera (Arthropoda: Insecta) and genome descriptions of six new beetles. Mol Biol Evol. 2008;25(11):2499–509. https://doi.org/10.1093/molbev/msn198.
Lee W, Park J, Lee GS, Lee S, Akimoto S. Taxonomic status of the Bemisia tabaci complex (Hemiptera: Aleyrodidae) and reassessment of the number of its constituent species. PLoS One. 2013;8(5):e63817. https://doi.org/10.1371/journal.pone.0063817.
Wang J, Zhang L, Zhang QL, Zhou MQ, Wang XT, Yang XZ, et al. Comparative mitogenomic analysis of mirid bugs (Hemiptera: Miridae) and evaluation of potential DNA barcoding markers. Peer J. 2017;5:e3661. https://doi.org/10.7717/peerj.3661.
Zhu JY, Wei SJ, Li QW, Yang S, Li YH. Mitochondrial genome of the pine tip moth Rhyacionia leptotubula (Lepidoptera: Tortricidae). Mitochondrial DNA. 2012;23(5):376–8. https://doi.org/10.3109/19401736.2012.696632.
Zhao JL, Wu YP, Su TJ, Jiang GF, Wu CS, Zhu CD. The complete mitochondrial genome of Acleris fimbriana (Lepidoptera: Tortricidae). Mitochondrial DNA Part A. 2016;27:2200–2.
Wu YP, Zhao JL, Su TJ, Luo AR, Zhu CD. The complete mitochondrial genome of Choristoneura longicellana (Lepidoptera: Tortricidae) and phylogenetic analysis of Lepidoptera. Gene. 2016;591(1):161–76. https://doi.org/10.1016/j.gene.2016.07.003.
Song F, Li H, Jiang P, Zhou X, Liu J, Sun C, et al. Capturing the phylogeny of holometabola with mitochondrial genome data and Bayesian site-heterogeneous mixture models. Genome Biol Evol. 2016;8(5):1411–26. https://doi.org/10.1093/gbe/evw086.
Yang H, Li T, Dang K, Bu W. Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera). BMC Genomics. 2018;19(1):264. https://doi.org/10.1186/s12864-018-4650-9.
Tang P, Zhu J, Zheng B, Wei S, Sharkey M, Chen X, et al. Mitochondrial phylogenomics of the Hymenoptera. Mol Phylogenet Evol. 2019;131:8–18. https://doi.org/10.1016/j.ympev.2018.10.040.
Nie R, Andújar C, Gómez-Rodríguez C, Bai M, Xue HJ, Tang M, et al. The phylogeny of leaf beetles (Chrysomelidae) inferred from mitochondrial genomes. Syst Entomol. 2020;45(1):188–204. https://doi.org/10.1111/syen.12387.
Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41(3):353–8. https://doi.org/10.1007/BF01215182.
Wei SJ, Shi M, Chen XX, Sharkey MJ, van Achterberg C, Ye GY, et al. New views on strand asymmetry in insect mitochondrial genomes. PLoS One. 2010;5(9):e12708. https://doi.org/10.1371/journal.pone.0012708.
Cameron SL, Whiting MF. The complete mitochondrial genome of the tobacco hornworm, Manduca sexta (Insecta: Lepidoptera: Sphingidae), and an examination of mitochondrial gene variability within butterflies and moths. Gene. 2008;408(1-2):112–23. https://doi.org/10.1016/j.gene.2007.10.023.
Hershberg R, Petrov DA. Selection on codon Bias. Annu Rev Genet. 2008;42(1):287–99. https://doi.org/10.1146/annurev.genet.42.110807.091442.
Plotkin JB, Kudla G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2011;12(1):32–42. https://doi.org/10.1038/nrg2899.
Yuan ML, Zhang QL, Guo ZL, Wang J, Shen YY. Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications. BMC Genomics. 2015;16(1):460. https://doi.org/10.1186/s12864-015-1679-x.
Jia WZ, Yan HB, Guo AJ, Zhu XQ, Wang YC, Shi WG, et al. Complete mitochondrial genomes of Taenia multiceps, T. hydatigena and T. pisiformis: additional molecular markers for a tapeworm genus of human and animal health significance. BMC Genomics. 2010;11(1):447. https://doi.org/10.1186/1471-2164-11-447.
Ma LY, Liu FF, Chiba H, Yuan XQ. The mitochondrial genomes of three skippers: insights into the evolution of the family Hesperiidae (Lepidoptera). Genomics. 2020;112(1):432–41. https://doi.org/10.1016/j.ygeno.2019.03.006.
Garey JR, Wolstenholme DR. Platyhelminth mitochondrialDNA: evidence for early evolutionary origin of a tRNAserAGN that contains a dihydrouridine arm replacement loop, and of serine-specifying AGA and AGG codons. J Mol Evol. 1989;28(5):374–87. https://doi.org/10.1007/BF02603072.
Lavrov DV, Brown WM, Boore JL. A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus. Proc Natl Acad Sci U. S. A. 2000;97(25):13738–42. https://doi.org/10.1073/pnas.250402997.
Gong YJ, Shi BC, Kang ZJ, Zhang F, Wei SJ. The complete mitochondrial genome of the oriental fruit moth Grapholita moleata (Buck) (Lepidoptera: Tortricidae). Mol Biol Rep. 2012;39(3):2893–900. https://doi.org/10.1007/s11033-011-1049-y.
Wu YP, Zhao JL, Su TJ, Li J, Yu F, Chesters D, et al. The complete mitochondrial genome of Leucoptera malifoliella Costa (Lepidoptera: Lyonetiidae). DNA Cell Biol. 2012;31(10):1508–22. https://doi.org/10.1089/dna.2012.1642.
Gillespie JJ, Johnston JS, Cannone JJ, Gutell RR. Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements. Insect Mol Biol. 2006;15(5):657–86. https://doi.org/10.1111/j.1365-2583.2006.00689.x.
Su T, He B, Li K, Liang A. Comparative analysis of the mitochondrial genomes of oriental spittlebug tribe Cosmoscartini: insights into the relationships among closely related taxa. BMC Genomics. 2018;19(1):961. https://doi.org/10.1186/s12864-018-5365-7.
Song L, Shi Y, Zhang H, Wang Z, Liu X, Yang M. Complete mitochondrial genome of the hemp borer, Grapholita delineana (Lepidoptera: Tortricidae): gene variability and phylogeny among Grapholita. J Asia Pac Entomol. 2021;24(2):250–8. https://doi.org/10.1016/j.aspen.2021.02.002.
Nie RE, Breeschoten T, Timmermans MJTN, Nadein K, Xue HJ, Bai M, et al. The phylogeny of Galerucinae (Coleoptera: Chrysomelidae) and the performance of mitochondrial genomes in phylogenetic inference compared to nuclear rRNA genes. Cladistics. 2018;34(2):113–30. https://doi.org/10.1111/cla.12196.
Yang M, Zhang Y. Phylogenetic utility of ribosomal genes for reconstructing the phylogeny of five Chinese satyrine tribes (Lepidoptera, Nymphalidae). ZooKeys. 2015a;488:105–20. https://doi.org/10.3897/zookeys.488.9171.
Yang M, Zhang Y. Molecular phylogeny of the butterfly tribe Satyrini with emphasis on the utility of ribosomal mitochondrial genes 16s rDNA and nuclear 28s rDNA (Nymphalidae: Satyrinae). Zootaxa. 2015b;3985(1):125–41. https://doi.org/10.11646/zootaxa.3985.1.7.
Li H. Microlepidoptera of Qinling Moutains (Insecta: Lepidoptera). Beijing: Science Press; 2012.
Liu YQ, Li GW. Fauna Sinica. Vol. 27: Lepidoptera: Tortricidae. Beijing: Science Press; 2002.
Hebert PDN, Ratnasingham S, de Waard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc Royal Soc. 2003;270:96–9.
Schubert M, Lindgreen S, Orlando L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging, BMC res. Notes. 2016;9:88.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1(1):18. https://doi.org/10.1186/2047-217X-1-18.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9. https://doi.org/10.1093/bioinformatics/bts199.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9. https://doi.org/10.1093/molbev/msy096.
Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, et al. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69(2):313–9. https://doi.org/10.1016/j.ympev.2012.08.023.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64. https://doi.org/10.1093/nar/25.5.955.
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large datasets. Mol Biol Evol. 2017;34(12):3299–302. https://doi.org/10.1093/molbev/msx248.
Peden JF. Analysis of codon usage. Univ Nottingham. 2000;90:73–4.
Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38(suppl_2):7–13. https://doi.org/10.1093/nar/gkq291.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6. https://doi.org/10.1093/bib/bbx108.
Zhang D, Gao F, Jakovlić I, Zou H, Zhang J, Li WX, et al. PhyloSuite: an integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol Ecol Resour. 2020;20(1):348–55. https://doi.org/10.1111/1755-0998.13096.
Nguyen LT, Schmidt HAA, von Haeseler A, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74. https://doi.org/10.1093/molbev/msu300.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9. https://doi.org/10.1038/nmeth.4285.
Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35(2):518–22. https://doi.org/10.1093/molbev/msx281.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42. https://doi.org/10.1093/sysbio/sys029.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution formolecular and morphological phylogenetic analyses. Mol Biol Evol. 2017;34(3):772–3. https://doi.org/10.1093/molbev/msw260.
Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarisation in Bayesian phylogenetics using tracer 1.7. Syst Biol. 2018;67(5):901–4. https://doi.org/10.1093/sysbio/syy032.
Owen CL, Marshall DC, Hill KBR, Simon C. The phylogenetic utility of acetyltransferase (ARD1) and glutaminyl tRNA synthetase (QtRNA) for reconstructing Cenozoic relationships as exemplified by the large Australian cicada Pauropsalta generic complex. Mol Phylogenet Evol. 2014;83:258–77. https://doi.org/10.1016/j.ympev.2014.07.008.
Townsend JP. Profiling phylogenetic informativeness. Syst Biol. 2007;56(2):222–31. https://doi.org/10.1080/10635150701311362.
Lopez-Giraldez F, Townsend JP. PhyDesign: an online application for profiling phylogenetic informativeness. BMC Evol Biol. 2011;11(1):152. https://doi.org/10.1186/1471-2148-11-152.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73. https://doi.org/10.1093/molbev/mss075.
We sincerely appreciate the anonymous reviewers for comments on this manuscript.
This work was funded by the National Natural Science Foundation of China (31702046), Science and Technology Project of Henan Province (202102110223), Key Scientific Research Project of Henan Colleges and Universities (21A210030) and the Project of Scientific Research Innovation Fund for College Student (ZKNUD2020035 and ZKNUD2020036).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The tortricid samples used in phylogenetic analyses.
Annotation and comparison of mitochondrial genome organizations of 13 mitogenomes sequenced in this study.
Reannotation and comparison of mitochondrial genome organizations of 27 previous sequenced tortricid species.
GC-Skew of the mitogenomes in six tortricid tribes.
Codon usage pattern in six tortricid tribes analyzed in this study.
Putative secondary structures of tRNAs from three Archips mitogenomes. The tRNAs are labeled with the abbreviations of their corresponding amino acids. The tRNA arms are illustrated as for trnV. Dashes indicate the Watson-Crick base pairs; dots indicate the wobble GU pairs; and the other non-canonical pairs are not marked.
Phylogenetic tree inferred from maximum likelihood method based on P12 dataset. Numbers on nodes represent the bootstrap supports.
Phylogenetic tree inferred from Bayesian inference method based on P12 dataset. Numbers on nodes represent the posterior probabilities.
Phylogenetic tree inferred from maximum likelihood method based on P123 dataset. Numbers on nodes represent the bootstrap supports.
Phylogenetic tree inferred from Bayesian inference method based on P123 dataset. Numbers on nodes represent the posterior probabilities.
Phylogenetic tree inferred from maximum likelihood method based on P123R dataset. Numbers on nodes represent the bootstrap supports.
Phylogenetic tree inferred from Bayesian inference method based on P123R dataset. Numbers on nodes represent the posterior probabilities.
Information of samples sequenced in this study.
The partitioning schemes and corresponding substitution models determined by ModelFinder.
The partitioning schemes and corresponding substitution models determined by PartitionFinder.
About this article
Cite this article
Yang, M., Li, J., Su, S. et al. The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance. BMC Genomics 22, 755 (2021). https://doi.org/10.1186/s12864-021-08041-y
- Leaf roller moths