Habitat use,but not gene flow,is influenced by human activities in two ecotypes of Egyptian fruit bat (Rousettus aegyptiacus)

Understanding the ecological, behavioural and evolutionary response of organisms to changing environments is of primary importance in a human‐altered world. It is crucial to elucidate how human activities alter gene flow and what are the consequences for the genetic structure of a species. We studied two lineages of the Egyptian fruit bat (Rousettus aegyptiacus) throughout the contact zone between mesic and arid Ecozones in the Middle East to evaluate the species' response to the growing proportion of human‐altered habitats in the desert. We integrated population genetics, morphometrics and movement ecology to analyse population structure, morphological variation and habitat use from GPS‐ or radio‐tagged individuals from both desert and Mediterranean areas. We classified the spatial distribution and environmental stratification by describing physical–geographical conditions and land cover. We analysed this information to estimate patch occupancy and used an isolation‐by‐resistance approach to model gene flow patterns. Our results suggest that lineages from desert and Mediterranean habitats, despite their admixture, are isolated by environment and by adaptation supporting their classification as ecotypes. We found a positive effect of human‐altered habitats on patch occupancy and habitat use of fruit bats by increasing the availability of roosting and foraging areas. While this commensalism promotes the distribution of fruit bats throughout the Middle East, gene flow between colonies has not been altered by human activities. This discrepancy between habitat use and gene flow patterns may, therefore, be explained by the breeding system of the species and modifications of natal dispersal patterns.     

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.     

The influence of landscape configuration and environment on population genetic structure in a sedentary passerine: insights from loci located in different genomic regions

The study of the factors structuring genetic variation can help to infer the neutral and adaptive processes shaping the demographic and evolutionary trajectories of natural populations. Here, we analyse the role of isolation by distance (IBD), isolation by resistance (IBR, defined by landscape composition) and isolation by environment (IBE, estimated as habitat and elevation dissimilarity) in structuring genetic variation in 25 blue tit (Cyanistes caeruleus) populations. We typed 1385 individuals at 26 microsatellite loci classified into two groups by considering whether they are located into genomic regions that are actively (TL; 12 loci) or not (NTL; 14 loci) transcribed to RNA. Population genetic differentiation was mostly detected using the panel of NTL. Landscape genetic analyses showed a pattern of IBD for all loci and the panel of NTL, but genetic differentiation estimated at TL was only explained by IBR models considering high resistance for natural vegetation and low resistance for agricultural lands. Finally, the absence for IBE suggests a lack of divergent selection pressures associated with differences in habitat and elevation. Overall, our study shows that markers located in different genomic regions can yield contrasting inferences on landscape‐level patterns of realized gene flow in natural populations.     

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.     

Environmental selection is a main driver of divergence in house sparrows (Passer domesticus) in Romania and Bulgaria

Both neutral and adaptive evolutionary processes can cause population divergence, but their relative contributions remain unclear. We investigated the roles of these processes in population divergence in house sparrows (Passer domesticus) from Romania and Bulgaria, regions characterized by high landscape heterogeneity compared to Western Europe. We asked whether morphological divergence, complemented with genetic data in this human commensal species, was best explained by environmental variation, geographic distance, or landscape resistance—the effort it takes for an individual to disperse from one location to the other—caused by either natural or anthropogenic barriers. Using generalized dissimilarity modeling, a matrix regression technique that fits biotic beta diversity to both environmental predictors and geographic distance, we found that a small set of climate and vegetation variables explained up to ~30% of the observed divergence, whereas geographic and resistance distances played much lesser roles. Our results are consistent with signals of selection on morphological traits and of isolation by adaptation in genetic markers, suggesting that selection by natural environmental conditions shapes population divergence in house sparrows. Our study thus contributes to a growing body of evidence that adaptive evolution may be a major driver of diversification.     

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.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

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.     

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.     

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east‐west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.     

Living in two worlds: Evolutionary mechanisms act differently in the native and introduced ranges of an invasive plant

Identifying the factors that influence spatial genetic structure among populations can provide insights into the evolution of invasive plants. In this study, we used the common reed (Phragmites australis), a grass native in Europe and invading North America, to examine the relative importance of geographic, environmental (represented by climate here), and human effects on population genetic structure and its changes during invasion. We collected samples of P. australis from both the invaded North American and native European ranges and used molecular markers to investigate the population genetic structure within and between ranges. We used path analysis to identify the contributions of each of the three factors—geographic, environmental, and human‐related—to the formation of spatial genetic patterns. Genetic differentiation was observed between the introduced and native populations, and their genetic structure in the native and introduced ranges was different. There were strong effects of geography and environment on the genetic structure of populations in the native range, but the human‐related factors manifested through colonization of anthropogenic habitats in the introduced range counteracted the effects of environment. The between‐range genetic differences among populations were mainly explained by the heterogeneous environment between the ranges, with the coefficient 2.6 times higher for the environment than that explained by the geographic distance. Human activities were the primary contributor to the genetic structure of the introduced populations. The significant environmental divergence between ranges and the strong contribution of human activities to the genetic structure in the introduced range suggest that invasive populations of P. australis have evolved to adapt to a different climate and to human‐made habitats in North America.     

Is population genetic structure of vascular plants shaped more by ecological or geographic factors? A study case on the Mediterranean endemic Centaurea filiformis (Asteraceae)

• All known populations of the Sardinian endemic Centaurea filiformis Viv. (Asteraceae) were studied in order to understand the impact of both geographic and ecological factors on the genetic structuring of this species.
• Fourteen populations and 234 individuals were sampled. The demographic structure of the populations and the reproductive ecology were estimated in 28 plots. Population genetic analyses were based on SSR markers. Genetic structure was investigated by spatial Bayesian methods.
• Average densities of 0.51 individuals m^?2 were detected, with a prevalence of adults. Ten species of pollinators were identified; C. filiformis ability to self‐pollinate and myrmecochory were demonstrated experimentally. The populations displayed an average heterozygosity value of H _e = 0.576 and high genetic differentiation (overall F _ST = 0.218). Bayesian analysis suggests that five is the most probable number of gene pools of origin. A strong correlation between geographic distances and genetic distances among populations was highlighted.
• The demographic population structure of C. filiformis is dominated by adults, suggesting that it is a stable‐regressive or senile species, investing more in local persistence than colonisation ability. Despite the scattered distribution, the populations studied do not present evidence of genetic erosion. The analysis of genetic differentiation reveals very high differentiation levels among populations, thus indicating that effective barriers exist against gene flow. A general conclusion is that population distribution results in a clear genetic structure for the populations studied, and that geography and not ecology is shaping the present distribution of this species.

     

Genetic evidence for species cohesion,substructure and hybrids in spruce

The origin and history of species are shaped by various evolutionary dynamics, including their persistence in the face of potential gene flow from related taxa. In this study, we use broad geographical and taxonomic sampling (2,219 individuals) to establish the distribution of species, hybrids and cryptic genetic variation within the conifer genus Picea (spruce) across western North America. We demonstrate that the six species of spruce in this region are distinguishable based on their genetic composition, and that the more closely related Engelmann spruce (P. engelmannii) and white spruce (P. glauca) have generated numerous and widespread hybrids. These hybrids occur in the central Rocky Mountains, well to the south of the well‐established region of admixture in Canada. Additionally, we provide evidence for subdivision within Engelmann spruce, manifested as a southern Rocky Mountains form, and a northern Rocky Mountain and Cascade mountains (western) form. In the intervening central Rocky Mountains region (forests in Wyoming and adjacent states) we found primarily individuals with admixed ancestry. Following their origin, these species of spruce have interacted repeatedly and in different geographical contexts. Multiple pairs of species have been shown to hybridize, yet the species persist and retain distinguishable compositions. At the same time, large geographical areas exist where hybrids are pervasive. Consequently, spruce provide a case study for the maintenance of species boundaries, particularly for how widespread hybridization need not lead to the collapse and loss of species.     

Living in the city: urban environments shape the evolution of a native annual plant

Urban environments are warmer, have higher levels of atmospheric CO₂ and have altered patterns of disturbance and precipitation than nearby rural areas. These differences can be important for plant growth and are likely to create distinct selective environments. We planted a common garden experiment with seeds collected from natural populations of the native annual plant Lepidium virginicum, growing in five urban and nearby rural areas in the northern United States to determine whether and how urban populations differ from those from surrounding rural areas. When grown in a common environment, plants grown from seeds collected from urban areas bolted sooner, grew larger, had fewer leaves, had an extended time between bolting and flowering, and produced more seeds than plants grown from seeds collected from rural areas. Interestingly, the rural populations exhibited larger phenotypic differences from one another than urban populations. Surprisingly, genomic data revealed that the majority of individuals in each of the urban populations were more closely related to individuals from other urban populations than they were to geographically proximate rural areas – the one exception being urban and rural populations from New York which were nearly identical. Taken together, our results suggest that selection in urban environments favors different traits than selection in rural environments and that these differences can drive adaptation and shape population structure.     

Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.     

Different processes shape the patterns of divergence in the nuclear and chloroplast genomes of a relict tree species in East Asia

Isolation by spatial distance (IBD), environment (IBE), and historical climatic instability (IBI) are three common processes assessed in phylogeographic and/or landscape genetic studies. However, the relative contributions of these three processes with respect to spatial genetic patterns have seldom been compared. Moreover, whether the relative contribution differs in different regions or when assessed using different genetic markers has rarely been reported. Lindera obtusiloba has been found to have two sister genetic clades of chloroplast (cpDNA) and nuclear microsatellite (nSSR), both of which show discontinuous distribution in northern and southern East Asia. In this study, we used the Mantel test and multiple matrix regression with randomization (MMRR) to determine the relative contributions of IBD, IBE, and IBI with respect to L. obtusiloba populations. Independent Mantel tests and MMRR calculations were conducted for two genetic data sets (cpDNA and nSSR) and for different regions (the overall species range, and northern and southern subregions of the range). We found a significant IBI pattern in nSSR divergence for all assessed regions, whereas no clear IBI pattern was detected with respect to cpDNA. In contrast, significant (or marginal) divergent IBD patterns were detected for cpDNA in all regions, whereas although a significant IBE was apparent with respect to the overall range, the effect was not detected in the two subregions. The differences identified in nSSR and cpDNA population divergence may be related to differences in the heredity and ploidy of the markers. Compared with the southern region, the northern region showed less significant correlation patterns, which may be related to the shorter population history and restricted population range. The findings of this study serve to illustrate that comparing between markers or regions can contribute to gaining a better understanding the population histories of different genomes or within different regions of a species' range.     

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.