Mitochondrial genomes of two Barklice, Psococerastis albimaculata and Longivalvus hyalospilus (Psocoptera: Psocomorpha): contrasting rates in mitochondrial gene rearrangement between major lineages of Psocodea

Comparative analyses of mitogenomes in the social bees with insights into evolution of long inverted repeats in the Meliponini

The insect mitogenome is typically a compact circular molecule with highly conserved gene contents. Nonetheless, mitogenome structural variations have been reported in specific taxa, and gene rearrangements, usually the tRNAs, occur in different lineages. Because synapomorphies of mitogenome organizations can provide information for phylogenetic inferences, comparative analyses of mitogenomes have been given increasing attention. However, most studies use a very few species to represent the whole genus, tribe, family, or even order, overlooking potential variations at lower taxonomic levels, which might lead to some incorrect inferences. To provide new insights into mitogenome organizations and their implications for phylogenetic inference, this study conducted comparative analyses for mitogenomes of three social bee tribes (Meliponini, Bombini, and Apini) based on the phylogenetic framework with denser taxonomic sampling at the species and population levels. Comparative analyses revealed that mitogenomes of Apini and Bombini are the typical type, while those of Meliponini show diverse variations in mitogenome sizes and organizations. Large inverted repeats (IRs) cause significant gene rearrangements of protein coding genes (PCGs) and rRNAs in Indo-Malay/Australian stingless bee species. Molecular evolution analyses showed that the lineage with IRs have lower d N/ d S ratios for PCGs than lineages without IRs, indicating potential effects of IRs on the evolution of mitochondrial genes. The finding of IRs and different patterns of gene rearrangements suggested that Meliponini is a hotspot in mitogenome evolution. Unlike conserved PCGs and rRNAs whose rearrangements were found only in the mentioned lineages within Meliponini, tRNA rearrangements are common across all three tribes of social bees, and are significant even at the species level, indicating that comprehensive sampling is needed to fully understand the patterns of tRNA rearrangements, and their implications for phylogenetic inference.

Mitochondrial genome rearrangements and phylogenomics of the Hymenoptera (Insecta) using an expanded taxon sample

The order Hymenoptera is one of the most species-rich insect orders, with more than 150,000 described extant species. Many hymenopteran insects have very different mitochondrial genome (mitogenome) organizations compared to the putative ancestral organization of insects. In this study, we sequenced 18 mitogenomes of representatives in the order Hymenoptera to increase taxonomic sampling. A total of 475 species were used in phylogenetic analyses, including 18 new mitogenomes and 457 existing mitogenomes. Using a site-heterogeneous model, Bayesian's inference from amino acid data yielded more resolved relationships among Hymenoptera than maximum-likelihood analysis and coalescent-based species analyses. The monophyly of Symphyta was not supported. The Xyeloidea was the earliest branching clade in the Hymenoptera. The Orussoidea was closely related to Apocrita. Within Apocrita, the Parasitoida was non-monophyletic. The monophyly of most Parasitoida superfamilies received strong support. The Proctotrupomorpha clade was supported in Bayesian's analysis. The Apoidea was monophyletic when excluding Ampulex compressa from consideration. The superfamilies Vespoidea and Chrysidoidea were found to be non-monophyletic. Comparisons of mitochondrial gene order revealed a higher frequency of gene rearrangement among lineages with a parasitoid lifestyle, particularly prominent in Chalcidoidea. The degree of gene rearrangement ranked second in specific taxa of Cynipoidea and Ichneumonoidea.

Positive Correlation of the Gene Rearrangements and Evolutionary Rates in the Mitochondrial Genomes of Thrips (Insecta: Thysanoptera)

Simple Summary

Aeolothrips, commonly known as banded thrips, is the largest genus of the family Aeolothripidae (predatory thrips). In the current study, we sequenced the mitochondrial genome (mitogenome) of the banded thrip species Aeolothrips xinjiangensis. We found a novel gene arrangement in this mitogenome that has not been reported in Thysanoptera. By comparing the gene order and rearrangement patterns, we found seven identical gene blocks and three identical rearrangement events in two mitogenomes of banded thrips. There was marked variation in the mitochondrial gene order across thrip species, with only two conserved gene blocks shared by all 14 thrips. In addition, we found a positive correlation between the degree of gene rearrangement and evolutionary rate. Our results suggested that the mitogenomes of thrips have tended to be stable since their massive rearrangement.

Abstract

Extensive gene rearrangement is characteristic in the mitogenomes of thrips (Thysanoptera), but the historical process giving rise to the contemporary gene rearrangement pattern remains unclear. To better understand the evolutionary processes of gene rearrangement in the mitogenomes of thrips, we sequenced the mitogenome of the banded thrip species Aeolothrips xinjiangensis. First, we found a novel mitochondrial gene order in this species. This mitogenome is 16,947 bp in length and encodes the typical 37 coding genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) of insects. The gene arrangement was dramatically different from the putative ancestral mitogenome, with 26 genes being translocated, eight of which were inverted. Moreover, we found a novel, conserved gene block, trnC-trnY, which has not been previously reported in the mitogenomes of thrips. With this newly assembled mitogenome, we compared mitogenome sequences across Thysanoptera to assess the evolutionary processes giving rise to the current gene rearrangement pattern in thrips. Seven identical gene blocks were shared by two sequenced banded thrip mitogenomes, while the reversal of ND2 combined with TDRL events resulted in the different gene orders of these two species. In phylogenetic analysis, the monophyly of the suborders and families of Thysanoptera was well supported. Across the gene orders of 14 thrips, only two conserved gene blocks, ATP8-ATP6 and ND4-ND4L, could be found. Correlation analysis showed that the degree of gene rearrangement was positively correlated with the non-synonymous substitution rate in thrips. Our study suggests that the mitogenomes of thrips remain stable over long evolutionary timescales after massive rearrangement during early diversification.

Characterization of the complete mitochondrial genome of a barklouse, Lepinotus sp. (Psocodea: Trogiomorpha: Trogiidae)

Barklice in the genus Lepinotus (Psocoptera: Trogiidae) are small, soft-bodied stored-product pests that are difficult to control. We sequenced and annotated the mitochondrial (mt) genome of Lepinotus sp. The mt genome of Lepinotus sp. is 16,299 bp in size with 74.4% A + T content. The gene order was highly conserved in some of the Trogimorpha barklice. Two types of tandem repeat units were identified in CR of Lepinotus sp. The phylogenetic analysis showed that Trogiidae species was the sister group to Lepidopsocidae barklice, and the suborder Troctomorpha was polyphyletic.

Tanyptera (Tanyptera) hebeiensis Yang et Yang (Diptera: Tipulidae) newly recorded from Shandong, China: sequencing and phylogenetic analysis of the mitochondrial genome

The genus Tanyptera Latreille, 1804 is recorded from Shandong Province, China for the first time with T. (T.) hebeiensis Yang et Yang, 1988 found in Mount Kunyu, Shandong. In this study, we report the complete mitochondrial genome sequence of T. (T.) hebeiensis, representing the first mitochondrial genome of the subfamily Ctenophorinae (Diptera: Tipulidae), which is a circular molecule of 15,888 bp with an AT content of 77.6%. The mitochondrial genome contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding region. Gene overlaps are found at nine gene junctions, ranging from 1 to 8 bp in length. The canonical mitochondrial start codons for invertebrate mitochondrial genomes are found in 12 PCGs, except for COI which uses the uncanonical start codons TCG. Stop codons of 10 PCGs are invariably complete TAA and TAG, while COII, ND4, and ND5 end with a single thymine stop codon. Phylogenetic analysis reveals that the Pediciidae is a sister group to the remaining Tipuloidea, the Cylindrotomidae has a sister-group relationship with the Tipulidae, and the Limoniidae is not a monophyletic clade.

Sequence analysis of mitochondrial genome of the false and phantom crane-fly Ptychoptera qinggouensis Kang, Yao and Yang, 2013 (Diptera, Ptychopteridae)

The genus Ptychoptera Meigen, 1803 is the largest genus of the family Ptychopteridae with 78 known species. In this study, we report a nearly complete mitochondrial (mt) genome of this genus, which is a circular molecule of more than 15,028 bp. The mt genome contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a non-coding region. The overall base composition is A (38.1%), T (36.7%), C (14.9%), and G (10.4%), with an AT content of 74.8%. The AT content of N-strand genes (75.7%) is higher than that of the J-strand genes (71.7%). Phylogenetic analysis reveals that the monophyly of Ptychopteridae, Bibiomorpha, Tipulomorpha and Brachycera are strongly supported, and the sister group relationship between Tanyderidae and Ptychopteridae is not supported.

Mitochondrial Genomes Provide Insights into the Phylogeny of Culicomorpha (Insecta: Diptera)

Culicomorpha is a monophyletic group containing most bloodsucking lower dipterans, including many important vectors of pathogens. However, the higher-level phylogenetic relationships within Culicomorpha are largely unresolved, with multiple competing hypotheses based on molecular sequence data. Here we sequenced four nearly complete mitochondrial (mt) genomes representing four culicomorph families, and combined these new data with published mt genomes to reconstruct the phylogenetic relationships of all eight extant culicomorph families. We estimated phylogenies using four datasets and three methods. We also used four-cluster likelihood mapping to study potential incongruent topologies supported by the different datasets and phylogenetic questions generated by the previous studies. The results showed that a clade containing Ceratopogonidae, Thaumaleidae and Simuliidae was the sister group to all other Culicomorpha; in another clade, the Dixidae was basal to the remaining four families; Chaoboridae, Corethrellidae and Culicidae formed a monophyletic group and the Chironomidae was the sister group to this clade; Culicidae and Corethrellidae were sister groups in all trees. Our study provides novel mt genome data in Culicomorpha for three new family representatives, and the resulting mt phylogenomic analysis helps to resolve the phylogeny and taxonomy of Culicomorpha.

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.

Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera)

Background

Mitochondrial genome (mt-genome) data can potentially return artefactual relationships in the higher-level phylogenetic inference of insects due to the biases of accelerated substitution rates and compositional heterogeneity. Previous studies based on mt-genome data alone showed a paraphyly of Cimicomorpha (Insecta, Hemiptera) due to the positions of the families Tingidae and Reduviidae rather than the monophyly that was supported based on morphological characters, morphological and molecular combined data and large scale molecular datasets. Various strategies have been proposed to ameliorate the effects of potential mt-genome biases, including dense taxon sampling, removal of third codon positions or purine-pyrimidine coding and the use of site-heterogeneous models. In this study, we sequenced the mt-genomes of five additional Tingidae species and discussed the compositional and mutational rate heterogeneity in mt-genomes and its effect on the phylogenetic inferences of Cimicomorpha by implementing the bias-reduction strategies mentioned above.

Results

Heterogeneity in nucleotide composition and mutational biases were found in mt protein-coding genes, and the third codon exhibited high levels of saturation. Dense taxon sampling of Tingidae and Reduviidae and the other common strategies mentioned above were insufficient to recover the monophyly of the well-established group Cimicomorpha. When the sites with weak phylogenetic signals in the dataset were removed, the remaining dataset of mt-genomes can support the monophyly of Cimicomorpha; this support demonstrates that mt-genomes possess strong phylogenetic signals for the inference of higher-level phylogeny of this group. Comparison of the ratio of the removal of amino acids for each PCG showed that ATP8 has the highest ratio while CO1 has the lowest. This pattern is largely congruent with the evolutionary rate of 13 PCGs that ATP8 represents the highest evolutionary rate, whereas CO1 appears to be the lowest. Notably, the value of Ka/Ks ratios of all PCGs is less than 1, indicating that these genes are likely evolving under purifying selection.

Conclusions

Our results demonstrate that mt-genomes have sites with strong phylogenetic signals for the inference of higher-level phylogeny of Cimicomorpha. Consequently, bioinformatic approaches to removing sites with weak phylogenetic signals in mt-genome without relying on an a priori tree topology would greatly improve this field.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4650-9) contains supplementary material, which is available to authorized users.

The Highly Divergent Mitochondrial Genomes Indicate That the Booklouse, Liposcelis bostrychophila (Psocoptera: Liposcelididae) Is a Cryptic Species

The booklouse, Liposcelis bostrychophila is an important storage pest worldwide. The mitochondrial (mt) genome of an asexual strain (Beibei, China) of the L. bostrychophila comprises two chromosomes; each chromosome contains approximate half of the 37 genes typically found in bilateral animals. The mt genomes of two sexual strains of L. bostrychophila, however, comprise five and seven chromosomes, respectively; each chromosome contains one to six genes. To understand mt genome evolution in L. bostrychophila, and whether L. bostrychophila is a cryptic species, we sequenced the mt genomes of six strains of asexual L. bostrychophila collected from different locations in China, Croatia, and the United States. The mt genomes of all six asexual strains of L. bostrychophila have two chromosomes. Phylogenetic analysis of mt genome sequences divided nine strains of L. bostrychophila into four groups. Each group has a distinct mt genome organization and substantial sequence divergence (48.7–87.4%) from other groups. Furthermore, the seven asexual strains of L. bostrychophila, including the published Beibei strain, are more closely related to two other species of booklice, L. paeta and L. sculptilimacula, than to the sexual strains of L. bostrychophila. Our results revealed highly divergent mt genomes in the booklouse, L. bostrychophila, and indicate that L. bostrychophila is a cryptic species.

Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera

BACKGROUND

Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae.

RESULTS

We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny.

CONCLUSIONS

This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.

Nearly complete mitogenome of hairy sawfly, Corynis lateralis (Brullé, 1832) (Hymenoptera: Cimbicidae): rearrangements in the IQM and ARNS1EF gene clusters

The Cimbicidae is a small family of the primitive and relatively less diverse suborder Symphyta (Hymenoptera). Here, nearly complete mitochondrial genome (mitogenome) of hairy sawfly, Corynis lateralis (Hymenoptera: Cimbicidae) was sequenced using next generation sequencing and comparatively analysed with the mitogenome of Trichiosoma anthracinum. The sequenced length of C. lateralis mitogenome was 14,899 bp with an A+T content of 80.60%. All protein coding genes (PCGs) are initiated by ATN codons and all are terminated with TAR or T- stop codon. All tRNA genes preferred usual anticodons. Compared with the inferred insect ancestral mitogenome, two tRNA rearrangements were observed in the IQM and ARNS1EF gene clusters, representing a new event not previously reported in Symphyta. An illicit priming of replication and/or intra/inter-mitochondrial recombination and TDRL seem to be responsible mechanisms for the rearrangement events in these gene clusters. Phylogenetic analyses confirmed the position of Corynis within Cimbicidae and recovered a relationship of Tenthredinoidea + (Cephoidea + Orussoidea) in Symphyta.

New Mitogenomes of Two Chinese Stag Beetles (Coleoptera, Lucanidae) and Their Implications for Systematics

Although conspicuous and well-studied, stag beetles have been slow to join the genomic era. In this study, mitochondrial genomes of two stag beetles, Sinodendron yunnanense and Prosopocoilus confucius, are sequenced for the first time. Both of their genomes consisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. The mitogenome of S. yunnanense was 16,921 bp in length, and P. confucius was 16,951 bp. The location of the gene trnL(UUR), between the A + T-rich and control region in S. yunnanense, is the first observed in Lucanidae. In P. confucius, an unexpected noncoding region of 580 bp was discovered. Maximum likelihood and Bayesian inference on the 13 mitochondrial PCGs were used to infer the phylogenetic relationships among 12 representative stag beetles and three scarab beetles. The topology of the two phylogenetic trees was almost identical: S. yunnanense was recovered as the most basal Lucanid, and the genus Prosopocoilus was polyphyletic due to P. gracilis being recovered sister to the genera Dorcus and Hemisodorcus. The phylogenetic results, genetic distances and mitogenomic characteristics call into question the cohesion of the genus Prosopocoilus. The genetic resources and findings herein attempts to redress understudied systematics and mitogenomics of the stag beetles.

The first complete mitochondrial genome of Dacus longicornis (Diptera: Tephritidae) using next-generation sequencing and mitochondrial genome phylogeny of Dacini tribe

The genus Dacus is one of the most economically important tephritid fruit flies. The first complete mitochondrial genome (mitogenome) of Dacus species - D. longicornis was sequenced by next-generation sequencing in order to develop the mitogenome data for this genus. The circular 16,253 bp mitogenome is the typical set and arrangement of 37 genes present in the ancestral insect. The mitogenome data of D. longicornis was compared to all the published homologous sequences of other tephritid species. We discovered the subgenera Bactrocera, Daculus and Tetradacus differed from the subgenus Zeugodacus, the genera Dacus, Ceratitis and Procecidochares in the possession of TA instead of TAA stop codon for COI gene. There is a possibility that the TA stop codon in COI is the synapomorphy in Bactrocera group in the genus Bactrocera comparing with other Tephritidae species. Phylogenetic analyses based on the mitogenome data from Tephritidae were inferred by Bayesian and Maximum-likelihood methods, strongly supported the sister relationship between Zeugodacus and Dacus.

The mitochondrial genome of booklouse, Liposcelis sculptilis (Psocoptera: Liposcelididae) and the evolutionary timescale of Liposcelis

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. To date, the complete mt genome has only been determined for four species of Liposcelis, a genus with economic importance, including L. entomophila, L. decolor, L. bostrychophila, and L. paeta. They belong to A, B, or D group of Liposcelis, respectively. Unlike most bilateral animals, L. bostrychophila, L. entomophila and L. paeta have a bitipartite mt genome with genes on two chromosomes. However, the mt genome of L. decolor has the typical mt chromosome of bilateral animals. Here, we sequenced the mt genome of L. sculptilis, and identified 35 genes, which were on a single chromosome. The mt genome fragmentation is not shared by the D group of Liposcelis and the single chromosome of L. sculptilis differed from those of booklice known in gene content and gene arrangement. We inferred that different evolutionary patterns and rate existed in Liposcelis. Further, we reconstructed the evolutionary history of 21 psocodean taxa with phylogenetic analyses, which suggested that Liposcelididae and Phthiraptera have evolved 134 Ma and the sucking lice diversified in the Late Cretaceous.

Comparative Mt Genomics of the Tipuloidea (Diptera: Nematocera: Tipulomorpha) and Its Implications for the Phylogeny of the Tipulomorpha

A traditionally controversial taxon, the Tipulomorpha has been frequently discussed with respect to both its familial composition and relationships with other Nematocera. The interpretation of internal relationships within the Tipuloidea, which include the Tipulidae sensu stricto, Cylindrotomidae, Pediciidae and Limoniidae, is also problematic. We sequenced the first complete mitochondrial (mt) genome of Symplecta hybrida (Meigen, 1804), which belongs to the subfamily Chioneinae of family Limoniidae, and another five nearly complete mt genomes from the Tipuloidea. We did a comparative analysis of these mt genomics and used them, along with some other representatives of the Nematocera to construct phylogenetic trees. Trees inferred by Bayesian methods strongly support a sister-group relationship between Trichoceridae and Tipuloidea. Tipulomorpha are not supported as the earliest branch of the Diptera. Furthermore, phylogenetic trees indicate that the family Limoniidae is a paraphyletic group.

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.

Capturing the Phylogeny of Holometabola with Mitochondrial Genome Data and Bayesian Site-Heterogeneous Mixture Models

After decades of debate, a mostly satisfactory resolution of relationships among the 11 recognized holometabolan orders of insects has been reached based on nuclear genes, resolving one of the most substantial branches of the tree-of-life, but the relationships are still not well established with mitochondrial genome data. The main reasons have been the absence of sufficient data in several orders and lack of appropriate phylogenetic methods that avoid the systematic errors from compositional and mutational biases in insect mitochondrial genomes. In this study, we assembled the richest taxon sampling of Holometabola to date (199 species in 11 orders), and analyzed both nucleotide and amino acid data sets using several methods. We find the standard Bayesian inference and maximum-likelihood analyses were strongly affected by systematic biases, but the site-heterogeneous mixture model implemented in PhyloBayes avoided the false grouping of unrelated taxa exhibiting similar base composition and accelerated evolutionary rate. The inclusion of rRNA genes and removal of fast-evolving sites with the observed variability sorting method for identifying sites deviating from the mean rates improved the phylogenetic inferences under a site-heterogeneous model, correctly recovering most deep branches of the Holometabola phylogeny. We suggest that the use of mitochondrial genome data for resolving deep phylogenetic relationships requires an assessment of the potential impact of substitutional saturation and compositional biases through data deletion strategies and by using site-heterogeneous mixture models. Our study suggests a practical approach for how to use densely sampled mitochondrial genome data in phylogenetic analyses.

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.

Mitochondrial genomes of praying mantises (Dictyoptera, Mantodea): rearrangement, duplication, and reassignment of tRNA genes

Insect mitochondrial genomes (mitogenomes) contain a conserved set of 37 genes for an extensive diversity of lineages. Previously reported dictyopteran mitogenomes share this conserved mitochondrial gene arrangement, although surprisingly little is known about the mitogenome of Mantodea. We sequenced eight mantodean mitogenomes including the first representatives of two families: Hymenopodidae and Liturgusidae. Only two of these genomes retain the typical insect gene arrangement. In three Liturgusidae species, the trnM genes have translocated. Four species of mantis (Creobroter gemmata, Mantis religiosa, Statilia sp., and Theopompa sp.-HN) have multiple identical tandem duplication of trnR, and Statilia sp. additionally includes five extra duplicate trnW. These extra trnR and trnW in Statilia sp. are erratically arranged and form another novel gene order. Interestingly, the extra trnW is converted from trnR by the process of point mutation at anticodon, which is the first case of tRNA reassignment for an insect. Furthermore, no significant differences were observed amongst mantodean mitogenomes with variable copies of tRNA according to comparative analysis of codon usage. Combined with phylogenetic analysis, the characteristics of tRNA only possess limited phylogenetic information in this research. Nevertheless, these features of gene rearrangement, duplication, and reassignment provide valuable information toward understanding mitogenome evolution in insects.

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.

Mitochondrial genomes of four katydids (Orthoptera: Phaneropteridae): New gene rearrangements and their phylogenetic implications

Phaneropteridae is a family of Orthoptera that displays an amazing amount of diversity in terms of both forms and species. We sequenced the mitochondrial genomes (mitogenomes) of two bush katydids: Ruidocollaris obscura and Kuwayamaea brachyptera (Phaneropterinae), and two true katydids: Orophyllus montanus and Phyllomimus detersus (Pseudophyllinae), to obtain further insight into the characteristics of the katydid mitogenomes and to investigate the taxonomic status of subfamily Pseudophyllinae and the diversity of gene arrangements among Phaneropteridae. The following general genomic characteristics were observed in the four katydids: a longer length of the mitogenomes (16,007bp-16,667bp) compared with Caelifera, abundant intergenic spacers, and accepted atypical initiation codons (GTG and TTG, found in cox1, nad1 and nad2). A new orientation of the gene arrangement "trnM-trnI-trnQ" was identified in P. detersus, which is the first representative of Polyneoptera found to carry this gene cluster. Large identical fragments (492bp) were detected in control region 1 (CR1) and control region 2 (CR2) of R. obscura. The high similarity of the duplicated CRs is likely due to a recent gene duplication or concerted evolution. Analyses of the duplicated CRs revealed one conserved stem-loop (on the N-strand) located in the identical sequences of both CRs that might be linked to replication initiation. Phylogenetic analyses based on 13 protein-coding genes and 2 ribosomal RNA genes from 20 Ensiferan species yielded the identical topologies between two different methods (maximum likelihood and bayesian inference). The newly sequenced Pseudophyllinae species was placed as the sister group of Phaneropterinae, and Mecopodinae clustered with Pseudophyllinae+Phaneropterinae. Additionally, we speculate that the species in Ruidocollaris and Sinochlora, as well as their closely related genera, may have undergone numerous rearrangement events.

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.

The complete mitochondrial genome of Aix galericulata and Tadorna ferruginea: bearings on their phylogenetic position in the Anseriformes

Aix galericulata and Tadorna ferruginea are two Anatidae species representing different taxonomic groups of Anseriformes. We used a PCR-based method to determine the complete mtDNAs of both species, and estimated phylogenetic trees based on the complete mtDNA alignment of these and 14 other Anseriforme species, to clarify Anseriform phylogenetics. Phylogenetic trees were also estimated using a multiple sequence alignment of three mitochondrial genes (Cyt b, ND2, and COI) from 68 typical species in GenBank, to further clarify the phylogenetic relationships of several groups among the Anseriformes. The new mtDNAs are circular molecules, 16,651 bp (Aix galericulata) and 16,639 bp (Tadorna ferruginea) in length, containing the 37 typical genes, with an identical gene order and arrangement as those of other Anseriformes. Comparing the protein-coding genes among the mtDNAs of 16 Anseriforme species, ATG is generally the start codon, TAA is the most frequent stop codon, one of three, TAA, TAG, and T-, commonly observed. All tRNAs could be folded into canonical cloverleaf secondary structures except for tRNASer (AGY) and tRNALeu (CUN), which are missing the "DHU" arm.Phylogenetic relationships demonstrate that Aix galericula and Tadorna ferruginea are in the same group, the Tadorninae lineage, based on our analyses of complete mtDNAs and combined gene data. Molecular phylogenetic analysis suggests the 68 species of Anseriform birds be divided into three families: Anhimidae, Anatidae, and Anseranatidae. The results suggest Anatidae birds be divided into five subfamilies: Anatinae, Tadorninae, Anserinae, Oxyurinae, and Dendrocygninae. Oxyurinae and Dendrocygninae should not belong to Anserinae, but rather represent independent subfamilies. The Anatinae includes species from the tribes Mergini, Somaterini, Anatini, and Aythyini. The Anserinae includes species from the tribes Anserini and Cygnini.

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

BACKGROUND

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

RESULTS

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

CONCLUSIONS

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

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.

Complete mitochondrial genome of the flat bug Brachyrhynchus hsiaoi (Hemiptera: Aradidae)

The mitochondrial genome of a flat bug, Brachyrhynchus hsiaoi (Blöte), is a typical circular DNA molecule of 15,250 bp with 37 genes and 70.4% A + T content. The gene order is different from that of the putative ancestral arrangement of insects with a tRNA gene rearrangement, trnQ-trnI. This rearrangement has been found in other sequenced flat bugs and is likely synapomorphic for the Aradidae or some subgroup within this family. Ten protein-coding genes start with ATN codon and others use TTG. All the 22 tRNAs, ranging from 61 to 70 bp, have the clover-leaf structure except for the dihydrouridine (DHU) arm of trnS1 forms a simple loop. The sizes of the large and small ribosomal RNA genes are 1245 and 808 bp, respectively. The control region is located between rrnS and trnQ with 703 bp in length and 69.8% A + T content.

Substantial variation in the extent of mitochondrial genome fragmentation among blood-sucking lice of mammals

Blood-sucking lice of humans have extensively fragmented mitochondrial (mt) genomes. Human head louse and body louse have their 37 mt genes on 20 minichromosomes. In human pubic louse, the 34 mt genes known are on 14 minichromosomes. To understand the process of mt genome fragmentation in the blood-sucking lice of mammals, we sequenced the mt genomes of the domestic pig louse, Haematopinus suis, and the wild pig louse, H. apri, which diverged from human lice approximately 65 Ma. The 37 mt genes of the pig lice are on nine circular minichromosomes; each minichromosome is 3–4 kb in size. The pig lice have four genes per minichromosome on average, in contrast to two genes per minichromosome in the human lice. One minichromosome of the pig lice has eight genes and is the most gene-rich minichromosome found in the sucking lice. Our results indicate substantial variation in the rate and extent of mt genome fragmentation among different lineages of the sucking lice.

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.