A clever solution to a vexing problem

F_ST (as well as related measures such as G_ST) has long been used both as a measure of the relative amount of genetic variation between populations and as an indicator of the amount of gene flow among populations. Unfortunately, F_ST and its clones are also sensitive to mutation, particularly when the mutation rate per locus is greater than the migration rate among populations. Relatively high mutation rates cause estimates of F_ST and G_ST to be much lower than researchers sometimes expect, when migration rates are low in the studied species. Several recent suggestions for dealing with this problem have been unsatisfactory for one reason or another, and no general solution exists (if we are not to abandon these otherwise useful measures of differentiation). In an important article in this issue, Jinliang Wang (2015) shows that it is possible to identify whether the genetic markers in a given study are likely to give estimates of F_ST that are strongly affected by mutation. The proposed test is simple and elegant, and with it, molecular ecologists can determine whether the F_ST from their makers can be depended on for further inference about their species’ genome and the demographic forces which shaped its patterns.     

Hybridization in natural sympatric populations of Dermacentor ticks in northwestern North America

Hybridization in ticks has been described in a handful of species and mostly as a result of laboratory experiments. We used 148 AFLP loci to describe putative hybridization events between D. andersoni and D. variabilis in sympatric populations from northwestern North America. Recently, D. variabilis has expanded its range westward into the natural range of D. andersoni. Using a sample of 235 D. andersoni and 62 D. variabilis, we identified 31 individuals as putative hybrids: four F₂ individuals and 27 backcrosses to D. andersoni (as defined by NewHybrids ). We found no evidence of hybrids backcrossing into D. variabilis. Furthermore, all hybrids presented 16S mtDNA signatures characteristic of D. andersoni, which indicates the directionality of the hybrid crosses: female D. andersoni × male D. variabilis. We also discovered 13 species‐specific AFLP fragments for D. andersoni. These loci were found to have a decreased occurrence in the putative hybrids and were absent altogether in D. variabilis samples. AFLP profiles were also used to determine the levels of genetic population structure and gene flow among nine populations of D. andersoni and three of D. variabilis. Genetic structure exists in both species (D. andersoni, Φ_ST = 0.110; D. variabilis, Φ_ST = 0.304) as well as significant estimates of isolation by distance (D. andersoni, ρ = 0.066, P = 0.001; D. variabilis, ρ = 0.729, P = 0.001).     

Population genomic analyses from low‐coverage RAD‐Seq data: a case study on the non‐model cucurbit bottle gourd

Restriction site‐associated DNA sequencing (RAD‐Seq), a next‐generation sequencing‐based genome ‘complexity reduction’ protocol, has been useful in population genomics in species with a reference genome. However, the application of this protocol to natural populations of genomically underinvestigated species, particularly under low‐to‐medium sequencing depth, has not been well justified. In this study, a Bayesian method was developed for calling genotypes from an F₂ population of bottle gourd [Lagenaria siceraria (Mol.) Standl.] to construct a high‐density genetic map. Low‐depth genome shotgun sequencing allowed the assembly of scaffolds/contigs comprising approximately 50% of the estimated genome, of which 922 were anchored for identifying syntenic regions between species. RAD‐Seq genotyping of a natural population comprising 80 accessions identified 3226 single nuclear polymorphisms (SNPs), based on which two sub‐gene pools were suggested for association with fruit shape. The two sub‐gene pools were moderately differentiated, as reflected by the Hudson's F_ST value of 0.14, and they represent regions on LG7 with strikingly elevated F_ST values. Seven‐fold reduction in heterozygosity and two times increase in LD (r²) were observed in the same region for the round‐fruited sub‐gene pool. Outlier test suggested the locus LX3405 on LG7 to be a candidate site under selection. Comparative genomic analysis revealed that the cucumber genome region syntenic to the high F_ST island on LG7 harbors an ortholog of the tomato fruit shape gene OVATE. Our results point to a bright future of applying RAD‐Seq to population genomic studies for non‐model species even under low‐to‐medium sequencing efforts. The genomic resources provide valuable information for cucurbit genome research.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white‐nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292‐bp region of the mitochondrial genome. Analyses of F_ST, Φ_ST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F_ST = 0.001, P < 0.05, Global Φ_ST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male‐biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.     

QST–FST comparisons with unbalanced half‐sib designs

Q_ST, a measure of quantitative genetic differentiation among populations, is an index that can suggest local adaptation if Q_ST for a trait is sufficiently larger than the mean F_ST of neutral genetic markers. A previous method by Whitlock and Guillaume derived a simulation resampling approach to statistically test for a difference between Q_ST and F_ST, but that method is limited to balanced data sets with offspring related as half‐sibs through shared fathers. We extend this approach (i) to allow for a model more suitable for some plant populations or breeding designs in which offspring are related through mothers (assuming independent fathers for each offspring; half‐sibs by dam); and (ii) by explicitly allowing for unbalanced data sets. The resulting approach is made available through the R package QstFstComp.     

Pattern‐oriented modelling of population genetic structure

Although several statistical approaches can be used to describe patterns of genetic variation and infer stochastic differentiation, selective responses, or interruptions of gene flow due to physical or environmental barriers, it is worthwhile to note that similar processes, controlled by several parameters in theoretical models, frequently give rise to similar patterns. Here, we develop a Pattern‐Oriented Modelling (POM) approach that allows us to determine how a complex set of parameters potentially driving empirical genetic differentiation among populations generate alternative scenarios that can be fitted to observed data. We generated 10 000 random combinations of parameters related to population size, gene flow and response to gradients (both driven by dispersal and selection) in a spatially explicit model, and analysed simulated patterns with F_ST statistics and mean correlograms using Moran's I spatial autocorrelation coefficients. These statistics were compared with observed patterns for a tree species endemic to the Brazilian Cerrado. For a best match with observed F_ST (equal to 0.182), the important parameters driving simulated scenario are mainly related to population structure, including low population size with closed populations (low N_m), strong distance decay of gene flow, in addition to a strong effect of the initial variance of allele frequencies. These scenarios present a low autocorrelation of allele frequencies. Best matching of correlograms, on the other hand, appears in simulations with a large population size, high N_m and low population differentiation and F_ST (as well as more gene flow). Thus, targeting the two statistics (correlograms and F_ST) shows that best matches with empirical data with two distinct sets of parameters in the simulations, because observed patterns involve both a relatively high F_ST and significant autocorrelation. This conflict can be resolved by assuming that initial variance in allele frequencies can be interpreted as reflecting deep‐time historical variation and evolutionary dynamics of allele frequencies, creating a relatively high level of population differentiation, whereas current patterns in gene flow creates spatial autocorrelation. This make sense in terms of the previous knowledge on population differentiation in D. alata, especially if patterns are explained by a combination of isolation‐by‐distance and allelic surfing due to range expansion after the last glacial maximum. This reveals the potential for more complex applications of POM in population genetics. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 1152–1161.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage F_ST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

High gene flow despite opposite chirality in hybrid zones between enantiomorphic door snails

We studied differentiation and geneflow patterns between enantiomorphic door‐snail species in two hybrid zones in the Bucegi Mountains (Romania) to investigate the effects of intrinsic barriers (complications in copulation) and extrinsic selection by environmental factors. A mitochondrial gene tree confirmed the historical separation of the examined populations into the dextral Alopia livida and the sinistral Alopia straminicollis in accordance with the morphological classification, but also indicated gene flow between the species. By contrast, a network based on amplified fragment length polymorphisms (AFLP) markers revealed local groups of populations as units independent of their species affiliation. Admixture analyses based on AFLP data showed that the genomes of most individuals in the hybrid zones are composed of parts of the genomes of both parental taxa. The introgression patterns of a notable fraction of the examined markers deviated from neutral introgression. However, the patterns of most non‐neutral markers were not concordant between the two hybrid zones. There was also no concordance between non‐neutral markers in the two genomic clines and markers that were correlated with environmental variables or markers that were correlated with the proportion of dextral individuals in the populations. Neither extrinsic selection by environmental factors nor intrinsic barriers resulting from positive frequency‐dependent selection of the prevailing coiling direction were sufficient to maintain the distinctness of A. straminicollis and A. livida. Despite being historically separated units, we conclude that these taxa now merge where they come into contact.     

Positive selection of IL‐33 in adaptive immunity of domestic Chinese goats

The identification of the candidate genes that play key role in phenotypic variation in livestock populations can provide new information about evolution and positive selection. IL‐33 (71954) (Interleukin) gene is associated with the increased nematode resistance in small ruminants; however, the role of IL‐33 for the genetic control of different diseases in Chinese goat breeds is poorly described in scientific literature. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single‐nucleotide polymorphism in IL‐33 gene. Fixation Index (F_ST)‐based method was used for the outlier loci determination and found that IL‐33 was present in outlier area with the provisional combined allocation of mean heterozygosity and F_ST. Positively selected IL‐33 gene was significantly, that is, p(Simul F_ST < sample F_ST = 0.98*) present in corresponding positive selection area. Hence, our study provided novel information about the nucleotide variations in IL‐33 gene and found to be nonsynonymous which may helpful for the genetic control of diseases by enhancing the immune system in local Chinese goat breeds as well as in other analyzed vertebrate species.     

Development of genetic diversity,differentiation and structure over 500 years in four ponderosa pine populations

Population history plays an important role in shaping contemporary levels of genetic variation and geographic structure. This is especially true in small, isolated range‐margin populations, where effects of inbreeding, genetic drift and gene flow may be more pronounced than in large continuous populations. Effects of landscape fragmentation and isolation distance may have implications for persistence of range‐margin populations if they are demographic sinks. We studied four small, disjunct populations of ponderosa pine over a 500‐year period. We coupled demographic data obtained through dendroecological methods with microsatellite data to discern how and when contemporary levels of allelic diversity, among and within‐population levels of differentiation, and geographic structure, arose. Alleles accumulated rapidly following initial colonization, demonstrating proportionally high levels of gene flow into the populations. At population sizes of approximately 100 individuals, allele accumulation saturated. Levels of genetic differentiation among populations (F_ST and Jost's D_est) and diversity within populations (F_IS) remained stable through time. There was no evidence of geographic genetic structure at any time in the populations' history. Proportionally, high gene flow in the early stages of population growth resulted in rapid accumulation of alleles and quickly created relatively homogenous genetic patterns among populations. Our study demonstrates that contemporary levels of genetic diversity were formed quickly and early in population development. How contemporary genetic diversity accumulates over time is a key facet of understanding population growth and development. This is especially relevant given the extent and speed at which species ranges are predicted to shift in the coming century.     

Separating Effects of Gene Flow and Natural Selection along an Environmental Gradient

Genetic differentiation along environmental clines is often observed as a result of interplay between gene flow and natural selection. In order to understand the relative roles of these processes in shaping this differentiation, we designed a study in which we used two approaches that have not previously been combined, the Q _ST–F _ST comparison and crossbreeding. We examined (1) interpopulation phenotypic and genetic (AFLP) variation, and (2) performance of interpopulation hybrids in a common annual Senecio glaucus. Fitness of interpopulation hybrids (F1 and F2) was tested under simulated population natural conditions in terms of aridity and analyzed for a relationship with (1) spatial distance and (2) environmental differences (amount of annual rainfall). While phenotypic variation corresponded to the clinal changes in aridity along population locations, viz. narrower and longer leaves, longer leaf outgrowths and advanced flowering in more arid environments, the F _ST < 0.1 calculated from AFLP data suggested intensive interpopulation gene flow, with little if any contribution of genetic drift. Performance of hybrids in simulated natural environments revealed heterosis in F1, but a hybrid breakdown in F2 generation. These effects were related to both the spatial distance between hybrid parents and their population rainfall differences. The detected clinal phenotypic variation and outbreeding depression in F2 strongly support presence of aridity-induced clinal natural selection, which is matched by the observed Q _ST ≫ F _ST. From this we conclude that Q _ST–F _ST comparison can detect effect of diversifying selection when patterns of phenotypic variation across sampled locations can be reliably predicted from environmental variation.     

Historical vicariance and male-mediated gene flow in the toad-headed lizards Phrynocephalus przewalskii

Using mitochondrial and microsatellite DNA data and a population genetic approach, we tested male‐mediated gene flow in the toad‐headed lizards Phrynocephalus przewalskii. The mitochondrial DNA (ND2 gene), on the one hand, revealed two major lineages and a strong population genetic structure (F_ST = 0.692; F_ST′ = 0.995). The pairwise differences between the two lineages ranged from 2.1% to 6.4% and the geographical division of the two lineages coincided with a mountain chain consisting of the Helan and Yin Mountains, suggesting a historical vicariant pattern. On the other hand, the nuclear microsatellite DNA revealed a significant but small population genetic structure (F_ST = 0.017; F_ST′ = 0.372). The pairwise F_ST among the nine populations examined with seven microsatellite DNA loci ranged from 0.0062 to 0.0266; the assignment test failed to detect any naturally occurring population clusters. Furthermore, the populations demonstrated a weak isolation by distance and a northeast to southwest clinal variation, rather than a vicariant pattern. A historical vicariant event followed by male‐mediated gene flow appears to be the best explanation for the data. Approximately 2–5 Ma, climatic change may have created an uninhabitable zone along the Helan‐Yin mountain chain and initiated the divergence between the two mitochondrial lineages. With further climatic changes, males were able to disperse across the mountain chain, causing sufficient gene flow that eventually erased the vicariant pattern and drastically reduced the population genetic structure, while females remained philopatric and maintained the mitochondrial DNA (mtDNA) divergence. Although polygyny mating system and female philopatry may partially contribute to the reduced movement of females, other hypotheses, such as female intrasexual aggression, should also be explored.     

Isolated in an ocean of grass: low levels of gene flow between termite subpopulations

Habitat fragmentation is one of the most important causes of biodiversity loss, but many species are distributed in naturally patchy habitats. Such species are often organized in highly dynamic metapopulations or in patchy populations with high gene flow between subpopulations. Yet, there are also species that exist in stable patchy habitats with small subpopulations and presumably low dispersal rates. Here, we present population genetic data for the ‘magnetic’ termite Amitermes meridionalis, which show that short distances between subpopulations do not hinder exceptionally strong genetic differentiation (F_ST: 0.339; R_ST: 0.636). Despite the strong genetic differentiation between subpopulations, we did not find evidence for genetic impoverishment. We propose that loss of genetic diversity might be counteracted by a long colony life with low colony turnover. Indeed, we found evidence for the inheritance of colonies by so‐called ‘replacement reproductives’. Inhabiting a mound for several generations might result in loss of gene diversity within a colony but maintenance of gene diversity at the subpopulation level.     

Walk,swim or fly? Locomotor mode predicts genetic differentiation in vertebrates

Limited dispersal is commonly used to explain differences in diversification rates. An obvious but unexplored factor affecting dispersal is the mode of locomotion used by animals. Whether individuals walk, swim or fly can dictate the type and severity of geographical barriers to dispersal, and determine the general range over which genetic differentiation might occur. We collated information on locomotion mode and genetic differentiation (F_ST) among vertebrate populations from over 400 published articles. Our results showed that vertebrate species that walk tend to have higher genetic differentiation among populations than species that swim or fly. Within species that swim, vertebrates in freshwater systems have higher genetic differentiation than those in marine systems, which is consistent with the higher number of species in freshwater environments. These results show that locomotion mode can impact gene flow among populations, supporting at a broad‐scale what has previously been proposed at smaller taxonomical scales.     

Contemporary human‐altered landscapes and oceanic barriers reduce bumble bee gene flow

Much of the world's terrestrial landscapes are being altered by humans in the form of agriculture, urbanization and pastoral systems, with major implications for biodiversity. Bumble bees are one of the most effective pollinators in both natural and cultivated landscapes, but are often the first to be extirpated in human‐altered habitats. Yet, little is known about the role of natural and human‐altered habitats in promoting or limiting bumble bee gene flow. In this study, I closely examine the genetic structure of the yellow‐faced bumble bee, Bombus vosnesenskii, across the southwestern US coast and find strong evidence that natural oceanic barriers, as well as contemporary human‐altered habitats, limit bee gene flow. Heterozygosity and allelic richness were lower in island populations, while private allelic richness was higher in island populations compared to mainland populations. Genetic differentiation, measured for three indices across the 1000 km study region, was significantly greater than the null expectation (F_ST = 0.041, F’_ST = 0.044 and D_est = 0.155) and correlated with geographic distance. Furthermore, genetic differentiation patterns were most strongly correlated with contemporary (2011) not past (2006, 2001) resistance maps calibrated for high dispersal limitation over oceans, impervious habitat and croplands. Despite the incorporation of dramatic elevation gradients, the analyses reveal that oceans and contemporary human land use, not mountains, are the primary dispersal barriers for B. vosnesenskii gene flow. These findings reinforce the importance of maintaining corridors of suitable habitat across the distribution range of native pollinators to promote their persistence and safeguard their ability to provide essential pollination services.