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

Comparative genomic and phylogenetic analysis of the complete mitochondrial genome of Cricula trifenestrata (Helfer) among lepidopteran insects

Cricula trifenestrata Helfer (commonly known as Amphutukoni muga/Cricula silkworm), a wild sericigenous insect produces golden yellow silk similar to Antheraea assamensis (muga silkworm), with significant potential as a natural fiber and biomaterial. Cricula is considered as a pest as it competes for food with muga, which produces the prized golden silk. This study focuses on decoding the mitochondrial genome of C. trifenestrata using next-generation sequencing technology and includes comparative analysis with Bombycoids and other lepidopteran insects. We found that the Cricula mitogenome spans 15 425 bp and exhibits typical gene content and arrangement consistent with other Saturniids and lepidopterans. All protein-coding genes were found to undergo purifying selection, with the highest and lowest conservation observed in the cox1 and atp8 gene, respectively, indicating their potential role in future evolutionary events. We identified two types of mismatches: 23 "G-U" and 6 "U-U" pairs, similar to those found in Actias selene among the Saturniids. Additionally, our study uncovered the presence of two 33 bp repeat units and a "TTAGA" motif in the control region, in contrast to the typical "ATAGA" motif, suggesting functional similarity with evolving sequences. Furthermore, phylogenetic analysis supports the close relationship of Cricula with other species within the Saturniidae family.

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

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.

Comparative analysis of mitochondrial genomes reveals family-specific architectures and molecular features in scorpions (Arthropoda: Arachnida: Scorpiones)

Scorpions are a group of arachnids with great evolutionary success that comprise more than 2,000 described species. Mitochondrial genomes have been little studied in this clade. We describe and compare different scorpion mitochondrial genomes and analyze their architecture and molecular characteristics. We assembled eight new scorpion mitochondrial genomes from transcriptomic datasets, annotated them, predicted the secondary structures of tRNAs, and compared the nucleotide composition, codon usage, and relative synonymous codon usage of 16 complete scorpion mitochondrial genomes. Lastly, we provided a phylogeny based on all mitochondrial protein coding genes. We characterized the mitogenomes in detail and reported particularities such as dissimilar synteny in the family Buthidae compared to other scorpions, unusual tRNA secondary structures, and unconventional start and stop codons in all scorpions. Our comparative analysis revealed that scorpion mitochondrial genomes exhibit different architectures and features depending on taxonomic identity. We highlight the parvorder Buthida, particularly the family Buthidae, as it invariably exhibited different mitogenome features such as synteny, codon usage, and AT-skew compared to the parvorder Iurida that included the rest of the scorpion families we analyzed in this study. Our results provide a better understanding of the evolution of mitogenome features and phylogenetic relationships in scorpions.

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

Universal Mitochondrial Multi-Locus Sequence Analysis (mtMLSA) to Characterise Populations of Unanticipated Plant Pest Biosecurity Detections

Simple Summary

Agricultural and environmental sustainability requires effective biosecurity responses that prevent the establishment or spread of exotic insect pests. Understanding where new detections may have come from or if recurrent detections are connected contributes to this. Suitable population genetic markers use relatively rapidly evolving gene regions which render the PCR method species-specific at best. Because resource limitations mean these are pre-emptively developed for the highest risk species, populations of other exotic pests are unable to be characterised at the time. Here we have developed a generic method that is useful across species within the same taxonomic Order, including where there is little or no prior knowledge of their gene sequences. Markers are formed by concomitant sequencing of four gene regions. Sequence concatenation was shown to retrieve higher resolution signatures than standard DNA barcoding. The method is encouragingly universal, as illustrated across species in ten fly and 11 moth superfamilies. Although as-yet untested in a biosecurity situation, this relatively low-tech, off-the-shelf method makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for impromptu species-specific research and development.

Abstract

Biosecurity responses to post-border exotic pest detections are more effective with knowledge of where the species may have originated from or if recurrent detections are connected. Population genetic markers for this are typically species-specific and not available in advance for any but the highest risk species, leaving other less anticipated species difficult to assess at the time. Here, new degenerate PCR primer sets are designed for within the Lepidoptera and Diptera for the 3′ COI, ND3, ND6, and 3′ plus 5′ 16S gene regions. These are shown to be universal at the ordinal level amongst species of 14 and 15 families across 10 and 11 dipteran and lepidopteran superfamilies, respectively. Sequencing the ND3 amplicons as an example of all the loci confirmed detection of population-level variation. This supported finding multiple population haplotypes from the publicly available sequences. Concatenation of the sequences also confirmed that higher population resolution is achieved than for the individual genes. Although as-yet untested in a biosecurity situation, this method is a relatively simple, off-the-shelf means to characterise populations. This makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for unprepared species-specific research and development.

The complete mitochondrial genome of the smudged eighty-eight butterfly Diaethria gabaza eupepla (Salvin & Godman, 1868) (Insecta: Lepidoptera: Nymphalidae)

The smudged eighty-eight butterfly Diaethria gabaza eupepla (Salvin & Godman, 1868) (Nymphalidae) is a vividly colored aposematic butterfly from Central and South America. A complete circular mitochondrial genome (mitogenome) of 15,156 bp from D. gabaza eupepla was assembled from a genome skimming Illumina sequence library. The AT-rich (80.5% AT) mitogenome consists of 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a control region in the typical butterfly gene order. Diaethria gabaza eupepla COX1 begins with an atypical CGA start codon and ATP6, COX1, COX2, CYTB, ND1, ND4, ND4L, and ND5 mRNAs contain incomplete stop codons completed by the addition of 3' A residues. Phylogenetic reconstruction places Diaethria as the sister clade to Hamadryas within monophyletic nymphalid subfamily Biblidinae, consistent with previous phylogenetic hypotheses.

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 of the pirate butterfly Catacroptera cloanthe (Stoll, 1781) (Insecta: Lepidoptera: Nymphalidae: Kallimini)

The pirate butterfly Catacroptera cloanthe (Stoll, 1781) (Nymphalidae: Kallimini) is a monotypic genus of butterfly that occupies grassland and savanna habitats in Sub-Saharan Africa, and exhibits seasonal variation. Genome skimming by Illumina sequencing allowed the assembly of a 78.8% AT-rich complete circular mitogenome of 15,204 bp from C. cloanthe. The mitogenome has a typical butterfly gene order consisting of 13 protein-coding genes, two rRNAs, 22 tRNAs, and a control region. Catacroptera cloanthe COX1 begins with an atypical CGA start codon, while COX2, ND3, ND4, and ND5 end with incomplete T or TA stop codons, completed by the addition of the poly-A tail during mRNA processing. Bayesian phylogenetic reconstruction placed Catacroptera cloanthe as sister to Mallika jacksoni in the monophyletic tribe Kallimini, which was consistent with previous phylogenetic hypotheses.

Phylogenetic analysis of the complete mitochondrial genome of the Blomfild's Beauty butterfly Smyrna blomfildia (Fabricius 1781) (Insecta: Lepidoptera: Nymphalidae: Nymphalini)

The Blomfild's Beauty butterfly Smyrna blomfildia (Fabricius 1781) (Lepidoptera: Nymphalidae: Nymphalini) is a sexually dimorphic species found in Mexico, Central, and South America. Males are territorial and are more vibrantly colored than females. Genome skimming by Illumina sequencing allowed the assembly of a complete circular mitochondrial genome (mitogenome) of 15,149 bp from S. blomfildia consisting of 83.9% AT nucleotides, 13 protein-coding genes, 22 tRNAs, two rRNAs, and a control region in the typical butterfly gene order. The S. blomfilda COX1 gene features an atypical start codon (CGA) while ATP6, COX1, COX2, CYTB, ND1, ND3, ND4, and ND5 display partial stop codons completed by the addition of 3' A residues to the mRNA. Bayesian phylogenetic reconstruction places Smyrna as a member of the tribe Nymphalini and sister to a clade containing genera Araschnia, Vanessa, Polygonia, and Aglais, which differs from its classic taxonomic placement in tribe Coeini.

Phylogenetic analysis of the complete mitochondrial genome of the Japanese peacock butterfly Aglais io geisha (Stichel 1907) (Insecta: Lepidoptera: Nymphalidae)

The peacock butterfly Aglais io (Linnaeus, 1758) (Nymphalidae: Nymphalinae: Nymphalini) is a colorful and charismatic flagship butterfly species whose range spans from the British Isles and Europe through temperate Asia and the Far East. In Europe, it has been used as a model species for studying the effects of GMO maize pollen on caterpillar growth and survivorship. The Japanese subspecies, Aglais io geisha (Stichel 1907), is not as well studied as its European counterpart. Genome skimming by Illumina sequencing allowed the assembly of a complete circular mitochondrial genome (mitogenome) of 15,252 bp from A. io geisha consisting of 80.6% AT nucleotides, 13 protein-coding genes, 22 tRNAs, two rRNAs, and a control region in the gene order typical of butterfly species. Aglais io geisha COX1 gene features an atypical start codon (CGA) while COX1, COX2, CYTB, ND1, ND3, ND4, and ND5 display incomplete stop codons finished by the addition of 3' A residues to the mRNA. Bayesian phylogenetic reconstruction places A. io geisha within a clade with European A. io mitogenomes in the tribe Nymphalini, which is consistent with previous phylogenetic hypotheses.

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.

Complete mitochondrial genome of the highly fecund Bombyx mori linnaeus, 1758 (Lepidoptera: Bombycidae) strain Jam 146

To meet the increasing demands of the society in the current era, new strains of the domesticated silkworm Bombyx mori Linnaeus, 1758 (Lepidoptera: Bombycidae) are being continuously bred. Consequently, cataloging the genetic information of pure lines is essential. The strain Jam 146, whose larvae have atypical pale, crescent-shaped body markings, is an important breeding resource due to its excellent fecundity. In this study, we sequenced the mitochondrial genome (mitogenome) of this strain using next-generation sequencing. The complete genome of this strain has a gene arrangement typical of Lepidoptera. The length of the Jam 146 mitogenome (15,661 bp) is well within the range reported in other B. mori strains, i.e. between 15,629 (Baiyun strain, China) and 15,676 bp (Hukpyobeom strain, South Korea). However, the total length of protein-coding genes, 3,733 codons in Jam 146 and two other silkworm strains previously reported from South Korea, is 13 codons longer than that in other B. mori strains. Phylogenetic analysis of 22 silkworm strains from nine countries showed that the Jam 146 strain forms a strong cluster with three other strains from China, Japan, and South Korea, suggesting that after their split from a common ancestor, the evolutionary divergence among the silkworm strains in these countries has been limited.

The complete mitochondrial genome of Meleonoma mirabilis (Wang 2003) (Lepidoptera: Autostichidae)

The mitogenome of Meleonoma mirabilis was determined in this study. It was 15,268 bps long and strongly AT biased. It harbored 13 PCGs, 22 tRNAs, 2 rRNAs, and 1 non-coding control region (334 bps). BI phylogenetic analysis based on 13 concatenated PCGs from 24 moth species indicated that M. mirabilis was clustered in the family Autostichidae, which was consistent with the latest phylogenetic study.

Complete mitochondrial genome of Coleophora therinella Tengström, 1848 (Lepidoptera: Coleophoridae)

The mitochondrial genome (mitogenome) of Coleophora therinella Tengström, 1848 is the first report for the family Coleophoridae in Lepidoptera. The 15,539-bp long complete genome has an arrangement identical to that observed in most lepidopteran genomes. COI had the atypical CGA codon that is frequently found in the start region of the lepidopteran COI, and COII had the GTG codon found previously in Drosophila yakuba ND5 and Rattus norvegicus ND1. The 457-bp long A + T-rich region was the second largest, next to Blastobasis lacticolella, which belongs to Blastobasidae in the superfamily Gelechioidea. The A/T content of the whole mitogenome was 80.7%; however, it varied among the regions/genes as follows: A + T-rich region, 94.8%; srRNA, 85.0%; lrRNA, 84.3%; tRNAs, 81.5%; and PCGs, 78.9%. Phylogenetic analyses with concatenated sequences of the 13 PCGs and two RNA genes using the maximum likelihood method, placed Coleophoridae, represented only by C. therinella, as the most basal lineage of the Gelechioidea families consisted of Stathmopodidae, Scythrididae, Blastobasidae, Autostichidae, and Oecophoridae, but nodal support for this grouping was very low (27%). Currently, several families of Gelechioidea are represented by a single species. Thus, extended sampling is required for further reasonable inference for the relationships of these 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.

Complete mitochondrial genome of the Small-Branded Swift: Pelopidas mathias (Lepidoptera, Hesperiidae)

We report the complete mitochondrial genome of the Small-Branded Swift: Pelopidas mathias, which is an important pest of rice. The total length of the circular double-stranded mitogenome is 15,524 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs) and a non-coding AT-rich region with the nucleotide base composition of 40.07% A, 40.83% T, 11.59% C, and 7.51% G, showing a relatively strong AT bias. The gene order and organization are consistent with typical Lepidoptera species. This work will provide molecular data support for the study of the phylogeny and evolution in the family Hesperiidae.

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.

It's a moth! It's a butterfly! It's the complete mitochondrial genome of the American moth-butterfly Macrosoma conifera (Warren, 1897) (Insecta: Lepidoptera: Hedylidae)!

The taxonomic placement of the moth-butterfly, Macrosoma conifera (Warren 1897) (Lepidoptera: Hedylidae), has been controversial. The 15,344 bp complete M. conifera circular mitogenome, assembled by genome skimming, consists of 81.7% AT nucleotides, 22 tRNAs, 13 protein-coding genes, 2 rRNAs and a control region in the typical butterfly gene order. Macrosoma conifera COX1 features an atypical CGA start codon while ATP6, COX1, COX2, and ND5 exhibit incomplete stop codons completed by the post-transcriptional addition of 3' A residues. Phylogenetic reconstruction places M. conifera as sister to the skippers (Hesperiidae), which is consistent with several recent phylogenetic analyses.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

Complete mitochondrial genome of Scythris sinensis (Lepidoptera: Scythrididae)

The complete mitochondrial genome (mitogenome) of Scythris sinensis Felder & Rogenhofer, 1775 (Lepidoptera: Scythrididae) was 15,216 bp with a typical set of genes (13 protein-coding genes [PCGs], 2 rRNA genes, and 22 tRNA genes) and one non-coding region, with an arrangement identical to that observed in most lepidopteran genomes. Twelve PCGs had the typical ATN start codon, whereas COI had the atypical CGA codon that is frequently found in the start region of the lepidopteran COI. The 271-bp long A + T-rich region was the shortest among sequenced Gelechioidea, which ranged from 290 – 375 bp. Phylogenetic analyses with concatenated sequences of the 13 PCGs, two RNA genes, and 22 tRNA genes using the Bayesian inference (BI) method, placed S. sinensis in the Scythrididae, as a sister to the family Stathmopodidae. The nodal support for this sister relationship was the highest at Bayesian posterior probabilities = 1.

Complete mitochondrial genome of Amorophaga japonica Robinson, 1986 (Lepidoptera: Tineidae)

The complete mitochondrial genome (mitogenome) of Amorophaga japonica Robinson, 1986 (Lepidoptera: Tineidae), comprises 15,027 base pairs (bp) and contains a typical set of genes (13 protein-coding genes [PCGs], 2 rRNA genes, and 22 tRNA genes), and 1 non-coding region. The genome has an arrangement, trnW-trnY-trnC, instead of typical trnW- trnC-trnY at the ND2 and COI junction. This arrangement is unique in lepidopteran mitogenomes. Unlike most lepidopteran insects, which have CGA as the start codon for the COI gene sequence, A. japonica COI had a typical ATT codon. The A + T-rich region was unusually short, with only 199 bp. Phylogenetic analyses with concatenated sequences of the 13 PCGs and two rRNA genes using the Bayesian inference method placed A. japonica in Tineidae as a sister to the cofamilial species, Tineola bisselliella, with high nodal support (Bayesian posterior probability [BPP] = 0.99), presenting the superfamily Tineoidea in a monophyletic group with a BPP of 0.99. Gracillarioidea, represented by three species of Gracillariidae, formed a monophyletic group with the highest BPP, but the Leucoptera malifoliella in Yponomeutoidea was unusually grouped together with the Gracillarioidea with the highest nodal support. As more mitogenome sequences are available, further analysis to infer the relationships among superfamilies of Lepidoptera might be possible.

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.

The complete mitochondrial genome of Amesia sanguiflua (Lepidoptera, Zygaenidae)

Amesia sanguiflua (Lepidoptera, Zygaenidae) is found in northern India, Myanmar, Indochina, Malay Peninsula, Sumatra, Java, and China. In the present study, we sequenced the complete mitochondrial genome of A. sanguiflua. The mitochondrial genome was 15,203 bp in length, containing a typical set of 37 genes (13 protein-coding genes, 2 rRNA genes, 22 tRNA genes) and a 346bp non-coding A+T-rich region. Phylogenetic analysis using mitochondrial genomes of 40 species showed that A. sanguiflua formed a well-supported monophyletic group with other Zygaenidae species.

Complete mitochondrial genome of the silkworm strain, Chilseongjam Bombyx mori (Lepidoptera: Bombycidae), with a unique larval body marking

Recently, a new silkworm strain with a peculiar larval marking and rare cocoon colour was bred in Korea for educational learning and exhibition. In order to obtain the genetic information of the newly bred strain, Chilseongjam Bombyx mori (Lepidoptera: Bombycidae), its complete mitochondrial genome (mitogenome) was sequenced. The mitogenome is 15,660 bp in length, contains a typical set of genes, and has gene arrangement and composition typical of Lepidoptera. However, the Chilseongjam strain mitogenome is 4-36 bp longer than 19 other strains originating from other countries and 16 bp shorter than the whole genome of a Korean Hukpyobeom strain. In particular, the Chilseongjam strain has an intergenic spacer sequence that is shorter than that of the Hukpyobeom strain at the tRNAHis and ND4 junction as it has fewer microsatellite-like AT repeats. Phylogenetic analyses conducted using a total of 21 silkworm strains originating from nine countries revealed a few subgroups with moderate-to-high nodal support (80-94%). The Korean Chilseongjam strain formed a relatively strong subgroup (85%) with a Japanese strain (J106) instead of the Korean Hukpyobeom strain.

Mitochondrial Recombination Reveals Mito-Mito Epistasis in Yeast

Genetic variation in mitochondrial DNA (mtDNA) provides adaptive potential although the underlying genetic architecture of fitness components within mtDNAs is not known. To dissect functional variation within mtDNAs, we first identified naturally occurring mtDNAs that conferred high or low fitness in Saccharomyces cerevisiae by comparing growth in strains containing identical nuclear genotypes but different mtDNAs. During respiratory growth under temperature and oxidative stress conditions, mitotype effects were largely independent of nuclear genotypes even in the presence of mito-nuclear interactions. Recombinant mtDNAs were generated to determine fitness components within high- and low-fitness mtDNAs. Based on phenotypic distributions of isogenic strains containing recombinant mtDNAs, we found that multiple loci contributed to mitotype fitness differences. These mitochondrial loci interacted in epistatic, nonadditive ways in certain environmental conditions. Mito-mito epistasis (i.e., nonadditive interactions between mitochondrial loci) influenced fitness in progeny from four different crosses, suggesting that mito-mito epistasis is a widespread phenomenon in yeast and other systems with recombining mtDNAs. Furthermore, we found that interruption of coadapted mito-mito interactions produced recombinant mtDNAs with lower fitness. Our results demonstrate that mito-mito epistasis results in functional variation through mitochondrial recombination in fungi, providing modes for adaptive evolution and the generation of mito-mito incompatibilities.

Complete mitochondrial genome of Saturnia jonasii (Lepidoptera: Saturniidae): Genomic comparisons and phylogenetic inference among Bombycoidea

The complete mitochondrial genome (mitogenome) of Saturnia jonasii (Lepidoptera: Saturniidae) was sequenced and compared to those of 19 other bombycoid species. Furthermore, the mitogenome sequences were used to infer phylogenetic relationships among bombycoid species. The 15,261-bp Saturnia jonasii mitogenome contained the typical sets of genes and gene arrangements found in majority of Lepidoptera. All Bombycoidea species, including Saturnia jonasii, have a 15-33-bp spacer sequence at the trnS₂-ND1 junction. The phylogenetic reconstruction of bombycoid species consistently and strongly supported monophylies of the families, Saturniidae, Bombycidae, and Sphingidae, based on Bayesian inference (BI) and maximum-likelihood (ML) methods. Among these families, the Bombycidae and Sphingidae species consistently showed a sister relationship, regardless of data partitions; the BI method strongly supported this relationship, whereas it was moderately supported using the ML method.

The mitochondrial genome of Muga silkworm (Antheraea assamensis) and its comparative analysis with other lepidopteran insects

Muga (Antheraea assamensis) is an economically important silkmoth endemic to the states of Assam and Meghalaya in India and is the producer of the strongest known commercial silk. However, there is a scarcity of genomic and proteomic data for understanding the organism at a molecular level. Our present study is on decoding the complete mitochondrial genome (mitogenome) of A. assamensis using next generation sequencing technology and comparing it with other available lepidopteran mitogenomes. Mitogenome of A. assamensis is an AT rich circular molecule of 15,272 bp (A+T content ~80.2%). It contains 37 genes comprising of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA genes along with a 328 bp long control region. Its typical tRNA^Met-tRNA^Ile-tRNA^Gln arrangement differed from ancestral insects (tRNA^Ile-tRNA^Gln-tRNA^Met). Two PCGs cox1 and cox2 were found to have CGA and GTG as start codons, respectively as reported in some lepidopterans. Interestingly, nad4l gene showed higher transversion mutations at intra-species than inter-species level. All PCGs evolved under strong purifying selection with highest evolutionary rates observed for atp8 gene while lowest for cox1 gene. We observed the typical clover-leaf shaped secondary structures of tRNAs with a few exceptions in case of tRNA^Ser1 and tRNA^Tyr where stable DHU and TΨC loop were absent. A significant number of mismatches (35) were found to spread over 19 tRNA structures. The control region of mitogenome contained a six bp (CTTAGA/G) deletion atypical of other Antheraea species and lacked tandem repeats. Phylogenetic position of A. assamensis was consistent with the traditional taxonomic classification of Saturniidae. The complete annotated mitogenome is available in GenBank (Accession No. KU379695). To the best of our knowledge, this is the first report on complete mitogenome of A. assamensis.

Sequencing of the complete mitochondrial genome of a fish-parasitic flatworm Paratetraonchoides inermis (Platyhelminthes: Monogenea): tRNA gene arrangement reshuffling and implications for phylogeny

BACKGROUND

Paratetraonchoides inermis (Monogenea: Tetraonchoididae) is a flatworm parasitising the gills of uranoscopid fishes. Its morphological characteristics are ambiguous, and molecular data have never been used to study its phylogenetic relationships, which makes its taxonomic classification controversial. Also, several decades of unsuccessful attempts to resolve the relationships within the Monogenea present a strong indication that morphological datasets may not be robust enough to be used to infer evolutionary histories. As the use of molecular data is currently severely limited by their scarcity, we have sequenced and characterized the complete mitochondrial (mt) genome of P. inermis. To investigate its phylogenetic position, we performed phylogenetic analyses using Bayesian inference and maximum likelihood approaches using concatenated amino acid sequences of all 12 protein-coding genes on a dataset containing all available monogenean mt genomes.

RESULTS

The circular mt genome of P. inermis (14,654 bp) contains the standard 36 genes: 22 tRNAs, two rRNAs, 12 protein-encoding genes (PCGs; Atp8 is missing) and a major non-coding region (mNCR). All genes are transcribed from the same strand. The A + T content of the whole genome (82.6%), as well as its elements, is the highest reported among the monogeneans thus far. Three tRNA-like cloverleaf structures were found in mNCR. Several results of the phylogenomic analysis are in disagreement with previously proposed relationships: instead of being closely related to the Gyrodactylidea, Tetraonchidea exhibit a phylogenetic affinity with the Dactylogyridea + Capsalidea clade; and the order Capsalidea is neither basal within the subclass Monopisthocotylea, nor groups with the Gyrodactylidea, but instead forms a sister clade with the Dactylogyridea. The mt genome of P. inermis exhibits a unique gene order, with an extensive reorganization of tRNAs. Monogenea exhibit exceptional gene order plasticity within the Neodermata.

CONCLUSIONS

This study shows that gene order within monopisthocotylid mt genomes is evolving at uneven rates, which creates misleading evolutionary signals. Furthermore, our results indicate that all previous attempts to resolve the evolutionary history of the Monogenea may have produced at least partially erroneous relationships. This further corroborates the necessity to generate more molecular data for this group of parasitic animals.

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.

A transfer RNA gene rearrangement in the lepidopteran mitochondrial genome

Gene arrangements in the mitochondrial genomes (mitogenomes) of insects are conserved across the major lineages, but can be rearranged within derived groups and may provide valuable phylogenetic characters. In this study, we sequenced the entire mitogenome of Parasa consocia, a moth of the family Limacodidae (Lepidoptera: Zygaenoidea). Compared with other lepidopterans and ancestral insects, the P. consocia mitogenome features a transfer RNA gene arrangement novel among lepidopterans between the ND3 and ND5 genes: RANSEF (the underline signifies an inverted gene), which differs from the ARNSEF arrangement of ancestral insects. This rearrangement can be explained by the tandem duplication-random loss model. We inferred a phylogenetic hypothesis for the lepidopteran superfamily based on mitochondrial amino-acid sequences using the Bayesian-inference and maximum-likelihood methods. Our results showed that P. consocia belongs to the Zygaenoidea superfamily and supported the following phylogenetic relationship: Yponomeutoidea + (Tortricoidea + Zygaenoidea + (Papilionoidea + (Pyraloidea + (Noctuoidea + (Geometroidea + Bombycoidea)))))). Comparative analyses indicated that mitogenomes are a useful phylogenetic tool at the subfamily level within the order Lepidoptera. Our findings also suggest that mitogenomes are likely to represent a valuable tool for systematics in other groups of lepidopterans.

The mitochondrial genome of the lepidopteran host cadaver (Thitarodes sp.) of Ophiocordyceps sinensis and related phylogenetic analysis

To understand the phylogeny of the host insect (Thitarodes sp.) of the fungus Ophiocordyceps sinensis, we sequenced, annotated and characterized the complete mitochondrial (mt) genome of the host cadaver of a natural O. sinensis. Further, we compared the Thitarodes sp. mt genome with those of the other 7 sequenced Hepialidae and examined the phylogenetic relationships using a constructed Maximum Likelihood (ML) phylogenetic tree and mt genomic features (genetic distances and intergenic spacers). The mt genome is a circular molecule of 16,280bp in length with a high A+T content (81.20%) and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and an AT-rich region. The gene arrangement is identical to the ancestral arrangement but differs from those of other lepidopteran mt genomes because of the arrangement of tRNA genes. The tRNA region, which is located between the AT-rich region and nad2, is trnI/trnQ/trnM (IQM) in Thitarodes sp., rather than the trnM/trnI/trnQ (MIQ) of the Lepidoptera-specific rearrangement. All PCGs begin with the canonical start codons ATN or NTG, except for cox1, which starts with CGA. Most PCGs terminate with the typical stop codon TAA, although some have an incomplete stop codon (T). The 1473bp AT-rich region is located between the rrnS (12S rRNA) and trnI, which is the longest sequenced in a Thitarodes mt genome to date, containing nine 112bp copies and one partial copy of a 55bp sequence. The results derived from the phylogenetic tree, the genetic distances and the intergenic spacers of the mt genome show that the host insect of O. sinensis belongs to the Thitarodes, while Endoclita signifer and Napialus hunanensis form a relatively distinct lineage from Thitarodes. The sequence and full annotation of this moth mt genome will provide more molecular information about the Exoporia within the Lepidoptera, and the clarification of its phylogeny will improve the management of this insect resource and the conservation and sustainable use of this endangered medicinal species in China.