Population structure of Columbia spotted frogs (Rana luteiventris) is strongly affected by the landscape

The landscape genetics of yellow perch (Perca flavescens) in a large fluvial ecosystem

Landscape genetics is being increasingly applied to elucidate the role of environmental features on the population structure of terrestrial organisms. However, the potential of this framework has been little explored in aquatic ecosystems such as large rivers. Here, we used a landscape genetics approach in order to (i) document the population structure of the yellow perch (Perca flavescens) by means of genetic variation at microsatellite markers, (ii) assess to what extent the structure was explained by landscape heterogeneity, and (iii) interpret the relevance of interactions between genetics and landscape for management and conservation. Analysis of the genetic variation among 1715 individuals from 16 localities and distributed over 310 km in the freshwater section of the Saint Lawrence River (Québec, Canada) revealed a relatively modest level of genetic structuring (F(ST) = 0.039). Application of the Monmonier's algorithm combining geographical and genetic information identified three zones of restricted gene flow defining four distinct populations. Physical barriers played a more important role on gene flow and genetic structure than waterway geographical distance. We found correlations between genetic differentiation and presence of distinct water masses in the sector of Lake Saint-Louis (r = 0.7177, P = 0.0340) and with fragmentation of spawning habitats in the sector of Lake Saint-Pierre (r = 0.8578, P = 0.0095). Our results support the treatment of four distinct biological units, which is in contrast with the current basis for yellow perch management. Finally, this study showed that landscape genetics is a powerful means to identify environmental barriers to gene flow causing genetic discontinuities in apparently highly connected aquatic landscapes.

Postglacial range expansion from northern refugia by the wood frog, Rana sylvatica

Although the range dynamics of North American amphibians during the last glacial cycle are increasingly better understood, the recolonization history of the most northern regions and the impact of southern refugia on patterns of intraspecific genetic diversity and phenotypic variation in these regions are not well reconstructed. Here we present the phylogeographic history of a widespread and primarily northern frog, Rana sylvatica. We surveyed 551 individuals from 116 localities across the species' range for a 650-bp region of the NADH dehydrogenase subunit 2 and tRNA(TRP) mitochondrial genes. Our phylogenetic analyses revealed two distinct clades corresponding to eastern and western populations, as well as a Maritime subclade within the eastern lineage. Patterns of genetic diversity support multiple refugia. However, high-latitude refugia in the Appalachian highlands and modern-day Wisconsin appear to have had the biggest impact on northern populations. Clustering analyses based on morphology further support a distinction between eastern and western wood frogs and suggest that postglacial migration has played an important role in generating broad-scale patterns of phenotypic variation in this species.

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.

GEOLOGICAL AND CLIMATIC FORCES DRIVING SPECIATION IN THE CONTINENTALLY DISTRIBUTED TRILLING CHORUS FROGS (PSEUDACRIS)

Tertiary geological events and Quaternary climatic fluctuations have been proposed as important factors of speciation in the North American flora and fauna. Few studies, however, have rigorously tested hypotheses regarding the specific factors driving divergence of taxa. Here, we test explicit speciation hypotheses by correlating geologic events with divergence times among species in the continentally distributed trilling chorus frogs (Pseudacris). In particular, we ask whether marine inundation of the Mississippi Embayment, uplift of the Appalachian Mountains, or modification of the ancient Teays-Mahomet River system contributed to speciation. To examine the plausibility of ancient rivers causing divergence, we tested whether modern river systems inhibit gene flow. Additionally, we compared the effects of Quaternary climatic factors (glaciation and aridification) on levels of genetic variation. Divergence time estimates using penalized likelihood and coalescent approaches indicate that the major lineages of chorus frogs diversified during the Tertiary, and also exclude Quaternary climate change as a factor in speciation of chorus frogs. We show the first evidence that inundation of the Mississippi Embayment contributed to speciation. We reject the hypotheses that Cenozoic uplift of the Appalachians and that diversion of the Teays-Mahomet River contributed to speciation in this clade. We find that by reducing gene flow, rivers have the potential to cause divergence of lineages. Finally, we demonstrate that populations in areas affected by Quaternary glaciation and aridification have reduced levels of genetic variation compared to those from more equable regions, suggesting recent colonization.

Applying new inter-individual approaches to assess fine-scale population genetic diversity in a neotropical frog, Eleutherodactylus ockendeni

We assess patterns of genetic diversity of a neotropical leaflitter frog, Eleutherodactylus ockendeni, in the upper Amazon of Ecuador without a priori delineation of biological populations and with sufficiently intensive sampling to assess inter-individual patterns. We mapped the location of each collected frog across a 5.4 x 1 km landscape at the Jatun Sacha Biological Station, genotyped 185 individuals using five species-specific DNA microsatellite loci, and sequenced a fragment of mitochondrial cytochrome b for a subset of 51 individuals. The microsatellites were characterized by high allelic diversity and homozygote excess across all loci, suggesting that when pooled the sample is not a panmictic population. We conclude that the lack of panmixia is not attributable to the influence of null alleles or biased sampling of consanguineous family groups. Multiple methods of population cluster analysis, using both Bayesian and maximum likelihood approaches, failed to identify discrete genetic clusters across the sampled area. Using multivariate spatial autocorrelation, kinship coefficients and relatedness coefficients, we identify a continuous isolation by distance population structure, with a first patch size of ca. 260 m and apparently large population sizes. Analysis of mtDNA corroborates the observation of high genetic diversity at fine scales: there are multiple haplotypes, they are non-randomly distributed and a binary haplotype correlogram shows significant spatial genetic autocorrelation. We demonstrate the utility of inter-individual genetic methods and caution against making a priori assumptions about population genetic structure based simply on arbitrary or convenient patterns of sampling.

Responses of pond-breeding amphibians to wildfire: short-term patterns in occupancy and colonization

Wildland fires are expected to become more frequent and severe in many ecosystems, potentially posing a threat to many sensitive species. We evaluated the effects of a large, stand-replacement wildfire on three species of pond-breeding amphibians by estimating changes in occupancy of breeding sites during the three years before and after the fire burned 42 of 83 previously surveyed wetlands. Annual occupancy and colonization for each species was estimated using recently developed models that incorporate detection probabilities to provide unbiased parameter estimates. We did not find negative effects of the fire on the occupancy or colonization rates of the long-toed salamander (Ambystoma macrodactylum). Instead, its occupancy was higher across the study area after the fire, possibly in response to a large snowpack that may have facilitated colonization of unoccupied wetlands. Naive data (uncorrected for detection probability) for the Columbia spotted frog (Rana luteiventris) initially led to the conclusion of increased occupancy and colonization in wetlands that burned. After accounting for temporal and spatial variation in detection probabilities, however, it was evident that these parameters were relatively stable in both areas before and after the fire. We found a similar discrepancy between naive and estimated occupancy of A. macrodactylum that resulted from different detection probabilities in burned and control wetlands. The boreal toad (Bufo boreas) was not found breeding in the area prior to the fire but colonized several wetlands the year after they burned. Occupancy by B. boreas then declined during years 2 and 3 following the fire. Our study suggests that the amphibian populations we studied are resistant to wildfire and that B. boreas may experience short-term benefits from wildfire. Our data also illustrate how naive presence-non-detection data can provide misleading results.

The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

A primary goal of molecular ecology is to understand the influence of abiotic factors on the spatial distribution of genetic variation. Features including altitudinal clines, topography and landscape characteristics affect the proportion of suitable habitat, influence dispersal patterns, and ultimately structure genetic differentiation among populations. We studied the effects of altitude and topography on genetic variation of long-toed salamanders (Ambystoma macrodactylum), a geographically widespread amphibian species throughout northwestern North America. We focused on 10 low altitude sites (< 1200 m) and 11 high-altitude sites in northwestern Montana and determined multilocus genotypes for 549 individuals using seven microsatellite loci. We tested four hypotheses: (1) gene flow is limited between high- and low-altitude sites; and, (2) gene flow is limited among high-altitude sites due to harsh habitat and extreme topographical relief between sites; (3) low-altitude sites exhibit higher among-site gene flow due to frequent flooding events and low altitudinal relief; and (4) there is a negative correlation between altitude and genetic variation. Overall F(ST) values were moderate (0.08611; P < 0.001). Pairwise F(ST) estimates between high and low populations and a population graphing method supported the hypothesis that low-altitude and high-altitude sites, taken together, are genetically differentiated from each other. Also as predicted, gene flow is more prominent among low-altitude sites than high-altitude sites; low-altitude sites had a significantly lower F(ST) (0.03995; P < 0.001) than high altitude sites (F(ST) = 0.10271; P < 0.001). Use of Bayesian analysis of population structure (BAPS) resulted in delineation of 10 genetic groups, two among low-altitude populations and eight among high-altitude populations. In addition, within high altitude populations, basin-level genetic structuring was apparent. A nonequilibrium algorithm for detecting current migration rates supported these population distinctions. Finally, we also found a significant negative correlation between genetic diversity and altitude. These results are consistent with the hypothesis that topography and altitudinal gradients shape the spatial distribution of genetic variation in a species with a broad geographical range and diverse life history. Our study sheds light on which key factors limit dispersal and ultimately species' distributions.

Genetic diversity, but not hatching success, is jointly affected by postglacial colonization and isolation in the threatened frog, Rana latastei

Both postglacial colonization and habitat fragmentation can reduce the genetic diversity of populations, which in turn can affect fitness. However, since these processes occur at different spatial and temporal scales, the consequences of either process may differ. To disentangle the relative role of isolation and postglacial colonization in determining genetic diversity and fitness, we studied microsatellite diversity of 295 individuals from 10 populations and measured the hatch rate of 218 clutches from eight populations of a threatened frog, R. latastei. The populations that were affected by fragmentation to a greater extent suffered higher embryo mortality and reduced hatch rate, while no effects of distance from glacial refugium on hatch rate were detected. Altogether, distance from glacial refugium and isolation explained > 90% of variation in genetic diversity. We found that the genetic diversity was lowest in populations both isolated and far from the glacial refugium, and that distance from refugium seems to have the primary role in determining genetic diversity. The relationship between genetic diversity and hatch rate was not significant. However, the proportion of genetic diversity lost through recent isolation had a significant, negative effect on fitness. It is possible that selection at least partially purged the negative effects of the ancestral loss of genetic diversity.

Fine-scale spatial genetic structure and dispersal among spotted salamander (Ambystoma maculatum) breeding populations

We examined fine-scale genetic variation among breeding aggregations of the spotted salamander (Ambystoma maculatum) to quantify dispersal, interpopulation connectivity and population genetic structure. Spotted salamanders rely on temporary ponds or wetlands for aggregate breeding. Adequate breeding sites are relatively isolated from one another and field studies suggest considerable adult site fidelity; therefore, we expected to find population structure and differentiation at small spatial scales. We used microsatellites to estimate population structure and dispersal among 29 breeding aggregations in Tompkins County, New York, USA, an area encompassing 1272 km(2). Bayesian and frequency-based analyses revealed fine-scale genetic structure with two genetically defined demes: the North deme included seven breeding ponds, and the South deme included 13 ponds. Nine ponds showed evidence of admixture between these two genetic pools. Bayesian assignment tests for detection of interpopulation dispersal indicate that immigration among ponds is common within demes, and that certain populations serve as sources of immigrants to neighbouring ponds. Likewise, spatial genetic correlation analyses showed that populations < or = 4.8 km distant from each other show significant genetic correlation that is not evident at higher scales. Within-population levels of relatedness are consistently larger than expected if mating were completely random across ponds, and in the case of a few ponds, within-population processes such as inbreeding or reproductive skew contribute significantly to differentiation from neighbouring ponds. Our data underscore the importance of these within-population processes as a source of genetic diversity across the landscape, despite considerable population connectivity. Our data further suggest that spotted salamander breeding groups behave as metapopulations, with population clusters as functional units, but sufficient migration among demes to allow for potential rescue and recolonization. Amphibian habitats are becoming increasingly fragmented and a clear understanding of dispersal and patterns of population connectivity for taxa with different ecologies and life histories is crucial for their conservation.

High dispersal in a frog species suggests that it is vulnerable to habitat fragmentation

Global losses of amphibian populations are a major conservation concern and their causes have generated substantial debate. Habitat fragmentation is considered one important cause of amphibian decline. However, if fragmentation is to be invoked as a mechanism of amphibian decline, it must first be established that dispersal is prevalent among contiguous amphibian populations using formal movement estimators. In contrast, if dispersal is naturally low in amphibians, fragmentation can be disregarded as a cause of amphibian declines and conservation efforts can be focused elsewhere. We examined dispersal rates in Columbia spotted frogs (Rana luteiventris) using capture-recapture analysis of over 10,000 frogs in combination with genetic analysis of microsatellite loci in replicate basins. We found that frogs had exceptionally high juvenile dispersal rates (up to 62% annually) over long distances (>5km), large elevation gains (>750m) and steep inclines (36 degrees incline over 2km) that were corroborated by genetic data showing high gene flow. These findings show that dispersal is an important life-history feature of some amphibians and suggest that habitat fragmentation is a serious threat to amphibian persistence.

Amphibians as Models for Studying Environmental Change

The use of amphibians as models in ecological research has a rich history. From an early foundation in studies of amphibian natural history sprang generations of scientists who used amphibians as models to address fundamental questions in population and community ecology. More recently, in the wake of an environment that human disturbances rapidly altered, ecologists have adopted amphibians as models for studying applied ecological issues such as habitat loss, pollution, disease, and global climate change. Some of the characteristics of amphibians that make them useful models for studying these environmental problems are highlighted, including their trophic importance, environmental sensitivity, research tractability, and impending extinction. The article provides specific examples from the recent literature to illustrate how studies on amphibians have been instrumental in guiding scientific thought on a broad scale. Included are examples of how amphibian research has transformed scientific disciplines, generated new theories about global health, called into question widely accepted scientific paradigms, and raised awareness in the general public that our daily actions may have widespread repercussions. In addition, studies on amphibian declines have provided insight into the complexity in which multiple independent factors may interact with one another to produce catastrophic and sometimes unpredictable effects. Because of the complexity of these problems, amphibian ecologists have been among the strongest advocates for interdisciplinary research. Future studies of amphibians will be important not only for their conservation but also for the conservation of other species, critical habitats, and entire ecosystems.

Putting the ‘landscape’ in landscape genetics

Landscape genetics has emerged as a new research area that integrates population genetics, landscape ecology and spatial statistics. Researchers in this field can combine the high resolution of genetic markers with spatial data and a variety of statistical methods to evaluate the role that landscape variables play in shaping genetic diversity and population structure. While interest in this research area is growing rapidly, our ability to fully utilize landscape data, test explicit hypotheses and truly integrate these diverse disciplines has lagged behind. Part of the current challenge in the development of the field of landscape genetics is bridging the communication and knowledge gap between these highly specific and technical disciplines. The goal of this review is to help bridge this gap by exposing geneticists to terminology, sampling methods and analysis techniques widely used in landscape ecology and spatial statistics but rarely addressed in the genetics literature. We offer a definition for the term "landscape genetics", provide an overview of the landscape genetics literature, give guidelines for appropriate sampling design and useful analysis techniques, and discuss future directions in the field. We hope, this review will stimulate increased dialog and enhance interdisciplinary collaborations advancing this exciting new field.

Fine-scale population differentiation and gene flow in a terrestrial salamander (Plethodon cinereus) living in continuous habitat

Several recent studies have shown that amphibian populations may exhibit high genetic subdivision in areas with recent fragmentation and urban development. Less is known about the potential for genetic differentiation in continuous habitats. We studied genetic differentiation of red-backed salamanders (Plethodon cinereus) across a 2-km transect through continuous forest in Virginia, USA. Mark-recapture studies suggest very little dispersal for this species, whereas homing experiments and post-Pleistocene range expansion both suggest greater dispersal abilities. We used six microsatellite loci to examine genetic population structure and differentiation between eight subpopulations of red-backed salamanders at distances from 200 m to 2 km. We also used several methods to extrapolate dispersal frequencies and test for sex-biased dispersal. We found small, but detectable differentiation among populations, even at distances as small as 200 m. Differentiation was closely correlated with distance and both Mantel tests and assignment tests were consistent with an isolation-by-distance model for the population. Extrapolations of intergenerational variance in spatial position (sigma(2)<15 m(2)) and pair-wise dispersal frequencies (4 Nm < 25 for plots separated by 300 m) both suggest limited gene flow. Additionally, tests for sex-biased dispersal imply that dispersal frequency is similarly low for both sexes. We suggest that these low levels of gene flow and the infrequent dispersal observed in mark-recapture studies may be reconciled with homing ability and range expansion if dispersing animals rarely succeed in breeding in saturated habitats, if dispersal is flexible depending on the availability of habitat, or if dispersal frequency varies across the geographic range of red-backed salamanders.

Genetic structure is influenced by landscape features: empirical evidence from a roe deer population

The delimitation of population units is of primary importance in population management and conservation biology. Moreover, when coupled with landscape data, the description of population genetic structure can provide valuable knowledge about the permeability of landscape features, which is often difficult to assess by direct methods (e.g. telemetry). In this study, we investigated the genetic structuring of a roe deer population which recently recolonized a fragmented landscape. We sampled 1148 individuals from a 40 x 55-km area containing several putative barriers to deer movements, and hence to gene flow, namely a highway, rivers and several canals. In order to assess the effect of these landscape features on genetic structure, we implemented a spatial statistical model known as geneland which analyses genetic structure, explicitly taking into account the spatial nature of the problem. Two genetic units were inferred, exhibiting a very low level of differentiation (F(ST) = 0.008). The location of their boundaries suggested that there are no absolute barriers in this study area, but that the combination of several landscape features with low permeability can lead to population differentiation. Our analysis hence suggests that the landscape has a significant influence on the structuring of the population under study. It also illustrates the use of geneland as a powerful method to infer population structure, even in situations of young populations exhibiting low genetic differentiation.

Genetic structure in northeastern populations of the Alpine newt (Triturus alpestris): evidence for post-Pleistocene differentiation

Genetic variation in 13 populations of the Alpine newt, Triturus alpestris, was assessed at the northeastern margin of its range (southern Poland). Variation at six microsatellite loci was scored in 354 newts, and two mitochondrial DNA fragments (c. 2000 bp) were sequenced in a subset of 27 individuals. Significant differences in allele frequencies and the presence of private alleles determined genetic units corresponding to three separate mountain ranges, i.e. the Carpathian, Sudetes and Holy Cross Mountains. F(ST)'s were three times greater in among than in within mountain range pairwise comparisons. An assignment test and pairwise F(ST)'s suggested relatively high levels of gene flow at the local level, although the Sudetes populations revealed some subtle structuring. Genetic variation was lower in the Carpathians and Holy Cross Mountains. The geographic pattern of mitochondrial DNA variation indicated that these newt populations originated from a single glacial refugium/founder population, and that the colonization of southern Poland took place in an easterly direction. The data show that substantial neutral variation and between group divergence has accumulated relatively quickly in these low-vagility organisms. The Alpine newt case exemplifies species history as a factor determining patterns of genetic diversity in marginal populations.

Gene flow and functional connectivity in the natterjack toad

Functional connectivity is a key factor for the persistence of many specialist species in fragmented landscapes. However, connectivity estimates have rarely been validated by the observation of dispersal movements. In this study, we estimated functional connectivity of a real landscape by modelling dispersal for the endangered natterjack toad (Bufo calamita) using cost distance. Cost distance allows the evaluation of 'effective distances', which are distances corrected for the costs involved in moving between habitat patches in spatially explicit landscapes. We parameterized cost-distance models using the results of our previous experimental investigation of natterjack's movement behaviour. These model predictions (connectivity estimates from the GIS study) were then confronted to genetic-based dispersal rates between natterjack populations in the same landscape using Mantel tests. Dispersal rates between the populations were inferred from variation at six microsatellite loci. Based on these results, we conclude that matrix structure has a strong effect on dispersal rates. Moreover, we found that cost distances generated by habitat preferences explained dispersal rates better than did the Euclidian distances, or the connectivity estimate based on patch-specific resistances (patch viscosity). This study is a clear example of how landscape genetics can validate operational functional connectivity estimates.

Geographic patterns of microsatellite variation in Boechera stricta, a close relative of Arabidopsis

The genus Boechera is a widespread North American group with great potential for studies of ecology and evolution: Boechera is closely related to Arabidopsis and exhibits different ecological and reproductive strategies. Boechera stricta (previously Arabis drummondii) is a morphologically and genetically well-defined, perennial crucifer species. Fifteen natural populations of diploid individuals from the Rocky Mountains were analysed using 21 microsatellite loci. In accordance with our expectation for this predominately inbreeding species, a high F IS value (0.89) was observed. Furthermore, populations of B. stricta were highly differentiated, as indicated by F ST = 0.56. Three clusters were identified using structure- the majority of populations belonged to either the Northern or Southern cluster. Together, the north-south partitioning and evenness of genetic variation across the two clusters suggested multiple refugia for this perennial herb in the Rocky Mountains. Pleistocene glaciation, together with the topographically and climatologically heterogeneous cordillera, has profoundly influenced the genetic architecture of B. stricta. Genetic population structure was also influenced by relatively recent genome admixture at two levels: within species (involving individuals from the Northern and Southern clusters) and between species (with the hybridization of B. stricta and Boechera holboellii). This complexity of population structure at presumably neutral microsatellite loci located throughout the genome in B. stricta provides a baseline against which to test whether functional genetic variation is undergoing local adaptive evolution throughout the natural species range.

Landscape genetics of the blotched tiger salamander (Ambystoma tigrinum melanostictum)

The field of landscape genetics has great potential to identify habitat features that influence population genetic structure. To identify landscape correlates of genetic differentiation in a quantitative fashion, we developed a novel approach using geographical information systems analysis. We present data on blotched tiger salamanders (Ambystoma tigrinum melanostictum) from 10 sites across the northern range of Yellowstone National Park in Montana and Wyoming, USA. We used eight microsatellite loci to analyse population genetic structure. We tested whether landscape variables, including topographical distance, elevation, wetland likelihood, cover type and number of river and stream crossings, were correlated with genetic subdivision (F(ST)). We then compared five hypothetical dispersal routes with a straight-line distance model using two approaches: (i) partial Mantel tests using Akaike's information criterion scores to evaluate model robustness and (ii) the BIOENV procedure, which uses a Spearman rank correlation to determine the combination of environmental variables that best fits the genetic data. Overall, gene flow appears highly restricted among sites, with a global F(ST) of 0.24. While there is a significant isolation-by-distance pattern, incorporating landscape variables substantially improved the fit of the model (from an r2 of 0.3 to 0.8) explaining genetic differentiation. It appears that gene flow follows a straight-line topographic route, with river crossings and open shrub habitat correlated with lower F(ST) and thus, decreased differentiation, while distance and elevation difference appear to increase differentiation. This study demonstrates a general approach that can be used to determine the influence of landscape variables on population genetic structure.

Conservation genetics of amphibians

Amphibians are good models for investigating the genetics of wild animal populations because they are: (1) widely distributed in most ecosystems; (2) easy to sample in breeding assemblages; (3) often philopatric to breeding sites, generating high levels of population genetic structure; (4) amenable to controlled crossings in the laboratory; and (5) of major conservation concern. Neutral genetic markers, mostly microsatellites, have been used successfully in studies of amphibian effective population sizes and structures, and in assessing the consequences of hybridisation. Phylogeography has provided important insights into population histories and the fates of introductions. Quantitative genetic methods have demonstrated adaptive variation in life history traits of importance to fitness and therefore to population viability.