Analysis of the mitochondrial genome of Schistocerca gregaria gregaria (Orthoptera: Acrididae)

Phylomitogenomics of two Neotropical species of long-legged crickets Endecous Saussure, 1878 (Orthoptera: Phalangopsidae)

Mitochondrial genomes have provided significant insights into the evolution of several insects. A typical mitogenome contains 37 genes, and variations in gene order can indicate evolutionary relationships between species. In this study, we have assembled the first complete mitogenomes of Endecous chape and E. onthophagus and analyzed the phylogenetic implications for the Gryllidea infraorder. We performed DNA extractions and genome sequencing for both Endecous species. Subsequently, we searched for raw data in the Sequence Read Archive (SRA) in NCBI. Using the SRA data, we assembled the partial mitogenome of Dianemobius nigrofasciatus and annotated the protein-coding genes (PCGs) for nine species. Phylogenomic relationships were reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI), utilizing the PCGs from 49 Gryllidea species. The mitogenome lengths of E. chape and E. onthophagus are 16,266 bp and 16,023 bp, respectively, while D. nigrofasciatus has a length of 15,359 bp. Our results indicate that species within the infraorder exhibit four types of gene order arrangements that align with their phylogenetic relationships. Both phylogenomic trees displayed strong support, and the ML corroborated with the literature. Gryllidea species have significantly contributed to various fields, and studying their mitogenomes can provide valuable insights into this infraorder evolution.

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.

Understanding the molecular evolution of tiger diversity through DNA barcoding marker ND4 and NADH dehydrogenase complex using computational biology

BACKGROUND

Currently, Tigers (the top predator of an ecosystem) are on the list of endangered species. Thus the need is to understand the tiger's population genomics to design their conservation strategies.

OBJECTIVE

We analyzed the molecular evolution of tiger diversity using NADH dehydrogenase subunit 4 (ND4), a significant electron transport chain component.

METHODS

We have analyzed nucleotide composition and distribution pattern of ND genes, molecular evolution, evolutionary conservation pattern and conserved blocks of NADH, phylogenomics of ND4, and estimating species divergence, etc., using different bioinformatics tools and software, and MATLAB programming and computing environment.

RESULTS

The nucleotide composition and distribution pattern of ND genes in the tiger genome demonstrated an increase in the number of adenine (A) and a lower trend of A+T content in some place of the distribution analysis. However, the observed distributions were not significant (P > 0.05). Evolutionary conservation analysis showed three highly align blocks (186 to 198, 406 to 416, and 527 to 545). On mapping the molecular evolution of ND4 among model species (n = 30), we observed its presence in a broader range of species. ND4 based molecular evolution of tiger diversity and time divergence for a tiger (20 different other species) shows that genus Panthera originated more or less at a similar time.

CONCLUSIONS

The nucleotide composition and nucleotide distribution pattern of tiger ND genes showed the evolutionary pattern and origin of tiger and Panthera lineage concerning the molecular clock, which will help to understand their adaptive evolution.

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.

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.

Characteristics of the complete mitochondrial genome of Suhpalacsa longialata (Neuroptera, Ascalaphidae) and its phylogenetic implications

The owlflies (Family Ascalaphidae) belong to the Neuroptera but are often mistaken as dragonflies because of morphological characters. To date, only three mitochondrial genomes of Ascalaphidae, namely Libelloides macaronius; Ascaloptynx appendiculatus; Ascalohybris subjacens, are published in GenBank, meaning that they are greatly under-represented in comparison with the 430 described species reported in this family. In this study, we sequenced and described the complete mitochondrial genome of Suhpalacsa longialata (Neuroptera, Ascalaphidae). The total length of the S. longialata mitogenome was 15,911 bp, which is the longest known to date among the available family members of Ascalaphidae. However, the size of each gene was similar to the other three Ascalaphidae species. The S. longialata mitogenome included a transposition of tRNACys and tRNATrp genes and formed an unusual gene arrangement tRNACys-tRNATrp-tRNATyr (CWY). It is likely that the transposition occurred by a duplication of both genes followed by random loss of partial duplicated genes. The nucleotide composition of the S. longialata mitogenome was as follows: A = 41.0%, T = 33.8%, C = 15.5%, G = 9.7%. Both Bayesian inference and ML analyses strongly supported S. longialata as a sister clade to (Ascalohybris subjacens + L. macaronius), and indicated that Ascalaphidae is not monophyletic.

Characterization of the mitochondrial genome of the montane grasshopper, Qinlingacris elaeodes (Orthoptera: Catantopidae)

Qinlingacris elaeodes is the dominant grasshopper at an altitude of 3000 meters and above, and is a representative species of the genus Qinlingacris endemic to China. The sequenced mitochondrial genome of this grasshopper is 14,818 bp in length, including 13 protein-coding genes (ND1-6, COI-III, ATP6, ATP8, ND4L, CTYB), 21 transfer RNAs, and 2 ribosomal RNAs (12S and 16S). The orientation and gene order of these genes are identical to those found in the putative ancestral insect mitogenome. The 13 PCGs start with a typical ATN codon as their start codons. The usual TAA and TAG termination codons are found for 12 PCGs. However, the ND5 gene has an incomplete termination codon (T).

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.

Ancestral gene organization in the mitochondrial genome of Thyridosmylus langii (McLachlan, 1870) (Neuroptera: Osmylidae) and implications for lacewing evolution

The first complete mitochondrial genome of the lacewing family Osmylidae (Thyridosmylus langii (McLachlan, 1870)) (Neuroptera) was sequenced in this study. The genome is a circular molecule of 16,221 bp containing the typical 37 genes but is arranged in the same order as that of the putative ancestor of hexapod and lacks translocation of trnC as shared by all previously sequenced neuropteran mtDNAs. This reveals that trnC translocation does not represent an organizational synapomorphy in the mitochondrion for the entire Neuroptera clade. Comparative analysis of neuropteran tRNA genes reveals a relatively slow and conserved evolution of the mitochondrion throughout the order. Secondary structure models of the ribosomal RNA genes of T. langii largely agree with those proposed for other insect orders. Nevertheless, domain I of T. langii rrnL is consisted of nine helices rather than eight helices which is typical for neuropteran rrnL. Protein-coding genes have typical mitochondrial start codons, with the exception of COI, which uses the TCG start codon also found in Ithonidae and Chrysopidae. Like other neuropteran insects, the control region is the most AT-rich region and comparatively simple, with little evidence of conserved blocks or long tandem repeats. Considering the issues of base-compositional and branch length heterogeneity, we used a range of phylogenetic approaches to recover neuropteridan relationships and explored the effect of method choice on recovery of monophyly of Neuropterida: ((Neuroptera + Megaloptera) + Raphidioptera). The monophyly of Neuroptera and the more basal position of Osmylidae were also recovered by different datasets and phylogenetic methods.

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 Sequence of Acrida cinerea (Acrididae: Orthoptera) and Comparative Analysis of Mitochondrial Genomes in Orthoptera

The complete 15,599-bp mitogenome of Acrida cinerea was determined and compared with that of the other 20 orthopterans. It displays characteristic gene content, genome organization, nucleotide composition, and codon usage found in other Caelifera mitogenomes. Comparison of 21 orthopteran sequences revealed that the tRNAs encoded by the H-strand appear more conserved than those by the L-stand. All tRNAs form the typical clover-leaf structure except trnS (agn), and most of the size variation among tRNAs stemmed from the length variation in the arm and loop of TΨC and the loop of DHU. The derived secondary structure models of the rrnS and rrnL from 21 orthoptera species closely resemble those from other insects on CRW except a considerably enlarged loop of helix 1399 of rrnS in Caelifera, which is a potentially autapomorphy of Caelifera. In the A+T-rich region, tandem repeats are not only conserved in the closely related mitogenome but also share some conserved motifs in the same subfamily. A stem-loop structure, 16 bp or longer, is likely to be involved in replication initiation in Caelifera and Grylloidea. A long T-stretch (>17 bp) with conserved stem-loop structure next to rrnS on the H-strand, bounded by a purine at either end, exists in the three species from Tettigoniidae.