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

Complete Mitochondrial Genome and Phylogenetic Position of Nurudea zhengii Ren (Insecta, Hemiptera, Aphididae)

Nurudea zhengii Ren was identified by aphid morphological characteristics as well as the gall shape and host plant species, and placed in the tribe Fordini (Hemiptera, Aphididae, Eriosomatinae). Here, its whole genome was firstly sequenced by a genome-skimming method and its mitochondrial genome (mitogenome) was assembled to examine its genetic variation and phylogenetic position. The complete mitogenome of Nurudea zhengii is 15,392 bp in length, and consists of 13 protein-coding genes, 22 tRNAs, two rRNAs and one D-loop region. The gene order follows the mitogenomes of the other Rhus gall aphids, and similarly has an AT bias with the content of 83.9%. The majority strand is A-skewed and C-skewed, and shows opposite skewness for G-skewed in the minority strands. The ratios of nonsynonymous to synonymous substitution rates of protein-coding genes are lower than one except for ATP8, which indicated that ATP8 was undergoing positive selection. Phylogenetic analysis among the Rhus gall aphids based on 13 protein-coding genes and two rRNA genes showed that N. zhengii was sister to N. shiraii, and then clustered with N. yanoniella as a group with high support value. The two species, N. shiraii and N. yanoniella, share the same host plant Rhus chinensis, while the host of N. zhengii is R. hypoleuca. However, the phylogenetic relationship indicated that the taxa sharing the same host plant were not absolutely clustered as the closest taxa at least at species level.

The Complete Mitochondrial Genome of Mytilisepta virgata (Mollusca: Bivalvia), Novel Gene Rearrangements, and the Phylogenetic Relationships of Mytilidae

The circular mitochondrial genome of Mytilisepta virgata spans 14,713 bp, which contains 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. Analysis of the 13 PCGs reveals that the mitochondrial gene arrangement of Mytilisepta is relatively conserved at the genus level. The location of the atp8 gene in Mytilisepta keenae differs from that of other species. However, compared with the putative molluscan ancestral gene order, M. virgata exhibits a high level of rearrangement. We constructed phylogenetic trees based on concatenated 12 PCGs from Mytilidae. As a result, we found that M. virgata is in the same clade as other Mytilisepta spp. The result of estimated divergence times revealed that M. virgata and M. keenae diverged around the early Paleogene period, although the oldest Mytilisepta fossil was from the late or upper Eocene period. Our results provide robust statistical evidence for a sister-group relationship within Mytilida. The findings not only confirm previous results, but also provide valuable insights into the evolutionary history of Mytilidae.

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.

Characterization of Seventeen Complete Mitochondrial Genomes: Structural Features and Phylogenetic Implications of the Lepidopteran Insects

Lepidoptera (moths and butterflies) are widely distributed in the world, but high-level phylogeny in Lepidoptera remains uncertain. More mitochondrial genome (mitogenome) data can help to conduct comprehensive analysis and construct a robust phylogenetic tree. Here, we sequenced and annotated 17 complete moth mitogenomes and made comparative analysis with other moths. The gene order of trnM-trnI-trnQ in 17 moths was different from trnI-trnQ-trnM of ancestral insects. The number, type, and order of genes were consistent with reported moths. The length of newly sequenced complete mitogenomes ranged from 14,231 bp of Rhagastis albomarginatus to 15,756 bp of Numenes albofascia. These moth mitogenomes were typically with high A+T contents varied from 76.0% to 81.7% and exhibited negative GC skews. Among 13 protein coding genes (PCGs), some unusual initiations and terminations were found in part of newly sequenced moth mitogenomes. Three conserved gene-overlapping regions and one conserved intergenic region were detected among 17 mitogenomes. The phylogenetic relationship of major superfamilies in Macroheterocera was as follows: (Bombycoidea + Lasiocampoidea) + ((Drepanoidea + Geometroidea) + Noctuoidea)), which was different from previous studies. Moreover, the topology of Noctuoidea as (Notodontidae + (Erebidae + Noctuidae)) was supported by high Bayesian posterior probabilities (BPP = 1.0) and bootstrapping values (BSV = 100). This study greatly enriched the mitogenome database of moth and strengthened the high-level phylogenetic relationships of Lepidoptera.

Complete mitochondrial genome analysis and molecular phylogenetic implications of Kennelia xylinana (Lepidoptera: Tortricidae)

Kennelia is a small genus in Tortricidae that is distributed in the Oriental and Palaearctic regions, and its taxonomic position within the subfamily Olethreutinae is controversial. For a comprehensive understanding of the genus, we sequenced the mitogenome of Kennelia xylinana, the type species of Kennelia, and Ancylis unculana, a species of Enarmoniini; analyzed the mitogenome characteristics of K. xylinana; and explored its phylogenetic position. Similar to other members of Lepidoptera, the mitogenome of K. xylinana is 15,762-bp long and consists of 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and a noncoding control region. In particular, we found a structure (TATAATTAATAA)₁₁ in the middle of the AT-rich region. Based on the Bayesian inference and maximum likelihood analyses of the 13 PCGs of 40 tortricid species, representing 8 tribes of 2 subfamilies, K. xylinana was clustered with two members of Enarmoniini, A. unculana and Loboschiza koenigiana, and formed highly supported monophyly. The results indicate that Kennelia should be placed in the tribe Enarmoniini.

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

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.

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

BACKGROUND

Mitochondrial genomes (mitogenomes) have greatly improved our understanding of the backbone phylogeny of Lepidoptera, but few studies on comparative mitogenomics below the family level have been conducted. Here, we generated 13 mitogenomes of eight tortricid species, reannotated 27 previously reported mitogenomes, and systematically performed a comparative analysis of nucleotide composition, gene variation and phylogenetic performance.

RESULTS

The lengths of completely sequenced mitogenomes ranged from 15,440 bp to 15,778 bp, and the gene content and organization were conserved in Tortricidae and typical for Lepidoptera. Analyses of AT-skew and GC-skew, the effective number of codons and the codon bias index all show a base bias in Tortricidae, with little heterogeneity among the major tortricid groups. Variations in the divergence rates among 13 protein-coding genes of the same tortricid subgroup and of the same PCG among tortricid subgroups were detected. The secondary structures of 22 transfer RNA genes and two ribosomal RNA genes were predicted and comparatively illustrated, showing evolutionary heterogeneity among different RNAs or different regions of the same RNA. The phylogenetic uncertainty of Enarmoniini in Tortricidae was confirmed. The synonymy of Bactrini and Olethreutini was confirmed for the first time, with the representative Bactrini consistently nesting in the Olethreutini clade. Nad6 exhibits the highest phylogenetic informativeness from the root to the tip of the resulting tree, and the combination of the third coding positions of 13 protein-coding genes shows extremely high phylogenetic informativeness.

CONCLUSIONS

This study presents 13 mitogenomes of eight tortricid species and represents the first detailed comparative mitogenomics study of Tortricidae. The results further our understanding of the evolutionary architectures of tortricid mitogenomes and provide a basis for future studies of population genetics and phylogenetic investigations in this group.

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.

Alpha and beta diversity patterns of macro-moths reveal a breakpoint along a latitudinal gradient in Mongolia

Little is known about the diversity and distribution patterns of moths along latitudinal gradients. We studied macro-moths in Mongolia along an 860 km latitudinal climatic gradient to gain knowledge on community composition, alpha, beta, and gamma diversity as well as underlying factors, which can be used as baseline information for further studies related to climate change. We identified 236 species of moths of ten families. Our study shows that the diversity of moths increased with the latitude, i.e., low species richness in the south and higher richness in the north. Moth community composition changed along the gradient, and we revealed a breakpoint of beta diversity that divided grassland and desert communities. In the desert, beta diversity was driven by species loss (i.e., nestedness), and few tolerant species existed with high abundance. In contrast, in the grassland, beta diversity was driven by species replacement with more unique species, (i.e., species which occurred only in one site). We found the lowest species diversity in the transitional zones dominated by few generalist species such as Agrotis ripae and Anarta trifolii. Low precipitation and an increasing number of grazing goats are drivers of species loss. We suggest different conservation strategies regarding the contrasting patterns of beta diversity in desert and grassland.

The complete mitochondrial genome of Lyssa zampa (Lepidoptera: Uraniidae)

The complete mitochondrial genome (mitogenome) of Lyssa zampa was first reported. It is 15,314 bp in length (GenBank accession number: MW435592) and consists of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes. The nucleotide composition is A (41.5%), C (11.1%), G (7.4%), and T (40.0%). Based on the sequences of complete mitogenome from 12 geometroid species and three drepanoid species as ingroups, and two noctuoid species as outgroups, the phylogenetic tree was constructed. The results showed that the closest relationship between Uraniidae and Epicopeiidae was strongly supported by Bayesian posterior probabilities values of 0.99.

Complete mitochondrial genome of Spodoptera littoralis (Lepidoptera: Noctuidae) from Egypt

The polyphagous cotton leafworm (Spodoptera littoralis) is one of the most destructive herbivorous insects worldwide. The present study reports the complete mitochondrial genome of S. littoralis collected from Egypt. The circular-mapping mitogenome was 15,408 bp in length with an overall A + T content of 81.1%, encoding a common set of 37 genes, including 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes. Most PCGs were found to use conventional ATN as the start codon and TAN as the stop codon. The phylogenetic tree based on the nucleic acid sequences of 13 shared PCGs of 29 Noctuidae species revealed that S. littoralis and Spodoptera litura are sister species. The data in this study will be helpful to understand geographical genetic variations, phylogenetic relationships, and species identification of S. littoralis.

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

The first mitochondrial genome for Phaudidae (Lepidoptera) with phylogenetic analyses of Zygaenoidea

Phauda flammans Walker belongs to Phaudidae (Lepidoptera), which is a holometabolous and leaf-eating pest that harms trees. So far, there is no mitochondrial (mt) genome reported of Phaudidae. Herein, we sequenced and annotated the complete mt genome of P. flammans representing the first mt genome of Phaudidae and predicted the secondary structures of its RNAs in this study. This mt genome is 15470 bp long consisting of 13 protein coding genes (PCGs), 22 tRNAs, 2 rRNAs and the control region, which are usually conserved in insects. Most PCGs used the standard ATN start codons and complete TAA/TAG termination codons. Almost all of tRNA genes exhibited cloverleaf secondary structures except that the dihydorouridine (DHU) arm of tRNA^Ser(AGN) was absent. The phylogenetic analyses using both Bayesian inference (BI) and maximum likelihood (ML) methods all supported that Phaudidae was a single family being the sister group to Zygaenidae. More mt genomes are needed to better understand the phylogenetic relationships within Zygaenoidea in the future.

Characteristic and variability of five complete aphid mitochondrial genomes: Aphis fabae mordvilkoi, Aphis craccivora, Myzus persicae, Therioaphis tenera and Appendiseta robiniae (Hemiptera; Sternorrhyncha; Aphididae)

The complete mitochondrial genomes of aphids Aphis fabae mordvilkoi, A. craccivora, Myzus persicae from Aphidinae as well as Therioaphis tenera and Appendiseta robiniae from Calaphidinae were sequenced and compared with the genomes of other aphid species. A. fabae mordvilkoi, Th. tenera and A. robiniae mitogenomes were sequenced and analyzed for the first time. The annotation of A. craccivora and M. persicae were corrected compared to what was previously published. According to our data there is no translocation of tRNA-Tyr gene in A. craccivora mitogenome and this aphid species has an ancestral type of mitochondrial gene order. A + T content in all 5 mitogenomes was higher than 80%. A + T content in the Th. tenera CR was 59.5% which is untypically low. CRs of all 5 studied mitogenomes had 2 conserved motifs at their ends and extended G + C rich region. A. craccivora, M. persicae and Th. tenera had large tandem repeats inside the CRs. Detailed molecular analysis of all 5 aphid mitochondrial genomes showed the importance of a deep understanding of the molecular organization of all the functional regions of the mitochondrial DNA, which helps to avoid mistakes during genome annotation.

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.