Reducing bias in population and landscape genetic inferences: the effects of sampling related individuals and multiple life stages

Disease-driven top predator decline affects mesopredator population genomic structure

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

Practical application of the linkage disequilibrium method for estimating contemporary effective population size: A review

The method to estimate contemporary effective population size (Ne ) based on patterns of linkage disequilibrium (LD) at unlinked loci has been widely applied to natural and managed populations. The underlying model makes many simplifying assumptions, most of which have been evaluated in numerous studies published over the last two decades. Here, these performance evaluations are reviewed and summarized, with a focus on information that facilitates practical application to real populations in nature. Potential sources of bias that are discussed include calculation of r2 (a measure of LD), adjustments for sampling error, physical linkage, age structure, migration and spatial structure, mutation and selection, mating systems, changes in abundance, rare alleles, missing data, genotyping errors, data filtering choices and methods for combining multiple Ne estimates. Factors that affect precision are reviewed, including pseudoreplication that limits the information gained from large genomics datasets, constraints imposed by small samples of individuals, and the challenges in obtaining robust estimates for large populations. Topics that merit further research include the potential to weight r2 values by allele frequency, lump samples of individuals, use genotypic likelihoods rather than called genotypes, prune large LD values and apply the method to species practising partial monogamy.

Terrain Ruggedness and Canopy Height Predict Short-Range Dispersal in the Critically Endangered Black-and-White Ruffed Lemur

Dispersal is a fundamental aspect of primates' lives and influences both population and community structuring, as well as species evolution. Primates disperse within an environmental context, where both local and intervening environmental factors affect all phases of dispersal. To date, research has primarily focused on how the intervening landscape influences primate dispersal, with few assessing the effects of local habitat characteristics. Here, we use a landscape genetics approach to examine between- and within-site environmental drivers of short-range black-and-white ruffed lemur (Varecia variegata) dispersal in the Ranomafana region of southeastern Madagascar. We identified the most influential drivers of short-range ruffed lemur dispersal as being between-site terrain ruggedness and canopy height, more so than any within-site habitat characteristic evaluated. Our results suggest that ruffed lemurs disperse through the least rugged terrain that enables them to remain within their preferred tall-canopied forest habitat. Furthermore, we noted a scale-dependent environmental effect when comparing our results to earlier landscape characteristics identified as driving long-range ruffed lemur dispersal. We found that forest structure drives short-range dispersal events, whereas forest presence facilitates long-range dispersal and multigenerational gene flow. Together, our findings highlight the importance of retaining high-quality forests and forest continuity to facilitate dispersal and maintain functional connectivity in ruffed lemurs.

Life-history stage and the population genetics of the tiger mosquito Aedes albopictus at a fine spatial scale

As a widespread vector of disease with an expanding range, the mosquito Aedes albopictus Skuse (Diptera: Culicidae) is a high priority for research and management. A. albopictus has a complex life history with aquatic egg, larval and pupal stages, and a terrestrial adult stage. This requires targeted management strategies for each life stage, coordinated across time and space. Population genetics can aid in A. albopictus control by evaluating patterns of genetic diversity and dispersal. However, how life stage impacts population genetic characteristics is unknown. We examined whether patterns of A. albopictus genetic diversity and differentiation changed with life stage at a spatial scale relevant to management efforts. We first conducted a literature review of field-caught A. albopictus population genetic papers and identified 101 peer-reviewed publications, none of which compared results between life stages. Our study uniquely examines population genomic patterns of egg and adult A. albopictus at five sites in Wake County, North Carolina, USA, using 8425 single nucleotide polymorphisms. We found that the level of genetic diversity and connectivity between sites varied between adults and eggs. This warrants further study and is critical for research aimed at informing local management.

The first assessment of the genetic diversity and structure of the endangered West Indian manatee in Cuba

The coastal waters of Cuba are home to a small, endangered population of West Indian manatee, which would benefit from a comprehensive characterization of the population's genetic variation. We conducted the first genetic assessment of Cuban manatees to determine the extent of the population's genetic structure and characterize the neutral genetic diversity among regions within the archipelago. We genotyped 49 manatees at 18 microsatellite loci, a subset of 27 samples on 1703 single nucleotide polymorphisms (SNPs), and sequenced 59 manatees at the mitochondrial control region. The Cuba manatee population had low nuclear (microsatellites H_E = 0.44, and SNP H_E = 0.29) and mitochondrial genetic diversity (h = 0.068 and π = 0.00025), and displayed moderate departures from random mating (microsatellite F_IS = 0.12, SNP F_IS = 0.10). Our results suggest that the western portion of the archipelago undergoes periodic exchange of alleles based on the evidence of shared ancestry and low but significant differentiation. The southeast Guantanamo Bay region and the western portion of the archipelago were more differentiated than southwest and northwest manatees. The genetic distinctiveness observed in the southeast supports its recognition as a demographically independent unit for natural resource management regardless of whether it is due to historical isolation or isolation by distance. Estimates of the regional effective population sizes, with the microsatellite and SNP datasets, were small (all N_e < 60). Subsequent analyses using additional samples could better examine how the observed structure is masking simple isolation by distance patterns or whether ecological or biogeographic forces shape genetic patterns.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

Genomic assessment of global population structure in a highly migratory and habitat versatile apex predator, the tiger shark (Galeocerdo cuvier)

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST=0.125-0.144) and candidate outlier loci (FST=0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared to sharks from the Indian Ocean (FST=0.003-0.005, P<0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Evaluation of the Minimum Sampling Design for Population Genomic and Microsatellite Studies: An Analysis Based on Wild Maize

Massive parallel sequencing (MPS) is revolutionizing the field of molecular ecology by allowing us to understand better the evolutionary history of populations and species, and to detect genomic regions that could be under selection. However, the economic and computational resources needed generate a tradeoff between the amount of loci that can be obtained and the number of populations or individuals that can be sequenced. In this work, we analyzed and compared two simulated genomic datasets fitting a hierarchical structure, two extensive empirical genomic datasets, and a dataset comprising microsatellite information. For all datasets, we generated different subsampling designs by changing the number of loci, individuals, populations, and individuals per population to test for deviations in classic population genetics parameters (H_S, F_IS, F_ST). For the empirical datasets we also analyzed the effect of sampling design on landscape genetic tests (isolation by distance and environment, central abundance hypothesis). We also tested the effect of sampling a different number of populations in the detection of outlier SNPs. We found that the microsatellite dataset is very sensitive to the number of individuals sampled when obtaining summary statistics. F_IS was particularly sensitive to a low sampling of individuals in the simulated, genomic, and microsatellite datasets. For the empirical and simulated genomic datasets, we found that as long as many populations are sampled, few individuals and loci are needed. For the empirical datasets, we found that increasing the number of populations sampled was important in obtaining precise landscape genetic estimates. Finally, we corroborated that outlier tests are sensitive to the number of populations sampled. We conclude by proposing different sampling designs depending on the objectives.

Genetic consequences of social structure in the golden-crowned sifaka

Many species are structured in social groups (SGs) where individuals exhibit complex mating strategies. Yet, most population genetic studies ignore SGs either treating them as small random-mating units or focusing on a higher hierarchical level (the population). Empirical studies acknowledging SGs have found an overall excess of heterozygotes within SGs and usually invoke inbreeding avoidance strategies to explain this finding. However, there is a lack of null models against which ecological theories can be tested and inbreeding avoidance quantified. Here, we investigate inbreeding (deviation from random mating) in an endangered forest-dwelling pair-living lemur species (Propithecus tattersalli). In particular, we measure the inbreeding coefficient (F_IS) in empirical data at different scales: SGs, sampling sites and forest patches. We observe high excess of heterozygotes within SGs. The magnitude of this excess is highly dependent on the sampling scheme: while offspring are characterised by a high excess of heterozygotes (F_IS < 0), the reproductive pair does not show dramatic departures from Hardy-Weinberg expectations. Moreover, the heterozygosity excess disappears at larger geographic scales (sites and forests). We use a modelling framework that incorporates details of the sifaka mating system but does not include active inbreeding avoidance mechanisms. The simulated data show that, although apparent "random mating" or even inbreeding may occur at the "population" level, outbreeding is maintained within SGs. Altogether our results suggest that social structure leads to high levels of outbreeding without the need for active inbreeding avoidance mechanisms. Thus, demonstrating and measuring the existence of active inbreeding avoidance mechanisms may be more difficult than usually assumed.

Replicated Landscape Genomics Identifies Evidence of Local Adaptation to Urbanization in Wood Frogs

Native species that persist in urban environments may benefit from local adaptation to novel selection factors. We used double-digest restriction-side associated DNA (RAD) sequencing to evaluate shifts in genome-wide genetic diversity and investigate the presence of parallel evolution associated with urban-specific selection factors in wood frogs (Lithobates sylvaticus). Our replicated paired study design involved 12 individuals from each of 4 rural and urban populations to improve our confidence that detected signals of selection are indeed associated with urbanization. Genetic diversity measures were less for urban populations; however, the effect size was small, suggesting little biological consequence. Using an FST outlier approach, we identified 37 of 8344 genotyped single nucleotide polymorphisms with consistent evidence of directional selection across replicates. A genome-wide association study analysis detected modest support for an association between environment type and 12 of the 37 FST outlier loci. Discriminant analysis of principal components using the 37 FST outlier loci produced correct reassignment for 87.5% of rural samples and 93.8% of urban samples. Eighteen of the 37 FST outlier loci mapped to the American bullfrog (Rana [Lithobates] catesbeiana) genome, although none were in coding regions. This evidence of parallel evolution to urban environments provides a powerful example of the ability of urban landscapes to direct evolutionary processes.

Fine-scale population genetic structure of dengue mosquito vector, Aedes aegypti, in Metropolitan Manila, Philippines

Dengue is a highly endemic disease in Southeast Asia and is transmitted primarily by the mosquito, Aedes aegypti. The National Capital Region (NCR) of the Philippines, or Metropolitan Manila, is a highly urbanized area that is greatly affected by this arboviral disease. Urbanization has been shown to increase the dispersal of this mosquito vector. For this reason, we conducted a fine-scale population genetic study of Ae. aegypti in this region. We collected adult Ae. aegypti mosquitoes (n = 526 individuals) within the region (n = 21 study areas) and characterized the present population structure and the genetic relatedness among mosquito populations. We genotyped 11 microsatellite loci from all sampled mosquito individuals and analyzed their genetic diversity, differentiation and structure. The results revealed low genetic differentiation across mosquito populations which suggest high gene flow and/or weak genetic drift among mosquito populations. Bayesian analysis indicated multiple genetic structures (K = 3-6), with no clear genetically distinct population structures. This result implies the passive or long-distance dispersal capability nature Ae. aegypti possibly through human-mediated transportation. The constructed dendrogram in this study describes the potential passive dispersal patterns across Metropolitan Manila. Furthermore, spatial autocorrelation analysis showed the limited and active dispersal capability (<1km) of the mosquito vector. Our findings are consistent with previous studies that investigated the genetic structure and dual (active and passive) dispersal capability of Ae. aegypti in a fine-scale highly urbanized area.

Landscape genetics reveals unique and shared effects of urbanization for two sympatric pool-breeding amphibians

Metapopulation-structured species can be negatively affected when landscape fragmentation impairs connectivity. We investigated the effects of urbanization on genetic diversity and gene flow for two sympatric amphibian species, spotted salamanders (Ambystoma maculatum) and wood frogs (Lithobates sylvaticus), across a large (>35,000 km2) landscape in Maine, USA, containing numerous natural and anthropogenic gradients. Isolation-by-distance (IBD) patterns differed between the species. Spotted salamanders showed a linear and relatively high variance relationship between genetic and geographic distances (r = .057, p < .001), whereas wood frogs exhibited a strongly nonlinear and lower variance relationship (r = 0.429, p < .001). Scale dependence analysis of IBD found gene flow has its most predictable influence (strongest IBD correlations) at distances up to 9 km for spotted salamanders and up to 6 km for wood frogs. Estimated effective migration surfaces revealed contrasting patterns of high and low genetic diversity and gene flow between the two species. Population isolation, quantified as the mean IBD residuals for each population, was associated with local urbanization and less genetic diversity in both species. The influence of geographic proximity and urbanization on population connectivity was further supported by distance-based redundancy analysis and multiple matrix regression with randomization. Resistance surface modeling found interpopulation connectivity to be influenced by developed land cover, light roads, interstates, and topography for both species, plus secondary roads and rivers for wood frogs. Our results highlight the influence of anthropogenic landscape features within the context of natural features and broad spatial genetic patterns, in turn supporting the premise that while urbanization significantly restricts interpopulation connectivity for wood frogs and spotted salamanders, specific landscape elements have unique effects on these two sympatric species.

Post-Pleistocene differentiation in a Central Interior Highlands endemic salamander

AIM

For many endemic species with limited dispersal capacities, the relationship between landscape changes and species distributions is still unclear. We characterized the population structure of the endemic ringed salamander (Ambystoma annulatum) across its distribution in the Central Interior Highlands (CIH) of North America, an area of high species endemism, to infer the ecological and evolutionary history of the species.

METHODS

We sampled 498 individuals across the species distribution and characterized the population genetic structure using nuclear microsatellite and mitochondrial DNA (mtDNA) markers.

RESULTS

Ambystoma annulatum exist in two strongly supported nuclear genetic clusters across the CIH that correspond to a northern cluster that includes the Missouri Ozark populations and a southern cluster that includes the Arkansas and Oklahoma Ozarks and the Ouachita Mountains. Our demographic models estimated that these populations diverged approximately 2,700 years ago. Pairwise estimates of genetic differentiation at microsatellite and mtDNA markers indicated limited contemporary gene flow and suggest that genetic differentiation was primarily influenced by changes in the post-Pleistocene landscape of the CIH.

MAIN CONCLUSIONS

Both the geologic history and post-European settlement history of the CIH have influenced the population genetic structure of A. annulatum. The low mtDNA diversity suggests a retraction into and expansion out of refugial areas in the south-central Ozarks, during temperature fluctuations of the Pleistocene and Holocene epochs. Similarly, the estimated divergence time for the two nuclear clusters corresponds to changes in the post-Pleistocene landscape. More recently, decreased A. annulatum gene flow may be a result of increased habitat fragmentation and alteration post-European settlement.

Single Nucleotide Polymorphisms (SNPs) Reveal Sibship Among Founders of a Bangladeshi Rohu (Labeo rohita) Breeding Population

Rohu (Labeo rohita) is a significant freshwater aquaculture species with approximately 1.8 Mt produced annually. Fin clips obtained from the founders of a newly established Bangladesh-based breeding population (∼140 fish from each of the Halda, Jamuna, and Padma rivers) were used to identify 9157 SNPs and 14 411 silicoDArT markers using the Diversity Arrays Technology (DArT) genotyping-by-sequencing platform known as DArTseq. After quality control, 1985 SNPs were retained and used to examine population structure within and among river systems. Examination of genomic relationships revealed evidence of full- and half-sibling relationships among founders. Accordingly, sibship and dummy parents were assigned within each river population using a maximum likelihood approach with COLONY software. Founders that had no dummy parents in common were then identified for population genetic analyses. Only 40 unique dummy parents and 17 founders with no common dummy parents were identified from the Halda river, compared with 206 (96) from the Jamuna and 184 (83) from the Padma. Overall pairwise FST estimates among rivers were low (<0.005) and the optimum number of clusters using unsupervised K-means clustering was one, indicating little genetic divergence among the river populations in our SNPs. These results suggest that observed sibship among founders should be accounted for in future pedigree-based analyses and it cannot be assumed that fertilized spawn collections are representative samples of river populations.

Sampling related individuals within ponds biases estimates of population structure in a pond-breeding amphibian

Effective conservation and management of pond-breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding-site connectivity. Population-level studies of pond-breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within-pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome-wide SNPs generated using a double-digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within-pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.

Reducing cryptic relatedness in genomic data sets via a central node exclusion algorithm

Cryptic relatedness is a confounding factor in genetic diversity and genetic association studies. Development of strategies to reduce cryptic relatedness in a sample is a crucial step for downstream genetic analyses. This study uses a node selection algorithm, based on network degrees of centrality, to evaluate its applicability and impact on evaluation of genetic diversity and population stratification. 1,036 Guzerá (Bos indicus) females were genotyped using Illumina Bovine SNP50 v2 BeadChip. Four strategies were compared. The first and second strategies consist on a iterative exclusion of most related individuals based on PLINK kinship coefficient (φij) and VanRaden's φij, respectively. The third and fourth strategies were based on a node selection algorithm. The fourth strategy, Network G matrix, preserved the larger number of individuals with a better diversity and representation from the initial sample. Determining the most probable number of populations was directly affected by the kinship metric. Network G matrix was the better strategy for reducing relatedness due to producing a larger sample, with more distant individuals, a more similar distribution when compared with the full data set in the MDS plots and keeping a better representation of the population structure. Resampling strategies using VanRaden's φij as a relationship metric was better to infer the relationships among individuals. Moreover, the resampling strategies directly impact the genomic inflation values in genomewide association studies. The use of the node selection algorithm also implies better selection of the most central individuals to be removed, providing a more representative sample.

Amphibian population genetics in agricultural landscapes: does viniculture drive the population structuring of the European common frog (Rana temporaria)?

Amphibian populations have been declining globally over the past decades. The intensification of agriculture, habitat loss, fragmentation of populations and toxic substances in the environment are considered as driving factors for this decline. Today, about 50% of the area of Germany is used for agriculture and is inhabited by a diverse variety of 20 amphibian species. Of these, 19 are exhibiting declining populations. Due to the protection status of native amphibian species, it is important to evaluate the effect of land use and associated stressors (such as road mortality and pesticide toxicity) on the genetic population structure of amphibians in agricultural landscapes. We investigated the effects of viniculture on the genetic differentiation of European common frog (Rana temporaria) populations in Southern Palatinate (Germany). We analyzed microsatellite data of ten loci from ten breeding pond populations located within viniculture landscape and in the adjacent forest block and compared these results with a previously developed landscape permeability model. We tested for significant correlation of genetic population differentiation and landscape elements, including land use as well as roads and their associated traffic intensity, to explain the genetic structure in the study area. Genetic differentiation among forest populations was significantly lower (median pairwise F_ST = 0.0041 at 5.39 km to 0.0159 at 9.40 km distance) than between viniculture populations (median pairwise F_ST = 0.0215 at 2.34 km to 0.0987 at 2.39 km distance). Our analyses rejected isolation by distance based on roads and associated traffic intensity as the sole explanation of the genetic differentiation and suggest that the viniculture landscape has to be considered as a limiting barrier for R. temporaria migration, partially confirming the isolation of breeding ponds predicted by the landscape permeability model. Therefore, arable land may act as a sink habitat, inhibiting genetic exchange and causing genetic differentiation of pond populations in agricultural areas. In viniculture, pesticides could be a driving factor for the observed genetic impoverishment, since pesticides are more frequently applied than any other management measure and can be highly toxic for terrestrial life stages of amphibians.