Description of new mitochondrial genomes (Spodoptera litura, Noctuoidea and Cnaphalocrocis medinalis, Pyraloidea) and phylogenetic reconstruction of Lepidoptera with the comment on optimization schemes

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.

Characterization of the mitochondrial genome of the lawn cutworm Spodoptera depravata (Lepidoptera: Noctuidae)

The lawn cutworm, Spodoptera depravata, is one of the most important pests that causes economic damage to grass crops. This study reports the complete mitochondrial genome of an S. depravata sample collected in China. The genome is a circular molecule 15,460 bp in length with an overall A + T content of 81.6%. It contains 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. The gene content and organization of the mitogenome of S. depravata are identical to those of other Spodoptera species. Maximum-likelihood phylogenetic analysis based on mitogenomes showed a close evolutionary relationship between S. depravata and S. exempta. This study provides new molecular data for the identification and further phylogenetic analyses of Spodoptera species.

Complete mitochondrial genome of Idea leuconoe (Lepidoptera: Danaidae) and related phylogenetic analyses

In the present study, we first sequenced and determined the complete mitochondrial genome (mitogenome) of the giant Danaidae butterfly, Idea leuconoe (Lepidoptera: Danaidae). The mitogenome was a typical closed, circular, double-stranded DNA molecule of 15,278 bp length (GenBank accession number: KR815449), similar to the metazoan mitogenomes containing 37 genes and one A + T-rich region. All the protein-coding genes (PCGs) were initiated with a typical ATN codon. Seven genes (COII, ATP6, COIII, nad4, nad4L, cytb, and nad1) adopted the standard ATG start codon, but the remaining six genes were initiated with ATA. All the 13 PCGs harbored complete termination codons (TAA). The overlap nucleotides ATGATAA were conserved for the ATP8/ATP6 gene. The largest intergenic spacer was located between trnGln and nad2, a common finding in Lepidoptera butterflies. All the transfer RNA genes in the I. leuconoe mitogenome could be folded into typical clover-leaf secondary structures, except for trnSer (AGN) that lacked a dihydrouridine arm. The control region with 94.8% A + T content was 444 bp in length and located between rrnS and trnMet. Finally, the phylogenetic relationships obtained using the maximum likelihood and Bayesian methods confirmed two well-supported phylogenetic trees of Danaidae, Papilionidae, and Nymphalidae from the order Lepidoptera, which were consistent with the traditional morphological classification. Results provided additional information for butterfly phylogenetic analysis and insights into the evolution of genomes.

Nine Mitochondrial Genomes of the Pyraloidea and Their Phylogenetic Implications (Lepidoptera)

Simple Summary

The Pyraloidea is a large superfamily of Lepidoptera in species composition. To date, the higher-level phylogenetic relationships in this group remain unresolved, and many taxa, with taxonomic positions historically established by morphological characters, need to be confirmed through sequencing of DNA, including mitochondrial genome sequences (mitogenomes). Here, we newly generated nine complete mitogenomes for Pyraloidea that shared identical gene content, and arrangements that are typical of Lepidoptera. The current phylogenetic results confirmed previous multilocus studies, indicating the effectiveness of mitogenomes for inference of Pyraloidea higher-level relationships. Unexpectedly, Orybina Snellen was robustly placed as basal to the remaining Pyralidae taxa, rather than nested in the Pyralinae of Pyralidae as morphologically defined and placed. Our results bring a greater understanding to Pyraloidea phylogeny, and highlight the necessity of sequencing more pyraloid taxa to reevaluate their phylogenetic positions.

Abstract

The Pyraloidea is one of the species-rich superfamilies of Lepidoptera and contains numerous economically important pest species that cause great loss in crop production. Here, we sequenced and annotated nine complete mitogenomes for Pyraloidea, and further performed various phylogenetic analyses, to improve our understanding of mitogenomic evolution and phylogeny of this superfamily. The nine mitogenomes were circular, double-stranded molecules, with the lengths ranging from 15,214 bp to 15,422 bp, which are comparable to other reported pyraloid mitogenomes in size. Gene content and arrangement were highly conserved and are typical of Lepidoptera. Based on the hitherto most extensive mitogenomic sampling, our various resulting trees showed generally congruent topologies among pyraloid subfamilies, which are almost in accordance with previous multilocus studies, indicating the suitability of mitogenomes in inferring high-level relationships of Pyraloidea. However, nodes linking subfamilies in the “non-PS clade” were not completely resolved in terms of unstable topologies or low supports, and future investigations are needed with increased taxon sampling and molecular data. Unexpectedly, Orybina Snellen, represented in a molecular phylogenetic investigation for the first time, was robustly placed as basal to the remaining Pyralidae taxa across our analyses, rather than nested in Pyralinae of Pyralidae as morphologically defined. This novel finding highlights the need to reevaluate Orybina monophyly and its phylogenetic position by incorporating additional molecular and morphological evidence.

A Likelihood-Ratio Test for Lumpability of Phylogenetic Data: Is the Markovian Property of an Evolutionary Process retained in Recoded DNA?

In molecular phylogenetics, it is typically assumed that the evolutionary process for DNA can be approximated by independent and identically distributed Markovian processes at the variable sites and that these processes diverge over the edges of a rooted bifurcating tree. Sometimes the nucleotides are transformed from a 4-state alphabet to a 3- or 2-state alphabet by a procedure that is called recoding, lumping, or grouping of states. Here, we introduce a likelihood-ratio test for lumpability for DNA that has diverged under different Markovian conditions, which assesses the assumption that the Markovian property of the evolutionary process over each edge is retained after recoding of the nucleotides. The test is derived and validated numerically on simulated data. To demonstrate the insights that can be gained by using the test, we assessed two published data sets, one of mitochondrial DNA from a phylogenetic study of the ratites (Syst. Biol. 59:90-107 [2010]) and the other of nuclear DNA from a phylogenetic study of yeast (Mol. Biol. Evol. 21:1455-1458 [2004]). Our analysis of these data sets revealed that recoding of the DNA eliminated some of the compositional heterogeneity detected over the sequences. However, the Markovian property of the original evolutionary process was not retained by the recoding, leading to some significant distortions of edge lengths in reconstructed trees.

Complete mitochondrial genome of Ostrinia kasmirica (Lepidoptera: Crambidae)

The complete mitochondrial genome of Ostrinia kasmirica (Moore, 1888) was sequenced in this study. The circular mitogenome is 15,214 bp in length, containing 37 typical encoded genes and a non-coding control region. The gene organization and nucleotide composition are similar to those of most other sequenced Ostrinia species. All protein-coding genes (PCGs) initiate with ATN and terminate with TAN, except cox1 starts with CGA and cox1, cox2, nad5 terminate with an incomplete codon T. The control region of 308 bp contains three conserved features including the motif ‘TTAGA’ preceded a poly-T stretch, a microsatellite-like (TA)_n element, and a poly-A stretch upstream of trnM. Phylogenetic analysis based on mitogenome sequences revealed that the O. kasmirica (the second species group) was more closely related to the third species group of the genus and the first species group was not at the basal position of this genus as that Mutuura and Munroe indicated.

The First Mitogenomes of the Subfamily Odontiinae (Lepidoptera, Crambidae) and Phylogenetic Analysis of Pyraloidea

Simple Summary

The Odontiinae is a small group in the Pyraloidea comprised of 388 species in 88 genera, but externally, these moths are diverse, including heterogeneous maculation and a size range from 9 to 50 mm in total wingspan. The monophyly of Pyraloidea and the two families (Pyralidae and Crambidae) is well supported by phylogenetic analyses based on morphology and molecular data of multiple nuclear genes. However, only a few mito-phylogenetic analyses have been conducted and no mitogenome of Odontiinae species has been reported. Three complete mitogenomes of odontiine species were sequenced and analyzed for the first time herein. The results showed that Odontiinae mitogenomes shared similar genomic characters with other Pyraloidea. The phylogenetic analyses based on 13 PCGs of mitogenomes confirmed the monophyly of Odontiinae and its position within Crambidae.

Abstract

The complete mitochondrial genomes of three species of Odontiinae were newly sequenced: Dausara latiterminalis Yoshiyasu, Heortia vitessoides (Moore), and Pseudonoorda nigropunctalis (Hampson). These circular and double-stranded mitogenomes vary from 15,084 bp to 15,237 bp in size, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs) and an A + T-rich region. The nucleotide composition indicated a strong A/T bias. Most PCGs are initiated with an ATN codon and terminated by a codon of TAR. All tRNAs could be folded into the clover-leaf structure with the exception of trnS1 (AGN), in which the dihydrouridine (DHU) arm formed a simple loop, and the motif ‘ATAG’ and ‘ATTTA’ in the A + T-rich region was also founded. The phylogenomic analyses covering Odontiinae + 11 subfamilies of Pyraloidea were conducted. Similar topologies were generated from both Bayesian inference (BI) and maximum likelihood (ML) analyses based on the nucleotide and amino acid sequence data. There was some discrepancy in the sister-group relationship of Odontiinae and Glaphyriinae, and the relationships among the subfamilies in the ‘CAMMSS clade’ of the Crambidae. The results of this study suggest that mitogenomic data are useful for resolving the deep-level relationships of Pyraloidea and the topologies generated from amino acid data might be more realistic and reliable. Moreover, more mitogenomic taxon sampling and larger scale analyses with more genes or a combination of mitogenomic and nuclear genes are needed to reconstruct a comprehensive framework of the pyraloid phylogeny.

Complete mitochondrial genome of Actias dubernardi (Lepidoptera: Saturniidae)

The complete mitochondrial genome of Actias dubernardi (Lepidoptera: Saturniidae) is 15,270 bp in length, containing 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a putative control region. All of the protein-coding genes (PCGs) use the standard start codon ATN, except for cox1 which starts with CGA. The Bayesian phylogenetic analysis was performed using a dataset matrix containing 13 PCGs concatenated from the mitogenomes of 14 Saturniidae species. The monophyly of the five Actias species was highly supported and Antheraea was inferred as the sister group of Actias.

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.

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.

Comparative Mitochondrial Analysis of Cnaphalocrocis exigua (Lepidoptera: Crambidae) and Its Close Relative C. medinalis

Rice leaffolders are important pests on rice in Asia, Oceania, and Africa, causing serious loss to rice production. There are two main rice leaffolders in China, namely Cnaphalocrocis medinalis (Guenée) and C. exigua (Butler) with the former having the ability of long-distance migration. To reveal the differences in the mitochondrial genomes (mitogenome) between them, we compared the completed mitogenome of C. exigua with three C. medinalis individuals. Although phylogenetic analysis based on the mitogenomic data strongly supported the close relationship between these two species, many differences were still being revealed. The results showed that the mitogenome of C. exigua was shorter in length (15,262 bp) and slight lower in AT content than that of C. medinalis. Except for the different start codons of nad3 and nad6 gene, we also found the cox1 gene had a typical start codon 'ATG' which suggested that the starting position of this gene must be reconsidered in the entire superfamily Pyraloidea. All tRNAs have a typical clover-leaf structure, except for the dihydrouridine (DHU) stem losing of trnS1, which has the atypical anticondon 'TCT' instead of 'GCT' in C. medinalis and most Pyraloidea species. Two intergenic regions (between trnY and cox1, nad3 and trnA) featured by AT repeats were only found in C. medinalis and even rarely appeared in reported Pyraloidea species. Furthermore, regardless of interspecific comparison or intraspecific comparison of these two species, protein coding genes, especially the atp8 genes, had quite different evolutionary rates.

Complete Mitogenomic Structure and Phylogenetic Implications of the Genus Ostrinia (Lepidoptera: Crambidae)

To understand mitogenome characteristics and reveal phylogenetic relationships of the genus Ostrinia, including several notorious pests of great importance for crops, we sequenced the complete mitogenomes of four species: Ostrinia furnacalis (Guenée, 1854), Ostrinia nubilalis (Hübner, 1796), Ostrinia scapulalis (Walker, 1859) and Ostrinia zealis (Guenée, 1854). Results indicate that the four mitogenomes-O. furnacalis, O. nubilalis, O. scapulalis, and O. zealis-are 15,245, 15,248, 15,311, and 15,208 bp in size, respectively. All four mitogenomes are comprised of 37 encoded genes and a control region. All 13 protein-coding genes (PCGs) initiate with ATN and terminate with TAN, with the exception of cox1 that starts with CGA, and cox1, cox2, and nad5 that terminate with an incomplete codon T. All transfer RNA genes (tRNAs) present the typical clover-leaf secondary structure except for the trnS1 (AGN) gene. There are some conserved structural elements in the control region. Our analyses indicate that nad6 and atp6 exhibit higher evolution rates compared to other PCGs. Phylogenetic analyses based on mitogenomes using both maximum likelihood (ML) and Bayesian inference (BI) methods revealed the relationship (O. palustralis + (O. penitalis + (O. zealis + (O. furnacalis + (O. nubilalis + O. scapulalis))))) within Ostrinia.

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.

Mitochondrial DNA genomes of five major Helicoverpa pest species from the Old and New Worlds (Lepidoptera: Noctuidae)

Five species of noctuid moths, Helicoverpa armigera, H. punctigera, H. assulta, H. zea, and H. gelotopoeon, are major agricultural pests inhabiting various and often overlapping global distributions. Visual identification of these species requires a great deal of expertise and misidentification can have repercussions for pest management and agricultural biosecurity. Here, we report on the complete mitochondrial genomes of H. assulta assulta and H. assulta afra, H. gelotopoeon, H. punctigera, H. zea, and H. armigera armigera and H. armigera conferta' assembled from high-throughput sequencing data. This study significantly increases the mitogenome resources for these five agricultural pests with sequences assembled from across different continents, including an H. armigera individual collected from an invasive population in Brazil. We infer the phylogenetic relationships of these five Helicoverpa species based on the 13 mitochondrial DNA protein-coding genes (PCG's) and show that two publicly available mitogenomes of H. assulta (KP015198 and KR149448) have been misidentified or incorrectly assembled. We further consolidate existing PCR-RFLP methods to cover all five Helicoverpa pest species, providing an updated method that will contribute to species differentiation and to future monitoring efforts of Helicoverpa pest species across different continents. We discuss the value of Helicoverpa mitogenomes to assist with species identification in view of the context of the rapid spread of H. armigera in the New World. With this work, we provide the molecular resources necessary for future studies of the evolutionary history and ecology of these species.

The complete mitochondrial genome of Parantica sita sita (Lepidoptera: Nymphalidae: Danainae) revealing substantial genetic divergence from its sibling subspecies P. s. niphonica

Currently, there are two subspecies of the chestnut tiger butterfly (Parantica sita) recognized in China. P. s. sita is widely distributed in southwest China and P. s. niphonica in Taiwan. Periodically, Taiwan Island and the Chinese mainland have been connected and separated because of sea level changes caused by Pleistocene glaciations, most likely influencing the genetic structure among P. sita populations on both sides of the Taiwan Strait. Also, P. s. niphonica's well-documented long-distance migration may have influenced genetic differentiation within this species as well. So, investigation of the genetic differentiation of these two subspecies is well warranted. In this study, we sequenced the complete mitogenome (15,156 bp in length) of P. s. sita and its general characteristics agreed with general butterfly mitogenomic characteristics. However, when compared genetically with P. s. niphonica, P. s. sita diverged substantially. First, there were 4.1% variable sites between these two subspecies, with 5.3% at COI and 3.8% at COII, differences much greater than those expected of general interspecific divergences in butterflies. Second, there was a 59 bp fragment deletion in the A + T rich region of P. s. sita and, third, the relationships of these two subspecies and P. luzonensis could not be distinguished using Bayesian inference and P. s. niphonica first clustered with P. luzonensis, rather than P. s. sita, using maximum likelihood. Based on these results, we propose that P. s. sita and P. s. niphonica are independent species instead of subspecies. This proposal should be clarified through further research.

Comparative mitochondrial genome analysis of Grammodes geometrica and other noctuid insects reveals conserved mitochondrial genome organization and phylogeny

The mitochondrial genome (mitogenome) plays an important role in revealing molecular evolution. In this study, the complete mitogenome of Grammodes geometrica (G. geometrica) (Lepidoptera: Erebidae) was sequenced and characterized. The nucleotide composition of the genome is highly A + T biased, accounting for 80.49%. Most protein-coding genes (PCGs) are initiated by ATN codons except for the cytochrome oxidase subunit 1 (cox1) gene, which was initiated by CGA. The order and orientation of genes with the order trnM-trnI-trnQ-nad2 is a typical rearrangement compared with those ancestral insects in which trnM is located between trnQ and nad2. Most tRNA genes were folded into the typical cloverleaf structure except for trnS1 (AGN). The A + T-rich region contains the conserved motif "ATAGA" followed by a 19 bp poly-T stretch, which was also observed in other Noctuoidea species. In addition, we reconstructed phylogenetic trees among the nucleotide alignments of five families of Noctuoidea species except the Oenosandridae. Finally, we achieved a well-supported tree, which showed that G. geometrica belongs to the Erebidae family. Moreover, the relationships at the family-level can be displayed as follows: (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))).

Detection of sister-species in invasive populations of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) from Uganda

The fall armyworm (FAW) Spodoptera frugiperda (J. E. Smith) is a species native to the Americas. This polyphagous lepidopteran pest was first reported in Nigeria and the Democratic Republic of São Tomé and Principe in 2016, but its presence in eastern Africa has not been confirmed via molecular characterisation. In this study, FAW specimens from western and central Uganda were identified based on the partial mtDNA COI gene sequences, with mtDNA COI haplotypes matching those identified in Nigeria and São Tomé. In this study, we sequence an additional partial mtDNA Cyt b gene and also the partial mtDNA COIII gene in Ugandan FAW samples. We detected identical mitochondrial DNA haplotypes for both the mtDNA Cyt b and COI partial genes, while combining the mtDNA COI/Cyt b haplotypes and mtDNA COIII haplotypes enabled a new maternal lineage in the Ugandan corn-preferred FAW samples to be identified. Our results suggested that the African incursions of S. frugiperda involved at least three maternal lineages. Recent full genome, phylogenetic and microsatellite analyses provided evidence to support S. frugiperda as likely consisted of two sympatric sister species known as the corn-preferred and rice-preferred strains. In our Ugandan FAW populations, we identified the presence of mtDNA haplotypes representative of both sister species. It is not known if both FAW sister species were originally introduced together or separately, and whether they have since spread as a single population. Further analyses of additional specimens originally collected from São Tomé, Nigeria and throughout Africa would be required to clarify this issue. Importantly, our finding showed that the genetic diversity of the African corn-preferred FAW species is higher than previously reported. This potentially contributed to the success of FAW establishment in Africa. Furthermore, with the additional maternal lineages detected, there is likely an increase in paternal lineages, thereby increasing the diversity of the African FAW population. Knowledge of the FAW genetic diversity will be needed to assess the risks of introducing Bt-resistance traits and to understand the FAW incursion pathways into the Old World and its potential onward spread. The agricultural implications of the presence of two evolutionary divergent FAW lineages (the corn and the rice lineage) in the African continent are further considered and discussed.

Complete mitochondrial genome of Cnidocampa flavescens (Lepidoptera: Limacodidae)

Cnidocampa flavescens, lives in Nepal, Bhutan, China, Far East of Russia, Korea, and Japan, belongs to the Lepidoptera family Limacodidae. In this study, we describe the genomic features of the mitogenome sequences of the insects. The mitogenome of C. flavescens is 15,406 bp long consisting a typical set of genes (13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes) and one major 415 bp non-coding A + T-rich region. All PCGs of C. flavescens start with ATN codons and end with TAA codons. The gene arrangement of C. flavescens mitogenome is same to Monema flavescens while the intergenic spacers and overlaps are different. The 415 bp A + T-rich region contains a conserved ATAGA motif followed a poly-T stretch.

Complete mitochondrial genome of the lappet moth, Kunugia undans (Lepidoptera: Lasiocampidae): genomic comparisons among macroheteroceran superfamilies

The mitochondrial genome (mitogenome) characteristics of the monotypic Lasiocampoidea are largely unknown, because only limited number of mitogenomes is available from this superfamily. In this study, we sequenced the complete mitogenome of the lappet moth, Kunugia undans (Lepidoptera: Lasiocampidae) and compared it to those of Lasiocampoidea and macroheteroceran superfamilies (59 species in six superfamilies). The 15,570-bp K. undans genome had one additional trnR that was located between trnA and trnN loci and this feature was unique in Macroheterocera, including Lasiocampoidea. Considering that the two trnR copies are located in tandem with proper secondary structures and identical anticodons, a gene duplication event might be responsible for the presence of the two tRNAs. Nearly all macroheteroceran species, excluding Lasiocampoidea, have a spacer sequence (1-34 bp) at the trnS2 and ND1 junction, but most lasiocampid species, including K. undans, have an overlap at the trnS2 and ND1 junction, which represents a different genomic feature in Lasiocampoidea. Nevertheless, a TTAGTAT motif, which is typically detected in Macroheterocera at the trnS2 and ND1 junction, was also detected in all Lasiocampoidea. In summary, the general mitogenome characteristics of Lasiocampoidea did not differ greatly from the remaining macroheteroceran superfamilies, but it did exhibit some unique features.

Complete mitochondrial genomes of two gelechioids, Mesophleps albilinella and Dichomeris ustalella (Lepidoptera: Gelechiidae), with a description of gene rearrangement in Lepidoptera

We sequenced the entire mitochondrial genome (mitogenome) of two gelechioids, Mesophleps albilinella and Dichomeris ustalella, and compared their genome organization and sequence composition to those of available gelechioid mitogenomes for an enhanced understanding of Gelechioidea genomic characteristics. We compared all available lepidopteran mitogenome arrangements, including that of M. albilinella, which is unique in Gelechioidea, to comprehend the extensiveness and mechanisms of gene rearrangement in Lepidoptera. The genomes of M. albilinella and D. ustalella are 15,274 and 15,410 bp in size, respectively, with the typical sets of mitochondrial (mt) genes. The COI gene begins with CGA (arginine) in all sequenced gelechioids, including M. albilinella and D. ustalella, reinforcing the feature as a synapomorphic trait, at least in the Gelechioidea. Each 353- and 321-bp long A + T-rich region of M. albilinella and D. ustalella contains one (D. ustalella) or two (M. albilinella) tRNA-like structures. The M. albilinella mitogenome has a unique gene arrangement among the Gelechioidea: ARNESF (the underline signifies an inverted gene) at the ND3 and ND5 junction, as opposed to the ARNSEF that is found in ancestral insects. An extensive search of available lepidopteran mitogenomes, including that of M. albilinella, turned up six rearrangements that differ from those of ancestral insects. Most of the rearrangements can be explained by the tandem duplication-random loss model, but inversion, which requires recombination, is also found in two cases, including M. albilinella. Excluding the MIQ rearrangement at the A + T-rich region and ND2 junction, which is found in nearly all Ditrysia, most of the remaining rearrangements found in Lepidoptera appear to be independently derived in that they are automorphic at several taxonomic scales, although current mitogenomic data are limited, particularly for congeneric data.

Reassessment of Species Diversity of the Subfamily Denticollinae (Coleoptera: Elateridae) through DNA Barcoding

The subfamily Denticollinae is a taxonomically diverse group in the family Elateridae. Denticollinae includes many morphologically similar species and crop pests, as well as many undescribed species at each local fauna. To construct a rapid and reliable identification system for this subfamily, the effectiveness of molecular species identification was assessed based on 421 cytochrome c oxidase subunit I (COI) sequences of 84 morphologically identified species. Among the 84 morphospecies, molecular species identification of 60 species (71.4%) was consistent with their morphological identifications. Six cryptic and/or pseudocryptic species with large genetic divergence (>5%) were confirmed by their sympatric or allopatric distributions. However, 18 species, including a subspecies, had ambiguous genetic distances and shared overlapping intra- and interspecific genetic distances (range: 2.12%–3.67%) suggesting incomplete lineage sorting, introgression of mitochondrial genome, or affection by endosymbionts, such as Wolbachia infection, between species and simple genetic variation within species. In this study, we propose a conservative threshold of 3.6% for convenient molecular operational taxonomic unit (MOTU) identification in the subfamily Denticollinae based on the results of pairwise genetic distances analyses using neighbor-joining, mothur, Automatic Barcode Gap Discovery analysis, and tree-based species delimitation by Poisson Tree Processes analysis. Using the 3.6% threshold, we identified 87 MOTUs and found 8 MOTUs in the interval between 2.5% to 3.5%. Evaluation of MOTUs identified in this range requires integrative species delimitation, including review of morphological and ecological differences as well as sensitive genetic markers. From this study, we confirmed that COI sequence is useful for reassessing species diversity for polymorphic and polytypic species occurring in sympatric and allopatric distributions, and for a single species having an extensively large habitat.

Extent and divergence of heteroplasmy of the DNA barcoding region in Anapodisma miramae (Orthoptera: Acrididae)

A partial sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene is widely used as a molecular marker for species identification in animals, also termed a DNA barcode. However, the presence of more than one sequence type in a single individual, also known as heteroplasmy, is one of the shortcomings of barcode identification. In this study, we examined the extent and divergence of COI heteroplasmy, including nuclear-encoded mitochondrial pseudogenes (NUMTs), at the genomic-DNA level from 13 insect species including orthopteran Anapodisma miramae, and a long fragment of mitochondrial DNA and cDNA from A. miramae as templates. When multiple numbers of clones originated from genomic DNA were sequenced, heteroplasmy was prevalent in all species and NUMTs were observed in five species. Long fragment DNA (∼13.5 kb) also is a source of heteroplasmic amplification, but the divergent haplotypes and NUMTs obtained from genomic DNA were not detected in A. miramae. On the other hand, cDNA was relatively heteroplasmy-free. Consistently, one dominant haplotype was always obtained from the genomic DNA-origin clones in all species and also from the long fragment- and cDNA-origin clones in the two tested individuals of A. miramae. Furthermore, the dominant haplotype was identical in sequence, regardless of the DNA source in A. miramae. Thus, one possible solution to avoid the barcoding problem in relationship to heteroplasmy could be the acquisition of multiple numbers of barcoding sequences to determine a dominant haplotype that can be assigned as barcoding sequence for a given species.

Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A+T composition (79.3%-80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA(Gln) and ND2 (64-87 bp), conserved ATACTAA motif between tRNA(Ser (UCN)) and ND1, and characteristic features of the A+T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA(Asn)/tRNA(Ser (AGN)), instead of the rearrangement tRNA(Ser (AGN))/tRNA(Asn), found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea+Noctuoidea+Geometroidea)+Pyraloidea)+Papilionoidea)+Tortricoidea)+Yponomeutoidea)+Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae+(Hesperiidae+(Pieridae+((Lycaenidae+Riodinidae)+Nymphalidae)))) were strongly supported in all analyses (0.98-1 by BI and 96-100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.

Complete mitochondrial genome of the Mormon metalmark butterfly, Apodemia mormo (Lepidoptera: Riodinidae)

We report the complete mitochondrial genome (mitogenome) of Apodemia mormo in the lepidopteran family Riodinidae. The 15,262 bp long complete genome was comprised of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and 1 major non-coding A + T-rich region, with the arrangement typically found in the majority of Lepidoptera. The genes of A. mormo were interleaved with a total of 168 bp, which were spread over 16 regions and overlapped 58 bp at eight locations. All tRNAs of the A. mormo mitogenome formed a typical cloverleaf structure, except tRNA(Ser(AGN)), which formed a truncated dihydrouridine arm. The COI gene started with CGA as seen in other Lepidoptera, instead of canonical ATN. The 349 bp long A + T-rich region harbored conserved sequence blocks, such as an ATAGA motif, a poly-T stretch, a conserved ATTTA sequence, and a microsatellite A/T repeat typically found in Lepidoptera, but absent from the tRNA-like pseudogene.

The complete mitochondrial genome of the moon moth, Actias aliena (Lepidoptera: Saturniidae)

In this study, we describe the complete mitochondrial genome (mitogenome) sequence of the Actias aliena belonging to the lepidopteran family Saturniidae in terms of general genomic features and composition. The 15,243 bp long genome consisted of a typical set of genes (13 protein-coding genes, 2 rRNA genes and 22 tRNA genes) and 1 major non-coding A+T-rich region. The A. aliena mitogenome harbored the gene order tRNA(Met), tRNA(Ile) and tRNA(Gln) between the A+T-rich region and ND2, as shown in most lepidopteran species. The COI gene possessed the CGA initiator, which is found in nearly all lepidopterans lacking a canonical ATN initiator. Twenty-one tRNAs formed the cloverleaf secondary structures but tRNA(Ser)(AGN) formed a simple loop in the DHU arm. The 328 bp long A+T-rich region, which was located between the S rRNA and tRNA(Met) genes, contained several Lepidoptera-specific sequences, such as the ATAGA motif, a poly-T stretch, an AT repeat, and a poly-A stretch, along with an unusual tRNA(Phe)-like structure.