Unraveling the habitat preferences of two closely related bumble bee species in Eastern Europe.

Co-occurrence of closely related species is often explained through resource partitioning, where key morphological or life-history traits evolve under strong divergent selection. In bumble bees (genus ), differences in tongue lengths, nest sites, and several life-history traits are the principal factors in resource partitioning. However, the buff-tailed and white-tailed bumble bee ( and respectively) are very similar in morphology and life history, but their ranges nevertheless partly overlap, raising the question how they are ecologically divergent. What little is known about the environmental factors determining their distributions stems from studies in Central and Western Europe, but even less information is available about their distributions in Eastern Europe, where different subspecies occur. Here, we aimed to disentangle the broad habitat requirements and associated distributions of these species in Romania and Bulgaria. First, we genetically identified sampled individuals from many sites across the study area. We then not only computed species distributions based on presence-only data, but also expanded on these models using relative abundance data. We found that is a more generalist species than previously thought, but that is restricted to forested areas with colder and wetter climates, which in our study area are primarily found at higher elevations. Both vegetation parameters such as annual mean Leaf Area Index and canopy height, as well as climatic conditions, were important in explaining their distributions. Although our models based on presence-only data suggest a large overlap in their respective distributions, results on their relative abundance suggest that the two species replace one another across an environmental gradient correlated to elevation. The inclusion of abundance enhances our understanding of the distribution of these species, supporting the emerging recognition of the importance of abundance data in species distribution modeling.