Population genomics of two congeneric Palaearctic shorebirds reveals differential impacts of Quaternary climate oscillations across habitats types

Intracontinental biotic divisions across the vast Palaearctic region are not well-characterized. Past research has revealed patterns ranging from a lack of population structure to deep divergences along varied lines of separation. Here we compared biogeographic patterns of two Palaearctic shorebirds with different habitat preferences, Whimbrel ( Numenius phaeopus ) and Eurasian curlew ( N. arquata ). Using genome-wide markers from populations across the Palaearctic, we applied a multitude of population genomic and phylogenomic approaches to elucidate population structure. Most importantly, we tested for isolation by distance and visualized barriers and corridors to gene flow. We found shallow Palaearctic population structure in subpolar bog and tundra-breeding whimbrels, consistent with other species breeding at a similarly high latitude, indicating connectivity across the tundra belt, both presently and during southward shifts in periods of global cooling. In contrast, the temperate grassland-breeding Eurasian curlew emerged in three distinct clades corresponding to glacial refugia. Barriers to gene flow coincided with areas of topographic relief in the central Palaearctic for whimbrels and further east for Eurasian curlews. Our findings highlight the interplay of historic and ecological factors in influencing present-day population structure of Palaearctic biota.