The complete mitochondrial genome of Dryophytes versicolor: Phylogenetic relationship among Hylidae and mitochondrial protein-coding gene expression in response to freezing and anoxia

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

Differential Mitochondrial Genome Expression of Three Sympatric Lizards in Response to Low-Temperature Stress

Ecological factors related to climate extremes have a significant influence on the adaptability of organisms, especially for ectotherms such as reptiles that are sensitive to temperature change. Climate extremes can seriously affect the survival and internal physiology of lizards, sometimes even resulting in the loss of local populations or even complete extinction. Indeed, studies have shown that the expression levels of the nuclear genes and mitochondrial genomes of reptiles change under low-temperature stress. At present, the temperature adaptability of reptiles has rarely been studied at the mitochondrial genome level. In the present study, the mitochondrial genomes of three species of lizards, Calotes versicolor, Ateuchosaurus chinensis, and Hemidactylus bowringii, which live in regions of sympatry, were sequenced. We used RT-qPCR to explore the level of mitochondrial gene expression under low-temperature stress, as compared to a control temperature. Among the 13 protein-coding genes (PCGs), the steady-state transcript levels of ND4L, ND1, ATP6, and COII were reduced to levels of 0.61 ± 0.06, 0.50 ± 0.08, 0.44 ± 0.16, and 0.41 ± 0.09 in C. versicolor, respectively, compared with controls. The transcript levels of the ND3 and ND6 genes fell to levels of just 0.72 ± 0.05 and 0.67 ± 0.05 in H. bowringii, compared with controls. However, the transcript levels of ND3, ND5, ND6, ATP6, ATP8, Cytb, and COIII genes increased to 1.97 ± 0.15, 2.94 ± 0.43, 1.66 ± 0.07, 1.59 ± 0.17, 1.46 ± 0.04, 1.70 ± 0.16, and 1.83 ± 0.07 in A. chinensis. Therefore, the differences in mitochondrial gene expression may be internally related to the adaptative strategy of the three species under low-temperature stress, indicating that low-temperature environments have a greater impact on A. chinensis, with a small distribution area. In extreme environments, the regulatory trend of mitochondrial gene expression in reptiles is associated with their ability to adapt to extreme climates, which means differential mitochondrial genome expression can be used to explore the response of different lizards in the same region to low temperatures. Our experiment aims to provide one new research method to evaluate the potential extinction of reptile species in warm winter climates.

Development of molecular identification methods for Dryophytes suweonensis and D. japonicus, and their hybrids

Background

As hybridization can reduce biodiversity or cause extinction, it is important to identify both purebred parental species and their hybrids prior to conserving them. The Suwon tree frog, Dryophytes suweonensis, is an endangered wildlife species in Korea that shares its habitat and often hybridizes with the Japanese tree frog, D. japonicus. In particular, D. suweonensis, D. japonicus, and their hybrids often have abnormal ovaries and gonads, which are known causes that could threaten their existence.

Methods

We collected 57 individuals from six localities where D. suweonensis is known to be present. High-resolution melting curve (HRM) analysis of the mitochondrial 12S ribosomal RNA gene was performed to determine the maternal species. Thereafter, the DNA sequences of five nuclear genes (SIAH, TYR, POMC, RAG1, and C-MYC) were analyzed to determine their parental species and hybrid status.

Results

The HRM analysis showed that the melting temperature of D. suweonensis was in the range of 79.0-79.3 °C, and that of D. japonicus was 77.7-78.0 °C, which clearly distinguished the two tree frog species. DNA sequencing of the five nuclear genes revealed 37 single-nucleotide polymorphism (SNP) sites, and STRUCTURE analysis showed a two-group structure as the most likely grouping solution. No heterozygous position in the purebred parental sequences with Q values ≥ 0.995 were found, which clearly distinguished the two treefrog species from their hybrids; 11 individuals were found to be D. suweonensis, eight were found to be D. japonicus, and the remaining 38 individuals were found to be hybrids.

Conclusion

Thus, it was possible to unambiguously identify the parental species and their hybrids using HRM analysis and DNA sequencing methods. This study provided fundamental information for D. suweonensis conservation and restoration research.

How Does Mitochondrial Protein-Coding Gene Expression in Fejervarya kawamurai (Anura: Dicroglossidae) Respond to Extreme Temperatures?

Unusual climates can lead to extreme temperatures. Fejervarya kawamurai, one of the most prevalent anurans in the paddy fields of tropical and subtropical regions in Asia, is sensitive to climate change. The present study focuses primarily on a single question: how do the 13 mitochondrial protein-coding genes (PCGs) respond to extreme temperature change compared with 25 °C controls? Thirty-eight genes including an extra tRNA-Met gene were identified and sequenced from the mitochondrial genome of F. kawamurai. Evolutionary relationships were assessed within the Dicroglossidae and showed that Dicroglossinae is monophyletic and F. kawamurai is a sister group to the clade of (F. multistriata + F. limnocharis). Transcript levels of mitochondrial genes in liver were also evaluated to assess responses to 24 h exposure to low (2 °C and 4 °C) or high (40 °C) temperatures. Under 2 °C, seven genes showed significant changes in liver transcript levels, among which transcript levels of ATP8, ND1, ND2, ND3, ND4, and Cytb increased, respectively, and ND5 decreased. However, exposure to 4 °C for 24 h was very different in that the expressions of ten mitochondrial protein-coding genes, except ND1, ND3, and Cytb, were significantly downregulated. Among them, the transcript level of ND5 was most significantly downregulated, decreasing by 0.28-fold. Exposure to a hot environment at 40 °C for 24 h resulted in a marked difference in transcript responses with strong upregulation of eight genes, ranging from a 1.52-fold increase in ND4L to a 2.18-fold rise in Cytb transcript levels, although COI and ND5 were reduced to 0.56 and 0.67, respectively, compared with the controls. Overall, these results suggest that at 4 °C, F. kawamurai appears to have entered a hypometabolic state of hibernation, whereas its mitochondrial oxidative phosphorylation was affected at both 2 °C and 40 °C. The majority of mitochondrial PCGs exhibited substantial changes at all three temperatures, indicating that frogs such as F. kawamurai that inhabit tropical or subtropical regions are susceptible to ambient temperature changes and can quickly employ compensating adjustments to proteins involved in the mitochondrial electron transport chain.

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.

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.

Mitochondrial gene expression in different organs of Hoplobatrachus rugulosus from China and Thailand under low-temperature stress

Background

Hoplobatrachus rugulosus (Anura: Dicroglossidae) is distributed in China and Thailand and the former can survive substantially lower temperatures than the latter. The mitochondrial genomes of the two subspecies also differ: Chinese tiger frogs (CT frogs) display two identical ND5 genes whereas Thai tiger frogs (TT frogs) have two different ND5 genes. Metabolism of ectotherms is very sensitive to temperature change and different organs have different demands on energy metabolism at low temperatures. Therefore, we conducted studies to understand: (1) the differences in mitochondrial gene expression of tiger frogs from China (CT frogs) versus Thailand (TT frogs); (2) the differences in mitochondrial gene expression of tiger frogs (CT and TT frogs) under short term 24 h hypothermia exposure at 25 °C and 8 °C; (3) the differences in mitochondrial gene expression in three organs (brain, liver and kidney) of CT and TT frogs.

Results

Utilizing RT-qPCR and comparing control groups at 25 °C with low temperature groups at 8 °C, we came to the following results. (1) At the same temperature, mitochondrial gene expression was significantly different in two subspecies. The transcript levels of two identical ND5 of CT frogs were observed to decrease significantly at low temperatures (P < 0.05) whereas the two different copies of ND5 in TT frogs were not. (2) Under low temperature stress, most of the genes in the brain, liver and kidney were down-regulated (except for COI and ATP6 measured in brain and COI measured in liver of CT frogs). (3) For both CT and TT frogs, the changes in overall pattern of mitochondrial gene expression in different organs under low temperature and normal temperature was brain > liver > kidney.

Conclusions

We mainly drew the following conclusions: (1) The differences in the structure and expression of the ND5 gene between CT and TT frogs could result in the different tolerances to low temperature stress. (2) At low temperatures, the transcript levels of most of mitochondrial protein-encoding genes were down-regulated, which could have a significant effect in reducing metabolic rate and supporting long term survival at low temperatures. (3) The expression pattern of mitochondrial genes in different organs was related to mitochondrial activity and mtDNA replication in different organs.

Cryptic Species Exist in Vietnamella sinensis Hsu, 1936 (Insecta: Ephemeroptera) from Studies of Complete Mitochondrial Genomes

Ephemeroptera (Insecta: Pterygota) are widely distributed all over the world with more than 3500 species. During the last decade, the phylogenetic relationships within Ephemeroptera have been a hot topic of research, especially regarding the phylogenetic relationships among Vietnamellidae. In this study, three mitochondrial genomes from three populations of Vienamella sinensis collected from Tonglu (V. sinensis TL), Chun’an (V. sinensis CN), and Qingyuan (V. sinensis QY) in Zhejiang Province, China were compared to discuss the potential existence of cryptic species. We also established their phylogenetic relationship by combining the mt genomes of 69 Ephemeroptera downloaded from NCBI. The mt genomes of V. sinensis TL, V. sinensis CN, and V. sinensis QY showed the same gene arrangement with lengths of 15,674 bp, 15,674 bp, and 15,610 bp, respectively. Comprehensive analyses of these three mt genomes revealed significant differences in mt genome organization, genetic distance, and divergence time. Our results showed that the specimens collected from Chun’an and Tonglu in Zhejiang Province, China belonged to V. sinensis, and the specimens collected from Qingyuan, Zhejiang Province, China were a cryptic species of V. sinensis. In maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees, the monophyly of the family Vietnamellidae was supported and Vietnamellidae has a close relationship with Ephemerellidae.

The complete mitochondrial genome of Choroterpes (Euthralus) yixingensis (Ephemeroptera: Leptophlebiidae) and its mitochondrial protein-coding gene expression under imidacloprid stress

As one of the most common benthic invertebrates in freshwater, mayflies are very sensitive to changes in water quality and have high requirements for the water environment to allow their nymphs to successfully live and grow. Neonicotinoids, such as imidacloprid, can enter fresh water and pollute the aquatic environment. The present study had two goals: (1) investigate imidacloprid effects on mayfly larvae Choroterpes (Euthralus) yixingensis, and (2) contribute to the phylogenetic status of Ephemeroptera that has always been controversial. Nymphs were collected from Jinhua, China and exposed to different concentrations imidacloprid (5, 10, 20, and 40 μg/L) in the laboratory. Survival of C. yixingensis nymphs decreased as a function of time and imidacloprid concentration with only ~ 55% survival after 72 h exposure to 40 μg/L imidacloprid. After culture under 40 μg/L imidacloprid for 24 h, the steady state transcript levels of mitochondrial COX3, ND4 and ND4L genes were reduced to just 0.07 ± 0.11, 0.30 ± 0.16, and 0.28 ± 0.13 as compared with respective control values (P < 0.01). Steady state transcript levels of ND4 and ND4L were also significantly reduced in a dose-dependent manner (P < 0.05), suggesting that the steady state transcript pattern of these genes in mayfly nymphs can change in response to different levels of environmental contamination. Hence, the mitochondrial protein-coding genes of mayflies could potentially be developed as biomarkers for water ecotoxicity monitoring in the future. In addition, we used the mitochondrial genome sequence of C. yixingensis for an assessment of the phylogenetic tree of Ephemeroptera. The monophyly of Leptophlebiidae was supported and showed that Leptophlebiidae was a sister group to the clade (Baetidae + Caenidae).

Mitogenome Analysis of Four Lamiinae Species (Coleoptera: Cerambycidae) and Gene Expression Responses by Monochamus alternatus When Infected with the Parasitic Nematode, Bursaphelenchus mucronatus

We determined the mitochondrial gene sequence of Monochamus alternatus and three other mitogenomes of Lamiinae (Insect: Coleoptera: Cerambycidae) belonging to three genera (Aulaconotus, Apriona and Paraglenea) to enrich the mitochondrial genome database of Lamiinae and further explore the phylogenetic relationships within the subfamily. Phylogenetic trees of the Lamiinae were built using the Bayesian inference (BI) and maximum likelihood (ML) methods and the monophyly of Monochamus, Anoplophora, and Batocera genera was supported. Anoplophora chinensis, An. glabripennis and Aristobia reticulator were closely related, suggesting they may also be potential vectors for the transmission of the pine wood pathogenic nematode (Bursaphelenchus xylophilus) in addition to M. alternatus, a well-known vector of pine wilt disease. There is a special symbiotic relationship between M. alternatus and Bursaphelenchus xylophilus. As the native sympatric sibling species of B. xylophilus, B. mucronatus also has a specific relationship that is often overlooked. The analysis of mitochondrial gene expression aimed to explore the effect of B. mucronatus on the energy metabolism of the respiratory chain of M. alternatus adults. Using RT-qPCR, we determined and analyzed the expression of eight mitochondrial protein-coding genes (COI, COII, COIII, ND1, ND4, ND5, ATP6, and Cty b) between M. alternatus infected by B. mucronatus and M. alternatus without the nematode. Expression of all the eight mitochondrial genes were up-regulated, particularly the ND4 and ND5 gene, which were up-regulated by 4-5-fold (p < 0.01). Since longicorn beetles have immune responses to nematodes, we believe that their relationship should not be viewed as symbiotic, but classed as parasitic.

Unique Evolution of Hyla hallowellii Among Amphibians of the Central Ryukyus, Japan (Anura: Hylidae)

Fauna of the Central Ryukyus includes a high percentage of endemic species, and Hallowell's tree frog Hyla hallowellii Thompson, 1912 is one of such elements, occurring in a total of eight islands in the Amami and Okinawa Island groups. Using samples representing all of these eight island populations, we studied variations in morphology, karyotype, allozyme, and mtDNA, to clarify the pattern of geographic differentiation of H. hallowellii and consider factors for its formation. We could not clearly discriminate one population from another in morphology, nor could we find any interpopulation difference in karyotype. From genetic analyses, using allozymes and cyt b, we found low overall differentiations among populations. However, the southern populations from Okinawajima and Yoronjima were genetically nearly identical with the northern Amamioshima population. From that group the geographically intermediate Tokunoshima and Kakeromajima populations showed prominent differentiations. These patterns of geographical differentiation greatly differ from those known in other amphibian species of the Central Ryukyus, and suggest that H. hallowellii has evolutionary history unique to the species.

Mitochondria and the Frozen Frog

The wood frog, Rana sylvatica, is the best-studied of a small group of amphibian species that survive whole body freezing during the winter months. These frogs endure the freezing of 65-70% of their total body water in extracellular ice masses. They have implemented multiple adaptations that manage ice formation, deal with freeze-induced ischemia/reperfusion stress, limit cell volume reduction with the production of small molecule cryoprotectants (glucose, urea) and adjust a wide variety of metabolic pathways for prolonged life in a frozen state. All organs, tissues, cells and intracellular organelles are affected by freeze/thaw and its consequences. This article explores mitochondria in the frozen frog with a focus on both the consequences of freezing (e.g., anoxia/ischemia, cell volume reduction) and mitigating defenses (e.g., antioxidants, chaperone proteins, upregulation of mitochondria-encoded genes, enzyme regulation, etc.) in order to identify adaptive strategies that defend and adapt mitochondria in animals that can be frozen for six months or more every year. A particular focus is placed on freeze-responsive genes in wood frogs that are encoded on the mitochondrial genome including ATP6/8, ND4 and 16S RNA. These were strongly up-regulated during whole body freezing (24 h at -2.5 °C) in the liver and brain but showed opposing responses to two component stresses: strong upregulation in response to anoxia but no response to dehydration stress. This indicates that freeze-responsive upregulation of mitochondria-encoded genes is triggered by declining oxygen and likely has an adaptive function in supporting cellular energetics under indeterminate lengths of whole body freezing.

Yellow sea mediated segregation between North East Asian Dryophytes species

While comparatively few amphibian species have been described on the North East Asian mainland in the last decades, several species have been the subject of taxonomical debates in relation to the Yellow sea. Here, we sampled Dryophytes sp. treefrogs from the Republic of Korea, the Democratic People's Republic of Korea and the People's Republic of China to clarify the status of this clade around the Yellow sea and determine the impact of sea level change on treefrogs' phylogenetic relationships. Based on genetics, call properties, adult morphology, tadpole morphology and niche modelling, we determined the segregated status species of D. suweonensis and D. immaculatus. We then proceeded to describe a new treefrog species, D. flaviventris sp. nov., from the central lowlands of the Republic of Korea. The new species is geographically segregated from D. suweonensis by the Chilgap mountain range and known to occur only in the area of Buyeo, Nonsan and Iksan in the Republic of Korea. While the Yellow sea is the principal element to the current isolation of the three clades, the paleorivers of the Yellow sea basin are likely to have been the major factor for the divergences within this clade. We recommend conducting rapid conservation assessments as these species are present on very narrow and declining ranges.

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.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

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