The complete mitochondrial genomes of two band-winged grasshoppers, Gastrimargus marmoratus and Oedaleus asiaticus

The complete mitochondrial genome of the grasshopper Chorthippus dubius (Zub.) (Orthoptera: Acrididae: Gomphocerinae): detailed characterization and phylogenetic position

Chorthippus dubius (Zub.) is one of the dominant grasshopper species. The limited data on Ch. dubius (Zub.) has impeded further understanding of its genetic characteristics and molecular detection. In this study, we analyzed the mitogenome of Ch. dubius (Zub.), which was 15,561 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and an AT-rich region. The entire mtDNA exhibited a strong AT bias, with an overall A+T content of 74.8%. The relative synonymous codon usage (RSCU) analysis revealed UUA (L) as the most frequently used codon. All the PCGs evolved under purifying selection (Ka/Ks <0.5), with ATP8 gene exhibited the highest Ka/Ks ratio. Maximum likelihood (ML) and Bayesian inference (BI) analyses reconstructed two topologically similar phylogenetic trees, and supported the monophly of the six subfamilies in Acrididae. Our results indicated two stable clades of the six subfamilies, with Oedipodinae emerging as the ancestral taxon and being sister group to the remaining taxa. In the alternative phylogenetic lineage, the remaining five subfamilies clustered the following relationship: Gomphocernae + (Acridinae + (Calliptaminae + (Melanoplinae + Oxyinae))). Both phylogenetic trees exhibited a closer relationship between Chorthippus dubius (Zub.) and Chorthippus aethalinus, members from the same genus.

A New Species of Scymnus (Coleoptera, Coccinellidae) from Pakistan with Mitochondrial Genome and Its Phylogenetic Implications

In this study, a new species of the subgenus Pullus belonging to the Scymnus genus from Pakistan, Scymnus (Pullus) cardi sp. nov., was described and illustrated, with information on its distribution, host plants, and prey. Additionally, the completed mitochondrial genome (mitogenome) of the new species using high-throughput sequencing technology was obtained. The genome contains the typical 37 genes (13 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs) and a non-coding control region, and is arranged in the same order as that of the putative ancestor of beetles. The AT content of the mitogenome is approximately 85.1%, with AT skew and GC skew of 0.05 and -0.43, respectively. The calculated values of relative synonymous codon usage (RSCU) determine that the codon UUA (L) has the highest frequency. Furthermore, we explored the phylogenetic relationship among 59 representatives of the Coccinellidae using Bayesian inference and maximum likelihood methods, the results of which strongly support the monophyly of Coccinellinae. The phylogenetic results positioned Scymnus (Pullus) cardi in a well-supported clade with Scymnus (Pullus) loewii and Scymnus (Pullus) rubricaudus within the genus Scymnus and the tribe Scymnini. The mitochondrial sequence of S. (P.) cardi will contribute to the mitochondrial genome database and provide helpful information for the identification and phylogeny of Coccinellidae.

Comparative analysis of mitogenomes among three species of grasshoppers (Orthoptera: Acridoidea: Gomphocerinae) and their phylogenetic implications

Whole mitochondrial genomes have been widely used in phylogenetic analysis, population genetics and biogeography studies. This study sequenced and characterized three complete mitochondrial genomes (Dasyhippus peipingensis, Myrmeleotettix palpalis, Aeropedellus prominemarginis) and determined their phylogenetic position in Acrididae. The length of the mitochondrial genomes ranged from 15,621-15,629 bp and composed of 13 PCGs, 2 rRNA, 22 tRNA genes and an AT control region. The arrangement and structure of the mitochondrial genomes were similar to those of other invertebrates. Comparative genomics revealed that the three mitochondrial genomes were highly conserved in terms of gene size, structure, and codon usage, all PCGs were purified selections with an ATN start codon and a TAN stop codon. All tRNAs could be folded into the typical clover-leaf structure, except tRNA Ser (AGN) that lacked a dihydrouridine (DHU) arm. Phylogenetic analysis based on 13 PCGs of 34 Acrididae species and seven outgroup species revealed that differences in the shape of antennae within the family Acrididae should be given less weight as a taxonomic character for higher-level classification. Moreover, the divergence time estimates indicates that in Gomphocerinae, the species with clubbed antennae were formed within the nearest 18 Mya, and Pacris xizangensis is more ancient.

An Analysis of the Possible Migration Routes of Oedaleus decorus asiaticus Bey-Bienko (Orthoptera: Acrididae) from Mongolia to China

Simple Summary

Airflow is very important for the long-distance migration of O. decorus asiaticus, and wind shear, in particular, is the main factor related to forced landing. Analyzing the weather records, we found that the northwest wind prevailed when the population invaded. Specifically, from July to August, a large number of emerging adults appeared in the source areas of Mongolia, and the large-scale northwest wind was the decisive condition for the successful long-distance migration of O. decorus asiaticus. The species has a strong migratory ability, flying along the airflow for several nights. If the northwest air current meets the southwest warm current going north, a large number of O. decorus asiaticus will drop due to wind shear, and then a major outbreak will occur. Analysis of the source of the insects shows that the O. decorus asiaticus break outs in China may have originated from Mongolia. They were brought into China by the southerly airflow at night, and they likely made a forced landing in Beijing due to wind shear, sinking airflow, rainfall and other reasons. In summary, through analysis of the insect’s prevalence and the meteorological conditions in Mongolia, we can provide a basis for predicting the occurrence of O. decorus asiaticus in China.

Abstract

Oedaleus decorus asiaticus (Bey-Bienko) is a destructive pest in grasslands and adjacent farmland in northern China, Mongolia, and other countries in Asia. It has been supposed that this insect pest can migrate a long distance and then induce huge damages, however, the migration mechanism is still unrevealed. The current study uses insect light trap data from Yanqing (Beijing), together with regional meteorological data to determine how air flow contributes to the long-distance migration of O. decorus asiaticus. Our results indicate that sinking airflow is the main factor leading to the insects’ forced landing, and the prevailing northwest wind was associated with O. decorus asiaticus taking off in the northwest and moving southward with the airflow from July to September. Meanwhile, the insects have a strong migratory ability, flying along the airflow for several nights. Thus, when the airflow from the northwest met the northward-moving warm current from the southwest, a large number of insects were dropped due to sinking airflow, resulting in a large outbreak. Our simulations suggest that the source of the grasshoppers involved in these outbreaks during early 2000s in northern China probably is in Mongolia, and all evidence indicates that there are two important immigrant routes for O. decorus asiaticus migration from Mongolia to Beijing. These findings improves our understanding of the factors guiding O. decorus asiaticus migration, providing valuable information to reduce outbreaks in China that have origins from outside the country.

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

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

Comparative Analysis of Mitogenomes among Five Species of Filchnerella (Orthoptera: Acridoidea: Pamphagidae) and Their Phylogenetic and Taxonomic Implications

Mitogenomes have been widely used for exploring phylogenetic analysis and taxonomic diagnosis. In this study, the complete mitogenomes of five species of Filchnerella were sequenced, annotated and analyzed. Then, combined with other seven mitogenomes of Filchnerella and four of Pamphagidae, the phylogenetic relationships were reconstructed by maximum likelihood (ML) and Bayesian (BI) methods based on PCGs+rRNAs. The sizes of the five complete mitogenomes are Filchnerella sunanensis 15,656 bp, Filchnerella amplivertica 15,657 bp, Filchnerella nigritibia 15,661 bp, Filchnerella pamphagoides 15,661 bp and Filchnerella dingxiensis 15,666 bp. The nucleotide composition of mitogenomes is biased toward A+T. All tRNAs could be folded into the typical clover-leaf structure, except that tRNA Ser (AGN) lacked a dihydrouridine (DHU) arm. The phylogenetic relationships of Filchnerella species based on mitogenome data revealed a general pattern of wing evolution from long wing to increasingly shortened wing.

Increasing 28 mitogenomes of Ephemeroptera, Odonata and Plecoptera support the Chiastomyaria hypothesis with three different outgroup combinations

BACKGROUND

The phylogenetic relationships of Odonata (dragonflies and damselflies) and Ephemeroptera (mayflies) remain unresolved. Different researchers have supported one of three hypotheses (Palaeoptera, Chiastomyaria or Metapterygota) based on data from different morphological characters and molecular markers, sometimes even re-assessing the same transcriptomes or mitochondrial genomes. The appropriate choice of outgroups and more taxon sampling is thought to eliminate artificial phylogenetic relationships and obtain an accurate phylogeny. Hence, in the current study, we sequenced 28 mt genomes from Ephemeroptera, Odonata and Plecoptera to further investigate phylogenetic relationships, the probability of each of the three hypotheses, and to examine mt gene arrangements in these species. We selected three different combinations of outgroups to analyze how outgroup choice affected the phylogenetic relationships of Odonata and Ephemeroptera.

METHODS

Mitochondrial genomes from 28 species of mayflies, dragonflies, damselflies and stoneflies were sequenced. We used Bayesian inference (BI) and Maximum likelihood (ML) analyses for each dataset to reconstruct an accurate phylogeny of these winged insect orders. The effect of outgroup choice was assessed by separate analyses using three outgroups combinations: (a) four bristletails and three silverfish as outgroups, (b) five bristletails and three silverfish as outgroups, or (c) five diplurans as outgroups.

RESULTS

Among these sequenced mitogenomes we found the gene arrangement IMQM in Heptageniidae (Ephemeroptera), and an inverted and translocated tRNA-Ile between the 12S RNA gene and the control region in Ephemerellidae (Ephemeroptera). The IMQM gene arrangement in Heptageniidae (Ephemeroptera) can be explained via the tandem-duplication and random loss model, and the transposition and inversion of tRNA-Ile genes in Ephemerellidae can be explained through the recombination and tandem duplication-random loss (TDRL) model. Our phylogenetic analysis strongly supported the Chiastomyaria hypothesis in three different outgroup combinations in BI analyses. The results also show that suitable outgroups are very important to determining phylogenetic relationships in the rapid evolution of insects especially among Ephemeroptera and Odonata. The mt genome is a suitable marker to investigate the phylogeny of inter-order and inter-family relationships of insects but outgroup choice is very important for deriving these relationships among winged insects. Hence, we must carefully choose the correct outgroup in order to discuss the relationships of Ephemeroptera and Odonata.

Mitogenome and phylogenetic analysis of typhlocybine leafhoppers (Hemiptera: Cicadellidae)

Mitogenomes have been widely used to estimate phylogenetic relationships among insects and provide data useful for augmenting traditional morphological characters in delimiting species. Here, complete mitogenome sequences of two closely related typhlocybine leafhoppers, Cassianeuracassiae (Ahmed, 1970) and C. bimaculata Dworakowska, 1984, were obtained and found to be 15,423 bp and 14,597 bp in length, respectively. The gene order was found to be similar to other published leafhopper mitogenomes, but the control region of C.bimaculata is the shortest among known leafhoppers and lacks tandem repeats. Phylogenetic analysis of 13 protein-coding genes (PCGs), the first and second codons of 13 PCGs, 13 PCGs and two rRNAs formed three well-supported tree topologies. The topologies of phylogenetic trees inferred from three datasets were almost identical, which was consistent with previous molecular phylogenies of this group. Comparative morphological study of the ovipositors revealed several characters potentially useful for diagnosing genera and resolving their phylogenetic relationships. Phylogenetic analysis of these and other morphological characters yielded a tree that is mostly consistent with the tree obtained from analysis of mitogenome sequences. In both molecular and morphological phylogenenies, Typhlocybini and Zyginellini clustered into one clade, but neither was recovered as monophyletic.

Comparative mitochondrial genomes of four species of Sinopodisma and phylogenetic implications (Orthoptera, Melanoplinae)

In this study, the whole mitochondrial genomes (mitogenomes) from four species were sequenced. The complete mitochondrial genomes of Sinopodisma pieli, S. houshana, S. qinlingensis, and S. wulingshanensis are 15,857 bp, 15,818 bp, 15,843 bp, and 15,872 bp in size, respectively. The 13 protein-coding genes (PCGs) begin with typical ATN codons, except for COXI in S. qinlingensis, which begins with ACC. The highest A+T content in all the sequenced orthopteran mitogenomes is 76.8% (S. qinlingensis), followed by 76.5% (S. wulingshanensis), 76.4% (S. pieli) and 76.4% (S. houshana) (measured on the major strand). The long polythymine stretches (T-stretch) in the A+T-rich region of the four species are not adjacent to the trnI locus but are inside the stem-loop sequences on the major strand. Moreover, several repeated elements are found in the A+T-rich region of the four species. Phylogenetic analysis based on 53 mitochondrial genomes using Bayesian Inference (BI) and Maximum Likelihood (ML) revealed that Melanoplinae (Podismini) was a monophyletic group; however, the monophyly of Sinopodisma was not supported. These data will provide important information for a better understanding of the phylogenetic relationship of Melanoplinae.

Comparative analysis of mitogenomes among six species of grasshoppers (Orthoptera: Acridoidea: Catantopidae) and their phylogenetic implications in wing-type evolution

The degree of wing development has a close relationship with insects' movement ability and range, and it should also be closely related to mitochondrial-related genes. The complete mitochondrial genomes of six species of Catantopidae were sequenced, annotated and analyzed. Then, combined with 37 mitogenomes of grasshoppers, the ratio of nonsynonymous substitution to synonymous substitution (Ka/Ks) of the combined sequences of protein coding genes (PCGs) was calculated by DnaSP5, and the phylogenetic relationships were reconstructed by maximum likelihood (ML) and Bayesian (BI) methods based on PCGs+rRNAs. The results showed that the sizes of the six complete mitogenomes are Stenocatantops mistshenkoi Willemse F., 1968, 15,573 bp; Traulia lofaoshana Tinkham, 1940, 15,645 bp; Sinopodisma rostellocerca You, 1980, 15,622 bp; Anapodisma miramae Dovnar-Zapolskij, 1932, 15,189 bp; Qinlingacris elaeodes Yin & Chou, 1979, 15,221 bp; and Eozubovskya planicaudata Zhang & Jin, 1985, 15,830 bp; their structures are the same as those of Acridoidea. The AT bias of the wing-degenerated group (lobiform and apterous) is higher than that of the longipennate group, and more nonsynonymous substitutions accumulated in the wing-degenerated group than in the longipennate group (P = 0.000), which indicates that the wing-degenerated group has undergone weaker evolutionary selection than the longipinnate group. The phylogenetic tree shows that the wing-degenerated group in the Catantopidae are multiorigin and present parallel evolution.

Comparative mitochondrial genome analysis of Dendrolimus houi (Lepidoptera: Lasiocampidae) and phylogenetic relationship among Lasiocampidae species

Dendrolimus houi is one of the most common caterpillars infesting Gymnosperm trees, and widely distributed in several countries in Southeast Asia, and exists soley or coexists with several congeners and some Lasiocampidae species in various forest habitats. However, natural hybrids occasionally occur among some closely related species in the same habitat, and host preference, extreme climate stress, and geographic isolation probably lead to their uncertain taxonomic consensus. The mitochondrial DNA (mtDNA) of D. houi was extracted and sequenced by using high-throughput technology, and the mitogenome composition and characteristics were compared and analyzed of these species, then the phylogenetic relationship was constructed using the maximum likelihood method (ML) and the Bayesian method (BI) based on their 13 protein-coding genes (PCGs) dataset, which were combined and made available to download which were combined and made available to download among global Lasiocampidae species data. Mitogenome of D. houi was 15,373 bp in length, with 37 genes, including 13 PCGs, 22 tRNA genes (tRNAs) and 2 rRNA genes (rRNAs). The positions and sequences of genes were consistent with those of most known Lasiocampidae species. The nucleotide composition was highly A+T biased, accounting for ~80% of the whole mitogenome. All start codons of PCGs belonged to typical start codons ATN except for COI which used CGA, and most stop codons ended with standard TAA or TAG, while COI, COII, ND4 ended with incomplete T. Only tRNASer (AGN) lacked DHU arm, while the remainder formed a typical "clover-shaped" secondary structure. For Lasiocampidae species, their complete mitochondrial genomes ranged from 15,281 to 15,570 bp in length, and all first genes started from trnM in the same direction. And base composition was biased toward A and T. Finally, both two methods (ML and BI) separately revealed that the same phylogenetic relationship of D. spp. as ((((D. punctatus + D. tabulaeformis) + D. spectabilis) + D. superans) + (D. kikuchii of Hunan population + D. houi) as in previous research, but results were different in that D. kikuchii from a Yunnan population was included, indicating that different geographical populations of insects have differentiated. And the phylogenetic relationship among Lasiocampidae species was ((((Dendrolimus) + Kunugia) + Euthrix) + Trabala). This provides a better theoretical basis for Lasiocampidae evolution and classification for future research directions.

Complete Mitogenome of a Leaf-Mining Buprestid Beetle, Trachys auricollis, and Its Phylogenetic Implications

A complete mitogenome of Trachys auricollis is reported, and a mitogenome-based phylogenetic tree of Elateriformia with all protein-coding genes (PCGs), rRNAs, and tRNAs is presented for the first time. The complete mitochondrial genome of T. auricollis is 16,429 bp in size and contains 13 PCGs, two rRNA genes, 22 tRNA genes, and an A + T-rich region. The A + T content of the entire genome is approximately 71.1%, and the AT skew and GC skew are 0.10 and -0.20, respectively. According to the the nonsynonymous substitution rate to synonymous substitution rates (Ka/Ks) of all PCGs, the highest and lowest evolutionary rates were observed for atp8 and cox1, respectively, which is a common finding among animals. The start codons of all PCGs are the typical ATN. Ten PCGs have complete stop codons, but three have incomplete stop codons with T or TA. As calculated based on the relative synonymous codon usage (RSCU) values, UUA(L) is the codon with the highest frequency. Except for trnS1, all 22 tRNA genes exhibit typical cloverleaf structures. The A + T-rich region of T. auricollis is located between rrnS and the trnI-trnG-trnM gene cluster, with six 72-bp tandem repeats. Both maximum likelihood (ML) and Bayesian (BI) trees suggest that Buprestoidea is close to Byrrhoidea and that Buprestoidea and Byrrhoidea are sister groups of Elateroidea, but the position of Psephenidae is undetermined. The inclusion of tRNAs might help to resolve the phylogeny of Coleoptera.

Positive Selection Drove the Adaptation of Mitochondrial Genes to the Demands of Flight and High-Altitude Environments in Grasshoppers

The molecular evolution of mitochondrial genes responds to changes in energy requirements and to high altitude adaptation in animals, but this has not been fully explored in invertebrates. The evolution of atmospheric oxygen content from high to low necessarily affects the energy requirements of insect movement. We examined 13 mitochondrial protein-coding genes (PCGs) of grasshoppers to test whether the adaptive evolution of genes involved in energy metabolism occurs in changes in atmospheric oxygen content and high altitude adaptation. Our molecular evolutionary analysis of the 13 PCGs in 15 species of flying grasshoppers and 13 related flightless grasshoppers indicated that, similar to previous studies, flightless grasshoppers have experienced relaxed selection. We found evidence of significant positive selection in the genes ATP8, COX3, ND2, ND4, ND4L, ND5, and ND6 in flying lineages. This results suggested that episodic positive selection allowed the mitochondrial genes of flying grasshoppers to adapt to increased energy demands during the continuous reduction of atmospheric oxygen content. Our analysis of five grasshopper endemic to the Tibetan Plateau and 13 non-Tibetan grasshoppers indicated that, due to positive selection, more non-synonymous nucleotide substitutions accumulated in Tibetan grasshoppers than in non-Tibetan grasshoppers. We also found evidence for significant positive selection in the genes ATP6, ND2, ND3, ND4, and ND5 in Tibetan lineages. Our results thus strongly suggest that, in grasshoppers, positive selection drives mitochondrial genes to better adapt both to the energy requirements of flight and to the high altitude of the Tibetan Plateau.

Comparative analysis of mitochondrial genomes of the superfamily Grylloidea (Insecta, Orthoptera) reveals phylogenetic distribution of gene rearrangements

To further characterize mitochondrial genome (mitogenome) features of the superfamily Grylloidea (Insecta, Orthoptera), mitogenomes of Cacoplistes rogenhoferi and Meloimorpha japonica representing the family Mogoplistidae and three Ornebius species of Phalangopsidae were sequenced. A repeat-containing control region (CR) and 37 genes were present in these mitogenomes. Unusual start codons (TCG, CCG, and CTG) of cox1 and, in Ornebius, a partial stop codon (T) of nad1 followed by a 15-17-bp intergenic spacer were proposed based on transcript information and sequence alignments. The mitogenome-based phylogenetic trees suggest strongly the familial relationships as (((Phalangopsidae + Gryllidae) + Trigonidiidae) + Mogoplistidae). The exclusive occurrence of the trnE-trnS1-trnN rearrangement in Phalangopsidae, Gryllidae, and Trigonidiidae is suggestive of its appearance in the common ancestor of these families after the separation of Mogoplistidae. The trnV transposition in O. bimaculatus and formerly sequenced Trigonidium sjostedti (Trigonidiidae) indicates a potential consequence of parallel evolution. This study offers novel insights into mitogenome evolution, especially gene rearrangements, of Grylloidea.

Complete mitochondrial genome of bamboo grasshopper, Ceracris fasciata, and the phylogenetic analyses and divergence time estimation of Caelifera (Orthoptera)

The bamboo grasshopper Ceracris fasciata is regarded as a major pest species because of the damage it causes to bamboo, and its classification within the families and subfamilies of the suborder Caelifera remains unclear. Thus, we attempted to resolve these questions using molecular biology methods and analyses. Our results are as follows: (1) the complete mitochondrial genome of C. fasciata is 15,569 bp in length. The mitochondrial genome contains a standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region in the same order as those of the other analysed Caeliferan species. The putative start codon for the COX1 gene in C. fasciata is ACC, although it is not defined in other genes. The presence of tandem repeats of different sizes in the A + T-rich region may lead to size differences in other mitochondrial genomes. The mitochondrial genome of C. fasciata harbours the typical 37 genes and an A + T-rich region, and it shows similar characteristics to those of other grasshopper species. Characterization of the mitochondrial genome has enriched our knowledge of the mitochondrial genomes of Orthoptera around the world. Therefore, the phylogenetic relationships in Orthoptera can be re-examined. (2) In phylogenetic analyses, the monophyly of Orthoptera and its two suborders (Caelifera and Ensifera) has been consistently recovered based on most of the datasets selected, regardless of the optimal criteria. Our results do not support the monophyly of the subfamily Oedipodinae of Caelifera. We found that Phlaeoba albonema of the Acridinae is sorted into a clade with Ceracris in all our phylogenetic trees, and field experiments show that Phlaeoba always lives with Ceracris in the same ecotopes. Therefore, we suggest that Phlaeoba should be classified as a member of the Oedipodinae. We found that C. fasciata always clustered with Ceracris kiangsu, and both were sisters to Ceracris versicolor. Therefore, the genetic relationship between C. fasciata and C. kiangsu is closer than that between C. fasciata and C. versicolor. (3) The oldest estimated time of divergence of Ensifera in this context was determined to be 146.16 million years ago (Mya), or around the late Jurassic or early Cretaceous. We estimated that katydids (Grylloidea) likely diverged from other groups in the early Cretaceous. According to our divergence time analyses, we concluded that the ancestral Acrididae probably originated in the early Paleogene, and it is likely that the major diversification events happened at the middle Paleogene, well into the next geologic time. We estimated that crickets (Tettigoniidae) likely diverged from other groups in the early Cretaceous. Acrididae and Romaleinae group, Pyrgacrididae and Ommexechidae group, the youngest two clades we observed, were estimated to have diverged 58.79 Mya, between the middle and early Paleogene. C. versicolor is a sister to the group containing C. kiangsu and C. fasciata. First, C. versicolor diverged from the sister group (C. kiangsu + C. fasciata) around 44.81 Mya, and then the C. kiangsu and C. fasciata group separated at 43.04 Mya.

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.

Deep-level phylogeny of Cicadomorpha inferred from mitochondrial genomes sequenced by NGS

Recent development and advancement of next-generation sequencing (NGS) technologies have enabled the determination of mitochondrial genome (mitogenome) at extremely efficiency. In this study, complete or partial mitogenomes for 19 cicadomorphan species and six fulgoroid species were reconstructed by using the method of high-throughput sequencing from pooled DNA samples. Annotation analyses showed that the mitogenomes obtained have the typical insect mitogenomic content and structure. Combined with the existing hemipteran mitogenomes, a series of datasets with all 37 mitochondrial genes (up to 14,381 nt total) under different coding schemes were compiled to test previous hypotheses of deep-level phylogeny of Cicadomorpha. Thirty-seven species representing Cicadomorpha constituted the ingroup. A taxon sampling with nine species from Fulgoroidea and six from Heteroptera comprised the outgroup. The phylogenetic reconstructions congruently recovered the monophyly of each superfamily within Cicadomorpha. Furthermore, the hypothesis (Membracoidea + (Cicadoidea + Cercopoidea)) was strongly supported under the heterogeneous CAT model.

Characterization of the complete mitochondrial genome of Acanthoscelides obtectus (Coleoptera: Chrysomelidae: Bruchinae) with phylogenetic analysis

Acanthoscelides obtectus is a common species of the subfamily Bruchinae and a worldwide-distributed seed-feeding beetle. The complete mitochondrial genome of A. obtectus is 16,130 bp in length with an A + T content of 76.4%. It contains a positive AT skew and a negative GC skew. The mitogenome of A. obtectus contains 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a non-coding region (D-loop). All PCGs start with an ATN codon, and seven (ND3, ATP6, COIII, ND3, ND4L, ND6, and Cytb) of them terminate with TAA, while the remaining five (COI, COII, ND1, ND4, and ND5) terminate with a single T, ATP8 terminates with TGA. Except tRNA ^Ser , the secondary structures of 21 tRNAs that can be folded into a typical clover-leaf structure were identified. The secondary structures of lrRNA and srRNA were also predicted in this study. There are six domains with 48 helices in lrRNA and three domains with 32 helices in srRNA. The control region of A. obtectus is 1354 bp in size with the highest A + T content (83.5%) in a mitochondrial gene. Thirteen PCGs in 19 species have been used to infer their phylogenetic relationships. Our results show that A. obtectus belongs to the family Chrysomelidae (subfamily-Bruchinae). This is the first study on phylogenetic analyses involving the mitochondrial genes of A. obtectus and could provide basic data for future studies of mitochondrial genome diversities and the evolution of related insect lineages.

Nutritional imbalance suppresses migratory phenotypes of the Mongolian locust (Oedaleus asiaticus)

For many species, migration evolves to allow organisms to access better resources. However, the proximate factors that trigger these developmental changes, and how and why these vary across species, remain poorly understood. One prominent hypothesis is that poor-quality food promotes development of migratory phenotypes and this has been clearly shown for some polyphenic insects. In other animals, particularly long-distance bird migrants, it is clear that high-quality food is required to prepare animals for a successful migration. We tested the effect of diet quality on the flight behaviour and morphology of the Mongolian locust, Oedaleus asiaticus. Locusts reared at high population density and fed low-N grass (performance-enhancing for this species) had enhanced migratory morphology relative to locusts fed high-N grass. Furthermore, locusts fed synthetic diets with an optimal 1 : 2 protein : carbohydrate ratio flew for longer times than locusts fed diets with lower or higher protein : carbohydrate ratios. In contrast to the hypothesis that performance-degrading food should enhance migration, our results support the more nuanced hypothesis that high-quality diets promote development of migratory characteristics when migration is physiologically challenging.

Molecular phylogeny of Polyneoptera (Insecta) inferred from expanded mitogenomic data

The Polyneoptera represents one of the earliest insect radiations, comprising the majority of hemimetabolous orders, in which many species have great economic importance. Here, we sequenced eleven mitochondrial genomes of the polyneopteran insects by using high throughput pooled sequencing technology, and presented a phylogenetic reconstruction for this group based on expanded mitochondrial genome data. Our analyses included 189 taxa, of which 139 species represent all the major polyneopteran lineages. Multiple results support the monophyly of Polyneoptera, the monophyly of Dictyoptera, and the monophyly of Orthoptera. Sister taxon relationships Plecoptera + Dermaptera, and Zoraptera + Embioptera are also supported by most analyses. Within Dictyoptera, the Blattodea is consistently retrieved as paraphyly due to the sister group relationship of Cryptocercus with Isoptera. In addition, the results demonstrate that model selection, data treatment, and outgroup choice can have significant effects on the reconstructed phylogenetic relationships of Polyneoptera.

The complete mitochondrial genome of Calliptamus abbreviatus Ikovnnikov (Orthoptera: acridoidea)

Calliptamus abbreviatus Ikovnnikov (Orthoptera: Acridoidea) is one of the important pests in the grasslands in northern China. The complete mitochondrial genome of this insect was sequenced. This genome is 16,615 bp long, with an AT content of 73.3%, containing 37 typical animal mitochondrial genes and an AT-rich region. All 13 PCGs share the start codon ATN, and the usual termination codons (TAA) are found from 13 protein-coding genes. All of the 22 typical animal tRNA genes were found in C. abbreviatus mt-genome, and most of the tRNAs could be folded into the classic cloverleaf secondary structure except for tRNA-Ser (AGN), which lacks the dihydrouracil (DHU) stem. The sizes of the large and small ribosomal RNA genes are 1555 and 799 bp long, respectively. The AT content of the AT-rich region is 87.0%. Phylogenetic analysis supports that the coleopteran insects from the same family cluster in the same group.

Application of RNA-seq for mitogenome reconstruction, and reconsideration of long-branch artifacts in Hemiptera phylogeny

Hemiptera make up the largest nonholometabolan insect assemblage. Despite previous efforts to elucidate phylogeny within this group, relationships among the major sub-lineages remain uncertain. In particular, mitochondrial genome (mitogenome) data are still sparse for many important hemipteran insect groups. Recent mitogenomic analyses of Hemiptera have usually included no more than 50 species, with conflicting hypotheses presented. Here, we determined the nearly complete nucleotide sequence of the mitogenome for the aphid species of Rhopalosiphum padi using RNA-seq plus gap filling. The 15,205 bp mitogenome included all mitochondrial genes except for trnF. The mitogenome organization and size for R. padi are similar to previously reported aphid species. In addition, the phylogenetic relationships for Hemiptera were examined using a mitogenomic dataset which included sequences from 103 ingroup species and 19 outgroup species. Our results showed that the seven species representing the Aleyrodidae exhibit extremely long branches, and always cluster with long-branched outgroups. This lead to the failure of recovering a monophyletic Hemiptera in most analyses. The data treatment of Degen-coding for protein-coding genes and the site-heterogeneous CAT model show improved suppression of the long-branch effect. Under these conditions, the Sternorrhyncha was often recovered as the most basal clade in Hemiptera.

Complete Mitochondrial Genome of the Citrus Spiny Whitefly Aleurocanthus spiniferus (Quaintance) (Hemiptera: Aleyrodidae): Implications for the Phylogeny of Whiteflies

In this study, we sequenced the complete mitochondrial genome (15,220 bp) of the citrus spiny whitefly, Aleurocanthus spiniferus (Quaintance), a well-known pest from the superfamily Aleyrodidae. The A. spiniferus mitogenome contains 36 genes, including 13 protein-coding genes (PCGs), 21 transfer RNAs (tRNA), two ribosomal RNAs (rRNA) and a large non-coding region (control region, CR). Like most whiteflies, the A. spiniferus mitogenome had a large degree of rearrangement due to translocation of the nad3-trnG-cox3 gene cluster. The 13 PCGs initiated with ATN and generally terminated with TAA, although some used TAG or T as stop codons; atp6 showed the highest evolutionary rate, whereas cox2 appeared to have the lowest rate. The A. spiniferus mitogenome had 21 tRNAs with a typical cloverleaf secondary structure composed of four arms. Modeling of the two rRNA genes indicated that their secondary structure was similar to that of other insects. The CR of A. spiniferus was 920 bp and mapped between the nad3-trnG-cox3 and trnI-trnM gene clusters. One potential stem-loop structure and five tandem repeats were identified in the CR. Phylogenetic relationships of 11 species from the Aleyrodidae were analyzed based on the deduced amino acid sequences of the 13 PCGs and evolutionary characteristics were explored. Species with more genetic rearrangements were generally more evolved within the Aleyrodidae.

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.

The Complete Mitochondrial Genome and Song Evolution of the Monotypic Genus U. Tarbinsky, 1932 (Orthoptera: Tettigoniidae)

The insect Uvarovites inflatus Uvarov is highly appreciated in China. It is known for its distinctive songs and horn-like forewings and is raised commercially for insect lovers. U. inflatus was previously categorized as part of the monotypic genus Uvarovites; however, there was little molecular evidence to support this taxonomic classification. This study obtained and investigated the mitogenome of U. inflatus, and its songs were characterized and compared with other Ensifera species whose mitogenomes are available. By performing the mitochondrial analysis, we were able to assess the phylogenetic relationships between these species and discuss the evolution of Ensifera calling songs. The mitogenome of U. inflatus is 15,956 bp in length and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region. The organization and orientation of the U. inflatus mitogenome are similar to those of other Tettigonioidea species. Phylogenetic analysis based on 13 protein-coding genes showed that the superfamily Tettigonioidea is monophyletic, as are the other six tested subfamilies from Tettigonioidea. Our results also indicated that Grylloidea is monophyletic. A Bayesian relaxed clock analysis showed that the differentiation of U. inflatus and Gampsocleis gratiosa Brunner occurred in the middle Miocene, suggesting that their speciation occurred over a long evolutionary period. The results provide significant support for the establishment of the monotypic genus Uvarovites. Calling song analysis showed that at least two discrete steps of independent evolution occurred during the change from pure tone to broadband noise, and that the ancestor of existing Ensifera was more likely to have emitted pure-tone songs than broadband signals. Together, the mitogenome, molecular clock, and acoustic data allowed us to clarify the taxonomic state of U. inflatus and propose a timeline for the evolution of Ensifera songs.

The Complete Mitochondrial Genome of two Tetragnatha Spiders (Araneae: Tetragnathidae): Severe Truncation of tRNAs and Novel Gene Rearrangements in Araneae

Mitogenomes can provide information for phylogenetic analysis and evolutionary biology. The Araneae is one of the largest orders of Arachnida with great economic importance. In order to develop mitogenome data for this significant group, we determined the complete mitogenomes of two long jawed spiders Tetragnatha maxillosa and T. nitens and performed the comparative analysis with previously published spider mitogenomes. The circular mitogenomes are 14578 bp long with A+T content of 74.5% in T. maxillosa and 14639 bp long with A+T content of 74.3% in T. nitens, respectively. Both the mitogenomes contain a standard set of 37 genes and an A+T-rich region with the same gene orientation as the other spider mitogenomes, with the exception of the different gene order by the rearrangement of two tRNAs (trnW and trnG). Most of the tRNAs lose TΨC arm stems and have unpaired amino acid acceptor arms. As interesting features, both trnS(AGN) and trnS(UCN) lack the dihydrouracil (DHU) arm and long tandem repeat units are presented in the A+T-rich region of both the spider mitogenomes. The phylogenetic relationships of 23 spider mitogenomes based on the concatenated nucleotides sequences of 13 protein-coding genes indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Araneae. The molecular information acquired from the results of this study should be very useful for future researches on mitogenomic evolution and genetic diversities in spiders.

Phylogenetic relationships of Hemiptera inferred from mitochondrial and nuclear genes

Here, we reconstructed the Hemiptera phylogeny based on the expanded mitochondrial protein-coding genes and the nuclear 18S rRNA gene, separately. The differential rates of change across lineages may associate with long-branch attraction (LBA) effect and result in conflicting estimates of phylogeny from different types of data. To reduce the potential effects of systematic biases on inferences of topology, various data coding schemes, site removal method, and different algorithms were utilized in phylogenetic reconstruction. We show that the outgroups Phthiraptera, Thysanoptera, and the ingroup Sternorrhyncha share similar base composition, and exhibit "long branches" relative to other hemipterans. Thus, the long-branch attraction between these groups is suspected to cause the failure of recovering Hemiptera under the homogeneous model. In contrast, a monophyletic Hemiptera is supported when heterogeneous model is utilized in the analysis. Although higher level phylogenetic relationships within Hemiptera remain to be answered, consensus between analyses is beginning to converge on a stable phylogeny.

Extent and divergence of heteroplasmy of the DNA barcoding region in Anapodisma miramae (Orthoptera: Acrididae)

A partial sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene is widely used as a molecular marker for species identification in animals, also termed a DNA barcode. However, the presence of more than one sequence type in a single individual, also known as heteroplasmy, is one of the shortcomings of barcode identification. In this study, we examined the extent and divergence of COI heteroplasmy, including nuclear-encoded mitochondrial pseudogenes (NUMTs), at the genomic-DNA level from 13 insect species including orthopteran Anapodisma miramae, and a long fragment of mitochondrial DNA and cDNA from A. miramae as templates. When multiple numbers of clones originated from genomic DNA were sequenced, heteroplasmy was prevalent in all species and NUMTs were observed in five species. Long fragment DNA (∼13.5 kb) also is a source of heteroplasmic amplification, but the divergent haplotypes and NUMTs obtained from genomic DNA were not detected in A. miramae. On the other hand, cDNA was relatively heteroplasmy-free. Consistently, one dominant haplotype was always obtained from the genomic DNA-origin clones in all species and also from the long fragment- and cDNA-origin clones in the two tested individuals of A. miramae. Furthermore, the dominant haplotype was identical in sequence, regardless of the DNA source in A. miramae. Thus, one possible solution to avoid the barcoding problem in relationship to heteroplasmy could be the acquisition of multiple numbers of barcoding sequences to determine a dominant haplotype that can be assigned as barcoding sequence for a given species.

Novel odorant-binding proteins and their expression patterns in grasshopper, Oedaleus asiaticus

Insects use olfaction to detect exogenous odors and adapt to environments. In their olfaction systems, odorant-binding proteins (OBPs) are believed to be a key component. The unique OBP system of each species reflects the evolution of chemosensation of insects with habits. Here, we for the first time identified 15 OBPs, OasiOBP1-15, of a grasshopper, Oedaleus asiaticus, that lives in the grasslands of Northern China and is closely related to the locust, Locusta migratoria. OasiOBP9 and OasiOBP10 are specifically expressed in the antennae. Other OBPs are expressed in the antennae as well as other chemosensory organs, such as the mouthparts and wings. Significantly more OasiOBP7 was detected in male than female antennae, but there are 9 OBPs that were more expressed in female than male antennae by quantitative real-time PCR. Phylogenetic analysis indicated that most of the O. asiaticus OBPs are similar to those of L. migratoria, but some are substantially different. This indicates that the OBPs originally evolved in a common ancestor, but their unique chemosensory systems are adapted to different ecosystems.

The compact mitochondrial genome of Zorotypus medoensis provides insights into phylogenetic position of Zoraptera

BACKGROUND

Zoraptera, generally regarded as a member of Polyneoptera, represents one of the most enigmatic insect orders. Although phylogenetic analyses based on a wide array of morphological and/or nuclear data have been performed, the position of Zoraptera is still under debate. Mitochondrial genome (mitogenome) information is commonly considered to be preferable to reconstruct phylogenetic relationships, but no efforts have been made to incorporate it in Zorapteran phylogeny. To characterize Zoraptera mitogenome features and provide insights into its phylogenetic placement, here we sequenced, for the first time, one complete mitogenome of Zoraptera and reconstructed the phylogeny of Polyneoptera.

RESULTS

The mitogenome of Zorotypus medoensis with an A+T content of 72.50% is composed of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding A+T-rich region. The gene content and arrangement are identical to those considered ancestral for insects. This mitogenome shows a number of very unusual features. First, it is very compact, comprising 14,572 bp, and is the smallest among all known polyneopteran mitogenomes. Second, both noncoding sequences and coding genes exhibit a significant decrease in size compared with those of other polyneopterans. Third, Z. medoensis mitogenome has experienced an accelerated substitution rate. Fourth, truncated secondary structures of tRNA genes occur with loss of dihydrouridine (DHU) arm in trnC, trnR, and trnS(AGN) and loss of TΨC arm in trnH and trnT. The phylogenetic analyses based on the mitogenome sequence information indicate that Zoraptera, represented by Z. medoensis, is recovered as sister to Embioptera. However, both Zoraptera and Embioptera exhibit very long branches in phylogenetic trees.

CONCLUSIONS

Characterization of Z. medoensis mitogenome contributes to our understanding of the enigmatic Zoraptera. Mitogenome data demonstrate an overall strong resolution of deep-level phylogenies of Polyneoptera but not Insecta. It is preferable to expand taxon sampling of Zoraptera and other poorly represented orders in future to break up long branches.

The mitochondrial genome of Dastarcus helophoroides (Coleoptera: Bothrideridae) and related phylogenetic analyses

The complete mitochondrial genome of Dastarcus helophoroides (Coleoptera: Bothrideridae) which consists of 13 PCGs, 22 tRNA genes, two rRNA genes and a non-coding region (D-loop), is sequenced for its nucleotide sequence of 15,878 bp (GenBank: KF811054.1). The genome has a typical gene order which is identical to other Coleoptera species. Except for COI gene generally starts with non-canonical initial codon, all protein-coding genes start with ATN codon and terminate with the stop codon TA(A) or TAG. The secondary structure of rrnL and rrnS consists of 48 helices (contains four newly proposed helices) and 35 helices (contains two newly proposed helices) respectively. All 22 tRNAs in D. helophoroides are predicted to fold into typical cloverleaf secondary structure, except trnS1 (AGN), in which the dihydrouracil arm (DHU arm) could not form stable stem-loop structure. Thirteen protein-coding genes (nucleotide dataset and nucleic acid dataset) of the available species (29 taxa) have been used to infer the phylogenetic relationships among these orders. Tenebrionoidea and Cucujoidea form a sister group, and D. helophoroides is classified into Cucujoidea (Bothrideridae). The study first research on the phylogenetic analyses involving to the D. helophoroides mitogenome, and the results strongly bolster the current morphology-based hypothesis.

The complete mitochondrial genome of Angaracris rhodopa Fischer & Walheim (Orthoptera: Acridoidea)

Angaracris rhodopa Fischer & Walheim (Orthoptera: Acridoidea) is one of the important pests in the grasslands in northern China. The complete mitochondrial genome of this insect was sequenced. This genome is 15,930 bp long, with an AT content of 75.4%, containing 37 typical animal mitochondrial genes and an AT-rich region. All 13 PCGs share the start codon ATN, and the usual termination codons (TAA) are found from 13 protein-coding genes, except for ND2, COII, ND3 (T). All of the 22 typical animal tRNA genes were found in A. rhodopa mt-genome, and most of the tRNAs could be folded into the classic cloverleaf secondary structure except for tRNA-Ser (AGN), which lacks the dihydrouracil (DHU) stem. The sizes of the large and small ribosomal RNA genes are 1319 and 830 bp long, respectively. The AT content of the AT-rich region is 85.3%.

The complete mitochondrial genome of Angaracris barabensis Pallas (Orthoptera: Acridoidea)

Angaracris barabensis Pallas (Orthoptera: Acridoidea) is one of important pests in the grasslands in northern China. The complete mitochondrial genome of this insect was sequenced. This genome is 15,930 bp long, with an AT content of 75.5%, containing 37 typical animal mitochondrial genes and a AT-rich region. All genes were arranged in the same order as most of other Acridoidea. All 13 mitochondrial PCGs share the start codon ATN, and the usual termination codons (TAA) are found from 13 protein-coding genes, except for ND2, COII, ND3 (T). All of the 22 typical animal tRNA genes were found in A. barabensis mt-genome, and most of the tRNAs could be folded into the classic cloverleaf secondary structure except for tRNA-Ser (AGN), which lacks the dihydrouracil (DHU) stem. The sizes of the large and small ribosomal RNA genes are 1319 and 830 bp, respectively. The AT content of the AT-rich region is 85.6%.

Mitochondrial genomes of two Sinochlora species (Orthoptera): novel genome rearrangements and recognition sequence of replication origin

Background

Orthoptera, the largest polyneopteran insect order, contains 2 suborders and 235 subfamilies. Orthoptera mitochondrial genomes (mitogenomes) follow the ancestral insect gene order, with the exception of a trnD-trnK rearrangement in Acridomorphs and rare tRNA inversions. A question still remains regarding whether a long thymine-nucleotide stretch (T-stretch) involved in the recognition of the replication origin exists in the control region (CR) of Orthoptera mitochondrial DNA (mtDNA). Herein, we completed the sequencing of whole mitogenomes of two congeners (Sinochlora longifissa and S. retrolateralis), which possess overlapping distribution areas. Additionally, we performed comparative mitogenomic analysis to depict evolutionary trends of Orthoptera mitogenomes.

Results

Both Sinochlora mitogenomes possess 37 genes and one CR, a common gene orientation, normal structures of transfer RNA and ribosomal RNA genes, rather low A+T bias, and significant C skew in the majority strand (J-strand), resembling all the other sequenced ensiferans. Both mitogenomes are characterized by (1) a large size resulting from multiple copies of an approximately 175 bp GC-rich tandem repeat within CR; (2) a novel gene order (rrnS-trnI-trnM-nad2-CR-trnQ-trnW), compared to the ancestral order (rrnS-CR-trnI-trnQ-trnM-nad2-trnW); and (3) redundant trnS(UCN) pseudogenes located between trnS(UCN) and nad1. Multiple independent duplication events followed by random and/or non-random loss occurred during Sinochlora mtDNA evolution. The Orthoptera mtDNA recognition sequence of the replication origin may be one of two kinds: a long T-stretch situated in or adjacent to a possible stem-loop structure or a variant of a long T-stretch located within a potential stem-loop structure.

Conclusions

The unique Sinochlora mitogenomes reveal that the mtDNA architecture within Orthoptera is more variable than previously thought, enriching our knowledge on mitogenomic genetic diversities. The novel genome rearrangements shed light on mtDNA evolutionary patterns. The two kinds of recognition sequences of replication origin suggest that the regulatory sequences involved in the replication initiation process of mtDNA have diverged through Orthoptera evolution.

The mitochondrial genome of the quiet-calling katydids, Xizicus fascipes (Orthoptera: Tettigoniidae: Meconematinae)

To help determine whether the typical arthropod arrangement was a synapomorphy for the whole Tettigoniidae, we sequenced the mitochondrial genome (mitogenome) of the quiet-calling katydids, Xizicus fascipes (Orthoptera: Tettigoniidae: Meconematinae). The 16,166-bp nucleotide sequences of X. fascipes mitogenome contains the typical gene content, gene order, base composition, and codon usage found in arthropod mitogenomes. As a whole, the X. fascipes mitogenome contains a lower A+T content (70.2%) found in the complete orthopteran mitogenomes determined to date. All protein-coding genes started with a typical ATN codon. Ten of the 13 protein-coding genes have a complete termination codon, but the remaining three genes (COIII, ND5 and ND4) terminate with incomplete T. All tRNAs have the typical clover-leaf structure of mitogenome tRNA, except for tRNA(Ser(AGN)), in which lengthened anticodon stem (9 bp) with a bulged nuleotide in the middle, an unusual T-stem (6 bp in constrast to the normal 5 bp), a mini DHU arm (2 bp) and no connector nucleotides. In the A+T-rich region, two (TA)n conserved blocks that were previously described in Ensifera and two 150-bp tandem repeats plus a partial copy of the composed at 61 bp of the beginning were present. Phylogenetic analysis found: i) the monophyly of Conocephalinae was interrupted by Elimaea cheni from Phaneropterinae; and ii) Meconematinae was the most basal group among these five subfamilies.

A molecular phylogeny of Hemiptera inferred from mitochondrial genome sequences

Classically, Hemiptera is comprised of two suborders: Homoptera and Heteroptera. Homoptera includes Cicadomorpha, Fulgoromorpha and Sternorrhyncha. However, according to previous molecular phylogenetic studies based on 18S rDNA, Fulgoromorpha has a closer relationship to Heteroptera than to other hemipterans, leaving Homoptera as paraphyletic. Therefore, the position of Fulgoromorpha is important for studying phylogenetic structure of Hemiptera. We inferred the evolutionary affiliations of twenty-five superfamilies of Hemiptera using mitochondrial protein-coding genes and rRNAs. We sequenced three mitogenomes, from Pyrops candelaria, Lycorma delicatula and Ricania marginalis, representing two additional families in Fulgoromorpha. Pyrops and Lycorma are representatives of an additional major family Fulgoridae in Fulgoromorpha, whereas Ricania is a second representative of the highly derived clade Ricaniidae. The organization and size of these mitogenomes are similar to those of the sequenced fulgoroid species. Our consensus phylogeny of Hemiptera largely supported the relationships (((Fulgoromorpha,Sternorrhyncha),Cicadomorpha),Heteroptera), and thus supported the classic phylogeny of Hemiptera. Selection of optimal evolutionary models (exclusion and inclusion of two rRNA genes or of third codon positions of protein-coding genes) demonstrated that rapidly evolving and saturated sites should be removed from the analyses.

Insect phylogenomics: results, problems and the impact of matrix composition

In this study, we investigated the relationships among insect orders with a main focus on Polyneoptera (lower Neoptera: roaches, mantids, earwigs, grasshoppers, etc.), and Paraneoptera (thrips, lice, bugs in the wide sense). The relationships between and within these groups of insects are difficult to resolve because only few informative molecular and morphological characters are available. Here, we provide the first phylogenomic expressed sequence tags data ('EST': short sub-sequences from a c(opy) DNA sequence encoding for proteins) for stick insects (Phasmatodea) and webspinners (Embioptera) to complete published EST data. As recent EST datasets are characterized by a heterogeneous distribution of available genes across taxa, we use different rationales to optimize the data matrix composition. Our results suggest a monophyletic origin of Polyneoptera and Eumetabola (Paraneoptera + Holometabola). However, we identified artefacts of tree reconstruction (human louse Pediculus humanus assigned to Odonata (damselflies and dragonflies) or Holometabola (insects with a complete metamorphosis); mayfly genus Baetis nested within Neoptera), which were most probably rooted in a data matrix composition bias due to the inclusion of sequence data of entire proteomes. Until entire proteomes are available for each species in phylogenomic analyses, this potential pitfall should be carefully considered.

Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ∼2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.

Deep phylogeographic divergence and cytonuclear discordance in the grasshopper Oedaleus decorus

The grasshopper Oedaleus decorus is a thermophilic insect with a large, mostly south-Palaearctic distribution range, stretching from the Mediterranean regions in Europe to Central-Asia and China. In this study, we analyzed the extent of phylogenetic divergence and the recent evolutionary history of the species based on 274 specimens from 26 localities across the distribution range in Europe. Phylogenetic relationships were determined using sequences of two mitochondrial loci (ctr, ND2) with neighbour-joining and Bayesian methods. Additionally, genetic differentiation was analyzed based on mitochondrial DNA and 11 microsatellite markers using F-statistics, model-free multivariate and model-based Bayesian clustering approaches. Phylogenetic analyses detected consistently two highly divergent, allopatrically distributed lineages within O. decorus. The divergence among these Western and Eastern lineages meeting in the region of the Alps was similar to the divergence of each lineage to the sister species O. asiaticus. Genetic differentiation for ctr was extremely high between Western and Eastern grasshopper populations (F(ct)=0.95). Microsatellite markers detected much lower but nevertheless very significant genetic structure among population samples. The nuclear data also demonstrated a case of cytonuclear discordance because the affiliation with mitochondrial lineages was incongruent in Northern Italy. Taken together these results provide evidence of an ancient separation within Oedaleus and either historical introgression of mtDNA among lineages and/or ongoing sex-specific gene flow in this grasshopper. Our study stresses the importance of multilocus approaches for unravelling the history and status of taxa of uncertain evolutionary divergence.

Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust

The migratory locust, Locusta migratoria, is the most widely distributed grasshopper species in the world. However, its global genetic structure and phylogeographic relationships have not been investigated. In this study, we explored the worldwide genetic structure and phylogeography of the locust populations based on the sequence information of 65 complete mitochondrial genomes and three mitochondrial genes of 263 individuals from 53 sampling sites. Although this locust can migrate over long distances, our results revealed high genetic differentiation among the geographic populations. The populations can be divided into two different lineages: the Northern lineage, which includes individuals from the temperate regions of the Eurasian continent, and the Southern lineage, which includes individuals from Africa, southern Europe, the Arabian region, India, southern China, South-east Asia and Australia. An analysis of population genetic diversity indicated that the locust species originated from Africa. Ancestral populations likely separated into Northern and Southern lineages 895 000 years ago by vicariance events associated with Pleistocene glaciations. These two lineages evolved in allopatry and occupied their current distributions in the world via distinct southern and northern dispersal routes. Genetic differences, caused by the long-term independent diversification of the two lineages, along with other factors, such as geographic barriers and temperature limitations, may play important roles in maintaining the present phylogeographic patterns. Our phylogeographic evidence challenged the long-held view of multiple subspecies in the locust species and tentatively divided it into two subspecies, L. m. migratoria and L. m. migratorioides.

The complete mitochondrial genomes of two ghost moths, Thitarodes renzhiensis and Thitarodes yunnanensis: the ancestral gene arrangement in Lepidoptera

Background

Lepidoptera encompasses more than 160,000 described species that have been classified into 45–48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to the non-ditrysian lineage Hepialoidea and conducted a comparative mitogenomic analysis across Lepidoptera.

Results

The mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,816 bp long with an A + T content of 81.28 % and 82.34 %, respectively. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the A + T-rich region. Different tandem repeats in the A + T-rich region mainly account for the size difference between the two mitogenomes. All the protein-coding genes start with typical mitochondrial initiation codons, except for cox1 (CGA) and nad1 (TTG) in both mitogenomes. The anticodon of trnS(AGN) in T. renzhiensis and T. yunnanensis is UCU instead of the mostly used GCU in other sequenced Lepidoptera mitogenomes. The 1,584-bp sequence from rrnS to nad2 was also determined for an unspecified ghost moth (Thitarodes sp.), which has no repetitive sequence in the A + T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A + T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif ‘ATAGA’ and poly-T stretch in the A + T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes.

Conclusion

The mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event(s) likely occurred after Hepialoidea diverged from other lepidopteran lineages. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.

The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae)

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