Table S1

Microsatellite loci in Bombus terrestris typed in multiplex reactions in this study. Table contains list of primers for respective loci as well as information on the labelling dye used

Table S4

Table S4 – STRUCTURE assignments of population probability to membership of K = 3 clusters and K = 7 clusters

Paup Tree File

Paup output file

Figure S2

Figure S2– A) Estimated likelihood of the posterior probability, LnP(X|K), for values of K (the number of genetic clusters) ranging from 1 to 12 as calculated using STRUCTURE. The mean LnP(X| K) for each K over 20 permutations is represented by solid squares. Bars correspond to standard deviations over 20 permutations. B) modal value quantity of ∆K based on the second order rate change is calculated according to Evanno et al (2005)

MrBayes Tree File

MrBayes output file

COI aligned sequences for all individuals_PAUP

COI aligned sequences for all individuals. PAUP input file

Table 2 - GPS locations of all sampled populations

GPS locations of all sampled population

Table S5

Table S5 – Correlations between genetic and geographic distance, as measured by Rousset’s distance (Fst/1-Fst) and the natural logarithm of geographical separation for Bombus terrestri

Table 1

Abundance of bumblebee species (Bombus s.str) in assemblages sampled throughout Europe. Species were separated according to the RFLP-based method of Murray et al. (2008)

Code used for Bayesian analysis

Code used for Bayesian analysi