Explore the vast range of titles available.

We studied the relationship between the variability and contemporary distribution of pelage phenotypes in one of most widely distributed felid species and an array of environmental and demographic conditions. We collected 672 photographic georeferenced records of the Eurasian lynx throughout Eurasia. We assigned each lynx coat to one of five phenotypes. Then we fitted the coat patterns to different environmental and anthropogenic variables, as well as the effective geographic distances from inferred glacial refugia. A majority of lynx were either of the large spotted (41.5%) or unspotted (uniform, 36.2%) phenotype. The remaining patterns (rosettes, small spots and pseudo-rosettes) were represented in 11.0%, 7.4%, and 3.9% of samples, respectively. Although various environmental variables greatly affected lynx distribution and habitat suitability, it was the effect of least-cost distances from locations of the inferred refugia during the Last Glacial Maximum that explained the distribution of lynx coat patterns the best. Whereas the occurrence of lynx phenotypes with large spots was explained by the proximity to refugia located in the Caucasus/Middle East, the uniform phenotype was associated with refugia in the Far East and Central Asia. Despite the widely accepted hypothesis of adaptive functionality of coat patterns in mammals and exceptionally high phenotypic polymorphism in Eurasian lynx, we did not find well-defined signs of habitat matching in the coat pattern of this species. Instead, we showed how the global patterns of morphological variability in this large mammal and its environmental adaptations may have been shaped by past climatic change.