A comparative analysis of mitochondrial genomes in Orthoptera (Arthropoda: Insecta) and genome descriptions of three grasshopper species

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.

The complete mitochondrial genome of a grasshopper endemic to the Qinghai-Tibet plateau, Uvaroviola multispinosa (Acrididae: Oedipodinae)

The grasshopper Uvaroviola multispinosa (Acrididae: Oedipodinae), is an agricultural pest to pasture and limitedly distributed in Qinghai-Tibet plateau of China. The complete mitochondrial genome of U. multispinosa is 15,620 bp long, which comprises of 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA (tRNA) genes and a putative non-coding control region (GenBank accession ID: MK829651). These genes are unequally distributed on different DNA chains, which 23 are located on the majority Chain. The nucleotide composition shows evidently bias (A, C, G, and T was 43.8, 14.8, 10.1, and 31.3%, respectively) with an overall AT content of 75.2%. All PCGs are initiated by ATN codons, among them ATG is the most preferred. Eleven PCGs use a common stop codon of TAA or TAG, whereas the remaining two were terminated with single T as incomplete stop codon. The phylogenetic relationships based on Bayes method showed that U. multispinosa is closely related to Compsorhipis davidiana, which is in accordance with its traditional morphological classification.

Characterization of the complete mitochondrial genome of the Yunnan endemic grasshopper Longchuanacris curvifurculus (Insecta: Orthoptera: Catantopidae)

Longchuanacris curvifurculus (L. curvifurculus) was once a dominating grasshopper in the Yunnan province (People's Republic of China) that occupy important ecological niche. However, its population has severely declined because of the deterioration of ecological environment. Identifying the species and source of L. curvifurculus is important for biodiversity conservation and ecological/environmental preservation. In the study, the complete mitochondrial genome of L. curvifurculus was assembled from high-coverage (36.8×) Illumina MiSeq sequencing data. The circular genome is 15,450 bp in length, harboring 37 typical mitochondrial genes and one control region. The nucleotide composition is asymmetric (43.0% A, 14.3% C, 10.5% G, and 32.2% T), with an overall A + T content of 75.2%. All the protein-coding genes (PCGs) are initiated with typical ATN start codons and terminated by the typical TAA codons or the incomplete T(aa) codon. The control region has a remarkably high A + T content (84.9%) and is located between genes rrnS and trnV.

Mitochondrial genomes of three Tetrigoidea species and phylogeny of Tetrigoidea

The mitochondrial genomes (mitogenomes) of Formosatettix qinlingensis, Coptotettix longjiangensis and Thoradonta obtusilobata (Orthoptera: Caelifera: Tetrigoidea) were sequenced in this study, and almost the entire mitogenomes of these species were determined. The mitogenome sequences obtained for the three species were 15,180, 14,495 and 14,538 bp in length, respectively, and each sequence included 13 protein-coding genes (PCGs), partial sequences of rRNA genes (rRNAs), tRNA genes (tRNAs) and a A + T-rich region. The order and orientation of the gene arrangement pattern were identical to that of most Tetrigoidea species. Some conserved spacer sequences between trnS(UCN) and nad1 were useful to identify Tetrigoidea and Acridoidea. The Ka/Ks value of atp8 between Trachytettix bufo and other four Tetrigoidea species indicated that some varied sites in this gene might be related with the evolution of T. bufo. The three Tetrigoidea species were compared with other Caelifera. At the superfamily level, conserved sequences were observed in intergenic spacers, which can be used for superfamily level identification between Tetrigoidea and Acridoidea. Furthermore, a phylogenomic analysis was conducted based on the concatenated data sets from mitogenome sequences of 24 species of Orthoptera in the superorders Caelifera and Ensifera. Both maximum likelihood and bayesian inference analyses strongly supported Acridoidea and Tetrigoidea as forming monophyletic groups. The relationships among six Tetrigoidea species were (((((Tetrix japonica, Alulatettix yunnanensis), Formosatettix qinlingensis), Coptotettix longjiangensis), Trachytettix bufo), Thoradonta obtusilobata).

The complete mitochondrial genome of Anopheles minimus (Diptera: Culicidae) and the phylogenetics of known Anopheles mitogenomes

Anopheles minimus is an important vector of human malaria in southern China and Southeast Asia. The phylogenetics of mosquitoes has not been well resolved, and the mitochondrial genome (mtgenome) has proven to be an important marker in the study of evolutionary biology. In this study, the complete mtgenome of An. minimus was sequenced for the first time. It is 15 395 bp long and encodes 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and a non-coding region. The gene organization is consistent with those of known Anopheles mtgenomes. The mtgenome performs a clear bias in nucleotide composition with a positive AT-skew and a negative GC-skew. All 13 PCGs prefer to use the codon UUA (Leu), ATN as initiation codon but cytochrome-oxidase subunit 1 (COI) and ND5, with TCG and GTG, and TAA as termination codon, but COI, COII, COIII and ND4, all with the incomplete T. tRNAs have the typical clover-leaf structure, but tRNA(Ser(AGN)) is consistent with known Anopheles mtgenomes. The control region includes a conserved T-stretch and a (TA)n stretch, and has the highest A+T content at 93.1%. The phylogenetics of An. minimus with 18 other Anopheles species was constructed by maximum likelihood and Bayesian inference, based on concatenated PCG sequences. The subgenera, Cellia and Anopheles, and Nyssorhynchus and Kerteszia have mutually close relationships, respectively. The Punctulatus group and Leucosphyrus group of Neomyzomyia Series, and the Albitarsis group of Albitarsis Series were suggested to be monophyletic. The monophyletic status of the subgenera, Cellia, Anopheles, Nyssorhynchus and Kerteszia need to be further elucidated.

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.

The complete mitochondrial genome of Mekongiella kingdoni (Uvarov, 1937) (Orthoptera: Acridoidea: Chrotogonidae)

The complete mitochondrial genome of Mekongiella kingdoni, which was collected from the Tibetan Plateau, is determined here. It is 15,932 bp in length and contains 73.1% AT. All M. kingdoni protein-coding sequences except for the cox1 start with a typical ATN codon. Instead, CCG, which is a rare but possible initiation codon, is located at the initiation context of cox1. The usual termination codons (TAA and TAG) and incomplete stop codons (TA or T) are found from 13 PCGs. All tRNA genes can be folded into the typical cloverleaf secondary structure, except trnS(AGN) due to the lacking of dihydrouridine arm. The sizes of the large and small ribosomal RNA genes are 1313 and 837 bp, respectively. The AT content of the A + T-rich region is 78.4%.

The phylogeny of the Orthoptera (Insecta) as deduced from mitogenomic gene sequences

Background

The phylogeny of the Orthoptera was analyzed based on 6 datasets from 47 orthopteran mitochondrial genomes (mitogenomes). The phylogenetic signals in the mitogenomes were rigorously examined under analytical regimens of maximum likelihood (ML) and Bayesian inference (BI), along with how gene types and different partitioning schemes influenced the phylogenetic reconstruction within the Orthoptera. The monophyly of the Orthoptera and its two suborders (Caelifera and Ensifera) was consistently recovered in the analyses based on most of the datasets we selected, regardless of the optimality criteria.

Results

When the seven NADH dehydrogenase subunits were concatenated into a single alignment (NADH) and were analyzed; a near-identical topology to the traditional morphological analysis was recovered, especially for BI_NADH. In both the concatenated cytochrome oxidase (COX) subunits and COX + cytochrome b (Cyt b) datasets, the small extent of sequence divergence seemed to be helpful for resolving relationships among major Orthoptera lineages (between suborders or among superfamilies). The conserved and variable domains of ribosomal (r)RNAs performed poorly when respectively analyzed but provided signals at some taxonomic levels.

Conclusions

Our findings suggest that the best phylogenetic inferences can be made when moderately divergent nucleotide data from mitogenomes are analyzed, and that the NADH dataset was suited for studying orthopteran phylogenetic relationships at different taxonomic levels, which may have been due to the larger amount of DNA sequence data and the larger number of phylogenetically informative sites.

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.

The complete mitochondrial genomes of three grasshoppers, Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis (Orthoptera: Pamphagidae)

The complete mitogenomes of Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis are 15,660 bp, 15,657 bp and 15,661 bp in size, respectively. All three mitogenomes contain a standard set of 13 protein - coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A+T-rich region in the same order as those of the other analysed caeliferan species, including the rearrangement of trnAsp and trnLys. The putative initiation codon for the cox1 gene in the three species is CCG. The long polythymine stretch (T-stretch) in the A+T-rich region of the three species is not adjacent to the trnIle but inside the stem-loop sequence in the majority strand. The mitogenomes of F. helanshanensis and P. rubimarginis have higher overall similarities. The characterization of the three mitogenomes will enrich our knowledge on the Pamphagidae mitogenome. The phylogenetic analyses indicated that within the Caelifera, Pyrgomorphoidea is a sister group to Acridoidea. The species from the Pamphagidae form a monophyletic group, as is the case for Acrididae. Furthermore, the two families cluster as sister groups, supporting the monophyly of Acridoidea. The relationships among eight acridid subfamilies were (Cyrtacanthacridinae+(Calliptaminae+(Catantopinae+(Oxyinae+(Melanopline+(Acridinae+(Oedipodinae+Gomphocerinae).

Complete mitochondrial genome of Locusta migratoria migratoria (Orthoptera: Oedipodidae): three tRNA-like sequences on the N-strand

The complete 16053 bp mitochondrial genome (mitogenome) sequence of Locusta migratoria migratoria has been determined. This mitogenome contains the base compositional biases and codon usage typical of metazoans, and the RSCU values indicate a negative correlation with the C and G contents in codon. The orientation and gene order of the L. migratoria migratoria is identical to Locusta migratoria migratoiodes. An unusual feature of the L. migratoria migratoria mitogenome is the presence of three tRNA-like structures on the N-strand: one tRNA(Ile)-like and two tRNA(Leu(CUN))-like sequences. The tRNA-like sequences have proper folding structures and anticodons sequences. Two repeated DNA sequences, Rpt I and Rpt II, were found in the A+T-rich region of the L. migratoria migratoria mitogenome. Both repeated sequences have various features. In the 5' region of Rpt I, a 51 bp fragment is localized in the srRNA gene; and there are two tandemly sub-repeated DNA sequences (sub-Rpts), Rpt 1-4, within Rpt I and Rpt II. One stem-loop structure on the N-strand that may be involved in the N-strand replication initiation was found in the A+T-rich region.