Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs

Enhanced dynamicity: evolutionary insights into amphibian mitogenomes architecture

Mitogenomes are known for their structural dynamics and the complexity of their rearrangement patterns. However, their utility in metazoan comparative biology has not been fully exploited. Vertebrate mitogenomes are now sufficiently representative to allow the development of more advanced methods for comparing genome architecture. Furthermore, the relatively robust phylogeny of vertebrates at higher taxonomic ranks allows us to infer the patterns of genome evolution accordingly. In this study, using amphibians as an example, we performed data cleaning and manual annotation on 1777 samples from the NCBI and identified 88 rearrangement types, most of which were clade-specific. In addition, we quantified genomic changes in an evolutionary framework and obtained stepwise growth curves of the architectural changes. This study provides new perspectives for understanding the evolution of the mitogenomes in amphibians and is expected to facilitate the qualitative and quantitative development of mitogenomes research.

Three new species of torrent treefrogs (Anura, Hylidae) of the Hyloscirtusbogotensis group from the eastern Andean slopes and the biogeographic history of the genus

The Hyloscirtusbogotensis group contains 17 species of treefrogs from the tropical Andes and Central America. A taxonomic review of the Amazonian clades of this group is presented based on DNA sequences of nuclear and mitochondrial DNA and a preliminary phylogenomic analysis of ultraconserved elements, as well as morphological, bioacoustic, and environmental characters. Additionally, the role of the Andes in the diversification of the genus Hyloscirtus is explored by reconstructing their ancestral basin (Amazon, Pacific, Caribbean). Our integrative analysis indicates the existence of eight undescribed candidate species within the group. Three of those species are described, previously masked within H.albopunctulatus, H.phyllognathus, and H.torrenticola. A lectotype is also designated for Hylaalbopunctulata. The new evidence suggests that neither Hyloscirtusphyllognathus nor H.torrenticola occur in Ecuador. The new species, H.elbakyanae sp. nov., H.dispersus sp. nov., and Hyloscirtusmaycu sp. nov. differ from other members of the group in bioacoustics and external morphology. The most useful diagnostic characters among species were advertisement calls. In contrast, skin coloration is highly variable intraspecifically and, as a result, of low diagnostic value. High variation in color is partly a result of phenotypic plasticity. Our biogeographic reconstructions indicate that the Andean barrier influenced the diversification of Hyloscirtus. Since the early Oligocene, there have been only four colonization events across de Andes, between the Pacific and Amazon basins. Two of those events occurred more than 14 Mya, when most of the tropical Andes were below 3000 m. Species in the highland H.larinopygion group are younger, suggesting recent diversification as high montane forests and paramo habitats emerged.

First Complete Mitochondrial Genome Analysis of Tree Frog, Dryophytes flaviventris and Comparison with Dryophytes suweonensis

Mitochondrial genomes (mitogenomes) play a key role in species identification and phylogenetic studies due to their stable gene arrangements and evolutionary insights. Dryophytes flaviventris, classified in 2020 and closely related to D. suweonensis, lacks mitochondrial DNA data for differentiation. This gap hinders accurate species identification, highlighting the need for further genomic studies. The complete mitogenome size of two D. flaviventris were 18,616-18,617 bp and those for two D. suweonensis were 18,610-18,616 bp, the mitogenomes of the two species consisting of 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA (tRNA) genes, and a D-loop. Phylogenetic analysis confirmed that the mitochondrial DNA of all four individuals formed a monophyletic group, showing no genetic differentiation. As a result, the two species do not form distinct clades, and mitogenomes could not differentiate them, suggesting they are not reciprocally monophyletic. This study presents the first mitogenome data for D. flaviventris and provides valuable insights into tree frog taxonomy.

Digging out another cryptic species: A new sand frog (Anura: Pyxicephalidae: Tomopterna) from Mozambique

A new species of Tomopterna is described from Banhine National Park, southern Mozambique. It differs from all other sand frogs by an uncorrected p-distance of 4.9-8.7% for 16S rRNA. The new species can be distinguished from all other described species of sand frogs by: males 37 mm SVL or less, undivided subarticular tubercles, more than three phalanges free of webbing on the fourth toe, a continuous glandular ridge below the tympanum, a distinct tympanic membrane, the nostrils closer to the snout tip than the eye, the rounded canthus rostralis, a round tympanum, the absence of an outer metatarsal tubercle, an interrupted inter-orbital bar, small dorsal warts, the absence (or weakness) of vomerine teeth, the solid pigmented throat in mature males, the distinct palmar tubercles, the smooth ventral surface, and barring on the limbs. This description adds a further amphibian to the list of Mozambican frogs, which has increased significantly in recent years.

The systematic codon usage bias has an important effect on genetic adaption in native species

Random mutations increase genetic variety and natural selection enhances adaption over generations. Codon usage biases (CUB) provide clues about the genome adaptation mechanisms of native species and extremophile species. Significant numbers of gene (CDS) of nine classes of endangered, native species, including extremophiles and mesophiles were utilised to compute CUB. Codon usage patterns differ among the lineages of endangered and extremophiles with native species. Polymorphic usage of nucleotides with codon burial suggests parallelism of native species within relatively confined taxonomic groups. Utilizing the deviation pattern of CUB of endangered and native species, I present a calculation parameter to estimate the extinction risk of endangered species. Species diversity and extinction risk are both positively associated with the propensity of random mutation in CDS (Coding DNA sequence). Codon bias tenet profoundly selected and it governs to adaptive evolution of native species.

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.

The status of Tomopterna pulchra Boulenger, 1896 (Anura: Pyxicephalidae)

Tomopterna pulchra (Boulenger, 1896) was described from Lake Tanganyika based on a single specimen. It was synonymised with Tomopterna tuberculosa (Boulenger, 1882) by Loveridge (1957) who noted that it did not fully agree with the description of T. tuberculosa, especially in the difference in dorsal colouration. Genetic differences show that two species are confused, although their advertisement calls and morphology are similar, except for the presence of a pale vertebral stripe mostly in eastern populations. We consider the advertisement calls, morphology, and genetic differences, and remove Tomopterna pulchra (Boulenger, 1896) from the synonymy of Tomopterna tuberculosa. Divergence dating shows that these two species separated in the mid-Miocene. Tomopterna pulchra breeds in rocky habitats, whereas T. tuberculosa breeds in sandy areas, typical of the genus. Based on genetically confirmed specimens, presently T. tuberculosa is known in the west (Angola and Namibia) while T. pulchra is only known from the east (Zambia and Tanzania). The examination of specimens for the presence or absence of a vertebral stripe infers that these species are widespread. Further work is required to determine if the two species are sympatric, and the extent of their ranges.

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.

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.

Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90β1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIβ1 sequences. The G90β1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

Crossing the Rift valley: using complete mitogenomes to infer the diversification and biogeographic history of ethiopian highlands Ptychadena (anura: Ptychadenidae)

The Ethiopian Highlands are considered a biodiversity hotspot, harboring a high number of endemic species. Some of the endemic species probably diversified in situ; this is, for example, the case of a monophyletic clade containing 12 known species of grass frogs of the genus Ptychadena. The different species occur at elevations ranging from 1,500 to above 3,400 m and constitute excellent models to study the process of diversification in the highlands as well as adaptations to high elevations. In this study, we sampled 294 specimens across the distribution of this clade and used complete mitogenomes and genome-wide SNP data to better understand how landscape features influenced the population structure and dispersal of these grass frogs across time and space. Using phylogenetic inference, population structure analyses, and biogeographic reconstructions, we found that the species complex probably first diversified on the south-east side of the Great Rift Valley. Later on, species dispersed to the north-west side, where more recent diversification occurred. We further demonstrate that Ptychadena species have dispersed across the Great Rift Valley at different times. Our analyses allowed for a more complete understanding of the contribution of geological events, biogeographic barriers and climatic changes as drivers of species diversification and adaptation in this important biogeographic region.

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.

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

The conservation paradox of an introduced population of a threatened species: spadefoot toads in the coastal dunes of the Netherlands

Species that are threatened in their native range may actually prosper as introduced populations. To investigate how such introduced populations were established involves determining from where within the natural range the founder individuals originated. This can be accomplished through mtDNA barcoding. The common spadefoot toad (Pelobates fuscus) naturally occurs in the south and east of the Netherlands and has shown a rapid decline. Yet, a flourishing introduced population was recently discovered in the coastal dunes in the west of the country. We use mtDNA barcoding to determine the provenance of the introduced population. We sampled both native and introduced populations from the Netherlands and compared our sequences to haplotypes from across the entire distribution range. The mtDNA haplotypes found in the introduced population are distinct from those naturally occurring in the Netherlands and point towards an origin in the Pannonian Basin, on the boundary between Central and Southeastern Europe. Paradoxically, the thriving P. fuscus population in the Dutch coastal dunes should be considered a conservation risk to local biodiversity, even though within the native range in the Netherlands the species is severely threatened. Our study illustrates the complicated conservation questions associated with species that are both native and invasive.

Draft genome and multi-tissue transcriptome assemblies of the Neotropical leaf-frog Phyllomedusa bahiana

Amphibians are increasingly threatened worldwide, but the availability of genomic resources that could be crucial for implementing informed conservation practices lags well behind that for other vertebrate groups. Here, we describe draft de novo genome, mitogenome, and transcriptome assemblies for the Neotropical leaf-frog Phyllomedusa bahiana native to the Brazilian Atlantic Forest and Caatinga. We used a combination of PacBio long reads and Illumina sequencing to produce a 4.74-Gbp contig-level genome assembly, which has a contiguity comparable to other recent nonchromosome level assemblies. The assembled mitogenome comprises 16,239 bp and the gene content and arrangement are similar to other Neobratrachia. RNA-sequencing from 8 tissues resulted in a highly complete (86.3%) reference transcriptome. We further use whole-genome resequencing data from P. bahiana and from its sister species Phyllomedusa burmeisteri, to demonstrate how our assembly can be used as a backbone for population genomics studies within the P. burmeisteri species group. Our assemblies thus represent important additions to the catalog of genomic resources available from amphibians.

Historical biogeography highlights the role of Miocene landscape changes on the diversification of a clade of Amazonian tree frogs

The diversification processes underlying why Amazonia hosts the most species-rich vertebrate fauna on earth remain poorly understood. We studied the spatio-temporal diversification of a tree frog clade distributed throughout Amazonia (Anura: Hylidae: Osteocephalus, Tepuihyla, and Dryaderces) and tested the hypothesis that Miocene mega wetlands located in western and central Amazonia impacted connectivity among major biogeographic areas during extensive periods. We assessed the group’s diversity through DNA-based (16S rRNA) species delimitation to identify Operational Taxonomic Units (OTUs) from 557 individuals. We then selected one terminal for each OTU (n = 50) and assembled a mitogenomic matrix (~14,100 bp; complete for 17 terminals) to reconstruct a Bayesian, time-calibrated phylogeny encompassing nearly all described species. Ancestral area reconstruction indicates that each genus was restricted to one of the major Amazonian biogeographic areas (western Amazonia, Guiana Shield and Brazilian Shield, respectively) between ~10 and 20 Mya, suggesting that they diverged and diversified in isolation during this period around the Pebas mega wetland. After 10 Mya and the transition to the modern configuration of the Amazon River watershed, most speciation within each genus continued to occur within each area. In Osteocephalus, only three species expanded widely across Amazonia (< 6 Mya), and all were pond-breeders. Species with other breeding modes remained mostly restricted to narrow ranges. The spectacular radiation of Osteocephalus was probably driven by climatic stability, habitat diversity and the acquisition of new reproductive modes along the Andean foothills and western Amazonia. Our findings add evidence to the importance of major hydrological changes during the Miocene on biotic diversification in Amazonia.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

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.

Complete Mitochondrial Genome of Two Ectoparasitic Capsalids (Platyhelminthes: Monogenea: Monopisthocotylea): Gene Content, Composition, and Rearrangement

The capsalid monogeneans are important pathogens that generally infect marine fishes and have a substantial impact on fish welfare in aquaculture systems worldwide. However, the current mitogenome information on capsalids has received little attention, limiting the understanding of their evolution and phylogenetic relationships with other monogeneans. This paper reports the complete mitochondrial genomes of Capsala katsuwoni and Capsala martinieri for the first time, which we obtained using a next-generation sequencing method. The mitogenomes of C. katsuwoni and C. martinieri are 13,265 and 13,984 bp in length, respectively. Both species contain the typical 12 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The genome compositions show a moderate A+T bias (66.5% and 63.9% for C. katsuwoni and C. martinieri, respectively) and exhibit a negative AT skew but a positive GC skew in both species. One gene block rearrangement was found in C. katsuwoni in comparison with other capsalid species. Instead of being basal to the Gyrodactylidea and Dactylogyridea or being clustered with Dactylogyridea, all species of Capsalidea are grouped into a monophyletic clade. Our results clarify the gene rearrangement process and evolutionary status of Capsalidae and lay a foundation for further phylogenetic studies of monogeneans.

Systematic position of the Clicking Frog (Kassinula Laurent, 1940), the problem of chimeric sequences and the revised classification of the family Hyperoliidae

The systematics of the African frog family Hyperoliidae has undergone turbulent changes in last decades. Representatives of several genera have not been genetically investigated or with only limited data, and their phylogenetic positions are thus still not reliably known. This is the case of the De Witte's Clicking Frog (Kassinula wittei) which belongs to a monotypic genus. This miniature frog occurs in a poorly studied region, southeastern Democratic Republic of the Congo, northern Zambia, Angola. So far it is not settled whether this genus belongs to the subfamily Kassininae as a relative of the genus Kassina, or to the subfamily Hyperoliinae as a relative of the genus Afrixalus. Here we present for the first time a multilocus phylogenetic reconstruction (using five nuclear and one mitochondrial marker) of the family Hyperoliidae, including Kassinula. We demonstrate with high confidence that Kassinula is a member of Hyperoliinae belonging to a clade also containing Afrixalus (sub-Saharan Africa), Heterixalus (Madagascar) and Tachycnemis (Seychelles). We find that Kassinula represents a divergent lineage (17-25 Mya), which supports its separate genus-level status, but its exact systematic position remains uncertain. We propose to name the clade to which the above four genera belong as the tribe Tachycnemini Channing, 1989. A new taxonomy of the family Hyperoliidae was recently proposed by Dubois et al. (2021: Megataxa 5, 1-738). We demonstrate here that the new taxonomy was based on a partially erroneous phylogenetic reconstruction resulting from a supermatrix analysis of chimeric DNA sequences combining data from two families, Hyperoliidae and Arthroleptidae (the case of Cryptothylax). We therefore correct the erroneous part and propose a new, revised suprageneric taxonomy of the family Hyperoliidae. We also emphasize the importance of inspecting individual genetic markers before their concatenation or coalescent-based tree reconstructions to avoid analyses of chimeric DNA sequences producing incorrect phylogenetic reconstructions. Especially when phylogenetic reconstructions are used to propose taxonomies and systematic classifications.

Comparative Mitogenome Analysis of Gerbils and the Mitogenome Phylogeny of Gerbillinae (Rodentia: Muridae)

To enrich the mitogenomic database of Gerbillinae (Rodentia: Muridae), mitogenomes of three gerbils from different genera, Meriones tamariscinus (16,393 bp), Brachiones przewalskii (16,357 bp), and Rhombomys opimus (16,352 bp), were elaborated and compared with those of other gerbils in the present study. The three gerbil mitogenomes consisted of 2 ribosomal RNA genes, 13 protein-coding genes (PCGs), 22 transfer RNA genes, and one control region. Here, gerbil mitogenomes have shown unique characteristics in terms of base composition, codon usage, non-coding region, and the replication origin of the light strand. There was no significant correlation between the nucleotide percentage of G + C and the phylogenetic status in gerbils, and between the GC content of PCGs and the leucine count. Phylogenetic relationships of the subfamily Gerbillinae were reconstructed by 7 gerbils that represented four genera based on concatenated mitochondrial DNA data using both Bayesian Inference and Maximum Likelihood. The phylogenetic analysis indicated that M. tamariscinus was phylogenetically distant from the genus Meriones, but has a close relationship with R. opimus. B. przewalskii was closely related to the genus Meriones rather than that of R. opimus.

A new toad species of the genus Brachytarsophrys Tian & Hu, 1983 (Anura, Megophryidae) from Guizhou Province, China

Background

The toads of the genus Brachytarsophrys Tian & Hu, 1983 are distributed in southern China, Myanmar, Vietnam, Laos and northern Thailand. Seven species of the genus have been recognised, of which five of them are known from China so far.

New information

Brachytarsophrysqiannanensis sp. nov., a new species of the short-legged toad genus is here described from southern Guizhou Province, China. Diagnostic characters of the new species are illustrated and comparisons with its congeners are provided. Its validity is also affirmed by its distinct mitochondrial gene sequence divergence with all congeners and its monophyly recovered in the mitochondrial gene-based phylogenetic analyses.

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.

Initial Phylotranscriptomic Confirmation of Homoplastic Evolution of the Conspicuous Coloration and Bufoniform Morphology of Pumpkin-Toadlets in the Genus Brachycephalus

The genus Brachycephalus is a fascinating group of miniaturized anurans from the Brazilian Atlantic Forest, comprising the conspicuous, brightly colored pumpkin-toadlets and the cryptic flea-toads. Pumpkin-toadlets are known to contain tetrodotoxins and therefore, their bright colors may perform an aposematic function. Previous studies based on a limited number of mitochondrial and nuclear-encoded markers supported the existence of two clades containing species of pumpkin-toadlet phenotype, but deep nodes remained largely unresolved or conflicting between data sets. We use new RNAseq data of 17 individuals from nine Brachycephalus species to infer their evolutionary relationships from a phylogenomic perspective. Analyses of almost 5300 nuclear-encoded ortholog protein-coding genes and full mitochondrial genomes confirmed the existence of two separate pumpkin-toadlet clades, suggesting the convergent evolution (or multiple reversals) of the bufoniform morphology, conspicuous coloration, and probably toxicity. In addition, the study of the mitochondrial gene order revealed that three species (B. hermogenesi, B. pitanga, and B. rotenbergae) display translocations of different tRNAs (NCY and CYA) from the WANCY tRNA cluster to a position between the genes ATP6 and COIII, showing a new mitochondrial gene order arrangement for vertebrates. The newly clarified phylogeny suggests that Brachycephalus has the potential to become a promising model taxon to understand the evolution of coloration, body plan and toxicity. Given that toxicity information is available for only few species of Brachycephalus, without data for any flea-toad species, we also emphasize the need for a wider screening of toxicity across species, together with more in-depth functional and ecological study of their phenotypes.

Germline competent mesoderm: the substrate for vertebrate germline and somatic stem cells?

In vitro production of tissue-specific stem cells [e.g. haematopoietic stem cells (HSCs)] is a key goal of regenerative medicine. However, recent efforts to produce fully functional tissue-specific stem cells have fallen short. One possible cause of shortcomings may be that model organisms used to characterize basic vertebrate embryology (Xenopus, zebrafish, chick) may employ molecular mechanisms for stem cell specification that are not conserved in humans, a prominent example being the specification of primordial germ cells (PGCs). Germ plasm irreversibly specifies PGCs in many models; however, it is not conserved in humans, which produce PGCs from tissue termed germline-competent mesoderm (GLCM). GLCM is not conserved in organisms containing germ plasm, or even in mice, but understanding its developmental potential could unlock successful production of other stem cell types. GLCM was first discovered in embryos from the axolotl and its conservation has since been demonstrated in pigs, which develop from a flat-disc embryo like humans. Together these findings suggest that GLCM is a conserved basal trait of vertebrate embryos. Moreover, the immortal nature of germ cells suggests that immortality is retained during GLCM specification; here we suggest that the demonstrated pluripotency of GLCM accounts for retention of immortality in somatic stem cell types as well.

This article has an associated Future Leaders to Watch interview with the author of the paper.

Summary: Recent findings that germline and stem cell specification may differ between species may have important implications for regenerative medicine and the future of stem cell biology.

Unexpectedly high levels of lineage diversity in Sundaland puddle frogs (Dicroglossidae: Occidozyga Kuhl and van Hasselt, 1822)

One of the most urgent contemporary tasks for taxonomists and evolutionary biologists is to estimate the number of species on earth. Recording alpha diversity is crucial for protecting biodiversity, especially in areas of elevated species richness, which coincide geographically with increased anthropogenic environmental pressures - the world's so-called biodiversity hotspots. Although the distribution of Puddle frogs of the genus Occidozyga in South and Southeast Asia includes five biodiversity hotspots, the available data on phylogeny, species diversity, and biogeography are surprisingly patchy. Samples analyzed in this study were collected throughout Southeast Asia, with a primary focus on Sundaland and the Philippines. A mitochondrial gene region comprising ~ 2000 bp of 12S and 16S rRNA with intervening tRNA Valine and three nuclear loci (BDNF, NTF3, POMC) were analyzed to obtain a robust, time-calibrated phylogenetic hypothesis. We found a surprisingly high level of genetic diversity within Occidozyga, based on uncorrected p-distance values corroborated by species delimitation analyses. This extensive genetic diversity revealed 29 evolutionary lineages, defined by the > 5% uncorrected p-distance criterion for the 16S rRNA gene, suggesting that species diversity in this clade of phenotypically homogeneous forms probably has been underestimated. The comparison with results of other anuran groups leads to the assumption that anuran species diversity could still be substantially underestimated in Southeast Asia in general. Many genetically divergent lineages of frogs are phenotypically similar, indicating a tendency towards extensive morphological conservatism. We present a biogeographic reconstruction of the colonization of Sundaland and nearby islands which, together with our temporal framework, suggests that lineage diversification centered on the landmasses of the northern Sunda Shelf. This remarkably genetically structured group of amphibians could represent an exceptional case for future studies of geographical structure and diversification in a widespread anuran clade spanning some of the most pronounced geographical barriers on the planet (e.g., Wallace's Line). Studies considering gene flow, morphology, ecological and bioacoustic data are needed to answer these questions and to test whether observed diversity of Puddle frog lineages warrants taxonomic recognition.

The complete mitogenome of the Mountain chicken frog, Leptodactylus fallax

The mountain chicken frog (Leptodactylus fallax) is a critically endangered frog native to the Caribbean islands of Dominica and Montserrat. Over the past 25 years their populations have declined by over 85%, largely due to a chytridiomycosis outbreak that nearly wiped out the Montserratian population. Within the context of developing tools that can aid in the conservation of the mountain chicken frog, we assembled its complete mitochondrial genome, contributing the first complete mitogenome of the genus Leptodactylus (Genbank Accession number MW260634). The circular genome is 18,669 bp long and contains 37 genes. A phylogenetic analysis reveals that L. fallax forms a clade with Leptodactylus melanonotus, highlighting the close relationship of Leptodactylus spp. relative to other species from the superfamily Hyloidea included in the analysis.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

Nanopore Sequencing Resolves Elusive Long Tandem-Repeat Regions in Mitochondrial Genomes

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

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.

Three new species of frogs of the genus Pristimantis (Anura, Strabomantidae) with a redefinition of the P. lacrimosus species group

A new phylogeny for the Pristimantis lacrimosus species group is presented, its species content reviewed, and three new species described from the eastern slopes of the Ecuadorian Andes. Our phylogeny includes, for the first time, samples of P. aureolineatus, P. bromeliaceus, and P. lacrimosus. The morphology of hyperdistal subarticular tubercles is also assessed among 21 species of Pristimantis. The P. lacrimosus species group is composed of 36 species distributed in the Chocó, Guiana, and Amazon regions of tropical South America with a single species reaching Central America. Ancestral area reconstruction indicates that, despite its high diversity in the Amazon region, the P. lacrimosus group originated in the Pacific basin, Chocó region of Ecuador and Colombia. Pristimantis amaguanae sp. nov. is most closely related to P. bromeliaceus. It differs from P. bromeliaceus by being smaller, having transversal dark bands in the hindlimbs (absent or faint in P. bromeliaceus) and the absence of discoidal fold (present in P. bromeliaceus). Pristimantis nankints sp. nov. and P. romeroae sp. nov. are part of a clade of predominantly light-green frogs that includes P. acuminatus, P. enigmaticus, P. limoncochensis, and P. omeviridis. Pristimantis nankints sp. nov. and P. romeroae sp. nov. can be distinguished from all of them by the presence of a dark dorsolateral stripe that borders a light green band on a green background. Hyperdistal tubercles are present in all examined species of the P. lacrimosus species group and its sister clade. Species with hyperdistal tubercles are characterized by having relatively long terminal phalanges and narrow T-shaped expansion at the end of the terminal phalange. We discuss the phylogenetic distribution of these characters and their potential diagnostic significance.

Genetic structure and demographic history of Indirana semipalmata, an endemic frog species of the Western Ghats, India

The evolutionary potential of a species mainly depends on the level of genetic variation in their populations. Maintenance of gene variation enables populations to adapt more quickly to environmental changes. The geographical gaps also influence the distribution and evolutionary history of many mountain frogs in the world. Hence, a sound knowledge in population genetic structure of a species will help understand its population dynamics and develop conservation strategies. In the context of facing threats to the amphibian fauna of Western Ghats due to habitat loss, we used both mitochondrial and nuclear DNA markers to investigate the genetic structure of an endemic frog species of the Western Ghats (Indirana semipalmata) with restricted distribution. The present study showed the importance of mountain gaps in shaping the species' structuring in the Western Ghats. Though a high genetic diversity was observed for the species when considering a single unit in the southern Western Ghats, the restricted gene flow on/between either side of the Shencottah gap with genetic clustering of the sampled populations may warrant a unique management plan for the species. The habitat fragmentation of the Western Ghats through anthropogenic activities may result in severe setbacks to the survival of the species in the future.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Review of the genus Brachytarsophrys (Anura: Megophryidae), with revalidation of Brachytarsophrys platyparietus and description of a new species from China

The genus-level recognition of monophyletic short-legged toads (Brachytarsophrys) has been recently implicated in the taxonomic debate of Megophrys sensu lato. In the present study, Brachytarsophrys is reasonably regarded as a distinct genus based on significant morphological differentiations and recent molecular analyses. Furthermore, a comprehensive review of this genus is performed, with two species groups proposed based on morphological differences and phylogenetic relationships. Particularly, Brachytarsophrys platyparietus is removed as a synonym of Brachytarsophrys carinense and considered a valid species due to significant genetic divergence and distinct morphological differences. In addition, a new species, Brachytarsophrys orientalissp. nov., is described based on a series of specimens collected from southeastern China. This work takes the member species of the genus Brachytarsophrys to seven, suggesting that the diversity of Brachytarsophrys is underestimated. In addition, the genus levels of other monophyletic groups within the subfamily Megophryinae are discussed.

A new species of the toothed toad Oreolalax (Anura, Megophryidae) from Sichuan Province, China

The toad genus Oreolalax is widely distributed in southwest China and northern Vietnam. A new species of the genus is described from Sichuan Province, China. Phylogenetic analyses based on the mitochondrial 12S rRNA and 16S rRNA gene sequences supported the new species as an independent clade clustered into the clade also containing O. nanjiangensis and O. chuanbeiensis. The new species can be distinguished from its congeners by a combination of the following characters: body size moderate (SVL 51.2–64.2 mm in males); head broad; tympanum hidden; interorbital region with dark triangular pattern; belly with marbling; male lacking spines on lip margin; spiny patches on chest small with thick sparse spines in male; nuptial spines thick and sparse; tibio-tarsal articulation reaching beyond nostril when leg stretched forward; toe webbing at base.

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.

Exceptional Enlargement of the Mitochondrial Genome Results from Distinct Causes in Different Rain Frogs (Anura: Brevicipitidae: Breviceps)

The mitochondrial (mt) genome of the bushveld rain frog (Breviceps adspersus, Brevicipitidae, Afrobatrachia) is the largest (28.8 kbp) among the vertebrates investigated to date. The major cause of genome size enlargement in this species is the duplication of multiple genomic regions. To investigate the evolutionary lineage, timing, and process of mt genome enlargement, we sequenced the complete mtDNAs of two congeneric rain frogs, B. mossambicus and B. poweri. The mt genomic organization, gene content, and gene arrangements of these two rain frogs are very similar to each other but differ from those of B. adspersus. The B. mossambicus mt genome (22.5 kbp) does not differ significantly from that of most other afrobatrachians. In contrast, the B. poweri mtDNA (28.1 kbp) is considerably larger: currently the second largest among vertebrates, after B. adspersus. The main causes of genome enlargement differ among Breviceps species. Unusual elongation (12.5 kbp) of the control region (CR), a single major noncoding region of the vertebrate mt genome, is responsible for the extremely large mt genome in B. poweri. Based on the current Breviceps phylogeny and estimated divergence age, it can be concluded that the genome enlargements occurred independently in each species lineage within relatively short periods. Furthermore, a high nucleotide substitution rate and relaxation of selective pressures, which are considered to be involved in changes in genome size, were also detected in afrobatrachian lineages. Our results suggest that these factors were not direct causes but may have indirectly affected mt genome enlargements in Breviceps.

Sperm collection and storage for the sustainable management of amphibian biodiversity

Current rates of biodiversity loss pose an unprecedented challenge to the conservation community, particularly with amphibians and freshwater fish as the most threatened vertebrates. An increasing number of environmental challenges, including habitat loss, pathogens, and global warming, demand a global response toward the sustainable management of ecosystems and their biodiversity. Conservation Breeding Programs (CBPs) are needed for the sustainable management of amphibian species threatened with extinction. CBPs support species survival while increasing public awareness and political influence. Current CBPs only cater for 10% of the almost 500 amphibian species in need. However, the use of sperm storage to increase efficiency and reliability, along with an increased number of CBPs, offer the potential to significantly reduce species loss. The establishment and refinement of techniques over the last two decades, for the collection and storage of amphibian spermatozoa, gives confidence for their use in CBPs and other biotechnical applications. Cryopreserved spermatozoa has produced breeding pairs of frogs and salamanders and the stage is set for Lifecycle Proof of Concept Programs that use cryopreserved sperm in CBPs along with repopulation, supplementation, and translocation programs. The application of cryopreserved sperm in CBPs, is complimentary to but separate from archival gene banking and general cell and tissue storage. However, where appropriate amphibian sperm banking should be integrated into other global biobanking projects, especially those for fish, and those that include the use of cryopreserved material for genomics and other research. Research over a broader range of amphibian species, and more uniformity in experimental methodology, is needed to inform both theory and application. Genomics is revolutionising our understanding of biological processes and increasingly guiding species conservation through the identification of evolutionary significant units as the conservation focus, and through revealing the intimate relationship between evolutionary history and sperm physiology that ultimately affects the amenability of sperm to refrigerated or frozen storage. In the present review we provide a nascent phylogenetic framework for integration with other research lines to further the potential of amphibian sperm banking.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Multi-Tissue Transcriptomes Yield Information on High-Altitude Adaptation and Sex-Determination in Scutiger cf. sikimmensis

The Himalayas are one of earth’s hotspots of biodiversity. Among its many cryptic and undiscovered organisms, including vertebrates, this complex high-mountain ecosystem is expected to harbour many species with adaptations to life in high altitudes. However, modern evolutionary genomic studies in Himalayan vertebrates are still at the beginning. Moreover, in organisms, like most amphibians with relatively high DNA content, whole genome sequencing remains bioinformatically challenging and no complete nuclear genomes are available for Himalayan amphibians. Here, we present the first well-annotated multi-tissue transcriptome of a Greater Himalayan species, the lazy toad Scutiger cf. sikimmensis (Anura: Megophryidae). Applying Illumina NextSeq 500 RNAseq to six tissues, we obtained 41.32 Gb of sequences, assembled to ~111,000 unigenes, translating into 54362 known genes as annotated in seven functional databases. We tested 19 genes, known to play roles in anuran and reptile adaptation to high elevations, and potentially detected diversifying selection for two (TGS1, SENP5) in Scutiger. Of a list of 37 genes, we also identify 27 candidate genes for sex determination or sexual development, all of which providing the first such data for this non-model megophryid species. These transcriptomes will serve as a valuable resource for further studies on amphibian evolution in the Greater Himalaya as a biodiversity hotspot.

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.

Temporin L and aurein 2.5 have identical conformations but subtly distinct membrane and antibacterial activities

Frogs such as Rana temporaria and Litoria aurea secrete numerous closely related antimicrobial peptides (AMPs) as an effective chemical dermal defence. Damage or penetration of the bacterial plasma membrane is considered essential for AMP activity and such properties are commonly ascribed to their ability to form secondary amphipathic, α-helix conformations in membrane mimicking milieu. Nevertheless, despite the high similarity in physical properties and preference for adopting such conformations, the spectrum of activity and potency of AMPs often varies considerably. Hence distinguishing apparently similar AMPs according to their behaviour in, and effects on, model membranes will inform understanding of primary-sequence-specific antimicrobial mechanisms. Here we use a combination of molecular dynamics simulations, circular dichroism and patch-clamp to investigate the basis for differing anti-bacterial activities in representative AMPs from each species; temporin L and aurein 2.5. Despite adopting near identical, α-helix conformations in the steady-state in a variety of membrane models, these two AMPs can be distinguished both in vitro and in silico based on their dynamic interactions with model membranes, notably their differing conformational flexibility at the N-terminus, ability to form higher order aggregates and the characteristics of induced ion conductance. Taken together, these differences provide an explanation of the greater potency and broader antibacterial spectrum of activity of temporin L over aurein 2.5. Consequently, while the secondary amphipathic, α-helix conformation is a key determinant of the ability of a cationic AMP to penetrate and disrupt the bacterial plasma membrane, the exact mechanism, potency and spectrum of activity is determined by precise structural and dynamic contributions from specific residues in each AMP sequence.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

Identification of leopard frogs (Anura: Ranidae: Lithobates) distributed in some localities of the Southern Mexican Plateau using mitochondrial DNA sequences

Species of the genus Lithobates are difficult to identify, especially on the 'Rana pipiens' complex, because the morphological differences among the species are often subtle. In this study, we used 12S ribosomal RNA gene partial sequences to identify specimens of leopard frogs from some new localities (not sampled on previous studies) of the Southern Mexican Plateau, to carry out a phylogenetic analysis, and also a morphometric analysis of some morphological features were analyzed to evaluate their morphological variation. A phylogenetic analysis using partial sequences of 12S rDNA mitochondrial gene for 31 samples was carried out using Bayesian Inference, Maximum Parsimony, and Maximum Likelihood. In addition, 20 morphological lineal measurements from 97 specimens were analyzed by principal component analysis (PCA) to evaluate if the frogs studied are conspecific. Partial sequences of the 12S rDNA obtained from frogs distributed in the Southern Mexican Plateau show two haplotypes with a divergence of 0.4%. Phylogenetic hypothesis shows an exclusive group with a previous sequence of Lithobates montezumae. The PCA indicates that variables are not linearly correlated and specimens belong to a single group. Evidence found, let us consider that specimens studied belong to one species of the Lithobates montezumae subgroup, distributed principally in the Southern Mexican Plateau.