The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

The Complete Mitochondrial Genome and Phylogenetic Analysis of Rhagastis binoculata (Matsumura, 1909) (Lepidoptera: Sphingidae)

The mitochondrial genome (mitogenome) Rhagastis binoculata (Matsumura, 1909), an endemic moth species in Taiwan, was sequenced and analyzed. The complete circular mitogenome of R. binoculata is 15,303 bp and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region. The mitogenome has an overall nucleotide composition of 41.2% A, 11.9% C, 7.5% G, and 39.4% T, with an AT content of 80.6%. Of the protein-coding genes (PCGs), 12 start with ATG, ATT, and ATC, and COX1 starts with a "CGA" codon. All of the stop codons are "TAA, TAG, or T". Our phylogenetic analysis of 21 species of Sphingidae insects suggests that R. binoculata is clustered with Rhagastis mongoliana, which belongs to the subfamily Macroglossinae.

The complete mitochondrial genome of Chibiraga houshuaii (Lepidoptera, Limacodidae) and its phylogenetic implications

Chibiraga is a mall East Asian genus in the family Limacodidae (slug-moths). The latter includes many agricultural pests. Mitochondrial genome analysis is an important tool for studying insect molecular identification and phylogenetics. However, there are very few mitogenome sequences available for Limacodidae species, and none for the genus Chibiraga at all. To explore the mitogenome features of Chibiraga and verify its phylogenetic position, the complete mitogenome of Chibiraga houshuaii was sequenced and annotated. The complete 15,487 bp genome encoded 37 mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region (CR). Most of the PCGs had typical ATN start codons and terminated with TAA or a single T residue. UUA (Leu2), AUU (Ile), UUU (Phe), AUA (Met) and AAU (Asn) were the five most frequently used codons. All tRNAs were folded into cloverleaf secondary structure, except for trnS1, which lacked the DHU arm. Phylogenetic analyses within the superfamily Zygaenoidea were performed based on multiple datasets from mitochondrial genes. The results showed that the families Phaudidae, Limacodidae and Zygaenidae were respectively recovered as monophyly; C. houshuaii was clustered in a clade with nettle type larvae in Limacodidae.

The Complete Mitochondrial Genome of Chilo infuscatellus (Lepidoptera: Pyralidae), and Related Phylogenetic Analysis

The Chilo infuscatellus (Lepidoptera: Pyralidae) is a significant pest of sugarcane in China. The genome-level characteristics of this pest are important genetic resources for identification, phylogenetic analysis, and even management. In the present study, the complete mitogenome of C. infuscatellus was sequenced and characterized. The assembled mitochondrial genome is 15,252 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and an A + T-rich region. Except for the CGA codon for the cox1 gene, the PCGs are initiated with ATN codons (ATG, ATT, and ATA). These PCGs are terminated with TAA or an incomplete termination codon of a single T. Except for the loss of the "DHU" arm for trnS1, the tRNA genes were folded into the typical cloverleaf structure. The A + T-rich region has a high AT content of 96.19% and contains the motifs "ATAGA" and "ATTTA", as well as a 19 bp poly-T stretch and microsatellite regions. The C. infuscatellus mitogenome exhibits a conserved gene order among lepidopteran insects, with a rearrangement of the trnM gene compared to the ancestral insect gene order. Phylogenetic analysis based on the 13 PCGs using Bayesian inference (BI) and maximum likelihood (ML) methods confirmed the monophyly of Pyralidae and Crambidae within Pyraloidea. The relationships between subfamilies in Pyralidae can be described as (Galleriinae + (Phycitinae + (Pyralinae + Epipaschiinae))). The "PS clade" and "non-PS clade" were formed within the family Crambidae. These findings provide valuable genetic resources for the identification, phylogenetic analysis, and management of sugarcane borers, contributing significantly to our understanding of the phylogeny of Pyraloidea insects and their evolution.

Life history and mitochondrial genomes of Salassinae and Agliinae (Insecta, Lepidoptera): New insights into the loss of cocooning behaviour and phylogeny of Saturniidae

The subfamilies Salassinae and Agliinae are two monogeneric groups of the family Saturniidae. They were regarded as the non-cocooning saturniids in Asia. Since very little information on their life history and mitogenome has been reported, their origin and evolution are still poorly understood. In this study, nature-imitated rearing is used to record the life history of two Aglia and five Salassa species. In addition, four complete mitogenomes are presented, which are the first ones of these two subfamilies. The results show that both Salassinae and Agliinae have lost their cocooning. Moreover, the phylogenetic analysis demonstrates that the subfamily Saturniinae is not monophyletic due to the inclusion of Agliinae and Salassinae.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.

Variants in the mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) infected with Oryctes rhinoceros nudivirus in oil palm and coconut plantations

The CRB (coconut rhinoceros beetle) haplotype was classified into CRB-S and CRB-G, based on the presence of single nucleotide polymorphisms (SNPs) in the mitochondrial cox1 gene. Mitochondrial genomes (mitogenomes) are the most widely used genetic resources for molecular evolution, phylogenetics, and population genetics in relation to insects. This study presents the mitogenome CRB-G and CRB-S which were collected in Johor, Malaysia. The mitogenome of CRB-G collected from oil palm plantations in 2020 and 2021, and wild coconut palms in 2021 was 15,315 bp, 15,475 bp, and 17,275 bp, respectively. The CRB-S was discovered in coconut and oil palms in 2021, and its mitogenome was 15,484 bp and 17,142 bp, respectively. All the mitogenomes have 37 genes with more than 99% nucleotide sequence homology, except the CRB-G haplotype collected from oil palm in 2021 with 89.24% nucleotide sequence homology. The mitogenome of Johor CRBs was variable in the natural population due to its elevated mutation rate. Substitutions and indels in cox1, cox2, nad2 and atp6 genes were able to distinguish the Johor CRBs into two haplotypes. The mitogenome data generated in the present study may provide baseline information to study the infection and relationship between the two haplotypes of Johor CRB and OrNV in the field. This study is the first report on the mitogenomes of mixed haplotypes of CRB in the field.

Determining the complete mitochondrial genome of Tethea albicostata (Lepidoptera: Drepanidae) and phylogenetic analysis

Tethea albicostata is a widely distributed insect species in northern and central China. To date, few studies have been conducted on this species, with the exception of morphological taxonomy studies. Here, we report the complete mitochondrial genome of T. albicostata collected in China. The circular-mapping mitogenome is 15,308 bp in length, with an overall A + T content of 80.52%, encoding 2 ribosomal RNA genes, 22 transfer RNA genes, and 13 protein-coding genes. The gene arrangement and components of T. albicostata are identical to those of most other Lepidopteran insects. Phylogenetic analysis based on mitogenomes showed that T. albicostata is grouped with Drepana pallida, which belongs to the same family as Drepanidae. The family Drepanidae formed a separate branch from other families in the phylogenetic tree. This study determined the second mitochondrial genome of the Drepanidae species.

Mitogenomes of Eight Nymphalidae Butterfly Species and Reconstructed Phylogeny of Nymphalidae (Nymphalidae: Lepidoptera)

The Nymphalidae family of cosmopolitan butterflies (Lepidoptera) comprises approximately 7200 species found on all continents and in all habitats. However, debate persists regarding the phylogenetic relationships within this family. In this study, we assembled and annotated eight mitogenomes of Nymphalidae, constituting the first report of complete mitogenomes for this family. Comparative analysis of 105 mitochondrial genomes revealed that the gene compositions and orders were identical to the ancestral insect mitogenome, except for Callerebia polyphemus trnV being before trnL and Limenitis homeyeri having two trnL genes. The results regarding length variation, AT bias, and codon usage were consistent with previous reports on butterfly mitogenomes. Our analysis indicated that the subfamilies Limenitinae, Nymphalinae, Apaturinae, Satyrinae, Charaxinae, Heliconiinae, and Danainae are monophyletic, while the subfamily the subfamily Cyrestinae is polyphyletic. Danainae is the base of the phylogenetic tree. At the tribe level, Euthaliini in Limenitinae; Melitaeini and Kallimini in Nymphalinae; Pseudergolini in Cyrestinae; Mycalesini, Coenonymphini, Ypthimini, Satyrini, and Melanitini in Satyrinae; and Charaxini in Charaxinae are regarded as monophyletic groups. However, the tribe Lethini in Satyrinae is paraphyletic, while the tribes Limenitini and Neptini in Limenitinae, Nymphalini and Hypolimni in Nymphalinae, and Danaini and Euploeini in Danainae are polyphyletic. This study is the first to report the gene features and phylogenetic relationships of the Nymphalidae family based on mitogenome analysis, providing a foundation for future studies of population genetics and phylogenetic relationships within this family.

Comparative genome and phylogenetic analysis revealed the complex mitochondrial genome and phylogenetic position of Conopomorpha sinensis Bradley

Conopomorpha sinensis Bradley is a destructive pest that causes severe economic damage to litchi and longan. Previous C. sinensis research has focused on population life tables, oviposition selectivity, pest population prediction, and control technology. However, there are few studies on its mitogenome and phylogenetic evolution. In this study, we sequenced the whole mitogenome of C. sinensis by the third-generation sequencing, and analyzed the characteristics of its mitogenome by comparative genome. The complete mitogenome of C. sinensis is a typical circular and double-stranded structure. The ENC-plot analyses revealed that natural selection could affect the information of codon bias of the protein-coding genes in the mitogenome of C. sinensis in the evolutionary process. Compared with 12 other Tineoidea species, the trnA-trnF gene cluster of tRNA in the C. sinensis mitogenome appears to have a new arrangement pattern. This new arrangement has not been found in other Tineoidea or other Lepidoptera, which needs further exploration. Meanwhile, a long AT repeated sequence was inserted between trnR and trnA, trnE and trnF, ND1 and trnS in the mitogenome of C. sinensis, and the reason for this sequence remains to be further studied. Furthermore, the results of phylogenetic analysis showed that the litchi fruit borer belonged to Gracillariidae, and Gracillariidae was monophyletic. The results will contribute to an improved understanding of the complex mitogenome and phylogeny of C. sinensis. It also will provide a molecular basis for further research on the genetic diversity and population differentiation of C. sinensis.

The Complete Mitochondrial Genome of Box Tree Moth Cydalima perspectalis and Insights into Phylogenetics in Pyraloidea

To resolve and reconstruct phylogenetic relationships within Pyraloidea based on molecular data, the mitochondrial genome (mitogenome) was widely applied to understand phylogenetic relations at different taxonomic levels. In this research, a complete mitogenome of Cydalima perspectalis was recorded, and the phylogenetic position of C. perspectalis was inferred based on the sequence in combination with other available sequence data. According to the research, the circular mitochondrial genome is 15,180 bp in length. It contains 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), 13 typical protein-coding genes (PCGs), and a non-coding control region. The arrangement of a gene of the C. perspectalis mitogenome is not the same as the putative ancestral arthropod mitogenome. All of the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is undertaken by CGA. Five genes have incomplete stop codons that contain only ‘T’. All tRNA genes display a typical clover–leaf structure of mitochondrial tRNA, except for trnS1 (AGN). The control region contained an ‘ATAGG(A)’-like motif followed by a poly-T stretch. Based on the mitochondrial data, phylogenetic analysis within Pyraloidea was carried out using Bayesian inference (BI) and maximum likelihood (ML) analyses. Phylogenetic analysis showed that C. perspectalis is more closely related to Pygospila tyres within Spilomelinae than those of Crambidae and Pyraloidea.

Evolution of mitogenomic gene order in Orthoptera

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

The mitochondrial genomes of the Geometroidea (Lepidoptera) and their phylogenetic implications

The Geometroidea is a large superfamily of Lepidoptera in species composition and contains numerous economically important pest species that cause great loss in crop and forest production. However, understanding of mitogenomes remains limited due to relatively fewer mitogenomes previously reported for this megadiverse group. Here, we sequenced and annotated nine mitogenomes for Geometridae and further analyzed the mitogenomic evolution and phylogeny of the whole superfamily. All nine mitogenomes contained 37 mitochondrial genes typical in insects, and gene organization was conserved except for Somatina indicataria. In S. indicataria, the positions of two tRNAs were rearranged. The trnR was located before trnA instead of after trnA typical in Lepidoptera, whereas the trnE was detected rarely on the minority strand (N-strand). This trnR-trnA-trnN-trnS1-trnE-trnF newly recognized in S. indicataria represents the first gene rearrangement reported for Geometroidea and is also unique in Lepidoptera. Besides, nucleotide composition analyses showed little heterogeneity among the four geometrid subfamilies involved herein, and overall, nad6 and atp8 have higher nucleotide diversity and Ka/Ks rate in Geometridae. In addition, the taxonomic assignments of the nine species, historically defined by morphological studies, were confirmed by various phylogenetic analyses based on the hitherto most extensive mitogenomic sampling in Geometroidea.

Evolution of Gene Arrangements in the Mitogenomes of Ensifera and Characterization of the Complete Mitogenome of Schizodactylus jimo

Gene arrangement (relative location of genes) is another evolutionary marker of the mitogenome that can provide extensive information on the evolutionary mechanism. To explore the evolution of gene arrangements in the mitogenome of diversified Ensifera, we sequenced the mitogenome of the unique dune cricket species found in China and used it for phylogenetic analysis, in combination with 84 known Ensiferan mitogenomes. The mitogenome of Schizodactylus jimo is a 16,428-bp circular molecule that contains 37 genes. We identified eight types of gene arrangement in the 85 ensiferan mitogenomes. The gene location changes (i.e., gene translocation and duplication) were in three gene blocks: I-Q-M-ND2, rrnl-rns-V, and ND3-A-R-N-S-E-F. From the phylogenetic tree, we found that Schizodactylus jimo and most other species share a typical and ancient gene arrangement type (Type I), while Grylloidea has two types (Types II and III), and the other five types are rare and scattered in the phylogenetic tree. We deduced that the tandem replication–random loss model is the evolutionary mechanism of gene arrangements in Ensifera. Selection pressure analysis revealed that purifying selection dominated the evolution of the ensiferan mitochondrial genome. This study suggests that most gene rearrangements in the ensiferan mitogenome are rare accidental events.

Enhanced Resolution of Evolution and Phylogeny of the Moths Inferred from Nineteen Mitochondrial Genomes

The vast majority (approximately 90%) of Lepidoptera species belong to moths whose phylogeny has been widely discussed and highly controversial. For the further understanding of phylogenetic relationships of moths, nineteen nearly complete mitochondrial genomes (mitogenomes) of moths involved in six major lineages were sequenced and characterized. These mitogenomes ranged from 15,177 bp (Cyclidia fractifasciata) to 15,749 bp (Ophthalmitis albosignaria) in length, comprising of the core 37 mitochondrial genes (13 protein-coding genes (PCGs) + 22 tRNAs + two rRNAs) and an incomplete control region. The order and orientation of genes showed the same pattern and the gene order of trnM-trnI-trnQ showed a typical rearrangement of Lepidoptera compared with the ancestral order of trnI-trnQ-trnM. Among these 13 PCGs, ATP8 exhibited the fastest evolutionary rate, and Drepanidae showed the highest average evolutionary rate among six families involved in 66 species. The phylogenetic analyses based on the dataset of 13 PCGs suggested the relationship of (Notodontidae + (Noctuidae + Erebidae)) + (Geometridae + (Sphingidae + Drepanidae)), which suggested a slightly different pattern from previous studies. Most groups were well defined in the subfamily level except Erebidae, which was not fully consistent across bayesian and maximum likelihood methods. Several formerly unassigned tribes of Geometridae were suggested based on mitogenome sequences despite a not very strong support in partial nodes. The study of mitogenomes of these moths can provide fundamental information of mitogenome architecture, and the phylogenetic position of moths, and contributes to further phylogeographical studies and the biological control of pests.

The first complete mitochondrial genome in the family Attevidae (Attevaaurea) of the order Lepidoptera

The superfamily Yponomeutoidea, one of the early-derived groups in the order Lepidoptera, consists of 11 families. However, mitochondrial genome (mitogenome) sequences, popularly used for phylogeny and evolutionary tracing, are available for only seven species across six genera and five families. Thus, a larger variety of mitogenome sequences in Yponomeutoidea are required to improve our understanding of lepidopteran phylogeny and genomic evolution. In this study, we present the complete mitogenome of Attevaaurea (Fitch, 1856), the first species in the family Attevidae (superfamily Yponomeutoidea, order Lepidoptera) to be sequenced. The complete mitogenome comprises 16,329 bp and contains a typical set of genes and one non-coding region. Within Yponomeutoidea, the mitogenome of A.aurea has a unique trnI-trnM-trnQ arrangement at the A + T-rich region and ND2 junction and trnA-ND3 arrangement at the trnG and trnR junction. Twelve of the 13 protein-coding genes (PCGs) of A.aurea have a typical ATN starting codon, whereas COI has the atypical CGA codon, which is frequently found in the starting region of lepidopteran COI. Phylogenetic analyses, based on the concatenated sequences of 13 PCGs and two rRNA genes, using the Maximum Likelihood method, revealed a sister relationship between Attevidae and Praydidae with moderately low nodal support (bootstrap support = 64%).

Comparative mitochondrial genome analysis and phylogenetic relationship among lepidopteran species

Lepidoptera has rich species including many agricultural pests and economical insects around the world. The mitochondrial genomes (mitogenomes) were utilized to explore the phylogenetic relationships between difference taxonomic levels in Lepidoptera. However, the knowledge of mitogenomic characteristics and phylogenetic position about superfamily-level in this order is unresolved. In this study, we integrated 794 mitogenomes consisting of 37 genes and a noncoding control region, which covered 26 lepidopteran superfamilies from newly sequenced and publicly available genomes for comparative genomic and phylogenetic analysis. In primitive taxon, putative start codon of cox1 gene was ATA or ATT instead of CGA, but stop codon of that showed four types, namely TAA, TAG, TA and T. The 7-bp overlap between atp8 and atp6 presented as "ATGATAA". Moreover, the most frequently utilized amino acids were leucine (UUA) in 13 PCGs. Phylogenetic analysis showed that the main backbone relationship in Lepidoptera was (Hepialoidea + (Nepticuloidea + (Adeloidea + (Tischerioidea + (Tineoidea + (Yponomeutoidea + (Gracillarioidea + (Papilionoidea + ((Zygaenoidea + Tortricoidea) + (Gelechioidea + (Pyraloidea + ((Geometroidea + Noctuoidea) + (Lasiocampoidea + Bombycoidea))))))))))))).

Mitogenomics and phylogenetics of twelve species of African Saturniidae (Lepidoptera)

African Saturniidae (Lepidoptera) include numerous species consumed at the caterpillar stage throughout the continent, and their importance to local communities as a source of nutrition and seasonal income cannot be overestimated. However, baseline genetic data with utility for the characterization of their diversity, phylogeography and phylogenetic relationships have remained scarce compared to their Asian counterparts. To bridge this gap, we sequenced the mitochondrial genomes of 12 species found in southern Africa for comparative mitogenomics and phylogenetic reconstruction of the family, including the first representatives of the tribes Eochroini and Micragonini. Mitochondrial gene content and organization were conserved across all Saturniidae included in the analyses. The phylogenetic positions of the 12 species were assessed in the context of publicly available mitogenomes using Bayesian inference and maximum likelihood (ML) methods. The monophyly of the tribes Saturniini, Attacini, Bunaeini and Micragonini, the sister relationship between Saturniini and Attacini, and the placement of Eochroa trimenii and Rhodinia fugax in the tribes Eochroini and Attacini, respectively, were strongly supported. These results contribute to significantly expanding genetic data available for African Saturniidae and allow for the development of new mitochondrial markers in future studies.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

Comparative mitogenomic analysis of two earwigs (Insecta, Dermaptera) and the preliminary phylogenetic implications

The phylogenetic position and inner relationships of Dermaptera remain unresolved despite the numerous efforts using morphological and molecular data. To facilitate the resolution of problems, this study sequenced the complete mitogenome of Apachyusfeae de Bormans, 1894 (Apachyidae) and the nearly complete mitogenome of Diplatysflavicollis Shiraki, 1907 (Diplatyidae). The 19,029-bp long mitogenome of A.feae exhibited an extra trnV gene and two control regions in addition to the typical set of 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and two ribosomal RNA (rRNA) genes. The 12,950-bp long partially sequenced mitogenome of D.flavicollis was composed of 10 and a partial fragment of PCGs, 18 tRNA genes, two rRNA genes, and a control region. Comparative analysis of available earwig mitogenomes revealed variable mitogenomic structure and extensive gene rearrangements in Dermaptera. The preliminary phylogenetic analyses using Bayesian inference and maximum likelihood methods showed identical results, but the limited sampling and different types of molecular data lead to an apparent incongruence with previous phylogenetic studies.

The complete mitochondrial genome data of the Common Rose butterfly, Pachliopta aristolochiae (Lepidoptera, Papilionoidea, Papilionidae) from Malaysia

Here, we present the complete mitochondrial genome of Pachliopta aristolochiae, a Common Rose butterfly from Malaysia. The sequence was generated using Illumina NovaSeq 6000 sequencing platform. The mitogenome is 15,235bp long, consisting of 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and two D-loop regions. The total base composition was (81.6%), with A (39.3%), T (42.3%), C (11.0%) and G (7.3%). The gene order of the three tRNAs was trnM-trnI-trnQ, which differs from the ancestral insect gene order trnI-trnQ-trnM. Phylogenetic tree analysis revealed that the sequenced Pachliopta aristolochiae in this data is closely related to Losaria neptunus (NC 037868), with highly supported ML and BI analysis. The data presented in this work can provide useful resources for other researchers to study deeper into the phylogenetic relationships of Lepidoptera and the diversification of the Pachliopta species. Also, as one of the bioindicator species, this data can be used to assess environmental changes in the terrestrial and aquatic ecosystem via enviromental DNA approahes. The mitogenome of Pachliopta aristolochiae is available in GenBank under the accession number MZ781228.

Complete mitochondrial genome of the recently discovered multivoltine Graphium (Pazala) confucius Hu, Duan & Cotton, 2018 (Lepidoptera: Papilionidae)

Graphium (Pazala) confucius Hu, Duan & Cotton, 2018 is a recently discovered, wide ranging, multivoltine swordtail butterfly in China and Vietnam. The present study reports the complete mitochondrial genome of this butterfly, which is the fifth mitochondrial genome record for subgenus Pazala Moore, 1888. The mitochondrial genome of G. (P.) confucius is circular and 15,212 bp in length, and consists of 37 genes, including 13 PCGs, 22 tRNAs, and two rRNAs. The Bayesian phylogenetic tree containing the focal species and 33 other Papilioninae members clusters G. (P.) confucius with other Pazala taxa inside tribe Leptocircini, which agrees with its taxonomic position. The findings of this study added data to the complex subgenus Pazala and are beneficial to future understanding and conservation planning of butterfly diversity.

Interspecific Hybridization and Complete Mitochondrial Genome Analysis of Two Ghost Moth Species

Simple Summary

The Chinese cordyceps is a valuable parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex. In view of culturing this complex, a method for the artificial rearing of the Thitarodes/Hepialus ghost moth hosts was established. Deterioration of the host insect population and low mummification rates in infected larvae constrain effective cultivation. Hybridization of Thitarodes/Hepialus populations may overcome this problem. Thitarodes shambalaensis and Thitarodes sp. were inbred or hybridized, and the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations were investigated. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. One hybrid population (T. shambalaensis females mated with Thitarodes sp. males) showed increased population growth as compared with the parental Thitarodes sp. population. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis differed. The complete mitochondrial genomes of T. shambalaensis, Thitarodes sp. and the hybrid population were 15,612 bp, 15,389 bp and 15,496 bp in length, respectively. A + T-rich regions were variable in sizes and repetitive sequences. The hybrid population was located in the same clade with T. shambalaensis, implying the maternal inheritance of mitochondrial DNA.

Abstract

The Chinese cordyceps, a parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex, is a valuable biological resource endemic to the Tibetan Plateau. Protection of the Plateau environment and huge market demand make it necessary to culture this complex in an artificial system. A method for the large-scale artificial rearing of the Thitarodes/Hepialus insect host has been established. However, the deterioration of the insect rearing population and low mummification of the infected larvae by the fungus constrain effective commercial cultivation. Hybridization of Thitarodes/Hepialus populations may be needed to overcome this problem. The species T. shambalaensis (GG♂ × GG♀) and an undescribed Thitarodes species (SD♂ × SD♀) were inbred or hybridized to evaluate the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations. The two parental Thitarodes species exhibited significant differences in adult fresh weights and body lengths but not in pupal emergence rates. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. The SD♂ × GG♀ population showed a higher population trend index than the SD♂ × SD♀ population, implying increased population growth compared with the male parent. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis also differed. This provides possibilities to create Thitarodes/Hepialus populations with increased growth potential for the improved artificial production of the insect hosts. The mitochondrial genomes of GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀ were 15,612 bp, 15,389 bp and 15,496 bp in length, with an A + T content of 80.92%, 82.35% and 80.87%, respectively. The A + T-rich region contains 787 bp with two 114 bp repetitive sequences, 554 bp without repetitive sequences and 673 bp without repetitive sequences in GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀, respectively. The hybrid population (SD♂ × GG♀) was located in the same clade with GG♂ × GG♀, based on the phylogenetic tree constructed by 13 PCGs, implying the maternal inheritance of mitochondrial DNA.

Comparative Analysis of Eight Mitogenomes of Bark Beetles and Their Phylogenetic Implications

Simple Summary

Many bark beetles are destructive pests in coniferous forests and cause extensive ecological and economic losses worldwide. Comparative studies of the structural characteristics of mitogenomes and phylogenetic relationships of bark beetles can improve our understanding of mitogenome evolution. In this study, we sequenced eight mitogenomes of bark beetles. Our results show that the use of start and stop codons, the abundance of amino acids, and the relative frequency of codon use are conserved among the eight bark beetles. Different regions of tRNA exhibit different degrees of conservatism. Together with the analysis of evolutionary rates and genetic distance among bark beetle species, our results reveal phylogenetic relationships among bark beetles of the subfamily Scolytinae.

Abstract

Many bark beetles of the subfamily Scolytinae are the most economically important insect pests of coniferous forests worldwide. In this study, we sequenced the mitochondrial genomes of eight bark beetle species, including Dendroctonus micans, Orthotomicus erosus, Polygraphus poligraphus, Dryocoetes hectographus, Ips nitidus, Ips typographus, Ips subelongatus, and Ips hauseri, to examine their structural characteristics and determine their phylogenetic relationships. We also used previously published mitochondrial genome sequence data from other Scolytinae species to identify and localize the eight species studied within the bark beetle phylogeny. Their gene arrangement matched the presumed ancestral pattern of these bark beetles. Start and stop codon usage, amino acid abundance, and the relative codon usage frequencies were conserved among bark beetles. Genetic distances between species ranged from 0.037 to 0.418, and evolutionary rates of protein-coding genes ranged from 0.07 for COI to 0.69 for ND2. Our results shed light on the phylogenetic relationships and taxonomic status of several bark beetles in the subfamily Scolytinae and highlight the need for further sequencing analyses and taxonomic revisions in additional bark beetle species.

The complete mitochondrial genome of Dudusa sphingiformis (Lepidoptera: Notodontidae) and phylogenetic analysis

Dudusa sphingiformis is an important lepidopteran pest widely distributed in tropical and subtropical zones of Asia. In this paper, the complete mitochondrial genome (mitogenome) of D. sphingiformis was determined by next-generation sequencing. The mitogenome was 15,806 bp in length, comprising 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and an AT-rich control region (D-loop). The gene arrangement of this mitogenome was identical to that of the previous studies of Notodontidae moths. Almost all the PCGs initiated with typical ATN codons, except for cox1 with CGA. Among them, nine PCGs terminated with TAA or TAG, while other four PCGs (cox1, cox2, nad5, and nad4) with incomplete stop codon T. All the 22 tRNAs had the typical cloverleaf structure, except for trnS1, whose dihydrouridine (DHU) arm forms a simple loop. Phylogenetic analysis based on the concatenated nucleotide sequences of 13 PCGs indicated that D. sphingiformis was more closely related to other species of family Notodontidae, forming a monophyletic group, with well-resolved relationships among five family of Noctuoidea.

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Saturnia japonica

The complete mitochondrial genome (mitogenome) of Saturnia japonica (Lepidoptera: Saturniidae) was sequenced and annotated. It is a circular molecule of 15, 376 bp, composed of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNA), and an adenine (A) + thymine (T)-rich region. All protein-coding genes (PCGs) are initiated by the ATN codon except for cytochrome c oxidase subunit 1 (cox1) gene that is seemingly initiated by the CGA codon. Except for cox2 and nad4, which were terminated by incomplete stop codon T or TA, the rest were terminated by canonical stop codon TAA. The A + T-rich region is high conservative, including 'ATAGA' motif followed by a 19 bp poly-T stretch, a microsatellite-like element (AT)₉ and also a poly-A element, with a total length of 332 bp. The Asn codon was the most frequently used codon, followed by Ile, Leu2, Lys, Met, Phe, and Tyr, while Cys was the least frequently used codon. Phylogenetic relationships analysis based on the 13 PCGs by using maximum likelihood (ML) and neighbor Joining (NJ) revealed that S. japonica belongs to the Saturniidae family. In this study, the annotation and characteristics of the mitogenome of S. japonica were resolved for the first time, which laid a foundation for species classification and the molecular evolution of Lepidoptera: Saturniidae.

The complete mitochondrial genome of Periacma orthiodes Meyrick, 1894 (Lepidoptera: Autostichidae)

The complete mitogenome of Periacma orthiodes Meyrick, 1894 was determined in this study. It was 15,306 bps long and strongly AT biased. It consisted of 13 PCGs, 22 tRNAs, 2 rRNAs and 1 non-coding control region (371 bps). Most PCGs used the typical ATN start codon, except for cox1. Four genes (cox1, cox2, nad4 and nad5) used truncated stop codons (a single T or TA) rather than the commonly used TAA or TAG. All tRNAs, excluding TrnS1, folded into the iconic cloverleaf structure. ML phylogenetic tree built on 13 PCGs from P. orthiodes and another 28 species in Gelechioidea demonstrated that the genus Periacma was a member of the family Autostichidae, which was consistent with the newest phylogenetic study.

Mitochondrial Genomes of Hestina persimilis and Hestinalis nama (Lepidoptera, Nymphalidae): Genome Description and Phylogenetic Implications

Simple Summary

In this study, the mitogenomes of Hestina persimilis and Hestinalis nama were obtained via sanger sequencing. Compared with other mitogenomes of Apaturinae butterflies, conclusions can be made that the mitogenomes of Hestina persimilis and Hestinalis nama are highly conservative. The phylogenetic trees build upon mitogenomic data showing that the relationships among Nymphalidae are similar to previous studies. Hestinalisnama is apart from Hestina, and closely related to Apatura, forming a monophyletic clade.

Abstract

In this study, the complete mitochondrial genomes (mitogenomes) of Hestina persimilis and Hestinalis nama (Nymphalidae: Apaturinae) were acquired. The mitogenomes of H. persimilis and H. nama are 15,252 bp and 15,208 bp in length, respectively. These two mitogenomes have the typical composition, including 37 genes and a control region. The start codons of the protein-coding genes (PCGs) in the two mitogenomes are the typical codon pattern ATN, except CGA in the cox1 gene. Twenty-one tRNA genes show a typical clover leaf structure, however, trnS1(AGN) lacks the dihydrouridine (DHU) stem. The secondary structures of rrnL and rrnS of two species were predicted, and there are several new stem loops near the 5′ of rrnL secondary structure. Based on comparative genomic analysis, four similar conservative structures can be found in the control regions of these two mitogenomes. The phylogenetic analyses were performed on mitogenomes of Nymphalidae. The phylogenetic trees show that the relationships among Nymphalidae are generally identical to previous studies, as follows: Libytheinae\Danainae + ((Calinaginae + Satyrinae) + Danainae\Libytheinae + ((Heliconiinae + Limenitidinae) + (Nymphalinae + (Apaturinae + Biblidinae)))). Hestinalisnama is apart from Hestina, and closely related to Apatura, forming monophyly.

The complete mitochondrial genome of Monopis longella Walker, 1863 (Lepidoptera: Tineidae)

The complete mitochondrial genome (mitogenome) of Monopis longella Walker, 1863 (Lepidoptera: Tineidae) comprises 15,541 bp and contains a typical set of genes and one non-coding region. The gene arrangement of M. longella is unique for Lepidoptera in that it has a trnI-trnM-trnQ sequence in the A + T-rich region and ND2 junction. Unlike most other lepidopteran insects, in which the COI gene has CGA as the start codon, M. longella COI has an ATT codon. Phylogenetic analyses based on the concatenated sequences of 13 protein-coding genes and two rRNA genes, using the Bayesian inference (BI) method, placed M. longella in the Tineidae, sister in position to the cofamilial species, Tineola bisselliella, with the highest nodal support. Tineidae, represented by three species including M. longella, formed a monophyletic group with high support (Bayesian posterior probability = 0.99). Within Tineoidea the sister relationship between Tineidae and Meessiidae was obtained with the highest support, leaving Psychidae occupying the basal lineage of the two families.

The complete mitochondrial genome of Ripeacma umbellata Wang, 2009 (Lepidoptera: Autostichidae)

The mitogenome of Ripeacma umbellata Wang, 2009 was reported in this study. It was 15,486 bps long and strongly AT biased, consisting of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and 1 non-coding control region (351 bp). Most PCGs used the conventional ATN start codon, except for cox1 and cox2. Four genes used single T residue as stop codon rather than the routinely used TAA or TAG. All tRNAs, except for TrnS1, could fold into the cloverleaf secondary structure. Bayesian inference phylogenetic tree built on 13 PCGs from R. umbellata and another 21 species in Gelechioidea demonstrated that genus Ripeacma was a member in Autostichidae, which was consistent with the latest phylogenetic study.

Comparative Analysis of Mitochondrial Genome Features among Four Clonostachys Species and Insight into Their Systematic Positions in the Order Hypocreales

The mycoparasite fungi of Clonostachys have contributed to the biological control of plant fungal disease and nematodes. The Clonostachys fungi strains were isolated from Ophiocordyceps highlandensis, Ophiocordycepsnigrolla and soil, which identified as Clonostachyscompactiuscula, Clonostachysrogersoniana, Clonostachyssolani and Clonostachys sp. To explore the evolutionary relationship between the mentioned species, the mitochondrial genomes of four Clonostachys species were sequenced and assembled. The four mitogenomes consisted of complete circular DNA molecules, with the total sizes ranging from 27,410 bp to 42,075 bp. The GC contents, GC skews and AT skews of the mitogenomes varied considerably. Mitogenomic synteny analysis indicated that these mitogenomes underwent gene rearrangements. Among the 15 protein-coding genes within the mitogenomes, the nad4L gene exhibited the least genetic distance, demonstrating a high degree of conservation. The selection pressure analysis of these 15 PCGs were all below 1, indicating that PCGs were subject to purifying selection. Based on protein-coding gene calculation of the significantly supported topologies, the four Clonostachys species were divided into a group in the phylogenetic tree. The results supplemented the database of mitogenomes in Hypocreales order, which might be a useful research tool to conduct a phylogenetic analysis of Clonostachys. Additionally, the suitable molecular marker was significant to study phylogenetic relationships in the Bionectriaceae family.

Unique Duplication of trnN in Odontoptilum angulatum (Lepidoptera: Pyrginae) and Phylogeny within Hesperiidae

To explore the variation and relationship between gene rearrangement and phylogenetic effectiveness of mitogenomes among lineages of the diversification of the tribe Tagiadini in the subfamily Pyrginae, we sequenced the complete mitogenome of Odontoptilum angulatum. The genome is 15,361 bp with the typical 37 genes, a large AT-rich region and an additional trnN (trnN2), which is completely identical to trnN (sequence similarity: 100%). The gene order differs from the typical Lepidoptera-specific arrangement and is unique to Hesperiidae. The presence of a "pseudo-trnS1" in the non-coding region between trnN1 and trnN2 supports the hypothesis that the presence of an extra trnN can be explained by the tandem duplication-random loss (TDRL) model. Regarding the phylogenetic analyses, we found that the dataset comprising all 37 genes produced the highest node support, as well as a monophyly of Pyrginae, indicating that the inclusion of RNAs improves the phylogenetic signal. Relationships among the subfamilies in Hesperiidae were also in general agreement with the results of previous studies. The monophyly of Tagiadini is strongly supported. Our study provides a new orientation for application of compositional and mutational biases of mitogenomes in phylogenetic analysis of Tagiadini and even all Hesperiidae based on larger taxon sampling in the future.

Complete mitochondrial genome of Acanthopsyche nigraplaga (Lepidoptera: Psychidae)

We report the mitochondrial genome (mitogenome) of a case-making moth Acanthopsyche nigraplaga Wileman, 1911 (Lepidoptera: Psychidae). The 15,704 bp long complete mitogenome comprises a typical set of genes [13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes] and one major non-coding, A + T-rich region, with an arrangement identical to that observed in most lepidopteran mitogenomes. Twelve of the 13 PCGs of the A. nigraplaga mitogenome initiate with a typical ATN start codon, however COI contains the atypical CGA start codon that is common for lepidopteran COI genes. A phylogenetic analysis using concatenated nucleotide sequences of the 13 PCGs and 2 rRNA genes using the Bayesian inference method fully resolved A. nigraplaga in a monophyletic clade within the Psychidae. Acanthopsyche nigraplaga was situated in a sister position to Eumeta variegata and Mahasena oolona with high nodal support. As more mitogenome sequences are available further scrutinized analysis for the superfamily Tineoidea including Psychidae will be possible.

Characterization of the mitochondrial genome of Spodoptera exempta (Lepidoptera: Noctuidae) from South Africa

The African armyworm, Spodoptera exempta, is an episodic migratory crop pest with an expanding distribution worldwide. This is the first report of the circular mitochondrial genome of S. exempta, with a length of 15,457 bp and an A + T content of 81.7%. It encoded a common set of 37 genes, including 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes, and contained a putative control region of 379 bp (94.7% in A + T proportion). The maximum-likelihood phylogenetic tree based on the complete mitogenome demonstrated that five species belonging to the Spodoptera genus formed one clade, in which S. exempta was the most isolated branch, followed by Spodoptera exigua. This data will contribute for the identification and phylogenetic analyses of S. exempta, providing useful information for its comprehensive control.

The complete mitochondrial genome of the Indian leafwing butterfly Kallima paralekta (insecta: Lepidoptera: Nymphalidae)

The Indian leafwing butterfly Kallima paralekta (Horsfield, 1829) (Nymphalidae) is an Asian forest-dwelling, leaf-mimic. Genome skimming by Illumina sequencing permitted assembly of a complete circular mitogenome of 15,200 bp from K. paralekta consisting of 79.5% AT nucleotides, 22 tRNAs, 13 protein-coding genes, two rRNAs and a control region in the typical butterfly gene order. Kallima paralekta COX1 features an atypical CGA start codon, while ATP6, COX1, COX2, ND4, ND4L, and ND5 exhibit incomplete stop codons completed by 3' A residues added to the mRNA. Phylogenetic reconstruction places K. paraleckta within the monophyletic genus Kallima, sister to Mallika in the subfamily Nymphalinae. These data support the monophyly of tribe Kallimini and contribute to the evolutionary systematics of the Nymphalidae.

The complete mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) based on long-read nanopore sequencing

BACKGROUND

The coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is a severe and invasive pest of coconut and other palms throughout Asia and the Pacific. The biocontrol agent, Oryctes rhinoceros nudivirus (OrNV), has successfully suppressed O. rhinoceros populations for decades but new CRB invasions started appearing after 2007. A single-SNP variant within the mitochondrial cox1 gene is used to distinguish the recently-invading CRB-G lineage from other haplotypes, but the lack of mitogenome sequence for this species hinders further development of a molecular toolset for biosecurity and management programmes against CRB. Here we report the complete circular sequence and annotation for CRB mitogenome, generated to support such efforts.

METHODS

Sequencing data were generated using long-read Nanopore technology from genomic DNA isolated from a CRB-G female. The mitogenome was assembled with Flye v.2.5, using the short-read Illumina sequences to remove homopolymers with Pilon, and annotated with MITOS. Independently-generated transcriptome data were used to assess the O. rhinoceros mitogenome annotation and transcription. The aligned sequences of 13 protein-coding genes (PCGs) (with degenerate third codon position) from O. rhinoceros, 13 other Scarabaeidae taxa and two outgroup taxa were used for the phylogenetic reconstruction with the Maximum likelihood (ML) approach in IQ-TREE and Bayesian (BI) approach in MrBayes.

RESULTS

The complete circular mitogenome of O. rhinoceros is 20,898 bp in length, with a gene content canonical for insects (13 PCGs, two rRNA genes, and 22 tRNA genes), as well as one structural variation (rearrangement of trnQ and trnI) and a long control region (6,204 bp). Transcription was detected across all 37 genes, and interestingly, within three domains in the control region. ML and BI phylogenies had the same topology, correctly grouping O. rhinoceros with one other Dynastinae taxon, and recovering the previously reported relationship among lineages in the Scarabaeidae. In silico PCR-RFLP analysis recovered the correct fragment set that is diagnostic for the CRB-G haplogroup. These results validate the high-quality of the O. rhinoceros mitogenome sequence and annotation.

Complete mitochondrial genome of Pterodecta felderi (Lepidoptera: Callidulidae)

We report the mitochondrial genome (mitogenome) of Pterodecta felderi (Callidulidae: Lepidoptera), which is the first mitogenome sequences in the family Callidulidae, a monotypic family in the superfamily Calliduloidea. The 15,340-bp long complete mitogenome consists of a typical set of genes (13 protein-coding genes [PCGs], 2 rRNA genes, and 22 tRNA genes) and 1 major non-coding A + T-rich region, which are arranged in a way that is frequently observed in Lepidoptera. Of the 13 PCGs, 12 P. felderi start with ATN, except for COI, which starts with CGA. The P. felderi mitogenome consists of 210-bp long intergenic-spacer sequences and 27-bp long overlaps. Phylogenetic analysis of superfamilial relationships in the lepidopteran clade Obtectomera with concatenated sequences of the 13 PCGs and 2 rRNA genes using the Bayesian inference method showed that Calliduloidea, which is only represented by P. felderi, was placed as the most basal lineage about Macroheterocera (Lasiocampoidea, Bombycoidea, Mimallonoidea, Noctuoidea, and Drepanoidea), Papilionoidea, and Pyraloidea.

Phylogenetic analysis of the complete mitochondrial genome of the white peacock butterfly Anartia jatrophae saturata (Insecta: Lepidoptera: Nymphalidae)

The white peacock butterfly Anartia jatrophae saturata Staudinger, 1884 (Nymphalidae: Nymphalinae: Victorini), lives in the neotropics. Genome skimming with Illumina sequencing of A. jatrophae saturata allowed the assembly of a complete circular mitogenome of 15,297 bp, consisting of 81.4% AT nucleotides, 22 tRNAs, 13 protein-coding genes, two rRNAs, and a control region. Anartia jatrophae COX1 features an atypical start codon (CGA); ATP6, COX1, ND1, ND4, ND4L, ND5, and ND6 exhibit incomplete stop codons completed in the mRNA by the addition of 3' A residues. Contrary to previous phylogenetic hypotheses, phylogenetic reconstruction places A. jatrophae as sister to nymphalid tribe Nymphalini.

The complete mitochondrial genome of the Jackson's leaf butterfly Mallika jacksoni (Insecta: Lepidoptera: Nymphalidae)

The Jackson's leaf butterfly Mallika jacksoni (Sharpe 1896), is a leaf-mimicking species from tropical East Africa. Genome skimming by Illumina sequencing permitted the assembly of the complete circular M. jacksoni 15,183 bp mitogenome. It consists of 79.4% AT nucleotides, 22 tRNAs, 13 protein-coding genes, 2 rRNAs, and a control region in the typical butterfly gene order. Mallika jacksoni COX1 has a CGA start codon while ATP6, COX1, COX2, ND3, ND4, and ND5 exhibit partial stop codons completed by 3'-A residues added to the mRNA. Phylogenetic reconstruction places M. jacksoni as sister to Kallima within nymphalid tribe Kallimini.

Complete mitochondrial genome of Spilosoma lubricipedum (Noctuoidea: Erebidae) and implications for phylogeny of noctuid insects

Mitochondrial genomes (mitogenomes) have been widely used for studies on phylogenetic relationships and molecular evolutionary biology. Here, the complete mitogenome sequence of Spilosoma lubricipedum (Noctuoidea: Erebidae: Arctiinae) was determined (total length 15,375 bp) and phylogenetic analyses S. lubricipedum were inferred from available noctuid sequence data. The mitogenome of S. lubricipedum was found to be highly A + T-biased (81.39%) and exhibited negative AT- and GC-skews. All 13 protein-coding genes (PCGs) were initiated by ATN codons, except for cox1 with CGA. All tRNAs exhibited typical clover-leaf secondary structures, except for trnS1. The gene order of the S. lubricipedum mitogenome was trnM-trnI-trnQ-nad2. The A + T-rich region of S. lubricipedum contained several conservative features common to noctuid insects. Phylogenetic analysis within Noctuoidea was carried out based on mitochondrial data. Results showed that S. lubricipedum belonged to Erebidae and the Noctuoidea insects could be divided into five well-supported families (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))).

Evolution of tRNA gene rearrangement in the mitochondrial genome of ichneumonoid wasps (Hymenoptera: Ichneumonoidea)

Gene rearrangements in the mitochondrial genome (mt genome) are common in certain insect groups and can be an informative character for phylogenetic reconstruction. However, knowledge of the mechanism and biases of gene rearrangement in insect mt genomes is still limited. With an accelerated rate of gene rearrangements, Hymenoptera is an important group for mt genome rearrangements diversity and for understanding the gene rearrangement evolution in mt genomes. Here, we sequenced the complete mt genome of Aphidius gifuensis and analyzed the evolution of tRNA gene rearrangements in the mt genomes of ichneumonoid wasps. Two control regions were detected in A. gifuensis and most of the tRNA rearrangement events occurred around these control regions. tRNA gene rearrangements occurred in almost all of the sequenced mt genomes of Ichneumonoidea and the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was the main hot spot of gene rearrangement. Mapped over the backbone phylogeny of Ichneumonoidea, we found that the inversion and translocation of both trnI and trnM is likely a synapomorphic rearrangement in Braconidae. Our study also demonstrated that the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was important for inferring the gene rearrangement dynamics in Ichneumonoidea.

Comparative mitochondrial genome analysis of Dendrolimus houi (Lepidoptera: Lasiocampidae) and phylogenetic relationship among Lasiocampidae species

Dendrolimus houi is one of the most common caterpillars infesting Gymnosperm trees, and widely distributed in several countries in Southeast Asia, and exists soley or coexists with several congeners and some Lasiocampidae species in various forest habitats. However, natural hybrids occasionally occur among some closely related species in the same habitat, and host preference, extreme climate stress, and geographic isolation probably lead to their uncertain taxonomic consensus. The mitochondrial DNA (mtDNA) of D. houi was extracted and sequenced by using high-throughput technology, and the mitogenome composition and characteristics were compared and analyzed of these species, then the phylogenetic relationship was constructed using the maximum likelihood method (ML) and the Bayesian method (BI) based on their 13 protein-coding genes (PCGs) dataset, which were combined and made available to download which were combined and made available to download among global Lasiocampidae species data. Mitogenome of D. houi was 15,373 bp in length, with 37 genes, including 13 PCGs, 22 tRNA genes (tRNAs) and 2 rRNA genes (rRNAs). The positions and sequences of genes were consistent with those of most known Lasiocampidae species. The nucleotide composition was highly A+T biased, accounting for ~80% of the whole mitogenome. All start codons of PCGs belonged to typical start codons ATN except for COI which used CGA, and most stop codons ended with standard TAA or TAG, while COI, COII, ND4 ended with incomplete T. Only tRNASer (AGN) lacked DHU arm, while the remainder formed a typical "clover-shaped" secondary structure. For Lasiocampidae species, their complete mitochondrial genomes ranged from 15,281 to 15,570 bp in length, and all first genes started from trnM in the same direction. And base composition was biased toward A and T. Finally, both two methods (ML and BI) separately revealed that the same phylogenetic relationship of D. spp. as ((((D. punctatus + D. tabulaeformis) + D. spectabilis) + D. superans) + (D. kikuchii of Hunan population + D. houi) as in previous research, but results were different in that D. kikuchii from a Yunnan population was included, indicating that different geographical populations of insects have differentiated. And the phylogenetic relationship among Lasiocampidae species was ((((Dendrolimus) + Kunugia) + Euthrix) + Trabala). This provides a better theoretical basis for Lasiocampidae evolution and classification for future research directions.

Phylogenetic relationships of Limacodidae and insights into the higher phylogeny of Lepidoptera

To determine the systematic status of family Limacodidae within Lepidoptera, the complete mitochondrial genome (mitogenome) of Thosea sinensis (Lepidoptera: Zygaenoidea: Limacodidae) was sequenced. The genome is 15,544 base pairs (bp), including 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and an AT-rich region. These characteristics are similar to of other lepidopterans. The gene order of T. sinensis is identical to that of Ditrysia lepidopterans. The nucleotide composition of the T. sinensis mitochondrial genome is highly biased toward A + T nucleotides (81.1%) and exhibits negative AT and GC skew. All the other 13 PCGs except cox1 are initiated by ATN codons. All tRNA genes are folded into the typical cloverleaf secondary structure, except for trnS1, which lacked the dihydrouridine (DHU) stem. There are 20 intergenic spacer regions ranging from 1 to 56 bp in length, and two gene overlap regions throughout the entire genome. The AT-rich region includes the ATAGA motif, followed by a 19-bp poly T stretch, a microsatellite-like (AT)₁₀, and a poly-A element. Analysis of phylogenetic relationships indicated that T. sinensis belongs to the Limacodidae, and the monophyly of each lepidopteran family was well supported.

First mitochondrial genome from Yponomeutidae (Lepidoptera, Yponomeutoidea) and the phylogenetic analysis for Lepidoptera

The complete mitochondrial genome (mitogenome) of Yponomeuta montanatus is sequenced and compared with other published yponomeutoid mitogenomes. The mitogenome is circular, 15,349 bp long, and includes the typical metazoan mitochondrial genes (13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes) and an A + T-rich region. All 13 protein-coding genes use a typical start codon ATN, the one exception being cox1, which uses CGA across yponomeutoid mitogenomes. Comparative analyses further show that the secondary structures of tRNAs are conserved, including loss of the Dihydorouidine (DHU) arm in trnS1 (AGN), but remarkable nucleotide variation has occurred mainly in the DHU arms and pseudouridine (TψC) loops. A + T-rich regions exhibit substantial length variation among yponomeutoid mitogenomes, and conserved sequence blocks are recognized but some of them are not present in all species. Multiple phylogenetic analyses confirm the position of Y. montanatus in Yponomeutoidea. However, the superfamily-level relationships in the Macroheterocera clade in Lepidoptera recovered herein show considerable difference with that recovered in previous mitogenomic studies, raising the necessity of extensive phylogenetic investigation when more mitogenomes become available for this clade.

Complete mitochondrial genome of two Thitarodes species (Lepidoptera, Hepialidae), the host moths of Ophiocordyceps sinensis and phylogenetic implications

Thitarodes (Lepidoptera, Hepialidae) is the only genus that hosts to the Ophiocordyceps sinensis, a traditional Chinese medicine considered as a powerful medicinal supplement. In this study, the complete mitochondrial genomes (mitogenomes) of two species, T. damxungensis and T. pui, have been sequenced, which are 15,928 bp and 15,362 bp in size respectively, and both contain 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs and an AT-rich region. Like other hepialoids, the gene arrangement of the mitogenomes of T. damxungensis and T. pui is identical to the ancestral arrangement but differs from those of other lepidopteran species on account of the different arrangements of trnM, trnI, and trnQ. The size of AT-rich region is 545 bp in T. damxungensis and 1030 bp in T. pui. Tandem repetition in the AT-rich region is responsible for the length difference of the A + T-rich region in both species. In Hepialidae, the phylogenetic study based on the dataset of the sequences that combined the protein-coding genes and RNA genes suggested that the species T. yunnanensis should still belong to the genus Thitarodes rather than Ahamns, which is different from the results based on the traditional phylogeny.

The first mitochondrial genome of the family Epicopeiidae and higher-level phylogeny of Macroheterocera (Lepidoptera: Ditrysia)

The Macroheterocera clade contains most of the Lepidoptera species. However, extensive comparative and phylogenetic analyses of this group using complete mitochondrial genomes (mitogenome) are limited particularly in the context of increasing macroheteroceran mitogenomes reported to date. In this study, complete mitogenome of the Epicopeia hainesii is determined as the first Epicopeiidae species with mitogenome available. The whole mitogenome is circular with 15,395 bp long, and is highly biased toward A + T nucleotides (80.6%) in nucleotide composition. Comparative analyses show that gene content and arrangement of macroheteroceran mitogenomes are generally conservative and are typical of Lepidoptera but exceptions exist. In newly sequenced mitogenome, the motif "ATACTAA" is putatively located at the end of gene nad1, rather than in intergenic sequences between trnS2 and nad1 genes routinely observed in Lepidoptera. Interestingly, multiple phylogenetic analyses recover the six macroheteroceran superfamilies as Mimallonoidea + (Drepanoidea + ((Bombycoidea + Lasiocampoidea) + (Noctuoidea + Geometroidea))), providing supports for a large-scale transcriptomic study rather than various mitogenome- and multiple-gene-based investigations. In addition, our analyses consistently place the Epicopeiidae as sister group with Geometroidea, firstly demonstrating that this family is closer with Geometroidea than Drepanoidea based on mitogenome data.

Phylogeographic structures of the host insects of Ophiocordyceps sinensis

A fungus-insect complex, known as DongChong XiaCao, is formed from the infection of the hepialid larvae by the fungus Ophiocordyceps sinensis, which is endemic to the Qinghai-Tibetan Plateau (QTP). Due to previously limited sample collection size, the data about the diversity and structure of the host insect was insufficient and lacked details. The purpose of this study was aimed to discuss the diversity and phylogeography of the host insects of O. sinensis with a large-scale sampling. The mitochondrial cytochrome oxidase I gene (cox1) was sequenced and analyzed among 710 samples representing 88 geographic locations. 205 haplotypes of cox1 were identified from all the 710 samples and 4 phylogenetic clades with 12 subclades were identified. Instead of the single latitude-based divergence suggested previously, three distribution patterns were deduced to correspond to the phylogeographic structures, including but not limited to the co-existence of a wide and specific local phylogeographic distribution structures. Two separate genetic diversity and differentiation centers, namely the northwestern Yunnan and the southeastern Tibet were identified. Dating analyses from three calibrations supported that the divergence of the 4 clades occurred in the Oligocene-Miocene period (30.54-13.66 million years ago) (Ma), which were connected with the second and third geological movements of the QTP (17-25, 8-13 Ma). Our results provide a more detailed understanding of the divergence and distribution patterns of the host insects of O. sinensis.