Landscape and climate-associated selection in the native and widespread bumblebee, Bombus terrestris

Anthropogenic land-use and climate change pose novel selection pressures on bees, yet their evolutionary responses in terms of morphological or physiological adaptations remain unclear. While adaptive responses are expected, these may be constrained by gene flow when changes in selection pressures are spatially heterogeneous. The buff-tailed bumblebee (Bombus terrestris) is a widespread species that copes well with anthropogenic land-use and climate change, suggesting high adaptive capacity or phenotypic plasticity. Here, we genotyped populations of native B. terrestris in south and central Sweden using RADseq to investigate genetic structure and local adaptation across a paired design of agricultural landscapes with high and low land-use complexity along a geographic climate gradient. We expected to find genetic structure reflective of regional barriers to gene flow, and molecular evidence for local adaptation to differing landscape and climate conditions. We found genetic structure separating southern Sweden from more northern regions, with a negative Tajima's D indicating a potential population expansion, likely northwards and inland into forested areas, consistent with observational data indicating a range shift. We found weak but significant evidence for local adaptation to climate and land use, specifically to agricultural land cover, including genes under putative selection linked to insecticide resistance. Signatures of selection were also identified in relation to latitude, temperature, and urban land cover, with other candidate SNPs associated with olfaction and immune response. Our results suggest that B. terrestris successfully responded to anthropogenic land-use and climate changes, likely due to its generalist traits, enabling phenotypic adaptation to changing environments.