Skip to main content

The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance



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.

Peer Review reports


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 [3]. 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 (; 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 [22] suggests that Bactrini and Endothenini should be synonymized with Olethreutini of Olethreutinae. Using more taxon samplings than Regier et al. [22], Fagua et al. [23] 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 [26]. 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. [17] 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 [36]. 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 [2]. The reannotation information of 27 previously reported tortricid mitogenomes is summarized in Additional file 3: Table S3.

Nucleotide composition

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 [3]. 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 [38]. 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 [43].

Fig. 1
figure 1

Nucleotide compositions of the tortricid mitogenomes. a A + T content; b AT-skew; c The averaged effective number of codons (ENC) of six tortricid tribes; d The averaged codon bias index (CBI) of six tortricid tribes; e Scatter plot of the GC content of 3rd codon sites versus ENC; f Scatter plot of the GC content of 3rd codon sites versus CBI

Synonymous codons are generally used with different frequencies [44]. 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. [17] and in other insect groups such as the Miridae of Hemiptera [33].

Fig. 2
figure 2

Mitochondrial gene variation of Tortricidae. a Mean genetic distances (Kimura-2-parameter) of mitochondrial gene sequences among Tortricidae for different taxonomic ranks; b Sliding window analysis shows the value of nucleotide diversity. Gene names and the values of averaged nucleotide diversity are attached

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 [55].

Fig. 3
figure 3

Putative secondary structures of rrnS from three Archips mitogenomes. Dashes indicate the Watson-Crick base pairs; dots indicate the wobble GU pairs; and the other noncanonical pairs are not marked. Roman numerals denote the conserved domain structure

Fig. 4
figure 4

Putative secondary structures of rrnL from three Archips mitogenomes. Dashes indicate the Watson-Crick base pairs; dots indicate the wobble GU pairs; and the other noncanonical pairs are not marked. Roman numerals denote the conserved domain structure

Phylogenetic informativeness

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 712: 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. [22] and Fagua et al. [23] based on multilocus data. The present study consistently showed Grapholita as nonmonophyletic, confirming our recent study [22] 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 [22], 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 [55], 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 [18]. Thus, the relationships among the six tortricid tribes recovered herein definitely need further confirmation with increased sampling.

Fig. 5
figure 5

Phylogenetic trees inferred from maximum likelihood and Bayesian inference methods based on the P123RT dataset. Numbers separated by a slash (/) on a node represent the bootstrap support and posterior probability, respectively. The species or individual samples with mitogenomes sequenced in this study are marked in bold

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. [56], 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.

Fig. 6
figure 6

Phylogenetic informative profiles for mitochondrial data partitions. The ultrametric tree was constructed using the P123RT 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 [61]. 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. [22].

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 ( was used for quality control of the raw sequences. Clean paired reads were obtained using AdapterRemoval version 2 [62] and SOAPdenovo version. 2.01 [63]. Geneious R11 [64] 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) [65]. After iteration up to 100 times with custom sensitivity, a target sequence with high coverage contigs was generated. Then, MEGA X [66] 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 [67]. 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 [68]. 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] [41]. 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 [69]. 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 [70].

Phylogenetic analyses

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 [71]. Two rRNAs and 22 tRNAs were independently aligned with the Q-INS-i algorithm as implemented in the MAFFT online platform [72]. MEGA X was used to check all alignments, and PhyloSuite version 1.2.1 [73] 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 [74] under the partitioning schemes and corresponding substitution models (Additional file 14: Table S7) determined by ModelFinder [75]. Branch supports were calculated using 1000 ultrafast bootstrap replicates [76]. Bayesian inference (BI) analyses were performed with MrBayes version 3.2.6 [77] with the partitioned models (Additional file 15: Table S8) determined by PartitionFinder version 2.1.1 [78]. 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 [79] and the potential scale reduction factor (PSRF) approached 1.0 [77]. 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

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 [80]. 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 [83].

Availability of data and materials

The data that support the findings of this study are openly available in GenBank of NCBI at, with the reference numbers MH013482, MW936632–MW936633, MW924656–MW924665.



Mitochondrial genome

atp6 and atp8 :

ATPase subunits 6 and 8

cob :

Cytochrome b

cox1cox3 :

Cytochrome c oxidase subunits I-III

nad1nad6 and nad4l :

NADH dehydrogenase subunits 1–6 and 4 L


Protein-coding gene

rrnL and rrnS :

Large and small subunit ribosomal RNA (rRNA)


Transfer RNA


Relative synonymous codon usage


Bayesian inference


Maximum likelihood


  1. Curole JP, Kocher TD. Mitogenomics: digging deeper with complete mitochondrial genomes. Trends Ecol Evol. 1999;14(10):394–8.

    Article  CAS  PubMed  Google Scholar 

  2. Cameron SL. Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014;59(1):95–117.

    Article  CAS  PubMed  Google Scholar 

  3. Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27(8):1767–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Timmermans MJTN, Lees DC, Simonsen TJ. Towards a mitogenomic phylogeny of Lepidoptera. Mol Phylogenet Evol. 2014;79:169–78.

    Article  PubMed  Google Scholar 

  5. 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.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Song N, Zhang H, Zhao T. Insights into the phylogeny of Hemiptera from increased mitogenomic taxon sampling. Mol Phylogenet Evol. 2019;137:239–49.

    Article  Google Scholar 

  8. 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.

    Article  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. Mitter C, Davis DR, Cummings MP. Phylogeny and evolution of Lepidoptera. Annu Rev Entomol. 2017;62(1):265–83.

    Article  CAS  PubMed  Google Scholar 

  11. Triant DA, Cinel SD, Kawahara AY. Lepidoptera genomes: current knowledge, gaps and future directions. Curr Opin Insect Sci. 2018;25:99–105.

    Article  PubMed  Google Scholar 

  12. 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.

    Article  PubMed  Google Scholar 

  13. 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.

    Article  CAS  PubMed  Google Scholar 

  14. 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.

    Article  Google Scholar 

  15. 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.

    Article  CAS  PubMed  Google Scholar 

  16. 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.

    Article  Google Scholar 

  17. 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.

    Article  PubMed  Google Scholar 

  18. Gilligan TM, Baixeras J, Brown JW. T@RTS: Online World Catalogue of the Tortricidae (Version 4.0). 2018.

    Google Scholar 

  19. 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.

    Google Scholar 

  20. Horak M. Olethreutine moths of Australia (Lepidoptera: Tortricidae), Monographs on Australian Lepidoptera; 2006.

    Book  Google Scholar 

  21. Brown JW. World catalogue of insects, Vol. 5: Tortricidae (Lepidoptera). Stenstrup: Apollo Books; 2005.

    Book  Google Scholar 

  22. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  PubMed  Google Scholar 

  24. van der Geest LPS, Evinhuis HH. Tortricid pests, their biology, natural enemies and control. Amsterdam: Elsevier Science Publishers B.V; 1991.

    Google Scholar 

  25. 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.

    Article  Google Scholar 

  26. 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.

    Article  CAS  PubMed  Google Scholar 

  27. 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.

    Article  Google Scholar 

  28. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 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.

    Article  PubMed  Google Scholar 

  30. 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.

    Google Scholar 

  31. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    CAS  Google Scholar 

  36. 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.

    Article  CAS  PubMed  Google Scholar 

  37. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  CAS  PubMed  Google Scholar 

  40. 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.

    Article  Google Scholar 

  41. Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41(3):353–8.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. 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.

    Article  CAS  PubMed  Google Scholar 

  44. Hershberg R, Petrov DA. Selection on codon Bias. Annu Rev Genet. 2008;42(1):287–99.

    Article  CAS  PubMed  Google Scholar 

  45. Plotkin JB, Kudla G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2011;12(1):32–42.

    Article  CAS  PubMed  Google Scholar 

  46. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. 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.

    Article  CAS  PubMed  Google Scholar 

  49. 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.

    Article  CAS  PubMed  Google Scholar 

  50. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    Article  Google Scholar 

  56. 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.

    Article  PubMed  Google Scholar 

  57. 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.

    Article  Google Scholar 

  58. 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.

    Article  PubMed  Google Scholar 

  59. Li H. Microlepidoptera of Qinling Moutains (Insecta: Lepidoptera). Beijing: Science Press; 2012.

    Google Scholar 

  60. Liu YQ, Li GW. Fauna Sinica. Vol. 27: Lepidoptera: Tortricidae. Beijing: Science Press; 2002.

    Google Scholar 

  61. 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.

    Google Scholar 

  62. Schubert M, Lindgreen S, Orlando L. AdapterRemoval v2: rapid adapter trimming, identification, and read merging, BMC res. Notes. 2016;9:88.

    Google Scholar 

  63. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  65. 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.

    CAS  PubMed  Google Scholar 

  66. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. 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.

    Article  PubMed  Google Scholar 

  68. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. 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.

    Article  CAS  PubMed  Google Scholar 

  70. Peden JF. Analysis of codon usage. Univ Nottingham. 2000;90:73–4.

    Google Scholar 

  71. 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.

    Article  CAS  Google Scholar 

  72. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6.

    Article  CAS  PubMed  Google Scholar 

  73. 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.

    Article  PubMed  Google Scholar 

  74. 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.

    Article  CAS  PubMed  Google Scholar 

  75. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. 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.

    Article  CAS  PubMed  Google Scholar 

  77. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  78. 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.

    Article  CAS  PubMed  Google Scholar 

  79. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. 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.

    Article  CAS  PubMed  Google Scholar 

  81. Townsend JP. Profiling phylogenetic informativeness. Syst Biol. 2007;56(2):222–31.

    Article  CAS  PubMed  Google Scholar 

  82. Lopez-Giraldez F, Townsend JP. PhyDesign: an online application for profiling phylogenetic informativeness. BMC Evol Biol. 2011;11(1):152.

    Article  PubMed  PubMed Central  Google Scholar 

  83. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


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).

Author information

Authors and Affiliations



MY designed the research, performed the data analysis and wrote the manuscript. JL, SS and WD analyzed the data. HZ and ZW collected the samples. LL designed the research and improved the manuscript writing. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Mingsheng Yang or Lili Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

The tortricid samples used in phylogenetic analyses.

Additional file 2: Table S2.

Annotation and comparison of mitochondrial genome organizations of 13 mitogenomes sequenced in this study.

Additional file 3: Table S3.

Reannotation and comparison of mitochondrial genome organizations of 27 previous sequenced tortricid species.

Additional file 4: Table S4.

GC-Skew of the mitogenomes in six tortricid tribes.

Additional file 5: Table S5.

Codon usage pattern in six tortricid tribes analyzed in this study.

Additional file 6: Figure S1.

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.

Additional file 7: Figure S2.

Phylogenetic tree inferred from maximum likelihood method based on P12 dataset. Numbers on nodes represent the bootstrap supports.

Additional file 8: Figure S3.

Phylogenetic tree inferred from Bayesian inference method based on P12 dataset. Numbers on nodes represent the posterior probabilities.

Additional file 9: Figure S4.

Phylogenetic tree inferred from maximum likelihood method based on P123 dataset. Numbers on nodes represent the bootstrap supports.

Additional file 10: Figure S5.

Phylogenetic tree inferred from Bayesian inference method based on P123 dataset. Numbers on nodes represent the posterior probabilities.

Additional file 11: Figure S6.

Phylogenetic tree inferred from maximum likelihood method based on P123R dataset. Numbers on nodes represent the bootstrap supports.

Additional file 12: Figure S7.

Phylogenetic tree inferred from Bayesian inference method based on P123R dataset. Numbers on nodes represent the posterior probabilities.

Additional file 13: Table S6.

Information of samples sequenced in this study.

Additional file 14: Table S7.

The partitioning schemes and corresponding substitution models determined by ModelFinder.

Additional file 15: Table S8.

The partitioning schemes and corresponding substitution models determined by PartitionFinder.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: