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

Phylogenetic and divergence analysis of Pentatomidae, with a comparison of the mitochondrial genomes of two related species (Hemiptera, Pentatomidae)

Pentatomidae, the most diverse family of Pentatomoidea, is found worldwide. Currently, the phylogenetic relationships among Pentatomidae tribes remain unstable, and subfamily divergence has not been estimated. Here, we sequenced and analyzed the complete mitochondrial genomes of two species of Lelia, and studied the phylogenetic relationships among Pentatominae tribes. We also selected three available fossil as the calibration points in the family, and preliminarily discussed the divergence time of Pentatomidae. Trees of Pentatomidae were reconstructed using the Bayesian inference method. Divergence times of Pentatominae were estimated based on the nucleotide sequences of protein-coding genes with a relaxed clock log-normal model in BEASTv.1.8.2. The results showed that the gene arrangements, nucleotide composition, and codon preferences were highly conserved in Lelia. Further, a phylogenetic analysis recovered Eysarcorini, Strachiini, Phyllocephalini, and Menidini as monophyletic with strong support, however, the monophyly of Antestiini, Nezarini, Carpocorini, Pentatomini and Cappaeini were rejected. Moreover, Pentatominae diverged from Pentatomidae soon after the origin of the Cretaceous Period, at approximately 110.38 Ma. This study enriches the mitochondrial genome database of Pentatomidae and provides a reference for further phylogenetic studies, and provides a more accurate estimate of divergence time.

The complete mitochondrial genome of Taiwanaptera montana (Hemiptera: Aradidae)

Taiwanaptera montana is a apterous flat bug found in Yunnan Province, China. This study is the first to sequence, assemble, and annotate the complete mitochondrial genome of T. montana. The mitochondrial genome of T. montana has a total length of 15615 bp, which is a typical circular DNA, containing 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and a control region, with A + T content of 69.4%. The phenomenon of gene rearrangement was found when the mitochondrial structure was compared with that of other Aradidae. The phylogenetic tree based on 37 mitochondrial genes showed that T. montana was most closely related to Libiocoris heissi. Aneurinae, Carventinae and Mezirinae are monophyletic groups. In addition, the results also confirmed that Aradinae and Calisiinae diverged earliest and were relatively primitive groups.

Mitochondrial Genome Characteristics and Phylogenetic Analysis of Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae)

Sugarcane thrips, Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae), is a common foliar pest that infests sugarcane and is found throughout tropical and subtropical countries. In this study, we obtained and analyzed the complete mitochondrial genome of F. serrata for the first time and explored the phylogenetic relationships of the higher-order elements of Thysanoptera members at the mitochondrial level. The complete mitochondrial genome of F. serrata is 16,596 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 1 noncoding control region. A+T accounted for 75% of the total bases in the mitochondrial genome of F. serrata, revealing an obvious AT bias. Among the 13 PCGs, except for nad5, which had a start codon of TTG, the remaining genes had ATNs typical of insects (ATA, ATT, ATC, and ATG); nad1, nad2, nad3, and atp8 had incomplete termination codons of TA or T. The remaining nine PCGs were complete with the termination codon TAA. Of the 22 tRNA secondary structures, all were typical cloverleaf secondary structures except for trnS1, which was missing the DHU arm. Compared with the hypothetical ancestral gene arrangement of arthropods, F. serrata presented extensive gene rearrangement, with 23 translocated genes, 8 inverted genes, and 5 shuffled genes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies: ((Thripidae + (Stenurothripidae + Aeolothripidae)) + Phlaeothripidae), with Thripidae, Aeolothripidae and Phlaeothripidae being monophyletic groups, whereas F. serrata is closely related to Thrips palmi, and the two are sister groups.

Differences between phytophagous and predatory species in Pentatomidae based on the mitochondrial genome

Pentatomidae includes many species of significant economic value as plant pests and biological control agents. The feeding habits of Pentatomidae are closely related to their energy metabolism and ecological adaptations. In this study, we sequenced the mitochondrial genomes of 12 Asopinae species using the next-generation sequencing to explore the effect of dietary changes on mitochondrial genome evolution. Notably, all sequences were double-stranded circular DNA molecules containing 37 genes and one control region. We then compared and analyzed the mitochondrial genome characteristics of phytophagous and predatory bugs. Notably, no significant difference was observed in the length of the mitochondrial genomes between the predatory and phytophagous bugs. However, the AT content was higher in the mitochondrial genomes of phytophagous bugs than that of predatory bugs. Moreover, phytophagous bugs prefer codon usage patterns ending in A/T compared with predatory bugs. The evolution rate of predatory bugs was lower than that of phytophagous bugs. The phylogenetic relationships across phytophagous bugs' lineages were largely consistent at depth nodes based on different datasets and tree-reconstructing methods, and strongly supported the monophyly of predatory bugs. Additionally, the estimated divergence times indicated that Pentatomidae explosively radiated in the Early Cretaceous. Subsequently, the subfamily Asopinae and the genus Menida diverged in the Late Cretaceous. Our research results provide data supporting for the evolutionary patterns and classification of Pentatomidae.

Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae

The mitochondrial genomes of three species of Yangiella were sequenced, annotated, and analyzed. The genome length of the three species of the genus is 15,070-15,202 bp, with a typical gene number, including a control region, 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). It was found that the mitochondrial genome of Yangiella had AT bias. Except for the lack of a DHU arm of the trnS1 gene, the other tRNAs had a typical cloverleaf structure, and the codon usage preferences of the three species exhibited high similarity. In addition, tRNA gene rearrangements were observed among the three subfamilies of Aradidae (Mezirinae, Calisiinae, Aradinae), and it was found that codon usage preferences appeared to be less affected by base mutation and more by natural selection. The Pi and Ka/Ks values indicated that cox1 was the most conserved gene in the mitochondrial genome of Aradidae, while atp8 and nad6 were rapidly evolved genes. Substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, suggesting the rationality of the phylogenetic analysis data. Bayesian and maximum likelihood methods were used to analyze the phylogeny of 16 species of Hemiptera insects, which supported the monophyly of Aneurinae, Carventinae, and Mezirinae, as well as the monophyly of Yangiella. Based on fossils and previous studies, the differentiation time was inferred, indicating that Yangiella diverged about 57 million years ago.

The complete mitochondrial genome of Neuroctenus taiwanicus (Hemiptera: Aradidae)

Neuroctenus taiwanicus Kormilev, 1955 is a flat-bodied and enigmatic bug that was first discovered on the island of Taiwan, China. In this study, the whole mitochondrial genome of N. taiwanicus was sequenced and annotated for the first time, and its genomic data were uploaded to Genbank feedback number OR675057. The mitochondrial genome of N. taiwanicus is 15,340 bp in length, a typical circular DNA encoding 37 genes and a control region with 68.4% A + T content. The phylogeny reveals the taxonomic status of N. taiwanicus, which is most closely related to N. yunnanensis, and demonstrates the sister relationship among Neuroctenus, Mezira, and Brachyrhynchus. In addition, the results also confirm that Aradinae and Calisiinae are the earliest branching and more primitive in the family Aradidae, which is consistent with the analysis of the traditional classification.

Describe the morphology and mitochondrial genome of Mecidea indica Dallas, 1851 (Hemiptera, Pentatomidae), with its phylogenetic position

We here describe the external morphology and complete mitochondrial genome characteristics of Mecidea indica Dallas, 1851, and clarify the evolutionary rate and divergence time. The M. indica mitochondrial genome length is 15,670 bp, and it exhibits a typical high A+T-skew (76.31%). The sequence shows strong synteny with the original gene arrangement of Drosophila yakuba Burla, 1954 without rearrangement. The M. indica mitochondrial genome characteristics were analyzed, and phylogenetic trees of Pentatomidae were reconstructed using Bayesian methods based on different datasets of the mitochondrial genome datasets. Phylogenetic analysis shows that M. indica belongs to Pentaotominae and form a sister-group with Anaxilaus musgravei Gross, 1976, and Asopinae is highly supported as monophyletic. Molecular clock analysis estimates a divergence time of Pentatomidae of 122.75 Mya (95% HPD: 98.76-145.43 Mya), within the Mesozoic Cretaceous; the divergence time of M. indica and A. musgravii was no later than 50.50 Mya (95% HPD: 37.20-64.80 Mya). In addition, the divergence time of Asopinae was 62.32 Mya (95% HPD: 47.08-78.23 Mya), which was in the Paleogene of the Cenozoic era. This study is of great significance for reconstructing the phylogeny of Pentatomidae and providing insights into its evolutionary history.

Complete mitochondrial genome of Cyclopelta obscura (Lepeletier & Serville, 1825) (Hemiptera: Pentatomoidea: Dinidoridae) and phylogenetic analysis of Pentatomoidea species

Cyclopelta obscura is a crop pest, which mainly damages legumes, especially Robinia pseudoacacia and Cercis chinensis. In recent years, many mitochondrial and nuclear DNA sequences of C. obscura have been sequenced and used for phylogenetic inference. However, the complete mitogenome has not been reported yet and studies on the phylogenetic relationships within Dinidoridae are rare. In this study, we sequenced the mitogenome of C. obscura and conducted comparative mitogenomic analyses of seven Dinidoridae species based on several different factors. The length of the mitogenome is 15,426 bp, which includes 37 typical mitochondrial genes (13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs) and a control region (796 bp long), as well as 13 intergenic spacers and 8 overlapping regions. Most PCGs of C. obscura began with the classical start codon ATN, while cox1 and nad4l used TTG, and nad1 used GTG. The Dinidoridae mitogenomes are highly conserved in terms of nucleotide composition, the codon usage of PCGs, and the secondary structure of tRNA. Phylogenetic analysis based on four datasets with two methods recovered the Dinidoridae as a monophyletic group with strong support values. All results indicate that Dinidoridae formed a sister group to Tessaratomidae, and (Tessaratomidae + Dinidoridae) formed a sister group to Cydnidae in most of the phylogenetic trees. Additionally, seven species within the Dinidoridae, we observed the following relationship: (Eumenotes sp. + (Cyclopelta parva + C. obscura)) + ((Megymenum gracilicorne + Megymenum brevicorne) + (Coridius chinensis + Coridius brunneu)).

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

Scanning Electron Microscopy of Antennae and Mouthparts of Mezira yunnana Hsiao (Hemiptera: Aradidae): Specialized Microstructures Reflecting Adaptation to Mycetophagy

Many species of the family Aradidae (also known as flat bugs) feed on fungal mycelia and fruiting bodies. In order to better understand the morphological adaptation to this unique feeding habit, we examined the microstructure of antennae and mouthparts of an aradid species, Mezira yunnana Hsiao, using scanning electron microscope, and documented the fungal feeding process under laboratory conditions. The antennal sensilla include three subtypes of sensilla trichodea, three subtypes of sensilla basiconica, two subtypes of sensilla chaetica, sensilla campaniformia, and sensilla styloconica. The apex of the second segment of flagellum has a large number of various sensilla forming a sensilla cluster. The labial tip is distally constricted, which is rarely observed in other Pentatomomorpha species. The labial sensilla include three subtypes of sensilla trichodea, three subtypes of sensilla basiconica, and a sensilla campaniformia. The tip of the labium has only three pairs of sensilla basiconica III and small comb-shaped cuticular processes. The external surface of the mandibular apex has 8-10 ridge-like central teeth. A series of key morphological structures associated with mycetophagous feeding habit were identified, which will facilitate future studies on adaptive evolution of species in Pentatomomorpha as well as in other heteropteran lineages.

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.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

The complete mitochondrial genome of Brachyrhynchus triangulus (Hemiptera: Aradoidea: Aradidae)

The complete mitochondrial genome of Brachyrhynchus triangulus Bergroth, 1889 was sequenced and annotated in the present study. It was a typical circular DNA molecule of 15,170 bp, with 37 genes and a control region. The percentages of A, C, G, and T nucleotides in the genome of B. triangulus were 41.1%, 17.4%, 11.9%, and 29.4%, respectively. Thirteen protein-coding genes (PCGs) start with a ATN codon or a TTG codon and terminate with a TAA codon or a TAG codon or a single T residue. With the exception of tRNA^Ser(AGN) , each of the 22 tRNA genes had a clover-leaf structure and ranged in length from 62 to 69 bp. The length of lrRNA and srRNA was 1241 bp and 828 bp, respectively. The control region had a length of 708 bp and an A + T content of 74.6%. The sister relationship between B. triangulus and Brachyrhynchus hsiaoi is supported by the phylogenetic tree. Additionally, it proved the sister relationship between Mezirinae and Aneurinae, supporting the classical taxonomy of the Aradidae.

A Rearrangement of the Mitochondrial Genes of Centipedes (Arthropoda, Myriapoda) with a Phylogenetic Analysis

Due to the limitations of taxon sampling and differences in results from the available data, the phylogenetic relationships of the Myriapoda remain contentious. Therefore, we try to reconstruct and analyze the phylogenetic relationships within the Myriapoda by examining mitochondrial genomes (the mitogenome). In this study, typical circular mitogenomes of Mecistocephalus marmoratus and Scolopendra subspinipes were sequenced by Sanger sequencing; they were 15,279 bp and 14,637 bp in length, respectively, and a control region and 37 typical mitochondrial genes were annotated in the sequences. The results showed that all 13 PCGs started with ATN codons and ended with TAR codons or a single T; what is interesting is that the gene orders of M. marmoratus have been extensively rearranged compared with most Myriapoda. Thus, we propose a simple duplication/loss model to explain the extensively rearranged genes of M. marmoratus, hoping to provide insights into mitogenome rearrangement events in Myriapoda. In addition, our mitogenomic phylogenetic analyses showed that the main myriapod groups are monophyletic and supported the combination of the Pauropoda and Diplopoda to form the Dignatha. Within the Chilopoda, we suggest that Scutigeromorpha is a sister group to the Lithobiomorpha, Geophilomorpha, and Scolopendromorpha. We also identified a close relationship between the Lithobiomorpha and Geophilomorpha. The results also indicate that the mitogenome can be used as an effective mechanism to understand the phylogenetic relationships within Myriapoda.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

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

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

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

External structures of the metathoracic scent gland efferent system in the true bug superfamily Pyrrhocoroidea (Hemiptera: Heteroptera: Pentatomomorpha)

Pyrrhocoroidea represents an important group of true bugs (Insecta: Hemiptera: Heteroptera) which includes fire bugs, cotton stainers and other taxa widely used in experimental studies or known as pests. However, the morphology and phylogeny of Pyrrhocoroidea have been only poorly studied so far. Here, structures of the external scent efferent system of the metathoracic scent glands are examined in 64 out of 71 currently valid genera of Pyrrhocoroidea and scanning electron micrographs are provided for most taxa. Several characters are revealed which define each of the three higher taxa within Pyrrhocoroidea: Larginae (small auriculate peritreme lacking manubrium and median furrow; metathoracic spiracle never surrounded by evaporatorium), Physopeltinae (large, widely open ostiole; large peritremal disc with manubrium [new term], lacking median furrow; mace-like mycoid filter processes of equal shape and size on both anterior and posterior margins of metathoracic spiracle), and Pyrrhocoridae (elongate auriculate peritreme with deep median furrow). Within Pyrrhocoridae, three main types (A, B and C) of the external scent efferent system are distinguished, differring in the amount of reductions. The findings are interpreted in the context of phylogenetic hypotheses available for Pyrrhocoroidea and their close relatives, Coreoidea and Lygaeoidea. An updated identification key to the families and subfamilies of Pyrrhocoroidea applicable for both sexes is provided.

Complete mitochondrial genome sequence for the Thrips hawaiiensis (Thysanoptera: Thripidae)

Thrips hawaiiensis (Morgan) (Thysanoptera: Thripidae) is a common Thysanoptera insect widely distributed in Asia and the Pacific, it damages various plants. In this study the complete mitochondrial genome of T. hawaiiensis was sequenced and characterized by using next-generation sequencing technique. The total length of the complete genome is 15,357 bp and A + T content of 77.8% (GeneBank accession No. MW582621). The T. hawaiiensis mitochondrial genome consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes (tRNAs) and 2 non-coding control regions (Dloop region). According to previous studies, only a few complete mitochondrial genomes from Order Thysanoptera have been reported. Thus, T. hawaiiensis complete mitochondrial genome sequence reported will provide molecular information for mitochondrial genome research on Thysanoptera.

The complete mitochondrial genome of Eclipophleps carinata (Orthoptera: Acridoidea: Gomphoceridae)

The complete Eclipophleps carinata mitogenome was sequenced through Illumina HiSeq 2500 platform and the resulting data were analyzed in this paper. The mitochondrial genome of E. carinata is a typical circular DNA molecule of 15,781 bp with 37 genes and 74.5% A + T content, which encoded 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and the control region. The E. carinata mitochondrial genome and 27 mitochondrial genomes (downloaded from NCBI) were employed to construct phylogenetic tree, in which Ruspolia dubia and Teleogryllus occipitails were the outgroups. Phylogenetic reconstruction validated the taxonomic status of E. carinata, which was placed in the monophyletic Gomphocerinae in Acrididae.

The complete mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) based on long-read nanopore sequencing

BACKGROUND

The coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is a severe and invasive pest of coconut and other palms throughout Asia and the Pacific. The biocontrol agent, Oryctes rhinoceros nudivirus (OrNV), has successfully suppressed O. rhinoceros populations for decades but new CRB invasions started appearing after 2007. A single-SNP variant within the mitochondrial cox1 gene is used to distinguish the recently-invading CRB-G lineage from other haplotypes, but the lack of mitogenome sequence for this species hinders further development of a molecular toolset for biosecurity and management programmes against CRB. Here we report the complete circular sequence and annotation for CRB mitogenome, generated to support such efforts.

METHODS

Sequencing data were generated using long-read Nanopore technology from genomic DNA isolated from a CRB-G female. The mitogenome was assembled with Flye v.2.5, using the short-read Illumina sequences to remove homopolymers with Pilon, and annotated with MITOS. Independently-generated transcriptome data were used to assess the O. rhinoceros mitogenome annotation and transcription. The aligned sequences of 13 protein-coding genes (PCGs) (with degenerate third codon position) from O. rhinoceros, 13 other Scarabaeidae taxa and two outgroup taxa were used for the phylogenetic reconstruction with the Maximum likelihood (ML) approach in IQ-TREE and Bayesian (BI) approach in MrBayes.

RESULTS

The complete circular mitogenome of O. rhinoceros is 20,898 bp in length, with a gene content canonical for insects (13 PCGs, two rRNA genes, and 22 tRNA genes), as well as one structural variation (rearrangement of trnQ and trnI) and a long control region (6,204 bp). Transcription was detected across all 37 genes, and interestingly, within three domains in the control region. ML and BI phylogenies had the same topology, correctly grouping O. rhinoceros with one other Dynastinae taxon, and recovering the previously reported relationship among lineages in the Scarabaeidae. In silico PCR-RFLP analysis recovered the correct fragment set that is diagnostic for the CRB-G haplogroup. These results validate the high-quality of the O. rhinoceros mitogenome sequence and annotation.

Comparative Mitogenomic Analysis of Two Cuckoo Bees (Apoidea: Anthophila: Megachilidae) with Phylogenetic Implications

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.

The mitochondrial genome of the assassin bug Sclomina erinacea (Hemiptera: Reduviidae)

The complete mitochondrial genome (mitogenome) of the assassin bug, Sclomina erinacea Stål, was determined in the present study. The sequenced mitogenome is a typical circular DNA molecule which is 15,828 bp in length, containing 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a putative control region. All protein-coding genes are initiated by ATN codons and terminated by TAA or TAG codons except COII, COIII and ND4 which use a single T residue as the stop codon. All tRNAs have the cloverleaf structure with the exception of tRNA^Ser(AGN) and the length of them range from 62 to 70 bp. The control region is 818 bp long with an A + T content of 67.1%. The phylogenetic analysis supports the monophyly of Harpactorinae and Sclomina erinacea is the closest relative to Macracanthopsis nodipes.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

Characterization of the Complete Mitochondrial Genome of Ostertagia trifurcata of Small Ruminants and its Phylogenetic Associations for the Trichostrongyloidea Superfamily

The complete mitochondrial (mt) genome of Ostertagia trifurcata, a parasitic nematode of small ruminants, has been sequenced and its phylogenetic relationship with selected members from the superfamily Trichostrongyloidea was investigated on the basis of deduced datasets of mt amino acid sequences. The entire mt genome of Ostertagia trifurcata is circular and 14,151 bp in length. It consists of a total of 36 genes comprising 12 genes coding for proteins (PCGs), 2 genes for ribosomal RNA (rRNA), 22 transfer RNA (tRNA) genes and 2 non-coding regions, since all genes are transcribed in the same direction. The phylogenetic analysis based on the concatenated datasets of predicted amino acid sequences of the 12 protein coding genes supported monophylies of the Haemonchidae, Dictyocaulidae and Molineidae families, but rejected monophylies of the Trichostrongylidae family. The complete characterization and provision of the mtDNA sequence of Ostertagia trifurcata provides novel genetic markers for molecular epidemiological investigations, systematics, diagnostics and population genetics of Ostertagia trifurcata and its correspondents.

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.

The Complete Mitochondrial Genome of the Plant Bug Lygus pratensis Linnaeus (Hemiptera: Miridae)

Lygus pratensis is a phytophagous pest responsible for yield losses in Bt alfalfa and other economic crops in Northwestern China. To better characterize Miridae at the genomic level, the complete mitochondrial (mt) genome of L. pratensis was sequenced and analyzed in this study. The mt genome was amplified via the polymerase chain reaction to generate overlapping fragments. These fragments were then sequenced, spliced, and analyzed to include the examination of nucleotide composition, codon usage, compositional biases, protein-coding genes (PCGs), and RNA secondary structures. Phylogenetic relationships between L. pratensis and other species in different Heteroptera families were also examined. The mt genome was found to be a typical circular genome with a length of 16,591 bp and a total AT content of 75.1%, encoded for 13 PCGs, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (lrRNA and srRNA), and a noncoding control region. The nucleotide composition of the entire mt genome was heavily biased toward A and T. All of the tRNAs were predicted to have classic clover leaf structures, but three of the tRNAs (tRNAAsn, tRNAHis, tRNAHis) were missing the TΨC loop. The control region (2,017 bp), which was found to be located between 12S and tRNAIle, contained three tandem repeat elements. Phylogenetic analyses showed that L. pratensis is closely related to the other three examined Lygus bugs, and that it is a sister group to Apolygus and Adelphocoris. This study confirms the usability of the mt genome in phylogenesis studies pertaining to the Lygus genus, within Miridae.

Characterization of the complete mitochondrial genome of Marshallagia marshalli and phylogenetic implications for the superfamily Trichostrongyloidea

Marshallagia marshalli (Nematoda: Trichostrongylidae) infection can lead to serious parasitic gastroenteritis in sheep, goat, and wild ruminant, causing significant socioeconomic losses worldwide. Up to now, the study concerning the molecular biology of M. marshalli is limited. Herein, we sequenced the complete mitochondrial (mt) genome of M. marshalli and examined its phylogenetic relationship with selected members of the superfamily Trichostrongyloidea using Bayesian inference (BI) based on concatenated mt amino acid sequence datasets. The complete mt genome sequence of M. marshalli is 13,891 bp, including 12 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. All protein-coding genes are transcribed in the same direction. Phylogenetic analyses based on concatenated amino acid sequences of the 12 protein-coding genes supported the monophylies of the families Haemonchidae, Molineidae, and Dictyocaulidae with strong statistical support, but rejected the monophyly of the family Trichostrongylidae. The determination of the complete mt genome sequence of M. marshalli provides novel genetic markers for studying the systematics, population genetics, and molecular epidemiology of M. marshalli and its congeners.

Novel insights into mitochondrial gene rearrangement in thrips (Insecta: Thysanoptera) from the grass thrips, Anaphothrips obscurus

We sequenced the mitochondrial (mt) genome of the grass thrips, Anaphothrips obscurus, which is highly rearranged and differs from the four thrips species reported previously in the arrangement of both tRNA genes and a protein-coding gene, nad3, and in the copy number of the control region (CR). We reconstructed the phylogeny of the thrips with mt genome sequences, and used it as a framework to gain insights into mt genome evolution in thrips. It is evident that A. obscurus is less rearranged in mt genome organization than the other four known thrips. nad3 is in its ancestral location in A. obscurus but was translocated in other four thrips. Also, A. obscurus has one CR, which is ancestral to hexapods whereas other thrips have two or three CRs. All of the five thrips whose mt genomes have been sequenced to date are from the subfamily Thripinae, which represents about a quarter of the species richness in the order Thysanoptera. The high variation in mt genome organization observed in a subfamily challenges our knowledge about animal mt genomes. It remains to be investigated why mt genomes evolved so fast in the subfamily Thripinae and how mt genomes evolved in other lineages of thrips.

Complete mitochondrial genomes of two flat-backed millipedes by next-generation sequencing (Diplopoda, Polydesmida)

A lack of mitochondrial genome data from myriapods is hampering progress across genetic, systematic, phylogenetic and evolutionary studies. Here, the complete mitochondrial genomes of two millipedes, Asiomorphacoarctata Saussure, 1860 (Diplopoda: Polydesmida: Paradoxosomatidae) and Xystodesmus sp. (Diplopoda: Polydesmida: Xystodesmidae) were assembled with high coverage using Illumina sequencing data. The mitochondrial genomes of the two newly sequenced species are circular molecules of 15,644 bp and 15,791 bp, within which the typical mitochondrial genome complement of 13 protein-coding genes, 22 tRNAs and two ribosomal RNA genes could be identified. The mitochondrial genome of Asiomorphacoarctata is the first complete sequence in the family Paradoxosomatidae (Diplopoda: Polydesmida) and the gene order of the two flat-backed millipedes is novel among known myriapod mitochondrial genomes. Unique translocations have occurred, including inversion of one half of the two genomes with respect to other millipede genomes. Inversion of the entire side of a genome (trnF-nad5-trnH-nad4-nad4L, trnP, nad1-trnL2-trnL1-rrnL-trnV-rrnS, trnQ, trnC and trnY) could constitute a common event in the order Polydesmida. Last, our phylogenetic analyses recovered the monophyletic Progoneata, subphylum Myriapoda and four internal classes.

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.

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.