Evolutionary rates of and selective constraints on the mitochondrial genomes of Orthoptera insects with different wing types

The Fast Evolution of the Stenobothrini Grasshoppers (Orthoptera, Acrididae, and Gomphocerinae) Revealed by an Analysis of the Control Region of mtDNA, with an Emphasis on the Stenobothrus eurasius Group

The two cryptic grasshopper species of the genus Stenobothrus, S. eurasius and S. hyalosuperficies, demonstrate different acoustic behavior despite a strong similarity in morphology. A hybridization between these species is possible in the contact zone; however, there are little molecular data about the relationships of these species. The analysis of the mtDNA control region (CR) reveals that haplotypes of S. hyalosuperficies have more in common with the more distant Stenobothrus species than with the closely related S. eurasius. In the contact zone, S. eurasius has mt-haplotypes shared with S. hyalosuperficies, which might indicate an introgression of mtDNA from S. hyalosuperficies to the S. eurasius gene pool. We also analyze the structure and evolutionary rate of the mtDNA CR for the Stenobothrus genus and estimate the time of divergence of the species within the genus. The phylogenetic tree of the tribe Stenobothrini reconstructed with either the CR or COI gave the same four groups. The phylogenetic tree of the Stenobothrus genus has a star-like topology with each mtDNA haplotype found in any analyzed species, except S. eurasius, which forms a separate branch. The maximum degree of incomplete lineage sorting can demonstrate either ancestral polymorphism or introgression.

Repetitive element expansions contribute to genome size gigantism in Pamphagidae: A comparative study (Orthoptera, Acridoidea)

Pamphagidae is a family of Acridoidea that inhabits the desert steppes of Eurasia and Africa. This study employed flow cytometry to estimate the genome size of eight species in the Pamphagidae. The results indicate that the genome size of the eight species ranged from 13.88 pg to 14.66 pg, with an average of 14.26 pg. This is the largest average genome size recorded for the Orthoptera families, as well as for the entire Insecta. Furthermore, the study explored the role of repetitive sequences in the genome, including their evolutionary dynamics and activity, using low-coverage next-generation sequencing data. The genome is composed of 14 different types of repetitive sequences, which collectively make up between 59.9% and 68.17% of the total genome. The Pamphagidae family displays high levels of transposable element (TE) activity, with the number of TEs increasing and accumulating since the family's emergence. The study found that the types of repetitive sequences contributing to the TE outburst events are similar across species. Additionally, the study identified unique repetitive elements for each species. The differences in repetitive sequences among the eight Pamphagidae species correspond to their phylogenetic relationships. The study sheds new light on genome gigantism in the Pamphagidae and provides insight into the correlation between genome size and repetitive sequences within the family.

Comparative Analysis of the Mitochondrial Genomes of Three Species of Yangiella (Hemiptera: Aradidae) and the Phylogenetic Implications of Aradidae

The mitochondrial genomes of three species of Yangiella were sequenced, annotated, and analyzed. The genome length of the three species of the genus is 15,070-15,202 bp, with a typical gene number, including a control region, 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). It was found that the mitochondrial genome of Yangiella had AT bias. Except for the lack of a DHU arm of the trnS1 gene, the other tRNAs had a typical cloverleaf structure, and the codon usage preferences of the three species exhibited high similarity. In addition, tRNA gene rearrangements were observed among the three subfamilies of Aradidae (Mezirinae, Calisiinae, Aradinae), and it was found that codon usage preferences appeared to be less affected by base mutation and more by natural selection. The Pi and Ka/Ks values indicated that cox1 was the most conserved gene in the mitochondrial genome of Aradidae, while atp8 and nad6 were rapidly evolved genes. Substitution saturation level analysis showed that the nucleic acid sequence of mitochondrial protein-coding genes in Aradidae did not reach saturation, suggesting the rationality of the phylogenetic analysis data. Bayesian and maximum likelihood methods were used to analyze the phylogeny of 16 species of Hemiptera insects, which supported the monophyly of Aneurinae, Carventinae, and Mezirinae, as well as the monophyly of Yangiella. Based on fossils and previous studies, the differentiation time was inferred, indicating that Yangiella diverged about 57 million years ago.

Comparative Mitogenomics and Phylogenetic Implications for Nine Species of the Subfamily Meconematinae (Orthoptera: Tettigoniidae)

Currently, the subfamily Meconematinae encompasses 1029 species, but whole-mitochondrial-genome assemblies have only been made available for 13. In this study, the whole mitochondrial genomes (mitogenomes) of nine additional species in the subfamily Meconematinae were sequenced. The size ranged from 15,627 bp to 17,461 bp, indicating double-stranded circular structures. The length of the control region was the main cause of the difference in mitochondrial genome length among the nine species. All the mitogenomes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (CR). The majority strand encoded 23 genes, and the minority strand encoded 14 genes. A phylogenetic analysis reaffirmed the monophyletic status of each subfamily, but the monophysitism of Xizicus, Xiphidiopsis and Phlugiolopsis was not supported.

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.

The complete mitochondrial genome of Phymateus saxosus (Coquerel, 1861) (Orthoptera: Pyrgomorphidae) and phylogenetic analysis

Phymateus saxosus is a member of the family Pyrgomorphidae, Orthoptera. In this study, the complete mitochondrial genome (mitogenome) of P. saxosus was determined and analyzed. Assembled mitogenome sequence of P. saxosus is 15,672 bp in size, containing 37 genes and a control region. The gene orientation and arrangement of P. saxosus are identical to other species in the Pyrgomorphoidea family. The overall nucleotide composition is as follows: A (43.6%) > T (30.2%) > C (16.1%) > G (10.1%). Phylogenetic analysis suggested that P. saxosus forms sister groups with P. morbillosus, and the monophyly of Pyrgomorphidae is supported. In general, this study provided valuable genetic information for P. saxosus and explored the phylogenetic relationships in the family Pyrgomorphidae.

Phylogeny and biogeography of the wingless orthopteran family Rhaphidophoridae

Cave crickets (Rhaphidophoridae) are insects of an ancient and wingless lineage within Orthoptera that are distributed worldwide except in Antarctica, and each subfamily has a high level of endemicity. Here, we show the comprehensive phylogeny of cave crickets using multi-gene datasets from mitochondrial and nuclear loci, including all extant subfamilies for the first time. We reveal phylogenetic relationships between subfamilies, including the sister relationship between Anoplophilinae and Gammarotettiginae, based on which we suggest new synapomorphies. Through biogeographic analyses based on divergence time estimations and ancestral range reconstruction, we propose novel hypotheses regarding the biogeographic history of cave crickets. We suggest that Gammarotettiginae in California originated from the Asian lineage when Asia and the Americas were connected by the Bering land bridge, and the opening of the western interior seaway affected the division of Ceuthophilinae from Tropidischiinae in North America. We estimate that Rhaphidophoridae originated at 138 Mya throughout Pangea. We further hypothesize that the loss of wings in Rhaphidophoridae could be the result of their adaptation to low temperatures in the Mesozoic era.

Multiple independent origins of duplicated mitochondrial control regions indicate an apomorphy in the Thysanoptera (Insecta)

The mitochondrial genome (mitogenome) of thrips is characterized by the presence of control region (CR) duplication. However, the evolution pattern of duplicated CRs in thrips is still unclear. In this study, the multiple independent origins of duplicated CR indicated that the CR duplication was not an ancestral state for Thysanoptera. The macroevolutionary pattern suggested that the earliest CR duplication event occurred in the middle Cretaceous (94.85 Ma) coincided with rearrangement events forming the ancestors of Aeolothripidae, but much later than that forming the ancestors of the suborder Terebrantia. The mitogenome with duplicated CRs showed a higher rate of gene rearrangement. The sequence similarity of the CR copies and divergence time were negatively correlated, indicating age-related deterioration of mitochondrial function. No significant differences were found in the mitochondrial DNA, the P123 and P4FD between the single and multiple-CR charactered mitogenomes, which suggested that the duplicated CRs may not affect the replication process in thrip mitogenome. The mitogenomes with duplicated CRs (mean: 0.0088 subs/s/my) show a significantly increased evolutionary rate than that with a single one (mean: 0.0058 subs/s/my). However, it seems that this higher evolutionary rate did not have adaptive mechanisms in Terebrantia. We speculated that the duplicated CRs may cause a more intense production of energy by mitochondria, and an accelerated mutation and substitution rate is expected in such mitogenomes. Our study provided new insights into the presence of CR duplications and their evolution in the mitogenomes of thrips.

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.

Adaptive evidence of mitochondrial genes in Pteromalidae and Eulophidae (Hymenoptera: Chalcidoidea)

Pteromalidae and Eulophidae are predominant and abundant taxa within Chalcidoidea (Hymenoptera: Apocrita). These taxa are found in diverse ecosystems, ranging from basin deserts (200 m) to alpine grasslands (4500 m). Mitochondria, cellular powerhouses responsible for energy production via oxidative phosphorylation, are sensitive to various environmental factors such as extreme cold, hypoxia, and intense ultraviolet radiation characteristic of alpine regions. Whether the molecular evolution of mitochondrial genes in these parasitoids corresponds to changes in the energy requirements and alpine environmental adaptations remains unknown. In this study, we performed a comparative analysis of mitochondrial protein-coding genes from 11 alpine species of Pteromalidae and Eulophidae, along with 18 lowland relatives, including 16 newly sequenced species. We further examined the codon usage preferences (RSCU, ENC-GC3s, neutrality, and PR2 bias plot) in these mitochondrial protein-coding sequences and conducted positive selection analysis based on their Bayesian phylogenetic relationships, and identified positive selection sites in the ATP6, ATP8, COX1, COX3, and CYTB genes, emphasizing the crucial role of mitochondrial gene adaptive evolution in the adaptation of Pteromalidae and Eulophidae to alpine environments. The phylogenetically independent contrast (PIC) analysis results verified the ω ratio of 13 PCGs from Pteromalidae and Eulophidae increased with elevation, and results from generalized linear model confirm that ATP6, ATP8, COX3, and ND1 are closely correlated with temperature-related environmental factors. This research not only enriched the molecular data of endemic alpine species but also underscores the significance of mitochondrial genes in facilitating the adaptation of these minor parasitoids to plateau habitats.

Mitogenome of Xya pfaendleri (Orthoptera: Caelifera): Its comparative description and phylogenetic position within Tridactylidea

We report the comparative examination of the complete mitochondrial genome of the pygmy mole cricket Xya pfaendleri (Orthoptera: Caelifera: Tridactylidae). The mitogenome consists of 13 protein-coding regions, 22 tRNAs, two rRNAs, and one control region, following the gene order of the ancestral pancrustacean mitogenome. The length of the mitogenome in Xya pfaendleri is 15352 bp. The start and stop codons of the protein-coding genes exhibit the general pattern observed in orthopterans. The data indicate that the pattern of gene overlapping/intergenic sequences exhibits a significant phylogenetic signal. A phylogenetic tree inferred using 12 mitogenomes (seven belonging to Tridactylidea, three to Acrididea, and two to Ensifera) confirms the sister group relationship of Acrididea and Tridactylidea. The relationship among the families of Tridactylidea is Cylindrachetidae + (Ripipterygidae + Tridactylidae). The mitogenome sequences of Xya and Tridactylus constitute a single clade, sharing a last common ancestor 94 million years ago, and rendering the first genus paraphyletic. The present preliminary data suggest that we still have much to learn about the evolution and diversity of Tridactylidea.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Two New Mitogenomes of Bibionidae and Their Comparison within the Infraorder Bibionomorpha (Diptera)

Despite the worldwide distribution and rich diversity of the infraorder Bibionomorpha in Diptera, the characteristics of mitochondrial genomes (mitogenomes) are still little-known, and the phylogenetics and evolution of the infraorder remains controversial. In the present study, we report complete and annotated mitogenome sequences of Penthetria simplioipes and Plecia hardyi representing Bibionidae. This is the first report of the complete mitogenomes for the superfamily Bibionoidea. There are 37 genes in each of the complete mitogenomes of all 20 studied species from eight families of four superfamilies within infraorder Bibionomorpha. The Ka/Ks analysis suggests that all 13 PCGs have undergone purifying selection. The gene rearrangement events exist in some families (Keroplatidae, Sciaridae, and Cecidomyiidae) but not in Mycetophilidae in Sciaroidea and also in Scatopsoidea, Anisopodoidea, and Bibionoidea, which suggests that these rearrangement events are derived in the late period in the evolution of the Bibionomorpha. The phylogenetic analysis suggests the phylogenetic relationships of Scatopsoidea + (Anisopodoidea + (Bibionoidea + Sciaroidea)) in Bibionomorpha. The divergence time analysis suggests that Bibionomorpha originated in the Triassic, Scatopsoidea and Anisopodoidea in the late Triassic, Bibionoidea in the Jurassic, and Sciaroidea in the Jurassic to the Cretaceous. The work lays a base for the study of mitogenomes in Bibionomorpha but further work and broader taxon sampling are necessary for a better understanding of the phylogenetics and evolution of the infraorder.

Comparative Mitogenomics of Jumping Spiders with First Complete Mitochondrial Genomes of Euophryini (Araneae: Salticidae)

Salticidae is the most species-rich family of spiders with diverse morphology, ecology and behavior. However, the characteristics of the mitogenomes within this group are poorly understood with relatively few well-characterized complete mitochondrial genomes. In this study, we provide completely annotated mitogenomes for Corythalia opima and Parabathippus shelfordi, which represent the first complete mitogenomes of the tribe Euophryini of Salticidae. The features and characteristics of the mitochondrial genomes are elucidated for Salticidae by thoroughly comparing the known well-characterized mitogenomes. The gene rearrangement between trnL2 and trnN was found in two jumping spider species, Corythalia opima and Heliophanus lineiventris Simon, 1868. Additionally, the rearrangement of nad1 to between trnE and trnF found in Asemonea sichuanensis Song & Chai, 1992 is the first protein-coding gene rearrangement in Salticidae, which may have an important phylogenetic implication for the family. Tandem repeats of various copy numbers and lengths were discovered in three jumping spider species. The codon usage analyses showed that the evolution of codon usage bias in salticid mitogenomes was affected by both selection and mutational pressure, but selection may have played a more important role. The phylogenetic analyses provided insight into the taxonomy of Colopsus longipalpis (Żabka, 1985). The data presented in this study will improve our understanding of the evolution of mitochondrial genomes within Salticidae.

Diversity, Form, and Postembryonic Development of Paleozoic Insects

While Mesozoic, Paleogene, and Neogene insect faunas greatly resemble the modern one, the Paleozoic fauna provides unique insights into key innovations in insect evolution, such as the origin of wings and modifications of postembryonic development including holometaboly. Deep-divergence estimates suggest that the majority of contemporary insect orders originated in the Late Paleozoic, but these estimates reflect divergences between stem groups of each lineage rather than the later appearance of the crown groups. The fossil record shows the initial radiations of the extant hyperdiverse clades during the Early Permian, as well as the specialized fauna present before the End Permian mass extinction. This review summarizes the recent discoveries related to the documented diversity of Paleozoic hexapods, as well as current knowledge about what has actually been verified from fossil evidence as it relates to postembryonic development and the morphology of different body parts.

Mitochondrial genomes of soft scales (Hemiptera: Coccidae): features, structures and significance

BACKGROUND

Soft scales (Hemiptera: Coccidae), including important agricultural and forestry pests, are difficult to identify directly by morphological characters. Mitochondrial genomes (mitogenomes) have been widely used in species identification and phylogenetic research. However, only three complete mitogenomes, and very few mitochondrial genes of scale insects (Hemiptera: Coccoidea) can be searched in GenBank. Mitogenome comparisons between scale insects or between scale insects and other hemipteran species have not yet been reported.

RESULTS

In this study, detailed annotation of three new mitogenomes and comparative analysis of scale insects were completed, as well as comparative analysis of the gene composition, gene arrangement, codon usage and evolutionary forces between scale insects and 488 other hemipteran species for the first time. We found that high A + T content, gene rearrangement and truncated tRNAs are common phenomena in soft scales. The average A + T content and codon usage bias of scale insects are higher and stronger than those of other hemipteran insects, respectively. The atp8 gene of Hemiptera and nine other protein-coding genes of scale insects are under positive selection with higher evolutionary rates.

CONCLUSIONS

The study revealed the particularity of the scale insect mitogenomes, which will provide a good reference for future research on insect phylogenetic relationships, insect pest control, biogeography and identification.

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.

Complete mitochondrial genome of Rectorisluxiensis (Teleostei, Cyprinidae): characterisation and phylogenetic implications

Mitochondrial genomes (mitogenomes) are widely used in scientific studies on phylogenetic relationships, molecular evolution and population genetics. Here, we sequenced and analysed the mitogenome of Rectorisluxiensis, a Yangtze River drainage endemic, but threatened cyprinid fish of Labeoninae. The complete mitogenome of R.luxiensis was 16,592 bp in length, encoding 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region. The mitogenome showed a high A+T content (58.2%) and a positive AT-skew (0.10) and negative GC-skew (-0.25) base composition pattern. All the 13 PCGs were found to start with ATG codons, except for the COXI, in which GTG was the start codon. The ratio of non-synonymous and synonymous substitutions (Ka/Ks) of all the 13 PCGs were less than 1, indicating negative or purifying selection evolved in these genes. Comparatively speaking, the evolutionary rate of ATP8 was the fastest and ND4L was the slowest. All tRNAs could fold into a typical cloverleaf secondary structure, except tRNA^Ser1 that lacked a dihydrouridine arm. Phylogenetic relationships, based on the PCGs dataset of 91 mitogenomes of Labeoninae, showed that R.luxiensis grouped with Rectorisposehensis and they formed a monophyletic Rectoris. However, many non-monophyletic genera were revealed in labeoninae fishes, such as Cirrhinus, Decorus, Garra, Labeo and Pseudocrossocheilus, which indicated that the validities of some traditional genera required a further check. This study reported the complete mitogenome of R.luxiensis for the first time, which provided valuable data for future molecular evolution and conservation related studies of Rectoris and other species in Labeoninae.

Evolution of Gene Arrangements in the Mitogenomes of Ensifera and Characterization of the Complete Mitogenome of Schizodactylus jimo

Gene arrangement (relative location of genes) is another evolutionary marker of the mitogenome that can provide extensive information on the evolutionary mechanism. To explore the evolution of gene arrangements in the mitogenome of diversified Ensifera, we sequenced the mitogenome of the unique dune cricket species found in China and used it for phylogenetic analysis, in combination with 84 known Ensiferan mitogenomes. The mitogenome of Schizodactylus jimo is a 16,428-bp circular molecule that contains 37 genes. We identified eight types of gene arrangement in the 85 ensiferan mitogenomes. The gene location changes (i.e., gene translocation and duplication) were in three gene blocks: I-Q-M-ND2, rrnl-rns-V, and ND3-A-R-N-S-E-F. From the phylogenetic tree, we found that Schizodactylus jimo and most other species share a typical and ancient gene arrangement type (Type I), while Grylloidea has two types (Types II and III), and the other five types are rare and scattered in the phylogenetic tree. We deduced that the tandem replication–random loss model is the evolutionary mechanism of gene arrangements in Ensifera. Selection pressure analysis revealed that purifying selection dominated the evolution of the ensiferan mitochondrial genome. This study suggests that most gene rearrangements in the ensiferan mitogenome are rare accidental events.

Mitogenomic Comparison of the Mole Crickets Gryllotalpidae with the Phylogenetic Implications (Orthoptera: Ensifera)

Owing to limited molecular data, the phylogenetic position of the family Gryllotalpidae is still controversial in the infraorder Gryllidea. Mitochondrial genome (mitogenome) plays a crucial role in reconstructing phylogenetic relationships and revealing the molecular evolution of insects. However, only four mitogenomes have been reported in Gryllotalpidae to date. Herein, we obtained the first mitogenomes of Gryllotalpa henana Cai & Niu, 1998 and the Chinese G. orientalis Burmeister, 1838, made a detailed comparison of all mitogenomes available in Gryllotalpidae and reconstructed the phylogeny of Gryllidea based on mitogenomes using Bayesian inference (BI) and maximum likelihood (ML) methods. The results show that the complete mitogenome sequences of G. henana (15,504 bp) and G. orientalis (15,497 bp) are conserved, both exhibiting the double-stranded circular structure, typical gene content and the ancestral insect gene arrangement. The complete mitogenome of G.henana exhibits the lowest average AT content ever detected in Gryllotalpidae, and even Gryllidea. The gene nad2 of both species has atypical initiation codon GTG. All tRNAs exhibit typical clover-leaf structure, except for trnS1 lacking the dihydrouridine (DHU) arm. A potential stem-loop structure, containing a (T)_n(TC)₂(T)_n sequence, is detected in the control region of all gryllotalpids investigated and is likely related to the replication initiation of the minority strand. The phylogenetic analyses recover the six families of Gryllidea as Gryllotalpidae + (Myrmecophilidae + (Mogoplistidae + (Trigonidiidae + (Phalangopsidae + Gryllidae)))), similar to the trees based on transcriptomic and mitogenomic data. However, the trees are slightly different from the multilocus phylogenies, which show the sister-group relationship of Gryllotalpidae and Myrmecophilidae. The contradictions between mitogenomic and multilocus trees are briefly discussed.

Comparative mitogenomic analyses provide evolutionary insights into the retrolateral tibial apophysis clade (Araneae: Entelegynae)

The retrolateral tibial apophysis (RTA) clade is the largest spider lineage within Araneae. To better understand the diversity and evolution, we newly determined mitogenomes of ten RTA species from six families and performed a comparative mitogenomics analysis by combining them with 40 sequenced RTA mitogenomes available on GenBank. The ten mitogenomes encoded 37 typical mitochondrial genes and included a large non-coding region (putative control region). Nucleotide composition and codon usage were well conserved within the RTA clade, whereas diversity in sequence length and structural features was observed in control region. A reversal of strand asymmetry in nucleotide composition, i.e., negative AT-skews and positive GC-skews, was observed in each RTA species, likely resulting from mitochondrial gene rearrangements. All protein-coding genes were evolving under purifying selection, except for atp8 whose Ka/Ks was larger than 1, possibly due to positive selection or selection relaxation. Both mutation pressure and natural selection might contribute to codon usage bias of 13 protein-coding genes in the RTA lineage. Phylogenetic analyses based on mitogenomic data recovered a family-level phylogeny within the RTA; {[(Oval calamistrum clade, Dionycha), Marronoid clade], Sparassidae}. This study characterized RTA mitogenomes and provided some new insights into the phylogeny and evolution of the RTA clade.

Complete mitochondrial genomes of Thyreophagus entomophagus and Acarus siro (Sarcoptiformes: Astigmatina) provide insight into mitogenome features, evolution, and phylogeny among Acaroidea mites

Mites from the Acaroidea (Sarcoptiformes: Astigmatina) are important pests of various stored products, posing potential threats to preserved foods. In addition, mites can cause allergic diseases. Complete mitochondrial genomes (mitogenomes) are valuable resources for different research fields, including comparative genomics, molecular evolutionary analysis, and phylogenetic inference. We sequenced and annotated the complete mitogenomes of Thyreophagus entomophagus and Acarus siro. A comparative analysis was made between mitogenomic sequences from 10 species representing nine genera within Acaroidea. The mitogenomes of T. entomophagus and A. siro contained 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region. In Acaroidea species, mitogenomes have highly conserved gene size and order, and codon usage. Among Acaroidea mites, most PCGs were found to be under purifying selection, implying that most PCGs might have evolved slowly. Our findings showed that nad4 evolved most rapidly, whereas cox1 and cox3 evolved most slowly. The evolutionary rates of Acaroidea vary considerably across families. In addition, selection analyses were also performed in 23 astigmatid mite species, and the evolutionary rate of the same genes in different superfamilies exhibited large differences. Phylogenetic results are mostly consistent with those identified by previous phylogenetic studies on astigmatid mites. The monophyly of Acaroidea was rejected, and the Suidasiidae and Lardoglyphidae appeared to deviate from the Acaroidea branch. Our research proposed a review of the current Acaroidea classification system.

The Complete Mitochondrial Genomes of Four Species in the Subfamily Limenitidinae (Lepidoptera, Nymphalidae) and a Phylogenetic Analysis

Simple Summary

As it is currently delineated, the subfamily Limenitidinae (Lepidoptera: Nymphalidae) is comprised of 50 genera with approximately 1100 species. The classification of this subfamily has always been unstable. There are tribes and genera whose status is doubtful. Their phylogenetic relationships are far from being clarified and the monophyly of some of them is under question. To provide further insight into the relationships among included tribes, four newly-completed mitochondrial genomes of Limenitidinae (Neptis thisbe, Athyma zeroca, and Aldania raddei) have been sequenced and analyzed. Results indicate that the gene orientation and arrangement are similar to typical mitogenomes in Lepidoptera. The inferred phylogenetic analysis shows that tribe levels are well-supported monophyletic groups. Taken together, this work will provide a well-resolved framework for future study of this subfamily.

Abstract

The complete mitogenomes of four species, Neptis thisbe, Neptis obscurior, Athyma zeroca, and Aldania raddei, were sequenced with sizes ranging from 15,172 bp (N. obscurior) to 16,348 bp (Al. raddei). All four mitogenomes display similar nucleotide content and codon usage of protein-coding genes (PCGs). Typical cloverleaf secondary structures are identified in 21 tRNA genes, while trnS1 (AGN) lacks the dihydrouridine (DHC) arm. The gene orientation and arrangement of the four mitogenomes are similar to that of other typical mitogenomes of Lepidoptera. The Ka/Ks ratio of 13 PCGs among 58 Limenitidinae species reveals that cox1 had the slowest evolutionary rate, while atp8 and nad6 exhibited a higher evolutionary rate. The phylogenetic analysis reveals that tribe-levels are well-supported monophyletic groups. Additionally, Maximum Likelihood analysis recovered the relationship (Parthenini + ((Chalingini + (Cymothoini + Neptini)) + (Adoliadini + Limenitidini))). However, a Bayesian analysis based on the same dataset recovered the relationship (Parthenini + (Adoliadini + ((Cymothoini + Neptini) + (Chalingini + Limenitidini)))). These results will offer valuable data for the future study of the phylogenetic relationships for Limenitidinae.

Trade-off between flight capability and reproduction in Acridoidea (Insecta: Orthoptera)

In many insect taxa, there is a well-established trade-off between flight capability and reproduction. The wing types of Acridoidea exhibit extremely variability from full length to complete loss in many groups, thus, provide a good model for studying the trade-off between flight and reproduction. In this study, we completed the sampling of 63 Acridoidea species, measured the body length, wing length, body weight, flight muscle weight, testis and ovary weight, and the relative wing length (RWL), relative flight muscle weight (RFW), and gonadosomatic index (GSI) of different species were statistically analyzed. The results showed that there were significant differences in RWL, RFW, and GSI among Acridoidea species with different wing types. RFW of long-winged species was significantly higher than that of short-winged and wingless species (p < .01), while GSI of wingless species was higher than that of long-winged and short-winged species. The RWL and RFW had a strong positive correlation in species with different wing types (correlation coefficient r = .8344 for male and .7269 for female, and p < .05), while RFW was strong negatively correlated with GSI (r = -.2649 for male and -.5024 for female, and p < .05). For Acridoidea species with wing dimorphism, males with relatively long wings had higher RFW than that of females with relatively short wings, while females had higher GSI. Phylogenetic comparative analysis showed that RWL, RFW, and GSI all had phylogenetic signals and phylogenetic dependence. These results revealed that long-winged individuals are flight capable at the expense of reproduction, while short-winged and wingless individuals cannot fly, but has greater reproductive output. The results support the trade-off between flight and reproduction in Acridoidea.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Phylogenetic Position of the Genus Alulacris (Orthoptera: Acrididae: Melanoplinae: Podismini) Revealed by Complete Mitogenome Evidence

Whole mitogenomes are a useful data source for a wide variety of research goals due to the vastly cheaper sequencing cost and the far less demanding high-quality templates. The mitogenome has demonstrated great potential in resolving phylogenetic questions in Orthoptera at different taxonomic scales as well as exploring patterns of molecular and morphological character evolutions. In this study, the complete mitogenome of Alulacrisshilinensis (Zheng, 1977) was sequenced using next-generation sequencing, the characteristics of the mitogenome are presented briefly, and the phylogeny of the Melanoplinae and Catantopinae was reconstructed using a selected dataset of mitogenome sequences under maximum likelihood and Bayesian inference frameworks. The results show that the genus was consistently assigned to the subfamily Melanoplinae rather than Catantopinae in all phylogenetic trees deduced from different datasets under different frameworks, and this finding is entirely consistent with its morphological characters. Therefore, it is more appropriate to place the genus Alulacris in Melanoplinae rather than in Catantopinae.

Comparative Analysis of Transposable Elements in Genus Calliptamus Grasshoppers Revealed That Satellite DNA Contributes to Genome Size Variation

Simple Summary

Calliptamus is a genus of grasshoppers belonging to the family Acrididae. The genus Calliptamus includes approximately 17 recognized species. Calliptamus abbreviatus, Calliptamus italicus, and Calliptamus barbarus are three species that are widely found in northern China. These species are polyphagous, feeding on a variety of wild plants as well as crops, particularly legumes. The genome sizes, phylogenetic position, and transcriptome analysis of the genus Calliptamus were already known previous to this research. The repeatome analysis of these species was missing, which is directly linked to the larger genome sizes of the grasshoppers. Here, we classified repetitive DNA sequences at the level of superfamilies and sub-families, and found that LINE, TcMar-Tc1 and Ty3-gypsy LTR retrotransposons dominated the repeatomes of all genomes, accounting for 16–34% of the total genomes of these species. Satellite DNA dynamic evolutionary changes in all three genomes played a role in genome size evolution. This study would be a valuable source for future genome assemblies.

Abstract

Transposable elements (TEs) play a significant role in both eukaryotes and prokaryotes genome size evolution, structural changes, duplication, and functional variabilities. However, the large number of different repetitive DNA has hindered the process of assembling reference genomes, and the genus level TEs diversification of the grasshopper massive genomes is still under investigation. The genus Calliptamus diverged from Peripolus around 17 mya and its species divergence dated back about 8.5 mya, but their genome size shows rather large differences. Here, we used low-coverage Illumina unassembled short reads to investigate the effects of evolutionary dynamics of satDNAs and TEs on genome size variations. The Repeatexplorer2 analysis with 0.5X data resulted in 52%, 56%, and 55% as repetitive elements in the genomes of Calliptamus barbarus, Calliptamus italicus, and Calliptamus abbreviatus, respectively. The LINE and Ty3-gypsy LTR retrotransposons and TcMar-Tc1 dominated the repeatomes of all genomes, accounting for 16–35% of the total genomes of these species. Comparative analysis unveiled that most of the transposable elements (TEs) except satDNAs were highly conserved across three genomes in the genus Calliptamus grasshoppers. Out of a total of 20 satDNA families, 17 satDNA families were commonly shared with minor variations in abundance and divergence between three genomes, and 3 were Calliptamus barbarus specific. Our findings suggest that there is a significant amplification or contraction of satDNAs at genus phylogeny which is the main cause that made genome size different.

Complete mitogenomes of Anopheles peditaeniatus and Anopheles nitidus and phylogenetic relationships within the genus Anopheles inferred from mitogenomes

BACKGROUND

Despite the medical importance of mosquitoes of the genus Anopheles in the transmission of malaria and other human diseases, its phylogenetic relationships are not settled, and the characteristics of mitochondrial genome (mitogenome) are not thoroughly understood.

METHODS

The present study sequenced and analyzed the complete mitogenomes of An. peditaeniatus and An. nitidus, investigated genome characteristics, and inferred the phylogenetic relationships of 76 Anopheles spp.

RESULTS

The complete mitogenomes of An. peditaeniatus and An. nitidus are 15,416 and 15,418 bp long, respectively, and both include 13 PCGs, 22 tRNAs, two tRNAs and one control region (CR). Mitogenomes of Anopheles spp. are similar to those of other insects in general characteristics; however, the trnR and trnA have been reversed to "trnR-trnA," as has been reported in other mosquito genera. Genome variations mainly occur in CR length (493-886 bp) with six repeat unit types identified for the first time that demonstrate an evolutionary signal. The subgenera Lophopodomyia, Stethomyia, Kerteszia, Nyssorhynchus, Anopheles and Cellia are inferred to be monophyletic, and the phylogenetic analyses support a new phylogenetic relationship among the six subgenera investigated, in that subgenus Lophopodomyia is the sister to all other five subgenera, and the remaining five subgenera are divided into two clades, one of which is a sister-taxon subgenera Stethomyia + Kerteszia, and the other consists of subgenus Nyssorhynchus as the sister to a sister-group subgenera Anopheles + Cellia. Four series (Neomyzomyia, Pyretophorus, Neocellia and Myzomyia) of the subgenus Cellia, and two series (Arribalzagia and Myzorhynchus) of the subgenus Anopheles were found to be monophyletic, whereas three sections (Myzorhynchella, Argyritarsis and Albimanus) and their subdivisions of the subgenus Nyssorhynchus were polyphyletic or paraphyletic.

CONCLUSIONS

The study comprehensively uncovered the characteristics of mitogenome and the phylogenetics based on mitogenomes in the genus Anopheles, and provided information for further study on the mitogenomes, phylogenetics and taxonomic revision of the genus.

Cuticular Structures in Micropterous Crickets (Orthoptera, Gryllidae, Petaloptilini, Gryllomorphini)

Orthoptera is a very diverse group that has colonized practically all terrestrial ecosystems on the planet. They have adapted to live in the endogenous environment as well as in caves so that some species exhibit troglomorphic characteristics. This group has been extensively studied due to its economic and social importance; however, many basic morphological and biological questions remain to be solved. In this study, a comparative morphological study by scanning electron microscopy of different structures of eight species of micropteran crickets of the tribes Gryllomorphini and Petaloptilini, whose tegmina had lost their flight and song functionality was carried out. Special emphasis was placed on the tegmina and their possible relationship to reproductive functions. In addition, to assess troglomorphism in the genus Petaloptila, the biometric parameters of six other species have been considered. Actualization of the lifestyle of the studied species has also been carried out. The results show structures not previously described in this group (gland openings, setae, pores, or group of campaniform sensilla). Structures not previously described in this group have been detected, and tegmina (glandular openings and devoured tegmina) seem to confer a role in reproduction. Troglomorphisms are only observed in species of the subgenus Zapetaloptila. Statistically, significant differences have been found in characters such as cephalic elongation, ocular reduction, greater length of appendages, and depigmentation.

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.

The complete mitochondrial genome of the intertidal spider (Desis jiaxiangi) provides novel insights into the adaptive evolution of the mitogenome and the evolution of spiders

BACKGROUND

Although almost all extant spider species live in terrestrial environments, a few species live fully submerged in freshwater or seawater. The intertidal spiders (genus Desis) built silk nests within coral crevices can survive submerged in high tides. The diving bell spider, Argyroneta aquatica, resides in a similar dynamic environment but exclusively in freshwater. Given the pivotal role played by mitochondria in supplying most energy for physiological activity via oxidative phosphorylation and the environment, herein we sequenced the complete mitogenome of Desis jiaxiangi to investigate the adaptive evolution of the aquatic spider mitogenomes and the evolution of spiders.

RESULTS

We assembled a complete mitogenome of the intertidal spider Desis jiaxiangi and performed comparative mitochondrial analyses of data set comprising of Desis jiaxiangi and other 45 previously published spider mitogenome sequences, including that of Argyroneta aquatica. We found a unique transposition of trnL2 and trnN genes in Desis jiaxiangi. Our robust phylogenetic topology clearly deciphered the evolutionary relationships between Desis jiaxiangi and Argyroneta aquatica as well as other spiders. We dated the divergence of Desis jiaxiangi and Argyroneta aquatica to the late Cretaceous at ~ 98 Ma. Our selection analyses detected a positive selection signal in the nd4 gene of the aquatic branch comprising both Desis jiaxiangi and Argyroneta aquatica. Surprisingly, Pirata subpiraticus, Hypochilus thorelli, and Argyroneta aquatica each had a higher Ka/Ks value in the 13 PCGs dataset among 46 taxa with complete mitogenomes, and these three species also showed positive selection signal in the nd6 gene.

CONCLUSIONS

Our finding of the unique transposition of trnL2 and trnN genes indicates that these genes may have experienced rearrangements in the history of intertidal spider evolution. The positive selection signals in the nd4 and nd6 genes might enable a better understanding of the spider metabolic adaptations in relation to different environments. Our construction of a novel mitogenome for the intertidal spider thus sheds light on the evolutionary history of spiders and their mitogenomes.

Molecular phylogeny and species delimitation of the genus Tonkinacris (Orthoptera, Acrididae, Melanoplinae) from China

Tonkinacris is a small group in Acrididae. While a few species were occasionally sampled in some previous molecular studies, there is no revisionary research devoted to the genus. In this study, we explored the phylogeny of and the relationships among Chinese species of the genus Tonkinacris using the mitochondrial COI barcode and the complete sequences of ITS1 and ITS2 of the nuclear ribosomal DNA. The phylogeny was reconstructed in maximum likelihood and Bayesian inference frameworks, respectively. The overlap range between intraspecific variation and interspecific divergence was assessed via K2P distances. Species boundaries were delimitated using phylogenetic species concept, NJ tree, K2P distance, the statistical parsimony network as well as the GMYC model. The results demonstrate that the Chinese Tonkinacris species is a monophyletic group and the phylogenetic relationship among them is (T. sinensis, (T. meridionalis, (T. decoratus, T. damingshanus))). While T. sinensis, T. meridionalis and T. decoratus were confirmed being good independent species strongly supported by both morphological and molecular evidences, the validity of T. damingshanus was not perfectly supported by molecular evidence in this study.

Phylogenetic position of the genera Caryandoides, Paratoacris, Fer and Longchuanacris (Orthoptera:Acrididae) revealed by complete mitogenome sequences

Many taxa in the Acrididae have controversial phylogenetic positions. A typical example of such controversies is the phylogenetic positions of the genera Caryandoides, Paratoacris, Fer and Longchuanacris, as well as some other related taxa, which were placed in Oxyinae by some authors, but were considered members of the subfamily Catantopinae by others. In this study, the complete mitogenomes of nine species were sequenced using next-generation sequencing, the characteristics of the newly sequenced mitogenomes are presented briefly, and the phylogeny of the Oxyinae and Catantopinae are reconstructed using a selected dataset of mitogenome sequences under maximum likelihood and Bayesian inference frameworks. The results show that the four controversial genera were consistently assigned to the subfamily Oxyinae rather than Catantopinae in all phylogenetic trees deduced from different datasets under different frameworks, and this finding is entirely consistent with their morphological characters. Therefore, it is more appropriate to place them in Oxyinae rather than Catantopinae. In addition, the results from our analysis also confirm the membership of the genus Apalacris in Coptacrinae rather than Catantopinae, and indicate the uncertainty in the phylogenetic position of the genus Traulia, and a more in-depth study is necessary to resolve the relationship of Traulia with other catantopine groups or Coptacrinae.

Genomics Reveals Widespread Ecological Speciation in Flightless Insects

Recent genomic analyses have highlighted parallel divergence in response to ecological gradients, but the extent to which altitude can underpin such repeated speciation remains unclear. Wing reduction and flight loss have apparently evolved repeatedly in montane insect assemblages and have been suggested as important drivers of hexapod diversification. We test this hypothesis using genomic analyses of a widespread wing-polymorphic stonefly species complex in New Zealand. We identified over 50,000 polymorphic genetic markers generated across almost 200 Zelandoperla fenestrata stonefly specimens using a newly generated plecopteran reference genome, to reveal widespread parallel speciation between sympatric full-winged and wing-reduced ecotypes. Rather than the existence of a single, widespread, flightless taxon (Zelandoperla pennulata), evolutionary genomic data reveal that wing-reduced upland lineages have speciated repeatedly and independently from full-winged Z. fenestrata. This repeated evolution of reproductive isolation between local ecotype pairs that lack mitochondrial DNA differentiation suggests that ecological speciation has evolved recently. A cluster of outlier single-nucleotide polymorphisms detected in independently wing-reduced lineages, tightly linked in an approximately 85 kb genomic region that includes the developmental “supergene” doublesex, suggests that this “island of divergence” may play a key role in rapid ecological speciation. [Ecological speciation; genome assembly; genomic island of differentiation; genotyping-by-sequencing; incipient species; plecoptera; wing reduction.]

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.

MtOrt: an empirical mitochondrial amino acid substitution model for evolutionary studies of Orthoptera insects

BACKGROUND

Amino acid substitution models play an important role in inferring phylogenies from proteins. Although different amino acid substitution models have been proposed, only a few were estimated from mitochondrial protein sequences for specific taxa such as the mtArt model for Arthropoda. The increasing of mitochondrial genome data from broad Orthoptera taxa provides an opportunity to estimate the Orthoptera-specific mitochondrial amino acid empirical model.

RESULTS

We sequenced complete mitochondrial genomes of 54 Orthoptera species, and estimated an amino acid substitution model (named mtOrt) by maximum likelihood method based on the 283 complete mitochondrial genomes available currently. The results indicated that there are obvious differences between mtOrt and the existing models, and the new model can better fit the Orthoptera mitochondrial protein datasets. Moreover, topologies of trees constructed using mtOrt and existing models are frequently different. MtOrt does indeed have an impact on likelihood improvement as well as tree topologies. The comparisons between the topologies of trees constructed using mtOrt and existing models show that the new model outperforms the existing models in inferring phylogenies from Orthoptera mitochondrial protein data.

CONCLUSIONS

The new mitochondrial amino acid substitution model of Orthoptera shows obvious differences from the existing models, and outperforms the existing models in inferring phylogenies from Orthoptera mitochondrial protein sequences.