Mitochondrial genomes of four katydids (Orthoptera: Phaneropteridae): New gene rearrangements and their phylogenetic implications

Mitogenome architecture supports the non-monophyly of the cosmopolitan parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) recovered by nuclear and mitochondrial phylogenomics

Mitochondrial DNA gene organisation is an important source of phylogenetic information for various metazoan taxa at different evolutionary timescales, though this has not been broadly tested for all insect groups nor within a phylogenetic context. The cosmopolitan subfamily Doryctinae is a highly diverse group of braconid wasps mainly represented by ectoparasitoids of xylophagous beetle larvae. Previous molecular studies based on Sanger and genome-wide (ultraconserved elements, UCE; and mitochondrial genomes) sequence data have recovered a non-monophyletic Doryctinae, though the relationships involved have always been weakly supported. We characterised doryctine mitogenomes and conducted separate phylogenetic analyses based on mitogenome and UCE sequence data of ~100 representative doryctine genera to assess the monophyly and higher-level classification of the subfamily. We identified rearrangements of mitochondrial transfer RNAs (tRNAs) that support a non-monophyletic Doryctinae consisting of two separate non-related clades with strong geographic structure ('New World' and 'Old World' clades). This geographic structure was also consistently supported by the phylogenetic analyses preformed with mitogenome and UCE sequence data. These results highlight the utility of the mitogenome gene rearrangements as a potential source of phylogenetic information at different evolutionary timescales.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

Mitochondrial genome of Poecilimon cretensis (Orthoptera: Tettigoniidae: Phaneropterinae): Strong phylogenetic signals in gene overlapping regions

We report the complete mitochondrial genome of the Cretan bush cricket Poecilimon cretensis. The mitogenome consists of 13 protein-coding regions, 22 tRNAs, two rRNAs, and one control region. The length of mitogenome in P. cretensis varies between15477 and 15631 bp, mainly due to variability in control region. The start and stop codons of protein coding genes exhibit the general pattern in Phaneropterinae. Phylogenetic tree constructed with the mitogenome obtained during this study and 12 mitogenomes of Phaneropterinae downloaded from GenBank, placed P. cretensis in Barbitistini as sister group to Poecilimon luschani. Data indicate that the gene overlapping pattern exhibit strong phylogenetic signals.  

Evolution of mitogenomic gene order in Orthoptera

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

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.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea)

Simple Summary

Papilionoidea is a superfamily of Lepidoptera encompassing about 19,000 species. In the present work, we study the evolution of the structure of the mitogenome of these lepidopterans. The mechanisms generating the eight arrangements known for Papilionoidea were investigated analysing the movements of different mitochondrial genes. Five newly sequenced/assembled mitogenomes were included in our analysis involving more than 600 genomes. We provide new findings that help to understand the evolution of the gene orders MIQGO, IMQGO, 2S1GO, ES1GO and S1NGO in different butterflies. We demonstrate that the evolution of the 2S1GO in Lycaenidae followed a complicated pathway with multiple events of duplication and loss of trnS1 and changes in anticodon. We describe two new gene orders 2FFGO and 4QGO for Ampittia subvittatus (Hesperiidae) and Bhutanitis thaidina (Papilionidae).

Abstract

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

The complete mitochondrial genome of Gyps coprotheres (Aves, Accipitridae, Accipitriformes): phylogenetic analysis of mitogenome among raptors

Three species of Old World vultures on the Asian peninsula are slowly recovering from the lethal consequences of diclofenac. At present the reason for species sensitivity to diclofenac is unknown. Furthermore, it has since been demonstrated that other Old World vultures like the Cape (Gyps coprotheres; CGV) and griffon (G. fulvus) vultures are also susceptible to diclofenac toxicity. Oddly, the New World Turkey vulture (Cathartes aura) and pied crow (Corvus albus) are not susceptible to diclofenac toxicity. As a result of the latter, we postulate an evolutionary link to toxicity. As a first step in understanding the susceptibility to diclofenac toxicity, we use the CGV as a model species for phylogenetic evaluations, by comparing the relatedness of various raptor species known to be susceptible, non-susceptible and suspected by their relationship to the Cape vulture mitogenome. This was achieved by next generation sequencing and assembly. The Cape vulture mitogenome had a genome size of 16,908 bp. The mitogenome phylogenetic analysis indicated a close evolutionary relationship between Old World vultures and other members of the Accipitridae as indicated by bootstrap value of 100% on the phylogenetic trees. Based on this, we postulate that the other species could also be sensitive to the toxic effects of diclofenac. This warrants further investigations.

Evolution of mitochondrial gene arrangements in Arcidae (Bivalvia: Arcida) and their phylogenetic implications

Arcidae is a diverse group of ark shells with over 260 described species. The phylogenetic relationships and the evolution of the mitochondrial genomes in this family were poorly understood. Comparisons of mitogenomes have been widely used to explore the phylogenetic relationship among animal taxa. We described the complete mitogenomes of Arca navicularis, Scapharca gubernaculum and one nearly complete mitogenome of Anadara consociata. The mitogenome of A. navicularis (18,103 bp) is currently the smallest known Arcidae mitogenome, while the mitogenomes of S. gubernaculum (45,697 bp) and A. consociata (44,034 bp) are relatively large. The mitochondrial gene orders of the three taxa were substantially different from each other, as well as the patterns found in other ark shells. The relationships among Arcidae species recovered from different mitochondrial characters (nucleotide sequence versus gene order) were in disagreement. The phylogeny based on nucleotide sequences did not support the monophyly of Arcidae, as Cucullaea labiata (Cucullaeidae) appeared as a subgroup within Arcinae, rather than sister group to the family Arcidae. In addition, we presented the first time-calibrated evolutionary tree of Arcidae based on mitochondrial DNA (mtDNA) sequences, which placed the deepest divergence within Arcidae at 342.36 million years ago (Mya), around the Carboniferous (360-300 Mya).

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.

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

Orthoptera is the most diverse order of polyneopterans, and the forewing and hindwing of its members exhibit extremely variability from full length to complete loss in many groups; thus, this order provides a good model for studying the effects of insect flight ability on the evolutionary constraints on and evolutionary rate of the mitochondrial genome. Based on a data set of mitochondrial genomes from 171 species, including 43 newly determined, we reconstructed Orthoptera phylogenetic relationships and estimated the divergence times of this group. The results supported Caelifera and Ensifera as two monophyletic groups, and revealed that Orthoptera originated in the Carboniferous (298.997 Mya). The date of divergence between the suborders Caelifera and Ensifera was 255.705 Mya, in the late Permian. The major lineages of Acrididae seemed to have radiated in the Cenozoic, and the six patterns of rearrangement of 171 Orthoptera mitogenomes mostly occurred in the Cretaceous and Cenozoic. Based on phylogenetic relationships and ancestral state reconstruction, we analysed the evolutionary selection pressure on and evolutionary rate of mitochondrial protein-coding genes (mPCGs). The results indicated that during approximately 300 Mya of evolution, these genes experienced purifying selection to maintain their function. Flightless orthopteran insects accumulated more non-synonymous mutations than flying species and experienced more relaxed evolutionary constraints. The different wing types had different evolutionary rates, and the mean evolutionary rate of Orthoptera mitochondrial mPCGs was 13.554 × 10^-9 subs/s/y. The differences in selection pressures and evolutionary rates observed between the mitochondrial genomes suggested that functional constraints due to locomotion play an important role in the evolution of mitochondrial DNA in orthopteran insects with different wing types.

The complete mitochondrial genomes of two vent squat lobsters, Munidopsis lauensis and M. verrilli: Novel gene arrangements and phylogenetic implications

Hydrothermal vents are considered as one of the most extremely harsh environments on the Earth. In this study, the complete mitogenomes of hydrothermal vent squat lobsters, Munidopsis lauensis and M. verrilli, were determined through Illumina sequencing and compared with other available mitogenomes of anomurans. The mitogenomes of M. lauensis (17,483 bp) and M. verrilli (17,636 bp) are the largest among all Anomura mitogenomes, while the A+T contents of M. lauensis (62.40%) and M. verrilli (63.99%) are the lowest. The mitogenomes of M. lauensis and M. verrilli display novel gene arrangements, which might be the result of three tandem duplication-random loss (tdrl) events from the ancestral pancrustacean pattern. The mitochondrial gene orders of M. lauensis and M. verrilli shared the most similarities with S. crosnieri. The phylogenetic analyses based on both gene order data and nucleotide sequences (PCGs and rRNAs) revealed that the two species were closely related to Shinkaia crosnieri. Positive selection analysis revealed that eighteen residues in seven genes (atp8, Cytb, nad3, nad4, nad4l, nad5, and nad6) of the hydrothermal vent anomurans were positively selected sites.

The complete mitochondrial genome of Xizicus (Haploxizicus) maculatus revealed by Next-Generation Sequencing and phylogenetic implication (Orthoptera, Meconematinae)

Xizicus Gorochov, 1993, the quiet-calling katydid, is a diverse genus with 68 species in world, which includes more than 45 species in China, has undergone numerous taxonomic revisions with contradicting conclusions. In this study the complete mitochondrial genome of Xizicus (Haploxizicus) maculatus collected from Hainan for the first time was sequenced using the Next-Generation Sequencing (NGS) technology. The length of whole mitogenome is 16,358 bp and contains the typical gene arrangement, base composition, and codon usage found in other related species. The overall base composition of the mitochondrial genome is 37.0 % A, 32.2 % T, 20.2 % C, and 10.6 % G. All 13 protein-coding genes (PCGs) began with typical ATN initiation codon. Nine of the 13 PCGs have a complete termination codon, but the remaining four genes (COI, COIII, ND5, and ND4) terminate with an incomplete T. Phylogenetic analyses are carried out based on the concatenated dataset of 13 PCGs and two rRNAs of Tettigoniidae species available in GenBank. Both Bayesian inference and Maximum Likelihood analyses recovered each subfamily as a monophyletic group. Regardless of the position of Lipotactinae, the relationships among the subfamilies of Tettigoniidae were as follows: ((((Tettigoniinae, Bradyporinae) Meconematinae) Conocephalinae) Hexacentrinae). The topological structure of the phylogeny trees showed that the Xizicus (Haploxizicus) maculatus is closer to Xizicus (Xizicus) fascipes than Xizicus (Eoxizicus) howardi.

Two new complete mitochondrial genomes of Dorcus stag beetles (Coleoptera, Lucanidae)

The systematics of Dorcus MacLeay has been a long-standing debate. Mitochondrial genomes were widely used to deeply understand the phylogeny of problematic taxa in virtue of their genetic importance and comprehensiveness. To provide more useful genetic data for resolving the systematic disputation of Dorcus stag beetles. The complete mitochondrial genomes of Dorcus hopei and Dorcus seguyi were obtained using the next generation sequencing. Characteristics of the two genomes are explicated through comparing their genome organization and base composition, protein-coding genes and codon usage, intergenic spacers and non-coding region, transfer and ribosomal RNA genes and control region. Phylogenetic relationships were reconstructed using Maximum likelihood and Bayesian inference analyses based on the concatenated nucleotide sequences of 13 PCGs from 9 stag beetles and 3 scarab beetles. The complete mitogenomes of D. hopei and D. seguyi was 16,026 bp/17,955 bp long, respectively. A tandem repeat with the length of 940 bp was presented in the A+T-rich region in D. hopei. An unexpected non-coding region of 332 bp was located between nad2 and trnW in D. seguyi. The phylogenetic analyses robustly supported that D. hopei formed a branch with the generic type of D. parallelipipedus. Whereas D. seguyi was not covered in the branch of (D. hopei + D. parallelipipedus), but was sister to them. The results indicated that D. hopei should be a good member of Dorcus MacLeay. The taxonomic status of D. seguyi remained to be studied furtherly.

The rearranged mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae), a parasitoid wasp of mantis

The complete sequence of the mitochondrial genome of Podagrion sp. (Hymenoptera: Torymidae) is described. The mitogenome was 15,845 bp in size, and contained typical sets of mitochondrial genes. The base composition of the Podagrion sp. mitogenome was also biased toward A + T bases (81.8%). The mitochondrial genome of Podagrion sp. has a weak AT skew (0.07) and a strong GC skew (-0.26). Podagrion sp. exhibits a novel rearrangement compared with the ancestral order, including six protein-coding genes (nad3, cox3, atp6, atp8, cox2 and cox1), which have inverted to the minor strand from the major strand. The A + T-rich region of Podagrion sp., which is located between trnN and trnI, have five tandem repeats. The apomorphic rearrangements, including the conserved block "cox3-atp6-atp8-cox2-cox1-nad5-nad4-nad4l-nad6-cob" and the special locations of trnV and trnA, were mapped onto the phylogeny of Proctotrupomorpha.

The highly rearranged mitochondrial genomes of the crabs Maja crispata and Maja squinado (Majidae) and gene order evolution in Brachyura

We sequenced the mitochondrial genomes of the spider crabs Maja crispata and Maja squinado (Majidae, Brachyura). Both genomes contain the whole set of 37 genes characteristic of Bilaterian genomes, encoded on both α- and β-strands. Both species exhibit the same gene order, which is unique among known animal genomes. In particular, all the genes located on the β-strand form a single block. This gene order was analysed together with the other nine gene orders known for the Brachyura. Our study confirms that the most widespread gene order (BraGO) represents the plesiomorphic condition for Brachyura and was established at the onset of this clade. All other gene orders are the result of transformational pathways originating from BraGO. The different gene orders exhibit variable levels of genes rearrangements, which involve only tRNAs or all types of genes. Local homoplastic arrangements were identified, while complete gene orders remain unique and represent signatures that can have a diagnostic value. Brachyura appear to be a hot-spot of gene order diversity within the phylum Arthropoda. Our analysis, allowed to track, for the first time, the fully evolutionary pathways producing the Brachyuran gene orders. This goal was achieved by coupling sophisticated bioinformatic tools with phylogenetic analysis.

Comparative Mitogenomic Analyses of Praying Mantises (Dictyoptera, Mantodea): Origin and Evolution of Unusual Intergenic Gaps

Praying mantises are a diverse group of predatory insects. Although some Mantodea mitogenomes have been reported, a comprehensive comparative and evolutionary genomic study is lacking for this group. In the present study, four new mitogenomes were sequenced, annotated, and compared to the previously published mitogenomes of other Mantodea species. Most Mantodea mitogenomes share a typical set of mitochondrial genes and a putative control region (CR). Additionally, and most intriguingly, another large non-coding region (LNC) was detected between trnM and ND2 in all six Paramantini mitogenomes examined. The main section in this common region of Paramantini may have initially originated from the corresponding control region for each species, whereas sequence differences between the LNCs and CRs and phylogenetic analyses indicate that LNC and CR are largely independently evolving. Namely, the LNC (the duplicated CR) may have subsequently degenerated during evolution. Furthermore, evidence suggests that special intergenic gaps have been introduced in some species through gene rearrangement and duplication. These gaps are actually the original abutting sequences of migrated or duplicated genes. Some gaps (G5 and G6) are homologous to the 5' and 3' surrounding regions of the duplicated gene in the original gene order, and another specific gap (G7) has tandem repeats. We analysed the phylogenetic relationships of fifteen Mantodea species using 37 concatenated mitochondrial genes and detected several synapomorphies unique to species in some clades.

Towards a higher-level Ensifera phylogeny inferred from mitogenome sequences

Although mitogenomes are useful tools for inferring evolutionary history, only a few representative ones can be used for most Ensifera lineages. Thirty-two ensiferan mitogenomes were determined using ABI Sanger sequencing and standard primer walking of 2-3 overlapping Long-PCR fragments, or Illumina® HiSeq2000 for "shotgun" sequenced long-PCR-amplified mitochondrial or total genomic DNA. Six patterns of gene arrangements, including the novel trnR-trnS^AGN-trnA-trnN-trnG-nad3 in Lipotactes tripyrga (Lipotactinae), were identified from 59 ensiferan mitogenomes. The results suggest that trnM-trnI-trnQ and trnA-trnR-trnE-trnS^AGN-trnN-trnF rearrangements might be a shared derived character in Pseudophyllinae and Gryllidae, respectively. We found base composition biases in our dataset, which potentially complicate the inference of higher-level ensiferan phylogeny. Site-heterogeneous Bayesian inference (BI) and site-homogeneous maximum likelihood (ML) analyses recovered all ensiferan superfamilies as monophyletic. The site-homogeneous BI analysis failed to recover the monophyly of Stenopelmatoidea. As Schizodactyloidea was only represented by Comicus campestris, its monophyly could not be tested. In the Triassic/Jurassic boundary, Ensifera diverged into grylloid and non-grylloid clades. All analyses confirmed Grylloidea and Gryllotalpoidea as sister groups. Site-heterogeneous BI analysis found Schizodactyloidea as the most basal lineage and sister to the clade formed by Grylloidea and Gryllotalpoidea, but the site-homogeneous analyses placed it basally to the non-grylloid clade and recovered a sister relationship between Tettigonioidea and (Hagloidea, Rhaphidophoroidea, Stenopelmatoidea), although this clade had a low support. The site-heterogeneous BI analysis found Tettigonioidea and Hagloidea were sister groups (posterior probability (PP)=0.99), Stenopelmatoidea was sister to (Tettigonioidea, Hagloidea) (PP>0.91), and Rhaphidophoroidea was basal to the non-grylloid clade. At a lower level, all analyses divided Tettigonioidea into Phaneropteridae and Tettigoniidae.

New Mitogenomes of Two Chinese Stag Beetles (Coleoptera, Lucanidae) and Their Implications for Systematics

Although conspicuous and well-studied, stag beetles have been slow to join the genomic era. In this study, mitochondrial genomes of two stag beetles, Sinodendron yunnanense and Prosopocoilus confucius, are sequenced for the first time. Both of their genomes consisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. The mitogenome of S. yunnanense was 16,921 bp in length, and P. confucius was 16,951 bp. The location of the gene trnL(UUR), between the A + T-rich and control region in S. yunnanense, is the first observed in Lucanidae. In P. confucius, an unexpected noncoding region of 580 bp was discovered. Maximum likelihood and Bayesian inference on the 13 mitochondrial PCGs were used to infer the phylogenetic relationships among 12 representative stag beetles and three scarab beetles. The topology of the two phylogenetic trees was almost identical: S. yunnanense was recovered as the most basal Lucanid, and the genus Prosopocoilus was polyphyletic due to P. gracilis being recovered sister to the genera Dorcus and Hemisodorcus. The phylogenetic results, genetic distances and mitogenomic characteristics call into question the cohesion of the genus Prosopocoilus. The genetic resources and findings herein attempts to redress understudied systematics and mitogenomics of the stag beetles.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.