The complete mitochondrial genome ofBabina adenopleura(Anura: Ranidae)

A new odorous frog species of Odorrana (Amphibia, Anura, Ranidae) from Guizhou Province, China

The frog genus Odorrana is distributed across east and southeastern Asia. Based on morphological differences and molecular phylogenetics, a new species of the genus occurring from Leigong Mountain in Guizhou Province, China is described. Phylogenetic analyses based on DNA sequences of the mitochondrial 12S rRNA, 16S rRNA, and ND2 genes supported the new species as an independent lineage. The uncorrected genetic distances between the 12S rRNA, 16S rRNA, and ND2 genes between the new species and its closest congener were 5.0%, 4.9%, and 16.3%, respectively. The new species is distinguished from its congeners by a combination of the following characters: body size moderate (SVL 39.1-49.4 mm in males, 49.7 mm in female); head width larger than head length; tympanum distinctly visible; small rounded granules scattered all over dorsal body and limbs; dorsolateral folds absent; heels overlapping when thighs are positioned at right angles to the body; tibiotarsal articulation reaching the level between eye to nostril when leg stretched forward; vocal sacs absent in male and nuptial pads present on the base of finger I.

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.

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.

Characterization of the mitochondrial genomes of two toads, Anaxyrus americanus (Anura: Bufonidae) and Bufotes pewzowi (Anura: Bufonidae), with phylogenetic and selection pressure analyses

Mitogenomes are useful in analyzing phylogenetic relationships and also appear to influence energy metabolism, thermoregulation and osmoregulation. Much evidence has accumulated for positive selection acting on mitochondrial genes associated with environmental adaptation. Hence, the mitogenome is a likely target for environmental selection. The family Bufonidae (true toads) has only nine complete and four partial mitogenomes published compared to the 610 known species of this family. More mitogenomes are needed in order to obtain a clearer understanding of the phylogenetic relationships within Bufonidae that are currently controversial. To date, no mitogenomes have been reported from the genera Anaxyrus and Bufotes. Anaxyrus americanus can live in low temperature environments and Bufotes pewzowi can live in high salinity environments. We sequenced the mitogenomes of these two species to discuss the phylogenetic relationships within Bufonidae and the selection pressures experienced by specimens living in low temperature or saline environments. Like other toads, the circular mitogenomes of both species contained the typical 37 genes. Anaxyrus americanus had the highest A+T content of the complete mitogenome among the Bufonidae. In addition, A. americanus showed a negative AT-skew in the control region, whereas Bufotes pewzowi showed a positive AT-skew. Additionally, both toad species had unique molecular features in common: an ND1 gene that uses TTG as the start codon, an extra unpaired adenine (A) in the anticodon arm of trnS (AGY), and the loss of the DHU loop in trnC. The monophyly of Bufonidae was corroborated by both BI and ML trees. An analysis of selective pressure based on the 13 protein coding genes was conducted using the EasyCodeML program. In the branch model analysis, we found two branches of A. americanus and Bufotes pewzowi that were under negative selection. Additionally, we found two positively selected sites (at positions 115 and 119, BEB value > 0.90) in the ND6 protein in the site model analysis. The residue D (119) was located only in A. americanus and may be related to adaptive evolution in low temperature environments. However, there was no evidence of a positively selected site in Bufotes pewzowi in this study.

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).

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.

DNA barcoding and molecular phylogeny in Ranidae

The family Ranidae has the widest distribution compared with other frog family. The phylogeny of the Ranidae is still a matter of dispute. In the present study, we analyzed the COI barcodes of 29 species from six genera belonging to the family Ranidae. Twenty-seven species (93.10% of all the species) were correctly identified by their DNA barcodes. Pelophylax lessonae and Pelophylax ridibundus shared the same one barcode sequence. Kimura two-parameter distances were calculated between barcodes. Pair-wise comparisons among-species were distributed from 0.16% (between Pelophylax lessonae and Pelophylax esculenta) to 29.13% (between Rana warszewitschii and Rana dybowskii). The average genetic distance between species was 28 times higher than the average genetic distance within species. The neighbor-joining method was used to construct a phylogenetic tree, which grouped all the genera into two divergent clades. The results indicated that some currently recognized genera of Ranidae may not be monophyletic. COI gene data supported the hypothesis of polyphyly for Rana, Amolops, Babina, and Hylarana. DNA barcoding is an effective molecular tool for Ranidae species identification and phylogenetic inference.

The complete mitochondrial genome of Hyla annectans (Anura: Hylidae)

The complete mitochondrial genome sequence of the Hyla annectans (Anura: Hylidae) was determined in this study. It is a circular molecule of 17,973 bp in length including 37 genes for 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs and a control region, showing a typical neobatrachian-type arrangement. The A + T content of H-strand and the control region is 60.7% and 64.6%, respectively. The control region is 2554 bp in length and possesses distinct repeat regions at both 5' and 3' ends.

Complete mitochondrial genomes of two brown frogs, Rana dybowskii and Rana cf. chensinensis (Anura: Ranidae)

We first determined complete mitochondrial genomes of Rana dybowskii and Rana cf. chensinensis (Anura: Ranidae). The mitogenomic lengths of R. dybowskii and R. cf. chensinensis were 18,864 and 18,808 bp, respectively. The two mitogenomes have similar gene compositions including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a control region. Rana dybowskii and R. cf. chensinensis mitogenomes displayed same gene order arrangements and similar base compositions with an A + T bias. Mitogenomic data of the two species contributed to provide molecular marker for their conservative genetics and clarified their phylogenetic position under mitogenome-based phylogeny of the genus Rana.

The complete mitochondrial genome of Lithobates catesbeianus (Anura: Ranidae)

The complete mitochondrial genome of Lithobates catesbeianus (Anura: Ranidae) is sequenced to analyze the gene arrangement. It is a circular molecule of 18,241 bp in length including 37 genes typically found in other frogs. The AT content of the overall base composition of L. catesbeianus is 59.9%. The length of control region is 2783 bp with 66.0% AT content. Protein-coding genes begin with ATG as start codon except except ND1 and ATP6 began with ATA, COI and ND4L with GTG, and ND2 with ATT. COI end with AGG as stop codon, COII and ND6 end with AGA, ND2 end with TAG, ATP8. ND4L. ND5 and Cytb end with TAA, and the other five PCGs end with a incomplete stop codon (a single stop nucleotide T).