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

Phylogenetic placement and comparative analysis of the mitochondrial genomes of Idiostoloidea (Hemiptera: Heteroptera)

The classification system and the higher level phylogenetic relationships of Pentatomomorpha, the second largest infraorder of Heteroptera (Insecta: Hemiptera), have been debated and remain controversial over decades. In particular, the placement and phylogenetic relationship of Idiostoloidea are not well resolved, which hampers a better understanding of the evolutionary history of Pentatomomorpha. In this study, for the first time, we reported the complete mitochondrial genome for two narrowly distributed families of Idiostoloidea (including Idiostolidae and Henicocoridae), respectively. The length of the mitochondrial genome of Monteithocoris hirsutus and Henicocoris sp. is 16,632 and 16,013 bp, respectively. The content of AT is ranging from 75.15% to 80.48%. The mitogenomic structure of Idiostoloidea is highly conservative and there are no gene arrangements. By using the Bayesian inference, maximum likelihood, and Bayesian site-heterogeneous mixture model, we inferred the phylogenetic relationships within Pentatomomorpha and estimated their divergence times based on concatenated mitogenomes and nuclear ribosomal genes. Our results support the classification system of six superfamilies within Pentatomomorpha and confirm the monophyletic groups of each superfamily, with the following phylogenetic relationships: (Aradoidea + (Pentatomoidea + (Idiostoloidea + (Coreoidea + (Pyrrhocoroidea + Lygaeoidea))))). Furthermore, estimated divergence times revealed that most pentatomomorphan superfamilies and families diverged during the Late Jurassic to Early Cretaceous, which coincides with the explosive radiation of angiosperms.

Characterization of the Complete Mitochondrial Genome of Schizothorax kozlovi (Cypriniformes, Cyprinidae, Schizothorax) and Insights into the Phylogenetic Relationships of Schizothorax

Schizothorax kozlovi is an endemic and vulnerable fish species found in the upper Yangtze River in China. Over the past few years, the population resources of S. kozlovi have been nearly completely depleted owing to multiple contributing threats. While the complete mitochondrial genomes serve as important molecular markers for phylogenetic and genetic studies, the mitochondrial genome of S. kozlovi has still received little attention. In this study, we analyzed the characterization of the mitochondrial genome of S. kozlovi and investigated the phylogenetic relationships of Schizothorax. The complete mitochondrial genome of S. kozlovi was 16,585 bp in length, which contained thirty-seven genes (thirteen protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), twenty-two transfer RNA genes (tRNAs)) and two non-coding regions for the origin of light strand (OL) and the control region (CR). There were nine overlapping regions and seventeen intergenic spacers regions in the mitochondrial genome. The genome also showed a bias towards A + T content (55.01%) and had a positive AT-skew (0.08) and a negative GC-skew (-0.20). All the PCGs employed the ATG or GTG as the start codon and TAA, TAG, or single T as the stop codon. Additionally, all of the tRNAs displayed a typical cloverleaf secondary structure, except trnS1 which lacked the D arm. The phylogenetic analysis, based on the maximum likelihood (ML) and Bayesian inference (BI) methods, revealed that the topologies of the phylogenetic tree divided the Schizothorax into four clades and did not support the classification of Schizothorax based on morphology. The phylogenetic status of S. kozlovi was closely related to that of S. chongi. The present study provides valuable genomic information for S. kozlovi and new insights in phylogenetic relationships of Schizothorax. These data could also offer fundamental references and guidelines for the management and conservation of S. kozlovi and other species of Schizothorax.

The first complete mitochondrial genome sequences for Ulidiidae and phylogenetic analysis of Diptera

BACKGROUND

Tetanops sintenisi is a pest that mainly damages the root of quinoa (Chenopodium quinoa) and it is first discovered in China in 2018.

METHODS AND RESULTS

Here, the complete mitochondrial genome (mitogenome) of T. sintenisi was sequenced and compared with the mitogenomes of other Diptera species. The results revealed that the mitogenome of T. sintenisi is 15,763 bp in length (GenBank accession number: MT795181) and is comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA genes, and a non-coding A + T-rich region (959 bp). The highly conserved gene arrangement of the mitogenome of T. sintenisi was identical to that of other Diptera insects. Twelve PCGs contained the typical insect start codon ATN, while cox1 had CGA as the start codon. The genes cox2, nad4, and nad1 contained an incomplete termination codon T; nad3, nad5, and cob contained the complete termination codon TAG; and the remaining seven PCGs contained the termination codon TAA. All tRNA genes were predicted to fold into the typical cloverleaf secondary structure. Phylogenetic analysis of 48 species based on the mitogenome sequence revealed that T. sintenisi clustered with the Tephritidae family, indicating that T. sintenisi and Tephritidae have a close phylogenetic relationship.

CONCLUSIONS

The phylogenetic relationship of T. sintenisi based on the mitogenome is consistent with the traditional morphological taxonomy, according to which T. sintenisi belongs to the family Otitidae, which is closely related to the family Muscidae.

Characterizing the Complete Mitochondrial Genomes of Three Bugs (Hemiptera: Heteroptera) Harming Bamboo

Herein, we report the mitochondrial genomic characteristics of three insect pests, Notobitus meleagris, Macropes harringtonae, and Homoeocerus bipunctatus, collected from bamboo plants in Guizhou Province, China. For the first time, the damaged conditions and life histories of M. harringtonae and H. bipunctatus are described in detail and digital photographs of all their life stages are provided. Simultaneously, the mitochondrial genome sequences of three bamboo pests were sequenced and analyzed. Idiocerus laurifoliae and Nilaparvata lugens were used as outgroups, and the phylogenetic trees were constructed. The mitochondrial genomes of the three bamboo pests contained 37 classical genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNAs (tRNAs), and a control region, with a total length of 16,199 bp, 15,314 bp, and 16,706 bp, respectively. The A+T values of the three bamboo pests were similar, and trnS1 was a cloverleaf structure with missing arms. The phylogenetic analyses, using the Bayesian inference (BI) and Maximum likelihood (ML), supported that N. meleagris and H. bipunctatus belonged to the Coreoidea family, whereas M. harringtonae belonged to the Lygaeoidea family with high support values. This study involves the first complete sequencing of the mitochondrial genomes of two bamboo pests. By adding these newly sequenced mitochondrial genome data and detailed descriptions of life histories, the database of bamboo pests is improved. These data also provide information for the development of bamboo pest control methods by quick identification techniques and the use of detailed photographs.

Full Mitochondrial Genomes Reveal Species Differences between the Venerid Clams Ruditapes philippinarum and R. variegatus

In natural sea areas along the coast of China, venerid clams Ruditapes philippinarum and R. variegatus exhibit similar adult shell forms and are especially difficult to distinguish as spat and juveniles. This study used comparative mitochondrial genome analysis to reveal differences between these species. The results showed that: (1) the mitochondrial genomes of R. philippinarum and R. variegatus share a large number of similar gene clusters arranged in consistent order, yet they also display noncommon genes, with both gene rearrangements and random losses found; (2) the 13 protein-coding genes in R. philippinarum as well as two-fold and four-fold degenerate sites in R. variegatus have an evident AT bias; (3) the Ka/Ks ratio of the mitochondrial ATP8 gene was significantly higher in R. philippinarum than in R. variegatus, and an analysis of selection pressure revealed that the mitochondrial NADH dehydrogenase subunit 2 gene and NADH dehydrogenase subunit 6 gene of R. variegatus were under great selective pressure during its evolution; and finally, (4) the two species clustered into one branch on a phylogenetic tree, further affirming their phylogenetic closeness. Based on these results, we speculate that the species differences between R. variegatus and R. philippinarum are largely attributable to adaptive evolution to the environment. The present findings provide a reference for the development of germplasm identification.

Mitogenomics of Chinch Bugs from China and Implications for Its Coevolutionary Relationship with Grasses

Blissidae (the Chinch bug) is a group with high species richness in Lygaeoidea, but there are only a few descriptions of mitochondrial genomes available. We obtained mitogenomes from 10 species of eight genera from Blissidae through second-generation sequencing technology. The length of the mitochondrial genome (excluding the control region) is between 14643 and 14385 bp; the content of AT is between 74.1% and 77.9%. The sequence of the evolution rate of protein coding genes was as follows: ND5 > ATP8 > ND6 > ND2 > ND4 > ND4L > ND1 > ATP6 > ND3 > COIII > COII > CYTB > COI. The mitogenomic structure of Blissidae is highly conservative. Gene rearrangement is only found in Pirkimeru japonicus (PiGXBS1), which is formed as the duplication of tRNA-H. The intergenic spacer between ND4 and tRNA-H, which form an obvious stem-and-loop structure, was found in all samples in this study. The phylogenetic trees generated by BI and ML indicated that Blissidae can be divided into three major clades: Clade A (only included Macropes); Clade B ((Pirkimerus + Bochrus) + Iphicrates); and Clade C ((Ischnodemus + Capodemus) + (Cavelerius + Dimorphopterus)). The divergence within the Blissidae began at about 56 Ma. At the genus level, the divergence was concentrated at 30−51 Ma, slightly later than the diversification of Poaceae. The consistency of divergence time between Blissidae and Poaceae might hint at the coevolutionary relationship between them, but further molecular and biological evidence is still needed to prove it.

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.

Comparative analysis of the complete mitochondrial genomes in two limpets from Lottiidae (Gastropoda: Patellogastropoda): rare irregular gene rearrangement within Gastropoda

To improve the systematics and taxonomy of Patellogastropoda within the evolution of gastropods, we determined the complete mitochondrial genome sequences of Lottia goshimai and Nipponacmea fuscoviridis in the family Lottiidae, which presented sizes of 18,192 bp and 18,720 bp, respectively. In addition to 37 common genes among metazoa, we observed duplication of the trnM gene in L. goshimai and the trnM and trnW genes in N. fuscoviridis. The highest A + T contents of the two species were found within protein-coding genes (59.95% and 54.55%), followed by rRNAs (56.50% and 52.44%) and tRNAs (56.42% and 52.41%). trnS1 and trnS2 could not form the canonical cloverleaf secondary structure due to the lack of a dihydrouracil arm in both species. The gene arrangements in all Patellogastropoda compared with those of ancestral gastropods showed different levels of gene rearrangement, including the shuffling, translocation and inversion of single genes or gene fragments. This kind of irregular rearrangement is particularly obvious in the Lottiidae family. The results of phylogenetic and gene rearrangement analyses showed that L. goshimai and Lottia digitalis clustered into one group, which in turn clustered with N. fuscoviridis in Patellogastropoda. This study demonstrates the significance of complete mitogenomes for phylogenetic analysis and enhances our understanding of the evolution of Patellogastropoda.

The complete mitochondrial genome of the Columbia lance nematode, Hoplolaimus columbus, a major agricultural pathogen in North America

BACKGROUND

The plant-parasitic nematode Hoplolaimus columbus is a pathogen that uses a wide range of hosts and causes substantial yield loss in agricultural fields in North America. This study describes, for the first time, the complete mitochondrial genome of H. columbus from South Carolina, USA.

METHODS

The mitogenome of H. columbus was assembled from Illumina 300 bp pair-end reads. It was annotated and compared to other published mitogenomes of plant-parasitic nematodes in the superfamily Tylenchoidea. The phylogenetic relationships between H. columbus and other 6 genera of plant-parasitic nematodes were examined using protein-coding genes (PCGs).

RESULTS

The mitogenome of H. columbus is a circular AT-rich DNA molecule 25,228 bp in length. The annotation result comprises 12 PCGs, 2 ribosomal RNA genes, and 19 transfer RNA genes. No atp8 gene was found in the mitogenome of H. columbus but long non-coding regions were observed in agreement to that reported for other plant-parasitic nematodes. The mitogenomic phylogeny of plant-parasitic nematodes in the superfamily Tylenchoidea agreed with previous molecular phylogenies. Mitochondrial gene synteny in H. columbus was unique but similar to that reported for other closely related species.

CONCLUSIONS

The mitogenome of H. columbus is unique within the superfamily Tylenchoidea but exhibits similarities in both gene content and synteny to other closely related nematodes. Among others, this new resource will facilitate population genomic studies in lance nematodes from North America and beyond.

The complete mitochondrial genome of the eusocial sponge-dwelling snapping shrimp Synalpheus microneptunus

In the marine realm, eusociality is only known to have evolved within a clade of sponge-dwelling snapping shrimps in the genus Synalpheus. Deciphering the genomic underpinnings of eusociality in these marine shrimps has been limited by the sparse genomic resources in this genus. Here, we report, for a eusocial shrimp Synalpheus microneptunus, a complete mitochondrial genome (22X coverage) assembled from short Illumina 150 bp paired-end reads. The 15,603 bp long mitochondrial genome of S. microneptunus is AT-rich and includes 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and an 834 bp intergenic region assumed to be the D-loop. The gene order is identical to that reported for most caridean shrimps and corresponds to the presumed Pancrustacean ground pattern. All PCGs showed signs of purifying selection, with K_A/K_S <<1 across the whole PCGs and most sliding windows within PCGs. Maximum-likelihood and Bayesian inference phylogenetic analyses of 13 PCGs and 68 terminals supports the monophyly of the Caridea and the family Alpheidae. The complete mitochondrial genome of the eusocial shrimp Synalpheus microneptunus will contribute to a better understanding of the selective pressures and rates of molecular evolution in marine eusocial animals.

Rearrangement and evolution of mitochondrial genomes in Thysanoptera (Insecta)

Prior to this study, complete mitochondrial genomes from Order Thysanoptera were restricted to a single family, the Thripidae, resulting in a biased view of their evolution. Here we present the sequences for the mitochondrial genomes of four additional thrips species, adding three extra families and an additional subfamily, thus greatly improving taxonomic coverage. Thrips mitochondrial genomes are marked by high rates of gene rearrangement, duplications of the control region and tRNA mutations. Derived features of mitochondrial tRNAs in thrips include gene duplications, anticodon mutations, loss of secondary structures and high gene translocation rates. Duplicated control regions are found in the Aeolothripidae and the 'core' Thripinae clade but do not appear to promote gene rearrangement as previously proposed. Phylogenetic analysis of thrips mitochondrial sequence data supports the monophyly of two suborders, a sister-group relationship between Stenurothripidae and Thripidae, and suggests a novel set of relationships between thripid genera. Ancestral state reconstructions indicate that genome rearrangements are common, with just eight gene blocks conserved between any thrips species and the ancestral insect mitochondrial genome. Conversely, 71 derived rearrangements are shared between at least two species, and 24 of these are unambiguous synapomorphies for clades identified by phylogenetic analysis. While the reconstructed sequence of genome rearrangements among the protein-coding and ribosomal RNA genes could be inferred across the phylogeny, direct inference of phylogeny from rearrangement data in MLGO resulted in a highly discordant set of relationships inconsistent with both sequence-based phylogenies and previous morphological analysis. Given the demonstrated rates of genomic evolution within thrips, extensive sampling is needed to fully understand these phenomena across the order.

The Complete Mitogenome of Pyrrhocoris tibialis (Hemiptera: Pyrrhocoridae) and Phylogenetic Implications

We determined the complete mitogenome of Pyrrhocoris tibialis (Hemiptera: Heteroptera: Pyrrhocoridae) to better understand the diversity and phylogeny within Pentatomomorpha, which is the second largest infra-order of Heteroptera. Gene content, gene arrangement, nucleotide composition, codon usage, ribosomal RNA (rRNA) structures, and sequences of the mitochondrial transcription termination factor were well conserved in Pyrrhocoroidea. Different protein-coding genes have been subject to different evolutionary rates correlated with the G + C content. The size of control regions (CRs) was highly variable among mitogenomes of three sequenced Pyrrhocoroidea species, with the P. tibialis CR being the largest. All the transfer RNA genes found in Pyrrhocoroidea had the typical clover leaf secondary structure, except for trnS1 (AGN), which lacked the dihydrouridine arm and possessed an unusual anticodon stem (9 bp vs. the normal 5 bp). A total of three different phylogenetic relationships among the five super-families of Pentatomomorpha were obtained using three analytical methods (MrBayes and RAxML under site-homogeneous models and PhyloBayes under a site-heterogeneous CAT + GTR model) and two mitogenomic datasets (nucleotides and amino acids). The tree topology test using seven methods statistically supported a phylogeny of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))) as the best topology, as recognized by both RAxML and MrBayes based on the two datasets.

The complete mitochondrial genome and description of a new cryptic species of Benedenia Diesing, 1858 (Monogenea: Capsalidae), a major pathogen infecting the yellowtail kingfish Seriola lalandi Valenciennes in the South-East Pacific

Background

The monogenean Benedenia seriolae parasitizes fishes belonging to the genus Seriola, represents a species complex, and causes substantial impact on fish welfare in aquaculture systems worldwide. This study reports, for the first time, the complete mitochondrial genome of B. humboldti n. sp., a new cryptic species from the South-East Pacific (SEP).

Methods

The mitogenome of B. humboldti n. sp. was assembled from short Illumina 150 bp pair-end reads. The phylogenetic position of B. humboldti n. sp. among other closely related congeneric and confamiliar capsalids was examined using mitochondrial protein-coding genes (PCGs). Morphology of B. humboldti n. sp. was examined based on fixed and stained specimens.

Results

The AT-rich mitochondrial genome of B. humboldti is 13,455 bp in length and comprises 12 PCGs (atp8 was absent as in other monogenean genomes), 2 ribosomal RNA genes, and 22 transfer RNA genes. All protein-coding, ribosomal RNA, and transfer RNA genes are encoded on the H-strand. The gene order observed in the mitochondrial genome of B. humboldti n. sp. was identical to that of B. seriolae from Japan but different from that of B. seriolae from Australia. The genetic distance between B. humboldti n. sp. and B. seriolae from Japan was high. Minor but reliable differences in the shape of the penis were observed between Benedenia humboldti n. sp. and congeneric species.

Conclusions

Phylogenetic analyses based on PCGs in association with differences in the shape of the penis permitted us to conclude that the material from the South-East Pacific represents a new species of Benedenia infecting S. lalandi off the coast of Chile. The discovery of this parasite represents the first step to improving our understanding of infestation dynamics and to develop control strategies for this pathogen infecting the farmed yellowtail kingfish, Seriola lalandi, in the South-East Pacific.

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.

Characterization of the complete mitochondrial genomes of two species of the genus Aphaena Guérin-Méneville (Hemiptera: Fulgoridae) and its phylogenetic implications

The complete mitochondrial genomes (mitogenomes) of Aphaena (Callidepsa) amabilis and Aphaena (Aphaena) discolor nigrotibiata were sequenced. The mitogenomes of these two species are 16,237 bp and 16,116 bp in length with an A + T content of 77.9% and 77.0%, respectively. Each contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (A + T-rich region). All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA or TAG except for atp6, where nad1 ends with an incomplete T codon. All tRNAs have the typical clover-leaf structure except for trnS1 and trnV which have a reduced DHU arm. Moreover, these two mitogenomes have trnL2, trnR and trnT with an unpaired base in the acceptor stem. The putative A + T-rich region includes multiple types of tandem repeat regions. These phylogenetic analyses are reconstructed based on 13 protein-coding genes of 25 auchenorrhynchan mitogenomes, with both maximum likelihood and Bayesian analyses yielding robust identical phylogenetic trees. These results support a monophyletic Auchenorrhyncha and the relationship (Pyrops + (Lycorma + Aphaena)) within Fulgoridae.

Comparative analysis of the complete mitochondrial genomes of five Achilidae species (Hemiptera: Fulgoroidea) and other Fulgoroidea reveals conserved mitochondrial genome organization

In the present study, the complete mitochondrial genomes (mitogenomes) of five Achilidae (Hemiptera: Fulgoroidea), Betatropis formosana, two new species (Magadhaideus luodiana sp. nov and Peltatavertexalis horizontalis sp. nov), Plectoderini sp. and Paracatonidia sp., were sequenced for the first time through next-generation sequencing. The five mitogenomes ranged from 15,214 to 16,216 bp in length, with the typical gene content and arrangement usually observed in Hexapods. The motif "ATGATAA" between atp8 and atp6 was found in all the analyzed species. An overlap "AAGCTTA" between trnW and trnC was observed in the mitogenomes of most Fulgoroidea. The structural and compositional analyses of 26 Fulgoroidea mitogenomes, including the gene rearrangement of five tRNAs (trnW, trnC and trnY; trnT and trnP), the A + T content and AT-skew of the whole mitogenomes, and the nuclear acid and amino acid compositions of the protein-coding genes (PCGs), revealed family-level differences between Delphacidae and other families (Achilidae, Flatidae, Fulgoridae, Issidae and Ricaniidae). Phylogenetic analyses of 13 protein-coding genes from 26 Fulgoroidea species by maximum likelihood and Bayesian Inference were consistent and well supported the basal position of Delphacidae, a close affinity among the families Flatidae, Issidae and Ricaniidae, and a close relationship between Achilidae and Fulgoridae.

The complete mitochondrial genome of Tetraphleps aterrimus (Hemiptera: Anthocoridae): Genomic comparisons and phylogenetic analysis of Cimicomorpha

The mitochondrial genome (mitogenome) provides important information for better understanding the phylogenetic relationships within heteropteran infraorder Cimicomorpha (Hemiptera: Heteroptera), but there are still limited representations at the family level of Anthocoridae. Here we sequenced the complete mitogenome of Tetraphleps aterrimus. It is 15,803 bp in size, and contains the expected 37 genes (13 PCGs, 22 tRNAs and 2 rRNAs) and control region. Gene order is identical to that of typical cimicomorphans. In comparison with other cimicomorphans, the ratios of Ka/Ks are increasing from 0.17 for COI to 0.85 for ATP8, which demonstrates COI shows the lowest evolutionary rate, while ATP8 appears to be the highest. The ratios of conserved sites of COI is the highest, while ATP8 is the lowest, suggesting that the evolutionary rate of ATP8 is higher than COI. The phylogenetic relationships based on mitogenomes using Bayesian inference (BI) and Maximum likelihood (ML) methods show that Tetraphleps aterrimus is sister to (Orius niger + Orius sauteri), suggesting that Tetraphleps aterrimus belongs to Anthocoridae. The monophyly of each superfamily is generally well supported and Reduvioidea is placed as basal branch in Cimicomorpha. The results support the remaining superfamily groupings (Miroidea + (Cimicoidea + (Velocipedoidea + Nabioidea))).

Comparative analysis of the mitochondrial genomes of oriental spittlebug trible Cosmoscartini: insights into the relationships among closely related taxa

Background

Cosmoscartini (Hemiptera: Cercopoidea: Cercopidae) is a large and brightly colored Old World tropical tribe, currently containing over 310 phytophagous species (including some economically important pests of eucalyptus in China) in approximately 17 genera. However, very limited information of Cosmoscartini is available except for some scattered taxonomic studies. Even less is known about its phylogenetic relationship, especially among closely related genera or species. In this study, the detailed comparative genomic and phylogenetic analyses were performed on nine newly sequenced mitochondrial genomes (mitogenomes) of Cosmoscartini, with the purpose of exploring the taxonomic status of the previously defined genus Okiscarta and some closely related species within the genus Cosmoscarta.

Results

Mitogenomes of Cosmoscartini display similar genomic characters in terms of gene arrangement, nucleotide composition, codon usage and overlapping regions. However, there are also many differences in intergenic spacers, mismatches of tRNAs, and the control region. Additionally, the secondary structures of rRNAs within Cercopidae are inferred for the first time.

Based on comparative genomic (especially for the substitution pattern of tRNA secondary structure) and phylogenetic analyses, the representative species of Okiscarta uchidae possesses similar structures with other Cosmoscarta species and is placed consistently in Cosmoscarta. Although Cosmoscarta bimacula is difficult to be distinguished from Cosmoscarta bispecularis by traditional morphological methods, evidence from mitogenomes highly support the relationships of (C. bimacula + Cosmoscarta rubroscutellata) + (C. bispecularis + Cosmoscarta sp.).

Conclusions

This study presents mitogenomes of nine Cosmoscartini species and represents the first detailed comparative genomic and phylogenetic analyses within Cercopidae. It is indicated that knowledge of mitogenomes can be effectively used to resolve phylogenetic relationships at low taxonomic levels. Sequencing more mitogenomes at various taxonomic levels will also improve our understanding of mitogenomic evolution and phylogeny in Cercopidae.

Electronic supplementary material

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

Comparative analysis of seven mitochondrial genomes of Phymatostetha (Hemiptera: Cercopidae) and phylogenetic implications

In this study, we present seven mitochondrial genomes (mitogenomes) of Phymatostetha. Each mitogenome contains the entire set of 37 genes, which arranged in the same order as the putative ancestral pattern of insects. The nucleotide composition of Phymatostetha mitogenomes is biased toward A/T, with rRNAs and PCG12 (i.e. the first and second codon positions of PCGs) exhibit the highest and lowest A + T content, respectively. Relative synonymous codon usage of PCGs also show that degenerate codons are biased to use more A/T than G/C. All tRNAs exhibit typical clover-leaf structure, with the exception of trnS1. Additionally, unpaired nucleotides are detected in trnS1 anticodon stem and trnR acceptor stem. Phylogenetic relationships, based on the dataset of 13 PCGs, 22tRNAs, and two rRNAs, are analyzed using both the Bayesian and maximum likelihood methods. Our results clearly revealed the systematic status of Phymatostetha species and robustly supported the monophyly of this genus, in which Phymatostetha semele is sister to other Phymatostetha species. It was demonstrated that mitogenome was an effective molecular marker to adequately resolve phylogeny at low taxonomic levels.

The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus

Panulirus argus is a keystone species and target of the most lucrative fishery in the Caribbean region. This study reports, for the first time, the complete mitochondrial genome of Panulirus argus (average coverage depth nucleotide-1 = 70×) assembled from short Illumina 150 bp PE reads. The AT-rich mitochondrial genome of Panulirus argus was 15 739 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 801 bp long intergenic space was assumed to be the D-loop. Most of the PCGs were encoded on the H-strand. The gene order observed in the mitochondrial genome of Panulirus argus corresponds to the presumed Pancrustacean ground pattern. KA/KS ratios calculated for all mitochondrial PCGs showed values < 1, indicating that all these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 154 arthropods) supported the monophyly of the Achelata and other infraorders within the Decapoda. Mitochondrial PCGs have enough phylogenetic informativeness to explore high-level genealogical relationships in the Pancrustacea. The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus will contribute to the better understanding of meta-population connectivity in this keystone overexploited species.

Comparative Characterization of the Complete Mitochondrial Genomes of the Three Apple Snails (Gastropoda: Ampullariidae) and the Phylogenetic Analyses

The apple snails Pomacea canaliculata, Pomacea diffusa and Pomacea maculate (Gastropoda: Caenogastropoda: Ampullariidae) are invasive pests causing massive economic losses and ecological damage. We sequenced and characterized the complete mitochondrial genomes of these snails to conduct phylogenetic analyses based on comparisons with the mitochondrial protein coding sequences of 47 Caenogastropoda species. The gene arrangements, distribution and content were canonically identical and consistent with typical Mollusca except for the tRNA-Gln absent in P. diffusa. An identifiable control region (d-loop) was absent. Bayesian phylogenetic analysis indicated that all the Ampullariidae species clustered on the same branch. The genus Pomacea clustered together and then with the genus Marisa. The orders Architaenioglossa and Sorbeoconcha clustered together and then with the order Hypsogastropoda. Furthermore, the intergenic and interspecific taxonomic positions were defined. Unexpectedly, Ceraesignum maximum, Dendropoma gregarium, Eualetes tulipa and Thylacodes squamigerus, traditionally classified in order Hypsogastropoda, were isolated from the order Hypsogastropoda in the most external branch of the Bayesian inference tree. The divergence times of the Caenogastropoda indicated that their evolutionary process covered four geological epochs that included the Quaternary, Neogene, Paleogene and Cretaceous periods. This study will facilitate further investigation of species identification to aid in the implementation of effective management and control strategies of these invasive species.

Characterization of the Complete Mitochondrial Genome Sequences of Three Croakers (Perciformes, Sciaenidae) and Novel Insights into the Phylogenetics

The three croakers (Nibea coibor, Protonibea diacanthus and Argyrosomus amoyensis, Perciformes, Sciaenidae) are important commercial species inhabiting the Eastern Indian Ocean and Western Pacific. Molecular data employed in previous research on phylogenetic reconstruction have not been adequate and complete, and systematic and comprehensive phylogenetic relationships for these fish are unresolved. We sequenced the complete mitochondrial genomes of the three croakers using next-generation sequencing for the first time. We analyzed the composition and phylogenies between 19 species in the family Sciaenidae using the mitochondrial protein coding sequences of 204 species in the Series Eupercaria. We present the characterization of the complete mitochondrial genome sequences of the three croakers. Gene arrangement and distribution of the three croakers are canonically identical and consistent with other vertebrates. We found that the family Sciaenidae is an independent branch that is isolated from the order Perciformes and does not belong to any extant classification. Therefore, this family is expected to belong to a new classification at the order level and needs further analysis. The evolution of Sciaenidae has lagged far behind the Perciformes differentiation. This study presents a novel insight into the phylogenetics of the family Sciaenidae from the order Perciformes and facilitates additional studies on the evolution and phylogeny of Series Eupercaria.

The Complete Mitochondrial Genome of Ugyops sp. (Hemiptera: Delphacidae)

The complete mitochondrial genome (mitogenome) of Ugyops sp. (Hemiptera: Delphacidae) was sequenced, making it the first determined mitogenome from the subfamily Asiracinae, the basal clade of the family Delphacidae. The mitogenome was 15,259 bp in length with A + T content of 77.65% and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. The gene order was identical with that of the ancestral insect. The nucleotide composition analysis indicated that the whole mitogenome was strongly A-skewed (0.288) and highly C-skewed (-0.270). For PCGs on the J-strand, the AT skew was positive, and the GC skew was negative. All PCGs started with canonical ATN codons, except for cox1 and nad5, which used CTG and GTG as start codon, respectively. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnS1 (AGN), in which the dihydrouridine arm was reduced to a simple loop. The control region included a poly-T stretch downstream of the small rRNA gene (rrnS), a subregion of higher A + T content and tandemly repeated sequence near trnI. The mitogenome of Ugyops sp. could be very helpful in exploring the diversity and evolution of mitogenomes in Delphacidae.

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.

Vectors of diversity: Genome wide diversity across the geographic range of the Chagas disease vector Triatoma dimidiata sensu lato (Hemiptera: Reduviidae)

To date, the phylogeny of Triatoma dimidiata sensu lato (s. l.) (Hemiptera: Reduviidae: Triatominae), the epidemiologically most important Chagas disease vector in Central America and a secondary vector in Mexico and northern South America, has only been investigated by one multi-copy nuclear gene (Internal Transcribed Spacer - 2) and a few mitochondrial genes. We examined 450 specimens sampled across most of its native range from Mexico to Ecuador using reduced representation next-generation sequencing encompassing over 16,000 single nucleotide polymorphisms (SNPs). Using a combined phylogenetic and species delimitation approach we uncovered two distinct species, as well as a well-defined third group that may contain multiple species. The findings are discussed with respect to possible drivers of diversification and the epidemiological importance of the distinct species and groups.

Compositional heterogeneity in true bug mitochondrial phylogenomics

Mitochondrial phylogenomics is often controversial, in particular for inferring deep relationships. The recent rapid increase of mitochondrial genome data provides opportunities for better phylogenetic estimates and assessment of potential biases resulting from heterogeneity in nucleotide composition and mutation rates. Here, we gathered 76 mitochondrial genome sequences for Heteroptera representing all seven infraorders, including 17 newly sequenced mitochondrial genomes. We found strong heterogeneity in base composition and contrasting evolutionary rates among heteropteran mitochondrial genomes, which affected analyses with various datasets and partitioning schemes under site-homogeneous models and produced false groupings of unrelated taxa exhibiting similar base composition and accelerated evolutionary rates. Bayesian analyses using a site-heterogeneous mixture CAT+GTR model showed high congruence of topologies with the currently accepted phylogeny of Heteroptera. The results confirm the monophyly of the six infraorders within Heteroptera, except for Cimicomorpha which was recovered as two paraphyletic clades. The monophyly of Terheteroptera (Cimicomorpha and Pentatomomorpha) and Panheteroptera (Nepomorpha, Leptopodomorpha and Terheteroptera) was recovered demonstrating a significant improvement over previous studies using mitochondrial genome data. Our study shows the power of the site-heterogeneous mixture models for resolving phylogenetic relationships with Heteroptera and provides one more case showing that model adequacy is critical for accurate tree reconstruction in mitochondrial phylogenomics.