The complete mitochondrial genomes of two rice planthoppers, Nilaparvata lugens and Laodelphax striatellus: conserved genome rearrangement in Delphacidae and discovery of new characteristics of atp8 and tRNA genes

The effect of migration on transmission of Wolbachia in Nilaparvata lugens

Brown planthopper Nilaparvata lugens, which can transmit rice ragged stunt virus, is a serious and damaging pest to rice plants. Rice plants can protect themselves from the associated diseases of N.lugens by either suppressing or replacing N.lugens by releasing N.lugens infected by a special strain of Wolbachia wStri. The long-distance migration habit of N.lugens is one of the important precursors leading up to the large-scale occurrence of N.lugens. To study the effect of migration on the transmission of Wolbachia in N.lugens, a Wolbachia spreading dynamics model with migration of N.lugens between two patches is put forward. The existence and local stability conditions of equilibrium points of the system and its subsystems are obtained. Moreover, the effects of migration on the dynamic properties and the control of N.lugens are analyzed; the results show that the system can exhibit a bistable phenomenon, and the migration can change the stability of equilibrium infected with wStri from stable to unstable. The quantitative control methods for the migration of the insect N.lugens are proposed, which provide a theoretical guidance for future field experiments. Lastly, we use the Markov chain Monte Carlo (MCMC) method to estimate the parameters of the wild N.lugens migration model based on limited observational data; the numerical simulation results show that migration can increase the quantity of N.lugens, which is consistent with the relevant experimental results.

Polymorphisms in the hypervariable control region of the mitochondrial DNA differentiate BPH populations

The brown planthopper (BPH; Nilaparvata lugens) is one of India's most destructive pests of rice. BPH, a monophagous migratory insect, reported from all major rice-growing ecosystems of the country, is capable of traversing large distances and causing massive crop loss. A crucial step for developing viable management strategies is understanding its population dynamics. Very few reliable markers are currently available to screen BPH populations for their diversity. In the current investigation, we developed a combinatorial approach using the polymorphism present within the mitochondrial Control Region of BPH and in the nuclear genome (genomic simple sequence repeats; gSSRs) to unravel the diversity present in BPH populations collected from various rice-growing regions of India. Using two specific primer pairs, the complete Control Region (1112 to 2612 bp) was PCR amplified as two overlapping fragments, cloned and sequenced from BPH individuals representing nine different populations. Results revealed extensive polymorphism within this region due to a variable number of tandem repeats. The three selected gSSR markers also exhibited population-specific amplification patterns. Overall genetic diversity between the nine populations was high (>5%). Further, in silico double-digestion of the consensus sequences of the Control Region, with HpyCH4IV and Tsp45I restriction enzymes, revealed unique restriction fragment length polymorphisms (digital-RFLPs; dRFLPs) that differentiated all the nine BPH populations. To the best of our knowledge, this is the first report of markers developed from the Control Region of the BPH mitogenome that can differentiate populations. Eventually, such reliable and rapid marker-based identification of BPH populations will pave the way for an efficient pest management strategy.

Genome Organization and Comparative Evolutionary Mitochondriomics of Brown Planthopper, Nilaparvata lugens Biotype 4 Using Next Generation Sequencing (NGS)

Nilaparvata lugens is the main rice pest in India. Until now, the Indian N. lugens mitochondrial genome has not been sequenced, which is a very important basis for population genetics and phylogenetic evolution studies. An attempt was made to sequence two examples of the whole mitochondrial genome of N. lugens biotype 4 from the Indian population for the first time. The mitogenomes of N. lugens are 16,072 and 16,081 bp long with 77.50% and 77.45% A + T contents, respectively, for both of the samples. The mitochondrial genome of N. lugens contains 37 genes, including 13 protein-coding genes (PCGs) (cox1-3, atp6, atp8, nad1-6, nad4l, and cob), 22 transfer RNA genes, and two ribosomal RNA (rrnS and rrnL) subunits genes, which are typical of metazoan mitogenomes. However, both samples of N. lugens mitogenome in the present study retained one extra copy of the trnC gene. Additionally, we also found 93 bp lengths for the atp8 gene in both of the samples, which were 60–70 bp less than that of the other sequenced mitogenomes of hemipteran insects. The phylogenetic analysis of the 19 delphacids mitogenome dataset yielded two identical topologies when rooted with Ugyops sp. in one clade, and the remaining species formed another clade with P. maidis and M. muiri being sisters to the remaining species. Further, the genus Nilaparvata formed a separate subclade with the other genera (Sogatella, Laodelphax, Changeondelphax, and Unkanodes) of Delphacidae. Additionally, the relationship among the biotypes of N. lugens was recovered as the present study samples (biotype-4) were separated from the three biotypes reported earlier. The present study provides the reference mitogenome for N. lugens biotype 4 that may be utilized for biotype differentiation and molecular-aspect-based future studies of N. lugens.

Insights into the Evolution of Aphid Mitogenome Features from New Data and Comparative Analysis

The complete mitochondrial genomes and their rearrangement patterns can provide useful information for inferring evolutionary history of organisms. Aphids are one of the insect groups with some unique mitogenome features. In this study, to examine whether some features in aphid mitogenomes are independent species-specific evolutionary events or clade-specific events at certain taxonomic levels, we sequenced three new aphid mitogenomes (Hormaphidinae: Ceratovacuna keduensis, Pseudoregma panicola; Lachninae: Nippolachnus piri) and compared them with all known aphid mitogenomes. The three mitogenomes are 16,059–17,033 bp in length, with a set of 37 typical mitochondrial genes, a non-coding control region and a tandem repeat region. The gene orders of them are all highly rearranged. Within the subfamily Hormaphidinae, the presence of repeat region and mitogenome rearrangement in Cerataphidini species but not in the other two tribes indicate that these may be Cerataphidini-specific features. The same gene rearrangement pattern in the two Lachninae species, N. piri (Tuberolachnini) and Stomaphis sinisalicis (Stomaphidini), supports that this feature should be at least derived from the common ancestor of two tribes. Overall, our data and analyses provide new insights into the evolutionary patterns of gene rearrangement and repeat region in aphid mitogenomes, and further corroborate the potential role of gene rearrangement in elucidating the evolutionary history of different insect lineages.

Characterizing the Complete Mitochondrial Genome of Arma custos and Picromerus lewisi (Hemiptera: Pentatomidae: Asopinae) and Conducting Phylogenetic Analysis

We characterized the mitochondrial genome (mitogenome) and conducted phylogenetic analyses of 48 Hemiptera species by sequencing and analyzing the mitogenome of Arma custos (Fabricius) and Picromerus lewisi (Scott). The complete mitogenomes of the two predators were 16,024 bp and 19,587 bp in length, respectively, and it contained 37 classical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), and a control region. Most PCGs in these predators use ATN as the start codon. This research revealed that the genes of the two natural enemy species have an A + T content of 75.40% and all tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. This is the first study to compare the mitochondrial genetic structure of two predatory insects; the mitochondrial genetic structure of individual predatory insects has been sequenced in previous studies. Here, phylogenetic analysis on the basis of amino acid and nucleotide sequences of 13 mitochondrial PCGs using Bayesian inference and maximum likelihood methods were conducted to generate similar tree topologies, which suggested that the two predators with close genetic relationships belong to Asopinae subfamily. Furthermore, the monophyly of the Pentatomoidea superfamily is well accepted despite limited taxon and species sampling. Finally, their complete mitogenome provided data to establish a predator-prey food web, which is the foundation of effective pest management. Our results also enhanced the database of natural enemy insects.

A New Mitochondrial Genome of Sogatella furcifera (Horváth) (Hemiptera: Delphacidae) and Mitogenome-Wide Investigation on Polymorphisms

Simple Summary

We completed one mitogenome of white-backed planthopper (WBPH), Sogatella furcifera (Horváth), with finding heteroplasmy phenomenon confirmed by PCR reaction and Sanger sequencing method. This heteroplasmy was not observed in WBPHs (n = 24) collected from the fields, suggesting that it may be uncommon in fields. We also analyzed single nucleotide polymorphisms, insertion and deletions, and simple sequence repeats among three currently available WBPH mitogenomes of Korea and China, suggesting that identified intraspecific variations could be potential candidates for developing markers to distinguish geographical populations of WBPH including Korean and Chinese. Phylogenetic analysis of 32 mitogenomes of Delphacidae including the three WBPH mitogenomes suggested that Delphacinae seems to be monophyletic and Sogatella species including WBPH are clearly formed as one clade.

Abstract

White-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the major sap-sucking rice pests in East Asia. We have determined a new complete mitochondrial genome of WBPH collected in the Korean peninsula using NGS technology. Its length and GC percentages are 16,613 bp and 23.8%, respectively. We observed one polymorphic site, a non-synonymous change, in the COX3 gene with confirmation heteroplasmy phenomenon within individuals of WBPH by PCR amplification and Sanger sequencing, the first report in this species. In addition, this heteroplasmy was not observed in wild WBPH populations, suggesting that it may be uncommon in fields. We analyzed single nucleotide polymorphisms, insertion, and deletions, and simple sequence repeats among the three WBPH mitogenomes from Korea and China and found diverse intraspecific variations, which could be potential candidates for developing markers to distinguish geographical populations. Phylogenetic analysis of 32 mitogenomes of Delphacidae including the three WBPH mitogenomes suggested that Delphacinae seems to be monophyletic and Sogatella species including WBPH are clearly formed as one clade. In the future, it is expected that complete mitogenomes of individuals of geographically dispersed WBPH populations will be used for further population genetic studies to understand the migration pathway of WBPH.

Comparative analysis of twelve mitogenomes of Caliscelidae (Hemiptera: Fulgoromorpha) and their phylogenetic implications

Here, the complete mitochondrial genomes (mitogenomes) of 12 Caliscelidae species, Augilina tetraina, Augilina triaina, Symplana brevistrata, Symplana lii, Neosymplana vittatum, Pseudosymplanella nigrifasciata, Symplanella brevicephala, Symplanella unipuncta, Augilodes binghami, Cylindratus longicephalus, Caliscelis shandongensis, and Peltonotellus sp., were determined and comparatively analyzed. The genomes varied from 15,424 to 16,746 bp in size, comprising 37 mitochondrial genes and an A+T-rich region. The typical gene content and arrangement were similar to those of most Fulgoroidea species. The nucleotide compositions of the mitogenomes were biased toward A/T. All protein-coding genes (PCGs) started with a canonical ATN or GTG codon and ended with TAN or an incomplete stop codon, single T. Among 13 PCGs in 16 reported Caliscelidae mitogenomes, cox1 and atp8 showed the lowest and highest nucleotide diversity, respectively. All PCGs evolved under purifying selection, with atp8 considered a comparatively fast-evolving gene. Phylogenetic relationships were reconstructed based on 13 PCGs in 16 Caliscelidae species and five outgroups using maximum likelihood and Bayesian inference analyses. All species of Caliscelidae formed a steadily monophyletic group with high support. Peltonotellini was present at the basal position of the phylogenetic tree. Augilini was the sister group to Caliscelini and Peltonotellini.

Characterization, Comparative Analysis and Phylogenetic Implications of Mitogenomes of Fulgoridae (Hemiptera: Fulgoromorpha)

The complete mitogenomes of nine fulgorid species were sequenced and annotated to explore their mitogenome diversity and the phylogenetics of Fulgoridae. All species are from China and belong to five genera: Dichoptera Spinola, 1839 (Dichoptera sp.); Neoalcathous Wang and Huang, 1989 (Neoalcathous huangshanana Wang and Huang, 1989); Limois Stål, 1863 (Limois sp.); Penthicodes Blanchard, 1840 (Penthicodes atomaria (Weber, 1801), Penthicodes caja (Walker, 1851), Penthicodes variegata (Guérin-Méneville, 1829)); Pyrops Spinola, 1839 (Pyrops clavatus (Westwood, 1839), Pyrops lathburii (Kirby, 1818), Pyrops spinolae (Westwood, 1842)). The nine mitogenomes were 15,803 to 16,510 bp in length with 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (A + T-rich region). Combined with previously reported fulgorid mitogenomes, all PCGs initiate with either the standard start codon of ATN or the nonstandard GTG. The TAA codon was used for termination more often than the TAG codon and the incomplete T codon. The nad1 and nad4 genes varied in length within the same genus. A high percentage of F residues were found in the nad4 and nad5 genes of all fulgorid mitogenomes. The DHU stem of trnV was absent in the mitogenomes of all fulgorids sequenced except Dichoptera sp. Moreover, in most fulgorid mitogenomes, the trnL2, trnR, and trnT genes had an unpaired base in the aminoacyl stem and trnS1 had an unpaired base in the anticodon stem. The similar tandem repeat regions of the control region were found in the same genus. Phylogenetic analyses were conducted based on 13 PCGs and two rRNA genes from 53 species of Fulgoroidea and seven outgroups. The Bayesian inference and maximum likelihood trees had a similar topological structure. The major results show that Fulgoroidea was divided into two groups: Delphacidae and ((Achilidae + (Lophopidae + (Issidae + (Flatidae + Ricaniidae)))) + Fulgoridae). Furthermore, the monophyly of Fulgoridae was robustly supported, and Aphaeninae was divided into Aphaenini and Pyropsini, which includes Neoalcathous, Pyrops, Datua Schmidt, 1911, and Saiva Distant, 1906. The genus Limois is recovered in the Aphaeninae, and the Limoisini needs further confirmation; Dichoptera sp. was the earliest branch in the Fulgoridae.

Structural Features and Phylogenetic Implications of 11 New Mitogenomes of Typhlocybinae (Hemiptera: Cicadellidae)

To explore the characteristics of mitogenomes and discuss the phylogenetic relationships and molecular evolution of the six tribes within Typhlocybinae, 11 complete mitogenomes are newly sequenced and comparatively analyzed. In all of these complete mitogenomes, the number and order of the genes are highly conserved in overall organization. The PCGs initiate with ATN/TTG/GTG and terminate with TAA/TAG/T. Almost all tRNAs are folded into the typical clover-leaf secondary structure. The control region is always variable in length and in numbers of multiple tandem repeat units. The atp8 and nad2 exhibits the highest evolution rate among all the PCGs. Phylogenetic analyses based on whole mitogenome sequences, with three different datasets, using both maximum likelihood and Bayesian methods, indicate the monophyly of Typhlocybinae and its inner tribes, respectively, except for Typhlocybini and Zyginellini that are paraphyletic. Finally, we confirm that Erythroneurini is a subtribe of Dikraneurini.

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.

Characterization of Three Complete Mitogenomes of Flatidae (Hemiptera: Fulgoroidea) and Compositional Heterogeneity Analysis in the Planthoppers' Mitochondrial Phylogenomics

Although sequences of mitogenomes have been widely used for investigating phylogenetic relationship, population genetics, and biogeography in many members of Fulgoroidea, only one complete mitogenome of a member of Flatidae has been sequenced. Here, the complete mitogenomes of Cerynia lineola, Cromna sinensis, and Zecheuna tonkinensis are sequenced. The gene arrangements of the three new mitogenomes are consistent with ancestral insect mitogenomes. The strategy of using mitogenomes in phylogenetics remains in dispute due to the heterogeneity in base composition and the possible variation in evolutionary rates. In this study, we found compositional heterogeneity and variable evolutionary rates among planthopper mitogenomes. Phylogenetic analysis based on site-homogeneous models showed that the families (Delphacidae and Derbidae) with high values of Ka/Ks and A + T content tended to fall together at a basal position on the trees. Using a site-heterogeneous mixture CAT + GTR model implemented in PhyloBayes yielded almost the same topology. Our results recovered the monophyly of Fulgoroidea. In this study, we apply the heterogeneous mixture model to the planthoppers’ phylogenetic analysis for the first time. Our study is based on a large sample and provides a methodological reference for future phylogenetic studies of Fulgoroidea.

The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications

To reveal mtgenome characterizations and reconstruct phylogenetic relationships of Hylicinae, the complete mtgenomes of four hylicine species, including Nacolus tuberculatus, Hylica paradoxa, Balala fujiana, and Kalasha nativa, were sequenced and comparatively analyzed for the first time. We also carried out the richest (11) subfamily sampling of Cicadellidae to date, and reconstructed phylogenetic relationships of Membracoidea among 61 species based on three datasets using maximum likelihood and Bayesian inference analyses. All new sequenced mtgenomes are molecules ranging from 14,918 to 16,221 bp in length and are double stranded, circular in shape. The gene composition and arrangement of these mtgenomes are consistent with members of Membracoidea. Among 13 protein-coding genes, most show typical ATN start codons and TAR (TAA/TAG) or an incomplete stop codon T-, and several genes start by TTG/GTG. Results of the analysis for sliding window, nucleotide diversity, and nonsynonymous substitution/synonymous substitution indicate cox1 is a comparatively slower-evolving gene while atp8 is the fastest gene. In line with previous researches, phylogenetic results indicate that treehopper families are paraphyletic with respect to family Cicadellidae and also support the monophyly of all involved subfamilies including Hylicinae. Relationships among the four hylicine genera were recovered as (Hylica + (Nacolus + (Balala + Kalasha))).

Characterization of the Complete Mitochondrial Genomes of Two Species with Preliminary Investigation on Phylogenetic Status of Zyginellini (Hemiptera: Cicadellidae: Typhlocybinae)

To explore the characteristics of mitogenomes and reveal phylogenetic relationships of the tribes of Zyginellini and Typhlocybini in Typhlocybinae, mitogenomes of two species of the Zyginellini, Parazyginella tiani and Limassolla sp., were sequenced. Mitogenomes of both species contain 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a large non-coding region (A + T-rich region). These characteristics are similar to other Membracoidea mitogenomes. All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA/G or with an incomplete T codon. All tRNAs have the typical clover-leaf structure, except trnS1 which has a reduced DHU arm and the acceptor stem of trnR is 5 or 6 bp in some species, an unusual feature here reported for the first time in Typhlocybinae. The A + T-rich region is highly variable in length and in numbers of tandem repeats present. Our analyses indicate that nad6 and atp6 exhibit higher evolutionary rates compared to other PCGs. Phylogenetic analyses by both maximum likelihood and Bayesian methods based on 13 protein-coding genes of 12 species of Typhlocybinae suggest that Zyginellini are paraphyletic with respect to Typhlocybini.

Stable Introduction of Plant-Virus-Inhibiting Wolbachia into Planthoppers for Rice Protection

Progress has been made in developing the maternally inherited endosymbiotic bacterium Wolbachia as a tool for protecting humans from mosquito-borne diseases. In contrast, Wolbachia-based approaches have not yet been developed for the protection of plants from insect pests and their associated diseases, with a major challenge being the establishment of artificial Wolbachia infections expressing desired characteristics in the hemipterans that transmit the majority of plant viruses. Here, we report stable introduction of Wolbachia into the brown planthopper, Nilaparvata lugens, the most destructive rice pest that annually destroys millions of hectares of staple crops. The Wolbachia strain wStri from the small brown planthopper, Laodelphax striatellus, was transferred to this new host, where it showed high levels of cytoplasmic incompatibility, enabling rapid invasion of laboratory populations. Furthermore, wStri inhibited infection and transmission of Rice ragged stunt virus and mitigated virus-induced symptoms in rice plants, opening up the development of Wolbachia-based strategies against major agricultural pests and their transmitted pathogens. VIDEO ABSTRACT.

The challenge of Coccidae (Hemiptera: Coccoidea) mitochondrial genomes: The case of Saissetia coffeae with novel truncated tRNAs and gene rearrangements

There have been few reports of complete mitochondrial genomes (mitogenomes) of scale insects, and it has been indicated that complex and novel structures in their mitogenomes may lead to difficulties in sequencing, assembly and annotation. Transfer RNAs (tRNAs) usually possess typical cloverleaf secondary structures, and truncated tRNAs are rarely found in insect mitogenomes. Here, we report a complete Saissetia coffeae mitogenome (15,389 bp) with high A+T content (84.7%) sequenced by next-generation sequencing (NGS) methods. Genes in the mitogenome were annotated, and nine tRNAs were not found using MITOS. Most of the detected tRNAs were significantly truncated without the dihydrouridine (DHU) arm or the TΨC (T) arm. In addition, the 9 "lost" tRNAs containing mismatched base pairs were retrieved based on the tRNA annotation workflow for Coccidae described in our study. The gene arrangement in the Saissetia coffeae mitogenome was significantly different from that in other hemipteran insects. Additionally, Bayesian and maximum likelihood trees based on the mitochondrial genes showed a long branch of the Saissetia lineage, indicating significant nonsynonymous substitutions or high evolutionary rates in the Saissetia lineage. We provide a reference mitogenome for the assembly and annotation of the Coccidae mitogenome and offer insights into the evolution of scale insects.

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.

Contribution to the mitogenome diversity in Delphacinae: Phylogenetic and ecological implications

We document the complete (or nearly complete) mitogenomes of 20 Delphacidae taxa, and together with 17 other delphacid mitogenomes currently in GenBank, to reconstruct the phylogeny of the Delphacinae and to investigate mitogenome differences among members of Delphacini, Tropidocephalini and Saccharosydnini. The mitogenomes of the 20 species encode the complete set of 37 genes usually found in animal mitogenomes. The length of complete mitogenomes in Delphacinae ranges from 15,531 to 16,231 bp. The gene order of all newly sequenced mitogenomes are identical, and the mitogenome gene order of Stenocranus matsumurai Metcalf in Stenocraninae has a transposition of tRNA^Thr. The two-clade system in Tropidocephalini was supported with high value (PP = 1, BS = 100), and the monophyly of Bambusiphaga was recovered in this study. Finally, we found that the host shift from plants with a C₃ to a C₄ photosynthetic pathway appears to have occurred independently in several clades.

The complete mitochondrial genome of Laodelphax striatellus (Fallén, 1826) (Hemiptera: Delphacidae) collected in a mid-Western part of Korean peninsula

Laodelphax striatellus (Fallén, 1826) is one of key rice pests in Northeast Asia. We have determined the mitochondrial genome of L. striatellus collected in a mid-western part of Korean peninsula. The circular mitogenome of L. striatellus is 16,359 bp long including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 1,972 bp. The base composition was AT-biased (77.3%). In comparison of the two Chinese L. striatellus mitogenomes with Korean mitogenome, total 140 and 40 single nucleotide polymorphisms and 166 and 118 insertions and deletions were identified, presenting intra-species variations based on geographical distribution.

The complete mitochondrial genome of Laodelphax striatellus (Fallén, 1826) (Hemiptera: Delphacidae) collected in a southern part of Korean peninsula

A 16,359 bp mitochondrial genome of Laodelphax striatellus collected in a southern part of Korean peninsula was completed and their intraspecies variations were compared with Korean and Chinese L. striatellus mitogenomes. The circular mitogenome of L. striatellus contains 13 protein-coding genes, two rRNA genes, 22 tRNAs, and a single large non-coding region of 1,972 bp. The base composition was AT-biased (77.3%). It has only one single nucleotide polymorphism (SNP) in AT-rich region compared to other Korean mitogenome, but total 41 and 141 SNPs and 118 and 166 insertions and deletions, respectively, compared to two Chinese mitogenomes, suggesting the possibility of tracing migration based on geographic genetic diversity.

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill-coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function-altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal-level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.

The first two mitochondrial genomes of the family Aphelinidae with novel gene orders and phylogenetic implications

Chalcidoidea is one of the most diverse group in Hymenoptera by possessing striking mitochondrial gene arrangement. By using next generation sequencing method, the first two nearly complete mitochondrial genomes in the family Aphelinidae (Insecta, Hymenopetra, Chalcidoidea) were obtained in this study. Almost all previously sequenced mitochondrial genome of Chalcidoidea species have a large inversion including six genes (atp6-atp8-trnD-trnK-cox2-trnL2-cox1) as compared with ancestral mitochondrial genome, but these two Encarsia mitochondrial genomes had a large inversion including nine genes (nad3-trnG-atp6-atp8-trnD-trnK-cox2-trnL2-cox1), which was only congruent with the species in the genus Nasonia. Moreover, we found that one shuffling changes (trnD and trnK) happened in the species E. obtusiclava but not in another species E. formosa within the same genus, of which such shuffling within the same genus at this region was only detected in Polisters within Insecta. Phylogenetic analysis displayed that different data matrix (13PCG+ 2 rRNA or 13 PCG) and inference methods (BI or ML) indicate the identical topology with high nodal supports that Aphelinidae formed a sister group with (Trichogrammatidae + Aganoidae) and the monophyly of Pteramalidae. Our results also indicated the validity of assembling and feasibility of next-generation technology to obtain the mitochondrial genomes of parasitic Hymenoptera.

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.

Sequencing and analysis of the complete mitochondrial genome of Changeondelphax velitchkovskyi (Hemiptera: Fulgoroidea)

We determined the complete mitochondrial genome sequence of Changeondelphax velitchkovskyi (Melichar 1913). The complete mitogenome sequence of C. velitchkovskyi was observed to be a circular molecule 16,449 bp long and consists of 13 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes and 22 transfer RNA (tRNA) genes (GenBank accession no. MG049916). This nucleotide composition is biased toward adenine and thymine (75.7% A + T). The A + T-rich region is found between rrnS and trnI, and this entire region was 1781 bp long.

The complete mitochondrial genome sequence of the corn planthopper, Peregrinus maidis (Hemiptera: Fulgoroidea)

In this study, we analyzed the complete mitochondrial genome sequence of the corn planthopper, Peregrinus maidis. The complete mitogenome sequence of P. maidis was observed to be a circular molecule 16,279 bp long and consisting of 13 protein-coding genes (PCG), 2 ribosomal RNA (rRNA) genes and 22 transfer RNA (tRNA) genes (GenBank accession no. MG049917). The nucleotide composition is biased toward adenine and thymine (77.8% A + T). The A + T-rich region was found between rrnS and trnI, and this entire region was 1596 bp long.

C. Costa M, A. Andow D, Paula DP. Mitogenome sequence accuracy using different elucidation methods

Mitogenome sequences are highly desired because they are used in several biological disciplines. Their elucidation has been facilitated through the development of massive parallel sequencing, accelerating their deposition in public databases. However, sequencing, assembly and annotation methods might induce variability in their quality, raising concerns about the accuracy of the sequences that have been deposited in public databases. In this work we show that different sequencing methods (number of species pooled in a library, insert size and platform) and assembly and annotation methods generated variable completeness and similarity of the resulting mitogenome sequences, using three species of predaceous ladybird beetles as models. The identity of the sequences varied considerably depending on the method used and ranged from 38.19 to 90.1% for Cycloneda sanguinea, 72.85 to 91.06% for Harmonia axyridis and 41.15 to 93.60% for Hippodamia convergens. Dissimilarities were frequently found in the non-coding A+T rich region, but were also common in coding regions, and were not associated with low coverage. Mitogenome completeness and sequence identity were affected by the sequencing and assembly/annotation methods, and high within-species variation was also found for other mitogenome depositions in GenBank. This indicates a need for methods to confirm sequence accuracy, and guidelines for verifying mitogenomes should be discussed and developed by the scientific community.

Comparative mitogenomic analyses of three North American stygobiont amphipods of the genus Stygobromus (Crustacea: Amphipoda)

The mitochondrial genomes of three North American stygobiont amphipods Stygobromus tenuis potomacus, S. foliatus and S. indentatus collected from Caroline County, VA, were sequenced using a shotgun sequencing approach on an Illumina NextSeq500 (Illumina Inc., San Diego, CA). All three mitogenomes displayed 13 protein-coding genes, 22 tRNAs and two rRNAs typical of metazoans. While S. tenuis and S. indentatus displayed identical gene orders similar to the pancrustacean ground pattern, S. foliatus displayed a transposition of the trnL2-cox2 genes to after atp8-atp6. In addition, a short atp8 gene, longer rrnL gene and large inverted repeat within the Control Region distinguished S. foliatus from S. tenuis potomacus and S. indentatus. Overall, it appears that gene order varies considerably among amphipods, and the addition of these Stygobromus mitogenomes to the existing sequenced amphipod mitogenomes will prove useful for characterizing evolutionary relationships among various amphipod taxa, as well as investigations of the evolutionary dynamics of the mitogenome in general.

All 37 Mitochondrial Genes of Aphid Aphis craccivora Obtained from Transcriptome Sequencing: Implications for the Evolution of Aphids

The availability of mitochondrial genome data for Aphididae, one of the economically important insect pest families, in public databases is limited. The advent of next generation sequencing technology provides the potential to generate mitochondrial genome data for many species timely and cost-effectively. In this report, we used transcriptome sequencing technology to determine all the 37 mitochondrial genes of the cowpea aphid, Aphis craccivora. This method avoids the necessity of finding suitable primers for long PCRs or primer-walking amplicons, and is proved to be effective in obtaining the whole set of mitochondrial gene data for insects with difficulty in sequencing mitochondrial genome by PCR-based strategies. Phylogenetic analyses of aphid mitochondrial genome data show clustering based on tribe level, and strongly support the monophyly of the family Aphididae. Within the monophyletic Aphidini, three samples from Aphis grouped together. In another major clade of Aphididae, Pterocomma pilosum was recovered as a potential sister-group of Cavariella salicicola, as part of Macrosiphini.

Mitochondrial genomes of two phlebotomine sand flies, Phlebotomus chinensis and Phlebotomus papatasi (Diptera: Nematocera), the first representatives from the family Psychodidae

BACKGROUND

Leishmaniasis is a worldwide but neglected disease of humans and animal transmitted by sand flies, vectors that also transmit other important diseases. Mitochondrial genomes contain abundant information for population genetic and phylogenetic studies, important in disease management. However, the available mitochondrial sequences of these crucial vectors are limited, emphasizing the need for developing more mitochondrial genetic markers.

METHODS

The complete mitochondrial genome of Phlebotomus chinensis was amplified in eight fragments and sequenced using primer walking. The mitochondrial genome of Phlebotomus papatasi was reconstructed from whole-genome sequencing data available on Genbank. The phylogenetic relationship of 24 selected representatives of Diptera was deduced from codon positions 1 and 2 for 13 protein coding genes, using Bayesian inference (BI) and maximum likelihood (ML) methods.

RESULTS

We provide the first Phlebotomus (P. chinensis and P. papatasi) mitochondrial genomes. Both genomes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region. The gene order of Phlebotomus mitochondrial genomes is identical with the ancestral gene order of insect. Phylogenetic analyses demonstrated that Psychodidae and Tanyderidae are sister taxa. Potential markers for population genetic study of Phlebotomus species were also revealed.

CONCLUSION

The generated mitochondrial genomes of P. chinensis and P. papatasi represent a useful resource for comparative genomic studies and provide valuable future markers for the population genetic study of these important Leishmania vectors. Our results also preliminary demonstrate the phylogenetic placement of Psychodidae based on their mitochondrial genomes.

Intraspecific and Interspecific Variations in the Mitochondrial Genomes of Nilaparvata (Hemiptera: Delphacidae)

Planthoppers in the genus Nilaparvata Distant are serious pests of rice and many other crops in tropical and temperate Asia, and northern Australia. In this study, the mitochondrial genomes of four Nilaparvata planthoppers were sequenced, three in Nilaparvata lugens Stål and one in Nilaparvata muiri China. Mitochondrial genome of Nilaparvata contain the standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region. The nucleotide composition of Nilaparvata mitochondrial sequence is biased toward adenine and thymine, and the amino acid composition is affected to a similar degree by the bias to AT. We compare the four mitochondrial genomes and find intra- and interspecific variation in gene length, base composition, nucleotide and amino acid substitutions, intergenic spacer length, and gene overlap. The intra- and interspecific variations reveal that nucleotide and amino acid substitutions in mitochondrial protein-coding genes make a contribution to the formation of various insect biotypes in one species. Furthermore, the accumulation of nonsynonymous substitutions in the mitochondrial protein-coding genes, as well as differences in start codons, the length of intergenic spacers, and gene overlap regions contribute to differences between the two species investigated here. In addition, cox is the most conserved gene family and nad4-nad4l cluster is variable in Nilaparvata mitochondrial genes for the intra- and interspecific variation.

Characterization of 12 Novel Microsatellite Markers of Sogatella furcifera (Hemiptera: Delphacidae) Identified From Next-Generation Sequence Data

The white-backed planthopper, Sogatella furcifera (Horváth) (Hemiptera: Delphacidae), is a major pest of rice and has long-range migratory behavior in Asia. Microsatellite markers (simple sequence repeats) have been widely used to determine the origins and genetic diversity of insect pests. We identified novel microsatellite loci for S. furcifera samples collected from Laos, Vietnam, and three localities in Bangladesh from next-generation Roche 454 pyrosequencing data. Size polymorphism at 12 microsatellite loci was verified for 40 adult individuals collected from Shinan, South Korea. The average number of alleles per locus was 7.92. The mean values of observed (H(o)) and expected heterozygosities (H(E)) were 0.615 and 0.757, respectively. These new microsatellite markers will be a resource for future ecological genetic studies of S. furcifera samples across more broad geographic regions in Asia and may assist in estimations of genetic differentiation and gene flow among populations for implementation of more effective management strategies to control this serious rice pest.

Hemipteran mitochondrial genomes: features, structures and implications for phylogeny

The study of Hemipteran mitochondrial genomes (mitogenomes) began with the Chagas disease vector, Triatoma dimidiata, in 2001. At present, 90 complete Hemipteran mitogenomes have been sequenced and annotated. This review examines the history of Hemipteran mitogenomes research and summarizes the main features of them including genome organization, nucleotide composition, protein-coding genes, tRNAs and rRNAs, and non-coding regions. Special attention is given to the comparative analysis of repeat regions. Gene rearrangements are an additional data type for a few families, and most mitogenomes are arranged in the same order to the proposed ancestral insect. We also discuss and provide insights on the phylogenetic analyses of a variety of taxonomic levels. This review is expected to further expand our understanding of research in this field and serve as a valuable reference resource.

Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): an unusual elongation of rrnS gene

Two nearly complete mitochondrial genomes (mitogenomes) of wheat stem borers, Cephus pygmeus and Cephus sareptanus (Hymenoptera: Cephidae), were sequenced, characterised and compared with the previously known mitogenome of Cephus cinctus. The gene orders are mostly conserved, except for translocation of trnM and swapped position of trnI and trnQ. An A+T bias was found, but a deviation from strand asymmetry was also detected on the J strand. All protein coding genes (PCGs) are initiated by ATN codons, except for nad1, nad2 and atp8, and all are terminated with TAA, TA- or T- as a stop codon. The predicted secondary structures of rrnS and rrnL genes are mostly consistent with reported hymenopteran species. However, an unusual elongation in rrnS, not know elsewhere in the order, was discovered in Cephus species. Three autonomous sequences detected in domains I and II are mainly responsible for the length expansions.

Is it an ant or a butterfly? Convergent evolution in the mitochondrial gene order of Hymenoptera and Lepidoptera

Insect mitochondrial genomes (mtDNA) are usually double helical and circular molecules containing 37 genes that are encoded on both strands. The arrangement of the genes is not constant for all species, and produces distinct gene orders (GOs) that have proven to be diagnostic in defining clades at different taxonomic levels. In general, it is believed that distinct taxa have a very low chance of sharing identically arranged GOs. However, examples of identical, homoplastic local rearrangements occurring in distinct taxa do exist. In this study, we sequenced the complete mtDNAs of the ants Formica fusca and Myrmica scabrinodis (Formicidae, Hymenoptera) and compared their GOs with those of other Insecta. The GO of F. fusca was found to be identical to the GO of Dytrisia (the largest clade of Lepidoptera). This finding is the first documented case of an identical GO shared by distinct groups of Insecta, and it is the oldest known event of GO convergent evolution in animals. Both Hymenoptera and Lepidoptera acquired this GO early in their evolution. Using a phylogenetic approach combined with new bioinformatic tools, the chronological order of the evolutionary events that produced the diversity of the hymenopteran GOs was determined. Additionally, new local homoplastic rearrangements shared by distinct groups of insects were identified. Our study showed that local and global homoplasies affecting the insect GOs are more widespread than previously thought. Homoplastic GOs can still be useful for characterizing the various clades, provided that they are appropriately considered in a phylogenetic and taxonomic context.

Identification and characterization of microRNAs in small brown planthopper (Laodephax striatellus) by next-generation sequencing

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level and are thought to play critical roles in many metabolic activities in eukaryotes. The small brown planthopper (Laodephax striatellus Fallén), one of the most destructive agricultural pests, causes great damage to crops including rice, wheat, and maize. However, information about the genome of L. striatellus is limited. In this study, a small RNA library was constructed from a mixed L. striatellus population and sequenced by Solexa sequencing technology. A total of 501 mature miRNAs were identified, including 227 conserved and 274 novel miRNAs belonging to 125 and 250 families, respectively. Sixty-nine conserved miRNAs that are included in 38 families are predicted to have an RNA secondary structure typically found in miRNAs. Many miRNAs were validated by stem-loop RT-PCR. Comparison with the miRNAs in 84 animal species from miRBase showed that the conserved miRNA families we identified are highly conserved in the Arthropoda phylum. Furthermore, miRanda predicted 2701 target genes for 378 miRNAs, which could be categorized into 52 functional groups annotated by gene ontology. The function of miRNA target genes was found to be very similar between conserved and novel miRNAs. This study of miRNAs in L. striatellus will provide new information and enhance the understanding of the role of miRNAs in the regulation of L. striatellus metabolism and development.

The complete mitochondrial genome sequence of Sogatella furcifera (Horváth) and a comparative mitogenomic analysis of three predominant rice planthoppers

The white-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the most destructive pests of rice crops in many Asian countries. Using long-PCR and shotgun library methods, we sequenced the entire mitochondrial genomes (mt-genomes) of two WBPH individuals. Total lengths of the mt-genome of the two WBPH individuals were 16,612 bp and 16,654 bp with an identical AT content of 76.19%. Among the 13 protein coding genes (PCGs), only nad5 used an atypical initiation codon GTG. Most of the tRNA genes had the typical cloverleaf secondary structure except that the dihydrouridine (DHU) arms in two trnS genes and the TΨC arm of trnG gene did not form a stable stem-loop structure. Similar to the brown planthopper (BPH), Nilaparvata lugens (Stål), and the small brown planthopper (SBPH), Laodelphax striatellus (Fallén), some extraordinary features were observed in the WBPH mt-genome. These include similar gene rearrangement pattern, unusually short length of the atp8 gene and variable numbers of tandem repeat (VNTR) structure in control region. Interestingly, the same tandem repeat unit with stable secondary structure appeared in two different planthoppers, WBPH and SBPH, which belong to two different genera of the Delphacidae. This peculiar feature provides a direct evidence for the close relationship between the two planthoppers and updates our understanding of the evolutionary characteristics of mitochondrial control region. Comparison with two other predominant rice planthoppers (BPH and SBPH) revealed that different PCGs of mitochondria exhibit different evolutionary patterns.