Genetic Divergences and Phylogenetic Relationships Among the Fejervarya limnocharis Complex in Thailand and Neighboring Countries Revealed by Mitochondrial and Nuclear Genes

An insight into molecular taxonomy of bufonids, microhylids, and dicroglossid frogs: First genetic records from Pakistan

The current study was focused on documentation of amphibian assemblage in North Punjab and Islamabad Capital Territory, Pakistan, by using mitochondrial gene sequences of 16S rRNA. Our study entailed 37% of the known amphibian species of the country. We provided a phylogenetic analysis based on 74 newly generated mitochondrial 16S rRNAs from nine species of genus Microlyla, Duttaphrynus, Allopaa, Nanorana, Sphaerotheca, Minervarya, Hoplobatrachus, and Euphlyctis. We employed the maximum-likelihood inference and Bayesian analysis to assess the taxonomic status of the samples obtained from Pakistan, with respect to other congeneric species from surrounding regions. Our findings confirmed the taxonomic status of South Asian anuran species Duttaphrynus stomaticus, Duttaphrynus melanostictus, Microhyla nilphamariensis, Allopaa hazarensis, Nanorana vicina, Sphaerotheca maskeyi (synonym: S. pashchima), Minervarya pierrei, Hoplobatrachus tigerinus, and Euphlyctis kalasgramensis in Pakistan. We have reported new country records of genus Minervarya ( M. pierrei). Minervarya pierrei was previously misidentified as Fejervarya limnocharis, due to dearth of genetic information. We provided the first genetic records of our endemic species N. vicina. The results revealed the taxonomic placement of N. vicina with respect to its congeners and validated the taxonomic status of N. vicina from its type locality (Murree) for the first time. The findings of the present study also indicated the paraphyletic relationship of A.- hazarensis with Nanorana species. So, based on our phylogenetic inferences, morphological characters, and habitat preferences, validity of generic status of A. hazarensis is undecided. As our data were not enough to resolve this issue, we suggest sequencing of additional mitochondrial and nuclear genes in the future studies to get a better resolution. We recommend carrying out extensive surveys throughout the country for proper scientific documentation of amphibians of Pakistan. Many new species, some of them might be endemic to Pakistan, are expected to be discovered, and taxonomic status of other species would be resolved.

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.

Species delimitation of crab-eating frogs (Fejervarya cancrivora complex) clarifies taxonomy and geographic distributions in mainland Southeast Asia

The taxonomy and geographic distributions of species of crab-eating frogs (Fejervaryacancrivora complex) in mainland Southeast Asia have been highly uncertain. Three taxonomic names are used in recent literature (F.cancrivora, F.raja, and F.moodiei) but the applications of these names to localities has been inconsistent, especially owing to the lack of available molecular data for F.raja. Morphometric and mitochondrial DNA variation was examined in these frogs, including name-bearing types and topotypes of all three species. Findings corroborate evidence for the existence of two species in coastal mainland Southeast Asia, with F.moodiei having a wide geographic distribution and F.cancrivora sensu stricto occurring only in extreme southern Thailand and peninsular Malaysia. Fejervaryaraja is shown to be only a large-bodied population of F.cancrivora sensu stricto and is synonymized with that species. Revised descriptions of F.moodiei and F.cancrivora sensu stricto are provided.

A new species of genus Fejervarya (Anura: Dicroglossidae) from northern Thailand

We describe a new species of frog in the dicroglossid genus Fejervarya from Ban Monjong, Omkoi District, Chiang Mai Province, northern Thailand. Analysis of DNA sequence data from the mitochondrial gene 16S, advertisement calls, and morphological distinctiveness support recognition of the new species. Matrilineal genealogy suggests that the new population from Chiang Mai is a sister taxon to the South Asian clade that includes F. syhadrensis, F. granosa, and F. pierrei. The new species, Fejervarya chiangmaiensissp. nov., differs morphologically from its congeners by its relatively small body size and proportions and the presence of dorsal warts and dermal ridges. Discovery of this new species indicates that the biodiversity of amphibians in this region remains underestimated.

A New Species of Frog (Anura: Dicroglossidae) Discovered from the Mega City of Dhaka

We describe a new species of frog of the genus Zakerana discovered from the urban core of Dhaka, Bangladesh, one of the most densely populated cities in the world. Although the new species is morphologically similar to the geographically proximate congeners in the Bangladeshi cricket frog group, we show that it can be distinguished from all congeners on the basis of morphological characters, advertisement calls and variation in two mitochondrial DNA genes (12S rRNA and 16S rRNA). Apart from several diagnostic differences in body proportions, the new species differs from other Zakerana species in having a flattened snout (from ventral view) projecting over the lower jaw, and diagnostic trapezoid-shaped red markings on the vocal sac in males. Molecular genetic analyses show that the new species is highly divergent (3.1–20.1% sequence divergence) from all congeneric species, and forms a well-supported clade with its sister species, Zakerana asmati. The discovery of a new amphibian species from the urban core of Dhaka together with several recent descriptions of new amphibian species from Bangladesh may indicate that more amphibian species remain to be discovered from this country.

Cryptic anuran biodiversity in Bangladesh revealed by mitochondrial 16S rRNA gene sequences

To survey the diversity of anuran species in Bangladesh, we compared mitochondrial 16S rRNA gene sequences (approximately 1.4 kbp) from 107 Bangladesh frog specimens. The results of genetic divergence and phylogenetic analyses incorporating data from related species revealed the occurrence of at least eight cryptic species. Hoplobatrachus tigerinus from two districts diverged considerably, indicating the involvement of a cryptic species. Two Fejervarya sp. (large and medium types) and Hylarana cf. taipehensis formed lineages distinct from related species and are probably new species. Microhyla cf. ornata differed from M. ornata with respect to type locality area and involved two distinct species. In addition, we found that Hylarana sp. and Microhyla sp. did not match congeners examined to date in either morphology or 16S rRNA sequence. The occurrence of M. fissipes was tentatively suggested. Consequently, at least, 19 species were found from Bangladesh in this study. These findings revealed a rich anuran biodiversity in Bangladesh, which is unexpected considering the rather simple topographic features of the country.

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.