Long-branch attraction and the phylogeny of true water bugs (Hemiptera: Nepomorpha) as estimated from mitochondrial genomes

The complete mitochondrial genome of Sigara lateralis (Leach, 1817) (Nepomorpha: Corixidae)

The Sigara lateralis (Leach, 1817) is a small aquatic insect belonging to the family Corixidae. The study aims to reveal the genomic data of the mitochondrial genome of S. lateralis. The length of its circular mitochondrial genome is 15,725 bp long with an A + T bias (75.96%). The mitogenome comprises 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two rRNA genes. The phylogenetic analyses showed that the S. lateralis is the closest to S. septemlineata. These findings will help the conservation of Corixidae from the perspective of genetic evolution.

Genome-scale approach to reconstructing the phylogenetic tree of psyllids (superfamily Psylloidea) with account of systematic bias

Psyllids (class Insecta: order Hemiptera: superfamily Psylloidea) are a taxonomically and phylogenetically challenging clade. Recent studies have largely advanced the phylogeny of this group, yet the family-level relationships among Aphalaridae, Carsidaridae, and others remain unresolved. Genome-scale phylogenetic analysis is known to provide a finer resolution for problems like that. However, such phylogenomics also introduces new problems: incorrect trees with high confidence yielded due to systematic error (bias). Here we addressed these issues using hundreds of single-copy orthologous (SCO) genes in psyllid transcriptomes and genomes. Our analyses revealed conflicts between the nucleotide-based and amino-acid-based phylogenetic trees. While the nucleotide-based phylogeny strongly supported the (Aphalaridae + Carsidaridae) + Others relationship, the amino-acid-based one recovered Aphalaridae + (Carsidaridae + Others) with 100% support. Further inspection revealed significant compositional heterogeneity in nucleotide sequences for 67% of SCO genes, but not in the corresponding translated amino acid sequences. We then used different strategies to combat this compositional bias, and found that using the RY-coding strategy (coding the standard nucleotides as purines and pyrimidines) the nucleotide-based phylogeny became consistent with the amino-acid-based one. We further applied RY-coding to a published concatenated nucleotide dataset and recovered the Aphalaridae monophyly (which is refuted by the original literature on non-recoded sequences) at the base of psyllid tree. Moreover, it was found that variations in evolutionary rate could lead to errors in nucleotide-based phylogeny. The fast-evolving Heteropsylla cubana (Psyllidae: Ciriacreminae) was incorrectly placed within the subfamily Psyllinae. This bias can be avoided by using data removal or RY-coding strategies. Together, our results strongly support the family relationship of Aphalaridae + (Carsidaridae + Others), and show that the amino-acid-based concatenation analysis is more robust than nucleotide-based one. Future phylogenomic analysis of psyllid nucleotide sequences should take into account methods such as the RY-coding scheme to address potential systematic biases arising from composition and rate heterogeneities.

Three New Species of Velvety Shore Bugs (Hemiptera: Heteroptera: Ochteroidea) From Mid-Cretaceous Kachin Amber Shed Light on the Evolution of Rostrum Length in Ochteroidea

Two new genera of velvety shore bugs, Arcochterus Zhang, Ren and Yao gen. nov. and Parvochterus Zhang, Ren and Yao gen. nov. are described, with three new species between them—Arcochterus fuscus Zhang, Ren and Yao sp. nov., Parvochterus reticulatus Ren and Yao sp. nov., and P. lanceolarus Zhang, Ren and Yao sp. nov. Based on the combination of fossil and extant taxa, a cladistic analysis is conducted to confirm the phylogenetic position of these species and allows reconstruction of the inter-genus relationships within the superfamily Ochteroidea. Major conclusions of the phylogenetic analysis: (1) these new species and Grimaldinia pronotalis belong to Ochteridae; (2) Ochteroidea is a monophyletic group, Ochteridae and Gelastocoridae are sister group and monophyletic, respectively. (3) The ancestral character state reconstruction (ACSR) shows that the length of the rostrum has occurred in at least three independent transitions during the evolution of the superfamily Ochteroidea.

Providing a phylogenetic framework for trait-based analyses in brown algae: Phylogenomic tree inferred from 32 nuclear protein-coding sequences

In the study of the evolution of biological complexity, a reliable phylogenetic framework is needed. Many attempts have been made to resolve phylogenetic relationships between higher groups (i.e., interordinal) of brown algae (Phaeophyceae) based on molecular evidence, but most of these relationships remain unclear. Analyses based on small multi-gene data (including chloroplast, mitochondrial and nuclear sequences) have yielded inconclusive and sometimes contradictory results. To address this problem, we have analyzed 32 nuclear protein-coding sequences in 39 Phaeophycean species belonging to eight orders. The resulting nuclear-based phylogenomic trees provide virtually full support for the phylogenetic relationships within the studied taxa, with few exceptions. The relationships largely confirm phylogenetic trees based on nuclear, chloroplast and mitochondrial sequences, except for the placement of the Sphacelariales with weak bootstrap support. Our study indicates that nuclear protein-coding sequences provide significant support to conclusively resolve phylogenetic relationships among Phaeophyceae, and may be a powerful approach to fully resolve interordinal relationships with increased taxon sampling.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

A mitochondrial genome phylogeny of voles and lemmings (Rodentia: Arvicolinae): Evolutionary and taxonomic implications

Arvicolinae is one of the most impressive placental radiations with over 150 extant and numerous extinct species that emerged since the Miocene in the Northern Hemisphere. The phylogeny of Arvicolinae has been studied intensively for several decades using morphological and genetic methods. Here, we sequenced 30 new mitochondrial genomes to better understand the evolutionary relationships among the major tribes and genera within the subfamily. The phylogenetic and molecular dating analyses based on 11,391 bp concatenated alignment of protein-coding mitochondrial genes confirmed the monophyly of the subfamily. While Bayesian analysis provided a high resolution across the entire tree, Maximum Likelihood tree reconstruction showed weak support for the ordering of divergence and interrelationships of tribal level taxa within the most ancient radiation. Both the interrelationships among tribes Lagurini, Ellobiusini and Arvicolini, comprising the largest radiation and the position of the genus Dinaromys within it also remained unresolved. For the first time complex relationships between genus level taxa within the species-rich tribe Arvicolini received full resolution. Particularly Lemmiscus was robustly placed as sister to the snow voles Chionomys in the tribe Arvicolini in contrast with a long-held belief of its affinity with Lagurini. Molecular dating of the origin of Arvicolinae and early divergences obtained from the mitogenome data were consistent with fossil records. The mtDNA estimates for putative ancestors of the most genera within Arvicolini appeared to be much older than it was previously proposed in paleontological studies.

Three Complete Mitochondrial Genomes of Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis (Insecta: Phasmatodea) and Their Phylogeny

Insects of the order Phasmatodea are mainly distributed in the tropics and subtropics and are best known for their remarkable camouflage as plants. In this study, we sequenced three complete mitochondrial genomes from three different families: Orestes guangxiensis, Peruphasma schultei, and Phryganistria guangxiensis. The lengths of the three mitochondrial genomes were 15,896 bp, 16,869 bp, and 17,005 bp, respectively, and the gene composition and structure of the three stick insects were identical to those of the most recent common ancestor of insects. The phylogenetic relationships among stick insects have been chaotic for a long time. In order to discuss the intra- and inter-ordinal relationship of Phasmatodea, we used the 13 protein-coding genes (PCGs) of 85 species for maximum likelihood (ML) and Bayesian inference (BI) analyses. Results showed that the internal topological structure of Phasmatodea had a few differences in both ML and BI trees and long-branch attraction (LBA) appeared between Embioptera and Zoraptera, which led to a non-monophyletic Phasmatodea. Consequently, after removal of the Embioptera and Zoraptera species, we re-performed ML and BI analyses with the remaining 81 species, which showed identical topology except for the position of Tectarchus ovobessus (Phasmatodea). We recovered the monophyly of Phasmatodea and the sister-group relationship between Phasmatodea and Mantophasmatodea. Our analyses also recovered the monophyly of Heteropterygidae and the paraphyly of Diapheromeridae, Phasmatidae, Lonchodidae, Lonchodinae, and Clitumninae. In this study, Peruphasma schultei (Pseudophasmatidae), Phraortes sp. YW-2014 (Lonchodidae), and species of Diapheromeridae clustered into the clade of Phasmatidae. Within Heteropterygidae, O. guangxiensis was the sister clade to O. mouhotii belonging to Dataminae, and the relationship of (Heteropteryginae + (Dataminae + Obriminae)) was recovered.

New Data, Old Story: Molecular Data Illuminate the Tribal Relationships among Rove Beetles of the Subfamily Staphylininae (Coleoptera: Staphylinidae)

The megadiverse subfamily Staphylininae traditionally belonged to the best-defined rove beetle taxa, but the advent of molecular phylogenetics in the last decade has brought turbulent changes to the group's classification. Here, we reevaluate the internal relationships among the tribes of Staphylininae by implementing tree inference methods that suppress common sources of systematic error. In congruence with morphological data, and in contrast to some previous phylogenetic studies, we unambiguously recover Staphylininae and Paederinae as monophyletic in the traditional sense. We show that the recently proposed subfamily Platyprosopinae (Arrowinus and Platyprosopus) is a phylogenetic artefact and reinstate Arrowinus as a member of Arrowinini stat. res. and Platyprosopus as a member of Platyprosopini stat. res. We show that several recent changes to the internal classification of the subfamily are phylogenetically unjustified and systematically unnecessary. We, therefore, reestablish Platyprosopini, Staphylinini, and Xantholinini as tribes within Staphylininae (all stat. res.) and recognize Coomaniini as a tribe (stat. nov.) rather than subfamily. Consequently, the traditional ranks of the subtribes Acylophorina, Afroquediina, Amblyopinina, Antimerina, †Baltognathina, Cyrtoquediina, Erichsoniina, Hyptiomina, Indoquediina, Quediina, and Tanygnathinina are restored (all stat. res.). We review the current classification of Staphylininae and discuss sources of incongruence in multigene phylogenies.

Phylogeny and diversification of the true water bugs (Insecta: Hemiptera: Heteroptera: Nepomorpha)

Climate fluctuations and tectonic reconfigurations associated with environmental changes play large roles in determining patterns of adaptation and diversification, but studies documenting how such drivers have shaped the evolutionary history and diversification dynamics of limnic organisms during the Mesozoic are scarce. Members of the heteropteran infraorder Nepomorpha, or aquatic bugs, are ideal for testing the effects of these determinants on their diversification pulses because most species are confined to aquatic environments during their entire life. The group has a relatively mature taxonomy and is well represented in the fossil record. We investigated the evolution of Nepomorpha based on phylogenetic analyses of morphological and molecular characters sampled from 115 taxa representing all 13 families and approximately 40% of recognized genera. Our results were largely congruent with the phylogenetic relationships inferred from morphology. A divergence dating analysis indicated that Nepomorpha began to diversify in the late Permian (approximately 263 Ma), and diversification analyses suggested that palaeoecological opportunities probably promoted lineage diversification in this group.

Characterization of the Complete Mitochondrial Genome of Harpalus sinicus and Its Implications for Phylogenetic Analyses

In this study, we report the complete mitochondrial genome of Harpalus sinicus (occasionally named as the Chinese ground beetle) which is the first mitochondrial genome for Harpalus. The mitogenome is 16,521 bp in length, comprising 37 genes, and a control region. The A + T content of the mitogenome is as high as 80.6%. A mitochondrial origins of light-strand replication (OL)-like region is found firstly in the insect mitogenome, which can form a stem-loop hairpin structure. Thirteen protein-coding genes (PCGs) share high homology, and all of them are under purifying selection. All tRNA genes (tRNAs) can be folded into the classic cloverleaf secondary structures except tRNA-Ser (GCU), which lacks a dihydrouridine (DHU) stem. The secondary structure of two ribosomal RNA genes (rRNAs) is predicted based on previous insect models. Twelve types of tandem repeats and two stem-loop structures are detected in the control region, and two stem-loop structures may be involved in the initiation of replication and transcription. Additionally, phylogenetic analyses based on mitogenomes suggest that Harpalus is an independent lineage in Carabidae, and is closely related to four genera (Abax, Amara, Stomis, and Pterostichus). In general, this study provides meaningful genetic information for Harpalus sinicus and new insights into the phylogenetic relationships within the Carabidae.

A mitochondrial genome of Micronectidae and implications for its phylogenetic position

The mitochondrial genome (mitogenome) has been extensively used to better understand the phylogenetic relationships within the heteropteran infraorder Nepomorpha (Hemiptera), but no mitogenome in Micronectidae has been sequenced to date. Here we describe the first complete mitogenome of Micronecta sahlbergii (Jakovlev, 1881). The mitogenome is 15,005 bp in size, containing 13 typical PCGs, 22 tRNAs, two rRNAs and a control region (CR). All genes are arranged in the same gene order as the most other known heteropteran mitogenome. The phylogenetic relationships based on mitogenomes using Bayesian inference and Maximum likelihood methods showed that Micronecta sahlbergii was sister to Sigara septemlineata, suggesting that Micronecta sahlbergii belongs to Corixoidea. Corixoidea was basal within Nepomorpha. The PCG12 and PCG12RT matrices of BI and ML analyses yielded the consistent topology, respectively. Whereas there was no consistent conclusions in PCG123 and PCG123RT matrices. Saturation tests showed that PCG12 and PCG12RT were free of saturation in evaluation of transition and transversion substitution, while PCG123 and PCG123RT exhibited a plateau revealing saturation of transition suggesting that the third codon positions of PCGs were not suitable for addressing relationships at the superfamily level in Nepomorpha. So our results supported the phylogenetic analysis of PCG12 and PCG12RT in Nepomorpha.

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.

Complete Mitochondrial Genome of Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) and Phylogenetic Analysis of Pentatomomorpha Species

Dinorhynchus dybowskyi (Hemiptera: Pentatomidae: Asopinae) is used as a biological control agent against various insect pests for its predatory. In the present study, the complete mitochondrial genome (mitogenome) of the species was sequenced using the next-generation sequencing technology. The results showed that the mitogenome is 15,952 bp long, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a control region. Furthermore, the gene order and orientation of this mitogenome are identical to those of most heteropterans. There are 21 intergenic spacers (of length 1-28 bp) and 13 overlapping regions (of length 1-23 bp) throughout the genome. The control region is 1,291 bp long. The start codon of the PCGs is ATN, except cox1 (TTG), and stop codon is TAA, except nad1 (TAG). The 22 tRNAs exhibit a typical cloverleaf secondary structure, except trnS1, which lacks a dihydrouridine (DHU) arm and trnV, where the DHU arm forms a simple loop. The analyses based on nucleotide sequences of the 13 PCGs by Bayesian Inference and maximum likelihood methods. The results support the monophyly of five superfamilies Aradoidea, Pentatomoidea, Pyrrhocoroidea, Lygaeoidea, and Coreoidea. Within Pentatomoidea, the relationship observed is as follows: (Plataspidae + Urostylididae) + (Pentatomidae + (Acanthosomatidae + (Cydnidae + (Scutelleridae + (Dinidoridae + Tessaratomidae))))), and D. dybowskyi was placed in Pentatomidae and close to Eurydema gebleri.

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.

The complete mitochondrial genome of Nepa hoffmanni (Hemiptera: Heteroptera: Nepidae)

The complete mitochondrial genome (mt-genome) of Nepa hoffmanni has been reported in this study. This mitochondrial genome is 15 774 bp long, with an A + T content of 72.04%, containing the typical 37 genes (13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes) and a control region. All genes are arranged in the same gene order as most other known heteropteran mt-genome. This is the second completely sequenced mt-genome from the family Nepidae of Nepomorpha. Bayesian analyses were performed using the mt-genome of Nepa hoffmanni and its relatives, including 17 taxa, showing a reasonable placement of Nepa hoffmanni.

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.

Complete mitochondrial genome of Cacopsylla coccinae (Hemiptera: Psyllidae)

In this study, the first complete mitochondrial genome (mitogenome) sequence of Cacopsylla coccinae was determined by long PCR and primer walking methods. The complete mitochondrial genome is 14,832 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes as well as a control region. The overall base composition of the genome is A (38.16%), T (33.88%), C (17.95%) and G (10.01%). Stop codon was incomplete for coxII gene and ND1 gene. The gene overlaps were suggested between 13 pairs of the contiguous genes in C. coccinae. The mitogenome would contribute to resolving phylogenetic position and interrelationships of Cacopsylla.

Long-branch attraction and the phylogeny of true water bugs (Hemiptera: Nepomorpha) as estimated from mitochondrial genomes

BACKGROUND

Most previous studies of morphological and molecular data have consistently supported the monophyly of the true water bugs (Hemiptera: Nepomorpha). An exception is a recent study by Hua et al. (BMC Evol Biol 9: 134, 2009) based on nine nepomorphan mitochondrial genomes. In the analysis of Hua et al. (BMC Evol Biol 9: 134, 2009), the water bugs in the group Pleoidea formed the sister group to a clade that consisted of Nepomorpha (the remaining true water bugs) + Leptopodomorpha (shore bugs) + Cimicomorpha (assassin bugs and relatives) + Pentatomomorpha (stink bugs and relatives), thereby suggesting that fully aquatic hemipterans evolved independently at least twice. Based on these results, Hua et al. (BMC Evol Biol 9: 134, 2009) elevated the Pleoidea to a new infraorder, the Plemorpha.

RESULTS

Our reanalysis suggests that the lack of support for the monophyly of the true water bugs (including Pleoidea) by Hua et al. (BMC Evol Biol 9: 134, 2009) likely resulted from inadequate taxon sampling. In particular, long-branch attraction (LBA) between the distant outgroup taxa and Pleoidea, as well as LBA among taxa in the ingroup, made Nepomorpha appear to be polyphyletic. We used three complementary strategies to test and alleviate the effects of LBA: (1) the removal of distant outgroups from the analysis; (2) the addition of closely related outgroups; and (3) the addition of a mitochondrial genome from a second family of Pleoidea. We also performed likelihood-ratio tests to examine the support for monophyly of Nepomorpha with different combinations of taxa included in the analysis. Furthermore, we found that specimens of Helotrephes sp. were misidentified as Paraplea frontalis (Fieber, 1844) by Hua et al. (BMC Evol Biol 9: 134, 2009).

CONCLUSIONS

All analyses that included the addition of more taxa significantly and consistently supported the placement of Pleoidea within the Nepomorpha (i.e., supported the monophyly of the traditional true water bugs). Our analyses further support a close relationship between Notonectoidea and Pleoidea within Nepomorpha, and the superfamilies Nepoidea, Ochteroidea, Naucoroidea, and Pleoidea are resolved as monophyletic in all trees with strong support. Our results also confirmed that monophyly of Nepomorpha clearly is not refuted by the mitochondrial genome data.