Population genetic analysis identifies source-sink dynamics for two sympatric garter snake species (Thamnophis elegans and Thamnophis sirtalis)

Population genetic structure can be shaped by multiple ecological and evolutionary factors, but the genetic consequences of these factors for multiple species inhabiting the same environment remain unexplored. We used microsatellite markers to examine the population structures of two coexisting species of garter snake, Thamnophis elegans and Thamnophis sirtalis, to determine if shared landscape and biology imposed similar population genetic structures. These snakes inhabit a series of ponds, lakes and flooded meadows in northern California and tend to converge on prey type wherever they coexist. Both garter snakes had comparable effective population sizes and bidirectional migration rates (estimated using a maximum-likelihood method based on the coalescent) with low but significant levels of genetic differentiation (F(ST) = 0.024 for T. elegans and 0.035 for T. sirtalis). Asymmetrical gene flow revealed large source populations for both species as well as potential sinks, suggesting frequent extinction-recolonization and metapopulation dynamics. In addition, we found a significant correlation between their genetic structures based on both pairwise F(ST)s for shared populations (P = 0.009) and for bidirectional migration rates (P = 0.024). Possible ecological and evolutionary factors influencing similarities and differences in genetic structure for the two species are discussed. Genetic measures of effective population size and migration rates obtained in this study are also compared with estimates obtained from mark-recapture data.     

Historical influences dominate the population genetic structure of a sedentary marine fish,Atlantic wolffish (Anarhichas lupus), across the North Atlantic Ocean

Genetic variation was assessed in Atlantic wolffish, Anarhichas lupus, across the North Atlantic Ocean using microsatellite and amplified fragment length polymorphism (AFLP) markers. Despite unusual life history attributes such as large benthic eggs, large larvae, a limited pelagic stage and relatively sedentary adults, which suggest potential for strong population structure, range‐wide F_ST values were comparable to other marine fishes (≤0.035). Nevertheless, both significant genetic differentiation among regions and isolation by distance were observed, suggesting limited dispersal in this species. AFLP loci, evaluated on a subset of samples, revealed slightly higher F_ST values, but similar patterns of differentiation and isolation‐by‐distance estimates, compared to microsatellites. The genetic structure of Atlantic wolffish has likely been shaped by its post‐glacial history of recolonization, North Atlantic current patterns and continuity of habitat on continental shelves.     

Low genetic diversity and limited genetic structure across the range of the critically endangered Mexican howler monkey (Alouatta palliata mexicana)

Genetic diversity provides populations with the possibility to persist in ever-changing environments, where selective regimes change over time. Therefore, the long-term survival of a population may be affected by its level of genetic diversity. The Mexican howler monkey (Alouatta palliata mexicana) is a critically endangered primate restricted to southeast Mexico. Here, we evaluate the genetic diversity and population structure of this subspecies based on 83 individuals from 31 groups sampled across the distribution range of the subspecies, using 29 microsatellite loci. Our results revealed extremely low genetic diversity (H_O = 0.21, H_E = 0.29) compared to studies of other A. palliata populations and to other Alouatta species. Principal component analysis, a Bayesian clustering method, and analyses of molecular variance did not detect strong signatures of genetic differentiation among geographic populations of this subspecies. Although we detect small but significant F_ST values between populations, they can be explained by a pattern of isolation by distance. These results and the presence of unique alleles in different populations highlight the importance of implementing conservation efforts in multiple populations across the distribution range of A. p. mexicana to preserve its already low genetic diversity. This is especially important given current levels of population isolation due to the extreme habitat fragmentation across the distribution range of this primate.     

Landscape genetic structure of Scirpus mariqueter reveals a putatively adaptive differentiation under strong gene flow in estuaries

Estuarine organisms grow in highly heterogeneous habitats, and their genetic differentiation is driven by selective and neutral processes as well as population colonization history. However, the relative importance of the processes that underlie genetic structure is still puzzling. Scirpus mariqueter is a perennial grass almost limited in the Changjiang River estuary and its adjacent Qiantang River estuary. Here, using amplified fragment length polymorphism (AFLP), a moderate‐high level of genetic differentiation among populations (range F_ST: 0.0310–0.3325) was showed despite large ongoing dispersal. FLOCK assigned all individuals to 13 clusters and revealed a complex genetic structure. Some genetic clusters were limited in peripheries compared with very mixing constitution in center populations, suggesting local adaptation was more likely to occur in peripheral populations. 21 candidate outliers under positive selection were detected, and further, the differentiation patterns correlated with geographic distance, salinity difference, and colonization history were analyzed with or without the outliers. Combined results of AMOVA and IBD based on different dataset, it was found that the effects of geographic distance and population colonization history on isolation seemed to be promoted by divergent selection. However, none‐liner IBE pattern indicates the effects of salinity were overwhelmed by spatial distance or other ecological processes in certain areas and also suggests that salinity was not the only selective factor driving population differentiation. These results together indicate that geographic distance, salinity difference, and colonization history co‐contributed in shaping the genetic structure of S. mariqueter and that their relative importance was correlated with spatial scale and environment gradient.     

Maintenance of genetic differentiation across a transition zone in the sea: discordance between nuclear and cytoplasmic markers

To investigate the origin and maintenance of the genetic discontinuity between Atlantic and Mediterranean populations of the common sea bass (Dicentrarchus labrax) we analysed the genetic variation at a fragment of mitochondrial cytochrome b sequence for 18 population samples. The result were also compared with new or previously published microsatellite data. Seven mitochondrial haplotypes and an average nucleotidic divergence of 0.02 between Atlantic and Mediterranean populations that matches a Pleistocene allopatric isolation were found. The frequency variation at the cytochrome b locus was many times greater between Atlantic and Mediterranean populations (theta(C) = 0.67) than at microsatellite loci (theta(N)= 0.02). The examination of the different evolutionary forces at play suggests that a sex-biased hybrid breakdown is a likely explanation for part of the observed discrepancy between mitochondrial and nuclear loci. In addition, an analysis is made of the correlation between microsatellite loci points towards the possible existence of a hybrid zone in samples from the Alboran Sea.     

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.     

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.     

Individual assignment test reveals differential restriction to dispersal between two salmonids despite no increase of genetic differences with distance

In many species genes move over limited distances, such that genetic differences among populations or individuals are expected to increase as a function of geographical distance. In other species, however, genes may move any distance over a single generation time, such that no increase of genetic differences is expected to occur with distance. Patterns of gene dispersal have been assessed typically using this theoretical property. In this study, this classical approach based on a Mantel test was compared to a new method using individual assignment to reveal contrasts in dispersal patterns between 15 populations of brook charr Salvelinus fontinalis and 10 populations of Atlantic salmon Salmo salar sampled in eastern Canada, where both species co-occur naturally. Based on the Mantel test, we found evidence for neither an increase of genetic differences with distance in either species nor a significant contrast between them. The individual-based method, in contrast, revealed that individual assignment in both species was non random, being significantly biased toward geographically proximate locations. Furthermore, brook charr were on average assigned to a closer river than were salmon, according to a priori expectations based on the dispersal behaviour of the two species. We thus propose that individual assignment methods might be a promising and more powerful alternative to Mantel tests when isolation by distance cannot be postulated a priori.     

Likelihood and approximate likelihood analyses of genetic structure in a linear habitat: performance and robustness to model mis-specification

We evaluate the performance of maximum likelihood (ML) analysis of allele frequency data in a linear array of populations. The parameters are a mutation rate and either the dispersal rate in a stepping stone model or a dispersal rate and a scale parameter in a geometric dispersal model. An approximate procedure known as maximum product of approximate conditional (PAC) likelihood is found to perform as well as ML. Mis-specification biases may occur because the importance sampling algorithm is formally defined in term of mutation and migration rates scaled by the total size of the population, and this size may differ widely in the statistical model and in reality. As could be expected, ML generally performs well when the statistical model is correctly specified. Otherwise, mutation rate estimates are much closer to mutation probability scaled by number of demes in the statistical model than scaled by number of demes in reality when mutation probability is high and dispersal is most limited. This mis-specification bias actually has practical benefits. However, opposite results are found in opposite conditions. Migration rate estimates show roughly similar trends, but they may not always be easily interpreted as low-bias estimates of dispersal rate under any scaling. Estimation of the dispersal scale parameter is also affected by mis-specification of the number of demes, and the different biases compensate each other in such a way that good estimation of the so-called neighborhood size (or more precisely the product of population density and mean-squared parent-offspring dispersal distance) is achieved. Results congruent with these findings are found in an application to a damselfly data set.     

Fine-scale genetic structure and dispersal in cooperatively breeding apostlebirds

In cooperatively breeding species, restricted dispersal of offspring leads to clustering of closely related individuals, increasing the potential both for indirect genetic benefits and inbreeding costs. In apostlebirds (Struthidea cinerea), philopatry by both sexes results in the formation of large (up to 17 birds), predominantly sedentary breeding groups that remain stable throughout the year. We examined patterns of relatedness and fine-scale genetic structure within a population of apostlebirds using six polymorphic microsatellite loci. We found evidence of fine-scale genetic structure within the study population that is consistent with behavioural observations of short-distance dispersal, natal philopatry by both sexes and restricted movement of breeding groups between seasons. Global F(ST) values among breeding groups were significantly positive, and the average level of pairwise relatedness was significantly higher for individuals within groups than between groups. For individuals from different breeding groups, geographical distance was negatively correlated with pairwise relatedness and positively correlated with pairwise F(ST). However, when each sex was examined separately, this pattern was significant only among males, suggesting that females may disperse over longer distances. We discuss the potential for kin selection to influence the evolution and maintenance of cooperative breeding in apostlebirds. Our results demonstrate that spatial genetic structural analysis offers a useful alternative to field observations in examining dispersal patterns of cooperative breeders.     

Geographic variation of the major histocompatibility complex in Eastern Atlantic grey seals (Halichoerus grypus)

Pathogen-driven balancing selection maintains high genetic diversity in many vertebrates, particularly in the major histocompatibility complex (MHC) immune system gene family, which is often associated with disease susceptibility. In large natural populations where subpopulations face different pathogen pressures, the MHC should show greater genetic differentiation within a species than neutral markers. We examined genetic diversity at the MHC-DQB locus and nine putatively neutral microsatellite markers in grey seals (Halichoerus grypus) from eight United Kingdom (UK) colonies, the Faeroe Islands and Sable Island, Canada. Five DQB alleles were identified in grey seals, which varied in prevalence across the grey seal range. Among the seal colonies, significant differences in DQB allele and haplotype frequencies and in average DQB heterozygosity were observed. Additionally, the DQB gene exhibited greater differentiation among colonies compared with neutral markers, yet a weaker pattern of isolation by distance (IBD). After correcting for the underlying IBD pattern, subpopulations breeding in similar habitats were more similar to one another in DQB allele frequencies than populations breeding in different habitats, but the same did not hold true for microsatellites, suggesting that habitat-specific pathogen pressure influences MHC evolution. Overall, the data are consistent with selection at MHC-DQB loci in grey seals with both varying selective pressures and geographic population structure appearing to influence the DQB genetic composition of breeding colonies.     

Microsatellite analysis of genetic variation among and within Alpine marmot populations in the French Alps

The genetic structure of the Alpine marmot, Marmota marmota, was studied by an analysis of five polymorphic microsatellite loci. Eight locations were sampled in the French Alps, one from Les Ecrins valley (n = 160), another from La Sassière valley (n = 289) and the six others from the Maurienne valley (n = 139). Information on social group structure was available for both Les Ecrins and La Sassière but not for the other samples. The high levels of genetic diversity observed are at odds with the results obtained using microsatellites, minisatellites and allozymes on Alpine marmots from Germany, Austria and Switzerland. Strong deficits in heterozygotes were found in Les Ecrins and La Sassière. They are caused by a Wahlund effect due to the family structure (i.e. differentiation between the family groups). The family groups exhibit excess of heterozygotes rather than deficits. This may be caused by outbreeding and this is compatible with recent results from the genetics of related social species when information on the social structure is taken into account. The observed outbreeding could be the result of females mating with transient males or males coming from neighbouring colonies. Both indicate that the species may not be as monogamous as is usually believed. The results are also compatible with a male-biased dispersal but do not allow us to exclude some female migration. We also found a significant correlation between geographical and genetic distance indicating that isolation by distance could be an issue in marmots. This study is the first that analysed populations of marmots taking into account the social structure within populations and assessing inbreeding at different levels (region, valley, population, and family groups). Our study clearly demonstrated that the sampling strategy and behavioural information can have dramatic effects on both the results and interpretation of the genetic data.     

Spatial population genetic structure and limited dispersal in a Rocky Mountain alpine stream insect

Using the mitochondrial cytochrome oxidase I (COI) gene, we assessed the phylogeographic structure of Prosimulium neomacropyga, a black fly (Simuliidae) whose distribution in the US Southern Rockies ecoregion is limited to alpine tundra streams. Given high habitat specificity, lack of hydrological connection between streams, and a terrestrial environment restrictive to insect flight, we hypothesized limited gene flow. A spatially nested sampling design showed that grouping populations according to high-elevation 'islands' of alpine tundra (which typically include headwater streams of > 1 watershed) explained a significant proportion of genetic variation while grouping streams according to major watershed (across islands) did not. Nested clade analysis and isolation-by-distance (IBD) relationships further implicated limited ongoing gene flow within but not among the isolated alpine islands. IBD was strong among five streams within an individual island using each of four alternative models of pairwise landscape connectivity for flying insects. Results of all landscape models were positively correlated, suggesting that straight-line distance is an acceptable surrogate for presumably more biologically meaningful connectivity measures in this system. IBD was significantly weaker across the entire study area, comprised of three separate islands. Overall, population structure was significant with F(ST) = 0.38, suggesting limited dispersal across a small spatial extent.     

The sampling scheme matters: Pan troglodytes troglodytes and P. t. schweinfurthii are characterized by clinal genetic variation rather than a strong subspecies break

Populations of an organism living in marked geographical or evolutionary isolation from other populations of the same species are often termed subspecies and expected to show some degree of genetic distinctiveness. The common chimpanzee (Pan troglodytes) is currently described as four geographically delimited subspecies: the western (P. t. verus), the nigerian‐cameroonian (P. t. ellioti), the central (P. t. troglodytes) and the eastern (P. t. schweinfurthii) chimpanzees. Although these taxa would be expected to be reciprocally monophyletic, studies have not always consistently resolved the central and eastern chimpanzee taxa. Most studies, however, used data from individuals of unknown or approximate geographic provenance. Thus, genetic data from samples of known origin may shed light on the evolutionary relationship of these subspecies. We generated microsatellite genotypes from noninvasively collected fecal samples of 185 central chimpanzees that were sampled across large parts of their range and analyzed them together with 283 published eastern chimpanzee genotypes from known localities. We observed a clear signal of isolation by distance across both subspecies. Further, we found that a large proportion of comparisons between groups taken from the same subspecies showed higher genetic differentiation than the least differentiated between‐subspecies comparison. This proportion decreased substantially when we simulated a more clumped sampling scheme by including fewer groups. Our results support the general concept that the distribution of the sampled individuals can dramatically affect the inference of genetic population structure. With regard to chimpanzees, our results emphasize the close relationship of equatorial chimpanzees from central and eastern equatorial Africa and the difficult nature of subspecies definitions. Am J Phys Anthropol 156:181–191, 2015. © 2014 Wiley Periodicals, Inc.     

Microsatellite DNA polymorphism confirms reproductive isolation and reveals differences in population genetic structure of cryptic pipistrelle bat species

Previous studies have indicated that the common European pipistrelle bat ( Pipistrellus pipistrellus ) comprises two cryptic species, P. pipistrellus and Pipistrellus pygmaeus , which differ in echolocation call frequency and mitochondrial DNA sequence. However, levels of divergence based on nuclear markers have not been examined, and hence the potential for male-mediated gene flow between the species cannot be discounted. Moreover, little is known about population structure and migration patterns in either species. Here, we describe the use of microsatellites to investigate nuclear DNA differentiation between, and the pattern of population genetic structure within, the two cryptic pipistrelle species. In total, 1300 individuals from 82 maternity colonies were sampled across the British Isles and Continental Europe. We show, using multivariate analyses, that colonies of the same species are generally genetically more similar to each other than to those from the other species regardless of geographical location. Our findings support the hypothesis that the species are reproductively isolated. Significant patterns of genetic isolation by distance were identified in both species, indicating that mating may occur before any long-distance autumnal migration. The presence of a sea channel does not confer higher levels of genetic differentiation among colonies over and above distance alone in either species. Differences in genetic population structure were identified between the species, with P. pipistrellus showing a wider range of levels of genetic differentiation among colonies and a stronger relationship between genetic and geographical distance than P. pygmaeus . Differences in dispersal, mating behaviour, colony size and/or postglacial colonization patterns could contribute to the differences observed.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 90 , 539–550.     

Evolutionary and demographic processes shaping geographic patterns of genetic diversity in a keystone species,the African forest elephant (Loxodonta cyclotis)

The past processes that have shaped geographic patterns of genetic diversity may be difficult to infer from current patterns. However, in species with sex differences in dispersal, differing phylogeographic patterns between mitochondrial (mt) and nuclear (nu) DNA may provide contrasting insights into past events. Forest elephants (Loxodonta cyclotis) were impacted by climate and habitat change during the Pleistocene, which likely shaped phylogeographic patterns in mitochondrial (mt) DNA that have persisted due to limited female dispersal. By contrast, the nuclear (nu) DNA phylogeography of forest elephants in Central Africa has not been determined. We therefore examined the population structure of Central African forest elephants by genotyping 94 individuals from six localities at 21 microsatellite loci. Between forest elephants in western and eastern Congolian forests, there was only modest genetic differentiation, a pattern highly discordant with that of mtDNA. Nuclear genetic patterns are consistent with isolation by distance. Alternatively, male‐mediated gene flow may have reduced the previous regional differentiation in Central Africa suggested by mtDNA patterns, which likely reflect forest fragmentation during the Pleistocene. In species like elephants, male‐mediated gene flow erases the nuclear genetic signatures of past climate and habitat changes, but these continue to persist as patterns in mtDNA because females do not disperse. Conservation implications of these results are discussed.     

Temperature and precipitation,but not geographic distance,explain genetic relatedness among populations in the perennial grass Festuca rubra

Aims Knowledge of genetic structure of natural populations and its determinants may provide key insights into the ability of species to adapt to novel environments. In many genetic studies, the effects of climate could not be disentangled from the effects of geographic proximity. We aimed to understand the effects of temperature and moisture on genetic diversity of populations and separate these effects from the effects of geographic distance. We also wanted to explore the patterns of distribution of genetic diversity in the system and assess the degree of clonality within the populations. We also checked for possible genome size variation in the system.     

Old fish in a young lake: stone loach (Pisces: Barbatula barbatula) populations in Lake Constance are genetically isolated by distance

The genetic structure of 10 populations (453 individuals) of stone loach (Barbatula barbatula L.), a small bottom-dwelling cyprinid fish, in the littoral zone of Lake Constance, central Europe, was investigated by analysing the mitochondrial control region sequences and five microsatellite loci. An unexpectedly high degree of genetic diversity (up to 0.36%) and old estimated age of these populations (> 150 000 years) based on mitochondrial DNA (mtDNA) was found. These findings contrast with the relatively young age of the lake, which could be colonized by fish only after the last ice age around 15 000 bp. Stone loach appears to be an old species in a young lake. Both types of molecular markers showed population genetic structure pronounced in mtDNA (overall F(ST) = 0.15) but moderate in microsatellites (F(ST) = 0.03). As predicted by its life history, philopatry, and limited capacity for dispersal, stone loach populations of Lake Constance show a clear pattern of isolation by distance. Geographic distances along the shores are the best explanation for the observed geographical distribution of genetic differentiation (r = 0.88), indicating that open water represents a barrier for the dispersal of the stone loach. The colonization of Lake Constance might have occurred initially at one location and then populations spread throughout the lake in a stepwise manner following the shoreline, and subsequently remained largely genetically isolated as suggested by the large observed differences among them.     

Environmental correlates of population differentiation in Atlantic herring

The marine environment is characterized by few physical barriers, and pelagic fishes commonly show high migratory potential and low, albeit in some cases statistically significant, levels of genetic divergence in neutral genetic marker analyses. However, it is not clear whether low levels of differentiation reflect spatially separated populations experiencing gene flow or shallow population histories coupled with limited random genetic drift in large, demographically isolated populations undergoing independent evolutionary processes. Using information for nine microsatellite loci in a total of 1951 fish, we analyzed genetic differentiation among Atlantic herring from eleven spawning locations distributed along a longitudinal gradient from the North Sea to the Western Baltic. Overall genetic differentiation was low (theta = 0.008) but statistically significant. The area is characterized by a dramatic shift in hydrography from the highly saline and temperature stable North Sea to the brackish Baltic Sea, where temperatures show high annual variation. We used two different methods, a novel computational geometric approach and partial Mantel correlation analysis coupled with detailed environmental information from spawning locations to show that patterns of reproductive isolation covaried with salinity differences among spawning locations, independent of their geographical distance. We show that reproductive isolation can be maintained in marine fish populations exhibiting substantial mixing during larval and adult life stages. Analyses incorporating genetic, spatial, and environmental parameters indicated that isolating mechanisms are associated with the specific salinity conditions on spawning locations.     

Widespread evidence for incipient ecological speciation: a meta‐analysis of isolation‐by‐ecology

Ecologically mediated selection has increasingly become recognised as an important driver of speciation. The correlation between neutral genetic differentiation and environmental or phenotypic divergence among populations, to which we collectively refer to as isolation‐by‐ecology (IBE), is an indicator of ecological speciation. In a meta‐analysis framework, we determined the strength and commonality of IBE in nature. On the basis of 106 studies, we calculated a mean effect size of IBE with and without controlling for spatial autocorrelation among populations. Effect sizes were 0.34 (95% CI 0.24–0.42) and 0.26 (95% CI 0.13–0.37), respectively, indicating that an average of 5% of the neutral genetic differentiation among populations was explained purely by ecological contrast. Importantly, spatial autocorrelation reduced IBE correlations for environmental variables, but not for phenotypes. Through simulation, we showed how the influence of isolation‐by‐distance and spatial autocorrelation of ecological variables can result in false positives or underestimated correlations if not accounted for in the IBE model. Collectively, this meta‐analysis showed that ecologically induced genetic divergence is pervasive across time‐scales and taxa, and largely independent of the choice of molecular marker. We discuss the importance of these results in the context of adaptation and ecological speciation and suggest future research avenues.