Results of generalized linear mixed models.

Tests were conducted at the intraspecific level using phylogroups (upper table) and at the interspecific level using nearest neighbour species (lower table) mtDNA sequence divergence data. For the model testing the presence of a breakpoint at glacier margin, we conducted a piece-wise regression involving a break at 46o of latitude which corresponds to the maximal extent of Pleistocene glaciations events.