Shaped by uneven Pleistocene climate: mitochondrial phylogeographic pattern and population history of white wagtail Motacilla alba

We studied the phylogeography and population history of the white wagtail Motacilla alba , which has a vast breeding range, covering areas with diff erent Pleistocene climatic histories. Th e mitochondrial NADH dehydrogenase subunit II gene (ND2) and Control Region (CR) were analyzed for 273 individuals from 45 localities. Our data comprised all nine subspecies of white wagtail. Four primary clades were inferred (M, N, SW and SE), with indications of M. grandis being nested within M. alba . Th e oldest split was between two haplotypes from the endemic Moroccan M. a. subpersonata (clade M) and the others, at 0.63 – 0.96 Mya; other divergences were at 0.31 – 0.38 Mya. Th e entire diff erentiation falls within the part of the Pleistocene characterized by Milankovitch cycles of large amplitudes and durations. Clade N was distributed across the northern Palearctic; clade SW in southwestern Asia plus the British Isles and was predicted by Ecological niche models (ENMs) to occur also in central and south Europe; and clade SE was distributed in central and east Asia. Th e deep divergence within M. a. subpersonata may refl ect retention of ancestral haplotypes. Regional diff erences in historical climates have had diff erent impacts on diff erent populations: clade N expanded after the last glacial maximum (LGM), whereas milder Pleistocene climate of east Asia allowed clade SE a longer expansion time (since MIS 5); clade SW expanded over a similarly long time as clade SE, which is untypical for European species. ENMs supported these conclusions in that the northern part of the Eurasian continent was unsuitable during the LGM, whereas southern parts remained suitable. Th e recent divergences and poor structure in the mitochondrial tree contrasts strongly with the pronounced, well defi ned phenotypical diff erentiation, indicating extremely fast plumage divergence.