Intraspecific Differentiation in the Japanese Brown Frog Rana japonica Inferred from Mitochondrial DNA Sequences of the CytochromebGene

Genetic diversity and demography of Bufo japonicus and B. torrenticola (Amphibia: Anura: Bufonidae) influenced by the Quaternary climate

The Quaternary climate affected the present species richness and geographic distribution patterns of amphibians by limiting their activities during the glacial period. The present study examined the phylogenetic relationships of Japanese toads (Bufo japonicus and B. torrenticola) and the demography of each lineage from the past to the present based on mitochondrial sequences and ecological niche models. Japanese toads are a monophyletic group with two main clades (clades A and B). Clade A represents B. j. formosus, including three clades (clades A1, A2, and A3). Clade B contains three clades, two of which corresponded to B. j. japonicus (clades B1 and B2) and the other to B. torrenticola. Clade B2 and B. torrenticola made a sister group, and, thus, B. j. japonicus is paraphyletic. Clades A and B diverged in the late Miocene 5.7 million years ago (Mya) during the period when the Japanese archipelago was constructed. The earliest divergence between the three clades of clade A was estimated at 1.8 Mya. Clades A1 and A2 may have diverged at 0.8 Mya, resulting from the isolation in the multiple different refugia; however, the effects of the glacial climate on the divergence events of clade A3 are unclear. Divergences within clade B occurred from the late Pliocene to the early Pleistocene (3.2-2.2 Mya). Niche similarity between the parapatric clade in clade B (clades B1 and B2) indicated their allopatric divergence. It was suggested that niche segregation between B. japonicus and B. torrenticola contributed to a rapid adaptation of B. torrenticola for lotic breeding. All clade of Japanese toads retreated to each refugium at a low elevation in the glacial period, and effective population sizes increased to construct the current populations after the Last Glacial Maximum. Furthermore, we highlight the areas of climate stability from the last glacial maximum to the present that have served as the refugia of Japanese toads and, thus, affected their present distribution patterns.

Diversity, distribution and molecular species delimitation in frogs and toads from the Eastern Palaearctic

Biodiversity analyses can greatly benefit from coherent species delimitation schemes and up-to-date distribution data. In this article, we have made the daring attempt to delimit and map described and undescribed lineages of anuran amphibians in the Eastern Palaearctic (EP) region in its broad sense. Through a literature review, we have evaluated the species status considering reproductive isolation and genetic divergence, combined with an extensive occurrence dataset (nearly 85k localities). Altogether 274 native species from 46 genera and ten families were retrieved, plus eight additional species introduced from other realms. Independent hotspots of species richness were concentrated in southern Tibet (Medog County), the circum-Sichuan Basin region, Taiwan, the Korean Peninsula and the main Japanese islands. Phylogeographic breaks responsible for recent in situ speciation events were shared around the Sichuan Mountains, across Honshu and between the Ryukyu Island groups, but not across shallow water bodies like the Yellow Sea and the Taiwan Strait. Anuran compositions suggested to restrict the zoogeographical limits of the EP to East Asia. In a rapidly evolving field, our study provides a checkpoint to appreciate patterns of species diversity in the EP under a single, spatially explicit, species delimitation framework that integrates phylogeographic data in taxonomic research.

Evidence from the primary structures of dermal antimicrobial peptides that Rana tagoi okiensis and Rana tagoi tagoi (Ranidae) are not conspecific subspecies

Morphological evidence and data from comparisons of nucleotide sequences of mitochondrial genes demonstrate considerable intraspecies variation among populations of the Japanese brown frog Rana tagoi Okada 1928 (Tago's brown frog). Five peptides with antimicrobial activity were isolated from an extract of the skins of specimens of Rana tagoi okiensis collected on the Oki Islands, Japan. Determination of their primary structures demonstrated that two peptides belong to the ranatuerin-2 family, two peptides to the temporin family, and one peptide to the brevinin-1 family. Ranatuerin-2 peptides were not previously identified in the skin of specimens of R. t. tagoi collected in Chiba Prefecture, Japan and the structures of the temporin peptides from R. t. okiensis (temporin-TOa: FLPILGKLLSGFL.NH(2) and temporin-TOb: FLPILGKLLSGLL.NH(2)) are different from temporin-TGa (FLPILGKLLSGIL.NH(2)) isolated from R. t. tagoi. Similarly, the acyclic C-terminally alpha-amidated brevinin-1 peptide from R. t. okiensis (Brevinin-1TOa, GIGSILGVIAKGLPTLISWIKNR.NH(2)) shows three amino acid substitutions (Gly(1)-->Ala, Val(8)-->Ala, Ile(9)-->Leu) compared to the ortholog from R. t. tagoi. In addition, bradykinin, identical to the mammalian peptide, is present in high concentration in the skin of R. t. okiensis but not R. t. tagoi. The data provide evidence to support the proposal that R. t. tagoi and R. t. okiensis should be regarded as separate species (R. tagoi and R. okiensis) rather than conspecific subspecies.

Phylogenetic relationships of the Chinese brown frogs (genus Rana) inferred from partial mitochondrial 12S and 16S rRNA gene sequences

Based on partial sequences of the 12S and 16S ribosomal RNA genes, we estimated phylogenetic relationships among brown frogs of the Rana temporaria group from China. From the phylogenetic trees obtained, we propose to include Rana zhengi in the brown frogs. Monophyly of the brown frogs was not unambiguously supported, but four well-supported clades (A, B, C, and D) always emerged, although relationships among them remained unresolved. Clade A contained brown frogs with 24 chromosomes and was split into two distinct subclades (Subclade A-1: R. chensinensis and R. huanrenensis; Subclade A-2: R. dybowskii). Polytomous relationships among populations of R. chensinensis and R. huanrenensis suggested the necessity of further taxonomic assessment. Rana kunyuensis proved to be the sister group to R. amurensis, and these two species formed Clade B. Clade C was composed of R. omeimontis and R. chaochiaoensis, and Clade D included R. sauteri, which has been placed in other ranid genera. These relationships did not change after adding published data, and monophyly of Subclade A-1, A-2, and other East Asian brown frogs with 24 chromosomes (R. pirica and R. ornativentris) was ascertained, though their relationships were unresolved. Clade C, together with R. japonica and R. longicrus, also formed a monophyletic group. Brown frogs related to Clades A and C were estimated to have dispersed from continental Asia to adjacent regions through multiple events.

Cytolytic peptides belonging to the brevinin-1 and brevinin-2 families isolated from the skin of the Japanese brown frog, Rana dybowskii

Peptidomic analysis of an extract of the skins of specimens of Dybowski's brown frog Rana dybowskii Gunther, 1876, collected on Tsushima Island, Japan led to the identification of 10 peptides with differential antibacterial and hemolytic activities. The primary structures of these peptides identified them as belonging to the brevinin-1 (5 peptides) and brevinin-2 (5 peptides) families of antimicrobial peptides. A peptide (FIGPIISALASLFG.NH(2)) with structural similarity to members of the temporin family was also isolated but this component lacked cytolytic activity. Phylogenetic relationships among the Japanese brown frogs (R. dybowskii, R. japonica, R. okinavana, R. ornativentris, R. pirica, R. sakuraii, R. tagoi, and R. tsushimensis) are only incompletely understood. Cladograms based upon maximum parsimony analyses of the brevinin-1 and brevinin-2 amino acid sequences provide strong support for a sister-group relationship between R. dybowskii and R. pirica and somewhat weaker support for a sister-group relationship between R. okinavana and R. tsushimensis. These conclusions are consistent with previous analyses based upon allozyme variations and comparisons of the nucleotide sequences of mitochondrial genes.

Molecular phylogenetic relationships of pond frogs distributed in the Palearctic region inferred from DNA sequences of mitochondrial 12S ribosomal RNA and cytochrome b genes

The evolutionary relationships of pond frogs distributed in the Far East and Europe were investigated by analyses of nucleotide sequences of mitochondrial 12S ribosomal RNA (12S rRNA) and cytochrome b (cyt b) genes. The nucleotide sequences of a 412-bp segment of the 12S rRNA gene and a 534-bp segment of the cyt b gene were determined by the PCR-direct sequencing method using 19 frogs belonging to six species and one subspecies distributed in the Palearctic region. Phylogenetic trees were constructed by the neighbor-joining and maximum-likelihood methods using Rana catesbeiana or Xenopus laevis as an outgroup. The 412-bp segment of the 12S rRNA gene contained 65 variable sites including gap sites, and the 534-bp segment of the cyt b gene contained 160 variable sites. The nucleotide sequence divergences of the 12S rRNA gene were 0.25-4.83% within the Far Eastern frogs, 0.25-6.22% within the European frogs, and 8.74-11.24% between the Far Eastern and the European frogs, whereas those of the cyt b gene were 3.64-14.73% within the Far Eastern frogs, 0.38-14.42% within the European frogs, and 16.53-23.58% between the Far Eastern and the European frogs. Although most nucleotide substitutions were at the third codon position of the cyt b gene and were silent mutations, 4 amino acid replacements occurred within the Far Eastern frogs, 4 within the European frogs, and 11 between the Far Eastern and the European frogs. The phylogenetic trees constructed from the nucleotide sequence divergences showed slightly different topologies for the 12S rRNA and cyt b genes. R. esculenta from Ukraine was closely related to R. lessonae from Luxembourg in both the 12S rRNA and the cyt b gene sequences.

Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.