Mitochondrial genome characteristics of two Sphingidae insects (Psilogramma increta and Macroglossum stellatarum) and implications for their phylogeny

Complete mitochondrial characteristics and phylogenetic analysis of Deilephila elpenor (Linnaeus, 1758) (Lepidoptera: Sphingidae)

Deilephila elpenor is widely distributed in countries of Asia and Central Europe, and the larva is recognized as significant agriculture pest. The complete mitochondrial genome of D. elpenor is 15,372 bp in length. It contains 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. TAA is utilized as the termination codon for most PCGs, however, cox1 and cox2 use the incomplete termination codon T, and nad3 uses TAG as the termination codon. UUA (Leu) is the most frequently used codon, GCG (Ala) and CCG (Pro) are the least frequently used codons. In addition, we selected 15 species sequences closest to this species from NCBI, and used Manduca quinquemaculata and Manduca sexta (Lepidoptera: Smerinthinae) as the outgroup. Phylogenetic analysis suggested that D. elpenor was the most closely related to genus Theretra, genus Rhagastis and Cechenena minor.

Comparative Mitochondrial Genomes between the Genera Amiota and Phortica (Diptera: Drosophilidae) with Evolutionary Insights into D-Loop Sequence Variability

To address the limited number of mitochondrial genomes (mitogenomes) in the subfamily Steganinae (Diptera: Drosophilidae), we assembled 12 complete mitogenomes for six representative species in the genus Amiota and six representative species in the genus Phortica. We performed a series of comparative and phylogenetic analyses for these 12 Steganinae mitogenomes, paying special attention to the commonalities and differences in the D-loop sequences. Primarily determined by the lengths of the D-loop regions, the sizes of the Amiota and Phortica mitogenomes ranged from 16,143-16,803 bp and 15,933-16,290 bp, respectively. Our results indicated that the sizes of genes and intergenic nucleotides (IGNs), codon usage and amino acid usage, compositional skewness levels, evolutionary rates of protein-coding genes (PCGs), and D-loop sequence variability all showed unambiguous genus-specific characteristics and provided novel insights into the evolutionary implications between and within Amiota and Phortica. Most of the consensus motifs were found downstream of the D-loop regions, and some of them showed distinct genus-specific patterns. In addition, the D-loop sequences were phylogenetically informative as the data sets of PCGs and/or rRNAs, especially within the genus Phortica.

Characterization of the Complete Mitochondrial Genome of Eight Diurnal Hawkmoths (Lepidoptera: Sphingidae): New Insights into the Origin and Evolution of Diurnalism in Sphingids

In this study, the mitochondrial genomes of 22 species from three subfamilies in the Sphingidae were sequenced, assembled, and annotated. Eight diurnal hawkmoths were included, of which six were newly sequenced (Hemaris radians, Macroglossum bombylans, M. fritzei, M. pyrrhosticta, Neogurelca himachala, and Sataspes xylocoparis) and two were previously published (Cephonodes hylas and Macroglossum stellatarum). The mitochondrial genomes of these eight diurnal hawkmoths were comparatively analyzed in terms of sequence length, nucleotide composition, relative synonymous codon usage, non-synonymous/synonymous substitution ratio, gene spacing, and repeat sequences. The mitogenomes of the eight species, ranging in length from 15,201 to 15,461 bp, encode the complete set of 37 genes usually found in animal mitogenomes. The base composition of the mitochondrial genomes showed A+T bias. The most commonly used codons were UUA (Leu), AUU (Ile), UUU (Phe), AUA (Met), and AAU (Asn), whereas GCG (Ala) and CCG (Pro) were rarely used. A phylogenetic tree of Sphingidae was constructed based on both maximum likelihood and Bayesian methods. We verified the monophyly of the four current subfamilies of Sphingidae, all of which had high support. In addition, we performed divergence time estimation and ancestral character reconstruction analyses. Diurnal behavior in hawkmoths originated 29.19 million years ago (Mya). It may have been influenced by the combination of herbaceous flourishing, which occurred 26-28 Mya, the uplift of the Tibetan Plateau, and the large-scale evolution of bats in the Oligocene to Pre-Miocene. Moreover, diurnalism in hawkmoths had multiple independent origins in Sphingidae.

Description of mitochon genome and phylogenetic considerations of Sabethes bipartipes, Sabethes cyaneus, Sabethes quasicyaneus, and Sabethes tarsopus (Diptera: Culicidae)

The genus Sabethes (Diptera: Culicidae) comprises species of great epidemiological relevance, particularly involved in transmission cycles of the Yellow fever virus in South America. Given the unavailability of information related to aspects of evolutionary biology and molecular taxonomy of species of this genus of mosquitoes, we report here the first sequencing of the mitochondrial genomes of Sabethes bipartipes, Sabethes cyaneus, Sabethes tarsopus, and Sabethes quasicyaneus. The sequences obtained showed an average length of 14,920 bp, comprising 37 functional genes (13 PCGs, 22 tRNA, and 02 rRNA). The phylogenies reconstructed by Maximum likelihood and Bayesian inference methods, based on the concatenated sequences of all 13 PCGs, produced similar topologies and strongly supported the monophyletic relationship between the Sabethes subgenera, corroborating the known taxonomic classification based on aspects of the external morphology of the taxa assessed. The data and information produced from the Sabethes species evaluated here may be useful for future taxonomic and evolutionary studies of the genus, as well as the Culicidae family.

Description of the mitogenome and phylogeny of Aedes spp. (Diptera: Culicidae) from the Amazon region

The genus Aedes (Diptera: Culicidae) includes species of great epidemiological relevance, particularly involved in transmission cycles of leading arboviruses in the Brazilian Amazon region, such as the Zika virus (ZIKV), Dengue virus (DENV), Yellow fever virus (YFV), and Chikungunya virus (CHIKV). We report here the first putatively complete sequencing of the mitochondrial genomes of Brazilian populations of the species Aedes albopictus, Aedes scapularis and Aedes serratus. The sequences obtained showed an average length of 14,947 bp, comprising 37 functional subunits, typical in animal mitochondria (13 PCGs, 22 tRNA, and 2 rRNA). The phylogeny reconstructed by Maximum likelihood method, based on the concatenated sequences of all 13 PCGs produced at least two non-directly related groupings, composed of representatives of the subgenus Ochlerotatus and Stegomyia of the genus Aedes. The data and information produced here may be useful for future taxonomic and evolutionary studies of the genus Aedes, as well as the Culicidae family.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region

The Culicidae family is distributed worldwide and comprises about 3587 species subdivided into the subfamilies Anophelinae and Culicinae. This is the first description of complete mitochondrial DNA sequences from Aedes fluviatilis, Aedeomyia squamipennis, Coquillettidia nigricans, Psorophora albipes, and Psorophora ferox. The mitogenomes showed an average length of 15,046 pb and 78.02% AT content, comprising 37 functional subunits (13 protein coding genes, 22 tRNAs, and two rRNAs). The most common start codons were ATT/ATG, and TAA was the stop codon for all PCGs. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogeny was inferred by analyzing the 13 PCGs concatenated nucleotide sequences of 48 mitogenomes. Maximum likelihood and Bayesian inference analysis placed Ps. albipes and Ps. ferox in the Janthinosoma group, like the accepted classification of Psorophora genus. Ae. fluviatilis was placed in the Aedini tribe, but was revealed to be more related to the Haemagogus genus, a result that may have been hampered by the poor sampling of Aedes sequences. Cq. nigricans clustered with Cq. chrysonotum, both related to Mansonia. Ae. squamipennis was placed as the most external lineage of the Culicinae subfamily. The yielded topology supports the concept of monophyly of all groups and ratifies the current taxonomic classification.

Complete mitochondrial genome of Hypomecis punctinalis Scopoli, 1763 and its phylogenetic position within family Geometridae

Hypomecis punctinalis Scopoli, 1763 belongs to the Lepidopteran family Geometridae. We sequenced the complete mitochondrial genome (mitogenome) of H. punctinalis. The mitogenome is 15,648 bp long and contains a typical set of genes (13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes) and a 484 bp AT-rich region. All PCGs start with ATN codons and stop at TAA codon except for cox1 using CGA as initiation codon and nad4 and nad5 using incomplete termination codon T. Within the mitogenome, 17 intergenic spacers and seven overlaps are founded. The intergenic nucleotides are 294 bp in total and two longest intervals locate between trnGln and nad2 as well as trnCys and trnTyr . The overlap nucleotides are 47 bp in total and the maximum overlap lies between cox2 and trnLys . The AT-rich region of the mitogenome contains an 'ATAGA + polyT' motif, three copies of 30-bp-repeat and a short polyA tail. The phylogenetic tree shows the relationships of four subfamilies of Geometridae are (((Ennominae + Geometrinae)+Larentiinae)+Sterrhinae)) and the relationships within subfamily Ennominae are ((((Erannis+Biston)+(Jankowskia+(Hypomecis+(Apocheima+Milionia))))+Ectropis)+Abraxas)+Phthonandria)+Celenna).

The complete mitogenome sequence of the hawk moth, Theretra latreillii subsp. lucasii (Lepidoptera: Sphingidae) from Zhejiang Province, China

The sphingid, Theretra latreillii subsp. lucasii is a common hawk moth distributed in southeast Asia and Australian regions. Although barcode analyses have been published, its complete mitogenome sequence has not been deciphered. In this study, the complete mitogenome of T. latreillii lucasii (GeneBank accession no. MW539688) was sequenced using Illumina HiSeq X Ten system for mitogenome-based phylogenetic analysis. The mitogenome was 15,354 bp in length and comprises 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 22 transfer RNAs (tRNAs) with the typical gene order and orientation of Sphingidae mitogenomes. The nucleotide composition of majority strand is 41.2% for A, 7.4% for G, 12.0% for C, and 39.4% for T, with an A + T content of 80.6%. Phylogenetic analysis using the 13 PCGs fully resolved T. latreillii lucasii in a clade with T. japonica, Macroglossum stellatarum, and Ampelophaga rubiginosa, with high nodal support both by Bayesian inference and maximum-likelihood methods, forming the Macroglossini monophyletic group.

First mitogenome of subfamily Langiinae (Lepidoptera: Sphingidae) with its phylogenetic implications

To date, a relatively complete classification of Sphingidae (Lepidoptera) has been generated, but the phylogeny of the family remains need to be fully resolved. Some phylogenetic relationships within Sphingidae still remains uncertain, especially the taxonomic status of the subfamily Langiinae and its sole included genus and species, Langia zenzeroides. To begin to address this problem, we generated nine new complete mitochondrial genomes, including that of Langia, and together with that of Theretra oldenlandiae from our previous study and 25 other Sphingidae mitogenomes downloaded from GenBank, analyzed the phylogenetic relationships of Sphingidae and investigated the mitogenomic differences among members of the Langiinae, Sphinginae, Smerinthinae and Macroglossinae. The mitogenomes of Sphingidae varied from 14995 bp to 15669 bp in length. The gene order of all newly sequenced mitogenomes was identical, containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and the A + T-rich region. Nucleotide composition was A + T biased, and all the protein-coding genes exhibited a positive AT-skew, which was reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region was comprised of nonrepetitive sequences, which contained regulatory elements related to the control of replication and transcription. We analyzed concatenated gene sequences, with third codon positions of protein coding genes and rRNAs excluded, using Maximum Likelihood and Bayesian Inference techniques. All four currently recognized subfamilies were recovered as monophyletic but in contrast to the most recent studies, our preferred tree placed Langiinae as the first subfamily to diverge within Sphingidae rather as sister to Smerinthinae + Sphinginae. Our results also support the removal of the genus Barbourion from the smerinthine tribe Ambulycini to an unresolved position in "Smerinthinae incertae sedis".

Complete mitochondrial genome of Cerace xanthocosma and its phylogenetic position in the family Tortricidae

Cerace xanthocosma Diakonoff, 1950 belongs to the Tortricidae of Lepidoptera, distributed in China and Japan. Its mitochondrial genome is 15,344 bp in size, containing 37 genes (13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes) and a non-coding A + T-rich region. In the 37 genes, 14 genes are located on the minority-strand (N-strand) with the remaining genes on the majority-strand (J-strand). The A + T-rich region has a poly-T stretch and a motif ATAGA. Phylogenetic analysis using Bayesian Inference method shows the sister relationship between Tortricinae and Olethreutinae with the placement of C. xanthocosma as the basal lineage of Tortricinae.

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 complete mitochondrial genome sequence of Scolopendra mutilans L. Koch, 1878 (Scolopendromorpha, Scolopendridae), with a comparative analysis of other centipede genomes

Scolopendramutilans L. Koch, 1878 is an important Chinese animal with thousands of years of medicinal history. However, the genomic information of this species is limited, which hinders its further application. Here, the complete mitochondrial genome (mitogenome) of S.mutilans was sequenced and assembled by next-generation sequencing. The genome is 15,011 bp in length, consisting of 13 protein-coding genes (PCGs), 14 tRNA genes, and two rRNA genes. Most PCGs start with the ATN initiation codon, and all PCGs have the conventional stop codons TAA and TAG. The S.mutilans mitogenome revealed nine simple sequence repeats (SSRs), and an obviously lower GC content compared with other seven centipede mitogenomes previously sequenced. After analysis of homologous regions between the eight centipede mitogenomes, the S.mutilans mitogenome further showed clear genomic rearrangements. The phylogenetic analysis of eight centipedes using 13 conserved PCG genes was finally performed. The phylogenetic reconstructions showed Scutigeromorpha as a separate group, and Scolopendromorpha in a sister-group relationship with Lithobiomorpha and Geophilomorpha. Collectively, the S.mutilans mitogenome provided new genomic resources, which will improve its medicinal research and applications in the future.

The complete mitochondrial genome of Dysgonia stuposa (Lepidoptera: Erebidae) and phylogenetic relationships within Noctuoidea

To determine the Dysgonia stuposa mitochondrial genome (mitogenome) structure and to clarify its phylogenetic position, the entire mitogenome of D. stuposa was sequenced and annotated. The D. stuposa mitogenome is 15,721 bp in size and contains 37 genes (protein-coding genes, transfer RNA genes, ribosomal RNA genes) usually found in lepidopteran mitogenomes. The newly sequenced mitogenome contained some common features reported in other Erebidae species, e.g., an A+T biased nucleotide composition and a non-canonical start codon for cox1 (CGA). Like other insect mitogenomes, the D. stuposa mitogenome had a conserved sequence 'ATACTAA' in an intergenic spacer between trnS2 and nad1, and a motif 'ATAGA' followed by a 20 bp poly-T stretch in the A+T rich region. Phylogenetic analyses supported D. stuposa as part of the Erebidae family and reconfirmed the monophyly of the subfamilies Arctiinae, Catocalinae and Lymantriinae within Erebidae.

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