No abstract

Hemiptera, the largest non-holometabolous order of insects, represents approximately 7% of metazoan diversity. With extraordinary life histories and highly specialized morphological adaptations, hemipterans have exploited diverse habitats and food sources through approximately 300 Myr of evolution. To elucidate the phylogeny and evolutionary history of Hemiptera, we carried out the most comprehensive mitogenomics analysis on the richest taxon sampling to date covering all the suborders and infraorders, including 34 newly sequenced and 94 published mitogenomes. With optimized branch length and sequence heterogeneity, Bayesian analyses using a site-heterogeneous mixture model resolved the higher-level hemipteran phylogeny as (Sternorrhyncha, (Auchenorrhyncha, (Coleorrhyncha, Heteroptera))). Ancestral character state reconstruction and divergence time estimation suggest that the success of true bugs (Heteroptera) is probably due to angiosperm coevolution, but key adaptive innovations (e.g. prognathous mouthpart, predatory behaviour, and haemelytron) facilitated multiple independent shifts among diverse feeding habits and multiple independent colonizations of aquatic habitats.

No abstract

Organelle genome fragmentation has been found in a wide range of eukaryotic lineages; however, its use in phylogenetic reconstruction has not been demonstrated. We explored the use of mitochondrial (mt) genome fragmentation in resolving the controversial suborder-level phylogeny of parasitic lice (order Phthiraptera). There are ∼5,000 species of parasitic lice in four suborders (Amblycera, Ischnocera, Rhynchophthirina and Anoplura), which infest mammals and birds. The phylogenetic relationships among these suborders are unresolved despite decades of studies. We sequenced the mt genomes of eight species of parasitic lice and compared them with 17 other species of parasitic lice sequenced previously. We found that the typical single-chromosome mt genome is retained in the lice of birds but fragmented into many minichromosomes in the lice of eutherian mammals. The shared derived feature of mt genome fragmentation unites the eutherian mammal lice of Ischnocera (family Trichodectidae) with Anoplura and Rhynchophthirina to the exclusion of the bird lice of Ischnocera (family Philopteridae). The novel clade, namely Mitodivisia, is also supported by phylogenetic analysis of mt genome and cox1 gene sequences. Our results demonstrate, for the first time, that organelle genome fragmentation is informative for resolving controversial high-level phylogenies.

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.

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.

Many of true bugs are important insect pests to cultivated crops and some are important vectors of human diseases, but few cladistic analyses have addressed relationships among the seven infraorders of Heteroptera. The Enicocephalomorpha and Nepomorpha are consider the basal groups of Heteroptera, but the basal-most lineage remains unresolved. Here we report the mitochondrial genome of the unique-headed bug Stenopirates sp., the first mitochondrial genome sequenced from Enicocephalomorpha. The Stenopirates sp. mitochondrial genome is a typical circular DNA molecule of 15, 384 bp in length, and contains 37 genes and a large non-coding fragment. The gene order differs substantially from other known insect mitochondrial genomes, with rearrangements of both tRNA genes and protein-coding genes. The overall AT content (82.5%) of Stenopirates sp. is the highest among all the known heteropteran mitochondrial genomes. The strand bias is consistent with other true bugs with negative GC-skew and positive AT-skew for the J-strand. The heteropteran mitochondrial atp8 exhibits the highest evolutionary rate, whereas cox1 appears to have the lowest rate. Furthermore, a negative correlation was observed between the variation of nucleotide substitutions and the GC content of each protein-coding gene. A microsatellite was identified in the putative control region. Finally, phylogenetic reconstruction suggests that Enicocephalomorpha is the sister group to all the remaining Heteroptera.

The higher‐level phylogeny of Pentatomomorpha, the second largest infraorder of true bugs (Hemiptera: Heteroptera), which includes many important agriculture and forestry pests, has been debated for decades. To investigate the phylogeny and evolutionary history of Pentatomomorpha, we assembled new mitochondrial genomes for 46 species through next‐generation sequencing of pooled genomic DNA. Based on a much broader taxon sampling than available previously, Bayesian analyses using a site‐heterogeneous mixture model (CAT+GTR) resolved the higher‐level phylogeny of Pentatomomorpha as (Aradoidea + (Pentatomoidea + (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)))). There was a transition from trnT/trnP to trnP/trnT in the common ancestor of Pyrrhocoroidea, which indicates that this gene rearrangement could be an autapomorphy for Pyrrhocoroidea. Divergence time analyses estimated that Pentatomomorpha originated c. 242 Ma in the Middle Triassic, and most of the recognized superfamilies originated during the Middle Jurassic to Early Cretaceous. The diversification of families within Pentatomomorpha largely coincided with the radiation of angiosperms during the Early Cretaceous.