Linking extinction-colonization dynamics to genetic structure in a salamander metapopulation

Genetic variability of the endangered fish lake minnow (Eupallasella percnurus) in populations newly established by translocation and those existing long term in Poland

The lake minnow Eupallasella percnurus is a small leuciscid fish. In Poland, this species has been in a continuous decline since the mid-20th century and is presently considered as a extremely endangered. According to Polish law, E. percnurus is a strictly protected species that requires active conservation measures. In Poland, one the most common and effective measure of active protection E. percnurus is initiation of new populations. For this purpose, in 2004-2012, juvenile individuals originating from aquaculture conditions were translocated to group of isolated water bodies not inhabited by this species. The juveniles were offspring of parental fish belonging to the same local population, which is extinct at present. Five of those attempts were successful. The aim of the present study was to assess the genetic variation in a group new populations and compare genetic variation indicators with 13 old populations that had existed for decades. The polymorphism of 13 microsatellite markers was investigated, significance of differences in the genetic variation indicators between the groups were tested using a one-way analysis of variance (ANOVA). The mean values of all summary statistics under study, i.e. observed heterozygosity, expected heterozygosity and the total number of alleles, were higher in the group of new populations compared to almost all old ones. A similar dependence was found for Garza-Williamson M values, where the mean for the group of new populations was higher than in almost all old populations. Our results indicate that all recently established E. percnurus populations have not yet experienced any extensive founder effects or bottlenecks. They have preserved a large part of the genetic variability typical of their maternal population, which might also have been relatively high. This feature of new populations, may give them a relatively high ability to adapt to changing environments in the future.

Phylogeography and phenotypic wing shape variation in a damselfly across populations in Europe

BACKGROUND

Describing geographical variation in morphology of organisms in combination with data on genetic differentiation and biogeography can provide important information on how natural selection shapes such variation. Here we study genetic structure using ddRAD seq and wing shape variation using geometric morphometrics in 14 populations of the damselfly Lestes sponsa along its latitudinal range in Europe.

RESULTS

The genetic analysis showed a significant, yet relatively weak population structure with high genetic heterozygosity and low inbreeding coefficients, indicating that neutral processes contributed very little to the observed wing shape differences. The genetic analysis also showed that some regions of the genome (about 10%) are putatively shaped by selection. The phylogenetic analysis showed that the Spanish and French populations were the ancestral ones with northern Swedish and Finnish populations being the most derived ones. We found that wing shape differed significantly among populations and showed a significant quadratic (but weak) relationship with latitude. This latitudinal relationship was largely attributed to allometric effects of wing size, but non-allometric variation also explained a portion of this relationship. However, wing shape showed no phylogenetic signal suggesting that lineage-specific variation did not contribute to the variation along the latitudinal gradient. In contrast, wing size, which is correlated with body size in L. sponsa, had a strong negative correlation with latitude.

CONCLUSION

Our results suggest a relatively weak population structure among the sampled populations across Europe, but a clear differentiation between south and north populations. The observed geographic phenotypic variation in wing shape may have been affected by different local selection pressures or environmental effects.

Genetic diversity, gene flow, and landscape resistance in a pond-breeding amphibian in agricultural and natural forested landscapes in Norway

Genetic diversity is a key part of biodiversity, threatened by human activities that lead to loss of gene flow and reduction of effective population sizes. Gene flow is a result of both landscape connectivity and demographic processes determining the number of dispersing individuals in space and time. Thus, the effect of human impact on processes determining the level of genetic diversity must be interpreted in the context of basic ecological conditions affecting survival and recruitment. When the intensity of human impact and habitat suitability correlate, the effect on genetic diversity and gene flow may be challenging to predict. We compared genetic diversity, gene flow and landscape resistance in two contrasting landscapes in Norway for the pond-breeding amphibian Triturus cristatus: a highly human-impacted, agricultural landscape with ecologically productive habitats, and a forested landscape with less productive habitats and lower levels of human impact. Our results show that genetic diversity was higher and gene flow lower within the forested landscape. Microclimatic moisture conditions and vegetation cover were important determinants of landscape resistance to gene flow within both landscapes. There were indications that landscape resistance was increased by minor roads in the forested landscape, which was not the case for the agricultural landscape, suggesting a higher vulnerability to human interference within the landscape matrix for the populations in less productive habitats. Our findings suggest that the effect of human impact on genetic diversity may not be straightforward but modulated by the ecological conditions underlying local demographic processes. Populations within both landscapes seem to be vulnerable to loss of genetic diversity, but due to different mechanisms. This has implications for the choice of relevant management actions, that is, increasing population stability may be more relevant within an agricultural landscape still permeable for dispersal, while conserving dispersal corridors may be more appropriate in the forested landscape, to avoid isolation and increased genetic drift.

Pool Age and Local Habitat Are Associated with Effective Number of Breeders in Spotted Salamander Populations Colonizing Created Vernal Pools

Population genetics can reveal whether colonization of created habitats has been successful and inform future strategies for habitat creation. We used genetic analysis to investigate spotted salamander (Ambystoma maculatum) colonization of created vernal pools and explored the impact of habitat characteristics on the genetic diversity and connectivity of the pools. Our first objective was to examine genetic structure, differentiation, diversity, and potential for a founder effect. Our second objective was to determine if habitat characteristics were associated with effective number of breeders, relatedness, or genetic diversity. We sampled spotted salamander larvae in 31 created vernal pools (1–5 years old) in Monongahela National Forest (WV) in May and June 2015 and 2016. The youngest pools exhibited genetic differentiation, a founder effect, and low effective number of breeders. Effective number of breeders was positively associated with pool age, vegetation cover, pool diameter, and sample size. Vegetation cover was also negatively associated with relatedness. Genetic diversity did not have strong environmental predictors. Our results indicated the effective number of breeders increased and genetic differentiation decreased within 4–5 years of pool creation, a sign of rapid colonization and potential population establishment. Our research also showed that higher vegetative cover within the pool and larger pool diameters could impact habitat quality and should be incorporated into future pool creation.

Genetic Drift Shapes the Evolution of a Highly Dynamic Metapopulation

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

Population connectivity in voles (Microtus sp.) as a gauge for tall grass prairie restoration in midwestern North America

Ecological restoration can promote biodiversity conservation in anthropogenically fragmented habitats, but effectiveness of these management efforts need to be statistically validated to determine ’success.’ One such approach is to gauge the extent of recolonization as a measure of landscape permeability and, in turn, population connectivity. In this context, we estimated dispersal and population connectivity in prairie vole (Microtus ochrogaster; N = 231) and meadow vole (M. pennsylvanicus; N = 83) within five tall-grass prairie restoration sites embedded within the agricultural matrix of midwestern North America. We predicted that vole dispersal would be constrained by the extent of agricultural land surrounding restored habitat patches, spatially isolating vole populations and resulting in significant genetic structure. We first employed genetic assignment tests based on 15 microsatellite DNA loci to validate field-derived species-designations, then tested reclassified samples with multivariate and Bayesian clustering to assay for spatial and temporal genetic structure. Population connectivity was further evaluated by calculating pairwise F_ST, then potential demographic effects explored by computing migration rates, effective population size (N_e), and average relatedness (r). Genetic species assignments reclassified 25% of initial field identifications (N = 11 M. ochrogaster; N = 67 M. pennsylvanicus). In M. ochrogaster population connectivity was high across the study area, reflected in little to no spatial or temporal genetic structure. In M. pennsylvanicus genetic structure was detected, but relatedness estimates identified it as kin-clustering instead, underscoring social behavior among populations rather than spatial isolation as the cause. Estimates of N_e and r were stable across years, reflecting high dispersal and demographic resilience. Combined, these metrics suggest the agricultural matrix is highly permeable for voles and does not impede dispersal. High connectivity observed confirms that the restored landscape is productive and permeable for specific management targets such as voles and also demonstrates population genetic assays as a tool to statistically evaluate effectiveness of conservation initiatives.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean F_ST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

Landscape genetics of northern crested newt Triturus cristatus populations in a contrasting natural and human-impacted boreal forest

Among vertebrates, amphibians currently have the highest proportion of threatened species worldwide, mainly through loss of habitat, leading to increased population isolation. Smaller amphibian populations may lose more genetic diversity, and become more dependent on immigration for survival. Investigations of landscape factors and patterns mediating migration and population genetic differentiation are fundamental for knowledge-based conservation. The pond-breeding northern crested newt (Triturus cristatus) populations are decreasing throughout Europe, and are a conservation concern. Using microsatellites, we studied the genetic structure of the northern crested newt in a boreal forest ecosystem containing two contrasting landscapes, one subject to recent change and habitat loss by clear-cutting and roadbuilding, and one with little anthropogenic disturbance. Newts from 12 breeding ponds were analyzed for 13 microsatellites and 7 landscape and spatial variables. With a Maximum-likelihood population-effects model we investigated important landscape factors potentially explaining genetic patterns. Results indicate that intervening landscape factors between breeding ponds, explain the genetic differentiation in addition to an isolation-by-distance effect. Geographic distance, gravel roads, and south/south-west facing slopes reduced landscape permeability and increased genetic differentiation for these newts. The effect was opposite for streams, presumably being more favorable for newt dispersal. Populations within or bordering on old growth forest had a higher allelic richness than populations in managed forest outside these areas. Old growth forest areas may be important source habitats in the conservation of northern crested newt populations.

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.

Genetic variability and structure of an isolated population of Ambystoma altamirani, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (F_ST = 0.005) and high levels of genetic variability (H_o = 0.883; H_e = 0.621); we also found a small population size (N_e = 8.8), the presence of historical (M = 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (F_IS = -0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.

River network architecture, genetic effective size and distributional patterns predict differences in genetic structure across species in a dryland stream fish community

Dendritic ecological network (DEN) architecture can be a strong predictor of spatial genetic patterns in theoretical and simulation studies. Yet, interspecific differences in dispersal capabilities and distribution within the network may equally affect species' genetic structuring. We characterized patterns of genetic variation from up to ten microsatellite loci for nine numerically dominant members of the upper Gila River fish community, New Mexico, USA. Using comparative landscape genetics, we evaluated the role of network architecture for structuring populations within species (pairwise F_ST ) while explicitly accounting for intraspecific demographic influences on effective population size (N_e ). Five species exhibited patterns of connectivity and/or genetic diversity gradients that were predicted by network structure. These species were generally considered to be small-bodied or habitat specialists. Spatial variation of N_e was a strong predictor of pairwise F_ST for two species, suggesting patterns of connectivity may also be influenced by genetic drift independent of network properties. Finally, two study species exhibited genetic patterns that were unexplained by network properties and appeared to be related to nonequilibrium processes. Properties of DENs shape community-wide genetic structure but effects are modified by intrinsic traits and nonequilibrium processes. Further theoretical development of the DEN framework should account for such cases.

Applying gene flow science to environmental policy needs: a boundary work perspective

One application of gene flow science is the policy arena. In this article, we describe two examples in which the topic of gene flow has entered into the U.S. national environmental policymaking process: regulation of genetically engineered crops and clarification of the jurisdictional scope of the Clean Water Act. We summarize both current scientific understanding and the legal context within which gene flow science has relevance. We also discuss the process by which scientific knowledge has been synthesized and communicated to decision-makers in these two contexts utilizing the concept of 'boundary work'. Boundary organizations, the work they engage in to bridge the worlds of science, policy, and practice, and the boundary objects they produce to translate scientific knowledge existed in both examples. However, the specific activities and attributes of the objects produced varied based on the needs of the decision-makers. We close with suggestions for how scientists can contribute to or engage in boundary work with policymakers.

The Value of Learning about Natural History in Biodiversity Markets

Markets for biodiversity have generated much controversy because of the often unstated and untested assumptions included in transactions rules. Simple trading rules are favored to reduce transaction costs, but others have argued that this leads to markets that favor development and erode biodiversity. Here, I describe how embracing complexity and uncertainty within a tradable credit system for the Red-cockaded Woodpecker (Picoides borealis) creates opportunities to achieve financial and conservation goals simultaneously. Reversing the effects of habitat fragmentation is one of the main reasons for developing markets. I include uncertainty in habitat fragmentation effects by evaluating market transactions using five alternative dispersal models that were able to approximate observed patterns of occupancy and movement. Further, because dispersal habitat is often not included in market transactions, I contrast how changes in breeding versus dispersal habitat affect credit values. I use an individually-based, spatially-explicit population model for the Red-cockaded Woodpecker (Picoides borealis) to predict spatial- and temporal- influences of landscape change on species occurrence and genetic diversity. Results indicated that the probability of no net loss of abundance and genetic diversity responded differently to the transient dynamics in breeding and dispersal habitat. Trades that do not violate the abundance cap may simultaneously violate the cap for the erosion of genetic diversity. To highlight how economic incentives may help reduce uncertainty, I demonstrate tradeoffs between the value of tradable credits and the value of information needed to predict the influence of habitat trades on population viability. For the trade with the greatest uncertainty regarding the change in habitat fragmentation, I estimate that the value of using 13-years of data to reduce uncertainty in dispersal behaviors is $6.2 million. Future guidance for biodiversity markets should at least encourage the use of spatially- and temporally-explicit techniques that include population genetic estimates and the influence of uncertainty.

Rapid genetic restoration of a keystone species exhibiting delayed demographic response

Genetic founder effects are often expected when animals colonize restored habitat in fragmented landscapes, but empirical data on genetic responses to restoration are limited. We examined the genetic response of banner-tailed kangaroo rats (Dipodomys spectabilis) to landscape-scale grassland restoration in the Chihuahuan Desert of New Mexico, USA. Dipodomys spectabilis is a grassland specialist and keystone species. At sites treated with herbicide to remove shrubs, colonization by D. spectabilis is slow and populations persist at low density for ≥10 years (≥6 generations). Persistence at low density and low gene flow may cause strong founder effects. We compared genetic structure of D. spectabilis populations between treated sites and remnant grasslands, and we examined how the genetic response to restoration depended on treatment age, area, and connectivity to source populations. Allelic richness and heterozygosity were similar between treated sites and remnant grasslands. Allelic richness at treated sites was greatest early in the restoration trajectory, and genetic divergence did not differ between recently colonized and established populations. These results indicated that founder effects during colonization of treated sites were weak or absent. Moreover, our results suggested founder effects were not mitigated by treatment area or connectivity. Dispersal is negatively density-dependent in D. spectabilis, and we hypothesize that high gene flow may occur early in the restoration trajectory when density is low. Our study shows genetic diversity can be recovered more rapidly than demographic components of populations after habitat restoration and that founder effects are not inevitable for animals colonizing restored habitat in fragmented landscapes.

How spatio-temporal habitat connectivity affects amphibian genetic structure

Heterogeneous landscapes and fluctuating environmental conditions can affect species dispersal, population genetics, and genetic structure, yet understanding how biotic and abiotic factors affect population dynamics in a fluctuating environment is critical for species management. We evaluated how spatio-temporal habitat connectivity influences dispersal and genetic structure in a population of boreal chorus frogs (Pseudacris maculata) using a landscape genetics approach. We developed gravity models to assess the contribution of various factors to the observed genetic distance as a measure of functional connectivity. We selected (a) wetland (within-site) and (b) landscape matrix (between-site) characteristics; and (c) wetland connectivity metrics using a unique methodology. Specifically, we developed three networks that quantify wetland connectivity based on: (i) P. maculata dispersal ability, (ii) temporal variation in wetland quality, and (iii) contribution of wetland stepping-stones to frog dispersal. We examined 18 wetlands in Colorado, and quantified 12 microsatellite loci from 322 individual frogs. We found that genetic connectivity was related to topographic complexity, within- and between-wetland differences in moisture, and wetland functional connectivity as contributed by stepping-stone wetlands. Our results highlight the role that dynamic environmental factors have on dispersal-limited species and illustrate how complex asynchronous interactions contribute to the structure of spatially-explicit metapopulations.

Determinants of genetic structure in a nonequilibrium metapopulation of the plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F_ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F_ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.

Drought, deluge and declines: the impact of precipitation extremes on amphibians in a changing climate

The Class Amphibia is one of the most severely impacted taxa in an on-going global biodiversity crisis. Because amphibian reproduction is tightly associated with the presence of water, climatic changes that affect water availability pose a particularly menacing threat to both aquatic and terrestrial-breeding amphibians. We explore the impacts that one facet of climate change—that of extreme variation in precipitation—may have on amphibians. This variation is manifested principally as increases in the incidence and severity of both drought and major storm events. We stress the need to consider not only total precipitation amounts but also the pattern and timing of rainfall events. Such rainfall “pulses” are likely to become increasingly more influential on amphibians, especially in relation to seasonal reproduction. Changes in reproductive phenology can strongly influence the outcome of competitive and predatory interactions, thus potentially altering community dynamics in assemblages of co-existing species. We present a conceptual model to illustrate possible landscape and metapopulation consequences of alternative climate change scenarios for pond-breeding amphibians, using the Mole Salamander, Ambystoma talpoideum, as an example. Although amphibians have evolved a variety of life history strategies that enable them to cope with environmental uncertainty, it is unclear whether adaptations can keep pace with the escalating rate of climate change. Climate change, especially in combination with other stressors, is a daunting challenge for the persistence of amphibians and, thus, the conservation of global biodiversity.