Isolation by environment in White‐breasted Nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach

Understanding landscape processes driving patterns of population genetic differentiation and diversity has been a long‐standing focus of ecology and evolutionary biology. Gene flow may be reduced by historical, ecological or geographic factors, resulting in patterns of isolation by distance (IBD) or isolation by environment (IBE). Although IBE has been found in many natural systems, most studies investigating patterns of IBD and IBE in nature have used anonymous neutral genetic markers, precluding inference of selection mechanisms or identification of genes potentially under selection. Using landscape genomics, the simultaneous study of genomic and ecological landscapes, we investigated the processes driving population genetic patterns of White‐breasted Nuthatches (Sitta carolinensis) in sky islands (montane forest habitat islands) of the Madrean Archipelago. Using more than 4000 single nucleotide polymorphisms and multiple tests to investigate the relationship between genetic differentiation and geographic or ecological distance, we identified IBE, and a lack of IBD, among sky island populations of S. carolinensis. Using three tests to identify selection, we found 79 loci putatively under selection; of these, seven matched CDS regions in the Zebra Finch. The loci under selection were highly associated with climate extremes (maximum temperature of warmest month and minimum precipitation of driest month). These results provide evidence for IBE – disentangled from IBD – in sky island vertebrates and identify potential adaptive genetic variation.     

Ecological connectivity assessment in a strongly structured fire salamander (Salamandra salamandra) population

Small populations are more prone to extinction if the dispersal among them is not adequately maintained by ecological connections. The degree of isolation between populations could be evaluated measuring their genetic distance, which depends on the respective geographic (isolation by distance, IBD) and/or ecological (isolation by resistance, IBR) distances. The aim of this study was to assess the ecological connectivity of fire salamander Salamandra salamandra populations by means of a landscape genetic approach. The species lives in broad‐leaved forest ecosystems and is particularly affected by fragmentation due to its habitat selectivity and low dispersal capability. We analyzed 477 biological samples collected in 47 sampling locations (SLs) in the mainly continuous populations of the Prealpine and Eastern foothill lowland (PEF) and 10 SLs in the fragmented populations of the Western foothill (WF) lowland of Lombardy (northern Italy). Pairwise genetic distances (Chord distance, D_C) were estimated from allele frequencies of 16 microsatellites loci. Ecological distances were calculated using one of the most promising methodology in landscape genetics studies, the circuit theory, applied to habitat suitability maps. We realized two habitat suitability models: one without barriers (EcoD) and a second one accounting for the possible barrier effect of main roads (EcoD_b). Mantel tests between distance matrices highlighted how the Log‐D_C in PEF populations was related to log‐transformed geographic distance (confirming a prevalence of IBD), while it was explained by the Log‐EcoD, and particularly by the Log‐EcoD_b, in WF populations, even when accounting for the confounding effect of geographic distance (highlighting a prevalence of IBR). Moreover, we also demonstrated how considering the overall population, the effect of Euclidean or ecological distances on genetic distances acting at the level of a single group (PEF or WF populations) could not be detected, when population are strongly structured.     

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high‐latitude B. huntii habitats exhibit low ENS when compared to low‐latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low‐latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Drivers of population genetic differentiation in the wild: isolation by dispersal limitation,isolation by adaptation and isolation by colonization

Empirical population genetic studies have been dominated by a neutralist view, according to which gene flow and drift are the main forces driving population genetic structure in nature. The neutralist view in essence describes a process of isolation by dispersal limitation (IBDL) that generally leads to a pattern of isolation by distance (IBD). Recently, however, conceptual frameworks have been put forward that view local genetic adaptation as an important driver of population genetic structure. Isolation by adaptation (IBA) and monopolization (M) posit that gene flow among natural populations is reduced as a consequence of local genetic adaptation. IBA stresses that effective gene flow is reduced among habitats that show dissimilar ecological characteristics, leading to a pattern of isolation by environment. In monopolization, local genetic adaptation of initial colonizing genotypes results in a reduction in gene flow that fosters the persistence of founder effects. Here, we relate these different processes driving landscape genetic structure to patterns of IBD and isolation by environment (IBE). We propose a method to detect whether IBDL, IBA and M shape genetic differentiation in natural landscapes by studying patterns of variation at neutral and non‐neutral markers as well as at ecologically relevant traits. Finally, we reinterpret a representative number of studies from the recent literature by associating patterns to processes and identify patterns associated with local genetic adaptation to be as common as IBDL in structuring regional genetic variation of populations in the wild. Our results point to the importance of quantifying environmental gradients and incorporating ecology in the analysis of population genetics.     

Isolation by resistance

Despite growing interest in the effects of landscape heterogeneity on genetic structuring, few tools are available to incorporate data on landscape composition into population genetic studies. Analyses of isolation by distance have typically either assumed spatial homogeneity for convenience or applied theoretically unjustified distance metrics to compensate for heterogeneity. Here I propose the isolation-by-resistance (IBR) model as an alternative for predicting equilibrium genetic structuring in complex landscapes. The model predicts a positive relationship between genetic differentiation and the resistance distance, a distance metric that exploits precise relationships between random walk times and effective resistances in electronic networks. As a predictor of genetic differentiation, the resistance distance is both more theoretically justified and more robust to spatial heterogeneity than Euclidean or least cost path-based distance measures. Moreover, the metric can be applied with a wide range of data inputs, including coarse-scale range maps, simple maps of habitat and nonhabitat within a species' range, or complex spatial datasets with habitats and barriers of differing qualities. The IBR model thus provides a flexible and efficient tool to account for habitat heterogeneity in studies of isolation by distance, improve understanding of how landscape characteristics affect genetic structuring, and predict genetic and evolutionary consequences of landscape change.     

Isolation by distance,resistance and/or clusters? Lessons learned from a forest‐dwelling carnivore inhabiting a heterogeneous landscape

Landscape genetics provides a valuable framework to understand how landscape features influence gene flow and to disentangle the factors that lead to discrete and/or clinal population structure. Here, we attempt to differentiate between these processes in a forest‐dwelling small carnivore [European pine marten (Martes martes)]. Specifically, we used complementary analytical approaches to quantify the spatially explicit genetic structure and diversity and analyse patterns of gene flow for 140 individuals genotyped at 15 microsatellite loci. We first used spatially explicit and nonspatial Bayesian clustering algorithms to partition the sample into discrete clusters and evaluate hypotheses of ‘isolation by barriers’ (IBB). We further characterized the relationships between genetic distance and geographical (‘isolation by distance’, IBD) and ecological distances (‘isolation by resistance’, IBR) obtained from optimized landscape models. Using a reciprocal causal modelling approach, we competed the IBD, IBR and IBB hypotheses with each other to unravel factors driving population genetic structure. Additionally, we further assessed spatially explicit indices of genetic diversity using sGD across potentially overlapping genetic neighbourhoods that matched the inferred population structure. Our results revealed a complex spatial genetic cline that appears to be driven jointly by IBD and partial barriers to gene flow (IBB) associated with poor habitat and interspecific competition. Habitat loss and fragmentation, in synergy with past overharvesting and possible interspecific competition with sympatric stone marten (Martes foina), are likely the main factors responsible for the spatial genetic structure we observed. These results emphasize the need for a more thorough evaluation of discrete and clinal hypotheses governing gene flow in landscape genetic studies, and the potential influence of different limiting factors affecting genetic structure at different spatial scales.     

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes can drive repeated patterns of evolution in different groups of populations experiencing similar ecological gradients. In this paper, we used a combination of nuclear and mitochondrial DNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological and genetic divergence of European minnows in two mountain ranges: the Pyrenees and the Alps. European minnows (Phoxinus phoxinus) are cyprinid fish inhabiting most freshwater bodies in Europe, including those in different mountain ranges that could act as major geographical barriers to gene flow. We explored patterns of P. phoxinus phenotypic and genetic diversification along a gradient of altitude common to the two mountain ranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation (IBA). The results indicated that populations from the Pyrenees and the Alps belong to two well differentiated, reciprocally monophyletic mtDNA lineages. Substantial genetic differentiation due to geographical isolation within and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPs markers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Also, morphology varied strongly with elevation and so did genetic differentiation to a lower extent. Despite moderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneity was found at the genetic level, suggesting that selection and population specific genetic drift act in concert to affect genetic divergence.     

Effect of landscape features on genetic structure of the goitered gazelle (<Emphasis Type="Italic">Gazella subgutturosa</Emphasis>) in Central Iran

The populations of goitered gazelle suffered significant decline due to natural and anthropogenic factors over the last century. Investigating the effects of barriers on gene flow among the remaining populations is vital for conservation planning. Here we adopted a landscape genetics approach to evaluate the genetic structure of the goitered gazelle in Central Iran and the effects of landscape features on gene flow using 15 polymorphic microsatellite loci. Spatial autocorrelation, isolation by distance (IBD) and isolation by resistance (IBR) models were used to elucidate the effects of landscape features on the genetic structure. Ecological modeling was used to construct landscape permeability and resistance map using 12 ecogeographical variables. Bayesian algorithms revealed three genetically homogeneous groups and restricted dispersal pattern in the six populations. The IBD and spatial autocorrelation revealed a pattern of decreasing relatedness with increasing distance. The distribution of potential habitats was strongly correlated with bioclimatic factors, vegetation type, and elevation. Resistance distances and graph theory were significantly related with variation in genetic structure, suggesting that gazelles are affected by landscape composition. The IBD showed greater impact on genetic structure than IBR. The Mantel and partial Mantel tests indicated low but non-significant effects of anthropogenic barriers on observed genetic structure. We concluded that a combination of geographic distance, landscape resistance, and anthropogenic factors are affecting the genetic structure and gene flow of populations. Future road construction might impede connectivity and gene exchange of populations. Conservation measures on this vulnerable species should consider some isolated population as separate management units.     

Adaptation with gene flow across the landscape in a dune sunflower

Isolation by adaptation increases divergence at neutral loci when natural selection against immigrants reduces the rate of gene flow between different habitats. This can occur early in the process of adaptive divergence and is a key feature of ecological speciation. Despite the ability of isolation by distance (IBD) and other forms of landscape resistance to produce similar patterns of neutral divergence within species, few studies have used landscape genetics to control for these other forces. We have studied the divergence of Helianthus petiolaris ecotypes living in active sand dunes and adjacent non-dune habitat, using landscape genetics approaches, such as circuit theory and multiple regression of distance matrices, in addition to coalescent modelling. Divergence between habitats was significant, but not strong, and was shaped by IBD. We expected that increased resistance owing to patchy and unfavourable habitat in the dunes would contribute to divergence. Instead, we found that landscape resistance models with lower resistance in the dunes performed well as predictors of genetic distances among subpopulations. Nevertheless, habitat class remained a strong predictor of genetic distance when controlling for isolation by resistance and IBD. We also measured environmental variables at each site and confirmed that specific variables, especially soil nitrogen and vegetation cover, explained a greater proportion of variance in genetic distance than did landscape or the habitat classification alone. Asymmetry in effective population sizes and numbers of migrants per generation was detected using coalescent modelling with Bayesian inference, which is consistent with incipient ecological speciation being driven by the dune habitat.     

Genomic footprints of adaptation in a cooperatively breeding tropical bird across a vegetation gradient

Identifying the genetic basis of phenotypic variation and its relationship with the environment is key to understanding how local adaptations evolve. Such patterns are especially interesting among populations distributed across habitat gradients, where genetic structure can be driven by isolation by distance (IBD) and/or isolation by environment (IBE). Here, we used variation in ~1,600 high‐quality SNPs derived from paired‐end sequencing of double‐digest restriction site‐associated DNA (ddRAD‐Seq) to test hypotheses related to IBD and IBE in the Yucatan jay (Cyanocorax yucatanicus), a tropical bird endemic to the Yucatán Peninsula. This peninsula is characterized by a precipitation and vegetation gradient—from dry to evergreen tropical forests—that is associated with morphological variation in this species. We found a moderate level of nucleotide diversity (π = .008) and little evidence for genetic differentiation among vegetation types. Analyses of neutral and putatively adaptive SNPs (identified by complementary genome‐scan approaches) indicate that IBD is the most reliable explanation to account for frequency distribution of the former, while IBE has to be invoked to explain those of the later. These results suggest that selective factors acting along a vegetation gradient can promote local adaptation in the presence of gene flow in a vagile, nonmigratory and geographically restricted species. The putative candidate SNPs identified here are located within or linked to a variety of genes that represent ideal targets for future genomic surveys.     

Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea

The relatively recent fields of terrestrial landscape and marine seascape genetics seek to identify the influence of biophysical habitat features on the spatial genetic structure of populations or individuals. Over the last few years, there has been accumulating evidence for the effect of environmental heterogeneity on patterns of gene flow and connectivity in marine systems. Here, we investigate the population genetic patterns of an anemonefish, Amphiprion bicinctus, along the Saudi Arabian coast of the Red Sea. We collected nearly one thousand samples from 19 locations, spanning approximately 1500 km, and genotyped them at 38 microsatellite loci. Patterns of gene flow appeared to follow a stepping‐stone model along the northern and central Red Sea, which was disrupted by a distinct genetic break at a latitude of approximately 19°N. The Red Sea is characterized by pronounced environmental gradients along its axis, roughly separating the northern and central from the southern basin. Using mean chlorophyll‐a concentrations as a proxy for this gradient, we ran tests of isolation by distance (IBD, R² = 0.52) and isolation by environment (IBE, R² = 0.64), as well as combined models using partial Mantel tests and multiple matrix regression with randomization (MMRR). We found that genetic structure across our sampling sites may be best explained by a combined model of IBD and IBE (Mantel: R² = 0.71, MMRR: R² = 0.86). Our results highlight the potential key role of environmental patchiness in shaping patterns of gene flow in species with pelagic larval dispersal. We support growing calls for the integration of biophysical habitat characteristics into future studies of population genetic structure.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Drivers of genetic differentiation in a generalist insect‐pollinated herb across spatial scales

The isolation‐by‐distance model (IBD) predicts that genetic differentiation among populations increases with geographic distance. Yet, empirical studies show that a variety of ecological, topographic and historical factors may override the effect of geographic distance on genetic variation. This may particularly apply to species with narrow but highly heterogeneous distribution ranges, such as those occurring along elevational gradients. Using nine SSR markers, we study the genetic differentiation of the montane pollination‐generalist herb, Erysimum mediohispanicum. Because the effects of any given factor may depend on the geographic scale considered, we investigate the contribution of different environmental and historical factors at three different spatial scales. We evaluate five competing models that put forward the role of geographic distance, local environmental factors [biotic interactions (IBEb) and climatic variables (IBEa)], landscape resistance (IBR) and phylogeographic patterns (IBP), respectively. We find significant IBD regardless of the spatial scale and the genetic distance estimator considered. However, IBEa and IBP also play a prominent role in shaping genetic differentiation patterns at the larger spatial scales, and IBR is significant at the fine spatial scale. Overall, our results highlight the importance of combining different estimators, statistical approaches and spatial scales to disentangle the relative importance of the various ecological factors contributing to the shaping of genetic divergence patterns in natural populations.     

Landscape genetics of a pollinator longhorn beetle [Typocerus v. velutinus (Olivier)] on a continuous habitat surface

Landscape connectivity, the degree to which the landscape structure facilitates or impedes organismal movement and gene flow, is increasingly important to conservationists and land managers. Metrics for describing the undulating shape of continuous habitat surfaces can expand the usefulness of continuous gradient surfaces that describe habitat and predict the flow of organisms and genes. We adopted a landscape gradient model of habitat and used surface metrics of connectivity to model the genetic continuity between populations of the banded longhorn beetle [Typocerus v. velutinus (Olivier)] collected at 17 sites across a fragmentation gradient in Indiana, USA. We tested the hypothesis that greater habitat connectivity facilitates gene flow between beetle populations against a null model of isolation by distance (IBD). We used next‐generation sequencing to develop 10 polymorphic microsatellite loci and genotype the individual beetles to assess the population genetic structure. Isolation by distance did not explain the population genetic structure. The surface metrics model of habitat connectivity explained the variance in genetic dissimilarities 30 times better than the IBD model. We conclude that surface metrology of habitat maps is a powerful extension of landscape genetics in heterogeneous landscapes.     

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.     

GENETIC ISOLATION BY ENVIRONMENT OR DISTANCE: WHICH PATTERN OF GENE FLOW IS MOST COMMON?

Gene flow among populations can enhance local adaptation if it introduces new genetic variants available for selection, but strong gene flow can also stall adaptation by swamping locally beneficial genes. These outcomes can depend on population size, genetic variation, and the environmental context. Gene flow patterns may align with geographic distance (IBD—isolation by distance), whereby immigration rates are inversely proportional to the distance between populations. Alternatively gene flow may follow patterns of isolation by environment (IBE), whereby gene flow rates are higher among similar environments. Finally, gene flow may be highest among dissimilar environments (counter‐gradient gene flow), the classic “gene‐swamping” scenario. Here we survey relevant studies to determine the prevalence of each pattern across environmental gradients. Of 70 studies, we found evidence of IBD in 20.0%, IBE in 37.1%, and both patterns in 37.1%. In addition, 10.0% of studies exhibited counter‐gradient gene flow. In total, 74.3% showed significant IBE patterns. This predominant IBE pattern of gene flow may have arisen directly through natural selection or reflect other adaptive and nonadaptive processes leading to nonrandom gene flow. It also precludes gene swamping as a widespread phenomenon. Implications for evolutionary processes and management under rapidly changing environments (e.g., climate change) are discussed.     

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour‐joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within‐prairie genetic diversity (92%). Using Bayenv 2, 14 top‐ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScan F_ST outliers were in common with Bayenv 2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.     

Effects of sample size,number of markers,and allelic richness on the detection of spatial genetic pattern

The influence of study design on the ability to detect the effects of landscape pattern on gene flow is one of the most pressing methodological gaps in landscape genetic research. To investigate the effect of study design on landscape genetics inference, we used a spatially‐explicit, individual‐based program to simulate gene flow in a spatially continuous population inhabiting a landscape with gradual spatial changes in resistance to movement. We simulated a wide range of combinations of number of loci, number of alleles per locus and number of individuals sampled from the population. We assessed how these three aspects of study design influenced the statistical power to successfully identify the generating process among competing hypotheses of isolation‐by‐distance, isolation‐by‐barrier, and isolation‐by‐landscape resistance using a causal modelling approach with partial Mantel tests. We modelled the statistical power to identify the generating process as a response surface for equilibrium and non‐equilibrium conditions after introduction of isolation‐by‐landscape resistance. All three variables (loci, alleles and sampled individuals) affect the power of causal modelling, but to different degrees. Stronger partial Mantel r correlations between landscape distances and genetic distances were found when more loci were used and when loci were more variable, which makes comparisons of effect size between studies difficult. Number of individuals did not affect the accuracy through mean equilibrium partial Mantel r, but larger samples decreased the uncertainty (increasing the precision) of equilibrium partial Mantel r estimates. We conclude that amplifying more (and more variable) loci is likely to increase the power of landscape genetic inferences more than increasing number of individuals.