- Research article
- Open Access
Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations
BMC Evolutionary Biologyvolume 18, Article number: 156 (2018)
Resolving the historical biogeography of the leopard (Panthera pardus) is a complex issue, because patterns inferred from fossils and from molecular data lack congruence. Fossil evidence supports an African origin, and suggests that leopards were already present in Eurasia during the Early Pleistocene. Analysis of DNA sequences however, suggests a more recent, Middle Pleistocene shared ancestry of Asian and African leopards. These contrasting patterns led researchers to propose a two-stage hypothesis of leopard dispersal out of Africa: an initial Early Pleistocene colonisation of Asia and a subsequent replacement by a second colonisation wave during the Middle Pleistocene. The status of Late Pleistocene European leopards within this scenario is unclear: were these populations remnants of the first dispersal, or do the last surviving European leopards share more recent ancestry with their African counterparts?
In this study, we generate and analyse mitogenome sequences from historical samples that span the entire modern leopard distribution, as well as from Late Pleistocene remains. We find a deep bifurcation between African and Eurasian mitochondrial lineages (~ 710 Ka), with the European ancient samples as sister to all Asian lineages (~ 483 Ka). The modern and historical mainland Asian lineages share a relatively recent common ancestor (~ 122 Ka), and we find one Javan sample nested within these.
The phylogenetic placement of the ancient European leopard as sister group to Asian leopards suggests that these populations originate from the same out-of-Africa dispersal which founded the Asian lineages. The coalescence time found for the mitochondrial lineages aligns well with the earliest undisputed fossils in Eurasia, and thus encourages a re-evaluation of the identification of the much older putative leopard fossils from the region. The relatively recent ancestry of all mainland Asian leopard lineages suggests that these populations underwent a severe population bottleneck during the Pleistocene. Finally, although only based on a single sample, the unexpected phylogenetic placement of the Javan leopard could be interpreted as evidence for exchange of mitochondrial lineages between Java and mainland Asia, calling for further investigation into the evolutionary history of this subspecies.
Achieving a comprehensive understanding of a species’ history is important for both evolutionary research and for conservation management. However, this may be impossible using data derived solely from living individuals – particularly for endangered species whose current genetic diversity is depauperate. A potential solution is to study DNA sequences obtained from historical or ancient samples , allowing extinct populations and lineages to be investigated and compared with those that exist today. In studies that use such samples, mitochondrial DNA (mtDNA) continues to be an important marker as the much higher copy number per cell compared to the nuclear genome generally results in a higher success rate for sequence recovery – particularly for species with low fossil abundance and/or poor biomolecular preservation (e.g. [2,3,4,5]). The leopard (Panthera pardus Linnaeus, 1758) is an example of a species that is currently distributed across only a fraction of its historical and ancient range (e.g. [6,7,8]). It is one of only a few large-bodied carnivore species that naturally occurs in a wide variety of habitats; from the Himalayan highlands to the Ethiopian desert, and from the Congo rainforests to the Amur taiga . Human persecution and hunting (e.g. [10,11,12]), habitat destruction (e.g. [13,14,15]) and reduced prey availability [16, 17] have severely impacted the distribution of this elusive predator, and leopards are now extinct in large parts of their historic Asian and African distribution (Fig. 1) [7, 18].
The fossil record and genetic data from modern and historical samples are generally interpreted as indicative of an African origin of the leopard. The oldest leopard fossils have been recovered in Eastern Africa, and genetic diversity (estimated using both mtDNA and nuclear microsatellites) in living African populations is higher than in Asian populations, which has been interpreted as evidence for an African origin [19, 20]. The oldest potential evidence for leopards outside of Africa is found in the Early Pleistocene fossil record from South Asia (Pakistan), suggesting an initial out-of-Africa expansion into Eurasia around 2.0 Ma (Mega annus; million years ago; ). Mitochondrial and microsatellite data, however, suggest that all current Asian populations share a much more recent common ancestor (approximately 622 Ka [Kilo annus; thousand years ago]; [19, 22]); more than a million years younger than the oldest fossils found in Asia. The apparent incongruence between the fossil record in Eurasia (which suggests occupation since the Early Pleistocene) and the relatively recent coalescence time of African and Asian modern leopard mitochondrial lineages has been interpreted as indication for two independent out-of-Africa dispersal events, the latter of which founded all modern Asian leopard lineages .
Previous studies that have examined intraspecific variation in leopards have been based on short mitochondrial sequences and microsatellites from recently collected samples [19, 23,24,25,26,27,28,29]. In this study, we present near-complete leopard mitogenomes from seven Late Pleistocene specimens (up to 45,000 years old) and from 15 historical samples collected up to 150 years ago, from a range of geographical locations encompassing its entire distribution (Fig. 1). We investigate this mitogenomic data in the context of the proposed evolutionary history and past population dynamics of the leopard, and the role of Pleistocene populations within this scenario. Furthermore, by combining mitochondrial data from previous studies [19, 22, 24, 27] with ancient and historical mtDNA, we present new insights into ancient and recent population dynamics of one of the most widespread extant felid species.
Using hybridisation capture and high-throughput sequencing, we retrieved near-complete mitogenomes from 22 leopard samples: one Late Pleistocene Caucasian leopard (~ 35 Ka), six Late Pleistocene Central European leopards (~ 45 Ka), and 15 historical leopards from across the entire historical distribution (Fig. 1; Table 1). The mitogenome sequences represent 17 distinct haplotypes, which were analysed in conjunction with three previously published modern leopard mitogenomes. Using a two-step approach (see Methods section for details), we estimated the coalescence times of this dataset in BEAST. Our results support the deep bifurcation between Asian and African leopards, that has previously been proposed based on short mtDNA sequences (Fig. 2; Additional file 1: Figure S1; [19, 22]). The fossil-calibrated Bayesian analysis of the basal divergence time for all leopard mitochondrial lineages was 710 Ka (95% credibility interval [CI]: 457–956 Ka), which is also consistent with the previous estimates (932 Ka; , 471–825 Ka; ). The Pleistocene European sequences form a well-supported clade consisting of three distinct haplotypes, which is sister to a clade containing all mitogenome sequences from Asian leopards. The two clades are estimated to have diverged approximately 483 Ka (95% CI: 305–677 Ka). The ancient Caucasian sequence is sister to all modern mainland Asian sequences, with high support (Bayesian posterior probability of 1.0, Fig. 2; Additional file 1: Figure S1). The estimated coalescence time of mainland Asian mitogenomes, including the ancient Caucasian individual, is 244 Ka (95% CI: 148–352 Ka).
The South-East Asian samples included one specimen collected in Sumatra (PP35), where currently no leopards live and thus is likely to represent a traded specimen imported from elsewhere. A second specimen in this clade had no associated geographical provenance (‘East Indies’; PP3). These individuals may represent the Javan leopard lineage, as their mitogenome haplotypes are sister to all Asian leopard sequences (Fig. 2), and they exhibit ND5 sequences identical to previously published sequences from Javan leopards ([19, 22]; Additional file 1: Figure S2). The coalescence time of this lineage and the mainland leopards is estimated at 375 Ka (95% CI: 230–524 Ka). In contrast to previous studies, our Javan sample (PP32) was placed as a sister lineage to a specimen from Thailand (PP15), nested within all mainland Asian leopards rather than sister to these, with an estimated coalescence time of 64 Ka (95% CI: 32–100 Ka; Fig. 2).
Using modern, historical and ancient mitogenomic data from leopard samples from across their current and former geographic range, we provide novel insights into the historical biogeography of the leopard. The resulting data clarify the relationship between current leopard populations and Late Pleistocene fossils, and provide additional evidence for the interpretation of the earliest putative leopard fossils.
Origin of the leopard
Based on the fossil record, the origin of the leopard has been placed in Eastern Africa. There are fossils that may belong to P. pardus dating to about 3.4–3.8 Ma from Laetoli (Olduvai; ), although some authors suggest that these may be assigned to a different Panthera species , or some other large-bodied felid . Although the fossil record between 2.0 and 3.8 Ma is sparse, unequivocally identified leopard fossils do confirm their presence in Eastern Africa around 2 Ma , which strongly suggests an African origin of the species (Fig. 3).
We estimate the coalescence time between African and Eurasian leopard mitochondrial lineages to be ~ 710 Ka, similar to the age found by previous molecular studies [19, 22]. Although the bifurcation between African and Eurasian leopards provides no confirmation for either an African or Eurasian origin for the leopard, the combined fossil and genetic evidence together support Africa as the most likely place of origin. The relatively recent coalescence time between African and Eurasian mitochondrial lineages is also consistent with previous suggestions that the oldest African fossils of 3.8 Ma may not in fact represent P. pardus [31, 32]. Within Africa, we found considerable genetic divergence between haplotypes occurring in East Africa. In particular, the haplotype sampled from the Burundi individual forms a divergent sister lineage to all other sampled African haplotypes (95% CI: 368–814 Ka), including those from other regions of East Africa. Thus, the mitochondrial divergence occurring within this region of Africa is equivalent to that of the entire African continent, indicating this region as the potential point of origin of the leopard, which has also been suggested based on short mtDNA sequences . Given that, based on the fossil record, leopards were present in East Africa at least around 2 Ma ago , this hypothesis requires that mitochondrial lineages established prior to ~ 600 Ka have either been lost from modern populations through genetic drift or have not yet been sampled. These results do suggest, however, that, as well as for the genus Homo and a number of other species , East Africa may be the point of origin of modern leopards, although additional sampling of individuals and of nuclear markers is desirable to more robustly test this hypothesis.
The first evidence for leopards outside of Africa is ambiguous. In Asia, the earliest occurrence of leopard fossil remains is in the Early Pleistocene of South Asia (Pakistan), suggesting an initial out-of-Africa expansion into Asia around 2.0 Ma . Whether or not these fossils should be assigned to P. pardus is subject to discussion, however [31, 32]. These ancient specimens from Pakistan may be better assigned to other medium-sized felids (e.g. Eurasian puma, jaguar or snow leopard), or indicate an earlier dispersal of a leopard-like Panthera taxon into Asia (Fig. 3). Unequivocal leopard fossils in Asia are much younger; 0.6–0.8 Ma [35, 36]. Furthermore, in Europe, the oldest findings of this species are dated to the early Middle Pleistocene (nearly 0.6 Ma: [37,38,39].
Inferring the timing of colonisation events from the mitochondrial gene tree is challenging, as the population divergence and lineage coalescence will be asynchronous. Lineage coalescence of founding and source populations represents the maximum limit for dispersal time, but could be an overestimate caused by deep coalescence within the source population. Assuming monophyly of founding population lineages, the radiation of the founding population represents the minimum limit on the dispersal time, but this will tend to be an underestimation due to the mutation-lag of the formation of new lineages, incomplete sampling of extant lineages, or by lineage extinction during population bottlenecks. Thus, the true colonisation event lies somewhere along the branch connecting these upper and lower limits. Following this reasoning, the mitogenome data suggest the colonisation of Eurasia (or at least the end of maternal gene flow) some time between 710 and 483 Ka ago. This range overlaps with the age of the younger, unequivocal leopard fossils from Asia but not with the older Early Pleistocene Asian fossils whose precise taxonomic assignment has been debated. Achieving a unified biogeographic hypothesis for the colonisation of Asia by leopards thus hinges on the classification of these older fossils. Under the assumption that the first leopard occurrence in Eurasia is represented by the unequivocally identified specimens dating to around 0.6–0.8 Ma, there is no contradiction between the fossil record and the molecular evidence. If, however, older Early Pleistocene specimens are to be classified as P. pardus, then multiple out-of-Africa dispersals and the replacement of all mitochondrial lineages present in Asia from earlier dispersal events is required to explain the observed patterns. Considering the ambiguity of the Early Pleistocene fossil record in Asia, we consider a single out-of-Africa dispersal for P. pardus more parsimonious given the currently available combined palaeontological and genetic evidence. Clarifying whether the African and Eurasian leopard populations have indeed been separated since this time, or if gene flow continued to some extent, and if the earlier Eurasian leopard-like populations played a role in the evolutionary history of the more recent populations, will require nuclear genomic data with a similar geographical coverage.
Ancient Eurasian leopards
The origin and classification of the European Pleistocene leopard is equivocal based on fossil remains from the region (e.g. [6, 40,41,42]). The fossil record suggests that leopards first dispersed into Europe during the Middle Pleistocene, as evidenced by fossils from Western Europe estimated to be around 600 Ka (e.g. [6, 30, 35, 37, 38]). Several different forms have been defined in European fossil assemblages based on morphology, which has been interpreted as evidence for several immigration waves from Africa or Central Asia . The Late Pleistocene form, representing the age of the specimens we retrieved mitogenomes from, has been referred to as the subspecies P. p. spelaea [6, 42]. These populations likely went extinct before the Last Glacial Maximum (LGM), around 25,000 years ago [6, 8, 40, 41, 43]. The disappearance of the European leopard populations coincides with the extinction of many other Pleistocene carnivores, such as the cave bear (Ursus spelaeus), cave hyena (Crocuta crocuta spelaea) and Homotherium [3, 44,45,46].
We recover three distinct haplotypes from the six leopard samples, and as there are no overlapping skeletal elements in the sample set, the number of individuals sampled is uncertain (minimum of three). We find that the Late Pleistocene European leopards share a more recent common ancestor with modern Asian leopards, than either does with African leopards (Fig. 2; Additional file 1: Figure S1). The estimated coalescence time between the European and Asian mitogenomes is 485 Ka, with a 95% CI of 301–675 Ka, which encompasses the age for the oldest fossil records of Europe. This suggests that the European leopards originate from the same out-of-Africa event as the Asian leopard lineages, rather than from an independent migration. The ancient Caucasian leopard lineage is sister to all modern samples from mainland Asia, with an estimated coalescence time of ~ 244 Ka (95% CI: 148–352 Ka).
Morphological studies have failed to differentiate Pleistocene Eurasian leopard fossils from their modern counterparts . Cranial and postcranial measurements from Middle Pleistocene remains from Poland (estimated to 569–125 Ka;  and from the Caucasus (less than 350 Ka; ) overlapped with the variation found in modern leopards, which has been interpreted as evidence for some degree of population continuity between those ancient and modern populations. The dates of these fossils are not older, and actually correspond fairly well to the estimated coalescence time between ancient European mitochondrial lineages and modern leopards (95% CI: 305–677 Ka) or between the ancient Caucasian individual and modern leopards (95% CI: 148–352 Ka). Thus, some degree of population continuity since the Late Pleistocene, as suggested by the morphology, is not in conflict with our molecular dating. Alleles of Pleistocene European leopards may thus survive in the modern Caucasian population. However, as our dataset does not include modern Caucasian individuals, nor ancient East Asian ones, a direct assessment of population continuity between the Middle Pleistocene and today can only be based on short mtDNA sequences [24, 27]; Additional file 1: Figure S2). The mtDNA haplotype network has only limited resolution as it is only based on 456 bp of mtDNA, but does not support that the ancient Caucasus leopard is particularly divergent from modern Persian leopards (P. p. saxicolor; “SAX”; Additional file 1: Figure S2). The discrepancy between the coalescent times reported for the divergence between Iranian and other Asian leopards (16–270 Ka; ), and that between our Caucasian sample and other lineages (148–352 Ka) could suggest that the ancient Caucasian leopard represents a distinct lineage that diverged from other Asian leopards prior to the Iranian leopards, although it should be noted that the credibility intervals overlap considerably. Tracing the potential contribution of the European leopard to modern Caucasian or other modern leopard populations will require the analysis of nuclear DNA from Pleistocene leopards. Unfortunately, the Pleistocene samples analysed in the present study were not suitable for the analysis of nuclear DNA due to their low endogenous DNA content, despite multiple re-sampling attempts .
Java and the Asian leopard
The Javan leopard has been assigned to the subspecies P. p. melas (“MEL”), based on both morphological and genetic distinctiveness of these island leopards from mainland Asian populations [19, 22, 36]. Previous studies found Javan leopards to represent a mitochondrial lineage that is sister to the entire diversity of mainland Asian leopards [19, 22]. Due to this position as sister lineage to all Asian leopards, it has been suggested that the Javan leopard represents a relict population founded from a Middle Pleistocene leopard population from the South East Asian mainland, after which it became isolated on the Sundaic islands [19, 22]. Deep divergence between mainland and Sundaic species or populations has also been found in a range of other animals including terrestrial mammals  such as tigers  and leopard cats , but also birds , bats , and frogs .
Under the assumption that both of our samples with unclear provenance (‘Sumatra’ and ‘South East Asia’, sample codes PP3 and PP35) do in fact represent the Javan population, the molecular dating performed with our mitogenome dataset suggests that the Javan leopards became established on the island some time after 375 Ka (95% CI: 230–524 Ka). This is in line with previous estimates (393–886 Ka; ). Lowered sea levels during glacial periods exposed portions of the Sunda shelf, creating land bridges connecting the mainland with the Sundaic islands , which would have allowed the colonisation of and subsequent isolation on the South East Asian islands [22, 36]. If there was any more recent gene flow between the islands and the mainland during Late Pleistocene periods of low sea levels , it is not expected to have been at high frequency, as there has been no evidence found in the mitochondrial data investigated until now. However, the sample in our dataset that was reported to be of Javan origin (sample PP32) was placed as a sister lineage to the mitogenome recovered from an individual from Thailand. Investigation of short mtDNA sequences confirmed that the mitochondrial haplotype from this putative Javan leopard is closely related to P. p. delacouri (“DEL”), the mainland South-East Asian leopard subspecies (Additional file 1: Figure S2). Trade routes between the South East Asian islands and the mainland have been reported to exist for almost 2000 years , so leopards could have been moved to Java via such human-mitigated routes – either by trading live animals or as hunting trophies. Alternatively, our results also fit a scenario in which leopards migrated to the South-East Asian islands from the mainland during the Late Pleistocene, supplementing the original locally surviving Javan mitochondrial haplotypes with mainland Asian haplotypes; the latter of which could have been lost during the drastic population declines during the past century . The coalescence time of this putative Javan leopard and the Thai individual is estimated to be 64 Ka (95% CI: 32–100 Ka), which could be consistent with a dispersal from the mainland to the Sundaic islands some time (e.g. during the LGM) after the Toba supervolcanic eruption that occurred 73 Ka . Additional sampling and/or the investigation of nuclear data is needed to further clarify the origin and relatedness between mainland and Javan leopards.
Our dataset further included a number of South and East Asian leopards; two Indian leopards (P. p. fusca), one Indochinese P. p. delacouri and three published North-Chinese leopards, P. p. japonensis (Table 1; “FUS”, “DEL” and “JAP”, respectively). The estimated coalescence time for these individuals is relatively recent (122 Ka; 95% CI: 73–178 Ka). This recent divergence time could point towards extensive population reductions on the Eurasian mainland during the Pleistocene, leading to a bottleneck and loss of divergent lineages for all Eurasian leopards except for Java and Europe. A loss of genetic diversity during the Pleistocene as suggested by recent mitogenome coalescence times has also been found for other carnivores (e.g. 151 Ka for leopard cats, 95% CI: 87–215; ), suggesting that this pattern may be detectable across other Asian carnivores. Considering the deeper coalescence time of all modern African mtDNA lineages (600 Ka), the African leopards are not likely to have experienced a bottleneck to the same extent as Asian leopards.
This study represents an example of how genetic data from ancient DNA can inform on the evolutionary history of species and thereby provide additional clues on identification of fossils. This may include those from which no genetic data can be recovered because they are beyond the expected range of DNA survival, such as the 2 Ma year old fossils from Pakistan. The molecular dating of mitochondrial lineages is in line with the age of the earliest unequivocal leopard fossils from Eurasia, and thus it encourages a re-evaluation of the older leopard-like fossils. Our results are therefore compatible with previous suggestions that older fossils should be assigned to other large-bodied felids. Furthermore, similar to many other Eurasian carnivores, our data suggests that with the extinction of European populations and the proposed Pleistocene population bottleneck in mainland Asia, leopards experienced not only a reduction of their geographical distribution but also a loss of mitochondrial lineages that, to date, have not been detected in the modern gene pool. Leopard dispersal is generally driven by males, and thus may not be detectable using mitochondrial DNA alone. Finally, although we do not observe any evidence in our data to suggest that ancient European mitochondrial lineages persist in modern Asian populations, it is possible that at least part of the genetic legacy of Pleistocene leopards survives in the modern nuclear genome, as was recently shown to be the case for archaic hominins [58, 59] and the extinct cave bear .
DNA extraction, library preparation and capture
All pre-PCR steps for ancient and historical bone samples were performed in dedicated cleanroom facilities of the Evolutionary Adaptive Genomics group at Potsdam University, with all standard ancient DNA precautions (e.g. decontamination procedures for all reagents and materials, negative controls for both extraction and library preparation). Extraction was performed following a protocol optimised for the retrieval of short DNA fragments . Illumina sequencing libraries were constructed following a double-stranded library preparation protocol for the historical samples [62, 63] and a single-stranded library protocol for the ancient samples, using UDG/Endonuclease VIII to excise uracils from the molecules . Optimal amplification cycle numbers for each library were estimated using qPCR and then used for dual-indexing PCR. Four historical samples yielded high endogenous content, allowing the mitogenome sequence to be recovered using shotgun sequencing (Table 1; Additional file 1: Table S1). For all remaining samples, mitogenome enrichment was performed using an in-solution capture approach with synthetic baits. A bait-set targeting felid mitochondrial DNA was designed by placing 52-mer probes (1 bp tiling) across the mitogenomes of five felids: leopard (Panthera pardus, EF551002 (RefSeq: NC_010641), lion (Panthera leo, KF776494), domestic cat (Felis catus, NC_001700), bobcat (Lynx rufus, NC_014456.1) and cheetah (Acinonyx jubatus, NC_005212). Bait sequences containing simple repeats longer than 24 bp were removed . Baits were synthesized, amplified and converted into biotinylated single-stranded DNA probes as described elsewhere . Capture was performed following the protocol in Horn (2012). For ancient samples, two serial captures were performed to improve the enrichment rates. For historical samples, either one or two captures were used (Table 1). Sequencing was performed on the Illumina NextSeq 500 platform using 75 bp paired-end sequencing, using custom sequencing primers for the single-stranded libraries [64, 67].
Raw sequences were trimmed and merged using SeqPrep with default parameters and a minimum length of 30 bp (available from https://github.com/jstjohn/SeqPrep). The trimmed and merged reads were then aligned to a leopard mitogenome sequence available from GenBank (Acc. Nr. KP202265) using the Burrows-Wheeler Aligner (BWA aln) v0.7.8 and samtools v1.19 [68, 69] with default parameters. Reads with low mapping quality (<Q30) were removed. Duplications were marked and removed taking both mapping coordinates into consideration, using MarkDupsByStartEnd.jar (http://github.com/dariober/Java-cafe/tree/master/MarkDupsByStartEnd). Summary statistics can be found in Additional file 1: Table S1. Mitogenome consensus sequences were retrieved using Geneious v7.0 , using a minimum sequence depth of 3× and a strict 90% majority rule for base calling. The resulting consensus sequences for each sample were combined with three mitogenome sequences available from GenBank at the time of analysis, and aligned using ClustalW with default parameters (; as implemented in Geneious). The control region, as well as any positions in the alignment that contained missing or ambiguous data, were removed to avoid any adverse influence of misalignments or numts. The resulting alignment (13,688 bp in length) was manually annotated in Geneious using the published sequence (Panthera pardus; Genbank Acc. Nr. KP202265) as reference.
PartitionFinder v1.1.1  was used to identify an optimal partitioning scheme from all possible combinations of tRNAs, ribosomal RNA genes and protein-coding genes, considering all substitution models available in BEAUti, using the Bayesian Information Criterion (BIC). The partitioned maximum likelihood tree was calculated using RAxML-HPC v8.2.4  black box version, with default substitution models for each partition, on the CIPRES Science Gateway (; Additional file 1: Figure S1). Molecular dating was performed using BEAST v1.8.2 . We adopted a two step strategy for the calibrated analyses: first, we generated a fossil-calibrated phylogeny using a total of seven fossil dates (Additional file 1: Table S2) under a Yule speciation model for 15 sequences from various Felidae species, including two of the most divergent leopard mitogenomes. The divergence estimates recovered for the two leopard mitogenomes (mean 0.81 Ma, standard deviation 0.12 Ma) was then applied as normal prior on the root height for the intraspecific analyses, and were run with a Bayesian Skyline coalescent population model. Preliminary analysis using a lognormal relaxed clock model failed to reject zero variation in substitution rates across branches of the phylogeny, and so a strict clock model was employed with an open uniform prior on the mean per-lineage substitution rate of 0 to 20% per million years. The MCMC chain was run for a sufficient number of generations to achieve convergence (to a maximum of 10 million, sampling every 10,000 states) and adequate posterior sampling of all parameters (ESS > 200), checked using Tracer v1.5 (available from http://beast.community/tracer). TreeAnnotator v1.8.2 was then used to remove the first 25% of trees as burnin (corresponding to 2500 trees) and extract the Maximum Clade Credibility (MCC) tree with nodes scaled to the median heights recovered by the posterior sample. The minimum-spanning network for the short mtDNA alignment (267 individuals, 456 bp in length) was generated using Popart (Additional file 1: Figure S2) .
Kilo annus; thousand years ago
Last Glacial Maximum
Mega annus; million years ago
Hofreiter M, Paijmans JLA, Goodchild H, Speller CF, Barlow A, Fortes GG, et al. The future of ancient DNA: technical advances and conceptual shifts. BioEssays. 2015;37:284–93.
Delsuc F, Gibb GC, Kuch M, Billet G, Hautier L, Southon J, et al. The phylogenetic affinities of the extinct glyptodonts. Curr Biol. 2016;26:R155–6.
Paijmans JLA, Barnett R, Gilbert MTP, Zepeda-Mendoza ML, Reumer JWF, de VJ, et al. Evolutionary history of saber-toothed cats based on ancient Mitogenomics. Curr Biol. 2017;27:3330–3336.e5.
Paijmans JLA, Gilbert MTP, Hofreiter M. mitogenomic analyses from ancient DNA. Mol Phylogenet Evol. 2013;69:404–16.
Westbury M, Baleka S, Barlow A, Hartmann S, JLA P, Kramarz A, et al. A mitogenomic timetree for Darwin’s enigmatic south American mammal Macrauchenia patachonica. Nat Commun. 2017;8 ncomms15951.
Diedrich CG. Late Pleistocene leopards across Europe – northernmost European German population, highest elevated records in the Swiss Alps, complete skeletons in the Bosnia Herzegowina Dinarids and comparison to the ice age cave art. Quat Sci Rev. 2013;76:167–93.
Jacobson AP, Gerngross P, Lemeris Jr. JR, Schoonover RF, Anco C, Breitenmoser-Würsten C, et al. Leopard (Panthera pardus) status, distribution, and the research efforts across its range. Peer J 2016;4:e1974.
Sommer RS, Benecke N. Late Pleistocene and Holocene development of the felid fauna (Felidae) of Europe: a review. J Zool. 2006;269:7–19.
Nowell K, Jackson P. Wild cats: status survey and conservation action plan. IUCN Gland; 1996.
Packer C, Kosmala M, Cooley HS, Brink H, Pintea L, Garshelis D, et al. Sport hunting, predator control and conservation of large carnivores. PLoS One. 2009;4:e5941.
Rostro-García S, Kamler JF, Ash E, Clements GR, Gibson L, Lynam AJ, et al. Endangered leopards: range collapse of the Indochinese leopard (Panthera pardus delacouri) in Southeast Asia. Biol Conserv. 2016;201:293–300.
Swanepoel LH, Somers MJ, Van Hoven W, Schiess-Meier M, Owen G, Snyman A, et al. Survival rates and causes of mortality of leopards Panthera pardus in southern Africa; 2015. https://doi.org/10.1017/S0030605313001282.
Ebrahimi A, Farashi A, Rashki A. Habitat suitability of Persian leopard Panthera pardus saxicolor in Iran in future. Environ Earth Sci. 2017;76:697.
Farashi A, Shariati M. Evaluation of the role of the national parks for Persian leopard Panthera pardus saxicolor habitat conservation (case study: Tandooreh National Park, Iran). Mammal Res. 2018:1–8.
Kittle AM, Watson AC, Cushman SA, Macdonald DW. Forest cover and level of protection influence the island-wide distribution of an apex carnivore and umbrella species, the Sri Lankan leopard Panthera pardus kotiya. Biodivers Conserv. 2018;27:235–63.
Sandom CJ, Faurby S, Svenning J-C, Burnham D, Dickman A, Hinks AE, et al. Learning from the past to prepare for the future: felids face continued threat from declining prey. Ecography. 2018;41:140–52.
Wolf C, Ripple WJ. Prey depletion as a threat to the world’s large carnivores. R Soc Open Sci. 2016;3:160252.
Stein AB, Athreya V, Gerngross P, Balme G, Henschel P, Karanth U, et al. Panthera pardus (errata version published in 2016). The IUCN Red List of Threatened Species 2016. 2016;e.T15954A102421779.
Uphyrkina O, Johnson WE, Quigley H, Miquelle D, Marker L, Bush M, et al. Phylogenetics, genome diversity and origin of modern leopard, Panthera pardus. Mol Ecol. 2001;10:2617–33.
Werdelin L, Lewis ME. Plio-Pleistocene Carnivora of eastern Africa: species richness and turnover patterns. Zool J Linnean Soc 2005;144:121–44.
Hemmer H. Fossil history of living Felidae. Carnivore. 1979;2:58–61.
Wilting A, Patel R, Pfestorf H, Kern C, Sultan K, Ario A, et al. Evolutionary history and conservation significance of the Javan leopard Panthera pardus melas. J Zool. 2016;:n/a-n/a.
Anco C, Kolokotronis S-O, Henschel P, Cunningham SW, Amato G, Hekkala E. Historical mitochondrial diversity in African leopards (Panthera pardus) revealed by archival museum specimens. Mitochondrial DNA Part A. 2017;0:1–19.
Farhadinia MS, Farahmand H, Gavashelishvili A, Kaboli M, Karami M, Khalili B, et al. Molecular and craniological analysis of leopard, Panthera pardus (Carnivora: Felidae) in Iran: support for a monophyletic clade in Western Asia. Biol J Linn Soc. 2015;114:721–36.
Miththapala S, Seidensticker J, O’Brien SJ. Phylogeographic subspecies recognition in leopards (Panthera pardus): molecular genetic variation. Conserv Biol. 1996;10:1115–32.
Ropiquet A, Knight AT, Born C, Martins Q, Balme G, Kirkendall L, et al. Implications of spatial genetic patterns for conserving African leopards. C R Biol. 2015;338:728–37.
Rozhnov VV, Lukarevskiy VS, Sorokin PA. Application of molecular genetic characteristics for reintroduction of the leopard (Panthera pardus L., 1758) in the Caucasus. Dokl Biol Sci. 2011;437:97–102.
Spong G, Johansson M, Björklund M. High genetic variation in leopards indicates large and long-term stable effective population size. Mol Ecol. 2000;9:1773–82.
Uphyrkina O, Miquelle D, Quigley H, Driscoll C, O’Brien SJ. Conservation genetics of the far eastern leopard (Panthera pardus orientalis). J Hered. 2002;93:303–11.
Turner A, Antón M. The big cats and their fossil relatives: an illustrated guide to their evolution and natural history. In: Columbia University press; 1997.
Werdelin L, Dehghani R. Carnivora. In: Paleontology and geology of Laetoli: human evolution in context. Dordrecht: Springer; 2011. p. 189–232. https://doi.org/10.1007/978-90-481-9962-4_8.
Hemmer H, Kahlke R-D, Vekua AK. The Old World puma - Puma pardoides (Owen, 1846) (Carnivora: Felidae) - in the lower Villafranchian (upper Pliocene) of Kvabebi (East Georgia, Transcaucasia) and its evolutionary and biogeographical significance. Neues Jahrb Geol Palaontol Abh. 2004;233:197–231.
Werdelin L, Yamaguchi N, Johnson WE, O’Brien SJ. Phylogeny and evolution of cats (Felidae). Biol Conserv Wild Felids Oxf. 2010:59–82.
Harrison T. Paleontology and geology of Laetoli: human evolution in context: volume 2: fossil hominins and the associated Fauna. Springer Science & Business Media; 2011.
Baryshnikov GF. Late Pleistocene Felidae remains (Mammalia, Carnivora) from Geographical Society Cave in the Russian Far East. Proc Zool Inst RAS. 2016;320:84–120.
Meijaard E. Biogeographic history of the Javan leopard Panthera pardus based on a Craniometric analysis. J Mammal. 2004;85:302–10.
Baryshnikov GF. Pleistocene Felidae (Mammalia, Carnivora) from the Kudaro paleolithic cave sites in the Caucasus. In: Proceedings of the Zoological Institute Russian Academy of Science; 2011. p. 19.
Kurtén B. Pleistocene mammals of Europe. Transaction Publishers; 1968.
Marciszak A, Krajcarz MT, Krajcarz M, Stefaniak K. The first record of leopard Panthera pardus Linnaeus, 1758 from the Pleistocene of Poland. Acta Zool Cracoviensia - Ser Vertebr. 2011;54:39–46.
Ghezzo E, Rook L. The remarkable Panthera pardus (Felidae, Mammalia) record from Equi (Massa, Italy): taphonomy, morphology, and paleoecology. Quat Sci Rev. 2015;110(Supplement C):131–51.
Nagel D. Panthera pardus vraonensis n. ssp., a new leopard from the Pleistocene of Vraona/Greece. (With 5 figures and 4 tables). Neues Jahrb Für Geol Palaontologie Monatshefte. 1999:129–50.
Sauqué V, Rabal-Garcés R, Cuenca-Bescós G. Carnivores from Los Rincones, a leopard den in the highest mountain of the Iberian range (Moncayo, Zaragoza, Spain). Hist Biol. 2016;28:479–506.
Sabol M, Persico D, Troco E. First fossil record of leopard-like felid (Panthera cf. pardus) from alluvial deposits of the Po River in northern Italy. Quat Int. 2017. https://doi.org/10.1016/j.quaint.2016.12.036.
Baca M, Popović D, Stefaniak K, Marciszak A, Urbanowski M, Nadachowski A, et al. Retreat and extinction of the Late Pleistocene cave bear (Ursus spelaeus sensu lato). Naturwissenschaften. 2016;103. https://doi.org/10.1007/s00114-016-1414-8.
Pacher M, Stuart AJ. Extinction chronology and palaeobiology of the cave bear (Ursus spelaeus). Boreas. 2009;38:189–206.
Stuart AJ, Lister AM. Patterns of Late Quaternary megafaunal extinctions in Europe and northern Asia. Cour-Forschungsinstitut Senckenberg. 2007;259:287.
Alberti F, Gonzalez J, Paijmans JLA, Basler N, Preick M, Henneberger K, et al. Optimised DNA sampling of ancient bones using computed tomography (CT) scans. Mol Ecol Resour. 2018. https://doi.org/10.1111/1755-0998.12911.
Woodruff DS, Turner LM. The Indochinese–Sundaic zoogeographic transition: a description and analysis of terrestrial mammal species distributions. J Biogeogr. 2009;36:803–21.
Wilting A, Courtiol A, Christiansen P, Niedballa J, Scharf AK, Orlando L, et al. Planning tiger recovery: understanding intraspecific variation for effective conservation. Sci Adv. 2015;1:e1400175.
Patel RP, Wutke S, Lenz D, Mukherjee S, Ramakrishnan U, Veron G, et al. Genetic structure and phylogeography of the leopard cat (Prionailurus bengalensis) inferred from mitochondrial genomes. J Hered. 2017;108:349–60.
Hughes JB, Round PD, Woodruff DS. The Indochinese–Sundaic faunal transition at the Isthmus of Kra: an analysis of resident forest bird species distributions. Journal of Biogeography. 2003;30:569–80.
Hughes AC, Satasook C, Bates PJJ, Bumrungsri S, Jones G. Explaining the causes of the zoogeographic transition around the Isthmus of Kra: using bats as a case study. Journal of Biogeography. 2011;38:2362–72.
Inger RF, Voris HK. The biogeographical relations of the frogs and snakes of Sundaland. Journal of Biogeography. 2001;28:863–91.
Lohman DJ, de Bruyn M, Page T, von Rintelen K, Hall R, Ng PKL, et al. Biogeography of the Indo-Australian archipelago. Annu Rev Ecol Evol Syst. 2011;42:205–26.
Voris HK. Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. J Biogeogr. 2000;27:1153–67.
Hall KR. A history of early Southeast Asia: maritime trade and societal development, 100–1500. Lanham, Md: Rowman & Littlefield; 2011.
Williams MAJ, Ambrose SH, van der Kaars S, Ruehlemann C, Chattopadhyaya U, Pal J, et al. Environmental impact of the 73ka Toba super-eruption in South Asia. Palaeogeogr Palaeoclimatol Palaeoecol. 2009;284:295–314.
Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, et al. A draft sequence of the Neandertal genome. Science. 2010;328:710–22.
Meyer M, Kircher M, Gansauge M-T, Li H, Racimo F, Mallick S, et al. A high-coverage genome sequence from an archaic Denisovan individual. Science. 2012;338:222–6.
Barlow A, Cahill JA, Hartmann S, Theunert C, Xenikoudakis G, Fortes GG, et al. Partial genomic survival of cave bears in living brown bears. Nat Ecol Evol. 2018. https://doi.org/10.1038/s41559-018-0654-8.
Dabney J, Knapp M, Glocke I, Gansauge M-T, Weihmann A, Nickel B, et al. Complete mitochondrial genome sequence of a middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc Natl Acad Sci. 2013;110:15758–63.
Fortes GG, Paijmans JLA. In: Kroneis T, editor. Analysis of whole mitogenomes from ancient samples. USA: Whole genome amplification. Humana Press; 2016.
Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010;2010. https://doi.org/10.1101/pdb.prot5448.
Gansauge M-T, Meyer M. Single-stranded DNA library preparation for the sequencing of ancient or damaged DNA. Nat Protoc. 2013;8:737–48.
Slon V, Glocke I, Barkai R, Gopher A, Hershkovitz I, Meyer M. Mammalian mitochondrial capture, a tool for rapid screening of DNA preservation in faunal and undiagnostic remains, and its application to Middle Pleistocene specimens from Qesem cave (Israel). Quat Int. 2015. https://doi.org/10.1016/j.quaint.2015.03.039.
Fu Q, Meyer M, Gao X, Stenzel U, Burbano HA, Kelso J, et al. DNA analysis of an early modern human from Tianyuan cave, China. Proc Natl Acad Sci. 2013;110:2223–7.
Paijmans JLA, Baleka S, Henneberger K, Taron UH, Trinks A, Westbury MV, et al. Sequencing single-stranded libraries on the Illumina NextSeq 500 platform. In: ArXiv171111004 Q-Bio; 2017. http://arxiv.org/abs/1711.11004. Accessed 5 Dec 2017.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25:1754–60.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.
Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29:1695–701.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014; btu033.
Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE) 2010; 2010. p. 1–8.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Leigh JW, Bryant D. Popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.
Joger U, Rosendahl W. The Rübeland Caves (Harz Mountains) - Historical Excavations and Modern Analyses. Braunschweiger Naturkundliche Schriften. 2012;11:55–68.
Lei W, XiaoBing W, Zhu L, and Jiang Z. Mitogenomic Analysis of the Genus Panthera. Science China Life Sciences. 2011;54(10):917–30.
Dou H, Feng L, Xiao W, Wang T. The Complete Mitochondrial Genome of the North Chinese Leopard (Panthera Pardus Japonensis)”. Mitochondrial DNA 2016,27(2):1167–68. https://doi.org/10.3109/19401736.2014.936421.
Gang L, Davis BW, Eizirik E, Murphy WJ. Phylogenomic Evidence for Ancient Hybridization in the Genomes of Living Cats (Felidae). Genome Research 2016,26(1):1–11. https://doi.org/10.1101/gr.186668.114.
We would like to acknowledge Peter Rask Møller and Hans Baagøe from the Natural History Museum of Denmark for providing access to samples from their zoological collection. We thank Mathias Stiller for help with the capture bait design. We also acknowledge Tom Gilbert for supervision of RWH, and for providing helpful comments on earlier drafts of the manuscript.
This work was supported by the European Research Council (consolidator grant GeneFlow #310763 to M.H.). The NVIDIA TITAN-X GPU used for BEAST analyses was kindly donated by the NVIDIA Corporation. The funding bodies had no role in study design, analysis and interpretation, or writing the manuscript.
Availability of data and materials
Leopard mitogenome consensus sequences will be made available on GenBank (NCBI Accession Numbers MH588611-MH588632). The alignment used for phylogenetic inference and the BEAST XML input file used for fossil calibration analyses will be made available upon request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Maximum likelihood phylogeny, with different Panthera species as outgroup. Figure S2: Minimum-spanning network of short (456 bp) mtDNA sequences (ND5 gene), including previously published data (based on 267 sequences in total). Number of substitutions are indicated as tick-marks on the branches connecting the haplotypes. Colours indicate the subspecies [19, following 25]. Table S1: Summarised sequence statistics for samples included in our study. Table S2: Fossil constraints and calibration priors used in the time-calibrated BEAST analysis performed for Felidae alignment. The resulting root age was then applied as calibration for the leopards-only phylogeny. Table S3: Radiocarbon dating (14C) information. (PDF 1852 kb)