The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

The complete mitochondrial genome of Eterusia aedea (Lepidoptera,Zygaenidae) and comparison with other zygaenid moths

Zygaenidae comprises >1036 species, including many folivorous pests in agriculture. In the present study, the complete mitochondrial genome (mitogenome) of a major pest of tea trees, Eterusia aedea was determined. The 15,196-bp circular genome contained the common set of 37 mitochondrial genes (including 13 protein-coding genes, two rRNA genes, and 22 tRNA genes) and exhibited the similar genomic features to reported Zygaenidae mitogenome. Comparative analyses of Zygaenidae mitogenomes showed a typical evolutionary trend of lepidopteran mitogenomes. In addition, we also investigated the gene order of lepidopteran mitogenomes and proposed that the novel gene order trnA-trnR-trnN-trnE-trnS-trnF from Zygaenidae and Gelechiidae and most other gene rearrangements of this tRNA cluster evolved independently. Finally, the mitogenomic phylogeny of Lepidoptera was reconstructed based on multiple mitochondrial datasets. And all the phylogenetic results revealed the sister relationships of Cossoidea and Zygaenoidea with both BI and ML methods, which is the first stable mitogenomic evidence for this clade.     

Towards a comprehensive picture of alloacceptor tRNA remolding in metazoan mitochondrial genomes

Remolding of tRNAs is a well-documented process in mitochondrial genomes that changes the identity of a tRNA. It involves a duplication of a tRNA gene, a mutation that changes the anticodon and the loss of the ancestral tRNA gene. The net effect is a functional tRNA that is more closely related to tRNAs of a different alloacceptor family than to tRNAs with the same anticodon in related species. Beyond being of interest for understanding mitochondrial tRNA function and evolution, tRNA remolding events can lead to artifacts in the annotation of mitogenomes and thus in studies of mitogenomic evolution. Therefore, it is important to identify and catalog these events. Here we describe novel methods to detect tRNA remolding in large-scale data sets and apply them to survey tRNA remolding throughout animal evolution. We identify several novel remolding events in addition to the ones previously mentioned in the literature. A detailed analysis of these remoldings showed that many of them are derived from ancestral events.     

Next-generation sequencing,phylogenetic signal and comparative mitogenomic analyses in Metacrangonyctidae (Amphipoda: Crustacea)

Background

Comparative mitochondrial genomic analyses are rare among crustaceans below the family or genus level. The obliged subterranean crustacean amphipods of the family Metacrangonyctidae, found from the Hispaniola (Antilles) to the Middle East, including the Canary Islands and the peri-Mediterranean region, have an evolutionary history and peculiar biogeography that can respond to Tethyan vicariance. Indeed, recent phylogenetic analysis using all protein-coding mitochondrial sequences and one nuclear ribosomal gene have lent support to this hypothesis (Bauzà-Ribot et al. 2012).

Results

We present the analyses of mitochondrial genome sequences of 21 metacrangonyctids in the genera Metacrangonyx and Longipodacrangonyx, covering the entire geographical range of the family. Most mitogenomes were attained by next-generation sequencing techniques using long-PCR fragments sequenced by Roche FLX/454 or GS Junior pyro-sequencing, obtaining a coverage depth per nucleotide of up to 281×. All mitogenomes were AT-rich and included the usual 37 genes of the metazoan mitochondrial genome, but showed a unique derived gene order not matched in any other amphipod mitogenome. We compare and discuss features such as strand bias, phylogenetic informativeness, non-synonymous/synonymous substitution rates and other mitogenomic characteristics, including ribosomal and transfer RNAs annotation and structure.

Conclusions

Next-generation sequencing of pooled long-PCR amplicons can help to rapidly generate mitogenomic information of a high number of related species to be used in phylogenetic and genomic evolutionary studies. The mitogenomes of the Metacrangonyctidae have the usual characteristics of the metazoan mitogenomes (circular molecules of 15,000-16,000 bp, coding for 13 protein genes, 22 tRNAs and two ribosomal genes) and show a conserved gene order with several rearrangements with respect to the presumed Pancrustacean ground pattern. Strand nucleotide bias appears to be reversed with respect to the condition displayed in the majority of crustacean mitogenomes since metacrangonyctids show a GC-skew at the (+) and (-) strands; this feature has been reported also in the few mitogenomes of Isopoda (Peracarida) known thus far. The features of the rRNAs, tRNAs and sequence motifs of the control region of the Metacrangonyctidae are similar to those of the few crustaceans studied at present.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-566) contains supplementary material, which is available to authorized users.     

Phylomitogenomics reveals mito-nuclear concordance in social wasps: The performance of mitochondrial markers and gene order for hymenopteran systematics

Mitochondrial (mtDNA) genes have served as widely utilised genetic loci for phylogenetic and phylogeographic studies of animals. However, the phylogenetic performance of many mtDNA genes has not been empirically evaluated across lineages within hymenopteran wasps. To address this question, we assembled and analysed mitogenomic data from social wasps, representing the four recognised tribes of Polistinae and all Epiponini genera. Additionally, we evaluated whether mtDNA gene order in Polistinae is congruent with its tribal classification. Using concatenation phylogenetic methods, we show phylogenetic congruence between mitogenomic and nuclear data. Statistically comparing the phylogenetic performance of individual mtDNA genes, we demonstrate that for social wasps the molecular markers COI, 16S, NAD5, and NAD2 perform best, while ATP6, COII, and 12S show the worst results. Finally, we verified that the tRNA cluster close to the noncoding region is a hotspot of genetic rearrangements in Vespidae and can be used as additional information for the systematics of this group. Together, these results indicate that mitogenomes contain robust phylogenetic signal to elucidate the evolutionary history of Vespidae. Moreover, our study identifies the best choice of mtDNA markers for systematic investigations of social wasps.     

Creating,curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information—especially sequence information that includes intraspecific variation—creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial ‘barcode’ regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.     

Mitogenomic exploration of higher teleostean phylogenies: a case study for moderate-scale evolutionary genomics with 38 newly determined complete mitochondrial DNA sequences

Although adequate resolution of higher-level relationships of organisms apparently requires longer DNA sequences than those currently being analyzed, limitations of time and resources present difficulties in obtaining such sequences from many taxa. For fishes, these difficulties have been overcome by the development of a PCR-based approach for sequencing the complete mitochondrial genome (mitogenome), which employs a long PCR technique and many fish-versatile PCR primers. In addition, recent studies have demonstrated that such mitogenomic data are useful and decisive in resolving persistent controversies over higher-level relationships of teleosts. As a first step toward resolution of higher teleostean relationships, which have been described as the "(unresolved) bush at the top of the tree," we investigated relationships using mitogenomic data from 48 purposefully chosen teleosts, of which those from 38 were newly determined during the present study (a total of 632,315 bp), using the above method. Maximum-parsimony and maximum-likelihood analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and third codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) from the 48 species. The resultant two trees from the two methods were well resolved and largely congruent, with many internal branches supported by high statistical values. The tree topologies themselves, however, exhibited considerable variation from the previous morphology-based cladistic hypotheses, with most of the latter being confidently rejected by the mitogenomic data. Such incongruence resulted largely from the phylogenetic positions or limits of long-standing problematic taxa, which were quite unexpected from previous morphological and molecular analyses. We concluded that the present study provided a basis of and guidelines for future investigations of teleostean evolutionary mitogenomics and that purposeful higher-density taxonomic sampling, subsequent sequencing efforts, and phylogenetic analyses of their mitogenomes may be decisive in resolving persistent controversies over higher-level relationships of teleosts, the most diversified group of all vertebrates, comprising over 23,500 extant species.     

Mitogenomes provide new insights of evolutionary history of Boreheptagyiini and Diamesini (Diptera: Chironomidae: Diamesinae)

Mitogenomes have been widely used for phylogenetic reconstruction of various Dipteran groups, but specifically for chironomid, they have not been carried out to resolve the relationships. Diamesinae (Diptera: Chironomidae) are important bioindicators for freshwater ecosystem monitoring, but its evolutionary history remains uncertain for lack of information. Here, coupled with one previously published and 30 new mitogenomes of Diamesinae, we carried out comparative mitogenomic analysis and phylogenetic analysis. Mitogenomes of Diamesinae were conserved in structure, and all genes arranged in the same order as the ancestral insect mitogenome. All protein‐coding genes in Diamesinae were under stronger purifying selection than those of other nonbiting midge species, which may exhibit signs of adaptation to life at cold living conditions. Phylogenetic analyses strongly supported the monophyly of Diamesinae, with Boreheptagyiini deeply nested within Diamesini. In addition, phylogenetic relationship of selected six genera was resolved, except Sympotthastia remained unstable. Our study revealed that the mitogenomes of Diamesinae are highly conserved, and they are practically useful for phylogenetic inference.     

Comparative Mitogenomic Analysis of Damsel Bugs Representing Three Tribes in the Family Nabidae (Insecta: Hemiptera)

Background

Nabidae, a family of predatory heteropterans, includes two subfamilies and five tribes. We previously reported the complete mitogenome of Alloeorhynchus bakeri, a representative of the tribe Prostemmatini in the subfamily Prostemmatinae. To gain a better understanding of architecture and evolution of mitogenome in Nabidae, mitogenomes of five species representing two tribes (Gorpini and Nabini) in the subfamily Nabinae were sequenced, and a comparative mitogenomic analysis of three nabid tribes in two subfamilies was carried out.

Methodology/Principal Findings

Nabid mitogenomes share a similar nucleotide composition and base bias, except for the control region, where differences are observed at the subfamily level. In addition, the pattern of codon usage is influenced by the GC content and consistent with the standard invertebrate mitochondrial genetic code and the preference for A+T-rich codons. The comparison among orthologous protein-coding genes shows that different genes have been subject to different rates of molecular evolution correlated with the GC content. The stems and anticodon loops of tRNAs are extremely conserved, and the nucleotide substitutions are largely restricted to TψC and DHU loops and extra arms, with insertion-deletion polymorphisms. Comparative analysis shows similar rates of substitution between the two rRNAs. Long non-coding regions are observed in most Gorpini and Nabini mtDNAs in-between trnI-trnQ and/or trnS2-nad1. The lone exception, Nabis apicalis, however, has lost three tRNAs. Overall, phylogenetic analysis using mitogenomic data is consistent with phylogenies constructed mainly form morphological traits.

Conclusions/Significance

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the family Nabidae. Nucleotide diversity and mitogenomic traits are phylogenetically informative at subfamily level. Furthermore, inclusion of a broader range of samples representing various taxonomic levels is critical for the understanding of mitogenomic evolution in damsel bugs.     

Mitogenomes provide new insights into the evolutionary history of Prodiamesinae (Diptera: Chironomidae)

Evolutionary analysis of Prodiamesinae has long been impeded by lack of information, and its phylogenetic relationship with Orthocladiinae remains questionable. Here, ten complete mitochondrial genomes (mitogenomes) of Orthocladiinae sensu lato were newly sequenced, including three Prodiamesinae species and seven Orthocladiinae species. Coupled with published mitogenomes, a total of 12 mitogenomes of Orthocladiinae sensu lato were selected for a comparative mitogenomic analysis and phylogenetic reconstruction. Mitogenomes of Orthocladiinae sensu lato are conserved in structure, and all genes arrange the same gene order as the ancestral insect mitogenome. Nucleotide composition is highly biased, and the control region displayed the highest A + T content. All protein-coding genes are under purifying selection, and the ATP8 evolves at the fastest rate. In addition, the mitogenomes of Orthocladiinae sensu lato are highly conserved, and they are practically useful for phylogenetic inference, suggesting a re-classification of Orthocladiinae by sinking Prodiamesinae as a subgroup of Orthocladiinae.     

Morphological Characters Are Compatible with Mitogenomic Data in Resolving the Phylogeny of Nymphalid Butterflies (Lepidoptera: Papilionoidea: Nymphalidae)

Nymphalidae is the largest family of butterflies with their phylogenetic relationships not adequately approached to date. The mitochondrial genomes (mitogenomes) of 11 new nymphalid species were reported and a comparative mitogenomic analysis was conducted together with other 22 available nymphalid mitogenomes. A phylogenetic analysis of the 33 species from all 13 currently recognized nymphalid subfamilies was done based on the mitogenomic data set with three Lycaenidae species as the outgroups. The mitogenome comparison showed that the eleven new mitogenomes were similar with those of other butterflies in gene content and order. The reconstructed phylogenetic trees reveal that the nymphalids are made up of five major clades (the nymphaline, heliconiine, satyrine, danaine and libytheine clades), with sister relationship between subfamilies Cyrestinae and Biblidinae, and most likely between subfamilies Morphinae and Satyrinae. This whole mitogenome-based phylogeny is generally congruent with those of former studies based on nuclear-gene and mitogenomic analyses, but differs considerably from the result of morphological cladistic analysis, such as the basal position of Libytheinae in morpho-phylogeny is not confirmed in molecular studies. However, we found that the mitogenomic phylogeny established herein is compatible with selected morphological characters (including developmental and adult morpho-characters).     

Evolutionary Analysis of Mitogenomes from Parasitic and Free-Living Flatworms

Mitochondrial genomes (mitogenomes) are useful and relatively accessible sources of molecular data to explore and understand the evolutionary history and relationships of eukaryotic organisms across diverse taxonomic levels. The availability of complete mitogenomes from Platyhelminthes is limited; of the 40 or so published most are from parasitic flatworms (Neodermata). Here, we present the mitogenomes of two free-living flatworms (Tricladida): the complete genome of the freshwater species Crenobia alpina (Planariidae) and a nearly complete genome of the land planarian Obama sp. (Geoplanidae). Moreover, we have reanotated the published mitogenome of the species Dugesia japonica (Dugesiidae). This contribution almost doubles the total number of mtDNAs published for Tricladida, a species-rich group including model organisms and economically important invasive species. We took the opportunity to conduct comparative mitogenomic analyses between available free-living and selected parasitic flatworms in order to gain insights into the putative effect of life cycle on nucleotide composition through mutation and natural selection. Unexpectedly, we did not find any molecular hallmark of a selective relaxation in mitogenomes of parasitic flatworms; on the contrary, three out of the four studied free-living triclad mitogenomes exhibit higher A+T content and selective relaxation levels. Additionally, we provide new and valuable molecular data to develop markers for future phylogenetic studies on planariids and geoplanids.     

Structural features of the mitogenome of the leafhopper genus Cladolidia (Hemiptera: Cicadellidae: Coelidiinae) and phylogenetic implications in Cicadellidae

The first two complete mitogenomes of the leafhopper genus Cladolidia (C. biungulata and C. robusta) were sequenced and annotated to further explore the phylogeny of Cladolidia. Both the newly sequenced mitogenomes have a typical circular structure, with lengths of 15,247 and 15,376 bp and A + T contents of 78.2% and 78%, respectively. We identified a highly conserved genome organization in the two Cladolidia spp. through comparative analysis that included the following assessments: genome content, gene order, nucleotide composition, codon usage, amino acid composition, and tRNA secondary structure. Moreover, we detected the base heterogeneity of Cicadellidae mitogenomic data and constructed phylogenetic trees using the nucleotide alignments of 12 subfamilies of 58 leafhopper species. We noted a weak heterogeneity in the base composition among the Cicadellidae mitogenomes. Phylogenetic analyses showed that the monophyly of each subfamily was generally well supported in the family Cicadellidae; the main topology was as follows: (Deltocephalinae + (Treehoppers + ((Megophthalminae + (Macropsinae + (Hylicinae + (Coelidiinae +Iassinae)) + (Idiocerinae + (Cicadellinae + (Typhlocybinae + (Mileewinae + (Evacanthinae +Ledrinae)))))))))). Within Coelidiinae, phylogenetic analyses revealed that C. biungulata and C. robusta belong to Coelidiinae and the monophyly of Cladolidia is well supported. In addition, on the basis of complete mitogenome phylogenetic analysis and the comparison of morphological characteristics, we further confirm the genus Olidiana as a paraphyletic group, suggesting that the genus may need taxonomic revisions.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-930) contains supplementary material, which is available to authorized users.     

Mitochondrial genomes of three Bostrichiformia species and phylogenetic analysis of Polyphaga (Insecta,Coleoptera)

Here we determined mitogenomes of three Bostrichiformia species. These data were combined with 51 previously sequenced Polyphaga mitogenomes to explore the higher-level relationships within Polyphaga by using four different mitogenomic datasets and three tree inference approaches. Among Polyphaga mitogenomes we observed heterogeneity in nucleotide composition and evolutionary rates, which may have affected phylogenetic inferences across the different mitogenomic datasets. Elateriformia, Cucujiformia, and Scarabaeiformia were each inferred to be monophyletic by all analyses, as was Bostrichiformia by most analyses based on two datasets with low heterogeneity. The large series Staphyliniformia was never recovered as monophyletic in our analyses. The Bayesian tree using a degenerated nucleotide dataset (P123_Degen) and a site-heterogeneous mixture model in PhyloBayes was supported as the best Polyphaga phylogeny: (Scirtiformia, (Elateriformia, ((Bostrichiformia, Cucujiformia), (Scarabaeiformia + Staphyliniformia)))). For Cucujiformia, the largest series, we inferred a superfamily-level phylogeny: ((Cleroidea, Coccinelloidea), (Tenebrionoidea, (Cucujoidea + Curculionoidea + Chrysomeloidea))).     

Ancient DNA Suggests Single Colonization and Within-Archipelago Diversification of Caribbean Caviomorph Rodents

Reconstructing the evolutionary history of island biotas is complicated by unusual morphological evolution in insular environments. However, past human-caused extinctions limit the use of molecular analyses to determine origins and affinities of enigmatic island taxa. The Caribbean formerly contained a morphologically diverse assemblage of caviomorph rodents (33 species in 19 genera), ranging from ∼0.1 to 200 kg and traditionally classified into three higher-order taxa (Capromyidae/Capromyinae, Heteropsomyinae, and Heptaxodontidae). Few species survive today, and the evolutionary affinities of living and extinct Caribbean caviomorphs to each other and to mainland taxa are unclear: Are they monophyletic, polyphyletic, or paraphyletic? We use ancient DNA techniques to present the first genetic data for extinct heteropsomyines and heptaxodontids, as well as for several extinct capromyids, and demonstrate through analysis of mitogenomic and nuclear data sets that all sampled Caribbean caviomorphs represent a well-supported monophyletic group. The remarkable morphological and ecological variation observed across living and extinct caviomorphs from Cuba, Hispaniola, Jamaica, Puerto Rico, and other islands was generated through within-archipelago evolutionary radiation following a single Early Miocene overwater colonization. This evolutionary pattern contrasts with the origination of diversity in many other Caribbean groups. All living and extinct Caribbean caviomorphs comprise a single biologically remarkable subfamily (Capromyinae) within the morphologically conservative living Neotropical family Echimyidae. Caribbean caviomorphs represent an important new example of insular mammalian adaptive radiation, where taxa retaining “ancestral-type” characteristics coexisted alongside taxa occupying novel island niches. Diversification was associated with the greatest insular body mass increase recorded in rodents and possibly the greatest for any mammal lineage.     

Three New Ranidae Mitogenomes and the Evolution of Mitochondrial Gene Rearrangements among Ranidae Species

Various types of gene rearrangements have been discovered in the mitogenoes of the frog family Ranidae.In this study,we determined the complete mitogenome sequence of three Rana frogs.By combining the available mitogenomic data sets from Gen Bank,we evaluated the phylogenetic relationships of Ranidae at the mitogenome level and analyzed mitogenome rearrangement cases within Ranidae.The three frogs shared an identical mitogenome organization that was extremely similar to the typical Neobatrachian-type arrangement.Except for the genus Babina,the monophyly of each genus was well supported.The genus Amnirana occupied the most basal position among the Ranidae.The[Lithobates+Rana]was the closest sister group of Odorrana.The diversity of mitochondrial gene arrangements in ranid species was unexpectedly high,with 47 mitogenomes from 40 ranids being classified into 10 different gene rearrangement types.Some taxa owned their unique gene rearrangement characteristics,which had significant implication for their phylogeny analysis.All rearrangement events discovered in the Ranidae mitogenomes can be explained by the duplication and random loss model.     

Phylomitogenomic analyses strongly support the sister relationship of the Chaetognatha and Protostomia

The phylum Chaetognatha (arrow worms) comprises a group of small marine predators that constitute a critical component of the zooplankton community throughout the world's oceans. Various phylogenetic affiliations have been proposed for the Chaetognatha, for which there are at least nine possible phylogenetic positions. Resolving the phylogenetic position of the chaetognaths is a key in understanding the fundamental developmental features of bilaterians. In comparison with the typical gene content of metazoan mitogenomes, two protein‐coding genes (atp6 and atp8) are absent from all chaetognaths. The two mitogenomes sequenced from Sagitta crassa and Zonosagitta nagae in this study nevertheless contain two and four tRNA genes, respectively, in contrast to those of the other five chaetognaths reported where only one tRNA gene (trnMet) is present, thus invalidating the view that all chaetognath mitogenomes have a single tRNA gene. A conserved major gene order shared by all chaetognaths could be partially identified in many protostome mitogenomes, but not in any ancestral mitogenome gene arrangement of the four deuterostome groups. Phylogenetic analysis of the deduced amino acid sequences of protein‐coding genes from 85 mitogenomes of 19 groups suggests the Chaetognatha to be a sister group to the protostomes, a result consistent with evidences from the developmental pattern and other molecular analyses.