Combining demography and genetic analysis to assess the population structure of an amphibian in a human-dominated landscape

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure, genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

Multi-Scale Ecological Connectivity Dynamics Associated With Hydropower Station: A Case Study in the Lancang River Valley

Dam construction is a major threat to terrestrial ecological processes that exhibit spatial variation. As an effective indicator of species movements at large scales, the dynamics of forest landscape connectivity for the dispersal abilities of local species such as primates both upstream and downstream of the Manwan hydropower station were compared in three periods: before, during, and after dam construction (in 1974, 1988, and 2004, respectively). The equivalent connected area (ECA) index which was modified from the probability of connectivity (PC) index based on the graph theory was applied in this study. We used this method to analyze the spatial and temporal changes in the overall forest landscape connectivity due to dam construction. The results showed that forest connectivity decreased after dam construction. The forest connectivity upstream was much greater than that downstream in 1974 but was lower in 1988 and 2004. The importance of connectivity increased with increasing dispersal distance. Only 8, 12, and 18% of forest patches of small area exhibited changes in 1974, 1988, and 2004, respectively. Most of the large habitat patches (i.e., those with dECA values >50%) remained stable, and the greatest patch changes were found at a dispersal distance of 400 m in all three periods. These large forest patch changes often occurred near the boundary of the study area. The Betweenness Centrality indicator, which identifies patches as stepping stones while accounting for ecological processes and biological flows at a larger scale, indicated that some habitat patches near Manwan Dam acted as stepping stones in maintaining the forest connectivity. Furthermore, rank correlations between the forest patch area and the three dPC (intra, flux, connector) fractions indicated that the dPCflux indicator can be used to measure the prioritization of habitat patches.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

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.

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.

Movement ecology of amphibians: from individual migratory behaviour to spatially structured populations in heterogeneous landscapes,

Both genetic cohesion among local populations of animals and range expansion depend on the frequency of dispersers moving at an interpatch scale. Animal movement has an individual component that reflects behaviour and an ecological component that reflects the spatial organization of populations. The total movement capacity of an individual describes maximum movement distance theoretically achievable during a lifetime, whereas its variation among the members of a local population determines the magnitude of interpatch movements and thus of gene flow between neighbouring patches within metapopulation or patchy population systems. Here, I review information on dispersal and migration as components of the movement capacity of juvenile and adult pond-breeding amphibians and discuss how these components inform the spatial structure of populations. Amphibians disperse as juveniles and adults, but movement distances detected in tracking or capture–mark–recapture studies are usually far below the corresponding estimates based on molecular gene-flow data. This discrepancy reflects the constraints of available tracking methods for free-ranging individuals leading to inappropriate surrogates of annual movement capacity, but can be resolved using probabilistic approaches based on dispersal functions. There is remarkable capacity for and plasticity in movements in amphibians. Annual within-patch movements (migrations) of individuals can be large and likely represent an underestimated capacity for movement at the interpatch scale. Landscape resistance may influence the paths of dispersing amphibians, but rarely impedes interpatch movements. Juveniles emigrating unpredictably far from the natal pond and adults switching from within-patch migrations to dispersal to another patch demonstrate the plasticity of individual movement behaviour. Three basic conclusions can be drawn with respect to the linkage of individual movement behaviour and spatial or genetic structure of local amphibian populations embedded in a heterogeneous landscape: (1) individual movements or consecutive short-term series of movements are misleading surrogate measures of total movement capacity; (2) probabilistic modelling of movement capacity is the best available behavioural predictor of interpatch gene flow; (3) connectivity of local populations in heterogeneous landscapes is less affected by landscape resistance than previously expected.

Optimizing the trade-off between spatial and genetic sampling efforts in patchy populations: towards a better assessment of functional connectivity using an individual-based sampling scheme

Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organisms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic distances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population-based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sampling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analyses on simulated and empirical genetic data sets to investigate the efficiency of an individual-based sampling scheme (ISS) in detecting isolation-by-distance and isolation-by-barrier patterns. Provided that pseudo-replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of interindividual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situations. The ISS proved to be a promising methodological alternative to the more conventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.

Effects of road mortality and mitigation measures on amphibian populations

Road mortality is a widely recognized but rarely quantified threat to the viability of amphibian populations. The global extent of the problem is substantial and factors affecting the number of animals killed on highways include life-history traits and landscape features. Secondary effects include genetic isolation due to roads acting as barriers to migration. Long-term effects of roads on population dynamics are often severe and mitigation methods include volunteer rescues and under-road tunnels. Despite the development of methods that reduce road kill in specific locations, especially under-road tunnels and culverts, there is scant evidence that such measures will protect populations over the long term. There also seems little likelihood that funding will be forthcoming to ameliorate the problem at the scale necessary to prevent further population declines.

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.

Using genetic variation to infer associations with climate in the common frog, Rana temporaria

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species-level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine-scale population structure without consideration of confounding historical variation. Second, we assessed fine-scale microsatellite-based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average FST=0.02), despite a mean annual temperature difference of 4.5 °C between low- and high-altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.

Genetic structure of a Japanese brown frog (Rana japonica) population implies severe restriction of gene flow caused by recent urbanization in a satoyama landscape

Rapid urbanization is one of the major pressures on amphibian species. Elucidating changes in genetic structure will be useful in evaluating the effects of urbanization on amphibian populations. Our study focused on Rana japonica, which is common in complex agricultural landscapes known as satoyama, which are also under intense development pressure. We conducted landscape genetic analyses based on mitochondrial DNA haplotype frequencies of 13 breeding sites in a rapidly urbanizing area of Japan. We found several breeding sites had significantly higher F(st) values, and we also identified the barriers to gene flow between these sites. Observation of past aerial photographs revealed that these barriers coincided with the construction of man-made structures in the last few decades, suggesting that urbanization has restricted gene flow in R. japonica. Our results show that landscape genetic approaches are useful in conservation planning where rapid habitat degradation has taken place.

Permeability of the landscape matrix between amphibian breeding sites

For organisms that reproduce in discrete habitat patches, land cover between patches (known as the matrix) is important for dispersal among breeding sites. Models of patchy populations often incorporate information on the permeability of the matrix to dispersal, sometimes based on expert opinion. I estimated the relative resistance to gene flow of land cover types and barriers using F(ST) calculated from microsatellite markers in two amphibians, within an 800-km(2) area in northern Switzerland. The species included a frog (Rana temporaria: 996 individuals, 48 populations, seven markers) and a newt (Triturus alpestris: 816 individuals, 41 populations, seven markers). Open fields and urban areas were more resistant to gene flow than forested land; roads and highways also reduced permeability. Results were similar for the two species. However, differences in resistance among matrix elements were relatively low: gene flow through urban areas was reduced by only 24-42% relative to forest; a divided highway reduced gene flow by 11-40% and was 7-8 times more resistant than a secondary road. These data offer an empirically based alternative to expert opinion for setting relative resistance values in landscape models.

Intra-specific variation in nitrate tolerance in tadpoles of the Natterjack toad

Anthropogenic sources of nitrogen that pollute bodies of water can have toxic and sub-lethal effects on amphibians. It has been hypothesized that such exposure may promote local adaptation, that is, selection for higher tolerance in individuals in populations exposed to pollutants. We tested this hypothesis with respect to the Natterjack toad (Bufo calamita Laurenti, 1768), by comparing the nitrate dose response of tadpoles from eight populations (doses: 0, 50, 100, 500 and 1000 mg/l nitrate) from relatively unpolluted and intensively farmed environments. We evaluated the effect of nitrate exposure by observing the behavior (movements) of tadpoles exposed to different concentrations of nitrates. Exposure to high nitrate levels did not cause tadpole mortality in the populations used in our experiments; however, we did observe changes in activity for all populations, with these changes being either dose-related responses (decreased activity after exposure to 500 or 1000 mg/l), or more complex responses (increased activity when exposed to 50 or 100 mg/l nitrate, followed by decreased activity at higher concentrations). Natterjack toad tadpoles exhibited variable behavioural responses among the tested populations. Although these populations were selected on the basis of their potential agrochemical contamination, the observed variation in population tolerance was not related to the parameters used to estimate this contamination in these breeding sites. Possible explanations for this apparent lack of local adaptation in B. calamita tadpoles include inadequate estimates of the toads' actual nitrate exposure in the field, and the biological characteristics of B. calamita, which may limit the effects of exposure or favor phenotypic plasticity.

Comparison of Bayesian clustering and edge detection methods for inferring boundaries in landscape genetics

Recently, techniques available for identifying clusters of individuals or boundaries between clusters using genetic data from natural populations have expanded rapidly. Consequently, there is a need to evaluate these different techniques. We used spatially-explicit simulation models to compare three spatial Bayesian clustering programs and two edge detection methods. Spatially-structured populations were simulated where a continuous population was subdivided by barriers. We evaluated the ability of each method to correctly identify boundary locations while varying: (i) time after divergence, (ii) strength of isolation by distance, (iii) level of genetic diversity, and (iv) amount of gene flow across barriers. To further evaluate the methods' effectiveness to detect genetic clusters in natural populations, we used previously published data on North American pumas and a European shrub. Our results show that with simulated and empirical data, the Bayesian spatial clustering algorithms outperformed direct edge detection methods. All methods incorrectly detected boundaries in the presence of strong patterns of isolation by distance. Based on this finding, we support the application of Bayesian spatial clustering algorithms for boundary detection in empirical datasets, with necessary tests for the influence of isolation by distance.