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

Characterization and Comparison of the Two Mitochondrial Genomes in the Genus Rana

The mitochondrial genome (mitogenome) possesses several invaluable attributes, including limited recombination, maternal inheritance, a fast evolutionary rate, compact size, and relatively conserved gene arrangement, all of which make it particularly useful for applications in phylogenetic reconstruction, population genetics, and evolutionary research. In this study, we aimed to determine the complete mitogenomes of two morphologically similar Rana species (Rana hanluica and Rana longicrus) using next-generation sequencing. The entire circular mitogenome was successfully identified, with a length of 19,395 bp for R. hanluica and 17,833 bp for R. longicrus. The mitogenomes of both species contained 37 genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and one control region; mitogenome size varied predominantly with the length of the control region. The two synonymous codon usages in 13 PCGs showed that T and A were used more frequently than G and C. The ratios of non-synonymous to synonymous substitutions of all 13 PCGs were <1 in the Rana species, indicating that the PCGs were under purifying selection. Finally, phylogenetic relationship analyses suggested that R. hanluica and R. longicrus were classified in the R. japonica group. Our study provides valuable reference material for the taxonomy of the genus Rana.

The complete mitochondrial genome of Rana coreana (Anura: Ranidae)

Rana coreana is a brown frog species native to the Korean Peninsula. We characterized the complete mitochondrial genome of the species. The mitochondrial genome sequence of R. coreana is 22,262 bp and comprises 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and two control regions (CRs). The CR duplication and gene organization were identical to those observed in Rana kunyuensis and Rana amurensis. A total of 13 protein-coding genes were used to examine the phylogenetic relationships between this species and the genus Rana. R. coreana living on the Korean Peninsula, formed a cluster with R. kunyuensis and R. amurensis, with R. coreana showing the closest phylogenetic affinity for R. kunyuensis.

Comparative Mitogenomics of True Frogs (Ranidae, Anura), and Its Implications for the Phylogeny and Evolutionary History of Rana

Simple Summary

The true frogs of the genus Rana are a complex and diverse group. Many new species have been discovered with the help of molecular markers and morphological traits. However, the evolutionary history in Rana were not well understood. In this study, we sequenced and annotated the complete mitochondrial genome of R. longicrus and R. zhenhaiensis. In 13 protein codon genes, the COI was the most conserved, and ATP8 had a fast rate of evolution. The Ka/Ks ratio analysis among Rana indicated the protein-coding genes were suffering purify selection. There were three kinds of gene arrangement patterns found. This study provides mitochondrial genetic information, improving our understanding of mitogenomic structure and evolution, and recognizes the phylogenetic relationship and taxonomy among Rana.

Abstract

The true frogs of the genus Rana are a complex and diverse group, containing approximately 60 species with wide distribution across Eurasia and the Americas. Recently, many new species have been discovered with the help of molecular markers and morphological traits. However, the evolutionary history in Rana was not well understood and might be limited by the absence of mitogenome information. In this study, we sequenced and annotated the complete mitochondrial genome of R. longicrus and R. zhenhaiensis, containing 22 tRNAs, 13 protein-coding genes, two ribosomal RNAs, and a non-coding region, with 17,502 bp and 18,006 bp in length, respectively. In 13 protein codon genes, the COI was the most conserved, and ATP8 had a fast rate of evolution. The Ka/Ks ratio analysis among Rana indicated the protein-coding genes were suffering purify selection. There were three kinds of gene arrangement patterns found. The mitochondrial gene arrangement was not related to species diversification, and several independent shifts happened in evolutionary history. Climate fluctuation and environmental change may have played an essential role in species diversification in Rana. This study provides mitochondrial genetic information, improving our understanding of mitogenomic structure and evolution, and recognizes the phylogenetic relationship and taxonomy among Rana.

The complete mitochondrial genome of Rana johnsi (Smith, 2009) (Anura: Ranidae) and its phylogeny

Rana johnsi (Smith 2009) firstly considered as the member of genus Pseudorana, has been moved into the genus Rana. In this study, we sequenced the complete mitochondrial (mt) genome of R. johnsi using the Sanger method. The circular mt genome was 17,873 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosome RNA genes, and one control region. The overall nucleotide composition in majority-strand was 28% A, 29% T, 29% C, and 14% G. We discussed the phylogenetic relationship of R. johnsi in genus Rana using ML and BI analyses based on 13 PCGs. Excluding the clade of subgenus Lithobates, Rana draytonii was the basal clade to all other Rana species, which included R. johnsi as the basal clade. The monophyly of genus Rana was supported, whereas Pseudorana was failed to support.