Variation in Bufadienolide Composition of Parotoid Gland Secretion From Three Taxa of Japanese Toads

Toads of the genus Bufo synthesize and accumulate bufadienolides (BDs) in their parotoid glands. BDs are cardiotonic steroids that play an important role in defense against the toads' predators. Three bufonid taxa occur in mainland Japan, Bufo japonicus formosus, B. j. japonicus, and B. torrenticola. The chemical structures of BDs isolated from B. j. formosus were studied several decades ago, but there is no further information on the toxic components of Japanese toads and their metabolism. In this study, we analyzed BDs of toads from throughout Japan and compared the BD profiles by liquid chromatography/mass spectrometry (LC/MS) and hierarchical cluster analysis (HCA). We observed BDs in three taxa of Japanese toads, and identified five of the most common BDs by nuclear magnetic resonance (NMR) analyses. Of the five BDs, only bufalin was detected in all individuals. HCA of individual BD profiles divided the three taxa into five primary clusters and several subclusters. This result indicates that BD profiles differ both among and within the taxa. The clustering pattern of BDs is generally concordant with a phylogenetic tree reconstructed from the mitochondrial cytochrome b gene of Japanese toads. Our results suggest that the BDs of Japanese toads have diversified not in response to specific selective pressures, but simply due to population structuring over evolutionary time.

Phytogeographic patterns and cryptic diversity in an aposematic toad from NW Argentina

The Yungas Redbelly Toad, Melanophryniscus rubriventris, is patchily distributed in Argentina, confined to the upland portion (1000-2000m above sea level) of the montane forests of northern and central regions of Salta, and in central-eastern and south-eastern Jujuy. This species is known for its striking aposematic color variation across its geographic distribution, and was once treated as a complex of three subspecies based on distinctive color patterns. Here we assess the geographical genetic variation within M. rubriventris and quantify divergence in color and pattern among individuals sampled from Northwestern Argentina. We compare multi-gene phylogeography of M. rubriventris to patterns of dorsal and ventral coloration to test whether evolutionary affinities predict variation in warning color. Our results reveal two well-supported species lineages: one confined to the extreme northern portion of our sampling area, and the other extending over most of the Argentine portion of the species' range, within which there are two populations. However, these well-supported evolutionary relationships do not mirror the marked variation in warning coloration. This discordance between DNA genealogy and warning color variation may reflect selection brought about by differences in local predation pressures, potentially coupled with effects of sexual selection and thermoregulation.

Phylogenetic relationships of diurnal, phytotelm-breeding Melanophryniscus (Anura: Bufonidae) based on mitogenomic data

Melanophryniscus is a bufonid frog genus with a broad geographic distribution over southeastern South America. In recent years, several new species of Melanophryniscus have been discovered in southern Brazil showing a distinctive life-history strategy for the genus - breeding in phytotelmata - as well as a strong association with high-altitude regions. In this study, we use mitogenomic data to infer the phylogenetic relationships among diurnal, phytotelm-breeding Melanophryniscus and to determine the timing of their divergence. We obtained the mitochondrial genomes (not including the control region) for eight individuals of Melanophryniscus representing all three described species (M. alipioi, M. milanoi, and M. xanthostomus), as well as some recently-discovered and potentially new species. Gene order was conserved in all species and corresponded to the general order found in bufonids. Although the phylogenetic relationships among the studied species was poorly supported, dating confirmed that they diverged during the Pleistocene, suggesting that phytotelm breeding could have arisen during drier periods in the glacial/interglacial cycles due to a decrease in the availability of permanent streams or ephemeral/temporary streams or ponds in which Melanophryniscus species commonly breed.

Polyphyly of Asian Tree Toads, Genus Pedostibes Günther, 1876 (Anura: Bufonidae), and the Description of a New Genus from Southeast Asia

The Asian Tree Toad genus Pedostibes, as currently understood, exhibits a conspicuously disjunct distribution, posing several immediate questions relating to the biogeography and taxonomy of this poorly known group. The type species, P. tuberculosus and P. kempi, are known only from India, whereas P. hosii, P. rugosus, and P. everetti are restricted to Southeast Asia. Several studies have shown that these allopatric groups are polyphyletic, with the Indian Pedostibes embedded within a primarily South Asian clade of toads, containing the genera Adenomus, Xanthophryne, and Duttaphrynus. Southeast Asian Pedostibes on the other hand, are nested within a Southeast Asian clade, which is the sister lineage to the Southeast Asian river toad genus Phrynoidis. We demonstrate that Indian and Southeast Asian Pedostibes are not only allopatric and polyphyletic, but also exhibit significant differences in morphology and reproductive mode, indicating that the Southeast Asian species’ are not congeneric with the true Pedostibes of India. As a taxonomic solution, we describe a new genus, Rentapiagen. nov. to accommodate the Southeast Asian species.

Three New Species of Phytotelm-Breeding Melanophryniscus from the Atlantic Rainforest of Southern Brazil (Anura: Bufonidae)

Three new species of Melanophryniscus are described from the Serra do Mar mountain range of the state of Santa Catarina, southern Brazil. All species are found at intermediate to high altitudes and share phytotelm-breeding as their reproductive strategy. The new species are distinguished from other phytotelm-breeding Melanophryniscus based on different combinations of the following traits: snout-vent length, presence of white and/or yellow spots on forearms, mouth, belly and cloaca, pattern and arrangement of warts, and presence and number of corneous spines. The discovery of these species in a rather restricted geographical area suggests that the diversity of phytotelm-breeding species of Melanophryniscus might be severely underestimated. The conservation status of these species is of particular concern, given that one of them is at risk of extinction not only due to its restricted habitat, but also because of anthropogenic disturbances.

[Comparison of chemical composition between fresh and processed Bufonis Venenum by UPLC-TQ-MS]

Toad venom is the Bufo bufo gargarizans or B. melanostictus after the ears of the gland secretion, used in the treatment of various cancers in recent years. Research shows that the main anti-tumor components in bufadienolide. Bufadienolide have free type structure and conjunct type structure. To identify and clarify the difference between bufogenin and bufotoxin contained in Bufonis Venenum, which was from B. bufo gargarizans, an UPLC-TQ-MS method has been established. UPLC-TQ-MS method was used to identify and quantify the major bufadienolides in Bufonis Venenum. UPLC-TQ-MS assay with positive ion mode was performed on a Waters ACQUITY UPLC BEH C, (2.1 mm x 100 mm, 1.7 µm) with the mobile phase consisting of 0. 1% aqueous formic and acidacetonitrile in gradient elution at a flow rate of 0.4 mL · min⁻¹ and the column temperature was set at 35 °C. By comparing their retention time and high resolution mass data of Bufonis Venenum extracts, 37 effective components were primarily identified by MS/MS analysis in positive ion mode. Twenty-six of them were free-type bufadienolides (bufogenin), 11 of them were conjugated bufadienolides. There were significant differences in the main composition between fresh and processed Bufonis Venenum. The study found that the chemical composition of toad venom through great changes after processing, conjunct type content is much less, free type content as well change.

Origin and genome evolution of polyploid green toads in Central Asia: evidence from microsatellite markers

Polyploidization, which is expected to trigger major genomic reorganizations, occurs much less commonly in animals than in plants, possibly because of constraints imposed by sex-determination systems. We investigated the origins and consequences of allopolyploidization in Palearctic green toads (Bufo viridis subgroup) from Central Asia, with three ploidy levels and different modes of genome transmission (sexual versus clonal), to (i) establish a topology for the reticulate phylogeny in a species-rich radiation involving several closely related lineages and (ii) explore processes of genomic reorganization that may follow polyploidization. Sibship analyses based on 30 cross-amplifying microsatellite markers substantiated the maternal origins and revealed the paternal origins and relationships of subgenomes in allopolyploids. Analyses of the synteny of linkage groups identified three markers affected by translocation events, which occurred only within the paternally inherited subgenomes of allopolyploid toads and exclusively affected the linkage group that determines sex in several diploid species of the green toad radiation. Recombination rates did not differ between diploid and polyploid toad species, and were overall much reduced in males, independent of linkage group and ploidy levels. Clonally transmitted subgenomes in allotriploid toads provided support for strong genetic drift, presumably resulting from recombination arrest. The Palearctic green toad radiation seems to offer unique opportunities to investigate the consequences of polyploidization and clonal transmission on the dynamics of genomes in vertebrates.

[Bufo gargarizans mcl-1 cloning and its prokaryotic recombinant protein expression]

MCL-1 is encoded by myeloid cell leukemia-1 gene (mcl-1), which is one of the anti-apoptotic members of bcl-2 cell apoptotic gene superfamily. ChanSu is made of dorsal secretions of several Bufo species, commonly used in the prescriptions of traditional Chinese medicine for treating many diseases including cancer. To clarify if mcl-1 is expressed in the dorsal skin of B. gargarizans, the PCR (polymerase chain reaction) was performed with its dorsal skin first strand cDNA as the template and a pair of specific primers of mcl-1, and PCR products were cloned into the pGM-T vector. DNA sequencing indicated that the ORF length was 639 bp encoding 212 amino acid residues, and the homology of 44%-95% with the MCL-1 of several other animals. For the further studies on MCL-1 biological functions during the oncogenesis and preparation of its antibody, the prokaryotic expression construct of pET-28b-mcl-1 was prepared which was confirmed by DNA sequencing, and its recombinant protein expression (0.02% wet weight) in E. coli BL21 (DE3) strain was confirmed by SDS-PAGE and Western blotting.

Surviving chytridiomycosis: differential anti-Batrachochytrium dendrobatidis activity in bacterial isolates from three lowland species of Atelopus

In the Neotropics, almost every species of the stream-dwelling harlequin toads (genus Atelopus) have experienced catastrophic declines. The persistence of lowland species of Atelopus could be explained by the lower growth rate of Batrachochytrium dendrobatidis (Bd) at temperatures above 25 °C. We tested the complementary hypothesis that the toads' skin bacterial microbiota acts as a protective barrier against the pathogen, perhaps delaying or impeding the symptomatic phase of chytridiomycosis. We isolated 148 cultivable bacterial strains from three lowland Atelopus species and quantified the anti-Bd activity through antagonism assays. Twenty-six percent (38 strains representing 12 species) of the bacteria inhibited Bd growth and just two of them were shared among the toad species sampled in different localities. Interestingly, the strongest anti-Bd activity was measured in bacteria isolated from A. elegans, the only species that tested positive for the pathogen. The cutaneous bacterial microbiota is thus likely a fitness-enhancing trait that may (adaptation) or not (exaptation) have appeared because of natural selection mediated by chytridiomycosis. Our findings reveal bacterial strains for development of local probiotic treatments against chytridiomycosis and also shed light on the mechanisms behind the frog-bacteria-pathogen interaction.

Comparison of toad venoms from different Bufo species by HPLC and LC-DAD-MS/MS

ETHNOPHARMACOLOGICAL RELEVANCE

Toad venom, called Chansu in China, has been widely used for the treatment of heart failure, sores, pains, and various cancers for a long time in clinic.

AIM OF THE STUDY

The aim of the study is to investigate the chemical differences among a variety of toad venoms from different geographic locations and related Bufo species.

MATERIALS AND METHODS

Ten batches of commercial toad venom collected from different regions in China, one batch of fresh toad venom obtained from Bufo bufo gargarizans, and six batches of related Bufo species were analyzed by HPLC and LC-DAD-MS/MS. Individual components were identified by comparison of retention times, UV spectra, and mass spectra with authentic compounds, standard addition, as well as summarized MS fragmentation rules. Based on the profile of identified constituents and the content of cinobufagin and resibufogenin, the chemical differences observed among different samples are discussed.

RESULTS

Overall, 43 compounds were identified in the methanolic extracts of the different samples of toad venom. Besides of suberoyl arginine, several free bufadienolides, bufadienolide sulfates, and suberoyl esters of bufadienolides were found. The total amounts of cinobufagin and resibufogenin, which are the only two control markers according to the current Chinese Pharmacopoeia, varied widely from 0.7% to 10.9% in the commercial Chansu samples collected in the different locations in China. Low levels of resibufogenin, but no cinobufagin was observed in the samples from Bufo melanosticus and Bufo marinus, and even neither of both compounds was found in the sample from Bufo viridis.

CONCLUSIONS

The chemical profiles of the different commercial and collected toad venoms from related Bufo species differed significantly, not only in the absolute and relative contents, but also in the number and type of the constituents. The main reason for this variation are species-specific differences, but additional factors, such as the harvest and post-harvest processing, and adaption to environmental factors in different geographic locations, also seem to contribute.

Post-Messinian evolutionary relationships across the Sicilian channel: mitochondrial and nuclear markers link a new green toad from Sicily to African relatives

BACKGROUND

Little attention has been paid to the consequences of the last landbridge between Africa and Sicily on Mediterranean biogeography. Previous paleontological and scarce molecular data suggest possible faunal exchange later than the well-documented landbridge in the Messinian (5.3 My); however, a possible African origin of recent terrestrial Sicilian fauna has not been thoroughly tested with molecular methods. To gain insight into the phylogeography of the region, we examine two mitochondrial and two nuclear markers (one is a newly adapted intron marker) in green toads (Bufo viridis subgroup) across that sea barrier, the Strait of Sicily.

RESULTS

Extensive sampling throughout the western Mediterranean and North Africa revealed a deep sister relationship between Sicilian (Bufo siculus n.sp.) and African green toads (B. boulengeri) on the mitochondrial and nuclear level. Divergence times estimated under a Bayesian-coalescence framework (mtDNA control region and 16S rRNA) range from the Middle Pliocene (3.6 My) to Pleistocene (0.16 My) with an average (1.83 to 2.0 My) around the Pliocene/Pleistocene boundary, suggesting possible land connections younger than the Messinian (5.3 My). We describe green toads from Sicily and some surrounding islands as a new endemic species (Bufo siculus). Bufo balearicus occurs on some western Mediterranean islands (Corsica, Sardinia, Mallorca, and Menorca) and the Apennine Peninsula, and is well differentiated on the mitochondrial and nuclear level from B. siculus as well as from B. viridis (Laurenti), whose haplotype group reaches northeastern Italy, north of the Po River. Detection of Calabrian B. balearicus haplotypes in northeastern Sicily suggests recent invasion. Our data agree with paleogeographic and fossil data, which suggest long Plio-Pleistocene isolation of Sicily and episodic Pleistocene faunal exchange across the Strait of Messina. It remains unknown whether both species (B. balearicus, B. siculus) occur in sympatry in northern Sicily.

CONCLUSION

Our findings on green toads give the first combined mitochondrial and nuclear sequence evidence for a phylogeographic connection across the Strait of Sicily in terrestrial vertebrates. These relationships may have implications for comparative phylogeographic research on other terrestrial animals co-occurring in North Africa and Sicily.

Population genetics and phylogeography of Bufo gargarizans in China

The Chinese toad Bufo gargarizans has an extensive range, covering nearly all of China. This study addresses the population genetic structure of this toad and the relation to geographic location. Partial sequences of the mtDNA control region were obtained from 14 populations, 29 haplotypes were defined, and 73 variable sites were found, including 66 parsimony informative sites. All population genetic analyses indicated no clear geographic pattern in the distribution of haplotypes. The genetic divergence between the populations was significant. Phylogeographic analyses suggested that past fragmentation and/or long-distance colonization seemed to have shaped the present-day distribution of the haplotypes of B. gargarizans.

Tests of phenotypic and genetic concordance and their application to the conservation of Panamanian golden frogs (Anura, Bufonidae)

Evolutionarily significant units (ESUs) differ in the extent to which they capture, or even consider, adaptive variation, and most such designations are based solely on neutral genetic differences that may not capture variation relevant to species' adaptabilities to changing environmental conditions. While concordant patterns of divergence among data sets (i.e. neutral and potentially non-neutral characters) can strengthen ESU designations, determining whether such criteria are met for highly variable taxa is especially challenging. This study tests whether previously defined ESUs for endangered Panamanian golden frogs (Atelopus varius and Atelopus zeteki) exhibit concordant variation among multiple phenotypic traits and mitochondrial DNA sequences, and the extent to which such divergence corresponds to environmental differences. Multivariate analyses identify phenotypic and genetic differentiation consistent with proposed ESUs and support the status of A. varius and A. zeteki as separate species. Moreover, the significant association detected between ESU co-membership and genetic similarity, which remained strong after removing the effect of geographic distance, also indicates that genetic differences are not simply due to isolation by distance. Two phenotypic characters (body size and the extent of dorsal black patterning) that differ among ESUs also co-vary with environmental differences, suggesting that to the extent that these phenotypic differences are heritable, variation may be associated with adaptive divergence. Lastly, discriminant function analyses show that the frogs can be correctly assigned to ESUs based on simultaneous analysis of multiple characters. The study confirms the merit of conserving the previously proposed golden frog ESUs as well as demonstrates the utility and feasibility of combined analyses of ecological, morphological and genetic variation in evaluating ESUs, especially for highly variable taxa.

Evolution of mitochondrial relationships and biogeography of Palearctic green toads (Bufo viridis subgroup) with insights in their genomic plasticity

Taxa involving three bisexually reproducing ploidy levels make green toads a unique amphibian system. We put a cytogenetic dataset from Central Asia in a molecular framework and apply phylogenetic and demographic methods to data from the entire Palearctic range. We study the mitochondrial relationships of diploids to infer their phylogeography and the maternal ancestry of polyploids. Control regions (and tRNAs between ND1 and ND2 in representatives) characterize a deeply branched assemblage of twelve haplotype groups, diverged since the Lower Miocene. Polyploidy has evolved several times: Central Asian tetraploids (B. oblongus, B. pewzowi) have at least two maternal origins. Intriguingly, the mitochondrial ancestor of morphologically distinctive, sexually reproducing triploid taxa (B. pseudoraddei) from Karakoram and Hindukush represents a different lineage. We report another potential case of bisexual triploid toads (B. zugmayeri). Identical d-loops in diploids and tetraploids from Iran and Turkmenistan, which differ in morphology, karyotypes and calls, suggest multiple origins and retained polymorphism and/or hybridization. A similar system involves diploids, triploids and tetraploids from Kyrgyzstan and Kazakhstan where green toads exemplify vertebrate genomic plasticity. A new form from Sicily and its African sister species (B. boulengeri) allow internal calibration and divergence time estimates for major clades. The subgroup may have originated in Eurasia rather than Africa since the earliest diverged lineages (B. latastii, B. surdus) and earliest fossils occur in Asia. We delineate ranges, contact and hybrid zones. Phylogeography, including one of the first non-avian datasets from Central Asian high mountains, reflects Quaternary climate and glaciation.

Developmental plasticity mirrors differences among taxa in spadefoot toads linking plasticity and diversity

Developmental plasticity is found in most organisms, but its role in evolution remains controversial. Environmentally induced phenotypic differences may be translated into adaptive divergence among lineages experiencing different environmental conditions through genetic accommodation. To examine this evolutionary mechanism, we studied the relationship between plasticity in larval development, postmetamorphic morphology, and morphological diversity in spadefoot toads, a group of closely related species that are highly divergent in the larval period and body shape and are distributed throughout temperate areas of both the New and the Old World. Previous studies showed that accelerated metamorphosis is adaptive for desert-dwelling spadefoot toads. We show that even under common garden conditions, spadefoot toad species show divergent reaction norms for the larval period. In addition, experimentally induced changes in the larval period caused correlated morphological changes in postmetamorphic individuals such that long larval periods resulted in relatively longer hindlimbs and snouts. A comparative analysis of morphological variation across spadefoot toad species also revealed a positive correlation between the larval period and limb and snout lengths, mirroring the effects of within-species plasticity at a higher taxonomic level. Indeed, after approximately 110 Ma of independent evolution, differences in the larval period explain 57% of the variance in relative limb length and 33% of snout length across species. Thus, morphological diversity across these species appears to have evolved as a correlated response to selection for a reduced larval period in desert-dwelling species, possibly diverging from ancestral plasticity through genetic accommodation.

Gene flow and functional connectivity in the natterjack toad

Functional connectivity is a key factor for the persistence of many specialist species in fragmented landscapes. However, connectivity estimates have rarely been validated by the observation of dispersal movements. In this study, we estimated functional connectivity of a real landscape by modelling dispersal for the endangered natterjack toad (Bufo calamita) using cost distance. Cost distance allows the evaluation of 'effective distances', which are distances corrected for the costs involved in moving between habitat patches in spatially explicit landscapes. We parameterized cost-distance models using the results of our previous experimental investigation of natterjack's movement behaviour. These model predictions (connectivity estimates from the GIS study) were then confronted to genetic-based dispersal rates between natterjack populations in the same landscape using Mantel tests. Dispersal rates between the populations were inferred from variation at six microsatellite loci. Based on these results, we conclude that matrix structure has a strong effect on dispersal rates. Moreover, we found that cost distances generated by habitat preferences explained dispersal rates better than did the Euclidian distances, or the connectivity estimate based on patch-specific resistances (patch viscosity). This study is a clear example of how landscape genetics can validate operational functional connectivity estimates.

Lusitania revisited: A phylogeographic analysis of the natterjack toad Bufo calamita across its entire biogeographical range

Attempts to understand the current distributions of plants and animals require both historical and ecological information. Phylogeography has proved highly effective in elucidating historical events such as postglacial colonisations in north temperate zones. However, interesting questions still await resolution. Lusitanian distributions of fauna and flora in western Europe, for example, have puzzled biogeographers for more than 150 years. Lusitanian species have highly disjunct distributions in Ireland and in Iberia, often with few or no other populations inbetween. Despite much debate, no agreed explanation for Lusitanian distributions has yet emerged. We investigated the phylogeographic structure of one Lusitanian species, the natterjack toad Bufo calamita, using mitochondrial DNA control region sequences and allelic variation at eight microsatellite loci. Our results show that this amphibian must have survived in north European refugia, as well as in Iberia, during and since the last (Weichselian) glacial maximum around 20,000 years before present (BP). Subsequent local recolonisation after the Younger Dryas cooling around 11,000 years BP best explains the Lusitanian aspect of natterjack toad distribution.

Phylogeographic analysis of the Bufo gargarizans species complex: A revisit

Using mtDNA sequencing and allozyme electrophoresis data, we tested the "vicariance followed by dispersal" hypothesis of the Bufo gargarizans species group and re-evaluated the species status in the general lineages species concept. A phylogenetic analysis suggested that dispersal, instead of vicariance, dominated the history of the species group. There was a general trend of west to east dispersal, while some lineages from the east subsequently returned to the west. The secondary admixture of those previously allopatric lineages produced substantial levels of sympatric genetic diversity, often as high as 7.0% pairwise difference within populations. The phylogenetic hypothesis does not support the current two species designation. Neither B. andrewsi nor B. gargarizans represents an independent evolutionary lineage, and monophyletic groups did not correspond to geographically discrete groups. Allozyme data also failed to reveal any fixed allelic difference among the populations. Therefore, we recommend regarding the complex as a single species, Bufo gargarizans, without subspecies division.

Multiple origins of tetraploid taxa in the Eurasian Bufo viridis subgroup

We used Q-banding and analyzed nucleolar organizing regions (NORs) to study the cytogenetic evolution of tetraploids within the Palearctic Bufo viridis subgroup, the only known amphibian complex comprising di-, tri- and tetraploid bisexually reproducing taxa. We examined three diploid (2n) nominal taxa (Bufo viridis viridis, B. v. turanensis, B. v. kermanensis) from five Eurasian localities and six tetraploid (4n) nominal taxa (B. oblongus, B. o. danatensis, B. pewzowi pewzowi, B. p. taxkorensis, B. p. unicolor, B. p. strauchi) from eight Central Asian localities. Homeologous chromosomes of 2n and 4n toads exhibit a similar morphology. Silver-staining and in situ hybridization revealed terminal NORs in the long arms of chromosomes 6 in all 2n but in only two out of four chromosomes 6 in all 4n taxa. Q-banding and a rapidly evolving mitochondrial marker suggest at least two origination events for Asian 4n toads: "Western Central Asian tetraploids" (B. oblongus Nikolsky, 1896) exhibit distinct differences within some chromosome quartets, which are divisible into pairs of chromosomes and may be allopolyploid. In contrast, "Central Asian tetraploids" (B. pewzowi Bedriaga, 1898) showed homogenous Q-banding patterns within each quartet, suggesting autopolyploidy. In Northeastern Iran, we discovered a zone of either common ancestry or hybridization of 2n and Western Central Asian 4n toads. This raises intriguing questions about how diploid and tetraploid taxa may evolve by exchanging genetic material.

Contrasting population structures in two sympatric anurans: implications for species conservation

A general prediction of the neutral theory of evolution is that genetic diversity should correlate positively with effective population size. We show here that diversity across eight microsatellite loci was consistently and substantially lower in one common amphibian (Bufo bufo) than in another with similar life history traits (Rana temporaria) despite B. bufo having the larger breeding assemblage sizes. However, B. bufo breeding assemblages were much more highly differentiated than those of R. temporaria according to both Fst and Rst estimators. These differences occurred in shared habitats across identical geographical distances. The patterns of genetic diversity and differentiation detected in these two species were probably a consequence of high gene flow in R. temporaria but much lower gene flow among the larger but more dispersed B. bufo assemblages. These observations highlight the difficulty of defining the boundaries of wild populations, and show how two broadly similar species can exhibit very different population dynamics.

Multiple, recurring origins of aposematism and diet specialization in poison frogs

Aposematism is the association, in a prey organism, of the presence of a warning signal with unprofitability to predators. The origin of aposematism is puzzling, because of its predicted low probability of establishment in a population due to the prey's increased conspicuousness. Aposematism is a widespread trait in invertebrate taxa, but, in vertebrates, it is mostly evident in amphibians, reptiles, and fishes. Poison frogs (Dendrobatidae) are one of the most well known examples of the co-occurrence of warning coloration and toxicity. This monophyletic group of mostly diurnal leaf-litter Neotropical anurans has both toxic/colorful and palatable/cryptic species. Previous studies suggested a single origin of toxicity and warning coloration, dividing the family in two discrete groups of primitively cryptic and more derived aposematic frogs. Recent molecular phylogenetic analyses using mostly aposematic taxa supported this conclusion and proposed a single tandem origin of toxicity and conspicuous warning coloration. By using expanded taxon and character sampling, we reexamined the phylogenetic correlation between the origins of toxicity and warning coloration. At least four or five independent origins of aposematism have occurred within poison frogs; by using simulations, we rejected hypotheses of one, two, or three origins of aposematism (P < 0.002). We also found that diet specialization is linked with the evolution of aposematism. Specialization on prey, such as ants and termites, may have evolved independently at least two times.

Sexual size dimorphism in anurans

Several hypotheses have been proposed to explain the direction and extent of sexual size dimorphism in anurans (in which males are usually smaller than females) as a result of sexual selection. Here, we present an analysis to test the hypothesis that sexual dimorphism in anurans is largely a function of differences between the sexes in life-history strategies. Morphological and demographic data for anurans were collected from the literature, and the mean size and age in each sex were calculated for 51 populations, across 30 species and eight genera. Comparisons across 14 Rana species, eight Bufo species and across the genera showed a highly significant relationship between size dimorphism, measured using the female-male size ratio, and mean female-male age difference. A comparison of a subset of 17 of these species for which phylogenetic information was available, using the method of independent contrasts, yielded a similar result. These results indicate that most of the variation in size dimorphism in the anura can be explained in terms of differences in the age structure between the sexes in breeding populations. If sexual selection has an effect on size dimorphism in anurans, it is likely to be only a secondary one.

Molecular systematics, hybridization, and phylogeography of the Bufo americanus complex in Eastern North America

We reconstruct phylogenetic relationships among a well-studied group of toads and find relationships that differ greatly from the current taxonomic understanding. We use mitochondrial sequences encoding ND1, tRNA(Leu(UUR)), and part of 16S to infer relationships among members of the Bufo americanus complex. Focusing on the four taxa that historically have been most problematic due to morphological similarity and hybridization in sympatry, we sample 150 individuals from multiple populations across each species' geographic range. Our evidence conflicts with previous taxonomic hypotheses that were based on ability to hybridize, geographic distribution, and call variation. First, sequences from B. fowleri do not comprise the sister clade to sequences of B. woodhousii; therefore the previous classifications of B. fowleri as sister species to, or eastern subspecies of, B. woodhousii are both called into question. Second, sequences from B. americanus are more closely related to those of B. woodhousii than to those of B. terrestris, indicating that similar advertisement call characteristics evolved independently. Third, sequences of B. fowleri are paraphyletic, with sequences of B. terrestris embedded within. Lastly, sequences from B. fowleri cluster into three distinct mitochondrial clades, with some divergences corresponding to greater than 2mya. These clades are somewhat geographically structured, suggesting divergence in allopatry during the Pleistocene. These mitochondrial divergences are not accompanied by known phenotypic differences, however, suggesting either evolutionary stasis in morphology and behavior, cryptic phenotypic evolution, or that hybridization in secondary contact has homogenized phenotypic differences that may have arisen in allopatry.

The effect of miniaturized body size on skeletal morphology in frogs

Miniaturization has evolved numerous times and reached impressive extremes in the Anura. I compared the skeletons of miniature frog species to those of closely related larger species to assess patterns of morphological change, sampling 129 species from 12 families. Two types of morphological data were examined: (1) qualitative data on bone presence and absence; and (2) thin-plate spline morphometric descriptions of skull structure and bone shape. Phylogenetic comparative methods were used to address the shared history of species. Miniature anurans were more likely to lose skull bones and phalangeal elements of the limbs. Their skulls also showed consistent differences compared to those of their larger relatives, including relatively larger braincases and sensory capsules, verticalization of lateral elements, rostral displacement of the jaw joint, and reduction of some skull elements. These features are explained by functional constraints and by paedomorphosis. Variation among lineages in the morphological response to miniaturization was also explored. Certain lineages appear to be unusually resistant to the morphological trends that characterize miniature frogs as a whole. This study represents the first large-scale examination of morphology and miniaturization across a major, diverse group of organisms conducted in a phylogenetic framework and with statistical rigor.

A bisexually reproducing all-triploid vertebrate

Green toads are common in the Palaearctic region, where they have differentiated into several taxa. The toads exist with variable amounts of ploidy, similar to other anuran species or reptiles. In vertebrate biology, the very rare occurrence of triploidy is coupled with infertility or unisexuality, or requires the coexistence of individuals of different ploidy in a reproductive community. The reproduction of naturally occurring triploids has been reported to occur only through parthenogenesis, gynogenesis or hybridogenesis. The bisexual reproduction of pure triploids has been considered to be impossible because of the problem of equally distributing three chromosome sets in meiosis. Here we report geographically isolated populations of green toads (Bufo viridis complex) that are all-triploid and reproduce bisexually.

Molecular phylogeny and biogeography of West Indian toads (Anura: Bufonidae)

Toads of the Bufo peltocephalus Group (Anura: Bufonidae) occur throughout the Greater Antilles (Cuba to the Virgin Islands), a geographic region of relatively high endemicity. Previous morphological and immunological studies suggested that the West Indian toads are a monophyletic lineage derived from Neotropical Bufo but were unable to clarify relationships within the group. We examined the evolutionary relationships and biogeography of this group of frogs by collecting approximately 2 kb of mitochondrial DNA sequence data from eight West Indian species and selected non-West Indian species from the New World and the Old World. Our analyses support the monophyly of native West Indian toads and a New World origin for the group. Relationships among the West Indian species are less certain, but a Cuban lineage is defined in most analyses.

Cytogenetics of two species of Paratelmatobius (Anura: Leptodactylidae), with phylogenetic comments

In this paper we provide a cytogenetic analysis of Paratelmatobius cardosoi and Paratelmatobius poecilogaster. The karyotypes of both species showed a diploid number of 24 chromosomes and shared some similarity in the morphology of some pairs. On the other hand, pairs 4 and 6 widely differed between these complements. These karyotypes also differed in their NOR number and location. Size heteromorphism was seen in all NOR-bearing chromosomes of the two karyotypes. In addition, both karyotypes showed small centromeric C-bands and a conspicuous heterochromatic band in the short arm of chromosome 1, although with a different size in each species. The P. cardosoi complement also showed other strongly stained non-centromeric C-bands, with no counterparts in the P. cardosoi karyotype. Chromosome staining with fluorochromes revealed heterogeneity in the base composition of two of the non-centromeric C-bands of P. cardosoi. Comparison of the chromosomal morphology of these Paratelmatobius karyotypes with that of P. lutzii showed that the P. poecilogaster karyotype is more similar to that of P. lutzii than P. cardosoi. These cytogenetic results agree with the proposed species arrangements in the P. cardosoi and P. lutzii groups based on morphological and ecological data.

Phylogeography and speciation of the morphologically variable, widespread species Bufo valliceps, based on molecular evidence from mtDNA

The widespread, lowland toad Bufo valliceps has an unusual distribution in North and Middle America that straddles two major biogeographic areas; previous morphological studies of this species suggested the existence of two species. We used mitochondrial DNA sequences to examine the phylogeography of this species and discovered the existence of two distinct clades. We recognize these as two species: B. valliceps and B. nebulifer. These molecular data support morphological data from previous studies. Our results show low levels of molecular variation in a morphologically uniform temperate species (B. nebulifer) and high levels of molecular variation in a morphologically variable tropical species (B. valliceps), providing an example of molecules matching morphology. Two biogeographic hypotheses are tested to explain the current distribution of these species, based on a calibrated rate of evolution and the percent sequence divergence between the two species. A more recent Pleistocene dispersal event, followed by vicariance associated with rising sea level, is rejected in favor of an earlier Miocene-Pliocene vicariant hypothesis associated with the formation of the Trans-Mexican Neovolcanic Belt.

Phylogeny of east asian bufonids inferred from mitochondrial DNA sequences (Anura: amphibia)

We investigated the relationships of Asian bufonids using partial sequences of mitochondrial DNA genes. Twenty-six samples representing 14 species of Bufo from China and Vietnam and 2 species of Torrentophryne from China were examined. Three samples of Bufo viridis from Armenia and Georgia were also sequenced to make a comparison to its sibling tetraploid species B. danatensis. Bufo americanus, from Canada, was used as the outgroup. Sequences from the 12S ribosomal RNA, 16S ribosomal RNA, cytochrome b, and the control region were analyzed using parsimony. East Asian bufonids were grouped into two major clades. One clade included B. andrewsi, B. bankorensis, B. gargarizans, B. tibetanus, B. tuberculatus, its sister clade B. cryptotympanicus, and the 2 species of Torrentophryne. The second clade consisted of B. galeatus, B. himalayanus, B. melanostictus, and a new species from Vietnam. The placement of three taxa (B. raddei, B. viridis, and its sister species, B. danatensis) was problematic. The genus Torrentophryne should be synonymized with Bufo to remove paraphyly. Because B. raddei does not belong to the clade that includes B. viridis and B. danatensis, it was removed from the viridis species group. The species status of B. bankorensis from Taiwan is evaluated.

Phylogeny of frogs of the Physalaemus pustulosus species group, with an examination of data incongruence

Characters derived from advertisement calls, morphology, allozymes, and the sequences of the small subunit of the mitochondrial ribosomal gene (12S) and the cytochrome oxidase I (COI) mitochondrial gene were used to estimate the phylogeny of frogs of the Physalaemus pustulosus group (Leptodactylidae). The combinability of these data partitions was assessed in several ways: measures of phylogenetic signal, character support for trees, congruence of tree topologies, compatibility of data partitions with suboptimal trees, and homogeneity of data partitions. Combined parsimony analysis of all data equally weighted yielded the same tree as the 12S partition analyzed under parsimony and maximum likelihood. The COI, allozyme, and morphology partitions were generally congruent and compatible with the tree derived from combined data. The call data were significantly different from all other partitions, whether considered in terms of tree topology alone, partition homogeneity, or compatibility of data with trees derived from other partitions. The lack of effect of the call data on the topology of the combined tree is probably due to the small number of call characters. The general incongruence of the call data with other data partitions is consistent with the idea that the advertisement calls of this group of frogs are under strong sexual selection.

Phylogenetic relationships among bufonoid frogs (Anura:Neobatrachia) inferred from mitochondrial DNA sequences

Nucleotide sequences of portions of the mitochondrial 12S and 16S ribosomal RNA genes were used to extend a recent study of anuran phylogeny (Hay et al., Mol. Biol. Evol. 12: 928-937, 1995) and to further evaluate phylogenetic relationships within the Neobatrachia. An analysis of almost 900 nucleotides from each of 8 new representatives of the Dendrobatidae, Hylidae, Leptodactyolidae, and Myobatrachidae, plus 14 available members of the Neobatrachia provides support for 2 major lineages (Bufonoidea and Ranoidea) within this anuran suborder. The neotropical Bufonoidea and their derivatives are monophyletic. There is an interesting association of the 2 Australian myobatrachids with the South African Heleophrynidae, and the Sooglossidae is one of the basal bufonoid lineages. Within the New World bufonoid frogs, a monophyletic Dendrobatidae is strongly supported. An Australian hylid (Pelodryadinae) shows close affinity with the South American hylid Phyllomedusinae. A group composed of Hylinae (Hyla and Smilisca), Centrolenidae, Bufonidae, and the hylid Hemiphractinae, with the latter two clustered, was supported significantly. The addition of new taxa has more clearly defined some relationships within the suborder Neobatrachia and has indicated that the families Hylidae, Leptodactylidae, and Myobatrachidae may not be monophyletic.

The phylogenetic utility of cytochrome b: lessons from bufonid frogs

The mitochondrial cytochrome b gene is widely used in systematic studies to resolve divergences of many different ages. To investigate phylogenetic relationships among frogs of the large family Bufonidae, and to explore the utility of cytochrome b for this purpose, approximately one-third of the gene was sequenced from representatives of this group. Samples were chosen to represent a range of divergence levels within Bufonidae: (1) deep (= old), among species from around the world; (2) middle, among North American species; and (3) shallow (= young), within a single species group (the Bufo boreas group). The inferred amino acid sequences of cytochrome b are highly similar in these frogs although most pairwise comparisons of the nucleotide sequences are 15-20% different. Consequently insufficient information is available to generate robust phylogenetic hypotheses for the older divergences; silent differences are saturated and yet almost no informative replacement differences exist. Among the younger divergences, silent differences are not saturated and some resolution is possible. These results show that (1) the amino acid sequence of cytochrome b evolves differently in Bufonidae than expected based on other vertebrates (2) it consequently provides surprisingly little information about old divergences in Bufonidae, and (3) phylogenetic studies applying particular genes to new groups should begin with preliminary surveys of exemplar taxa representing the range of divergence times within the group to estimate the likely phylogenetic utility of that gene in that group.