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

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.

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

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.

The mitochondrial genome and phylogenetic analysis of Rhacophorus rhodopus

Classification of the genus Rhacophorus has been problematic. In particular there has been considerable controversy surrounding the phylogenetic relationships among Rhacophorus rhodopus, R. bipunctatus, and R. reinwardtii. To examine the relationship among these Rhacophorus species, we assembled the complete mitochondrial genome sequence of R. rhodopus. The R. rhodopus genome is 15,789 bp in length with 12 protein-coding genes (PCGs) (losing ND5), two ribosomal genes, 22 transfer RNA genes, and a control region (D-loop). Base composition of the overall sequence was 60.86% for A + T content and 39.14% for C + G content. Most of the PCGs used ATG as a start codon, except for the COX I gene, which used the ATA start codon. COX I and ND6 used AGG and ATP8 stop codons respectively, while ND3 and ND4L used the TAA stop codon. For the remaining seven genes, the stop codons was incomplete. In addition, both 5' and 3' of the control areas had distinct repeating regions. Based on three datasets and two methods (Bayesian inference (BI) and maximum likelihood (ML)), we reconstructed three phylogenetic trees to explore the taxonomic status of the species and the phylogenetic relationship among R. rhodopus, R. bipunctatus and R. reinwardtii. Our results indicated that these three species are non-monophyletic; thus, the phylogenetic relationship among them is complex and difficult to determine. Further, R. rhodopus is divided into three lineages from different parts of China. The two Rhacophorus samples showed very close phylogenetic relationship with R. rhodopus. Our results add to the mitochondrial genome database of amphibians and will help to disentangle the phylogenetic relationships within the Rhacophoridae.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

Comparative analysis of the complete mitochondrial genomes of four cordyceps fungi

Cordyceps is a large group of entomogenous, medicinally important fungi. In this study, we sequenced, assembled, and annotated the entire mitochondrial genome of Ophiocordyceps xuefengensis, in addition to comparing it against other three complete cordyceps mitogenomes that were previously published. Comparative analysis indicated that the four complete mitogenomes are all composed of circular DNA molecules, although their sizes significantly differ due to high variability in intron and intergenic region sizes in the Ophiocordyceps sinensis and O. xuefengensis mitogenomes. All mitogenomes contain 14 conserved genes and two ribosomal RNA genes, but varying numbers of tRNA introns. The Ka/Ks ratios for all 14 PCGs and rps3 were all less than 1, indicating that these genes have been subject to purifying selection. Phylogenetic analysis was conducted using concatenated amino acid and nucleotide sequences of the 14 PCGs and rps3 using two different methods (Maximum Likelihood and Bayesian analysis), revealing highly supported relationships between O. xuefengensis and other Ophiocordyceps species, in addition to a close relationship with O. sinensis. Further, the analyses indicated that cox1 and rps3 play important roles in population differentiation. These mitogenomes will allow further study of the population genetics, taxonomy, and evolutionary biology of medicinally important cordyceps species.

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

Comparative mitogenomic and phylogenetic analysis of Apalone spinifera and Apalone ferox (Testudines: Trionychidae)

The soft-shell turtles Apalone spinifera (AS) and Apalone ferox (AF) are two important economic species. AF is found in the Yellow River of China, and is a confirmed member of the Trionychidae family. However, the classification of AS was in dispute. Mitochondrial genomes (mitogenomes) have been widely used for species identification, as well as population and phylogenetic analysis. In order to understand the phylogenetic and mitogenomic features of AS and AF, the complete mitogenomes were sequenced, annotated and analyzed in this study. The complete mitogenomes of AS and AF are 16,817 bp and 16,756 bp in length, respectively. Both mitogenomes contain 37 genes, seven short intergenic spacers and two long intergenic spacers. Comparative analysis showed that there are 1,137 variation sites (6.79%) between the two mitogenomes. AS and AF mitogenomes both show a usage preference in terms of nucleotides, codons and amino acids. In addition, the non-synonymous substitution rate/synonymous substitution rate indicates that all protein-coding genes (PCGs) have undergone a strong purifying selection. Phylogenetic trees constructed by 13 PCGs show a clear phylogenetic relationship of the soft-shell turtles and suggest that AS is a sister species to AF of the genus Apalone. The data could be useful for further research of species identification, population analysis and the mitogenomic features of soft-shell turtles.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.