Comparative mitogenomic analysis of damsel bugs representing three tribes in the family Nabidae (Insecta: Hemiptera)

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

The larval, pupal and mitogenomic characteristics of Agrilusadelphinus Kerremans, 1895 (Coleoptera, Buprestidae) from China

In this study, the larva and pupa of Agrilusadelphinus are described and illustrated. DNA barcoding (COI gene) was used to associate the larval and pupal stages with adults based on the maximum-likelihood method. In the resulting phylogenetic tree, species from the same species-group were found to be clustered on a branch with high support value. To better understand A.adelphinus, the complete mitochondrial genome of this species was also sequenced and annotated. Comparing this genome to the known mitogenomes of Agrilus species, the newly sequenced genome is shorter, with 15,732 bp. However, its whole mitogenome composition and gene orientation were consistent with that of most species of Buprestidae. In the mitogenome of A.adelphinus, the ATGATAG sequence was observed between ATP8 and ATP6, which is ATGATAA in other insect mitogenomes. Leu2, Phe, Ile, Gly, and Ser2 were the five most frequently encoded amino acids. The results further prove that DNA barcoding can remove the limitation of traditional taxonomy which cannot identify to species all developmental stages. This study also provides valuable molecular and morphological data for species identification and phylogenetic analyses of the genus Agrilus.

Comparative analysis of mitochondrial genomes among the family Peltoperlidae (Plecoptera: Systellognatha) and phylogenetic implications

Nowadays, the position of Peltoperlidae in Systellognatha has been resolved based on morphological analyses. However, there are different opinions based on molecular data. To date, only three peltoperlid mitogenomes are available, and more sampling is needed to obtain precise phylogenetic relationships. In this study, we obtained the complete mitogenomes of Cryptoperla kawasawai (15,832 bp) and Peltoperlopsis sagittata (15,756 bp). Our results show that gene content, gene order, DmTTF binding site, nucleotide composition, codon usage, ribonucleic acid (RNA) structure, and structural elements in the control region are highly conserved in peltoperlids. Heatmap analysis of codon usage shows that the AT-rich codons UUA, AUU, UUU, and AUA were commonly used codons in the Peltoperlidae. Evolutionary rate analyses of protein-coding genes reveal that different genes have been subject to different rates of molecular evolution correlated with the GC content. All tRNA genes in peltoperlid mitogenomes have a canonical cloverleaf secondary structure except for trnS1, whose dihydrouridine arm simply forms a loop. The control region of the family has several distinct structural characteristics and has the potential to serve as effective phylogenetic markers. Phylogenetic analyses support the monophyly of Perloidea, but the monophyly of Pteronarcyoidea is still not supported. The Peltoperlidae is placed as the earliest branch within the Systellognatha, and the estimated phylogenetic relationship is: Peltoperlidae + {(Styloperlidae + Pteronarcyidae) + [Perlidae + (Chloroperlidae + Perlodidae)]}. Our results provide new insight into the phylogeny of this group.

First Report of Complete Mitochondrial Genome in the Tribes Coomaniellini and Dicercini (Coleoptera: Buprestidae) and Phylogenetic Implications

The complete mitochondrial genomes (mitogenomes) of the tribes Coomaniellini and Dicercini were sequenced and described in this study, including Coomaniella copipes (16,196 bp), Coomaniella dentata (16,179 bp), and Dicerca corrugata (16,276 bp). These complete mitogenomes are very similar in length and encoded 37 typical mitochondrial genes, including 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and 13 protein-coding genes (PCGs). Most of PCGs had typical ATN start codons and terminated with TAR. Among these mitogenomes, Leu2 (L2), Ile (I), Ser2 (S2), and Phe (F) were the four most frequently encoded amino acids. Moreover, phylogenetic analyses were performed based on three kinds of nucleotide matrixes (13 PCGs, 2 rRNAs, and 13 PCGs + 2 rRNAs) among the available sequenced species of the family Buprestidae using Bayesian inference and Maximum-likelihood methods. The results showed that a Chrysochroninae species interspersed in Buprestinae, and Coomaniellini is more closely related to Dicercini than Melanophilini. Moreover, the clade of Buprestidae was well separated from outgroups and the monophyly of Agrilinae is confirmed again. Our whole mitogenome phylogenetic results support that the genus Dicerca can be transferred from Chrysochroinae to Buprestinae; whether Dicercini can be completely transferred remains to be further verified after enriching samples. Our results have produced new complete mitogenomic data, which will provide information for future phylogenetic and taxonomic research.

Characterizing the Complete Mitochondrial Genome of Arma custos and Picromerus lewisi (Hemiptera: Pentatomidae: Asopinae) and Conducting Phylogenetic Analysis

We characterized the mitochondrial genome (mitogenome) and conducted phylogenetic analyses of 48 Hemiptera species by sequencing and analyzing the mitogenome of Arma custos (Fabricius) and Picromerus lewisi (Scott). The complete mitogenomes of the two predators were 16,024 bp and 19,587 bp in length, respectively, and it contained 37 classical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), and a control region. Most PCGs in these predators use ATN as the start codon. This research revealed that the genes of the two natural enemy species have an A + T content of 75.40% and all tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. This is the first study to compare the mitochondrial genetic structure of two predatory insects; the mitochondrial genetic structure of individual predatory insects has been sequenced in previous studies. Here, phylogenetic analysis on the basis of amino acid and nucleotide sequences of 13 mitochondrial PCGs using Bayesian inference and maximum likelihood methods were conducted to generate similar tree topologies, which suggested that the two predators with close genetic relationships belong to Asopinae subfamily. Furthermore, the monophyly of the Pentatomoidea superfamily is well accepted despite limited taxon and species sampling. Finally, their complete mitogenome provided data to establish a predator-prey food web, which is the foundation of effective pest management. Our results also enhanced the database of natural enemy insects.

The complete mitochondrial genome of the jewel beetle, Anthaxia chinensis (Coleoptera: Buprestidae)

The mitogenome of Anthaxia chinensis is a typical circular DNA molecule of 15,881 bp. All protein-coding genes use ATN codons as the start codon except for ATP8 and COI. The stop codon TAA and TAG are assigned to all protein-coding genes. Most of tRNAs have the typical secondary clover-leaf structure except for tRNA-Ala, tRNA-Phe, tRNA-His, and tRNA-Ser (AGN). The length of tRNAs ranged from 63 bp to 70 bp. Genome organization and nucleotide composition of the mitogenome were also noted.

Mitochondrial Genomes from Two Specialized Subfamilies of Reduviidae (Insecta: Hemiptera) Reveal Novel Gene Rearrangements of True Bugs

Reduviidae, a hyper-diverse family, comprise 25 subfamilies with nearly 7000 species and include many natural enemies of crop pests and vectors of human disease. To date, 75 mitochondrial genomes (mitogenomes) of assassin bugs from only 11 subfamilies have been reported. The limited sampling of mitogenome at higher categories hinders a deep understanding of mitogenome evolution and reduviid phylogeny. In this study, the first mitogenomes of Holoptilinae (Ptilocnemus lemur) and Emesinae (Ischnobaenella hainana) were sequenced. Two novel gene orders were detected in the newly sequenced mitogenomes. Combined 421 heteropteran mitogenomes, we identified 21 different gene orders and six gene rearrangement units located in three gene blocks. Comparative analyses of the diversity of gene order for each unit reveal that the tRNA gene cluster trnI-trnQ-trnM is the hotspot of heteropteran gene rearrangement. Furthermore, combined analyses of the gene rearrangement richness of each unit and the whole mitogenome among heteropteran lineages confirm Reduviidae as a 'hot-spot group' of gene rearrangement in Heteroptera. The phylogenetic analyses corroborate the current view of phylogenetic relationships between basal groups of Reduviidae with high support values. Our study provides deeper insights into the evolution of mitochondrial gene arrangement in Heteroptera and the early divergence of reduviids.

The complete mitochondrial genome of the pyrophilous jewel beetle Melanophila acuminata (Coleoptera: Buprestidae)

The complete mitochondrial of genome Melanophila acuminata (DeGeer 1774) is a typical double-stranded circular molecule of 15,853 bp (GenBank accession number: MW287594). All tRNA genes, ranging from 62 to 72 bp, can be folded into typical clover-leaf secondary structure except for tRNA ^Ser(AGN) . The control region is 1,080 bp long with an A+T content of 87.5%. The phylogeny tree is monophyletic among 19 related species. The Melanophila acuminata cluster was more closely related to Chrysochroa fulgidissima. This mitochondrial genome can be used for further analyses of Buprestidae mitochondrial comparative genomics to improve the understanding of diverse coleopteran species.

Comparative mitochondrial genomes of four species of Sinopodisma and phylogenetic implications (Orthoptera, Melanoplinae)

In this study, the whole mitochondrial genomes (mitogenomes) from four species were sequenced. The complete mitochondrial genomes of Sinopodisma pieli, S. houshana, S. qinlingensis, and S. wulingshanensis are 15,857 bp, 15,818 bp, 15,843 bp, and 15,872 bp in size, respectively. The 13 protein-coding genes (PCGs) begin with typical ATN codons, except for COXI in S. qinlingensis, which begins with ACC. The highest A+T content in all the sequenced orthopteran mitogenomes is 76.8% (S. qinlingensis), followed by 76.5% (S. wulingshanensis), 76.4% (S. pieli) and 76.4% (S. houshana) (measured on the major strand). The long polythymine stretches (T-stretch) in the A+T-rich region of the four species are not adjacent to the trnI locus but are inside the stem-loop sequences on the major strand. Moreover, several repeated elements are found in the A+T-rich region of the four species. Phylogenetic analysis based on 53 mitochondrial genomes using Bayesian Inference (BI) and Maximum Likelihood (ML) revealed that Melanoplinae (Podismini) was a monophyletic group; however, the monophyly of Sinopodisma was not supported. These data will provide important information for a better understanding of the phylogenetic relationship of Melanoplinae.

The complete mitochondrial genome of the assassin bug Sycanus croceovittatus (Hemiptera: Reduviidae)

The complete mitochondrial genome (mitogenome) of the assassin bug, Sycanus croceovittatus, was sequenced and analyzed in the present study. This mitogenome spans 15,644 bp in size with a high A + T content (71.7%), containing 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a putative control region. All protein-coding genes are initiated by ATN codons expect ND1 use GTG as start codons and terminated with TAG or TAA, expect COX3 use a single T–– residue as the stop codon. All tRNAs have the typical clover-leaf like structures except for tRNA^Ser(AGN). A phylogenetic analysis of S. croceovittatus and 33 other assassin bugs is also presented using 13 protein-coding genes and 2 rRNA genes. The result supports the monophyly of Harpactorinae and the sister relationship between S. croceovittatus and Agriosphodrus dohrni.

The challenge of Coccidae (Hemiptera: Coccoidea) mitochondrial genomes: The case of Saissetia coffeae with novel truncated tRNAs and gene rearrangements

There have been few reports of complete mitochondrial genomes (mitogenomes) of scale insects, and it has been indicated that complex and novel structures in their mitogenomes may lead to difficulties in sequencing, assembly and annotation. Transfer RNAs (tRNAs) usually possess typical cloverleaf secondary structures, and truncated tRNAs are rarely found in insect mitogenomes. Here, we report a complete Saissetia coffeae mitogenome (15,389 bp) with high A+T content (84.7%) sequenced by next-generation sequencing (NGS) methods. Genes in the mitogenome were annotated, and nine tRNAs were not found using MITOS. Most of the detected tRNAs were significantly truncated without the dihydrouridine (DHU) arm or the TΨC (T) arm. In addition, the 9 "lost" tRNAs containing mismatched base pairs were retrieved based on the tRNA annotation workflow for Coccidae described in our study. The gene arrangement in the Saissetia coffeae mitogenome was significantly different from that in other hemipteran insects. Additionally, Bayesian and maximum likelihood trees based on the mitochondrial genes showed a long branch of the Saissetia lineage, indicating significant nonsynonymous substitutions or high evolutionary rates in the Saissetia lineage. We provide a reference mitogenome for the assembly and annotation of the Coccidae mitogenome and offer insights into the evolution of scale insects.

The complete mitochondrial genome of the soldier fly Ptecticus aurifer

The pig body was put in the wild area in summer for collecting sarcosaphagous insects. After 31-day (9 June 2015 to 9 July 2015), the pig body was to be mummification. Ptecticus aurifer (Walker 1854) was found as sarcosaphagous insect for the first time. The complete mitochondrial genome of P. aurifer (Walker 1854) was sequenced in this study. The complete mitochondrial genome is a typical double-stranded circular molecule of 15,775 bp (GenBank accession number: MN604259) containing 37 typical animal mitochondrial gene and an A + T-rich region. 11 of the 22 tRNAs, ranging from 63 to 72 bp, can be folded into classic clover-leaf secondary structure except for tRNA^Ser(AGN), in which the dihydrouridine (DHU) arm did not form a stable stem-loop structure. The control region is 954 bp long with an A + T content of 90.7%.

Mitogenome evolution in ladybirds: Potential association with dietary adaptation

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein-coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.

Complete Mitogenome of a Leaf-Mining Buprestid Beetle, Trachys auricollis, and Its Phylogenetic Implications

A complete mitogenome of Trachys auricollis is reported, and a mitogenome-based phylogenetic tree of Elateriformia with all protein-coding genes (PCGs), rRNAs, and tRNAs is presented for the first time. The complete mitochondrial genome of T. auricollis is 16,429 bp in size and contains 13 PCGs, two rRNA genes, 22 tRNA genes, and an A + T-rich region. The A + T content of the entire genome is approximately 71.1%, and the AT skew and GC skew are 0.10 and -0.20, respectively. According to the the nonsynonymous substitution rate to synonymous substitution rates (Ka/Ks) of all PCGs, the highest and lowest evolutionary rates were observed for atp8 and cox1, respectively, which is a common finding among animals. The start codons of all PCGs are the typical ATN. Ten PCGs have complete stop codons, but three have incomplete stop codons with T or TA. As calculated based on the relative synonymous codon usage (RSCU) values, UUA(L) is the codon with the highest frequency. Except for trnS1, all 22 tRNA genes exhibit typical cloverleaf structures. The A + T-rich region of T. auricollis is located between rrnS and the trnI-trnG-trnM gene cluster, with six 72-bp tandem repeats. Both maximum likelihood (ML) and Bayesian (BI) trees suggest that Buprestoidea is close to Byrrhoidea and that Buprestoidea and Byrrhoidea are sister groups of Elateroidea, but the position of Psephenidae is undetermined. The inclusion of tRNAs might help to resolve the phylogeny of Coleoptera.

The complete mitochondrial genome of the jewel beetle Trachys variolaris (Coleoptera: Buprestidae)

The complete mitochondrial genome (mitogenome) of the jewel beetle, Trachys variolaris, was sequenced and described. The mitogenome of T. variolaris is a typical circular DNA molecule with 16,771 bp in length. It contains the typical 37 mitochondrial genes (13 protein-coding genes, two rRNA genes, and 22 tRNA genes) and a long non-coding region called the control region. Twelve protein-coding genes initiate with ATN codons except ND1 uses TTG. Most of the protein-coding genes use TAA or TAG as the stop codon, but COII, COIII, ND4, and ND5 terminate with a single T-- or TA-. The length of the 22 tRNAs ranges from 61 to 70 bp and they all have the clover-leaf structure except for tRNASer(AGN) . The control region is 2155 bp long with the A + T content of 69.4%. The result of our phylogenetic analysis showed that Buprestoidea is monophyletic, and it is the sister group to (Byrrhoidea + Elateroidea).

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

The complete mitochondrial genome of the jewel beetle Coraebus cavifrons (Coleoptera: Buprestidae)

The complete mitochondrial genome (mitogenome) of the jewel beetle, Coraebus cavifrons, was sequenced and described in the present study. The mitogenome of Coraebus cavifrons is a typical circular DNA molecule of 15,686 bp. All of the 37 mitochondrial genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes) were annotated and a putative control region with 1135 bp in length was found between tRNAIle and srRNA. Protein-coding genes all initiate with ATN codons except for COII uses GTG. Most of the protein-coding genes use TAA or TAG as the stop codon, but COII, COIII, and ND5 terminate with a single T. All tRNAs have the clover-leaf structure except for tRNASer(AGN) and the length of them range from 60 to 71 bp. Genome organization and nucleotide composition of the mitogenome were also noted. Our phylogenetic analysis of Elateriformia supported the monophyly of Buprestoidea and the sister relationship between Buprestoidea and (Byrrhoidea + Elateroidea).

The complete mitochondrial genome of the assassin bug Reduvius gregoryi (Hemiptera: Reduviidae)

The complete mitochondrial genome (mitogenome) of the assassin bug, Reduvius gregoryi, was determined. The sequenced mitogenome is a typical circular DNA molecule of 16,477 bp, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a putative control region. Protein-coding genes all initiate with ATN codons and terminate with TAA codons except for ATP6, COI, COIII, ND4, and ND5 use a single T residue as the termination codon. All tRNAs have the clover-leaf structure except for the tRNA^Ser(AGN) and the length of them range from 61 to 70 bp. The control region is 1731 bp long with an A + T content of 72.3%. Our phylogenetic analysis supported the polyphyly of Reduviinae and the sister relationship between Reduvius gregoryi and Reduvius tenebrosus.

Characteristics of the complete mitochondrial genome of Suhpalacsa longialata (Neuroptera, Ascalaphidae) and its phylogenetic implications

The owlflies (Family Ascalaphidae) belong to the Neuroptera but are often mistaken as dragonflies because of morphological characters. To date, only three mitochondrial genomes of Ascalaphidae, namely Libelloides macaronius; Ascaloptynx appendiculatus; Ascalohybris subjacens, are published in GenBank, meaning that they are greatly under-represented in comparison with the 430 described species reported in this family. In this study, we sequenced and described the complete mitochondrial genome of Suhpalacsa longialata (Neuroptera, Ascalaphidae). The total length of the S. longialata mitogenome was 15,911 bp, which is the longest known to date among the available family members of Ascalaphidae. However, the size of each gene was similar to the other three Ascalaphidae species. The S. longialata mitogenome included a transposition of tRNACys and tRNATrp genes and formed an unusual gene arrangement tRNACys-tRNATrp-tRNATyr (CWY). It is likely that the transposition occurred by a duplication of both genes followed by random loss of partial duplicated genes. The nucleotide composition of the S. longialata mitogenome was as follows: A = 41.0%, T = 33.8%, C = 15.5%, G = 9.7%. Both Bayesian inference and ML analyses strongly supported S. longialata as a sister clade to (Ascalohybris subjacens + L. macaronius), and indicated that Ascalaphidae is not monophyletic.

Complete mitochondrial genome of the ambush bug, Amblythyreus gestroi (Hemiptera: Reduviidae)

The mitochondrial genome (mitogenome) of Amblythyreus gestroi is described in the present paper. The complete mitogenome is a 15,228 bp circular DNA molecule, containing 13 protein-coding genes, 22 tRNA genes, two rRNA genes and a control region. Genome organization, nucleotide composition and codon usage of the mitogenome are noted, secondary structures of all tRNAs are predicted. The monophyly of Phymatinae is highly supported by the phylogenetic tree. Phylogenetic implications of the A. gestroi mitogenome is briefly discussed.

The complete mitochondrial genome of Orancistrocerusaterrimusaterrimus and comparative analysis in the family Vespidae (Hymenoptera, Vespidae, Eumeninae)

To date, only one mitochondrial genome (mitogenome) in the Eumeninae has been reported in the world and this is the first report in China. The mitogenome of O.a.aterrimus is 17 972 bp long, and contains 38 genes, including 13 protein coding genes (PCGs), 23 tRNA genes, two rRNA genes, a long non-coding region (NCR), and a control region (CR). The mitogenome has 79.43% A + T content, its 13 PCGs use ATN as the initiation codon except for cox1 using TTG, and nine genes used complete translation termination TAA and four genes have incomplete stop codon T (cox2, cox3, nad4, and cytb). Twenty-two of 23 tRNAs can form the typical cloverleaf secondary structure except for trnS1. The CR is 1 078 bp long with 84.69% A+T content, comprising 28 bp tandem repeat sequences and 13 bp T-strech. There are two gene rearrangements which are an extra trnM2 located between trnQ and nad2 and the trnL2 in the upstream of nad1. Within all rearrangements of these mitogenomes reported in the family Vespidae, the translocation between trnS1 and trnE genes only appears in Vespinae, and the translocation of trnY in Polistinae and Vespinae. The absent codons of 13 PCGs in Polistinae are more than those both in Vespinae and Eumeninae in the family Vespidae. The study reports the complete mitogenome of O.a.aterrimus, compares the characteristics and construct phylogenetic relationships of the mitogenomes in the family Vespidae.

The mitochondrial genome of the ambush bug Carcinochelis bannaensis (Hemiptera: Reduviidae)

The complete mitochondrial genome (mitogenome) of the ambush bug, Carcinochelis bannaensis, was determined in this study. The sequenced mitogenome is a typical circular DNA molecule of 15,335 bp, containing 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a putative control region. All protein-coding genes initiate with ATN codons and terminate with TAA codons except for COII,COIII and ND5 use a single T residue as the stop codon. All tRNAs have the clover-leaf structure except for the tRNA^Ser(AGN) and the length of them range from 61 to 71 bp. The control region is 797 bp long with an A + T content of 66.3%. The phylogenetic analysis result supports the monophyly of Phymatinae.

The Mitochondrial Genome of the Guanaco Louse, Microthoracius praelongiceps: Insights into the Ancestral Mitochondrial Karyotype of Sucking Lice (Anoplura, Insecta)

Fragmented mitochondrial (mt) genomes have been reported in 11 species of sucking lice (suborder Anoplura) that infest humans, chimpanzees, pigs, horses, and rodents. There is substantial variation among these lice in mt karyotype: the number of minichromosomes of a species ranges from 9 to 20; the number of genes in a minichromosome ranges from 1 to 8; gene arrangement in a minichromosome differs between species, even in the same genus. We sequenced the mt genome of the guanaco louse, Microthoracius praelongiceps, to help establish the ancestral mt karyotype for sucking lice and understand how fragmented mt genomes evolved. The guanaco louse has 12 mt minichromosomes; each minichromosome has 2-5 genes and a non-coding region. The guanaco louse shares many features with rodent lice in mt karyotype, more than with other sucking lice. The guanaco louse, however, is more closely related phylogenetically to human lice, chimpanzee lice, pig lice, and horse lice than to rodent lice. By parsimony analysis of shared features in mt karyotype, we infer that the most recent common ancestor of sucking lice, which lived ∼75 Ma, had 11 minichromosomes; each minichromosome had 1-6 genes and a non-coding region. As sucking lice diverged, split of mt minichromosomes occurred many times in the lineages leading to the lice of humans, chimpanzees, and rodents whereas merger of minichromosomes occurred in the lineage leading to the lice of pigs and horses. Together, splits and mergers of minichromosomes created a very complex and dynamic mt genome organization in the sucking lice.

Hemiptera Mitochondrial Control Region: New Sights into the Structural Organization, Phylogenetic Utility, and Roles of Tandem Repetitions of the Noncoding Segment

As a major noncoding fragment, the control region (CR) of mtDNA is responsible for the initiation of mitogenome transcription and replication. Several structural features of CR sequences have been reported in many insects. However, comprehensive analyses on the structural organization and phylogenetic utility, as well as the role of tandem replications (TRs) on length variation, high A+T content, and shift of base skew of CR sequences are poorly investigated in hemipteran insects. In this study, we conducted a series of comparative analyses, using 116 samples covering all 11 infraorders of the five currently recognized monophyletic groups in the Hemiptera. Several structural elements (mononucleotide stretches containing conserved sequence blocks (CSBs), TRs, and GA-rich region) were identified in the mitochondrial control region in hemipteran insects, without showing a consistent location. The presence and absence of certain specific structural elements in CR sequences show the various structural organizations of that segment among the five monophyletic groups, which indicates the diversification of the control region’s structural organization in Hemiptera. Among the many groups within Hemiptera, eight monophyletic groups and three consistent phylogenetic trees were recovered, using CSBs datasets by maximum likelihood and Bayesian methods, which suggests the possible utility of CR sequences for phylogenetic reconstruction in certain groups of Hemiptera. Statistical analyses showed that TRs may contribute to the length variation, high AT content, and the shift of base skewing of CR sequences toward high AT content in the Hemiptera. Our findings enrich the knowledge of structural organization, phylogenetic utility, and roles of tandem replication of hemipteran CR, and provide a possible framework for mitochondrial control region analyses in hemimetabolous insects.

A Mitochondrial Genome of Rhyparochromidae (Hemiptera: Heteroptera) and a Comparative Analysis of Related Mitochondrial Genomes

The Rhyparochromidae, the largest family of Lygaeoidea, encompasses more than 1,850 described species, but no mitochondrial genome has been sequenced to date. Here we describe the first mitochondrial genome for Rhyparochromidae: a complete mitochondrial genome of Panaorus albomaculatus (Scott, 1874). This mitochondrial genome is comprised of 16,345 bp, and contains the expected 37 genes and control region. The majority of the control region is made up of a large tandem-repeat region, which has a novel pattern not previously observed in other insects. The tandem-repeats region of P. albomaculatus consists of 53 tandem duplications (including one partial repeat), which is the largest number of tandem repeats among all the known insect mitochondrial genomes. Slipped-strand mispairing during replication is likely to have generated this novel pattern of tandem repeats. Comparative analysis of tRNA gene families in sequenced Pentatomomorpha and Lygaeoidea species shows that the pattern of nucleotide conservation is markedly higher on the J-strand. Phylogenetic reconstruction based on mitochondrial genomes suggests that Rhyparochromidae is not the sister group to all the remaining Lygaeoidea, and supports the monophyly of Lygaeoidea.

Duplication and Remolding of tRNA Genes in the Mitochondrial Genome of Reduvius tenebrosus (Hemiptera: Reduviidae)

Most assassin bugs are predators that act as important natural enemies of insect pests. Mitochondrial (mt) genomes of these insects are double-strand circular DNAs that encode 37 genes. In the present study, we explore the duplication and rearrangement of tRNA genes in the mt genome of Reduvius tenebrosus, the first mt genome from the subfamily Reduviinae. The gene order rearranges from CR (control region)-trnI-trnQ-trnM-ND2 to CR-trnQ-trnI2-trnI1-trnM-ND2. We identified 23 tRNA genes, including 22 tRNAs commonly found in insects and an additional trnI (trnI2), which has high sequence similarity to trnM. We found several pseudo genes, such as pseudo-trnI, pseudo-CR, and pseudo-ND2, in the hotspot region of gene rearrangement (between the control region and ND2). These features provided evidence that this novel gene order could be explained by the tandem duplication/random loss (TDRL) model. The tRNA duplication/anticodon mutation mechanism further explains the presence of trnI2, which is remolded from a duplicated trnM in the TDRL process (through an anticodon mutation of CAT to GAT). Our study also raises new questions as to whether the two events proceed simultaneously and if the remolded tRNA gene is fully functional. Significantly, the duplicated tRNA gene in the mitochondrial genome has evolved independently at least two times within assassin bugs.

Rearrangement of mitochondrial tRNA genes in flat bugs (Hemiptera: Aradidae)

The typical insect mitochondrial (mt) genome organization, which contains a single chromosome with 37 genes, was found in the infraorder Pentatomomorpha (suborder Heteroptera). The arrangement of mt genes in these true bugs is usually the same as the ancestral mt gene arrangement of insects. Rearrangement of transfer RNA (tRNA) genes, however, has been found in two subfamilies of flat bugs (Mezirinae and Calisiinae, family Aradidae). In this study, we sequenced the complete mt genomes of four species from three other subfamilies (Aradinae, Carventinae and Aneurinae). We found tRNA gene rearrangement in all of these four species. All of the rearranged tRNA genes are located between the mitochondrial control region and cox1, indicating this region as a hotspot for gene rearrangement in flat bugs; the rearrangement is likely caused by events of tandem duplication and random deletion of genes. Furthermore, our phylogenetic and dating analyses indicated that the swap of positions between trnQ and trnI occurred ~162 million years ago (MYA) in the most recent common ancestor of the five subfamilies of flat bugs investigated to date, whereas the swap of positions between trnC and trnW occurred later in the lineage leading to Calisiinae, and the translocation of trnC and trnY occurred later than 134 MYA in the lineage leading to Aradinae.

The Complete Mitochondrial Genome of Mindarus keteleerifoliae (Insecta: Hemiptera: Aphididae) and Comparison with Other Aphididae Insects

The mitogenome of Mindarus keteleerifoliae Zhang (Hemiptera: Aphididae) is a 15,199 bp circular molecule. The gene order and orientation of M. keteleerifoliae is similarly arranged to that of the ancestral insect of other aphid mitogenomes, and, a tRNA isomerism event maybe identified in the mitogenome of M. keteleerifoliae. The tRNA-Trp gene is coded in the J-strand and the same sequence in the N-strand codes for the tRNA-Ser gene. A similar phenomenon was also found in the mitogenome of Eriosoma lanigerum. However, whether tRNA isomers in aphids exist requires further study. Phylogenetic analyses, using all available protein-coding genes, support Mindarinae as the basal position of Aphididae. Two tribes of Aphidinae were recovered with high statistical significance. Characteristics of the M. keteleerifoliae mitogenome revealed distinct mitogenome structures and provided abundant phylogenetic signals, thus advancing our understanding of insect mitogenomic architecture and evolution. But, because only eight complete aphid mitogenomes, including M. keteleerifoliae, were published, future studies with larger taxon sampling sizes are necessary.

Fragmented mitochondrial genomes in two suborders of parasitic lice of eutherian mammals (Anoplura and Rhynchophthirina, Insecta)

Parasitic lice (order Phthiraptera) infest birds and mammals. The typical animal mitochondrial (mt) genome organization, which consists of a single chromosome with 37 genes, was found in chewing lice in the suborders Amblycera and Ischnocera. The sucking lice (suborder Anoplura) known, however, have fragmented mt genomes with 9–20 minichromosomes. We sequenced the mt genome of the elephant louse, Haematomyzus elephantis – the first species of chewing lice investigated from the suborder Rhynchophthirina. We identified 33 mt genes in the elephant louse, which were on 10 minichromosomes. Each minichromosome is 3.5–4.2 kb in size and has 2–6 genes. Phylogenetic analyses of mt genome sequences confirm that the elephant louse is more closely related to sucking lice than to the chewing lice in the Amblycera and Ischnocera. Our results indicate that mt genome fragmentation is shared by the suborders Anoplura and Rhynchophthirina. Nine of the 10 mt minichromosomes of the elephant louse differ from those of the sucking lice (Anoplura) known in gene content and gene arrangement, indicating that distinct mt karyotypes have evolved in Anoplura and Rhynchophthirina since they diverged ~92 million years ago.

First complete mitochondrial genome from the brown lacewings (Neuroptera: Hemerobiidae)

The first complete mitochondrial genome of the family Hemerobiidae (Neuronema laminatum Tjeder, 1936) is sequenced in this study. The complete mitochondrial genome is a typical double-stranded circular molecule of 17,164 bp (GenBank accession number: KR078257) containing 37 typical animal mitochondrial gene and an A + T-rich region. The gene order is identical to that of the putative ancestral arrangement of insects and other lacewings. Thirteen protein-coding genes (PCGs) possessed common triplet initiation codons ATN and mostly terminate with TAN codons except for ND5 with a single T residue adjacent to a downstream tRNA gene. All the 22 tRNAs, ranging from 63 to 72 bp, can be folded into classic clover-leaf secondary structure except for tRNA(Ser(AGN)), in which the dihydrouridine (DHU) arm did not form a stable stem-loop structure. The control region is 2131 bp long with an A + T content of 87.4%. In the sampled families of Neuroptera, Osmylidae + the remaining families, Hemerobiidae + Mantispidae, Polystoechotidae + Rapismatidae, are recovered in phylogenetic analyses with high supports.

The mitochondrial genome of the butterfly Polyura schreiber (Lepidoptera: Nymphalidae)

The nearly complete mitochondrial genome (mitogenome) of the butterfly, Polyura schreiber, was determined. The sequenced mitogenome is a typical circular DNA molecule of 15 320 bp, containing 13 protein-coding genes, two rRNA genes, 21 tRNA genes, and a putative control region. tRNA(Phe) was failed to sequence, which was presumed to be located between tRNA(Glu) and ND5. Protein-coding genes all initiate with ATN codons and terminate with TAA codons except for COII and ND5 use a single T residue as the termination codon. All tRNAs have the clover-leaf structure except for the tRNA(Ser(AGN)) and the length of them range from 65 to 71 bp. The control region is 412 bp long with an A + T content of 90.5%. Our phylogenetic analysis recovered the sister-group relationship between Charaxinae and Satyrinae.

Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications

BACKGROUND

Insect mitochondrial genomes (mitogenomes) are the most extensively used genetic marker for evolutionary and population genetics studies of insects. The Pentatomoidea superfamily is economically important and the largest superfamily within Pentatomomorpha with over 7,000 species. To better understand the diversity and evolution of pentatomoid species, we sequenced and annotated the mitogenomes of Eurydema gebleri and Rubiconia intermedia, and present the first comparative analysis of the 11 pentatomoid mitogenomes that have been sequenced to date.

RESULTS

We obtained the complete mitogenome of Eurydema gebleri (16,005 bp) and a nearly complete mitogenome of Rubiconia intermedia (14,967 bp). Our results show that gene content, gene arrangement, base composition, codon usage, and mitochondrial transcription termination factor sequences are highly conserved in pentatomoid species, especially for species in the same family. Evolutionary rate analyses of protein-coding genes reveal that the highest and lowest rates are found in atp8 and cox1 and distinctive evolutionary patterns are significantly correlated with the G + C content of genes. We inferred the secondary structures for two rRNA genes for eleven pentatomoid species, and identify some conserved motifs of RNA structures in Pentatomidea. All tRNA genes in pentatomoid mitogenomes have a canonical cloverleaf secondary structure, except for two tRNAs (trnS1 and trnV) which appear to lack the dihydrouridine arm. Regions that are A + T-rich have several distinct characteristics (e.g. size variation and abundant tandem repeats), and have potential as species or population level molecular markers. Phylogenetic analyses based on mitogenomic data strongly support the monophyly of Pentatomoidea, and the estimated phylogenetic relationships are: (Urostylididae + (Plataspidae + (Pentatomidae + (Cydnidae + (Dinidoridae + Tessaratomidae))))).

CONCLUSIONS

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the superfamily Pentatomoidea. Mitogenomes can be effectively used to resolve phylogenetic relationships of pentatomomorphan insects at various taxonomic levels. Sequencing more mitogenomes at various taxonomic levels, particularly from closely related species, will improve the annotation accuracy of mitochondrial genes, as well as greatly enhance our understanding of mitogenomic evolution and phylogenetic relationships in pentatomoids.

The Complete Mitochondrial Genome of Corizus tetraspilus (Hemiptera: Rhopalidae) and Phylogenetic Analysis of Pentatomomorpha

Insect mitochondrial genome (mitogenome) are the most extensively used genetic information for molecular evolution, phylogenetics and population genetics. Pentatomomorpha (>14,000 species) is the second largest infraorder of Heteroptera and of great economic importance. To better understand the diversity and phylogeny within Pentatomomorpha, we sequenced and annotated the complete mitogenome of Corizus tetraspilus (Hemiptera: Rhopalidae), an important pest of alfalfa in China. We analyzed the main features of the C. tetraspilus mitogenome, and provided a comparative analysis with four other Coreoidea species. Our results reveal that gene content, gene arrangement, nucleotide composition, codon usage, rRNA structures and sequences of mitochondrial transcription termination factor are conserved in Coreoidea. Comparative analysis shows that different protein-coding genes have been subject to different evolutionary rates correlated with the G+C content. All the transfer RNA genes found in Coreoidea have the typical clover leaf secondary structure, except for trnS1 (AGN) which lacks the dihydrouridine (DHU) arm and possesses a unusual anticodon stem (9 bp vs. the normal 5 bp). The control regions (CRs) among Coreoidea are highly variable in size, of which the CR of C. tetraspilus is the smallest (440 bp), making the C. tetraspilus mitogenome the smallest (14,989 bp) within all completely sequenced Coreoidea mitogenomes. No conserved motifs are found in the CRs of Coreoidea. In addition, the A+T content (60.68%) of the CR of C. tetraspilus is much lower than that of the entire mitogenome (74.88%), and is lowest among Coreoidea. Phylogenetic analyses based on mitogenomic data support the monophyly of each superfamily within Pentatomomorpha, and recognize a phylogenetic relationship of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))).

Hemipteran mitochondrial genomes: features, structures and implications for phylogeny

The study of Hemipteran mitochondrial genomes (mitogenomes) began with the Chagas disease vector, Triatoma dimidiata, in 2001. At present, 90 complete Hemipteran mitogenomes have been sequenced and annotated. This review examines the history of Hemipteran mitogenomes research and summarizes the main features of them including genome organization, nucleotide composition, protein-coding genes, tRNAs and rRNAs, and non-coding regions. Special attention is given to the comparative analysis of repeat regions. Gene rearrangements are an additional data type for a few families, and most mitogenomes are arranged in the same order to the proposed ancestral insect. We also discuss and provide insights on the phylogenetic analyses of a variety of taxonomic levels. This review is expected to further expand our understanding of research in this field and serve as a valuable reference resource.

Higher-level phylogeny of paraneopteran insects inferred from mitochondrial genome sequences

Mitochondrial (mt) genome data have been proven to be informative for animal phylogenetic studies but may also suffer from systematic errors, due to the effects of accelerated substitution rate and compositional heterogeneity. We analyzed the mt genomes of 25 insect species from the four paraneopteran orders, aiming to better understand how accelerated substitution rate and compositional heterogeneity affect the inferences of the higher-level phylogeny of this diverse group of hemimetabolous insects. We found substantial heterogeneity in base composition and contrasting rates in nucleotide substitution among these paraneopteran insects, which complicate the inference of higher-level phylogeny. The phylogenies inferred with concatenated sequences of mt genes using maximum likelihood and Bayesian methods and homogeneous models failed to recover Psocodea and Hemiptera as monophyletic groups but grouped, instead, the taxa that had accelerated substitution rates together, including Sternorrhyncha (a suborder of Hemiptera), Thysanoptera, Phthiraptera and Liposcelididae (a family of Psocoptera). Bayesian inference with nucleotide sequences and heterogeneous models (CAT and CAT + GTR), however, recovered Psocodea, Thysanoptera and Hemiptera each as a monophyletic group. Within Psocodea, Liposcelididae is more closely related to Phthiraptera than to other species of Psocoptera. Furthermore, Thysanoptera was recovered as the sister group to Hemiptera.

Comparative mitogenomics of the assassin bug genus Peirates (Hemiptera: Reduviidae: Peiratinae) reveal conserved mitochondrial genome organization of P. atromaculatus, P. fulvescens and P. turpis

In this study, we sequenced four new mitochondrial genomes and presented comparative mitogenomic analyses of five species in the genus Peirates (Hemiptera: Reduviidae). Mitochondrial genomes of these five assassin bugs had a typical set of 37 genes and retained the ancestral gene arrangement of insects. The A+T content, AT- and GC-skews were similar to the common base composition biases of insect mtDNA. Genomic size ranges from 15,702 bp to 16,314 bp and most of the size variation was due to length and copy number of the repeat unit in the putative control region. All of the control region sequences included large tandem repeats present in two or more copies. Our result revealed similarity in mitochondrial genomes of P. atromaculatus, P. fulvescens and P. turpis, as well as the highly conserved genomic-level characteristics of these three species, e.g., the same start and stop codons of protein-coding genes, conserved secondary structure of tRNAs, identical location and length of non-coding and overlapping regions, and conservation of structural elements and tandem repeat unit in control region. Phylogenetic analyses also supported a close relationship between P. atromaculatus, P. fulvescens and P. turpis, which might be recently diverged species. The present study indicates that mitochondrial genome has important implications on phylogenetics, population genetics and speciation in the genus Peirates.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae)

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.

Characterization of the complete mitochondrial genomes of Nematodirus oiratianus and Nematodirus spathiger of small ruminants

BACKGROUND

Nematodirus spp. are among the most common nematodes of ruminants worldwide. N. oiratianus and N. spathiger are distributed worldwide as highly prevalent gastrointestinal nematodes, which cause emerging health problems and economic losses. Accurate identification of Nematodirus species is essential to develop effective control strategies for Nematodirus infection in ruminants. Mitochondrial DNA (mtDNA) could provide powerful genetic markers for identifying these closely related species and resolving phylogenetic relationships at different taxonomic levels.

METHODS

In the present study, the complete mitochondrial (mt) genomes of N. oiratianus and N. spathiger from small ruminants in China were obtained using Long-range PCR and sequencing.

RESULTS

The complete mt genomes of N. oiratianus and N. spathiger were 13,765 bp and 13,519 bp in length, respectively. Both mt genomes were circular and consisted of 36 genes, including 12 genes encoding proteins, 2 genes encoding rRNA, and 22 genes encoding tRNA. Phylogenetic analyses based on the concatenated amino acid sequence data of all 12 protein-coding genes by Bayesian inference (BI), Maximum likelihood (ML) and Maximum parsimony (MP) showed that the two Nematodirus species (Molineidae) were closely related to Dictyocaulidae.

CONCLUSIONS

The availability of the complete mtDNA sequences of N. oiratianus and N. spathiger not only provides new mtDNA sources for a better understanding of nematode mt genomics and phylogeny, but also provides novel and useful genetic markers for studying diagnosis, population genetics and molecular epidemiology of Nematodirus spp. in small ruminants.

Comparative mitogenomics of plant bugs (Hemiptera: Miridae): identifying the AGG codon reassignments between serine and lysine

Insect mitochondrial genomes are very important to understand the molecular evolution as well as for phylogenetic and phylogeographic studies of the insects. The Miridae are the largest family of Heteroptera encompassing more than 11,000 described species and of great economic importance. For better understanding the diversity and the evolution of plant bugs, we sequence five new mitochondrial genomes and present the first comparative analysis of nine mitochondrial genomes of mirids available to date. Our result showed that gene content, gene arrangement, base composition and sequences of mitochondrial transcription termination factor were conserved in plant bugs. Intra-genus species shared more conserved genomic characteristics, such as nucleotide and amino acid composition of protein-coding genes, secondary structure and anticodon mutations of tRNAs, and non-coding sequences. Control region possessed several distinct characteristics, including: variable size, abundant tandem repetitions, and intra-genus conservation; and was useful in evolutionary and population genetic studies. The AGG codon reassignments were investigated between serine and lysine in the genera Adelphocoris and other cimicomorphans. Our analysis revealed correlated evolution between reassignments of the AGG codon and specific point mutations at the antidocons of tRNALys and tRNASer(AGN). Phylogenetic analysis indicated that mitochondrial genome sequences were useful in resolving family level relationship of Cimicomorpha. Comparative evolutionary analysis of plant bug mitochondrial genomes allowed the identification of previously neglected coding genes or non-coding regions as potential molecular markers. The finding of the AGG codon reassignments between serine and lysine indicated the parallel evolution of the genetic code in Hemiptera mitochondrial genomes.

Mitogenomic sequences effectively recover relationships within brush-footed butterflies (Lepidoptera: Nymphalidae)

BACKGROUND

Mitogenomic phylogenies have revealed well-supported relationships for many eukaryote groups. In the order Lepidoptera, 113 species mitogenomes had been sequenced (May 14, 2014). However, these data are restricted to ten of the forty-three recognised superfamilies, while it has been challenging to recover large numbers of mitogenomes due to the time and cost required for primer design and sequencing. Nuclear rather than mitochondrial genes have been preferred to reconstruct deep-level lepidopteran phylogenies, without seriously evaluating the potential of entire mitogenomes. Next-generation sequencing methods remove these limitations by providing efficiently massive amounts of sequence data. In the present study, we simultaneously obtained a large number of nymphalid butterfly mitogenomes to evaluate the utility of mitogenomic phylogenies by comparing reconstructions to the now quite well established phylogeny of Nymphalidae.

RESULTS

We newly obtained 30 nymphalid mitogenomes via pyrosequencing on the Roche 454 GS Junior system, and combined these sequences with publicly accessible data to provide a 70-taxa dataset covering 37 genes for a 15,495 bp alignment. Polymorphic sites were not homogeneously distributed across the gene. Two gene regions, nad6 and 3' end of nad5, were most variable, whereas the cox1 and 5' ends of rrnL were most conserved. Phylogenetic relationships inferred by two likelihood methods were congruent and strongly supported (>0.95 posterior probability; ML bootstrap >85%), across the majority of nodes for multiple partitioning strategies and substitution models. Bayes factor results showed that the most highly partitioned dataset is the preferred strategy among different partitioning schemes. The most striking phylogenetic findings were that the subfamily Danainae not Libytheinae was sister of the remaining brush-footed butterflies and that, within Limenitidini, the genus Athyma was clearly polyphyletic. None of the single-gene phylogenies recovered the highly supported topologies generated on the basis of the whole mitogenomic data.

CONCLUSIONS

Thirty mitogenomes were assembled with 89% completeness from the contigs of pyrosequencing-derived reads. Entire mitogenomes or higher-quality sequences could be obtained by increasing pyrosequencing read coverage or by additional Sanger sequencing. Our mitogenomic phylogenies provide robust nodal support at a range of levels, demonstrating that mitogenomes are both accurate and efficient molecular markers for inferring butterfly phylogeny.

Complete mitochondrial genome sequence of the plant bug Adelphocoris fasciaticollis (Hemiptera: Heteroptera: Miridae)

The mitochondrial genome of Adelphocoris fasciaticollis Reuter, an important pest plant bug, is a typical circular DNA molecule of 15,434 bp with 37 genes and 77.4% A + T content. Twenty-three genes are located on the J-strand, the remaining being oriented on the N-strand. Gene order is identical to that of the typical arrangement of insects and other plant bugs. This genome is highly economized with 54 overlapped nucleotides between neighboring genes in 12 locations. All protein-coding genes initiate with ATN codons. All of the 22 tRNAs, ranging from 62 to 71 bp, have the clover-leaf structure. Two tRNAs, trnK and trnS1, use the non-standard anticodons TTT and TCT. The sizes of the large and small ribosomal RNA genes are 1229 and 793 bp, respectively. The control region is 922 bp in length with 78.2% A + T content and includes two 30-bp repeat units and two 35-bp repeat units.

The complete mitochondrial genome of an assassin bug Peirates arcuatus (Hemiptera: Reduviidae)

The complete mitochondrial genome of Peirates arcuatus (Stål) is a typical double-stranded circular molecule 16,176 bp long with 37 genes usually present in animal mitochondrial genomes and a control region. Gene order is identical to that of the putative ancestral arrangement of insects and other assassin bugs. Thirteen protein-coding genes initiate with ATN codons and mostly terminate with TAN codons except for COII and COIII use a single T residue as the termination codon. All tRNA sequences can be folded into classic clover-leaf secondary structure except that the dihydrouridine (DHU) arm of tRNA(Ser(AGN)) forms a simple loop, and their length range from 63 to 69 bp. The control region is 1552 bp long and includes two complete 118-bp tandem repeats and one partial copy of anterior repeat unit.

Complete mitochondrial genome of the assassin bug Oncocephalus breviscutum (Hemiptera: Reduviidae)

The mitochondrial genome of an assassin bug, Oncocephalus breviscutum Reuter, is a typical circular DNA molecule of 15,984 bp with 37 genes and a large control region. The gene order is identical to that of the putative ancestral arrangement of insects. Twelve protein-coding genes start with ATN codon and ND4L uses GTG. All of the 22 tRNAs, ranging from 61 to 70 bp, have the clover-leaf structure except for the dihydrouridine (DHU) arm of trnS2 forms a simple loop. The control region is 1345 bp in length and includes six tandem repeats of three 31-nt and three 145-nt units.

The complete mitochondrial genome sequence of Sogatella furcifera (Horváth) and a comparative mitogenomic analysis of three predominant rice planthoppers

The white-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the most destructive pests of rice crops in many Asian countries. Using long-PCR and shotgun library methods, we sequenced the entire mitochondrial genomes (mt-genomes) of two WBPH individuals. Total lengths of the mt-genome of the two WBPH individuals were 16,612 bp and 16,654 bp with an identical AT content of 76.19%. Among the 13 protein coding genes (PCGs), only nad5 used an atypical initiation codon GTG. Most of the tRNA genes had the typical cloverleaf secondary structure except that the dihydrouridine (DHU) arms in two trnS genes and the TΨC arm of trnG gene did not form a stable stem-loop structure. Similar to the brown planthopper (BPH), Nilaparvata lugens (Stål), and the small brown planthopper (SBPH), Laodelphax striatellus (Fallén), some extraordinary features were observed in the WBPH mt-genome. These include similar gene rearrangement pattern, unusually short length of the atp8 gene and variable numbers of tandem repeat (VNTR) structure in control region. Interestingly, the same tandem repeat unit with stable secondary structure appeared in two different planthoppers, WBPH and SBPH, which belong to two different genera of the Delphacidae. This peculiar feature provides a direct evidence for the close relationship between the two planthoppers and updates our understanding of the evolutionary characteristics of mitochondrial control region. Comparison with two other predominant rice planthoppers (BPH and SBPH) revealed that different PCGs of mitochondria exhibit different evolutionary patterns.