The mitochondrial genome of the Russian wheat aphid Diuraphis noxia: Large repetitive sequences between trnE and trnF in aphids

To characterize aphid mitochondrial genome (mitogenome) features, we sequenced the complete mitogenome of the Russian wheat aphid, Diuraphis noxia. The 15,784-bp mitogenome with a high A + T content (84.76%) and strong C skew (− 0.26) was arranged in the same gene order as that of the ancestral insect. Unlike typical insect mitogenomes, D. noxia possessed a large tandem repeat region (644 bp) located between trnE and trnF. Sequencing partial mitogenome of the cotton aphid (Aphis gossypii) further confirmed the presence of the large repeat region in aphids, but with different repeat length and copy number. Another motif (58 bp) tandemly repeated 2.3 times in the control region of D. noxia. All repeat units in D. noxia could be folded into stem-loop secondary structures, which could further promote an increase in copy numbers. Characterization of the D. noxia mitogenome revealed distinct mitogenome architectures, thus advancing our understanding of insect mitogenomic diversities and evolution.     

Characterization of the complete mitochondrial genome of the giant silkworm moth, Eriogyna pyretorum (Lepidoptera: Saturniidae)

The complete mitochondrial genome (mitogenome) of Eriogyna pyretorum (Lepidoptera: Saturniidae) was determined as being composed of 15,327 base pairs (bp), including 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a control region. The arrangement of the PCGs is the same as that found in the other sequenced lepidopteran. The AT skewness for the E. pyretorum mitogenome is slightly negative (-0.031), indicating the occurrence of more Ts than As. The nucleotide composition of the E. pyretorum mitogenome is also biased toward A + T nucleotides (80.82%). All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 and 2 (cox1 and cox2). Two of the 13 PCGs harbor the incomplete termination codon by T. All tRNA genes have a typical clover-leaf structure of mitochondrial tRNA, with the exception of trnS1(AGN) and trnS2(UCN). Phylogenetic analysis among the available lepidopteran species supports the current morphology-based hypothesis that Bombycoidea, Geometroidea, Notodontidea, Papilionoidea and Pyraloidea are monophyletic. As has been previously suggested, Bombycidae (Bombyx mori and Bombyx mandarina), Sphingoidae (Manduca sexta) and Saturniidae (Antheraea pernyi, Antheraea yamamai, E. pyretorum and Caligula boisduvalii) formed a group.     

Complete mitochondrial genome of Malaysian Mahseer (Tor tambroides)

This is the first documentation of the complete mitochondrial genome sequence of the Malaysian Mahseer, Tor tambroides. The 16,690 bp mitogenome with GenBank accession number JX444718 contains 13 protein genes, 22 tRNAs, two rRNAs, and a noncoding control region (D-loop) as is typical of most vertebrates. The phylogenomic reconstruction of this newly generated data with 21 Cypriniformes GenBank accession ID concurs with the recognized status of T. tambroides within the subfamily Cyprininae. This is in agreement with previous hypotheses based on morphological and partial mitochondrial analyses.     

Jumping and gliding rodents: Mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes + Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.     

The complete mitochondrial genome of the fall webworm,Hyphantria cunea (Lepidoptera: Arctiidae)

The complete mitochondrial genome (mitogenome) of the fall webworm, Hyphantria cunea (Lepidoptera: Arctiidae) was determined. The genome is a circular molecule 15 481 bp long. It presents a typical gene organization and order for completely sequenced lepidopteran mitogenomes, but differs from the insect ancestral type for the placement of tRNA^Met. The nucleotide composition of the genome is also highly A + T biased, accounting for 80.38%, with a slightly positive AT skewness (0.010), indicating the occurrence of more As than Ts, as found in the Noctuoidea species. All protein-coding genes (PCGs) are initiated by ATN codons, except for COI, which is tentatively designated by the CGA codon as observed in other lepidopterans. Four of 13 PCGs harbor the incomplete termination codon, T or TA. All tRNAs have a typical clover-leaf structure of mitochondrial tRNAs, except for tRNA^Ser(AGN), the DHU arm of which could not form a stable stem-loop structure. The intergenic spacer sequence between tRNA^Ser(AGN) and ND1 also contains the ATACTAA motif, which is conserved across the Lepidoptera order. The H. cunea A+T-rich region of 357 bp is comprised of non-repetitive sequences, but harbors several features common to the Lepidoptera insects, including the motif ATAGA followed by an 18 bp poly-T stretch, a microsatellite-like (AT)₈ element preceded by the ATTTA motif, an 11 bp poly-A present immediately upstream tRNA^Met. The phylogenetic analyses support the view that the H. cunea is closerly related to the Lymantria dispar than Ochrogaster lunifer, and support the hypothesis that Noctuoidea (H. cunea, L. dispar, and O. lunifer) and Geometroidea (Phthonandria atrilineata) are monophyletic. However, in the phylogenetic trees based on mitogenome sequences among the lepidopteran superfamilies, Papillonoidea (Artogeia melete, Acraea issoria, and Coreana raphaelis) joined basally within the monophyly of Lepidoptera, which is different to the traditional classification.     

Characterization of the complete mitochondrial genome of Diaphania pyloalis (Lepidoptera: Pyralididae)

The complete mitochondrial genome (mitogenome) of Diaphania pyloalis (Lepidoptera: Pyralididae) was determined to be 15,298 bp and has the typical gene organization of mitogenomes from lepidopteran insects. It consists of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The A + T content of this mitogenome is 80.83% and the AT skew is slightly positive. All PCGs are initiated by ATN codons, except for cytochrome c oxidase subunit 1 (cox1) gene which is initiated by CGA. Only the cox2 gene has an incomplete stop codon consisting of just a T. All the tRNA genes display a typical clover-leaf structure of mitochondrial tRNA. The A + T-rich region of the mitogenome is 332 bp in length, including several common features found in lepidopteran mitogenomes. Phylogenetic analysis showed that the D. pyloalis is close to Pyralididae.     

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.     

Moramonas marocensis gen. nov., sp. nov.: a jakobid flagellate isolated from desert soil with a bacteria-like,but bloated mitochondrial genome

A new jakobid genus has been isolated from Moroccan desert soil. The cyst-forming protist Moramonas marocensis gen. nov., sp. nov. has two anteriorly inserted flagella of which one points to the posterior cell pole accompanying the ventral feeding groove and is equipped with a dorsal vane—a feature typical for the Jakobida. It further shows a flagellar root system consisting of singlet microtubular root, left root (R1), right root (R2) and typical fibres associated with R1 and R2. The affiliation of M. marocensis to the Jakobida was confirmed by molecular phylogenetic analyses of the SSU rRNA gene, five nuclear genes and 66 mitochondrial protein-coding genes. The mitochondrial genome has the high number of genes typical for jakobids, and bacterial features, such as the four-subunit RNA polymerase and Shine–Dalgarno sequences upstream of the coding regions of several genes. The M. marocensis mitochondrial genome encodes a similar number of genes as other jakobids, but is unique in its very large genome size (greater than 264 kbp), which is three to four times higher than that of any other jakobid species investigated yet. This increase seems to be due to a massive expansion in non-coding DNA, creating a bloated genome like those of plant mitochondria.     

The complete mitochondrial genome sequence of the little egret (Egretta garzetta)

Many phylogenetic questions in the Ciconiiformes remain unresolved and complete mitogenome data are urgently needed for further molecular investigation. In this work, we determined the complete mitogenome sequence of the little egret (Egretta garzetta). The genome was 17,361 bp in length and the gene organization was typical of other avian mtDNA. In protein-coding genes (PCGs), a C insertion was found in ND3, and COIII and ND4 terminated with incomplete stop codons (T). tRNA-Val and tRNA-Ser (AGY) were unable to fold into canonical cloverleaf secondary structures because they had lost the DHU arms. Long repetitive sequences consisting of five types of tandem repeats were found at the 3′ end of Domain III in the control region. A phylogenetic analysis of 11 species of Ciconiiformes was done using complete mitogenome data and 12 PCGs. The tree topologies obtained with these two strategies were identical, which strongly confirmed the monophyly of Ardeidae, Threskiorothidae and Ciconiidae. The phylogenetic analysis also revealed that Egretta was more closely related to Ardea than to Nycticorax in the Ardeidae, and Platalea was more closely related to Threskiornis than to Nipponia in the Threskiornithidae. These findings contribute to our understanding of the phylogenetic relationships of Ciconiiformes based on complete mitogenome data.     

Complete mitochondrial genome sequence of bighead croaker Collichthys niveatus (Perciformes, Sciaenidae): a mitogenomic perspective on the phylogenetic relationships of Pseudosciaeniae

The monophyly and phylogenetic relationships of Pseudosciaeniae have long been controversial. Here we describe the mitochondrial genome (mitogenome) sequence of Collichthys niveatus. It is a circular double-stranded DNA molecule of 16,450 base pairs (bp) in length with a standard set of 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes as well as a non-coding control region. The mitogenome of C. niveatus shared common features with those of other bony fishes in terms of gene arrangement, base composition, and tRNA structures. The C. niveatus mitogenome exhibited pronounced strand-specific asymmetry in nucleotide composition, which was also reflected in the codon usage of genes oriented in opposite directions. Contrary to the typical structure of the control region, the central conserved blocks (CSB-D, -E, and -F) could not be detected in C. niveatus mitogenome. Phylogenetic analysis based on whole mitogenome sequences provided strong support for the monophyly of Pseudosciaeniae, and sister-group relationships of C. niveatus + Collichthys lucidus and Larimichthys crocea + Larimichthys polyactis, which was consistent with the traditional taxonomy. Unexpected divergence was found in two C. niveatus mitogenomes and several hypotheses were proposed to explain this observation including misidentification and introgressive hybridization between C. niveatus and L. polyactis, and polyphyletic origin of C. niveatus. We considered species misidentification to be the main hypothesis. However, additional data is essential to test these proposed hypotheses.     

The complete mitochondrial genome sequence of the world's largest fish,the whale shark (Rhincodon typus), and its comparison with those of related shark species

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a “vulnerable” species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species.     

Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla coelebs): New inferences about mitochondrial rate heterogeneity,neutral theory,and phylogenetic relationships within the order Passeriformes

We describe whole mitochondrial genome sequences from four subspecies of the common chaffinch (Fringilla coelebs), and compare them to 31 publicly available mitochondrial genome sequences from other Passeriformes. Rates and patterns of mitochondrial gene evolution are analyzed at different taxonomic levels within this avian order, and evidence is adduced for and against the nearly neutral theory of molecular evolution and the role of positive selection in shaping genetic variation of this small but critical genome. We find evidence of mitochondrial rate heterogeneity in birds as in other vertebrates, likely due to differences in mutational pressure across the genome. Unlike in gadine fish and some of the human mitochondrial work we do not observe strong support for the nearly neutral theory of molecular evolution; instead evidence from molecular clocks, distribution of dN/dS ratios at different levels of the taxonomic hierarchy and in different lineages, McDonald–Kreitman tests within Fringillidae, and site-specific tests of selection within Passeriformes, all point to a role for positive selection, especially for the complex I NADH dehydrogenase genes. The protein-coding mitogenome phylogeny of the order Passeriformes is broadly consistent with previously-reported molecular findings, but provides support for a sister relationship between the superfamilies Muscicapoidea and Passeroidea on a short basal internode of the Passerida where relationships have been difficult to resolve. An unexpected placement of the Paridae (represented by Hume's groundpecker) within the Muscicapoidea was observed. Consistent with other molecular studies the mtDNA phylogeny reveals paraphyly within the Muscicapoidea and a sister relationship of Fringilla with Carduelis rather than Emberiza.     

The complete mitochondrial genome of the sycamore lace bug Corythucha ciliata (Hemiptera: Tingidae)

The complete mitochondrial genome of the sycamore lace bug, Corythucha ciliata, was sequenced in this study. It represents the first sequenced mitogenome of family Tingidae in Heteroptera. The mitogenome of C. ciliata is 15,257 bp and contains 37 genes including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a large non-coding region. Gene arrangement, nucleotide content, codon usage, and amino acid composition and asymmetry indicate a high degree of conservation with six other species of Cimicomorpha. The 13 PCGs initiated with ATN as the start codon and terminated with TAA, TA or T as stop codon. The evolutionary rate of each PCG was different, among which ATP8 showed the highest rate while ATP6 indicated the lowest rate. The 22 tRNAs genes apparently fold into a typical cloverleaf structure; however, the anticodon (TTC) of trnSer (AGN) differs from other Heteropteran insects. Secondary structure modeling of rRNA genes revealed similarity to other insects, except for two incomplete helices (H1648 and H2735) in lrRNA. The predicted secondary structure of lrRNA indicates 45 helices in six domains, whereas srRNA has 27 helices in three domains. Three potential stem–loops and two tandem repeats (–TCTAAT–) were identified in the A+T-rich region. Phylogenetic analysis indicated that C. ciliata is a sister group to other Heteroptera species based on analysis of the 13 PCGs.     

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A + T content of the complete mitogenome is 80.83%. The AT skew value ([A − T] / [A + T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of “ATAGA” + 19-bp poly T. It also contains a microsatellite-like (TA)₂₆, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.     

Characterization of the complete mitochondrial genome and a set of polymorphic microsatellite markers through next-generation sequencing for the brown brocket deer Mazama gouazoubira

The complete mitochondrial genome of the brown brocket deer Mazama gouazoubira and a set of polymorphic microsatellite markers were identified by 454-pyrosequencing. De novo genome assembly recovered 98% of the mitochondrial genome with a mean coverage of 9-fold. The mitogenome consisted of 16,356 base pairs that included 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and the control region, as found in other deer. The genetic divergence between the mitogenome described here and a previously published report was ∼0.5%, with the control region and ND5 gene showing the highest intraspecific variation. Seven polymorphic loci were characterized using 15 unrelated individuals; there was moderate genetic variation across most loci (mean of 5.6 alleles/locus, mean expected heterozygosity = 0.70), with only one locus deviating significantly from Hardy-Weinberg equilibrium, probably because of null alleles. Marker independence was confirmed with tests for linkage disequilibrium. The genetic variation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the phylogeography and population genetic patterns in M. gouazoubira, particularly in the context of habitat fragmentation in South America.     

Complete mitochondrial DNA sequence of the endangered fish (Bahaba taipingensis): Mitogenome characterization and phylogenetic implications

To understand the systematic status of Bahaba taipingensis within Sciaenidae, the complete mitochondrial genome (mitogenome) sequence of Chinese bahaba has recently been determined by long PCR and primer walking methods. The complete mitochondrial genome is 16500 bp in length and contains 37 mitochondrial genes (13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes) as well as a control region (CR) as other bony fishes. Within the control region, we identified the extended termination associated sequence domain (ETAS), the central conserved sequence block domain (CSB-D, SCB-E and CSB-F) and the conserved sequence block domain (CSB-1, CSB-2 and CSB-3). Phylogenetic analyses revealed that Bahaba taipingensis is more closely related to Pseudosciaeniae than Argyrosominae and Sciaeninae. Additionally, Bahaba taipingensis is the sister taxon of Miichthys miiuy, and those two are sister to Collichthys plus Larimichthys.     

Complete mitochondrial exploration of Echinococcus multilocularis from French alveolar echinococcosis patients

Alveolar echinococcosis (AE) is a parasitosis that is expanding worldwide, including in Europe. The development of genotypic markers is essential to follow its spatiotemporal evolution. Sequencing of the commonly used mitochondrial genes cob, cox1, and nad2 shows low discriminatory power, and analysis of the microsatellite marker EmsB does not allow nucleotide sequence analysis. We aimed to develop a new method for the genotyping of Echinococcus multilocularis based on whole mitochondrial genome (mitogenome) sequencing, to determine the genetic diversity among 30 human visceral samples from French patients, and compare this method with those currently in use. Sequencing of the whole mitochondrial genome was carried out after amplification by PCR, using one uniplex and two multiplex reactions to cover the 13,738 bp of the mitogenome, combined with Illumina technology. Thirty complete mitogenome sequences were obtained from AE lesions. One showed strong identity with Asian genotypes (99.98% identity) in a patient who had travelled to China. The other 29 mitogenomes could be differentiated into 13 haplotypes, showing higher haplotype and nucleotide diversity than when using the cob, cox1, and nad2 gene sequences alone. The mitochondrial genotyping data and EmsB profiles did not overlap, probably because one method uses the mitochondrial genome and the other the nuclear genome. The pairwise fixation index (F_st) value between individuals living inside and those living outside the endemic area was high (F_st = 0.222, P = 0.002). This is consistent with the hypothesis of an expansion from historical endemic areas to peripheral regions.     

The complete mitochondrial genome of the intertidal copepod Tigriopus sp. (Copepoda, Harpactidae) from Korea and phylogenetic considerations

This paper reports on the basic characteristics of the mitochondrial genome of the intertidal copepod Tigriopus sp. from Korea, including its structural organization, base composition of rRNAs and protein-encoding genes, and the secondary structure of tRNAs. We amplified the complete mitochondrial DNA of the intertidal copepod Tigriopus sp. from Korea (sampling site Busan) by long-polymerase chain reaction (long-PCR) with conserved primers and sequenced this mitogenome by primer walking using flanking sequences as sequencing primers. The primer informations were obtained as expressed sequence tags (ESTs) from Tigriopus sp. The resultant Tigriopus sp. mitochondrial DNA sequence was 14,301 bp with a conserved structural organization, compared to that of T. japonicus from Japan with significant differences in several protein-coding regions including rRNAs, although the genomic organization of the mitochondrial genome was identical. In order to investigate biogeographic differences within the genus Tigriopus, we analyzed the CO1 gene by sequencing. This way, we compared several Tigriopus species from Korea, Japan, Hong Kong and Taiwan as well as other related species such as T. californicus, T. brevicornis and T. fulvus. The results further support the notion that the copepods display significantly different genomes within the same genus. These findings provide valuable genomic information for further studies on the population genetics and speciation processes within the genus Tigriopus.     

Characterisation of the complete mitochondrial genome and 13 microsatellite loci through next-generation sequencing for the New Caledonian spider-ant Leptomyrmex pallens

The complete mitochondrial genome and a set of polymorphic microsatellite markers were identified by 454 pyrosequencing (1/16th of a plate) for the New Caledonian rainforest spider-ant Leptomyrmex pallens. De novo genome assembly recovered the entire mitochondrial genome with mean coverage of 8.9-fold (range 1–27). The mitogenome consists of 15,591 base pairs including 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a non-coding AT-rich region. The genome arrangement is typical of insect taxa and very similar to the only other published ant mitogenome from the Solenopsis genus, with the main differences consisting of translocations and inversions of tRNAs. A total of 13 polymorphic loci were also characterized using 41 individuals from a single population in the Aoupinié region, corresponding to workers from 21 nests and 16 foraging workers. We observed moderate genetic variation across most loci (mean number of alleles per locus = 4.50; mean expected heterozygosity = 0.53) with evidence of only two loci deviating significantly from Hardy–Weinberg equilibrium due to null alleles. Marker independence was confirmed with tests for linkage disequilibrium. Most loci cross amplified for three additional Leptomyrmex species. The annotation of the mitogenome and characterization of microsatellite markers will provide useful tools for assessing the colony structure, population genetic patterns, and dispersal strategy of L. pallens in the context of rainforest fragmentation in New Caledonia. Furthermore, this paper confirms a recent line of evidence that comprehensive mitochondrial data can be obtained relatively easily from small next-generation sequencing analyses. Greater synthesis of next-generation sequencing data will play a significant role in expanding the taxonomic representation of mitochondrial genome sequences.