The sensitivity of genetic connectivity measures to unsampled and under-sampled sites

Limited effects of population age on the genetic structure of spatially isolated forest herb populations in temperate Europe

Due to multiple land-cover changes, forest herb populations residing in forest patches embedded in agricultural landscapes display different ages and, thus, experience differences in genetic exchange, mutation accumulation and genetic drift. The extent of divergence in present-day population genetic structure among these populations of different ages remains unclear, considering their diverse breeding systems and associated pollinators. Answering this question is essential to understand these species' persistence, maintenance of evolutionary potential and adaptability to changing environments. We applied a multi-landscape setup to compare the genetic structure of forest herb populations across forest patches of different ages (18-338 years). We studied the impact on three common slow-colonizer herb species with distinct breeding systems and associated pollinators: Polygonatum multiflorum (outcrossing, long-distance pollinators), Anemone nemorosa (outcrossing, short-distance pollinators) and Oxalis acetosella (mixed breeding). We aimed to assess if in general older populations displayed higher genetic diversity and lower differentiation than younger ones. We also anticipated that P. multiflorum would show the smallest while O. acetosella the largest difference, between old and young populations. We found that older populations had a higher observed heterozygosity (H o) but a similar level of allelic richness (A r) and expected heterozygosity (H e) as younger populations, except for A. nemorosa, which exhibited higher A r and H e in younger populations. As populations aged, their pairwise genetic differentiation measured by D PS decreased independent of species identity while the other two genetic differentiation measures showed either comparable levels between old and young populations (G" ST) or inconsistency among three species (cGD). The age difference of the two populations did not explain their genetic differentiation. Synthesis: We found restricted evidence that forest herb populations with different ages differ in their genetic structure, indicating that populations of different ages can reach a similar genetic structure within decades and thus persist in the long term after habitat disturbance. Despite their distinct breeding systems and associated pollinators, the three studied species exhibited partly similar genetic patterns, suggesting that their common characteristics, such as being slow colonizers or their ability to propagate vegetatively, are important in determining their long-term response to land-cover change.

Detection of Genetic Patterns in Endangered Marine Species Is Affected by Small Sample Sizes

Knowledge of Genetic diversity and its spatial distribution is crucial to improve conservation plans for endangered species. Genetic tools help ensure species' long-term persistence by unraveling connectivity patterns and evolutionary trajectories of populations. Here, microsatellite genotypes of individuals from populations of Patella ferruginea are used to assess the effect of sample size on metrics of within-and between-population genetic diversity by combining empirical and simulated data. Within-population metrics are slightly to moderately affected by small sample size, albeit the magnitude of the bias is proportional to the effective population size and gene flow. The power of detecting genetic differentiation among populations increases with sample size, albeit the gain of increasing the number of sampled individuals tends to be negligible between 30 and 50. Our results line up with those of previous studies and highlight that small sample sizes are not always a hindrance to investigating genetic patterns in endangered marine species. Caution is needed in interpreting genetic patterns based on small sample sizes when the observed genetic differentiation is weak. This study also highlights the importance of carrying out genetic monitoring in seemingly well-preserved but potentially isolated populations.

Applying novel connectivity networks to wood turtle populations to provide comprehensive conservation management strategies for species at risk

Genetic diversity within and among populations is frequently used in prioritization processes to rank populations based on their vulnerability or distinctiveness, however, connectivity and gene flow are rarely considered within these frameworks. Using a wood turtle (Glyptemys insculpta) population graph, we introduce BRIDES as a new tool to evaluate populations for conservation purpose without focusing solely on individual nodes. BRIDES characterizes different types of shortest paths among the nodes of a subgraph and compares the shortest paths among the same nodes in a complete network. The main objectives of this study were to (1) introduce a BRIDES selection process to assist conservation biologists in the prioritization of populations, and (2) use different centrality indices and node removal statistics to compare BRIDES results and assess gene flow among wood turtle populations. We constructed six population subgraphs and used a stepwise selection algorithm to choose the optimal number of additional nodes, representing different populations, required to maximize network connectivity under different weighting schemes. Our results demonstrate the robustness of the BRIDES selection process for a given scenario, while inconsistencies were observed among node-based metrics. Results showed repeated selection of certain wood turtle populations, which could have not been predicted following only genetic diversity and distinctiveness estimation, node-based metrics and node removal analysis. Contrary to centrality measures focusing on static networks, BRIDES allowed for the analysis of evolving networks. To our knowledge, this study is the first to apply graph theory for turtle conservation genetics. We show that population graphs can reveal complex gene flow dynamics and population resiliency to local extinction. As such, BRIDES offers an interesting complement to node-based metrics and node removal to better understand the global processes at play when addressing population prioritization frameworks.

Analysing landscape effects on dispersal networks and gene flow with genetic graphs

Graph-theoretic approaches have relevant applications in landscape genetic analyses. When species form populations in discrete habitat patches, genetic graphs can be used (a) to identify direct dispersal paths followed by propagules or (b) to quantify landscape effects on multi-generational gene flow. However, the influence of their construction parameters remains to be explored. Using a simulation approach, we constructed genetic graphs using several pruning methods (geographical distance thresholds, topological constraints, statistical inference) and genetic distances to weight graph links (F_ST , D_PS , Euclidean genetic distances). We then compared the capacity of these different graphs to (a) identify the precise topology of the dispersal network and (b) to infer landscape resistance to gene flow from the relationship between cost-distances and genetic distances. Although not always clear-cut, our results showed that methods based on geographical distance thresholds seem to better identify dispersal networks in most cases. More interestingly, our study demonstrates that a sub-selection of pairwise distances through graph pruning (thereby reducing the number of data points) can counter-intuitively lead to improved inferences of landscape effects on dispersal. Finally, we showed that genetic distances such as the D_PS or Euclidean genetic distances should be preferred over the F_ST for landscape effect inference as they respond faster to landscape changes.

Re-examination of population structure in Arctic ringed seals using DArTseq genotyping

Although Arctic ringed seals Phoca hispida hispida are currently abundant and broadly distributed, their numbers are projected to decline substantially by the year 2100 due to climate warming. While understanding population structure could provide insight into the impact of environmental changes on this subspecies, detecting demographically important levels of exchange can be difficult in taxa with high abundance. We used a next-generation sequencing approach (DArTseq) to genotype ~5700 single nucleotide polymorphisms in 79 seals from 4 Pacific Arctic regions. Comparison of the 2 most geographically separated strata (eastern Bering vs. northeastern Chukchi-Beaufort Seas) revealed a statistically significant level of genetic differentiation (FST = 0.001, p = 0.005) that, while small, was 1 to 2 orders of magnitude greater than expected based on divergence estimated for similarly sized populations connected by low (1% yr-1) dispersal. A relatively high proportion (72 to 88%) of individuals within these strata could be genetically assigned to their stratum of origin. These results indicate that demographically important structure may be present among Arctic ringed seals breeding in different areas, increasing the risk that declines in the number of seals breeding in areas most negatively affected by environmental warming could occur.

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.

Urban hubs of connectivity: contrasting patterns of gene flow within and among cities in the western black widow spider

As urbanization drastically alters the natural landscape and generates novel habitats within cities, the potential for changes to gene flow for urban-dwelling species increases. The western black widow spider (Latrodectus hesperus) is a medically relevant urban adapter pest species, for which we have previously identified population genetic signatures consistent with urbanization facilitating gene flow, likely due to human-mediated transport. Here, in an analysis of 1.9 million genome-wide SNPs, we contrast broad-scale geographical analyses of 10 urban and 11 non-urban locales with fine-scale within-city analyses including 30 urban locales across the western USA. These hierarchical datasets enable us to test hypotheses of how urbanization impacts multiple urban cities and their genetic connectivity at different spatial scales. Coupled fine-scale and broad-scale analyses reveal contrasting patterns of high and low genetic differentiation among locales within cities as a result of low and high genetic connectivity, respectively, of these cities to the overall population network. We discuss these results as they challenge the use of cities as replicates of urban eco-evolution, and have implications for conservation and human health in a rapidly growing urban habitat.

The genetic network of greater sage-grouse: Range-wide identification of keystone hubs of connectivity

Genetic networks can characterize complex genetic relationships among groups of individuals, which can be used to rank nodes most important to the overall connectivity of the system. Ranking allows scarce resources to be guided toward nodes integral to connectivity. The greater sage-grouse (Centrocercus urophasianus) is a species of conservation concern that breeds on spatially discrete leks that must remain connected by genetic exchange for population persistence. We genotyped 5,950 individuals from 1,200 greater sage-grouse leks distributed across the entire species' geographic range. We found a small-world network composed of 458 nodes connected by 14,481 edges. This network was composed of hubs-that is, nodes facilitating gene flow across the network-and spokes-that is, nodes where connectivity is served by hubs. It is within these hubs that the greatest genetic diversity was housed. Using indices of network centrality, we identified hub nodes of greatest conservation importance. We also identified keystone nodes with elevated centrality despite low local population size. Hub and keystone nodes were found across the entire species' contiguous range, although nodes with elevated importance to network-wide connectivity were found more central: especially in northeastern, central, and southwestern Wyoming and eastern Idaho. Nodes among which genes are most readily exchanged were mostly located in Montana and northern Wyoming, as well as Utah and eastern Nevada. The loss of hub or keystone nodes could lead to the disintegration of the network into smaller, isolated subnetworks. Protecting both hub nodes and keystone nodes will conserve genetic diversity and should maintain network connections to ensure a resilient and viable population over time. Our analysis shows that network models can be used to model gene flow, offering insights into its pattern and process, with application to prioritizing landscapes for conservation.

Consequences of population topology for studying gene flow using link-based landscape genetic methods

Many landscape genetic studies aim to determine the effect of landscape on gene flow between populations. These studies frequently employ link‐based methods that relate pairwise measures of historical gene flow to measures of the landscape and the geographical distance between populations. However, apart from landscape and distance, there is a third important factor that can influence historical gene flow, that is, population topology (i.e., the arrangement of populations throughout a landscape). As the population topology is determined in part by the landscape configuration, I argue that it should play a more prominent role in landscape genetics. Making use of existing literature and theoretical examples, I discuss how population topology can influence results in landscape genetic studies and how it can be taken into account to improve the accuracy of these results. In support of my arguments, I have performed a literature review of landscape genetic studies published during the first half of 2015 as well as several computer simulations of gene flow between populations. First, I argue why one should carefully consider which population pairs should be included in link‐based analyses. Second, I discuss several ways in which the population topology can be incorporated in response and explanatory variables. Third, I outline why it is important to sample populations in such a way that a good representation of the population topology is obtained. Fourth, I discuss how statistical testing for link‐based approaches could be influenced by the population topology. I conclude the article with six recommendations geared toward better incorporating population topology in link‐based landscape genetic studies.

Sex-specific graphs: Relating group-specific topology to demographic and landscape data

Sex-specific genetic structure is a commonly observed pattern among vertebrate species. Facing differential selective pressures, individuals may adopt sex-specific life history traits that ultimately shape genetic variation among populations. Although differential dispersal dynamics are commonly detected in the literature, few studies have used genetic structure to investigate sex-specific functional connectivity. The recent use of graph theoretic approaches in landscape genetics has demonstrated network capacities to describe complex system behaviours where network topology represents genetic interaction among subunits. Here, we partition the overall genetic structure into sex-specific graphs, revealing different male and female dispersal dynamics of a fisher (Pekania [Martes] pennanti) metapopulation in southern Ontario. Our analyses based on network topologies supported the hypothesis of male-biased dispersal. Furthermore, we demonstrated that the effect of the landscape, identified at the population level, could be partitioned among sex-specific strata. We found that female connectivity was negatively correlated with snow depth, whereas connectivity among males was not. Our findings underscore the potential of conducting sex-specific analysis by identifying landscape elements or configuration that differentially promotes or impedes functional connectivity between sexes, revealing processes that may otherwise remain cryptic. We propose that the sex-specific graph approach would be applicable to other vagile species where differential sex-specific processes are expected to occur.

Multi-species genetic connectivity in a terrestrial habitat network

BACKGROUND

Habitat fragmentation reduces genetic connectivity for multiple species, yet conservation efforts tend to rely heavily on single-species connectivity estimates to inform land-use planning. Such conservation activities may benefit from multi-species connectivity estimates, which provide a simple and practical means to mitigate the effects of habitat fragmentation for a larger number of species. To test the validity of a multi-species connectivity model, we used neutral microsatellite genetic datasets of Canada lynx (Lynx canadensis), American marten (Martes americana), fisher (Pekania pennanti), and southern flying squirrel (Glaucomys volans) to evaluate multi-species genetic connectivity across Ontario, Canada.

RESULTS

We used linear models to compare node-based estimates of genetic connectivity for each species to point-based estimates of landscape connectivity (current density) derived from circuit theory. To our knowledge, we are the first to evaluate current density as a measure of genetic connectivity. Our results depended on landscape context: habitat amount was more important than current density in explaining multi-species genetic connectivity in the northern part of our study area, where habitat was abundant and fragmentation was low. In the south however, where fragmentation was prevalent, genetic connectivity was correlated with current density. Contrary to our expectations however, locations with a high probability of movement as reflected by high current density were negatively associated with gene flow. Subsequent analyses of circuit theory outputs showed that high current density was also associated with high effective resistance, underscoring that the presence of pinch points is not necessarily indicative of gene flow.

CONCLUSIONS

Overall, our study appears to provide support for the hypothesis that landscape pattern is important when habitat amount is low. We also conclude that while current density is proportional to the probability of movement per unit area, this does not imply increased gene flow, since high current density tends to be a result of neighbouring pixels with high cost of movement (e.g., low habitat amount). In other words, pinch points with high current density appear to constrict gene flow.

Navigating the pitfalls and promise of landscape genetics

The field of landscape genetics has been evolving rapidly since its emergence in the early 2000s. New applications, techniques and criticisms of techniques appear like clockwork with each new journal issue. The developments are an encouraging, and at times bewildering, sign of progress in an exciting new field of study. However, we suggest that the rapid expansion of landscape genetics has belied important flaws in the development of the field, and we add an air of caution to this breakneck pace of expansion. Specifically, landscape genetic studies often lose sight of the fundamental principles and complex consequences of gene flow, instead favouring simplistic interpretations and broad inferences not necessarily warranted by the data. Here, we describe common pitfalls that characterize such studies, and provide practical guidance to improve landscape genetic investigation, with careful consideration of inferential limits, scale, replication, and the ecological and evolutionary context of spatial genetic patterns. Ultimately, the utility of landscape genetics will depend on translating the relationship between gene flow and landscape features into an understanding of long-term population outcomes. We hope the perspective presented here will steer landscape genetics down a more scientifically sound and productive path, garnering a field that is as informative in the future as it is popular now.

Node-based measures of connectivity in genetic networks

At-site environmental conditions can have strong influences on genetic connectivity, and in particular on the immigration and settlement phases of dispersal. However, at-site processes are rarely explored in landscape genetic analyses. Networks can facilitate the study of at-site processes, where network nodes are used to model site-level effects. We used simulated genetic networks to compare and contrast the performance of 7 node-based (as opposed to edge-based) genetic connectivity metrics. We simulated increasing node connectivity by varying migration in two ways: we increased the number of migrants moving between a focal node and a set number of recipient nodes, and we increased the number of recipient nodes receiving a set number of migrants. We found that two metrics in particular, the average edge weight and the average inverse edge weight, varied linearly with simulated connectivity. Conversely, node degree was not a good measure of connectivity. We demonstrated the use of average inverse edge weight to describe the influence of at-site habitat characteristics on genetic connectivity of 653 American martens (Martes americana) in Ontario, Canada. We found that highly connected nodes had high habitat quality for marten (deep snow and high proportions of coniferous and mature forest) and were farther from the range edge. We recommend the use of node-based genetic connectivity metrics, in particular, average edge weight or average inverse edge weight, to model the influences of at-site habitat conditions on the immigration and settlement phases of dispersal.

Genetic attributes of midwife toad (Alytes obstetricans) populations do not correlate with degree of species decline

Genetic diversity is crucial for long-term population persistence. Population loss and subsequent reduction in migration rate among the most important processes that are expected to lead to a reduction in genetic diversity and an increase in genetic differentiation. While the theory behind this is well-developed, empirical evidence from wild populations is inconsistent. Using microsatellite markers, we compared the genetic structure of populations of an amphibian species, the midwife toad (Alytes obstetricans), in four Swiss regions where the species has suffered variable levels of subpopulation extirpation. We also quantified the effects of several geographic factors on genetic structure and used a model selection approach to ascertain which of the variables were important for explaining genetic variation. Although subpopulation pairwise F ST-values were highly significant even over small geographic scales, neither any of the geographic variables nor loss of subpopulations were important factors for predicting spatial genetic structure. The absence of a signature of subpopulation loss on genetic differentiation may suggest that midwife toad subpopulations function as relatively independent units.