Countrywide screening supports model-based predictions of the distribution of Batrachochytrium dendrobatidis in Ukraine.

The chytrid Batrachochytrium dendrobatidis (Bd) is a widespread fungus causing amphibian declines across the globe. Although data on Bd occurrence in Eastern Europe are scarce, a recent species distribution model (SDM) for Bd reported that western and north-western parts of Ukraine are highly suitable to the pathogen. We verified the SDM-predicted range of Bd in Ukraine by sampling amphibians across the country and screening for Bd using qPCR. A total of 446 amphibian samples (tissue and skin swabs) from 11 species were collected from 36 localities. We obtained qPCR-positive results for 33 samples including waterfrogs (Pelophylax esculentus complex) and fire- and yellow-bellied toads (Bombina spp.) from 8 localities. We found that Bd-positive localities had significantly higher predicted Bd habitat suitability than sites that were pathogen-free. Amplification and sequencing of the internal transcribed spacer (ITS) region of samples with the highest Bd load revealed matches with ITS haplotypes of the globally distributed BdGPL strain, and a single case of the BdASIA-2/BdBRAZIL haplotype. We found that Bd was non-randomly distributed across Ukraine, with infections present in the western and north-central forested peripheries of the country with a relatively cool, moist climate. On the other hand, our results suggest that Bd is absent or present in low abundance in the more continental central, southern and eastern regions of Ukraine, corroborating the model-predicted distribution of chytrid fungus. These areas could potentially serve as climatic refugia for Bd-susceptible amphibian hosts.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates).

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana.

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.