Whole genome analyses reveal weak signatures of population structure and environmentally associated local adaptation in an important North American pollinator, the bumble bee Bombus vosnesenskii

Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history, and environmentally associated genomic variation in Bombus vosnesenskii, a common bumble bee in the western U.S.A., using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find that B. vosnesenskii exhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent (SMC) models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size (Ne) trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables, and environmental association analysis with latent factor mixed modeling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared to other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicate Bombus vosnesenskii is an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.