Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A + T composition (79.3%–80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA^Gln and ND2 (64–87 bp), conserved ATACTAA motif between tRNA^{Ser (UCN)} and ND1, and characteristic features of the A + T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA^Asn/tRNA^{Ser (AGN)}, instead of the rearrangement tRNA^{Ser (AGN)}/tRNA^Asn, found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea + Noctuoidea + Geometroidea) + Pyraloidea) + Papilionoidea) + Tortricoidea) + Yponomeutoidea) + Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae + (Hesperiidae + (Pieridae + ((Lycaenidae + Riodinidae) + Nymphalidae)))) were strongly supported in all analyses (0.98–1 by BI and 96–100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.     

Characterization of four mitochondrial genomes of Crambidae (Lepidoptera,Pyraloidea) and phylogenetic implications

Loxostege turbidalis, Loxostege aeruginalis, Pyrausta despicata, and Crambus perlellus belong to Crambidae, Pyraloidea. Their mitochondrial genomes (mitogenomes) were successfully sequenced. The mitogenomes of L. turbidalis, L. aeruginalis, P. despicata, and C. perlellus are 15 240 bp, 15 339 bp, 15 389 bp, and 15 440 bp. The four mitogenomes all have a typical insect mitochondrial gene order, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and one A + T rich region (control region). The PCGs are initiated by the typical ATN codons, except CGA for the cox1 gene. Most PCGs terminate with common codon TAA or TAG, the incomplete codon T is found as the stop codon for cox2, nad4, and nad5. Most tRNA genes exhibit typical cloverleaf structure, except trnS1 (AGN) lacking the dihydrouridine (DHU) arm. The secondary structure of rRNA of four mitogenomes were predicted. Poly-T structure and micro-satellite regions are conserved in control regions. The phylogenetic analyses based on 13 PCGs showed the relationships of subfamilies in Pyraloidea. Pyralidae, and Crambidae are monophyletic, respectively. Pyralidae comprises four subfamilies, which form the following topology with high support values: (Galleriinae + ((Pyralinae + Epipaschiinae)+ Phycitinae)). Crambidae includes seven subfamilies and is divided into two lineages. Pyraustinae and Spilomelinae are sister groups of each other, and form the “PS clade.” Other five subfamilies (Crambinae, Acentropinae, Scopariinae, Schoenobiinae, and Glaphyriinae) form the “non-PS clade” in the Bayesian inference tree. However, Schoenobiinae is not grouped with the other four subfamilies and located at the base of Crambidae in two maximum likelihood trees.     

The mitochondrial genomes of three skippers: Insights into the evolution of the family Hesperiidae (Lepidoptera)

We sequenced the mitogenomes of Astictopterus jama, Isoteinon lamprospilus and Notocrypta curvifascia to obtain further insight into the mitogenomic architecture evolution and performed phylogenetic reconstruction using 29 Hesperiidae mitogenome sequences. The complete mitogenome sequences of A. jama, I. lamprospilus and N. curvifascia are 15,430, 15,430 and 15,546 bp in size, respectively. All contain 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and an A + T-rich region. Nucleotide composition is A + T biased, and the majority of the protein-coding genes exhibit a negative AT-skew, which is reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region is comprised of nonrepetitive sequences, including the motif ATAGA followed by a poly-T stretch, a microsatellite-like element next to the ATTTA motif, and a poly-A adjacent to tRNAs. Although most genes evolve under a strong purifying selection, the entire nad gene family (especially nad6) exhibits somewhat relaxed purifying selection, and atp8, evolving under a highly relaxed selection, is an outlier in the family Hesperiidae. Several different approaches relatively consistently indicated that nad6, atp8 and nad4 are comparatively fast-evolving genes in this family, which may have implications for future phylogenetic, population genetics and species diagnostics studies. For phylogenetic analyses of Hesperiidae, we tested a few datasets, and found that the one comprising all 37 genes produced the highest node support, indicating that the inclusion of RNAs improves the phylogenetic signal. Results indicate that subfamilies Euschemoninae, Heteropterinae, and Coeliadinae are monophyletic with strong nodal support, but Pyrginae and Eudaminae are paraphyletic. Finally, we confirm that A. jama and I. lamprospilus are close relatives.     

Complete mitochondrial genome sequence of Bonasa sewerzowi(Galliformes: Phasianidae) and phylogenetic analysis

Bonasa sewerzowi, the smallest and most southerly distributed grouse species in the world, is a bird endemic to China. The population of B. sewerzowi had shown a declining trend, which made it to be the endangered species in the China Red Data Book and Category I of nationally protected animals. So far, however, most studies about this species were mainly focused on the morphological and ecological aspects. In order to further study the feature of B. sewerzowi, the complete mitochondrial genome(mitogenome) of B. sewerzowi was sequenced by Illumina Hiseq 2000 high-throughput sequencing. Then, we focused on comparative genomics of two Bonasa species to find their characteristics. Finally, phylogenetic position of Bonasa was made based on the mitogenome dataset. Our results revealed that:(1) the mitogenome of B. sewerzowi, consisting of 16 658 bp, displayed typical genome organization and gene order found in other previously determined Galliformes mitogenomes;(2) the structure and composition of mitogenomes were similar between B. sewerzowi and B. bonasia;(3) the monophyly of Bonasa was well supported, which had a closer phylogenetic relationship with Meleagris gallopavo.     

Complete mitochondrial genomes of three Oxya grasshoppers (Orthoptera) and their implications for phylogenetic reconstruction

Oxya is a genus of grasshoppers (Orthoptera: Acridoidea) attacking rice and other gramineous plants in Africa and Asia. In the present study, we characterized complete mitochondrial genomes (mitogenomes) of three species, Oxya japonica japonica (15,427 bp), Oxya hainanensis (15,443 bp) and Oxya agavisa robusta (15,552 bp) collected from China. The three mitogenomes contained a typical gene set of metazoan mitogenomes and shared the same gene order with other Acridid grasshoppers, including the rearrangement of tRNA^Asp and tRNA^Lys. Analyses of pairwise genetic distances showed that ATP8 was the least conserved gene, while COI the most conserved. To determine the position of Oxya grasshoppers in the phylogeny of Acrididae, we reconstructed phylogenetic trees among 64 species from across 11 subfamilies using nucleotide sequences of mitogenomes. While the tree confirms traditional classifications of Acrididae at major higher-levels, it suggests a few modifications for classifications at lower-levels.     

Evolutionary implications of analyses of complete mitochondrial genomes across order Zoantharia (Cnidaria: Hexacorallia)

Cnidarians are early-diverging metazoans, but evolutionary aspects of some taxa are still poorly understood, as in the order Zoantharia (Anthozoa: Hexacorallia). Zoantharians have been divided into two suborders based on the arrangement of the fifth septae as complete (Macrocnemina) or incomplete (Brachycnemina). Previous molecular phylogenetic analyses have indicated the need for re-evaluation as Macrocnemina has been found to be paraphyletic. Despite many phylogenetic studies, the recovery of complete mitochondrial genomes (mt-genomes) for systematic and evolutionary studies of zoantharians has been limited. The present study represents the first to sequence the complete mt-genomes of members of eight of nine zoantharian families. Although all examined mt-genomes had the same gene order arrangement, there were variations among mt-genomes' sizes, nucleotide substitution rates, and introns. Only two species did not have the cox1 intron, which harbors a gene coding a homing endonuclease of the LAGLIDADG type. Our mitogenomic analyses also showed relatively high nucleotide diversity in mt-DNA regions other than the standard regions traditionally considered for DNA barcoding of this group. Phylogenetic analyses using 13 mt-genome protein-coding genes recovered a fully resolved tree with clear separation between macrocnemic representatives. Ancestral state reconstruction analyses revealed three main transitions in arrangement of the marginal musculature through the evolutionary history of the order. An “early” transition from reticulate mesogleal to a cteniform endodermal arrangement was followed by transitions that occurred in the common ancestor of the Brachycnemina and family Hydrozoanthidae. Our results indicate the need for clarification of higher-level phylogeny and taxonomy of Zoantharia.     

Complete mitochondrial genomes of eleven extinct or possibly extinct bird species

Natural history museum collections represent a vast source of ancient and historical DNA samples from extinct taxa that can be utilized by high‐throughput sequencing tools to reveal novel genetic and phylogenetic information about them. Here, we report on the successful sequencing of complete mitochondrial genome sequences (mitogenomes) from eleven extinct bird species, using de novo assembly of short sequences derived from toepad samples of degraded DNA from museum specimens. For two species (the Passenger Pigeon Ectopistes migratorius and the South Island Piopio Turnagra capensis), whole mitogenomes were already available from recent studies, whereas for five others (the Great Auk Pinguinis impennis, the Imperial Woodpecker Campehilus imperialis, the Huia Heteralocha acutirostris, the Kauai Oo Moho braccathus and the South Island Kokako Callaeas cinereus), there were partial mitochondrial sequences available for comparison. For all seven species, we found sequence similarities of >98%. For the remaining four species (the Kamao Myadestes myadestinus, the Paradise Parrot Psephotellus pulcherrimus, the Ou Psittirostra psittacea and the Lesser Akialoa Akialoa obscura), there was no sequence information available for comparison, so we conducted blast searches and phylogenetic analyses to determine their phylogenetic positions and identify their closest extant relatives. These mitogenomes will be valuable for future analyses of avian phylogenetics and illustrate the importance of museum collections as repositories for genomics resources.     

Tsukiyamaia,a new genus of the tribe Baorini (Lepidoptera,Hesperiidae, Hesperiinae)

Skippers of the tribe Baorini are evidently a monophyletic group in the subfamily Hesperiinae. In this study, a new Baorini member Tsukiyamaia albimacula gen. n. et sp. n. is described from north Myanmar, southwest China and north Vietnam. Despite its peculiar and striking wing-pattern, this new genus has some important characters of Baorini, such as a broad and bifid uncus and a well-developed gnathos. Based on an analysis of male genitalia and the molecular phylogenies inferred from both mitochondrial and nuclear genes (28 taxa, total aligned length: 2968 bp), it is proposed that the genus Tsukiyamaia is closely related to the genus Polytremis, which has high species diversity in China. This study not only describes a new skipper but also highlights that Tsukiyamaia is important in clarifying phylogenetic relationship of Polytremis and its allies.     

Complete mitochondrial genomes of Conogethes punctiferalis and C. pinicolalis (Lepidoptera: Crambidae): Genomic comparison and phylogenetic inference in Pyraloidea

Conogethes punctiferalis Guenée, 1854 (Lepidoptera: Crambidae) is a serious polyphagous pest that attacks more than one hundred species of plants. Previously, C. punctiferalis was determined to be composed of two ecotypes; later, Conogethes pinicolalis, was described as a separate species. Due to the prolonged negligence of C. pinicolalis as an independent species, the genetic perspective of the two species is limited. Thus, in this study, 15,332 and 15,336 bp-long complete mitochondrial genome (mitogenome) of the two species were sequenced and compared to each other and to 54 available mitogenomes of Pyraloidea. The comparison of each protein-coding gene (PCG) and rRNA gene of the two congeneric species showed substantial sequence divergence, ranging from 3.13% (ATP8) to 8.3% (COIII), with an average of 5.92%. Phylogenetic analyses using concatenated sequences of 13 PCGs and 2 rRNAs (12,458 bp including gaps), both by maximum likelihood (ML) and Bayesian inference (BI) methods, consistently supported the monophyly of each family (Crambidae and Pyralidae) and subfamily, generally with the highest nodal supports. The subfamilial relationships of ((((Acentropinae + Schoenobiinae) + (Scopariinae + Crambinae)) + Evergestinae) + (Spilomelinae + Pyraustinae)) in Crambidae and the subfamilial relationships of ((((Pyralinae + Epipaschiinae) + Phycitinae) + Galleriilinae)) in Pyralidae were obtained in both analyses. However, nodal supports were substantially low in this study, mainly due to limited taxa.     

Complete mitochondrial genomes of three lizard species and the systematic position of the Lacertidae (Squamata)*

Complete mitochondrial (mt) genomes were sequenced from representatives of three lacertid lizards: Podarcis siculus, Podarcis muralis and Phoenicolacerta kulzeri. In all three genomes the arrangement of the 22 tRNAs, the two rRNAs and the 13 protein‐coding genes conforms to the common vertebrate arrangement. The phylogenetic position of Lacertidae within the order Squamata was determined through sequence analyses based on large sections of complete mt genomes. The number of nucleotide sites used for tree construction was 9234 when outgroup taxa were included, and 10 499 when only Squamata were compared. The phylogenetic analyses confirmed the sister group relationship between Lacertidae and Amphisbaenia as previously proposed on the basis of molecular data. Additionally, Bayesian analysis revealed a well supported clade comprising (Gekkonidae (Lacertidae + Amphisbaenia)), which is not in accordance with the traditional morphological view and most of the previous molecular studies. It confirms, however, the close relationship between Gekkonidae and Amphisbaenia as revealed in a recent study based on complete mt genomes from a smaller number of taxa. Intra‐ and intergeneric sequence comparisons of six commonly used marker genes showed rather high levels of divergence within the Lacertidae. In the intrageneric comparison the control region proved to be considerably more conserved than the protein coding genes.     

Complete mitochondrial genome sequence and phylogenetic position of the Amu Darya sturgeon,Pseudoscaphirhynchus kaufmanni (Acipenseriformes: Acipenseridae)

In this study, we successfully assembled the complete mitochondrial genome of the Amu Darya sturgeon Pseudoscaphirhynchus kaufmanni. Based on this mitochondrial genome and previously published mitochondrial genomes of members of the Acipenseridae family, we assessed the phylogenetic position of P. kaufmanni using maximum likelihood and Bayesian inference for phylogeny reconstruction. The resultant phylogenetic trees were well-resolved, with congruence between different phylogenetic methods. This robust phylogenetic analysis elucidated the relationship among the four acipenserid genera and strongly supported the division of the family into three main clades. Evaluation of molecular phylogeny using maximum likelihood and Bayesian analysis led to the following conclusions: (a) the most basal position within the Acipenseridae remains in the clade containing Acipenser oxyrinchus and Acipenser sturio; (b) the genus Scaphirhynchus belongs to the Atlantic clade and is a sister group of the remaining species of the clade; and (c) the close relationship between P. kaufmanni and Acipenser stellatus is well supported.     

Complete mitochondrial genome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae) and the phylogenetic relationship of Saturniidae species

Mitochondrial genome (mitogenome) can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. In this study, we present the complete mitogenome of the atlas moth, Attacus atlas (Lepidoptera: Saturniidae), a well-known silk-producing and ornamental insect with the largest wing surface area of all moths. The mitogenome of A. atlas is a circular molecule of 15,282 bp long, and its nucleotide composition shows heavily biased towards As and Ts, accounting for 79.30%. This genome comprises 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A + T-rich region. It is of note that this genome exhibits a slightly positive AT skew, which is different from the other known Saturniidae species. All PCGs are initiated by ATN codons, except for COI with CGA instead. Only six PCGs use a common stop codon of TAA or TAG, whereas the remaining seven use an incomplete termination codon T or TA. All tRNAs have the typical clover-leaf structure, with an exception for tRNA^Ser(AGN). The A. atlas A + T-rich region contains non-repetitive sequences, but harbors several features common to the Bombycoidea insects. The phylogenetic relationships based on Maximum Likelihood method provide a well-supported outline of Saturniidae, which is in accordance with the traditional morphological classification and recent molecular works.     

Antennal sensilla of male Lophocorona pediasia Common 1973 and their phylogenetic implications (Lepidoptera: Lophocoronidae)

The male antennae of Lophocorona pediasia reveal seven types of sensilla: multiporous sensilla trichodea, multiporous sensilla basiconica of subtypes 1 and 2, multiporous sensilla auricillica, aporous sensilla styloconica, multiporous sensilla coeloconica and aporous Böhm’s sensilla. All these sensillum types are present in the higher Lepidoptera. The existence of sensilla auricillica, hitherto denied by several authors, is demonstrated in this study. Typically, they are located on the distal edge of the flagellomeres as in most Lepidoptera. Therefore, their so-called absence must no longer be regarded as an autapomorphy of the Lophocoronidae. Their presence is a common characteristic of all «Dacnonypha», which include, apart from the Lophocoronidae, the Eriocraniidae and the Acanthopteroctetidae. Their number in L. pediasia is smaller than in these two families and more in accordance with the condition found in higher Lepidoptera, suggesting a more advanced position of the Lophocoronidae. The presence of very numerous sensilla basiconica is a possible autapomorphy of the family.     

基于线粒体基因组的华蝉族系统发育地位研究(半翅目: 蝉科)(英文)

【目的】本研究旨在明确无鼓膜发音器的华蝉族(Sinosenini)昆虫在蝉总科(Cicadoidea)的系统发育地位。【方法】依据在陕西宁陕采集的合哑蝉Karenia caelatata成虫标本,对华蝉族的合哑蝉K. caelatata线粒体基因组进行测序、注释和生物信息学分析;并与蝉总科其他类群的线粒体基因组进行了比较,然后利用最大似然法(ML)和贝叶斯法(BI)分别构建了蝉总科分子系统发育树。【结果】合哑蝉线粒体基因组长14 960 bp (GenBank登录号: MN922304),其基因组成、蛋白编码基因的核苷酸组成和密码子使用等,与蝉总科其他类群具相似特征。核苷酸多样性分析表明,atp8, nad6和nad2为易变基因,而cox1比较保守。非同义替换率和同义替换率比表明,蝉总科昆虫线粒体基因组进化处于高水平的纯化选择下。系统发育分析结果支持蝉次目(Cicadomorpha)的单系性,该次目3个总科的关系为：角蝉总科Membraciodea+(蝉总科Cicadoidea+沫蝉总科Cercopodidea)。无鼓膜发音器的哑蝉属与蝉亚科(Cicadinae)的蜩蝉族(Dundubiini)相关类群聚在一起,且与寒蝉属Meimuna关系最近;黑蝉族(Cicadatrini)的草蝉属Mogannia和音蝉属Vagitanus则与姬蝉亚科(Cicadettinae)相关类群聚在一起;日宁蝉属Yezoterpnosia并非一个单系群。【结论】华蝉族应从姬蝉亚科转移至蝉亚科并与蜩蝉族(Dundubiini)合并,而黑蝉族应从蝉亚科转移至姬蝉亚科。研究结果为进一步解析具有不同发声机制的蝉科昆虫系统演化提供了新信息。     

The complete mitochondrial genome of the mackerel icefish,Champsocephalus gunnari (Actinopterygii: Channichthyidae), with reference to the evolution of mitochondrial genomes in Antarctic notothenioids

The mackerel icefish (Champsocephalus gunnari Lönnberg, 1905) is a ray‐finned fish living in the Southern Ocean around Antarctica. We sequenced the complete mitochondrial (mt) genome of the mackerel icefish and a segment from cytochrome b to the control region (CR) in 32 individuals. The mt genome of the mackerel icefish was rearranged, containing two nicotinamide adenine dinucleotide (reduced form) dehydrogenase subunit 6 (ND6), two tRNA^Glu, and two CRs. However, variations in numbers of ND6 and tRNA^Glu were observed amongst individuals. These variations included type 1 (containing two ND6 and two tRNA^Glu), type 2 (containing one ND6, one incomplete ND6, and one tRNA^Glu), and type 3 (containing one ND6 and one tRNA^Glu). The gene orders of types 1 and 2, and variations in numbers of ND6 and tRNA^Glu were not previously found in any Antarctic notothenioids, whereas type 3 is the same as that of Racovitzia glacialis. Phylogenetic analyses of CR DNA sequences showed that duplicated CRs of the same species formed a monophyletic group, suggesting that duplication of CRs occurred in each species. The frequent duplication of mt genomes in Antarctic notothenioids is an unusual feature in vertebrates. We propose that interspecific hybridization and impairment of mismatch repair might account for the high frequency of gene duplications and rearrangement of mt genomes in Antarctic notothenioids.     

Phylogeny of Veneridae (Bivalvia) based on mitochondrial genomes

Veneridae is one of the most diverse families of bivalve molluscs. However, their phylogenetic relationships among subfamilies have been debated for years. To explore phylogenetic relationships of Veneridae, we sequenced 13 complete mitochondrial genome sequences from eight subfamilies and compared with available complete mitochondrial genome of other Veneridae taxa (18 previously reported sequences). Phylogenetic analyses using probabilistic methods recovered two highly supported clades. In addition, the protein‐coding gene order revealed a highly conserved pattern among the same subclade lineages. According to our molecular analyses, Tapetinae should be recognized as a valid subfamily, but the genera formed para‐polyphyletic clades. Chioninae was recovered not monophyletic that differs from a previously molecular phylogeny. Furthermore, the reconstructed chronogram calibrated with fossils recovered the Veneridae have originated during the early Permian (about 290 million years ago). Noticeably, programmed frameshift was found in the nad4 gene of Leukoma jedoensis, Anomalodiscus squamosus and Antigona lamellaris and cob gene of L. jedoensis. This is the first time that the presence of the programmed frameshift has been found in the protein‐coding genes of Heterodonta species. Our results improved the phylogenetic resolution within Veneridae, and a more taxonomic sampling analysis of the subfamily Chioninae is supposed to construct.     

Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous–Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India–Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.