Landscape heterogeneity predicts gene flow in a widespread polymorphic bumble bee, Bombus bifarius (Hymenoptera: Apidae)

Bombus bifarius is a widespread bumble bee that occurs in montane regions of western North America. This species has several major color pattern polymorphisms and shows evidence of genetic structuring among regional populations, and the taxonomic status of regional populations has repeatedly been debated. We test whether observed structure is evidence for discrete gene flow barriers that might indicate isolation or instead reflects clinal variation associated with spatially limited dispersal in a complex landscape. We first consider color pattern variation and identify geographical patterns of B. bifarius color variation using cluster analysis. We then use climate data and a comprehensive set of B. bifarius natural history records with an existing genetic data set to model the distribution of environmentally suitable habitat in western North America and predict pathways of potential gene flow using circuit theory. Resistance distances among populations that incorporate environmental suitability information predict patterns of genetic structure much better than geographic distance or Bayesian clustering alone. Results suggest that there may not be barriers to gene flow warranting further taxonomic considerations, but rather that the arrangement of suitable habitat at broad scales limits dispersal sufficiently to explain observed levels of population differentiation in B. bifarius.