Sequence and phylogenetic analysis of complete mitochondrial DNA genomes of two grasshopper speciesGomphocerus rufus(Linnaeus, 1758) andPrimnoa arctica(Zhang and Jin, 1985) (Orthoptera: Acridoidea)

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.

The Complete Mitogenomes of Three Grasshopper Species with Special Notes on the Phylogenetic Positions of Some Related Genera

Clarifying phylogenetic position and reconstructing robust phylogeny of groups using various evidences are an eternal theme for taxonomy and systematics. In this study, the complete mitogenomes of Longzhouacris mirabilis, Ranacris albicornis, and Conophyma zhaosuensis were sequenced using next-generation sequencing (NGS), and the characteristics of the mitogenomes are presented briefly. The mitogenomes of the three species are all circular molecules with total lengths of 16,164 bp, 15,720 bp, and 16,190 bp, respectively. The gene structures and orders, as well as the characteristics of the mitogenomes, are similar to those of other published mitogenomes in Caelifera. The phylogeny of the main subfamilies of Acrididae with prosternal process was reconstructed using a selected dataset of mitogenome sequences under maximum likelihood (ML) and Bayesian inference (BI) frameworks. The results showed that the genus Emeiacris consistently fell into the subfamily Melanoplinae rather than Oxyinae, and the genus Choroedocus had the closest relationship with Shirackiacris of the subfamily Eyprepocnemidinae in both phylogenetic trees deduced from mitogenome protein coding genes (PCGs). This finding is entirely consistent with the morphological characters, which indicate that Emeiacris belongs to Melanoplinae and Choroedocus belongs to Eyprepocnemidinae. In addition, the genera Conophymacris and Xiangelilacris, as well as Ranacris and Menglacris, are two pairs of the closest relatives, but their phylogenetic positions need further study to clarify.

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

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

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.

Characterization and analysis of full-length transcriptomes from two grasshoppers, Gomphocerus licenti and Mongolotettix japonicus

Acrididae are diverse in size, body shape, behavior, ecology and life history; widely distributed; easy to collect; and important to agriculture. They represent promising model candidates for functional genomics, but their extremely large genomes have hindered this research; establishing a reference transcriptome for a species is the primary means of obtaining genetic information. Here, two Acrididae species, Gomphocerus licenti and Mongolotettix japonicus, were selected for full-length (FL) PacBio transcriptome sequencing. For G. licenti and M. japonicus, respectively, 590,112 and 566,165 circular consensus sequences (CCS) were generated, which identified 458,131 and 428,979 full-length nonchimeric (FLNC) reads. After isoform-level clustering, next-generation sequencing (NGS) short sequences were used for error correction, and remove redundant sequences with CD-HIT, 17,970 and 16,766 unigenes were generated for G. licenti and M. japonicus. In addition, we obtained 17,495 and 16,373 coding sequences, 1,082 and 813 transcription factors, 11,840 and 10,814 simple sequence repeats, and 905 and 706 long noncoding RNAs by analyzing the transcriptomes of G. licenti and M. japonicus, respectively, and 15,803 and 14,846 unigenes were annotated in eight functional databases. This is the first study to sequence FL transcriptomes of G. licenti and M. japonicus, providing valuable genetic resources for further functional genomics research.

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.

DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety-three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.

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 Kingdonella bicollina (Orthoptera: Acridoidea: Catantopidae)

The complete mitochondrial genome of Kingdonella bicollina, which was collected from the Tibetan Plateau, is determined here. It is 15,630 bp in length and contains 75.6% AT. All 13 mitochondrial PCGs share the start codon ATN, except for cox1 (CCG start codon). The termination codons include TAA, TAG and T. All tRNA genes could be folded into the typical cloverleaf secondary structure, except trnS(AGN) lacking of dihydrouridine arm. The sizes of the large and small ribosomal RNA genes are 1312 and 859 bp, respectively. The A+T-rich region is 739 bp in length contains 79.7% AT content.

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 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).