Historical mitochondrial diversity in African leopards (Panthera pardus) revealed by archival museum specimens

Unravelling the maternal evolutionary history of the African leopard (Panthera pardus pardus)

The African leopard (Panthera pardus pardus) has lost a significant proportion of its historical range, notably in north-western Africa and South Africa. Recent studies have explored the genetic diversity and population structure of African leopards across the continent. A notable genetic observation is the presence of two divergent mitochondrial lineages, PAR-I and PAR-II. Both lineages appeared to be distributed widely, with PAR-II frequently found in southern Africa. Until now, no study has attempted to date the emergence of either lineage, assess haplotype distribution, or explore their evolutionary histories in any detail. To investigate these underappreciated questions, we compiled the largest and most geographically representative leopard data set of the mitochondrial NADH-5 gene to date. We combined samples (n = 33) collected in an altitudinal transect across the Mpumalanga province of South Africa, where two populations of leopard are known to be in genetic contact, with previously published sequences of African leopard (n = 211). We estimate that the maternal PAR-I and PAR-II lineages diverged approximately 0.7051 (0.4477-0.9632) million years ago (Ma). Through spatial and demographic analyses, we show that while PAR-I underwent a mid-Pleistocene population expansion resulting in several closely related haplotypes with little geographic structure across much of its range, PAR-II remained at constant size and may even have declined slightly in the last 0.1 Ma. The higher genetic drift experienced within PAR-II drove a greater degree of structure with little haplotype sharing and unique haplotypes in central Africa, the Cape, KwaZulu-Natal and the South African Highveld. The phylogeographic structure of PAR-II, with its increasing frequency southward and its exclusive occurrence in south-eastern South Africa, suggests that this lineage may have been isolated in South Africa during the mid-Pleistocene. This hypothesis is supported by historical changes in paleoclimate that promoted intense aridification around the Limpopo Basin between 1.0-0.6 Ma, potentially reducing gene flow and promoting genetic drift. Interestingly, we ascertained that the two nuclear DNA populations identified by a previous study as East and West Mpumalanga correspond to PAR-I and PAR-II, respectively, and that they have come into secondary contact in the Lowveld region of South Africa. Our results suggest a subdivision of African leopard mtDNA into two clades, with one occurring almost exclusively in South Africa, and we identify the potential environmental drivers of this observed structure. We caution that our results are based on a single mtDNA locus, but it nevertheless provides a hypothesis that can be further tested with a dense sample of nuclear DNA data, preferably whole genomes. If our interpretation holds true, it would provide the first genetic explanation for the smaller observed size of leopards at the southernmost end of their range in Africa.

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 unique single nucleotide polymorphism in Agouti Signalling Protein (ASIP) gene changes coat colour of Sri Lankan leopard (Panthera pardus kotiya) to dark black

The Sri Lankan leopard (Panthera pardus kotiya) is an endangered subspecies restricted to isolated and fragmented populations in Sri Lanka. Among them, melanistic leopards have been recorded on a few occasions. Literature suggests the evolution of melanism several times in the Felidae family, with three species having distinct mutations. Nevertheless, the mutations or other variations in the remaining species, including Sri Lankan melanistic leopard, are unknown. We used reference-based assembled nuclear genomes of Sri Lankan wild type and melanistic leopards and de novo assembled mitogenomes of the same to investigate the genetic basis, adaptive significance, and evolutionary history of the Sri Lankan melanistic leopard. Interestingly, we identified a single nucleotide polymorphism in exon-4 Sri Lankan melanistic leopard, which may completely ablate Agouti Signalling Protein (ASIP) function. The wild type leopards in Sri Lanka did not carry this mutation, suggesting the cause for the occurrence of melanistic leopords in the population. Comparative analysis of existing genomic data in the literature suggests it as a P. p. kotiya specific mutation and a novel mutation in the ASIP-gene of the Felidae family, contributing to naturally occurring colour polymorphism. Our data suggested the coalescence time of Sri Lankan leopards at ~0.5 million years, sisters to the Panthera pardus lineage. The genetic diversity was low in Sri Lankan leopards. Further, the P. p. kotiya melanistic leopard is a different morphotype of the P. p. kotiya wildtype leopard resulting from the mutation in the ASIP-gene. The ability of black leopards to camouflage, along with the likelihood of recurrence and transfer to future generations, suggests that this rare mutation could be environment-adaptable.

A Camera-Trap Home-Range Analysis of the Indian Leopard (Panthera pardus fusca) in Jaipur, India

The suitability of the camera trap-retrap method was explored for identifying territories and studying the spatial distribution of leopards (Panthera pardus fusca) in the Jhalana Reserve Forest, Jaipur, India. Data from two years (November 2017 to November 2019, N = 23,208 trap-hours) were used to provide estimates of minimum home-range size and overlap. We conducted home-range analysis and estimation, using the minimum convex polygon (MCP) method with geographic information system (GIS) tools. We are aware of the limitations and advantages of camera trapping for long-term monitoring. However, the limitations of the research permit allowed only the use of camera traps to estimate the home ranges. A total of 25 leopards were identified (male = 8, female = 17). No territorial exclusivity was observed in either of the sexes. However, for seven females, we observed familial home-range overlaps wherein daughters established home ranges adjacent to or overlapping their natal areas. The median home range, as calculated from the MCP, was 305.9 ha for males and 170.3 ha for females. The median percentage overlap between males was 10.33%, while that between females was 3.97%. We concluded that camera trapping is an effective technique to map the territories of leopards, to document inter- and intraspecific behaviors, and to elucidate how familial relationships affect dispersal.

Is there low maternal genetic variation in West Asian populations of leopard?

Persian leopards Panthera pardus saxicolor have been extirpated from over 84% of their historic range and are now limited to rugged landscapes of West Asia and the Caucasus. Understanding and maintaining genetic diversity and population connectivity is important for preventing inbreeding and genetic drift, both of which can threaten population viability. All previous analyses of intraspecific genetic variation of West Asian leopards based on the NADH dehydrogenase subunit 5 gene have reported low mitogenomic diversity. In the current study, we sequenced 959 bp of the mtDNA cytochrome b gene to describe the spatial genetic structure of 22 wild Persian leopards across Iran, which hosts most of the subspecies extant range. The findings based on phylogenetic trees and median-joining network indicated that leopards from Iran formed a distinct subclade, i.e., P. p. saxicolor. The AMOVA analysis showed significant differentiation (88.55%) between the subclades of Persian leopards and other Asian leopards. The lowest levels of haplotype (0.247) and nucleotide (0.00078) diversity were estimated in Persian leopards from Iran. Mitochondrial genome sequencing revealed only two closely related haplotypes. There was no evidence for recent sudden demographic expansion scenario in Persian leopards. The low diversity in cytochrome b gene could potentially be brought about by selective pressure on mitochondria to adapt to oxidative stress and higher metabolic rates in cold environments.

Shotgun sequencing decades-old lichen specimens to resolve phylogenomic placement of type material

Natural history collections, including name-bearing type specimens, are an important source of genetic information. These data can be critical for appropriate taxonomic revisions in cases where the phylogenetic position of name-bearing type specimens needs to be identified, including morphologically cryptic lichen-forming fungal species. Here, we use high-throughput metagenomic shotgun sequencing to generate genome-scale data from decades-old (i.e., more than 30 years old) isotype specimens representing three vagrant taxa in the lichen-forming fungal genusRhizoplaca, including one species and two subspecies. We also use data from high-throughput metagenomic shotgun sequencing to infer the phylogenetic position of an enigmatic collection, originally identified asR. haydenii, that failed to yield genetic data via Sanger sequencing. We were able to construct a 1.64 Mb alignment from over 1200 single-copy nuclear gene regions for theRhizoplacaspecimens. Phylogenomic reconstructions recovered an isotype representingRhizoplaca haydeniisubsp.arbusculawithin a clade comprising other specimens identified asRhizoplaca haydeniisubsp.arbuscula, while an isotype ofR. idahoensiswas recovered within a clade with substantial phylogenetic substructure comprisingRhizoplaca haydeniisubsp.haydeniiand other specimens. Based on these data and morphological differences,Rhizoplaca haydeniisubsp.arbusculais elevated to specific rank asRhizoplaca arbuscula. For the enigmatic collection, we were able to assemble the nearly complete nrDNA cistron and over 50 Mb of the mitochondrial genome. Using these data, we identified this specimen as a morphologically deviant form representingXanthoparmeliaaff.subcumberlandia. This study highlights the power of high-throughput metagenomic shotgun sequencing in generating larger and more comprehensive genetic data from taxonomically important herbarium specimens.

Assessing subspecies status of leopards (Panthera pardus) of northern Pakistan using mitochondrial DNA

Despite being classified as critically endangered, little work has been done on leopard protection in Pakistan. Once widely present throughout this region, leopards are now sparsely distributed, and possibly extinct from much of their previously recorded habitat. While leopards show morphological and genetic variation across their species range worldwide, resulting in the classification of nine different subspecies, the leopard genetic structure across Pakistan is unknown, with previous studies including only a very limited sampling. To clarify the genetic status of leopards in Pakistan we investigated the sequence variation in the subunit 5 of the mitochondrial gene NADH from 43 tissue samples and compared it with 238 sequences available from online databases. Phylogenetic analysis clearly separates the Pakistani leopards from the African and Arabian clades, confirming that leopards from Pakistan are members of the Asian clade. Furthermore, we identified two separate subspecies haplotypes within our dataset: P. p. fusca (N = 23) and P. p. saxicolor (N = 12).

Historical biogeography of the leopard (Panthera pardus) and its extinct Eurasian populations

Background

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?

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1268-0) contains supplementary material, which is available to authorized users.

No longer locally extinct? Tracing the origins of a lion (Panthera leo) living in Gabon

Lions (Panthera leo) are of particular conservation concern due to evidence of recent, widespread population declines in what has hitherto been seen as a common species, robust to anthropogenic disturbance. Here we use non-invasive methods to recover complete mitochondrial genomes from single hair samples collected in the field in order to explore the identity of the Gabonese Plateaux Batéké lion. Comparison of the mitogenomes against a comprehensive dataset of African lion sequences that includes relevant geographically proximate lion populations from both contemporary and ancient sources, enabled us to identify the Plateaux Batéké lion as a close maternal relative to now extirpated populations found in Gabon and nearby Congo during the twentieth century, and to extant populations of Southern Africa. Our study demonstrates the relevance of ancient DNA methods to field conservation work, and the ability of trace field samples to provide copious genetic information about free-ranging animals.