Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the bell-ring frog, Buergeria buergeri (family Rhacophoridae)

Rapid and collective determination of the complete "hot-spring frog" mitochondrial genome containing long repeat regions using Nanopore sequencing

The mitochondrial genome (mt-genome) is one of the promising molecular markers for phylogenetics and population genetics. Recently, various mt-genomes have been determined rapidly by using massively parallel sequencers. However, the control region (CR, also called D-loop) in mt-genomes remain difficult to precisely determine due to the presence of repeat regions. Here, using Nanopore sequencing, we succeeded in rapid and collective determination of complete mt-genome of the hot-spring frog, Buergeria japonica, and found that its mt-genome size was 22,274 bp including CR (6,929 bp) with two types of tandem repeat motifs forming repeat regions. Comparison of assembly strategies revealed that the long- and short-read data combined together enabled efficient determination of the CR, but the short-read data alone did not. The B. japonica CR was longer than that of a congenic species inhabiting cooler climate areas, Buergeria buergeri, because of the long repeat regions in the former. During the thermal adaptation of B. japonica, the longer repeat regions in its CR may have accumulated within a period after divergence from B. buergeri.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Characterization of Two Mitogenomes of Hyla sanchiangensis (Anura: Hylidae), with Phylogenetic Relationships and Selection Pressure Analyses of Hylidae

Hyla sanchiangensis (Anura: Hylidae) is endemic to China and is distributed across Anhui, Zhejiang, Fujian, Guangdong, Guangxi, Hunan, and Guizhou provinces. The mitogenomes of H. sanchiangensis from two different sites (Jinxiu, Guangxi, and Wencheng, Zhejiang) were sequenced. Phylogenetic analyses were conducted, including 38 mitogenomes of Hylidae from the NCBI database, and assessed the phylogenetic relationship of H. sanchiangensis within the analyzed dataset. Two mitogenomes of H. sanchiangensis showed the typical mitochondrial gene arrangement with 13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA (tRNA) genes, and one non-coding control region (D-loop). The lengths of the 12S rRNA and 16S rRNA genes from both samples (Jinxiu and Wencheng) were 933 bp and 1604 bp, respectively. The genetic distance (p-distance transformed into percent) on the basis of the mitogenomes (excluding the control region) of the two samples was calculated as 4.4%. Hyla sanchiangensis showed a close phylogenetic relationship with the clade of (H. annectans + H. tsinlingensis), which was supported by ML and BI analyses. In the branch-site model, five positive selection sites were found in the clade of Hyla and Dryophytes: Cytb protein (at position 316), ND3 protein (at position 85), and ND5 protein (at position 400) have one site, respectively, and two sites in ND4 protein (at positions 47 and 200). Based on the results, we hypothesized that the positive selection of Hyla and Dryophytes was due to their experience of cold stress in historical events, but more evidence is needed to support this conclusion.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

Simple Summary

Duplicated control regions have been reported several times in the tree frog family Rhacophoridae, and previous studies have mostly relied on sequence analysis to reconstruct their evolution. This is the first study to employ a phylogenetic method to demonstrate the existence of concerted and parallel evolution succinctly and intuitively in the duplicated control regions of the family Rhacophoridae. Phylogenetic relationships were also used to illustrate the parallel evolution of ATP8 loss of function in the genus Polypedates. In general, this study elucidated the evolutionary patterns and pathways of mitochondrial gene rearrangement of the family Rhacophoridae from a phylogenetic perspective, which aids in understanding the evolutionary history of this fascinating tree frog taxon from a molecular evolution standpoint.

Abstract

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6-trnE-Cytb-D-loop-trnL2 (CUN)-ND5-D-loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.

The complete mitochondrial genome of Pyxicephalus adspersus: high gene rearrangement and phylogenetics of one of the world's largest frogs

The family Pyxicephalidae including two subfamilies (Cacosterninae and Pyxicephalinae) is an ecologically important group of frogs distributed in sub-Saharan Africa. However, its phylogenetic position among the Anura has remained uncertain. The present study determined the complete mitochondrial genome sequence of Pyxicephalus adspersus, the first representative mitochondrial genome from the Pyxicephalinae, and reconstructed the phylogenetic relationships within Ranoidae using 10 mitochondrial protein-coding genes of 59 frog species. The P. adspersus mitochondrial genome showed major gene rearrangement and an exceptionally long length that is not shared with other Ranoidae species. The genome is 24,317 bp in length, and contains 15 protein-coding genes (including extra COX3 and Cyt b genes), four rRNA genes (including extra 12S rRNA and 16S rRNA genes), 29 tRNA genes (including extra tRNALeu (UAG), tRNALeu (UUR), tRNAThr , tRNAPro , tRNAPhe , tRNAVal , tRNAGln genes) and two control regions (CRs). The Dimer-Mitogenome and Tandem duplication and random loss models were used to explain these gene arrangements. Finally, both Bayesian inference and maximum likelihood analyses supported the conclusion that Pyxicephalidae was monophyletic and that Pyxicephalidae was the sister clade of (Petropedetidae + Ptychadenidae).

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).

Complete mitochondrial genomes of Nanorana taihangnica and N. yunnanensis (Anura: Dicroglossidae) with novel gene arrangements and phylogenetic relationship of Dicroglossidae

BACKGROUND

Complete mitochondrial (mt) genomes have been used extensively to test hypotheses about microevolution and to study population structure, phylogeography, and phylogenetic relationships of Anura at various taxonomic levels. Large-scale mt genomic reorganizations have been observed among many fork-tongued frogs (family Dicroglossidae). The relationships among Dicroglossidae and validation of the genus Feirana are still problematic. Hence, we sequenced the complete mt genomes of Nanorana taihangnica (=F. taihangnica) and N. yunnanensis as well as partial mt genomes of six Quasipaa species (dicroglossid taxa), two Odorrana and two Amolops species (Ranidae), and one Rhacophorus species (Rhacophoridae) in order to identify unknown mt gene rearrangements, to investigate the validity of the genus Feirana, and to test the phylogenetic relationship of Dicroglossidae.

RESULTS

In the mt genome of N. taihangnica two trnM genes, two trnP genes and two control regions were found. In addition, the trnA, trnN, trnC, and trnQ genes were translocated from their typical positions. In the mt genome of N. yunnanensis, three control regions were found and eight genes (ND6, trnP, trnQ, trnA, trnN, trnC, trnY and trnS genes) in the L-stand were translocated from their typical position and grouped together. We also found intraspecific rearrangement of the mitochondrial genomes in N. taihangnica and Quasipaa boulengeri. In phylogenetic trees, the genus Feirana nested deeply within the clade of genus Nanorana, indicating that the genus Feirana may be a synonym to Nanorana. Ranidae as a sister clade to Dicroglossidae and the clade of (Ranidae + Dicroglossidae) as a sister clade to (Mantellidae + Rhacophoridae) were well supported in BI analysis but low bootstrap in ML analysis.

CONCLUSIONS

We found that the gene arrangements of N. taihangnica and N. yunnanensis differed from other published dicroglossid mt genomes. The gene arrangements in N. taihangnica and N. yunnanensis could be explained by the Tandem Duplication and Random Loss (TDRL) and the Dimer-Mitogenome and Non-Random Loss (DMNR) models, respectively. The invalidation of the genus Feirana is supported in this study.

The complete mitochondrial genome sequence of the Sichuan Digging Frog, Kaloula rugifera (Anura: Microhylidae) and its phylogenetic implications

The Sichuan Digging Frog (Kaloula rugifera) belongs to the family Dicroglossidae, which is endemic to northeastern Sichuan and southernmost Gansu provinces, in southwestern China. In this study, the complete mitochondrial genome of K. rugifera was sequenced. The mitogenome was 17,074bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding control region. As in other vertebrates, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes which are encoded on the light strand. The overall base composition of the K. rugifera is 30.32% A, 25.76% C, 29.72% T, and 14.20% G, which is consistent with the lowest frequency for G content in typical amphibian animals' mitochondrial genomes. The alignment of the Kaloula species control regions exhibited high genetic variability and rich A+T content. Besides, 3 types of tandem repeat units were also identified in the control region. Phylogenetic tree demonstrated that K. rugifera was clustered together with K. borealis and K. verrucosa and they had a close relationship with each other. The complete mitogenome of K. rugifera can provide an important data for the studies on phylogenetic relationship to further explore the taxonomic status of Kaloula species.

Next-generation sequencing of mixed genomic DNA allows efficient assembly of rearranged mitochondrial genomes in Amolops chunganensis and Quasipaa boulengeri

Recent improvements in next-generation sequencing (NGS) technologies can facilitate the obtainment of mitochondrial genomes. However, it is not clear whether NGS could be effectively used to reconstruct the mitogenome with high gene rearrangement. These high rearrangements would cause amplification failure, and/or assembly and alignment errors. Here, we choose two frogs with rearranged gene order, Amolops chunganensis and Quasipaa boulengeri, to test whether gene rearrangements affect the mitogenome assembly and alignment by using NGS. The mitogenomes with gene rearrangements are sequenced through Illumina MiSeq genomic sequencing and assembled effectively by Trinity v2.1.0 and SOAPdenovo2. Gene order and contents in the mitogenome of A. chunganensis and Q. boulengeri are typical neobatrachian pattern except for rearrangements at the position of "WANCY" tRNA genes cluster. Further, the mitogenome of Q. boulengeri is characterized with a tandem duplication of trnM. Moreover, we utilize 13 protein-coding genes of A. chunganensis, Q. boulengeri and other neobatrachians to reconstruct the phylogenetic tree for evaluating mitochondrial sequence authenticity of A. chunganensis and Q. boulengeri. In this work, we provide nearly complete mitochondrial genomes of A. chunganensis and Q. boulengeri.

Concerted evolution in the mitochondrial control region of the Amazon small-bodied frog Pseudopaludicola canga (Anura, Leiuperidae)

This study presents evidence of concerted evolution in the mitochondrial control region of the frog Pseudopaludicola canga. Four repeat units of 88 bp (as well as a fifth, incomplete unit) were observed in the 5' domain, with the duplicated segments of the same specimen being more related to one another than to the equivalent regions in other specimens, as a result of concerted evolution. We highlight that drawing conclusions from phylogeographical analysis using the control region containing VNTRs must be interpreted with caution, because it violated a basic assumption of phylogeny, since the regions cannot be treated as independent characters.

Mitochondrial phylogeny of an Asian tree frog genus Theloderma (Anura: Rhacophoridae)

We assessed phylogenetic and systematic relationships among 17 out of 23 species of Theloderma and all three species of Nyctixalus from 2412bp sequences of the mitochondrial DNA genes of 12S rRNA, tRNA(val), and 16S rRNA using maximum likelihood and Bayesian inference methods. With the exception of T. moloch, Theloderma and Nyctixalus are confirmed to form a clade, in which each genus also forms a clade. Theloderma moloch is phylogenetically outside these clades and closer to samples from Chiromantis, Feihyla, Gracixalus, Kurixalus, Philautus, Polypedates, Raorchestes, and Rhacophorus. Within Theloderma, T. horridum and T. stellatum form the sister taxon to a clade comprising the remaining species. The basal split within the latter clade groups T. asperum, T. licin, T. petilum, and T. ryabovi as the sister to a clade comprising T. bicolor, T. chuyangsinense, T. corticale, T. gordoni, T. laeve, T. lateriticum, T. nebulosum, T. rhododiscus, and T. truongsonense. Our phylogenetic results indicate homoplastic evolution of four morphological characters: small vs. large body size, presence of vomerine teeth, presence of a vocal opening in males, and interdigital webbing on hands. The common ancestor of Theloderma and Nyctixalus is inferred to have arisen in the area including the current Sunda region.

Mitochondrial genomes of Japanese Babina frogs (Ranidae, Anura): unique gene arrangements and the phylogenetic position of genus Babina

Genus Babina is a member of Ranidae, a large family of frogs, currently comprising 10 species. Three of them are listed as endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered species, we determined the complete nucleotide sequences of the mt genomes of 3 endangered Japanese Babina frogs, B. holsti, B. okinavana, and B. subaspera and 1 ranid frog Lithobates catesbeianus. The genes of NADH dehydrogenase subunit 5 (nad5) and the control region (CR) were found to have high sequence divergences and to be usable for population genetics studies. At present, no consensus on the phylogenetic position of genus Babina has been reached. To resolve this problem, we performed molecular phylogenetic analyses with the largest dataset used to date (11,345 bp from 2 ribosomal RNA- and 13 protein-encoding genes) in studies dealing with Babina phylogeny. These analyses revealed monophyly of Babina and Odorrana. It is well known that mt gene rearrangements of animals can provide usable phylogenetic information. Thus, we also compared the mt gene arrangements among Babina species and other related genera. Of the surveyed species, only L. catesbeianus manifested typical neobatrachian-type mt gene organization. In the B. okinavana, an additional pseudogene of tRNA-His (trnH) was observed in the CR downstream region. Furthermore, in the B. holsti and B. subaspera, the trnH/nad5 block was translocated from its typical position to the CR downstream region, and the translocated trnH became a pseudogene. The position of the trnH pseudogene is consistent with the translocated trnH position reported in Odorrana. Consequently, the trnH rearrangement seems to be a common ancestry characteristic (synapomorphy) of Babina and Odorrana. Based on the "duplication and deletion" gene rearrangement model, a single genomic duplication event can explain the order of derived mt genes found in Babina and Odorrana.

Next generation sequencing and comparative analyses of Xenopus mitogenomes

Background

Mitochondrial genomes comprise a small but critical component of the total DNA in eukaryotic organisms. They encode several key proteins for the cell’s major energy producing apparatus, the mitochondrial respiratory chain. Additonally, their nucleotide and amino acid sequences are of great utility as markers for systematics, molecular ecology and forensics. Their characterization through nucleotide sequencing is a fundamental starting point in mitogenomics. Methods to amplify complete mitochondrial genomes rapidly and efficiently from microgram quantities of tissue of single individuals are, however, not always available. Here we validate two approaches, which combine long-PCR with Roche 454 pyrosequencing technology, to obtain two complete mitochondrial genomes from individual amphibian species.

Results

We obtained two new xenopus frogs (Xenopus borealis and X. victorianus) complete mitochondrial genome sequences by means of long-PCR followed by 454 of individual genomes (approach 1) or of multiple pooled genomes (approach 2), the mean depth of coverage per nucleotide was 9823 and 186, respectively. We also characterised and compared the new mitogenomes against their sister taxa; X. laevis and Silurana tropicalis, two of the most intensely studied amphibians. Our results demonstrate how our approaches can be used to obtain complete amphibian mitogenomes with depths of coverage that far surpass traditional primer-walking strategies, at either the same cost or less. Our results also demonstrate: that the size, gene content and order are the same among xenopus mitogenomes and that S. tropicalis form a separate clade to the other xenopus, among which X. laevis and X. victorianus were most closely related. Nucleotide and amino acid diversity was found to vary across the xenopus mitogenomes, with the greatest diversity observed in the Complex 1 gene nad4l and the least diversity observed in Complex 4 genes (cox1-3). All protein-coding genes were shown to be under strong negative (purifying selection), with genes under the strongest pressure (Complex 4) also being the most highly expressed, highlighting their potentially crucial functions in the mitochondrial respiratory chain.

Conclusions

Next generation sequencing of long-PCR amplicons using single taxon or multi-taxon approaches enabled two new species of Xenopus mtDNA to be fully characterized. We anticipate our complete mitochondrial genome amplification methods to be applicable to other amphibians, helpful for identifying the most appropriate markers for differentiating species, populations and resolving phylogenies, a pressing need since amphibians are undergoing drastic global decline. Our mtDNAs also provide templates for conserved primer design and the assembly of RNA and DNA reads following high throughput “omic” techniques such as RNA- and ChIP-seq. These could help us better understand how processes such mitochondrial replication and gene expression influence xenopus growth and development, as well as how they evolved and are regulated.

Mitogenomic analysis of the genus Panthera

The complete sequences of the mitochondrial DNA genomes of Panthera tigris, Panthera pardus, and Panthera uncia were determined using the polymerase chain reaction method. The lengths of the complete mitochondrial DNA sequences of the three species were 16990, 16964, and 16773 bp, respectively. Each of the three mitochondrial DNA genomes included 13 protein-coding genes, 22 tRNA, two rRNA, one O(L)R, and one control region. The structures of the genomes were highly similar to those of Felis catus, Acinonyx jubatus, and Neofelis nebulosa. The phylogenies of the genus Panthera were inferred from two combined mitochondrial sequence data sets and the complete mitochondrial genome sequences, by MP (maximum parsimony), ML (maximum likelihood), and Bayesian analysis. The results showed that Panthera was composed of Panthera leo, P. uncia, P. pardus, Panthera onca, P. tigris, and N. nebulosa, which was included as the most basal member. The phylogeny within Panthera genus was N. nebulosa (P. tigris (P. onca (P. pardus, (P. leo, P. uncia)))). The divergence times for Panthera genus were estimated based on the ML branch lengths and four well-established calibration points. The results showed that at about 11.3 MYA, the Panthera genus separated from other felid species and then evolved into the several species of the genus. In detail, N. nebulosa was estimated to be founded about 8.66 MYA, P. tigris about 6.55 MYA, P. uncia about 4.63 MYA, and P. pardus about 4.35 MYA. All these estimated times were older than those estimated from the fossil records. The divergence event, evolutionary process, speciation, and distribution pattern of P. uncia, a species endemic to the central Asia with core habitats on the Qinghai-Tibetan Plateau and surrounding highlands, mostly correlated with the geological tectonic events and intensive climate shifts that happened at 8, 3.6, 2.5, and 1.7 MYA on the plateau during the late Cenozoic period.

Mitochondrial genome organization and divergence in hybridizing central European waterfrogs of the Pelophylax esculentus complex (Anura, Ranidae)

Natural transfer of mitochondrial DNA has occurred between three western Palaearctic waterfrog taxa: Pelophylax lessonae, Pelophylax ridibundus and their hybridogenetic hybrid, Pelophylax kl. esculentus. The transfer is asymmetric with most P. kl. esculentus and approximately one third of all central European P. ridibundus having mtDNA derived from P. lessonae (L-mtDNA). We obtained complete nucleotide sequences of multiple mitochondrial genomes (15,376-78 bp without control regions) from all 3 taxa, including a P. ridibundus frog with introgressed L-mtDNA. The gene content and organization of the mitogenomes correspond to those typical of neobatrachians. Divergence between the mtDNAs of P. lessonae and P. ridibundus is high with an uncorrected p-distance of 11.9% across the entire mitogenome. However, the rate of nucleotide substitution depends on the degree of functional constraint with up to 30-fold differences in levels of divergence. In general, mitochondrial genes encoding the translational machinery evolve very slowly, whereas genes encoding polypeptides of the electron transport system, especially the ND genes, evolve rapidly. Only 25 of 211-213 observed amino acid replacements could be classified as radical and are therefore more likely to be exposed to selection. A disproportionately high number of amino acid substitutions has occurred in the ND4, ND4L and cytb genes of the P. lessonae lineage (including 36% of all radical changes). In contrast to the interspecific divergence, nucleotide polymorphism within L- and R-mtDNA is very low: L-mtDNA haplotypes differed on average by only 19 nucleotides, while there was no variation within two mtDNAs derived from P. ridibundus. This is an expected finding considering that we have sampled a post-glacial expansion area. Moreover, the introgressed L-mtDNA on a P. ridibundus background differed from other L-mtDNAs by only a few substitutions, indicative of a very recent introgression event. We discuss our findings in the context of natural selection acting on L-mtDNA and its potential significance in cytonuclear epistasis.

From Antarctica or Asia? New colonization scenario for Australian-New Guinean narrow mouth toads suggested from the findings on a mysterious genus Gastrophrynoides

Background

Microhylidae is a geographically widespread family of anurans. Although several extensive molecular analyses have attempted to elucidate their subfamilial relationships, and correlate these with Mesozoic and Cenozoic continental drifts, consensus has not been reached. Further, generic level relationships have not been well investigated in some microhylid subfamilies, and therefore subfamilial affiliations of some genera are still unclear. To elucidate the phylogenetic positions of two mysterious Asian genera, Gastrophrynoides and Phrynella, and to better understand the trans-continental distributions of microhylid taxa, we performed molecular phylogenetic and dating analyses using the largest molecular dataset applied to these taxa to date.

Results

Six nuclear and two mitochondrial genes (approx. 8 kbp) were sequenced from 22 microhylid frog species representing eight subfamilies. The maximum likelihood and Bayesian analyses could not fully elucidate the subfamilial relationships, suggesting a rapid radiation of these taxa between 85 and 66 million years ago. In contrast, generic relationships of Asian microhylines were generally well resolved.

Conclusion

Our results clearly showed that one of two problematic Asian genera, Phrynella, was nested in the clade of the Asian subfamily Microhylinae. By contrast, Gastrophrynoides occupied the most basal position of the Australian-New Guinean subfamily Asterophryinae. The estimated divergence of Gastrophrynoides from other asterophryine was unexpectedly around 48 million years ago. Although a colonization scenario via Antarctica to the Australian-New Guinean landmass has been suggested for Asterophryinae, our finding suggested a novel colonization route via Indo-Eurasia.

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5' and 3'-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.

Genetic divergence and evolutionary relationships in six species of genera Hoplobatrachus and Euphlyctis (Amphibia: Anura) from Bangladesh and other Asian countries revealed by mitochondrial gene sequences

To elucidate the species composition, genetic divergence, evolutionary relationships, and divergence time of Hoplobatrachus and Euphlyctis frogs (subfamily Dicroglossinae, family Ranidae) in Bangladesh and other Asian countries, we analyzed the mitochondrial Cyt b, 12S, and 16S rRNA genes of 252 specimens. Our phylogenetic analyses showed 13 major clades corresponding to several cryptic species as well as to nominal species in the two genera. The results suggested monophyly of Asian Hoplobatrachus species, but the position of African Hoplobatrachus occipitalis was not clarified. Nucleotide divergence and phylogenetic data suggested the presence of allopatric cryptic species allied to Euphlyctis hexadactylus in Sundarban, Bangladesh and several parapatric cryptic species in the Western Ghats, India. The presence of at least two allopatric cryptic species among diverged Euphlyctis cyanophlyctis in Bangladesh, India, and Sri Lanka was also suggested. In some cases, our estimated divergence times matched the paleogeological events of South and Southeast Asian regions that may have led to the divergence of Hoplobatrachus and Euphlyctis taxa. Especially, land formation at Bangladesh (15-10Ma) may have allowed the spread of these frog taxa to Southeast Asian areas, and the aridification of central India (5.1-1.6Ma) might have affected the gene flow of widely distributed species. The present study revealed prior underestimation of the richness of the amphibian fauna in this region, indicating the possible occurrence of many cryptic species among these groups.

Complete mitochondrial genomes of three neobatrachian anurans: a case study of divergence time estimation using different data and calibration settings

We sequenced the whole mitochondrial (mt) genomes of three neobatrachian species: Japanese tree frog Hyla japonica, Japanese common toad Bufo japonicus, and narrow-mouthed toad Microhyla okinavensis. The gene arrangements of these genomes diverged from that of basal anurans (suborder Archaeobatrachia), but are the same as that of the members of derived frogs (i.e., superfamily Hyloidae and Ranoidae in suborder Neobatrachia), suggesting the one-time occurrence of a gene rearrangement event in an ancestral lineage of derived anurans. Furthermore, several distinct repeat motifs including putative termination-associated sequences (TASs) and conserved sequence blocks (CSBs) were observed in the control regions (CRs) of B. japonicus and H. japonica, while no repeat motifs were found in that of M. okinavensis. Phylogenetic analyses using both nucleotide and amino acid data of mt genes support monophyly of neobatrachians. The estimated divergence time based on amino acid data with multiple reference points suggests that the three living amphibian orders may have originated in the Carboniferous period, and that the divergences of anurans had occurred between the Permian and Tertiary periods. We also checked the influence of the data types and the settings of reference times on divergence time estimation. The resultant divergence times estimated from several datasets and reference time settings suggest that the substitution saturation of nucleotide data may lead to overestimated (i.e., older) branching times, especially for early divergent taxa. We also found a highly accelerated substitution rate in neobatrachian mt genes, and fast substitution possibly resulted in overestimation. To correct this erroneous estimation, it is efficient to apply several reference points among neobatrachians.

Rearrangement of a mitochondrial tRNA gene of the concave-eared torrent frog, Amolops tormotus

In this study, the complete nucleotide sequence (17,962 bp) of the mitochondrial DNA of Amolops tormotus was determined using polymerase chain reaction (PCR). The gene content, base composition and codon usage of A. tormotus conformed to those of typical vertebrate patterns. Among 22 tRNAs, the novel position of the tRNA-His gene was in the D-loop region, which was a novel mtDNA gene rearrangement in amphibians. Phylogenetic analyses were based on a 885-bp sequence of 12S and 16S rRNA for species of Amolops and other related species, concatenated sequences of the 11 protein-encoding genes of 13 species.

Complete nucleotide sequence of the mitochondrial genome of Schlegel's tree frog Rhacophorus schlegelii (family Rhacophoridae): duplicated control regions and gene rearrangements

The complete nucleotide sequence (21,359 bp) of the mitochondrial DNA of the rhacophorid frog Rhacophorus schlegelii was determined. The gene content, nucleotide composition, and codon usage of this genome corresponded to those typical of vertebrates. However, the Rh. schlegelii genome was unusually large due to the inclusion of two control regions and the accumulation of lengthy repetitive sequences in these regions. The two control regions had 97% sequence similarity over 1,510 bp, suggesting the occurrence of concerted sequence evolution. Comparison of the gene organizations among anuran species revealed that the mitochondrial gene arrangement of Rh. schlegelii diverged from that of typical vertebrates but was similar to that of Buergeria buergeri. The positions of the tRNA-Leu(CUN) and tRNA-Thr genes were exchanged between Rh. schlegelii and B. buergeri. Based on parsimonious consideration and the basal phylogenetic position of B. buergeri, these genes seemed to have been rearranged in an ancestral lineage leading to Rh. schlegelii.

Complete nucleotide sequence of the mitochondrial genome of a Malagasy poison frog Mantella madagascariensis: Evolutionary implications on mitochondrial genomes of higher anuran groups

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of a Malagasy poison frog, Mantella madagascariensis (family Mantellidae), and partial sequences of two Mantella (M. baroni and M. bernhardi) and two additional mantellid species (Boophis madagascariensis and Mantidactylus cf. ulcerosus). The M. madagascariensis genome was shown to be the largest (23kbp) of all vertebrate mtDNAs investigated so far. Furthermore, the following unique features were revealed: (1) the positions of some genes and gene regions were rearranged compared to mitochondrial genomes typical for vertebrates and other anuran groups, (2) two distinct genes and a pseudogene corresponding to transfer RNA gene for methionine (tRNA-Met) were encoded, and (3) two control regions with very high sequence homology were present. These features were shared by the two other Mantella species but not the other mantellid species, indicating dynamic genome reorganization in a common ancestor linage before divergence of the Mantella genus. The reorganization pathway could be explained by a model of gene duplication and deletion. Duplication and deletion events also seem to have been responsible for concerted sequence evolution of the control regions in Mantella mt genomes. It is also suggested that the pseudo tRNA-Met gene sustained for a long time in Mantella mt genomes possibly functions as a punctuation marker for NADH dehydrogenase subunit (ND) 2 mRNA processing. Phylogenetic analyses employing a large sequence data set of mt genes supported the monophyly of Mantellidae and Rhacophoridae and other recent phylogenetic views for ranoid frogs. The resultant phylogenetic relationship also suggested parallel occurrence of two tRNA-Met genes, duplicated control regions, and ND5 gene translocation in independent ranoid lineages.