Intravital Dynamic and Correlative Imaging of Mouse Livers Reveals Diffusion-Dominated Canalicular and Flow-Augmented Ductular Bile Flux.

BACKGROUND AND AIMS: Small-molecule flux in tissue microdomains is essential for organ function, but knowledge of this process is scant due to the lack of suitable methods. We developed two independent techniques that allow the quantification of advection (flow) and diffusion in individual bile canaliculi and in interlobular bile ducts of intact livers in living mice, namely fluorescence loss after photoactivation and intravital arbitrary region image correlation spectroscopy. APPROACH AND RESULTS: The results challenge the prevailing "mechano-osmotic" theory of canalicular bile flow. After active transport across hepatocyte membranes, bile acids are transported in the canaliculi primarily by diffusion. Only in the interlobular ducts is diffusion augmented by regulatable advection. Photoactivation of fluorescein bis-(5-carboxymethoxy-2-nitrobenzyl)-ether in entire lobules demonstrated the establishment of diffusive gradients in the bile canalicular network and the sink function of interlobular ducts. In contrast to the bile canalicular network, vectorial transport was detected and quantified in the mesh of interlobular bile ducts. CONCLUSIONS: The liver consists of a diffusion-dominated canalicular domain, where hepatocytes secrete small molecules and generate a concentration gradient and a flow-augmented ductular domain, where regulated water influx creates unidirectional advection that augments the diffusive flux.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Tissue-specific epigenetic inheritance after paternal heat exposure in male wild guinea pigs.

External temperature change has been shown to modify epigenetic patterns, such as DNA methylation, which regulates gene expression. DNA methylation is heritable, and as such provides a mechanism to convey environmental information to subsequent generations. Studies on epigenetic response to temperature increase are still scarce in wild mammals, even more so studies that compare tissue-specific epigenetic responses. Here, we aim to address differential epigenetic responses on a gene and gene pathway level in two organs, liver and testis. We chose these organs, because the liver is the main metabolic and thermoregulation organ, and epigenetic modifications in testis are potentially transmitted to the F2 generation. We focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to an ambient temperature increase of 10 °C, and investigated differential methylated regions of sons sired before and after the paternal exposure using Reduced Representation Bisulfite Sequencing. We detected both a highly tissue-specific epigenetic response, reflected in genes involved in organ-specific metabolic pathways, and a more general regulation of single genes epigenetically modified in both organs. We conclude that genomes are context-specifically differentially epigenetically regulated in response to temperature increase. These findings emphasize the epigenetic relevance in cell differentiation, which is essential for the specific function(s) of complex organs, and is represented in a diverse molecular regulation of genes and gene pathways. The results also emphasize the paternal contribution to adaptive processes.

Phylogeography of red muntjacs reveals three distinct mitochondrial lineages.

BACKGROUND: The members of the genus Muntiacus are of particular interest to evolutionary biologists due to their extreme chromosomal rearrangements and the ongoing discussions about the number of living species. Red muntjacs have the largest distribution of all muntjacs and were formerly considered as one species. Karyotype differences led to the provisional split between the Southern Red Muntjac (Muntiacus muntjak) and the Northern Red Muntjac (M. vaginalis), but uncertainties remain as, so far, no phylogenetic study has been conducted. Here, we analysed whole mitochondrial genomes of 59 archival and 16 contemporaneous samples to resolve uncertainties about their taxonomy and used red muntjacs as model for understanding the evolutionary history of other species in Southeast Asia. RESULTS: We found three distinct matrilineal groups of red muntjacs: Sri Lankan red muntjacs (including the Western Ghats) diverged first from other muntjacs about 1.5 Mya; later northern red muntjacs (including North India and Indochina) and southern red muntjacs (Sundaland) split around 1.12 Mya. The diversification of red muntjacs into these three main lineages was likely promoted by two Pleistocene barriers: one through the Indian subcontinent and one separating the Indochinese and Sundaic red muntjacs. Interestingly, we found a high level of gene flow within the populations of northern and southern red muntjacs, indicating gene flow between populations in Indochina and dispersal of red muntjacs over the exposed Sunda Shelf during the Last Glacial Maximum. CONCLUSIONS: Our results provide new insights into the evolution of species in South and Southeast Asia as we found clear genetic differentiation in a widespread and generalist species, corresponding to two known biogeographical barriers: The Isthmus of Kra and the central Indian dry zone. In addition, our molecular data support either the delineation of three monotypic species or three subspecies, but more importantly these data highlight the conservation importance of the Sri Lankan/South Indian red muntjac.

Resistance to RHD virus in wild Australian rabbits: Comparison of susceptible and resistant individuals using a genomewide approach.

Deciphering the genes involved in disease resistance is essential if we are to understand host-pathogen coevolutionary processes. The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus). During the first outbreaks of the disease, RHDV caused mortality rates of up to 97%. Recently, however, increased genetic resistance to RHDV has been reported. Here, we have aimed to identify genomic differences between rabbits that survived a natural infection with RHDV and those that died in the field using a genomewide next-generation sequencing (NGS) approach. We detected 72 SNPs corresponding to 133 genes associated with survival of a RHD infection. Most of the identified genes have known functions in virus infections and replication, immune responses or apoptosis, or have previously been found to be regulated during RHD. Some of the genes identified in experimental studies, however, did not seem to play a role under natural selection regimes, highlighting the importance of field studies to complement the genomic background of wildlife diseases. Our study provides a set of candidate markers as a tool for the future scanning of wild rabbits for their resistance to RHDV. This is important both for wild rabbit populations in southern Europe where RHD is regarded as a serious problem decimating the prey of endangered predator species and for assessing the success of currently planned RHDV variant biocontrol releases in Australia.

Genetic Structure and Phylogeography of the Leopard Cat (Prionailurus bengalensis) Inferred from Mitochondrial Genomes.

The Leopard cat Prionailurus bengalensis is a habitat generalist that is widely distributed across Southeast Asia. Based on morphological traits, this species has been subdivided into 12 subspecies. Thus far, there have been few molecular studies investigating intraspecific variation, and those had been limited in geographic scope. For this reason, we aimed to study the genetic structure and evolutionary history of this species across its very large distribution range in Asia. We employed both PCR-based (short mtDNA fragments, 94 samples) and high throughput sequencing based methods (whole mitochondrial genomes, 52 samples) on archival, noninvasively collected and fresh samples to investigate the distribution of intraspecific genetic variation. Our comprehensive sampling coupled with the improved resolution of a mitochondrial genome analyses provided strong support for a deep split between Mainland and Sundaic Leopard cats. Although we identified multiple haplogroups within the species' distribution, we found no matrilineal evidence for the distinction of 12 subspecies. In the context of Leopard cat biogeography, we cautiously recommend a revision of the Prionailurus bengalensis subspecific taxonomy: namely, a reduction to 4 subspecies (2 mainland and 2 Sundaic forms).

Phylogeography of the Small Indian Civet and Origin of Introductions to Western Indian Ocean Islands.

The biogeographic dynamics affecting the Indian subcontinent, East and Southeast Asia during the Plio-Pleistocene has generated complex biodiversity patterns. We assessed the molecular biogeography of the small Indian civet (Viverricula indica) through mitogenome and cytochrome b + control region sequencing of 89 historical and modern samples to (1) establish a time-calibrated phylogeography across the species' native range and (2) test introduction scenarios to western Indian Ocean islands. Bayesian phylogenetic analyses identified 3 geographic lineages (East Asia, sister-group to Southeast Asia and the Indian subcontinent + northern Indochina) diverging 3.2-2.3 million years ago (Mya), with no clear signature of past demographic expansion. Within Southeast Asia, Balinese populations separated from the rest 2.6-1.3 Mya. Western Indian Ocean populations were assigned to the Indian subcontinent + northern Indochina lineage and had the lowest mitochondrial diversity. Approximate Bayesian computation did not distinguish between single versus multiple introduction scenarios. The early diversification of the small Indian civet was likely shaped by humid periods in the Late Pliocene-Early Pleistocene that created evergreen rainforest barriers, generating areas of intra-specific endemism in the Indian subcontinent, East, and Southeast Asia. Later, Pleistocene dispersals through drier conditions in South and Southeast Asia were likely, giving rise to the species' current natural distribution. Our molecular data supported the delineation of only 4 subspecies in V. indica, including an endemic Balinese lineage. Our study also highlighted the influence of prefirst millennium AD introductions to western Indian Ocean islands, with Indian and/or Arab traders probably introducing the species for its civet oil.

Paternal intergenerational epigenetic response to heat exposure in male Wild guinea pigs.

Epigenetic modifications, of which DNA methylation is the best studied one, can convey environmental information through generations via parental germ lines. Past studies have focused on the maternal transmission of epigenetic information to the offspring of isogenic mice and rats in response to external changes, whereas heterogeneous wild mammals as well as paternal epigenetic effects have been widely neglected. In most wild mammal species, males are the dispersing sex and have to cope with differing habitats and thermal changes. As temperature is a major environmental factor we investigated if genetically heterogeneous Wild guinea pig (Cavia aperea) males can adapt epigenetically to an increase in temperature and if that response will be transmitted to the next generation(s). Five adult male guinea pigs (F0) were exposed to an increased ambient temperature for 2 months, i.e. the duration of spermatogenesis. We studied the liver (as the main thermoregulatory organ) of F0 fathers and F1 sons, and testes of F1 sons for paternal transmission of epigenetic modifications across generation(s). Reduced representation bisulphite sequencing revealed shared differentially methylated regions in annotated areas between F0 livers before and after heat treatment, and their sons' livers and testes, which indicated a general response with ecological relevance. Thus, paternal exposure to a temporally limited increased ambient temperature led to an 'immediate' and 'heritable' epigenetic response that may even be transmitted to the F2 generation. In the context of globally rising temperatures epigenetic mechanisms may become increasingly relevant for the survival of species.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population.

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

A haplotype-resolved reference genome of a long-distance migratory bat, Pipistrellus nathusii (Keyserling & Blasius, 1839).

We present a complete, chromosome-scale reference genome for the long-distance migratory bat Pipistrellus nathusii. The genome encompasses both haplotypic sets of autosomes and the separation of both sex chromosomes by utilizing highly accurate long-reads and preserving long-range phasing information through the use of three-dimensional chromatin conformation capture sequencing (Hi-C). This genome, accompanied by a comprehensive protein-coding sequence annotation, provides a valuable genomic resource for future investigations into the genomic bases of long-distance migratory flight in bats as well as uncovering the genetic architecture, population structure and evolutionary history of Pipistrellus nathusii. The reference-quality genome presented here gives a fundamental resource to further our understanding of bat genetics and evolution, adding to the growing number of high-quality genetic resources in this field. Here, we demonstrate its use in the phylogenetic reconstruction of the order Chiroptera, and in particular, we present the resources to allow detailed investigations into the genetic drivers and adaptations related to long-distance migration.

Epigenetic signatures of social status in wild female spotted hyenas (Crocuta crocuta).

In mammalian societies, dominance hierarchies translate into inequalities in health, reproductive performance and survival. DNA methylation is thought to mediate the effects of social status on gene expression and phenotypic outcomes, yet a study of social status-specific DNA methylation profiles in different age classes in a wild social mammal is missing. We tested for social status signatures in DNA methylation profiles in wild female spotted hyenas (Crocuta crocuta), cubs and adults, using non-invasively collected gut epithelium samples. In spotted hyena clans, female social status influences access to resources, foraging behavior, health, reproductive performance and survival. We identified 149 differentially methylated regions between 42 high- and low-ranking female spotted hyenas (cubs and adults). Differentially methylated genes were associated with energy conversion, immune function, glutamate receptor signalling and ion transport. Our results provide evidence that socio-environmental inequalities are reflected at the molecular level in cubs and adults in a wild social mammal.

Coronavirus genotype diversity and prevalence of infection in wild carnivores in the Serengeti National Park, Tanzania.

Knowledge of coronaviruses in wild carnivores is limited. This report describes coronavirus genetic diversity, species specificity and infection prevalence in three wild African carnivores. Coronavirus RNA was recovered from fresh feces from spotted hyena and silver-backed jackal, but not bat-eared fox. Analysis of sequences of membrane (M) and spike (S) gene fragments revealed strains in the genus Alphacoronavirus, including three distinct strains in hyenas and one distinct strain in a jackal. Coronavirus RNA prevalence was higher in feces from younger (17 %) than older (3 %) hyenas, highlighting the importance of young animals for coronavirus transmission in wild carnivores.

First mitogenome phylogeny of the sun bear Helarctos malayanus reveals a deep split between Indochinese and Sundaic lineages.

The sun bear Helarctos malayanus is one of the most endangered ursids, and to date classification of sun bear populations has been based almost exclusively on geographic distribution and morphology. The very few molecular studies focussing on this species were limited in geographic scope. Using archival and non-invasively collected sample material, we have added a substantial number of complete or near-complete mitochondrial genome sequences from sun bears of several range countries of the species' distribution. We here report 32 new mitogenome sequences representing sun bears from Cambodia, Thailand, Peninsular Malaysia, Sumatra, and Borneo. Reconstruction of phylogenetic relationships revealed two matrilines that diverged ~295 thousand years ago: one restricted to portions of mainland Indochina (China, Cambodia, Thailand; "Mainland clade"), and one comprising bears from Borneo, Sumatra, Peninsular Malaysia but also Thailand ("Sunda clade"). Generally recent coalescence times in the mitochondrial phylogeny suggest that recent or historical demographic processes have resulted in a loss of mtDNA variation. Additionally, analysis of our data in conjunction with shorter mtDNA sequences revealed that the Bornean sun bear, classified as a distinct subspecies (H. m. euryspilus), does not harbor a distinctive matriline. Further molecular studies of H. malayanus are needed, which should ideally include data from nuclear loci.

First Steps towards the Development of Epigenetic Biomarkers in Female Cheetahs (Acinonyx jubatus).

Free-ranging cheetahs (Acinonyx jubatus) are generally healthy, whereas cheetahs under human care, such as those in zoological gardens, suffer from ill-defined infectious and degenerative pathologies. These differences are only partially explained by husbandry management programs because both groups share low genetic diversity. However, mounting evidence suggests that physiological differences between populations in different environments can be tracked down to differences in epigenetic signatures. Here, we identified differentially methylated regions (DMRs) between free-ranging cheetahs and conspecifics in zoological gardens and prospect putative links to pathways relevant to immunity, energy balance and homeostasis. Comparing epigenomic DNA methylation profiles obtained from peripheral blood mononuclear cells (PBMCs) from eight free-ranging female cheetahs from Namibia and seven female cheetahs living in zoological gardens within Europe, we identified DMRs of which 22 were hypermethylated and 23 hypomethylated. Hypermethylated regions in cheetahs under human care were located in the promoter region of a gene involved in host-pathogen interactions (KLC1) and in an intron of a transcription factor relevant for the development of pancreatic ß-cells, liver, and kidney (GLIS3). The most canonical mechanism of DNA methylation in promoter regions is assumed to repress gene transcription. Taken together, this could indicate that hypermethylation at the promoter region of KLC1 is involved in the reduced immunity in cheetahs under human care. This approach can be generalized to characterize DNA methylation profiles in larger cheetah populations under human care with a more granular longitudinal data collection, which, in the future, could be used to monitor the early onset of pathologies, and ultimately translate into the development of biomarkers with prophylactic and/or therapeutic potential.

Assessing Coagulation Parameters in Healthy Asian Elephants (Elephas maximus) from European and Thai Populations.

The Asian elephant population is continuously declining due to several extrinsic reasons in their range countries, but also due to diseases in captive populations worldwide. One of these diseases, the elephant endotheliotropic herpesvirus (EEHV) hemorrhagic disease, is very impactful because it particularly affects Asian elephant calves. It is commonly fatal and presents as an acute and generalized hemorrhagic syndrome. Therefore, having reference values of coagulation parameters, and obtaining such values for diseased animals in a very short time, is of great importance. We analyzed prothrombin time (PT), activated partial thromboplastin time (aPTT), and fibrinogen concentrations using a portable and fast point-of-care analyzer (VetScan Pro) in 127 Asian elephants from Thai camps and European captive herds. We found significantly different PT and aPTT coagulation times between elephants from the two regions, as well as clear differences in fibrinogen concentration. Nevertheless, these alterations were not expected to have biological or clinical implications. We have also sequenced the coagulation factor VII gene of 141 animals to assess the presence of a previously reported hereditary coagulation disorder in Asian elephants and to investigate the presence of other mutations. We did not find the previously reported mutation in our study population. Instead, we discovered the presence of several new single nucleotide polymorphisms, two of them being considered as deleterious by effect prediction software.

Genetic population structure defines wild boar as an urban exploiter species in Barcelona, Spain.

Urban wildlife ecology is gaining relevance as metropolitan areas grow throughout the world, reducing natural habitats and creating new ecological niches. However, knowledge is still scarce about the colonisation processes of such urban niches, the establishment of new communities, populations and/or species, and the related changes in behaviour and life histories of urban wildlife. Wild boar (Sus scrofa) has successfully colonised urban niches throughout Europe. The aim of this study is to unveil the processes driving the establishment and maintenance of an urban wild boar population by analysing its genetic structure. A set of 19 microsatellite loci was used to test whether urban wild boars in Barcelona, Spain, are an isolated population or if gene flow prevents genetic differentiation between rural and urban wild boars. This knowledge will contribute to the understanding of the effects of synurbisation and the associated management measures on the genetic change of large mammals in urban ecosystems. Despite the unidirectional gene flow from rural to urban areas, the urban wild boars in Barcelona form an island population genotypically differentiated from the surrounding rural ones. The comparison with previous genetic studies of urban wild boar populations suggests that forest patches act as suitable islands for wild boar genetic differentiation. Previous results and the genetic structure of the urban wild boar population in Barcelona classify wild boar as an urban exploiter species. These wild boar peri-urban island populations are responsible for conflict with humans and thus should be managed by reducing the attractiveness of urban areas. The management of peri-urban wild boar populations should aim at reducing migration into urban areas and preventing phenotypic changes (either genetic or plastic) causing habituation of wild boars to humans and urban environments.

Mating system, feeding type and ex situ conservation effort determine life expectancy in captive ruminants.

Zoo animal husbandry aims at constantly improving husbandry, reproductive success and ultimately animal welfare. Nevertheless, analyses to determine factors influencing husbandry of different species are rare. The relative life expectancy (rLE; life expectancy (LE) as proportion of longevity) describes husbandry success of captive populations. Correlating rLE with biological characteristics of different species, reasons for variation in rLE can be detected. We analysed data of 166 901 animals representing 78 ruminant species kept in 850 facilities. The rLE of females correlated with the percentage of grass in a species' natural diet, suggesting that needs of species adapted to grass can be more easily accommodated than the needs of those adapted to browse. Males of monogamous species demonstrate higher rLE than polygamous males, which matches observed differences of sexual bias in LE in free-living populations and thus supports the ecological theory that the mating system influences LE. The third interesting finding was that rLE was higher in species managed by international studbooks when compared with species not managed in this way. Our method facilitates the identification of biological characteristics of species that are relevant for their husbandry success, and they also support ecological theory. Translating these findings into feeding recommendations, our approach can help to improve animal husbandry.

Geographical variation in and evolutionary history of the Sunda clouded leopard (Neofelis diardi) (Mammalia: Carnivora: Felidae) with the description of a new subspecies from Borneo.

Recent morphological and molecular studies led to the recognition of two extant species of clouded leopards; Neofelis nebulosa from mainland southeast Asia and Neofelis diardi from the Sunda Islands of Borneo and Sumatra, including the Batu Islands. In addition to these new species-level distinctions, preliminary molecular data suggested a genetic substructure that separates Bornean and Sumatran clouded leopards, indicating the possibility of two subspecies of N. diardi. This suggestion was based on an analysis of only three Sumatran and seven Bornean individuals. Accordingly, in this study we re-evaluated this proposed subspecies differentiation using additional molecular (mainly historical) samples of eight Bornean and 13 Sumatran clouded leopards; a craniometric analysis of 28 specimens; and examination of pelage morphology of 20 museum specimens and of photographs of 12 wild camera-trapped animals. Molecular (mtDNA and microsatellite loci), craniomandibular and dental analyses strongly support the differentiation of Bornean and Sumatran clouded leopards, but pelage characteristics fail to separate them completely, most probably owing to small sample sizes, but it may also reflect habitat similarities between the two islands and their recent divergence. However, some provisional discriminating pelage characters are presented that need further testing. According to our estimates both populations diverged from each other during the Middle to Late Pleistocene (between 400 and 120 kyr). We present a discussion on the evolutionary history of Neofelis diardi sspp. on the Sunda Shelf, a revised taxonomy for the Sunda clouded leopard, N. diardi, and formally describe the Bornean subspecies, Neofelis diardi borneensis, including the designation of a holotype (BM.3.4.9.2 from Baram, Sarawak) in accordance with the rules of the International Code of Zoological Nomenclature.

Determinants of defence strategies of a hibernating European bat species towards the fungal pathogen Pseudogymnoascus destructans.

Pseudogymnoascus destructans (Pd), the causative agent of white-nose syndrome in North America, has decimated bat populations within a decade. The fungus impacts bats during hibernation when physiological functions, including immune responses, are down-regulated. Studies have shown that Pd is native to Europe, where it is not associated with mass mortalities. Moreover, genomic and proteomic studies indicated that European bats may have evolved an effective immune defence, which is lacking in North American bats. However, it is still unclear which defence strategy enables European bats to cope with the pathogen. Here, we analyzed selected physiological and immunological parameters in torpid, Pd infected European greater mouse-eared bats (Myotis myotis) showing three different levels of infection (asymptomatic, mild and severe symptoms). From a subset of the studied bats we tracked skin temperatures during one month of hibernation. Contrasting North American bats, arousal patterns remained unaffected by Pd infections in M. myotis. In general, heavier M. myotis aroused more often from hibernation and showed less severe disease symptoms than lean individuals; most likely because heavy bats were capable of reducing the Pd load more effectively than lean individuals. In the blood of severely infected bats, we found higher gene expression levels of an inflammatory cytokine (IL-1ß), but lower levels of an acute phase protein (haptoglobin), reactive oxygen metabolites (ROMs) and plasma non-enzymatic antioxidant capacity (OXY) compared to conspecifics with lower levels of infection. We conclude that M. myotis, and possibly also other European bat species, tolerate Pd infections during torpor by using selected acute phase response parameters at baseline levels, yet without arousing from torpor and without synthesizing additional immune molecules.

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

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 the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies' current range. L. l. dinniki harbors 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.