Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx

Phylogenetic and phylogeographic assessments of species should ideally include individuals of all known populations. However, this is difficult for species with large distribution areas. Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed that the subspecies L. l. isabellinus harbours the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbours the second-most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species. Genetic variation present in southern Palearctic populations of Lynx lynx is clearly important for elucidating the evolutionary history of this species.