Characterization of the complete mitochondrial genome of tea tussock moth, Euproctis pseudoconspersa (Lepidoptera: Lymantriidae) and its phylogenetic implications

Complete Mitochondrial Genome Characterization of Schrankia costaestrigalis (Insecta: Erebidae: Hypenodinae) and Its Phylogenetic Implication

The pinion-streaked snout Schrankia costaestrigalis is a new potato pest that has recently been recorded in China. In this study, we analyzed the complete mitochondrial genome of S. costaestrigalis. The results revealed the mitogenome (GenBank: OQ181231) to occur as a circular DNA molecule of 16,376 bp with 51.001% AT content, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 1 control region. Notably, the PCGs exhibited typical ATN (Met) start codons, including cox1, which deviated from the usual CGA start codon observed in other lepidopteran mitogenomes, and followed the conventional TAN stop codons. The 22 tRNA genes demonstrated the ability to form a cloverleaf structure, with the exception of trnS1-NCU, which lacked the DHU arm present in other Erebidae mitogenomes. Additionally, conserved motifs like "ATAGA + poly-T (19 bp) stretch" and five microsatellite-like elements (TA) were identified in the AT-rich region. The phylogenetic trees revealed that the Hypenodinae subfamily forms an independent lineage closely related to Erebinae and Catocalinae. The comprehensive mitogenome of S. costaestrigalis will greatly enhance future studies focused on the molecular classification and phylogenetic understanding of the Hypenodinae subfamily within the larger family Erebidae.

Characterization of four mitochondrial genomes of Crambidae (Lepidoptera, Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the "PS clade." Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the "non-PS clade" in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.

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

A supertree of Northern European macromoths

Ecological and life-history data on the Northern European macromoth (Lepidoptera: Macroheterocera) fauna is widely available and ideal for use in answering phylogeny-based research questions: for example, in comparative biology. However, phylogenetic information for such studies lags behind. Here, as a synthesis of all currently available phylogenetic information on the group, we produce a supertree of 114 Northern European macromoth genera (in four superfamilies, with Geometroidea considered separately), providing the most complete phylogenetic picture of this fauna available to date. In doing so, we assess those parts of the phylogeny that are well resolved and those that are uncertain. Furthermore, we identify those genera for which phylogenetic information is currently too poor to include in such a supertree, or entirely absent, as targets for future work. As an aid to studies involving these genera, we provide information on their likely positions within the macromoth tree. With phylogenies playing an ever more important role in the field, this supertree should be useful in informing future ecological and evolutionary studies.

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.

Characterization of the complete mitochondrial genome of Orthaga olivacea Warre (Lepidoptera Pyralidae) and comparison with other Lepidopteran insects

Orthaga olivacea Warre (Lepidoptera: Pyralidae) is an important agricultural pest of camphor trees (Cinnamomum camphora). To further supplement the known genome-level features of related species, the complete mitochondrial genome of Orthaga olivacea is amplified, sequenced, annotated, analyzed, and compared with 58 other species of Lepidopteran. The complete sequence is 15,174 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a putative control region. Base composition is biased toward adenine and thymine (79.02% A+T) and A+T skew are slightly negative. Twelve of the 13 PCGs use typical ATN start codons. The exception is cytochrome oxidase 1 (cox1) that utilizes a CGA initiation codon. Nine PCGs have standard termination codon (TAA); others have incomplete stop codons, a single T or TA nucleotide. All the tRNA genes have the typical clover-leaf secondary structure, except for trnS(AGN), in which dihydrouridine (DHU) arm fails to form a stable stem-loop structure. The A+T-rich region (293 bp) contains a typical Lepidopter motifs 'ATAGA' followed by a 17 bp poly-T stretch, and a microsatellite-like (AT)13 repeat. Codon usage analysis revealed that Asn, Ile, Leu2, Lys, Tyr and Phe were the most frequently used amino acids, while Cys was the least utilized. Phylogenetic analysis suggested that among sequenced lepidopteran mitochondrial genomes, Orthaga olivacea Warre was most closely related to Hypsopygia regina, and confirmed that Orthaga olivacea Warre belongs to the Pyralidae family.

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

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

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

Mitochondrial genome characteristics of Somena scintillans (Lepidoptera: Erebidae) and comparation with other Noctuoidea insects

In this study, mitogenome of Somena scintillans (Lepidoptera: Erebidae) were sequenced and compared with other Noctuoidea species. The mitogenome is 15,410 base pairs in length. All 13 protein-coding genes (PCGs) are initiated by ATN codons except cox1 with CGA and all of PCGs terminate with TAA except nad4 with TAG. The codons ACG and CGC are absent. All the tRNA genes could be folded into the typical cloverleaf secondary structure except the trnS1 which not only loses dihydrouridine (DHU) arm but also mutates its anticodon into TCT. In the AT-rich region of the mitogenome the motif 'ATAGA' mutates to 'ATATA' and two copies of 161 bp-tandem repeats and two 'TA' short tandem repeats are founded. Phylogenetic analyses showed that S. scintillans is clustered into subfamily Lymatriinae. The phylogenetic relationships within Noctuoidea is (((Nolidae + (Euteliidae + Noctuidae)) + Erebidae) + Notodontidae).

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

The complete mitochondrial genome of Euproctis similis (Lepidoptera: Noctuoidea: Erebidae) and phylogenetic analysis

The mitochondrial genome (mitogenome) can provide information for phylogenetic analyses and evolutionary biology. We sequenced, annotated, and characterized the mitogenome of Euproctis similis. The complete mitogenome is 15,437bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a control region (A+T-rich region). The A+T content in the mitogenome was 80.16%. All PCGs use standard ATN as a start codon, with the exception of cytochrome c coxidase 1 (cox1) with CGA. A gene rearrangement (trnM) was found. All transfer RNA (tRNA) genes have a typical clover-leaf structure except for trnS1 (AGN). Phylogenetic analysis was performed using Bayesian Inference and Maximum Likelihood based on the amino acid and nucleotide sequences of 13 mitochondrial PCGs. The well-supported phylogenetic relationships can be generally described as: Notodontidae+(Erebidae+(Nolidae+(Euteliidae+Noctuidae))). The tree support that E. similis shares a close ancestry with Erebidae insects. Our results indicate that Erebidae is a sister group to the other families (Euteliidae, Nolidae, and Noctuidae).

Characterization of the Complete Mitochondrial Genome of Leucoma salicis (Lepidoptera: Lymantriidae) and Comparison with Other Lepidopteran Insects

The complete mitochondrial genome (mitogenome) of Leucoma salicis (Lepidoptera: Lymantriidae) was sequenced and annotated. It is a circular molecule of 15,334 bp, containing the 37 genes usually present in insect mitogenomes. All protein-coding genes (PCGs) are initiated by ATN codons, other than cox1, which is initiated by CGA. Three of the 13 PCGs had an incomplete termination codon, T or TA, while the others terminated with TAA. The relative synonymous codon usage of the 13 protein-coding genes (PCGs) was consistent with those of published lepidopteran sequences. All tRNA genes had typical clover-leaf secondary structures, except for the tRNASer (AGN), in which the dihydrouridine (DHU) arm could not form a stable stem-loop structure. The A + T-rich region of 325 bp had several distinctive features, including the motif 'ATAGA' followed by an 18 bp poly-T stretch, a microsatellite-like (AT)7 element, and an 11-bp poly-A present immediately upstream of tRNAMet. Relationships among 32 insect species were determined using Maximum Likelihood (ML), Neighbor Joining (NJ) and Bayesian Inference (BI) phylogenetic methods. These analyses confirm that L. salicis belongs to the Lymantriidae; and that Lymantriidae is a member of Noctuoidea, and is a sister taxon to Erebidae, Nolidae and Noctuidae, most closely related to Erebidae.

The complete mitochondrial genome of Thitarodes sejilaensis (Lepidoptera: Hepialidae), a host insect of Ophiocordyceps sinensis and its implication in taxonomic revision of Hepialus adopted in China

The mitochondrial genome is widely used for phylogenetic analyses and evolutionary biology. The complete mitochondrial genome of Thitarodes sejilaensis (Lepidoptera: Hepialidae) was sequenced and analyzed in this study. This mitogenome is a typical circular molecule of 15,290bp, with the gene content, orientation and order identical to other insects in the family Hepialidae. The genome nucleotide composition is heavily biased towards As and Ts, accounting for 80.87% of total nucleotide content. The major strand shows a positive AT-skew and negative GC-skew. All 13 protein-coding genes (PCG) are initiated by the canonical putative start codons ATN, except for COI and ND1 that use the initiation codons CGA and TTG, respectively. Nine PCGs share the complete termination codons TAA, while the remaining PCGs use an incomplete termination codon T. Additionally, the codon distribution and Relative Synonymous Codon Usage of all PCGs in the T. sejilaensis mitogenome are consistent with other Hepialidae mitogenomes. Among 22 transfer RNAs, 21 have the typical clover-leaf structure, while tRNASer(AGN) does not possess the dihydrouridine (DHU) arm and could not form a stable stem-loop structure. The secondary structures of 2 ribosomal RNA genes broadly conform to the proposed models of these genes documented in other lepidopteran insects. T. sejilaensis AT-rich region exhibits three repetitive sequences of 118bp. Other regions contain 22-bp overlapping nucleotides and 72-bp intergenic nucleotides. The phylogenetic relationships were constructed by two datasets, the amino acid sequence derived from protein-coding genes and the nucleotide sequence of 13 PCGs and 2 rRNAs. Using Maximum Likelihood (ML), we reconstructed a phylogenetic tree which supported a more primitive taxa of Hepialoidea within Lepidoptera. Moreover, according to comparisons based on the CytB sequences and morphological characteristic, Hepialus species reported in China should be revised. Our taxonomic recommendations include assigning these species to the following genera: Thitarodes, Ahamus, Hepialus and Parahepialus.