Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria

Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian-Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area.

A validated protocol for eDNA-based monitoring of within-species genetic diversity in a pond-breeding amphibian

In light of the dramatic decline in amphibian biodiversity, new cost-efficient tools to rapidly monitor species abundance and population genetic diversity in space and time are urgently needed. It has been amply demonstrated that the use of environmental DNA (eDNA) for single-species detection and characterization of community composition can increase the precision of amphibian monitoring compared to traditional (observational) approaches. However, it has been suggested that the efficiency and accuracy of the eDNA approach could be further improved by more timely sampling; in addition, the quality of genetic diversity data derived from the same DNA has been confirmed in other vertebrate taxa, but not amphibians. Given the availability of previous tissue-based genetic data, here we use the common frog Rana temporaria Linnaeus, 1758 as our target species and an improved eDNA protocol to: (i) investigate differences in species detection between three developmental stages in various freshwater environments; and (ii) study the diversity of mitochondrial DNA (mtDNA) haplotypes detected in eDNA (water) samples, by amplifying a specific fragment of the COI gene (331 base pairs, bp) commonly used as a barcode. Our protocol proved to be a reliable tool for monitoring population genetic diversity of this species, and could be a valuable addition to amphibian conservation and wetland management.

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.

Microplastics Occurrence in the European Common Frog (Rana temporaria) from Cottian Alps (Northwest Italy)

Microplastics (MPs) pollution is arousing growing attention, yet knowledge about its occurrence in amphibians is scant to date. With this study, we aimed to determine whether plastic (>5000 μm) and MPs (10–5000 μm) could be detected in adult Rana temporaria from a high-mountain ecosystem (the Cottian Alps, northwest Italy). To do this, aquatic compartments and the digestive tract of adult R. temporaria were analyzed. Water, sediment, periphyton, aquatic macroinvertebrates, and tadpoles tested negative for plastic and MPs. Microplastics were detected in all the adult frogs (n = 5); all the identified items (one per specimen) were fibers (size range: 550.91–2355.51 µm). A statistically significant positive correlation between the particle length and frog size was recorded. The predominant fiber color was blue. The chemical composition was polyamide (60%), polyethylene (20%), and polyethylene terephthalate (20%). Since both the biotic and the abiotic freshwater compartments (tadpoles included) revealed the absence of MPs, it can be assumed that adult frogs ingest MPs from the surrounding terrestrial environment.

Slow and steady wins the race: contrasted phylogeographic signatures in two Alpine amphibians

A deeper phylogeographic structure is expected for slow-dispersing habitat specialists compared to widespread adaptable species, especially in topographically complex regions. We tested this classic assumption by comparing the genomic (RAD-sequencing) phylogeographies of two amphibians inhabiting the Swiss Alps: the mobile, cosmopolitan common frog (Rana temporaria) against the stationary, mountain endemic Alpine salamander (Salamandra atra). Our results ran opposite of predictions: the frog displayed significantly higher genetic divergences and lower within-population variation compared to the salamander. This implies a prominent role for their distinctive glacial histories in shaping intraspecific diversity and structure: diversification and recolonization from several circum-Alpine micro-refugia for the frog versus a single refugium for the salamander, potentially combined with better population connectivity and stability. These striking differences emphasize the great variability of phylogeographic responses to the Quaternary glaciations, hence the complexity to predict general patterns of genetic diversity at the regional scale, and the forces that underlie them.

Refugia within refugium of Geranium yesoense (Geraniaceae) in Japan were driven by recolonization into the southern interglacial refugium

Recent studies have found that geographical fragmentation and recurrent colonization result in complex genetic structures in refugial areas. This phenomenon, known as ‘refugia within refugium’, has been identified from many geographical locations. In Japan, the high-elevation mountains of central Honshu provided an interglacial refugium for alpine plants. Here we focused on the Geranium yesoense complex, which exhibits increased morphological variation in the refugial area, to determine whether this variation was shaped by recurrent colonization, range fragmentation or phenotypic changes independent of population history. We analysed single nucleotide polymorphism data and chloroplast genome sequences. Diversification in the G. yesoense species complex occurred in the mid-Pleistocene. The varieties are distinct entities and suggest the presence of a genetic cluster with highly disjunct distributions, occurring both in northern Japan and in southern refugial areas in central Honshu. Demographic analysis suggests that a single ancestral variety (var. nipponicum) evolved in the alpine region of central Honshu, and that subsequent migration from one of the two diverged northern varieties (var. pseudopratense) led to secondary contact with var. nipponicum during the last glacial period. Recolonization into refugial populations in central Honshu and hybridization between diverged populations have resulted in complex genetic structures among refugial populations.

Repeated hybridization increased diversity in the door snail complex Charpentieria itala in the Southern Alps

The door snail species complex Charpentieria itala is widely distributed in the Southern Alps and subdivided into several morphologically differentiated subspecies. Thus, it can be used as a model group for understanding migration and differentiation processes in the Southern Alps. We generated genome-wide double digest Restriction Site Associated DNA (ddRAD) sequencing data for 166 specimens from 36 populations of the door snail Charpentieria itala and for 8 specimens of the other three Charpentieria species to reconstruct their evolutionary history and phylogeography. Phylogenetic and structure analyses based on the ddRAD data indicated that the repeated separation of the populations in western and eastern groups by the Garda glacier during the glacials was the process that most strongly shaped the population structure of C. itala. This process may also explain a similar phylogeographic boundary in many other southern Alpine animal and plant species. Our study revealed that some populations that resemble Charpentieria stenzii morphologically and ecologically, the 'stenzioid' subspecies, originated by a hybridization event with Charpentieria stenzii. A further hybridization event between stenzioid populations that survived the glacials in mountain refuges and non-stenzioid populations that probably came into contact with stenzioid populations as a result of climate warming during an interglacial resulted in the origin of a hybrid subspecies that is adapted to intermediate altitudes. Our study demonstrated that the origin of new differentiated taxa by hybridization, is more frequent than previously assumed.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Deciphering the drivers of negative species-genetic diversity correlation in Alpine amphibians

The evolutionary and ecological importance of neutral and adaptive genetic diversity is widely recognized. Nevertheless, genetic diversity is rarely assessed for conservation planning, which often implicitly assumes a positive correlation between species and genetic diversity. Multiple drivers can cause the co-variation between the genetic diversity of one species and the richness of the whole communities, and explicit tests are needed to identify the processes that can determine species-genetic diversity correlations (SGDCs). Here, we tested whether intrapopulation genetic diversity (at neutral loci) and species richness co-vary in the amphibian communities of a southern Alpine region (Trentino, Italy), using the common frog (Rana temporaria) as focal species for the study of genetic diversity. We also analysed ecological similarity, niche overlap and interspecific interactions between the species, to unravel the processes determining SGDC. The neutral genetic diversity of common frogs was negatively related to species richness. The negative SGDC was probably due to an opposite influence of environmental gradients on the two levels of biodiversity, since the focal species and the other amphibians differ in ecological preferences, particularly in terms of thermal optimum. Conversely, we did not find evidence for a role of interspecific interactions in the negative SGDC. Our findings stress that species richness cannot be used as a universal proxy for genetic diversity, and only combining SGDC with analyses on the determinants of biodiversity can allow to identify the processes determining the relationships between genetic and species diversity.