Melanie A. Murphy

Landscape genetics: where are we now?

Landscape genetics has seen rapid growth in number of publications since the term was coined in 2003. An extensive literature search from 1998 to 2008 using keywords associated with landscape genetics yielded 655 articles encompassing a vast array of study organisms, study designs and methodology. These publications were screened to identify 174 studies that explicitly incorporated at least one landscape variable with genetic data. We systematically reviewed this set of papers to assess taxonomic and temporal trends in: (i) geographic regions studied; (ii) types of questions addressed; (iii) molecular markers used; (iv) statistical analyses used; and (v) types and nature of spatial data used. Overall, studies have occurred in geographic regions proximal to developed countries and more commonly in terrestrial vs. aquatic habitats. Questions most often focused on effects of barriers and/or landscape variables on gene flow. The most commonly used molecular markers were microsatellites and amplified fragment length polymorphism (AFLPs), with AFLPs used more frequently in plants than animals. Analysis methods were dominated by Mantel and assignment tests. We also assessed differences among journals to evaluate the uniformity of reporting and publication standards. Few studies presented an explicit study design or explicit descriptions of spatial extent. While some landscape variables such as topographic relief affected most species studied, effects were not universal, and some species appeared unaffected by the landscape. Effects of habitat fragmentation were mixed, with some species altering movement paths and others unaffected. Taken together, although some generalities emerged regarding effects of specific landscape variables, results varied, thereby reinforcing the need for species‐specific work. We conclude by: highlighting gaps in knowledge and methodology, providing guidelines to authors and reviewers of landscape genetics studies, and suggesting promising future directions of inquiry.

Quantifying the lag time to detect barriers in landscape genetics

Understanding how spatial genetic patterns respond to landscape change is crucial for advancing the emerging field of landscape genetics. We quantified the number of generations for new landscape barrier signatures to become detectable and for old signatures to disappear after barrier removal. We used spatially explicit, individual‐based simulations to examine the ability of an individual‐based statistic [Mantel’s r using the proportion of shared alleles’ statistic (Dps)] and population‐based statistic (FST) to detect barriers. We simulated a range of movement strategies including nearest neighbour dispersal, long‐distance dispersal and panmixia. The lag time for the signal of a new barrier to become established is short using Mantel’s r (1–15 generations). FST required approximately 200 generations to reach 50% of its equilibrium maximum, although G’ST performed much like Mantel’s r. In strong contrast, FST and Mantel’s r perform similarly following the removal of a barrier formerly dividing a population. Also, given neighbour mating and very short‐distance dispersal strategies, historical discontinuities from more than 100 generations ago might still be detectable with either method. This suggests that historical events and landscapes could have long‐term effects that confound inferences about the impacts of current landscape features on gene flow for species with very little long‐distance dispersal. Nonetheless, populations of organisms with relatively large dispersal distances will lose the signal of a former barrier within less than 15 generations, suggesting that individual‐based landscape genetic approaches can improve our ability to measure effects of existing landscape features on genetic structure and connectivity.

Utility of computer simulations in landscape genetics

Population genetics theory is primarily based on mathematical models in which spatial complexity and temporal variability are largely ignored. In contrast, the field of landscape genetics expressly focuses on how population genetic processes are affected by complex spatial and temporal environmental heterogeneity. It is spatially explicit and relates patterns to processes by combining complex and realistic life histories, behaviours, landscape features and genetic data. Central to landscape genetics is the connection of spatial patterns of genetic variation to the usually highly stochastic space–time processes that create them over both historical and contemporary time periods. The field should benefit from a shift to computer simulation approaches, which enable incorporation of demographic and environmental stochasticity. A key role of simulations is to show how demographic processes such as dispersal or reproduction interact with landscape features to affect probability of site occupancy, population size, and gene flow, which in turn determine spatial genetic structure. Simulations could also be used to compare various statistical methods and determine which have correct type I error or the highest statistical power to correctly identify spatio‐temporal and environmental effects. Simulations may also help in evaluating how specific spatial metrics may be used to project future genetic trends. This article summarizes some of the fundamental aspects of spatial–temporal population genetic processes. It discusses the potential use of simulations to determine how various spatial metrics can be rigorously employed to identify features of interest, including contrasting locus‐specific spatial patterns due to micro‐scale environmental selection.

Landscape genetics of high mountain frog metapopulations

Explaining functional connectivity among occupied habitats is crucial for understanding metapopulation dynamics and species ecology. Landscape genetics has primarily focused on elucidating how ecological features between observations influence gene flow. Functional connectivity, however, may be the result of both these between‐site (landscape resistance) landscape characteristics and at‐site (patch quality) landscape processes that can be captured using network based models. We test hypotheses of functional connectivity that include both between‐site and at‐site landscape processes in metapopulations of Columbia spotted frogs (Rana luteiventris) by employing a novel justification of gravity models for landscape genetics (eight microsatellite loci, 37 sites, n = 441). Primarily used in transportation and economic geography, gravity models are a unique approach as flow (e.g. gene flow) is explained as a function of three basic components: distance between sites, production/attraction (e.g. at‐site landscape process) and resistance (e.g. between‐site landscape process). The study system contains a network of nutrient poor high mountain lakes where we hypothesized a short growing season and complex topography between sites limit R. luteiventris gene flow. In addition, we hypothesized production of offspring is limited by breeding site characteristics such as the introduction of predatory fish and inherent site productivity. We found that R. luteiventris connectivity was negatively correlated with distance between sites, presence of predatory fish (at‐site) and topographic complexity (between‐site). Conversely, site productivity (as measured by heat load index, at‐site) and growing season (as measured by frost‐free period between‐sites) were positively correlated with gene flow. The negative effect of predation and positive effect of site productivity, in concert with bottleneck tests, support the presence of source–sink dynamics. In conclusion, gravity models provide a powerful new modelling approach for examining a wide range of both basic and applied questions in landscape genetics.

Critical considerations for the application of environmental DNA methods to detect aquatic species

Summary Species detection using environmental DNA (eDNA) has tremendous potential for contributing to the understanding of the ecology and conservation of aquatic species. Detecting species using eDNA methods, rather than directly sampling the organisms, can reduce impacts on sensitive species and increase the power of field surveys for rare and elusive species. The sensitivity of eDNA methods, however, requires a heightened awareness and attention to quality assurance and quality control protocols. Additionally, the interpretation of eDNA data demands careful consideration of multiple factors. As eDNA methods have grown in application, diverse approaches have been implemented to address these issues. With interest in eDNA continuing to expand, supportive guidelines for undertaking eDNA studies are greatly needed. Environmental DNA researchers from around the world have collaborated to produce this set of guidelines and considerations for implementing eDNA methods to detect aquatic macroorganisms. Critical considerations for study design include preventing contamination in the field and the laboratory, choosing appropriate sample analysis methods, validating assays, testing for sample inhibition and following minimum reporting guidelines. Critical considerations for inference include temporal and spatial processes, limits of correlation of eDNA with abundance, uncertainty of positive and negative results, and potential sources of allochthonous DNA. We present a synthesis of knowledge at this stage for application of this new and powerful detection method.

A simulation-based evaluation of methods for inferring linear barriers to gene flow

Different analytical techniques used on the same data set may lead to different conclusions about the existence and strength of genetic structure. Therefore, reliable interpretation of the results from different methods depends on the efficacy and reliability of different statistical methods. In this paper, we evaluated the performance of multiple analytical methods to detect the presence of a linear barrier dividing populations. We were specifically interested in determining if simulation conditions, such as dispersal ability and genetic equilibrium, affect the power of different analytical methods for detecting barriers. We evaluated two boundary detection methods (Monmonier’s algorithm and WOMBLING), two spatial Bayesian clustering methods (TESS and GENELAND), an aspatial clustering approach (STRUCTURE), and two recently developed, non‐Bayesian clustering methods [PSMIX and discriminant analysis of principal components (DAPC)]. We found that clustering methods had higher success rates than boundary detection methods and also detected the barrier more quickly. All methods detected the barrier more quickly when dispersal was long distance in comparison to short‐distance dispersal scenarios. Bayesian clustering methods performed best overall, both in terms of highest success rates and lowest time to barrier detection, with GENELAND showing the highest power. None of the methods suggested a continuous linear barrier when the data were generated under an isolation‐by‐distance (IBD) model. However, the clustering methods had higher potential for leading to incorrect barrier inferences under IBD unless strict criteria for successful barrier detection were implemented. Based on our findings and those of previous simulation studies, we discuss the utility of different methods for detecting linear barriers to gene flow.

Relationships between migration rates and landscape resistance assessed using individual-based simulations

Linking landscape effects on gene flow to processes such as dispersal and mating is essential to provide a conceptual foundation for landscape genetics. It is particularly important to determine how classical population genetic models relate to recent individual‐based landscape genetic models when assessing individual movement and its influence on population genetic structure. We used classical Wright–Fisher models and spatially explicit, individual‐based, landscape genetic models to simulate gene flow via dispersal and mating in a series of landscapes representing two patches of habitat separated by a barrier. We developed a mathematical formula that predicts the relationship between barrier strength (i.e., permeability) and the migration rate (m) across the barrier, thereby linking spatially explicit landscape genetics to classical population genetics theory. We then assessed the reliability of the function by obtaining population genetics parameters (m, FST) using simulations for both spatially explicit and Wright–Fisher simulation models for a range of gene flow rates. Next, we show that relaxing some of the assumptions of the Wright–Fisher model can substantially change population substructure (i.e., FST). For example, isolation by distance among individuals on each side of a barrier maintains an FST of ∼0.20 regardless of migration rate across the barrier, whereas panmixia on each side of the barrier results in an FST that changes with m as predicted by classical population genetics theory. We suggest that individual‐based, spatially explicit modelling provides a general framework to investigate how interactions between movement and landscape resistance drive population genetic patterns and connectivity across complex landscapes.

High-resolution landcover classification using Random Forest

Potential data sets for landcover classification, such as Landsat (or pre-processed data such as the National Land Cover Dataset (NLCD)), are often too coarse for fine-scale research needs or are cost-prohibitive (Quickbird, Ikonos and Geoeye). Repeated attempts at classifying high spatial resolution data, National Agricultural Imagery Program (NAIP) imagery, based on traditional techniques, such as a maximum likelihood supervised classification, have failed to produce a product with sufficient accuracy. We used the ensemble Random Forests (RFs) classifier to classify landcover at 1 m resolution using 2009 NAIP imagery in south-eastern Wyoming. We classified riparian areas within a 225 km2 study area, at 1 m spatial resolution, using RFs with emphasis on riparian corridors that yielded a land cover map with overall accuracy of 81% and a kappa coefficient of 79%. Users’ accuracy of important riparian vegetation species, aspen, riparian grasses and willow were 79%, 81% and 83%, respectively. Techniques presented in this paper, which exploit free NAIP imagery for landcover classification, represent an inexpensive and reliable alternative to purchasing commercial imagery when high spatial resolution landcover data is needed.

Effects of hydroperiod duration on survival, developmental rate, and size at metamorphosis in boreal chorus frog tadpoles ( Pseudacris maculata )

Abstract: Understanding the relationship between climate-driven habitat conditions and survival is key to preserving biodiversity in the face of rapid climate change. Hydroperiod—the length of time water is in a wetland—is a critical limiting habitat variable for amphibians as larvae must metamorphose before ponds dry. Changes in precipitation and temperature patterns are affecting hydroperiod globally, but the impact of these changes on amphibian persistence is poorly understood. We studied the responses of Boreal Chorus Frog (Pseudacris maculata) tadpoles to simulated hydroperiods (i.e., water level reductions) in the laboratory using individuals collected from ponds spanning a range of natural hydroperiods (Colorado Front Range, USA). To assess the effects of experimental hydroperiod reduction, we measured mortality, time to metamorphosis, and size at metamorphosis. We found that tadpoles grew at rates reflecting the hydroperiods of their native ponds, regardless of experimental treatment. Tadpoles from permanent ponds metamorphosed faster than those from ephemeral ponds across all experimental treatments, a pattern which may represent a predation selection gradient or countergradient variation in developmental rates. Size at metamorphosis did not vary across experimental treatments. Mortality was low overall but varied with pond of origin. Our results suggest that adaptation to local hydroperiod and/or predation and temperature conditions is important in P. maculata. Moreover, the lack of a plastic response to reduced hydroperiods suggests that P. maculata may not be able to metamorphose quickly enough to escape drying ponds. These results have important implications for amphibian persistence in ponds predicted to dry more quickly due to rapid climate change.

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.