Evolutionary relationships of Eastern Palearctic Brown Frogs, genus Rana: paraphyly of the 24-chromosome species group and the significance of chromosome number change

POSTGLACIAL RANGE FLUCTUATION, GENETIC SUBDIVISION AND SPECIATION IN THE WESTERN NORTH AMERICAN SPOTTED FROG COMPLEX, RANA PRETIOSA

The western North American complex of spotted frogs (Rana pretiosa) exhibits isolation-by-distance, genetic subdivision, and speciation in association with its extensive northward range shift in postglacial times. The southern relict populations of R. pretiosa species B existing at high altitudes or in desert springs have been subjected to restricted gene flow, high inbreeding, and bottlenecks to produce significant between-population genetic diversity. The more recently established northern populations, however, show genetic uniformity and isolation-by-distance, as estimated using Slatkin's (1993) statistic M̂. Middle latitude populations have higher heterozygosities than populations at either extreme. Fixed differences in allozyme variation separate 21 populations of species B from five populations of R. pretiosa species A found in southwest Washington State and the Cascades Mountains of Oregon. Morphological variation of 20 metric characters among 38 samples, examined using multiple discriminant function analysis, could partially resolve partitioning among populations but specimens from the vicinity of the type series of R. p. pretiosa could not be assigned to either species A or species B. Speciation in these frogs may not be correlated with morphological evolution since comparatively neutral allozyme changes may be established more rapidly than changes in morphology.

Molecular phylogenetic relationships among Anatolian-Hyrcanian brown frog taxa (Ranidae: Rana)

Although the phylogenetic relationship of Western Palearctic brown frogs has been repeatedly studied, the taxonomic status and phylogenetic relationship of Anatolian-Hyrcanian brown frogs is still not fully resolved. Here, we assess the phylogenetic status of these species among Western Palearctic brown frogs with special emphasize on Iranian populations based on two partial mitochondrial DNA sequences (16S rRNA and cytochrome b genes) and the application of a molecular clock. Our results clearly show that Western Palearctic brown frogs underwent a basal radiation in to two main monophyletic clades, the European brown frogs plus the Asian R. asiatica and the Anatolian-Hyrcanian brown frogs, during Early Miocene ca. 20.2 mya. The Hyrcanian (R. pseudodalmatina) and the Anatolian lineage diverged approximately 16.6 mya. The further diverged into two subclades, R. tavasensis and R. macrocnemis, during the Middle Miocene, 14.5 mya. Our results suggest that diversification within these lineages may be closely linked to the formation of Neotethys and Paratethys and the subsequent uplift of the Turkish-Iranian plateau during the Early Miocene which led to restricted gene flow among brown frogs in these regions. Contrary to previous studies, we conclude that the Plio-Pleistocene epoch seems to be not associated to further significant speciation events within Anatolian-Hyrcanian brown frogs.

Description of a new brown frog from Tsushima Island, Japan (Anura: Ranidae: Rana)

Because all available evidence from allozymes, mtDNA sequences, and artificial hybridization suggests presence of high genetic differentiation between populations of East Asian brown frogs currently assigned to Rana dybowskii Günther, 1876, I compared morphological characters between specimens from Tsushima Island of Japan and Maritime territory of Russia. The population from Tsushima is slightly, but significantly different from R. dybowskii from Russia, including the holotype. I therefore consider the Tsushima population to be specifically distinct, and describe it as a new species R. uenoi. The new species also occurs in the Korean Peninsula and adjacent islands, but the distributional relationships with R. dybowskii are unclear, as detailed distribution in northern Korea is lacking.

Identifying homomorphic sex chromosomes from wild-caught adults with limited genomic resources

We demonstrate a genotyping-by-sequencing approach to identify homomorphic sex chromosomes and their homolog in a distantly related reference genome, based on noninvasive sampling of wild-caught individuals, in the moor frog Rana arvalis. Double-digest RADseq libraries were generated using buccal swabs from 30 males and 21 females from the same population. Search for sex-limited markers from the unfiltered data set (411 446 RAD tags) was more successful than searches from a filtered data set (33 073 RAD tags) for markers showing sex differences in heterozygosity or in allele frequencies. Altogether, we obtained 292 putatively sex-linked RAD loci, 98% of which point to male heterogamety. We could map 15 of them to the Xenopus tropicalis genome, all but one on chromosome pair 1, which seems regularly co-opted for sex determination among amphibians. The most efficient mapping strategy was a three-step hierarchical approach, where R. arvalis reads were first mapped to a low-coverage genome of Rana temporaria (17 My divergence), then the R. temporaria scaffolds to the Nanorana parkeri genome (90 My divergence), and finally the N. parkeri scaffolds to the X. tropicalis genome (210 My). We validated our conclusions with PCR primers amplifying part of Dmrt1, a candidate sex determination gene mapping to chromosome 1: a sex-diagnostic allele was present in all 30 males but in none of the 21 females. Our approach is likely to be productive in many situations where biological samples and/or genomic resources are limited.

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.

A family of antimicrobial peptides related to japonicin-2 isolated from the skin of the chaochiao brown frog Rana chaochiaoensis

Four structurally-related peptides with antimicrobial activity were isolated from an extract of the skin of the Chinese brown frog, Rana chaochiaoensis Liu, 1946. Determination of their primary structures revealed that they are members of the japonicin-2 family, previously identified only in the skin of the Japanese brown frog, R. japonica. Japonicin-2CHa (FVLPLLGILPKELCIVLKKNC) represented the most abundant peptide in the extract but its growth-inhibitory potency against Staphylococcus aureus (minimum inhibitory concentration, MIC=100 microM) and Escherichia coli (MIC>100 microM) was appreciably less than that of the more cationic japonicin-2 (FGLPMLSILPKALCILLKRKC). The high degree of structural similarity of japonicin-2CHb (VVPAFVLLKKAICIMLKRNC) with japonicin-2CHc (K9 --> R), and japonicin-2CHd (L16 --> F) is suggestive of recent gene duplication events. The data indicate a close phylogenetic relationship between R. chaochiaoensis and R. japonica but demonstrate that the species are not conspecific.

Antimicrobial peptides from the skin of the Tsushima brown frog Rana tsushimensis

The Tsushima brown frog Rana tsushimensis Stejneger, 1907 exists in reproductive isolation on the island of Tsushima, Japan. Six peptides with antimicrobial activity were isolated in pure form from an extract of the skin of this species and their amino acid sequences identified them as members of the brevinin-1 (one peptide), brevinin-2 (one peptide) and temporin (four peptides) families. The C-terminally alpha-amidated brevinin-1 peptide (FLGSIVGALASALPSLISKIRN.NH2) lacks the cyclic heptapeptide domain Cys18-(Xaa)4-Lys-Cys24 at the COOH-terminus of the molecule that characterizes other members of that family. A structurally similar brevinin-1 peptide, also lacking the cyclic domain, was previously isolated from the skin of the Ryukyu brown frog Rana okinavana, indicative of a close phylogenetic relationship between the species. Brevinin-2TSa (GIMSLFKGVLKTAGKHVAGSLVDQLKCKITGGC) showed broad-spectrum growth inhibitory activity against a range of Gram-negative and Gram-positive bacteria (including methicillin-resistant Staphylococcus aureus) (minimum inhibitory concentrations< or =25 microM) and relatively low hemolytic activity against human erythrocytes (LD50=100 microM). The peptide therefore represents a candidate for drug development.

Multi-character perspectives on the evolution of intraspecific differentiation in a neotropical hylid frog

Background

Multi-character empirical studies are important contributions to our understanding of the process of speciation. The relatively conserved morphology of, and importance of the mate recognition system in anurans, combined with phylogenetic tools, provide an opportunity to address predictions about the relative role of each in the process of speciation. We examine the relationship among patterns of variation in morphology, call characters, and 16S gene sequences across seven populations of a neotropical hylid frog, Hyla leucophyllata, to infer their relative importance in predicting the early stages of population differentiation.

Results

Multivariate analyses demonstrate that both morphological and call characteristics were significantly variable among populations, characterized by significantly lower intra-population dispersion in call space than morphological space, and significantly greater among-population variation in call structure. We found lack of concordance between a 16S DNA phylogeny of Hyla leucophyllata and the significant population-level differentiation evident in both external morphology and male advertisement call. Comparisons of the reconstructed gene trees to simulated lineages support the notion that variation in call cannot be simply explained by population history.

Conclusion

Discordance among traits may reflect sampling biases (e.g. single genetic marker effects), or imply a decoupling of evolution of different suites of characters. Diagnostic differences among populations in call structure possibly reflect local selection pressures presented by different heterospecific calling assemblages and may serve as a precursor of species-wide differentiation. Differentiation among populations in morphology may be due to ecophenotypic variation or to diversifying selection on body size directly, or on frequency attributes of calls (mediated by female choice) that show a strong relationship to body size.

Climatic oscillations triggered post-Messinian speciation of Western Palearctic brown frogs (Amphibia, Ranidae)

Oscillating glacial cycles over the past 2.4 million years are proposed to have had a major impact on the diversity of contemporary species communities. We used mitochondrial and nuclear DNA sequence data to infer phylogenetic relationships within Western Palearctic brown frogs and to test the influence of Pliocene and Pleistocene climatic changes on their evolution. We sequenced 1976bp of the mitochondrial genes 16S rRNA and cytochrome b and of the nuclear rhodopsin gene for all current species and subspecies. Based on an established allozyme clock for Western Palearctic water frogs and substitution rate constancy among water frogs and brown frogs, we calibrated a molecular clock for 1425bp of the 16S and rhodopsin genes. We applied this clock to date speciation events among brown frogs. Western Palearctic brown frogs underwent a basal post-Messinian radiation about 4 million years ago (mya) into five major clades: three monotypic lineages (Rana dalmatina, Rana latastei, Rana graeca), an Anatolian lineage, and a lineage comprising Rana italica, Rana arvalis, and all Iberian taxa. Polytypic lineages radiated further in concordance with the onset of climatic oscillations ca. 3.2, 2.0, and 1.0-0.6 mya, respectively. The dated fossil record corroborates our paleobiogeographic scenario. We conclude that drastic climatic changes followed by successive temperature oscillations "trapped" most brown frog species in their southern European glacial refugia with enough time to speciate. Substantial dispersal was only possible during extensive interglacial periods of a constant subtropical climate.

Genetic relationships among Korean brown frog species (Anura, Ranidae), with special reference to evolutionary divergences between two allied species Rana dybowskii and R. huanrenensis

Allozyme analysis for 41 populations of brown frog species, Rana dybowskii, R. huanrenensis, and R. amurensis from Korea and three reference species (Chinese R. chensinensis and Japanese R. dybowskii and R. tsushimensis), were performed to clarify taxonomic status of Korean brown frogs. The level of average genetic differentiation (Nei's D) among local populations of each species in Korea was very low (D<0.01 2) and Korean and Japanese R. dybowskii also showed conspecific level of differentiation (D=0.070). Whereas, much larger, discrete genetic differences were detected in the interspecific comparisons (D>0.370). In the genetic relationships among five species examined, the 24 chromosome brown frogs (R. dybowskii, R. huanrenensis, and R. chensinensis) did not form a monophyletic group. Rana dybowskii with the chromosome number of 2n=24 was grouped together with R. amurensis with the chromosome number of 2n=26. The hypothesis of reversal change from 24 to 26 in Korean R. amurensis seems to better explain the phylogenetic relationships of east Asian brown frogs than the assumption of parallel reduction in chromosome number from 2n=26 to 24 in R. dybowskii and in the common ancestor of R. huanrenensis and R. chensinensis. The genetic, morphological, and reproductive divergences between Korean R. dybowskii and R. huanrenensis were compared.