Table S1 Pairwise FST values for the BTM populations based on the COI sequence dataset.

1 / 2 / 3 / 4 / 5 / 6 / 7 / 8 / 9 / 10 / 11 / 12 / 13
1:UK_Wig_Rub
2:UK_Wig_Pc / 0.015
3:US_Mai_Pru / 0.322 / 0.350
4:UK_Wig_Cm / 0.007 / 0.007 / 0.217
5:FR_Mon_Au / 0.364 / 0.384 / 0.815 / 0.428
6:SP_Gra_Au / 0.158 / 0.176 / 0.563 / 0.214 / 0.370
7:SP_Uix_Au / 0.600 / 0.619 / 0.951 / 0.666 / 0.121 / 0.584
8:SP_Men_Au / 0.171 / 0.228 / 0.708 / 0.302 / 0.518 / 0.280 / 0.768
9:SP_Gra_Qfb / 0.079 / 0.098 / 0.501 / 0.135 / 0.169 / 0.111 / 0.400 / 0.201
10:SP_Gra_Qf / 0.175 / 0.193 / 0.557 / 0.230 / 0.059 / 0.194 / 0.244 / 0.294 / 0.018
11:FR_Bri_Qr / 0.005 / 0.046 / 0.217 / -0.012 / 0.466 / 0.264 / 0.742 / 0.277 / 0.167 / 0.247
12:SP_Gra_St / 0.436 / 0.455 / 0.874 / 0.502 / -0.029 / 0.436 / 0.058 / 0.600 / 0.237 / 0.107 / 0.554
13:SP_Ter_Um / 0.240 / 0.260 / 0.708 / 0.302 / -0.005 / 0.256 / 0.241 / 0.382 / 0.060 / -0.009 / 0.326 / 0.032
14:RO_Bot_Rc / 0.043 / 0.096 / 0.479 / 0.122 / 0.446 / 0.230 / 0.682 / 0.067 / 0.151 / 0.245 / 0.068 / 0.519 / 0.320

Significant results (P<0.05) are in bold. Negative FST values indicate a greater genetic difference between individuals within a population than between individuals from different populations, i.e. more heterozygotes than expected.

Table S2 Pairwise FST values for the BTM populations based on the EF1a sequence dataset.

1 / 2 / 3 / 4 / 5 / 6 / 7
1:SP_Gra_Qfb
2:UK_Wig_Cm / -0.008
3:FR_Mon_Au / 0.033 / -0.022
4:SP_Gra_Au / -0.028 / 0.028 / 0.088
5:SP_Uix_Au / 0.095 / 0.030 / 0.000 / 0.174
6:US_Mai_Pru / -0.002 / -0.032 / -0.021 / 0.035 / 0.019
7:UK_Wig_Pc / 0.058 / 0.000 / 0.000 / 0.122 / 0.000 / -0.004
8:RO_Bot_Rc / -0.006 / 0.056 / 0.121 / -0.023 / 0.213 / 0.064 / 0.156

Significant results (P<0.05) are in bold. Negative FST values indicate a greater genetic difference between individuals within a population than between individuals from different populations, i.e. more heterozygotes than expected.

Table S3 Pairwise FST values for the BTM populations based on the WING sequence dataset.

1 / 2 / 3 / 4 / 5 / 6 / 7 / 8 / 9 / 10
1:UK_Wig_Pc
2:UK_Wig_Cm / 0.000
3:UK_Wig_Rub / 0.185 / 0.156
4:SP_Gra_Au / 0.077 / 0.042 / 0.228
5:SP_Gra_Qf / 0.047 / 0.012 / 0.055 / 0.045
6:SP_Gra_Qfb / 0.038 / -0.004 / 0.210 / 0.104 / 0.068
7:SP_Gra_St / 0.100 / 0.073 / 0.314 / -0.018 / 0.091 / 0.125
8:RO_Bot_Rc / 0.026 / 0.008 / 0.206 / 0.095 / 0.062 / 0.011 / 0.109
9:US_Mai_Pru / 0.484 / 0.493 / 0.689 / 0.242 / 0.429 / 0.510 / 0.234 / 0.537
10:SP_Uix_Au / 0.396 / 0.389 / 0.580 / 0.155 / 0.326 / 0.428 / 0.145 / 0.446 / 0.011
11:FR_Mon_Au / 0.025 / -0.009 / 0.210 / 0.094 / 0.056 / 0.015 / 0.117 / -0.001 / 0.559 / 0.461

Significant results (P<0.05) are in bold. Negative FST values indicate a greater genetic difference between individuals within a population than between individuals from different populations, i.e. more heterozygotes than expected.

Table S4 Models used in the phylogenetic reconstructions.

Gene / Position / Model / pi(A) / pi(C) / pi(G) / pi(T) / r(AC) / r(AG) / r(AT) / r(CG) / r(CT) / r(GT) / I / G
COI / 1st / GTR+G / 0.430 / 0.088 / 0.007 / 0.476 / 0.047 / 1.135 / 0.009 / 1.617 / 0.197 / 1.000 / - / 0.525
2nd / F81+I+G / 0.335 / 0.154 / 0.228 / 0.283 / 1.000 / 1.000 / 1.000 / 1.000 / 1.000 / 1.000 / 0.610 / 0.105
3rd / HKY+I / 0.179 / 0.230 / 0.175 / 0.416 / 1.000 / 2.436 / 1.000 / 1.000 / 2.436 / 1.000 / 0.873 / -
EF1α / 1st / K80 / 0.250 / 0.250 / 0.250 / 0.250 / 1.000 / 8.907 / 1.000 / 1.000 / 8.907 / 1.000 / - / -
2nd+3rd / JC+I / 0.250 / 0.250 / 0.250 / 0.250 / 1.000 / 1.000 / 1.000 / 1.000 / 1.000 / 1.000 / 0.925 / -
WING / 1st+2nd+3rd / SYM+I+G / 0.250 / 0.250 / 0.250 / 0.250 / 0.296 / 1.069 / 0.713 / 1.296 / 0.547 / 1.000 / 0.772 / 0.905

1st, 2nd and 3rd – codon positions; pi(A), pi(C), pi(G) and pi(T) – proportion of each nucleotide; r(AC), r(AG), r(AT), r(CG), r(CT) and r(GT) – rate of change from one nucleotide to another G - gamma-shape parameter; I - proportion of invariable sites.