How complex is the Bufo bufo species group?

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Differentiation of skull morphology and cranial kinesis in common toads

We examined the cranial morphology and cranial kinesis of the common toads Bufo bufo and B. spinosus with micro-computed tomography and geometric morphometrics and compared the results with published data for related species in a phylogenetic context. The species significantly diverge in skull shape. The skull of B. spinosus is shorter and higher, with a ventral arm of the squamosal bone and the jaw articulation point positioned perpendicular to the braincase, in comparison with a more lateral position in B. bufo. In either species, females have a shorter snout and a higher and wider skull at the jaw articulation point that is positioned more posteriorly, in comparison with conspecific males. High variation in the amount of bone ossification was recorded in both species, ranging from scarcely ossified and loosely connected bones to highly ossified and firmly connected bones. We also found that skull shape and inferred kinetic properties of the skull are highly variable across the Bufonini tribe. However, sample sizes are mostly small and intraspecific variation is high, which might compromise the analyses. Overall, the results suggest that developmental plasticity produces high variation in ossification and cranial kinesis, affecting individuals’ feeding performances. At the population level, this variation supports an efficient exploitation of the habitat and may promote morphological adaptation in a changing environment.

Cytogenetic Analysis in the Toad Species Bufo spinosus, Bufotes viridis and Epidalea calamita (Anura, Bufonidae) from the Mediterranean Area

Taxonomy in Bufonidae witnessed notable transformations. Bufotes viridis and Epidalea calamita, previously included in genus Bufo, were relocated in other genera, while the genus Bufo was restricted to members of the earlier Bufo bufo group. On the other hand, Bufo bufo sensu lato now includes four species: Bufo bufo, Bufo spinosus, Bufo verrucosissimus and Bufo eichwaldi. In this study, we examined three species of three Bufonidae genera (B. spinosus, B. viridis and E. calamita) by conventional (C-banding and Ag-NOR staining) and molecular (in situ hybridization with probes for telomeric repeats and rDNA loci, and genomic in situ hybridization (GISH)) cytogenetic methods. C-banding patterns are reported for the first time for B. spinosus and E. calamita populations from Iberian Peninsula and for B. viridis from Greece, and reveal several differences with the reported C-banded karyotypes described for other European populations of these species. Silver staining shows size heteromorphisms of the signals at the Nucleolar Organizing Region (NOR). By contrast, FISH with ribosomal probes only reveal size heteromorphism of rDNA sequences in E. calamita, suggesting that the differences observed after silver staining in B. spinosus and B. viridis should be attributed to differences in chromosomal condensation and/or gene activity rather than to differences in the copy number for ribosomal genes. Regarding telomeric repeats, E. calamita is the only species with interstitial telomeric sequences (ITS) located on centromeric regions, probably originated by accumulation of telomeric sequences in the centromeric heterochromatin. Finally, we analyzed the composition and distribution of repetitive sequences by genome in situ hybridization. These experiments reveal the accumulation of repetitive sequences in centromeric regions of the three species, although these sequences are not conserved when species from different genera are compared.

Identification and characterization of a new family of long satellite DNA, specific of true toads (Anura, Amphibia, Bufonidae)

Amphibians have some of the most variable genome sizes among vertebrates. Genome size variation has been attributed to repetitive and noncoding DNA, including satellite repeats, transposable elements, introns, and nuclear insertions of viral and organelle DNA. In vertebrates, satellite DNAs have been widely described in mammals, but few molecular studies have been carried out in amphibians. Here, we provide a detailed characterization of a new family of satellite DNA, present in all 15 examined species of the family Bufonidae. Southern-blot analysis and PCR reveal that this satellite is formed by monomers of 807 bp, is organized in tandem arrays, and has an AT-content of 57.4%. Phylogenetic analyses show that most clades exhibit species-specific variances, indicating that this satellite DNA has evolved by concerted evolution. The homogenization/fixation process is heterogeneous in Bufonidae, where the genera Bufo and Bufotes do not show species-specific differences, while populations from Rhinella marina exhibit population-specific changes. Additionally, variants of this satellite DNA have been identified in Duttaphrynus melanostictus and R. marina, supporting the 'library hypothesis' (a set, 'library', of satellite DNAs is shared by a species group). Physical mapping in Bufo bufo, Bufo spinosus, Epidalea calamita and Bufotes viridis provides evidence that this repetitive DNA is not dispersed in the karyotype, but accumulated in pericentromeric regions of some chromosomal pairs. This location, together with its presence in the transcriptomes of bufonids, could indicate a role in centromere function or heterochromatin formation and maintenance.

Phylogeography of Rhabdias spp. (Nematoda: Rhabdiasidae) collected from Bufo species in Honshu, Shikoku, and Kyushu, Japan including possible cryptic species

The genus Rhabdias Stiles & Hassall, 1905 comprises lung parasites of amphibians and reptiles worldwide. In Japan, 9 species have been recorded, including Rhabdias incerta Wilkie, 1930 which has been reported only in Bufo species. In this study, to assess the diversity of R. incerta, we performed molecular analyses of Rhabdias species sampled from three species/subspecies of Japanese toads namely Bufo japonicus, B. japonicus formosus, and B. torrenticola, collected in various regions of Honshu, Shikoku, and Kyushu, Japan. DNA sequence divergence was compared using mtDNA (COI) and nuclear DNA (28S) to identify possible cryptic species. Morphological analysis was performed through light microscopy and scanning electron microscopy (SEM). The results revealed that Bufo spp. serve as hosts for at least three Rhabdias species. Morphologically, most samples were identified as R. incerta but with a longer body and larger buccal cavity than originally described. Rhabdias incerta appears to be specific to the genus Bufo and is further subdivided into two or three phylogroups based on subspecies divisions and biogeography of their host. Some Rhabdias specimens collected in this study resemble R. tokyoensis Wilkie, 1930, parasitic in hosts from the order Caudata, which suggests host switching. Both molecular and morphological analyses suggested the presence of undescribed and cryptic Rhabdias species within toads collected in Japan. This study was the first to molecularly characterize Rhabdias species in Japan, including novel sequences of R. incerta and two undescribed species.

Asymmetrical sex reversal: Does the type of heterogamety predict propensity for sex reversal?

Sex reversal, a mismatch between phenotypic and genetic sex, can be induced by chemical and thermal insults in ectotherms. Therefore, climate change and environmental pollution may increase sex-reversal frequency in wild populations, with wide-ranging implications for sex ratios, population dynamics, and the evolution of sex determination. We propose that reconsidering the half-century old theory "Witschi's rule" should facilitate understanding the differences between species in sex-reversal propensity and thereby predicting their vulnerability to anthropogenic environmental change. The idea is that sex reversal should be asymmetrical: more likely to occur in the homogametic sex, assuming that sex-reversed heterogametic individuals would produce new genotypes with reduced fitness. A review of the existing evidence shows that while sex reversal can be induced in both homogametic and heterogametic individuals, the latter seem to require stronger stimuli in several cases. We provide guidelines for future studies on sex reversal to facilitate data comparability and reliability.

From Gondwana to the Yellow Sea, evolutionary diversifications of true toads Bufo sp. in the Eastern Palearctic and a revisit of species boundaries for Asian lineages

Taxa with vast distribution ranges often display unresolved phylogeographic structures and unclear taxonomic boundaries resulting in hidden diversity. This hypothesis-driven study reveals the evolutionary history of Bufonidae, covering the phylogeographic patterns found in Holarctic bufonids from the West Gondwana to the phylogenetic taxonomy of Asiatic true toads in the Eastern Palearctic. We used an integrative approach relying on fossilized birth-death calibrations, population dynamics, gene-flow, species distribution, and species delimitation modeling to resolve the biogeography of the clade and highlight cryptic lineages. We verified the near-simultaneous Miocene radiations within Western and Eastern Palearctic Bufo, c. 14.49-10.00 Mya, temporally matching with the maximum dust outflows in Central Asian deserts. Contrary to earlier studies, we demonstrated that the combined impacts of long dispersal and ice-age refugia equally contributed to the current genetic structure of Bufo in East Asia. Our findings reveal a climate-driven adaptation in septentrional Eastern Asian Bufo, explaining its range shifts toward northern latitudes. We resolve species boundaries within the Eastern Palearctic Bufo, and redefine the taxonomic and conservation units of the northeastern species: B. sachalinensis and its subspecies.

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene-environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

Hybridization and introgression between toads with different sex chromosome systems

The growing interest in the lability of sex determination in non‐model vertebrates such as amphibians and fishes has revealed high rates of sex chromosome turnovers among closely related species of the same clade. Can such lineages hybridize and admix with different sex‐determining systems, or could the changes have precipitated their speciation? We addressed these questions in incipient species of toads (Bufonidae), where we identified a heterogametic transition and characterized their hybrid zone with genome‐wide markers (RADseq). Adult and sibship data confirmed that the common toad B. bufo is female heterogametic (ZW), while its sister species the spined toad B. spinosus is male heterogametic (XY). Analysis of a fine scale transect across their parapatric ranges in southeastern France unveiled a narrow tension zone (∼10 km), with asymmetric mitochondrial and nuclear admixture over hundreds of kilometers southward and northward, respectively. The geographic extent of introgression is consistent with an expansion of B. spinosus across B. bufo’s former ranges in Mediterranean France, as also suggested by species distribution models. However, widespread cyto‐nuclear discordances (B. spinosus backrosses carrying B. bufo mtDNA) run against predictions from the dominance effects of Haldane's rule, perhaps because Y and W heterogametologs are not degenerated. Common and spined toads can thus successfully cross‐breed despite fundamental differences in their sex determination mechanisms, but remain partially separated by reproductive barriers. Whether and how the interactions of their XY and ZW genes contribute to these barriers shall provide novel insights on the debated role of labile sex chromosomes in speciation.

Impact of the Mid-Pleistocene Revolution and Anthropogenic Factors on the Dispersion of Asian Black-Spined Toads (Duttaphrynus melanostictus)

Divergence-time estimation critically improves the understanding of biogeography processes underlying the distribution of species, especially when fossil data is not available. We hypothesise that the Asian black-spined toad, Duttaphrynus melanostictus, expanded into the Eastern Indomalaya following the Quaternary glaciations with the subsequent colonisation of new landscapes during the Last Glacial Maximum. Divergence dating inferred from 364 sequences of mitochondrial tRNAGly ND3 supported the emergence of a common ancestor to the three D. melanostictus clades around 1.85 (±0.77) Ma, matching with the Lower to Mid-Pleistocene transition. Duttaphrynus melanostictus then dispersed into Southeast Asia from the central Indo-Pacific and became isolated in the Southern Sundaic and Wallacea regions 1.43 (±0.10) Ma through vicariance as a result of sea level oscillations. The clade on the Southeast Asian mainland then colonised the peninsula from Myanmar to Vietnam and expanded towards Southeastern China at the end of the Mid-Pleistocene Revolution 0.84 (±0.32) Ma. Population dynamics further highlight an expansion of the Southeast Asian mainland population towards Taiwan, the Northeastern edge of the species' range after the last interglacial, and during the emergence of the Holocene human settlements around 7000 BP. Thus, the current divergence of D. melanostictus into three segregated clades was mostly shaped by Quaternary glaciations, followed by natural dispersion events over land bridges and accelerated by anthropogenic activities.

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

The Caucasian Toad, Bufo verrucosissimus (Pallas, 1814) in the Levant: evidence from mitochondrial DNA

We present the first molecular analysis of the Bufo bufo (Linnaeus, 1758) complex population in the Levant. This complex contains four species distributed through the Western Palearctic. Based on the analysis of the widely used mitochondrial marker 16S rRNA, populations from the Levant belong to B. verrucosissimus (Pallas, 1814) and have a close affiliation to populations from the Transcaucasian region rather than southern Turkey. Such identified relationships may imply past rapid colonisation from north to south and support the need for further research.

Complex hybridization patterns in European pond turtles (Emys orbicularis) in the Pyrenean Region

Hybrid zones are natural laboratories allowing insights in genetic processes like lineage diversification, speciation and introgression. Using large sampling, 15 microsatellite loci and a mitochondrial marker, we examined the Pyrenean contact zone of three pond turtle taxa (Emys orbicularis orbicularis, E. o. galloitalica, E. o. occidentalis). The Pyrenees are a biogeographically important region separating many lineages endemic to the Iberian Peninsula from their Western European counterparts. We discovered limited admixture, reflecting a complex biogeographic scenario. Simulations using Approximate Bayesian Computing supported that E. o. orbicularis invaded the Iberian Peninsula in the Holocene, circumventing the Pyrenees along the Mediterranean coast, and hybridized in the northern peninsula with the local coastal subspecies galloitalica, and to a lesser extent, with occidentalis. While E. o. occidentalis, and in particular E. o. orbicularis, expanded their ranges considerably during Holocene warming, E. o. galloitalica remained largely confined to its former Iberian refuge. Admixture among the three taxa is surprisingly low, and a future taxonomic investigation that includes the unstudied subspecies of E. orbicularis from North Africa, Eastern Europe and Asia has to determine whether their current status properly reflects their evolutionary divergence or whether certain taxa should be regarded as full species.

A common toad hybrid zone that runs from the Atlantic to the Mediterranean

We document the distribution of the common toadBufo bufoand the spined toadB. spinosusat their contact zone across France with data from a mitochondrial DNA RFLP assay, complementing similar work including nuclear markers in the northwest and southeast of France and in Italy. We also reconstruct geographical clines across the species’ contact zone in central France.Bufo bufois found in the north-eastern half of France.Bufo spinosusis found in the south-western complement. The contact zone they form runs from the Atlantic coast near Caen, France, to the Mediterranean coast near Savona, Italy, and has a length of over 900 km. In central FranceB. bufoandB. spinosusengage in a hybrid zone with a unimodal genetic signature. Hybrid zone width is ca. 10 km at mitochondrial DNA and averages at 61 km for four nuclear loci. The hybrid zone is distinctly asymmetric with a signature ofB. spinosusinB. bufoand not the other way round. We attribute this observation toB. bufomoving southwards at the expense ofB. spinosus, with introgression in the direction of the advancing species. We noted substantial geographic variation in characters for species identification. Morphological species identification performs well in France, but breaks down in Italy. Mitochondrial DNA is inconclusive in south-eastern France and Italy. The nuclear genetic markers perform consistently well but have not yet been applied to the zone in full. Possible, but surely heterogeneous ecological correlates for the position of the hybrid zone are mountains and rivers.

Effects of Sample Size and Full Sibs on Genetic Diversity Characterization: A Case Study of Three Syntopic Iberian Pond-Breeding Amphibians

Accurate characterization of genetic diversity is essential for understanding population demography, predicting future trends and implementing efficient conservation policies. For that purpose, molecular markers are routinely developed for nonmodel species, but key questions regarding sampling design, such as calculation of minimum sample sizes or the effect of relatives in the sample, are often neglected. We used accumulation curves and sibship analyses to explore how these 2 factors affect marker performance in the characterization of genetic diversity. We illustrate this approach with the analysis of an empirical dataset including newly optimized microsatellite sets for 3 Iberian amphibian species: Hyla molleri, Epidalea calamita, and Pelophylax perezi. We studied 17-21 populations per species (total n = 547, 652, and 516 individuals, respectively), including a reference locality in which the effect of sample size was explored using larger samples (77-96 individuals). As expected, FIS and tests for Hardy-Weinberg equilibrium and linkage disequilibrium were affected by the presence of full sibs, and most initially inferred disequilibria were no longer statistically significant when full siblings were removed from the sample. We estimated that to obtain reliable estimates, the minimum sample size (potentially including full sibs) was close to 20 for expected heterozygosity, and between 50 and 80 for allelic richness. Our pilot study based on a reference population provided a rigorous assessment of marker properties and the effects of sample size and presence of full sibs in the sample. These examples illustrate the advantages of this approach to produce robust and reliable results for downstream analyses.

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long-distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present-day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.

Neogene amphibians and reptiles (Caudata, Anura, Gekkota, Lacertilia, and Testudines) from the south of Western Siberia, Russia, and Northeastern Kazakhstan

Background

The present-day amphibian and reptile fauna of Western Siberia are the least diverse of the Palaearctic Realm, as a consequence of the unfavourable climatic conditions that predominate in this region. The origin and emergence of these herpetofaunal groups are poorly understood. Aside from the better-explored European Neogene localities yielding amphibian and reptile fossil remains, the Neogene herpetofauna of Western Asia is understudied. The few available data need critical reviews and new interpretations, taking into account the more recent records of the European herpetofauna. The comparison of this previous data with that of European fossil records would provide data on palaeobiogeographic affiliations of the region as well as on the origin and emergence of the present-day fauna of Western Siberia. An overview of the earliest occurrences of certain amphibian lineages is still needed. In addition, studies that address such knowledge gaps can be useful for molecular biologists in their calibration of molecular clocks.

Methods and Results

In this study, we considered critically reviewed available data from amphibian and reptile fauna from over 40 Western Siberian, Russian and Northeastern Kazakhstan localities, ranging from the Middle Miocene to Early Pleistocene. Herein, we provided new interpretations that arose from our assessment of the previously published and new data. More than 50 amphibians and reptile taxa were identified belonging to families Hynobiidae, Cryptobranchidae, Salamandridae, Palaeobatrachidae, Bombinatoridae, Pelobatidae, Hylidae, Bufonidae, Ranidae, Gekkonidae, Lacertidae, and Emydidae. Palaeobiogeographic analyses were performed for these groups and palaeoprecipitation values were estimated for 12 localities, using the bioclimatic analysis of herpetofaunal assemblages.

Conclusion

The Neogene assemblage of Western Siberia was found to be dominated by groups of European affinities, such as Palaeobatrachidae, Bombina, Hyla, Bufo bufo, and a small part of this assemblage included Eastern Palaearctic taxa (e.g. Salamandrella, Tylototriton, Bufotes viridis). For several taxa (e.g. Mioproteus, Hyla, Bombina, Rana temporaria), the Western Siberian occurrences represented their most eastern Eurasian records. The most diverse collection of fossil remains was found in the Middle Miocene. Less diversity has been registered towards the Early Pleistocene, potentially due to the progressive cooling of the climate in the Northern Hemisphere. The results of our study showed higher-amplitude changes of precipitation development in Western Siberia from the Early Miocene to the Pliocene, than previously assumed.

Concordant morphological and molecular clines in a contact zone of the Common and Spined toad (Bufo bufo and B. spinosus) in the northwest of France

BACKGROUND

Hybrid zones are regions where individuals of two species meet and produce hybrid progeny, and are often regarded as natural laboratories to understand the process of species formation. Two microevolutionary processes can take place in hybrid zones, with opposing effects on population differentiation. Hybridization tends to produce genetic homogenization, reducing species differences, whereas the presence of mechanisms of reproductive isolation result in barriers to gene flow, maintaining or increasing differences between taxa.

RESULTS

Here we study a contact zone between two hybridizing toad species, Bufo bufo and B. spinosus, through a combination of molecular (12 polymorphic microsatellites, four nuclear and two mitochondrial SNP markers) and morphological data in a transect in the northwest of France. The results show largely concordant clines across markers, defining a narrow hybrid zone of ca. 30 km wide. Most hybrids in the centre of the contact zone are classified as F₂ or backcrossed individuals, with no individuals assigned to the F₁ hybrid class.

CONCLUSIONS

We discuss the implications of these results for our understanding of the evolutionary history of these species. We anticipate that the toad contact zone here described will become an important asset in the study of hybrid zone dynamics and evolutionary biology because of its easy access and the abundance of the species involved.

Turning into frogs: Asymmetry in forelimb emergence and escape direction in metamorphosing anurans

There is considerable debate about the pattern and origin of laterality in forelimb emergence and turning behaviour within amphibians, with the latter being poorly investigated in tadpoles around metamorphic climax. Using 6 species of metamorphosing anurans, we investigated the effect of asymmetrical spiracle location, and disturbance at the time of forelimb emergence, on the pattern of forelimb emergence. Turning behaviour was observed to assess whether motor lateralization occurred in non-neobatrachian anurans and was linked to patterns of forelimb emergence. Biases in forelimb emergence differed among species, supporting the hypothesis that asymmetrical spiracle position results in the same asymmetry in forelimb emergence. However, this pattern only occurred when individuals were undisturbed. Therefore, context at the time of the emergence of the forelimbs may be important, and might explain some discrepancies in the literature. Turning biases, unconnected to forelimb emergence, were found in Pipidae and Bombinatoridae, confirming the basal origin of lateralized behaviour among anurans. Turning direction in our metamorphs differed from the leftward bias commonly observed in tadpoles, but may be analogous to the prevalent right-"handedness" among adult anurans. Therefore, the transitions occurring during metamorphosis may affect lateralized behaviour and metamorphosis may be fruitful for understanding the development of lateralization.

High genetic diversity of common toad (Bufo bufo) populations under strong natural fragmentation on a Northern archipelago

The last decades have shown a surge in studies focusing on the interplay between fragmented habitats, genetic variation, and conservation. In the present study, we consider the case of a temperate pond-breeding anuran (the common toad Bufo bufo) inhabiting a naturally strongly fragmented habitat at the Northern fringe of the species' range: islands offshore the Norwegian coast. A total of 475 individuals from 19 populations (three mainland populations and 16 populations on seven adjacent islands) were genetically characterized using nine microsatellite markers. As expected for a highly fragmented habitat, genetic distances between populations were high (pairwise F st values ranging between 0.06 and 0.33), with however little differences between populations separated by ocean and populations separated by terrestrial habitat (mainland and on islands). Despite a distinct cline in genetic variation from mainland populations to peripheral islands, the study populations were characterized by overall high genetic variation, in line with effective population sizes derived from single-sample estimators which were on average about 20 individuals. Taken together, our results reinforce the notion that spatial and temporal scales of fragmentation need to be considered when studying the interplay between landscape fragmentation and genetic erosion.

Phylogenetic relationships of toads of the Rhinella granulosa group (Anura: Bufonidae): a molecular perspective with comments on hybridization and introgression

The Rhinella granulosa group consists of 13 species of toads distributed throughout open areas of South America and Panama. In this paper we perform a phylogenetic analysis considering all but one species of the group, employing five nuclear and four mitochondrial genes, for up to 7910 bp per specimen. Separate phylogenetic analyses under direct optimization (DO) of nuclear and mitochondrial sequences recovered the R. granulosa group as monophyletic and revealed topological incongruence that can be explained mainly by multiple events of hybridization and introgression, both mitochondrial and nuclear. The DO combined analysis, after the exclusion of putatively introgressed or heterozygous genomes, resulted in a phylogenetic hypothesis for the R. granulosa group in which most of the species are recovered as monophyletic, but with interspecific relationships poorly supported. The optimization of morphological (adult and larval), chromosomal, and behavioural characters resulted in 12 putative phenotypic synapomorphies for this species group and some other synapomorphies for internal clades. Our results indicate the need for additional population genetic studies on R. dorbignyi and R. fernandezae to corroborate the taxonomic status of both taxa. Finally, we discuss biological and genetic characteristics of Bufonidae, as possible explanations for the common occurrence of hybridization and introgression observed in some lineages of this family.

Phylogeography of the common toad (Bufo bufo, Lissamphibia: Anura) in Switzerland

WhileBufo bufois widespread in Europe, the closely relatedB. spinosusis parapatrically distributed in Western Europe and Northern Africa. The exact course of the contact zone between both taxa, however, is still unknown. 24 samples from Switzerland were analysed using mitochondrial and nuclear markers combined with 243 previously published samples from Europe and Asia. No haplotypes ofB. spinosuswere detected in Switzerland, only two mitochondrial haplotypes ofB. bufo: one north and south of the Alps, the other only north of the Alps. Both haplotypes are also widely distributed in Eastern and Central Europe. These results agree with a postglacial recolonisation of Central Europe from refugia on the Balkan Peninsula. The occurrence of one haplotype north and south of the Alps can be explained by colonisation either from the north by crossing the Alps or from the east along the southern edge of the Alps. The rapid postglacial recolonisation from Balkan refugia might have prevented the spread of southern haplotypes from Italy and ofB. spinosusfrom France.