The gray wolf (Canis lupus) is among the few large carnivores that survived the Late Pleistocene megafaunal extinctions. Thanks to their complex history of admixture and extensive geographic range, the number of gray wolf subspecies and their phylogenetic relationships remain poorly understood. Here, we perform whole-genome sequencing of a gray wolf collected from peninsular India that was phenotypically distinct from gray wolves outside India. Genomic analyses reveal that the Indian gray wolf is an evolutionarily distinct lineage that diverged from other extant gray wolf lineages ∼110 thousand years ago. Demographic analyses suggest that the Indian wolf population declined continuously decline since separating from other gray wolves and, today, has exceptionally low genetic diversity. We also find evidence for pervasive and mosaic gene flow between the Indian wolf and African canids including African wolf, Ethiopian wolf, and African wild dog despite their current geographical separation. Our results support the hypothesis that the Indian subcontinent was a Pleistocene refugium and center of diversification and further highlight the complex history of gene flow that characterized the evolution of gray wolves.

Genetic introgression not only provides material for adaptive evolution but also confounds our understanding of evolutionary history. This is particularly true for canids, a species complex in which genome sequencing and analysis has revealed a complex history of admixture and introgression. Here, we sequence 19 new whole genomes from high-altitude Tibetan and Himalayan wolves and dogs and combine these into a larger data set of 166 whole canid genomes. Using these data, we explore the evolutionary history and adaptation of these and other canid lineages. We find that Tibetan and Himalayan wolves are closely related to each other, and that ∼39% of their nuclear genome is derived from an as-yet-unrecognized wolf-like lineage that is deeply diverged from living Holarctic wolves and dogs. The EPAS1 haplotype, which is present at high frequencies in Tibetan dog breeds and wolves and confers an adaptive advantage to animals living at high altitudes, was probably derived from this ancient lineage. Our study underscores the complexity of canid evolution and demonstrates how admixture and introgression can shape the evolutionary trajectories of species.