Implications of spatial genetic patterns for conserving African leopards

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.

Leopard density and the ecological and anthropogenic factors influencing density in a mixed-use landscape in the Western Cape, South Africa

Large carnivores face numerous threats, including habitat loss and fragmentation, direct killing, and prey depletion, leading to significant global range and population declines. Despite such threats, leopards (Panthera pardus) persist outside protected areas throughout most of their range, occupying diverse habitat types and land uses, including peri-urban and rural areas. Understanding of leopard population dynamics in mixed-use landscapes is limited, especially in South Africa, where the majority of leopard research has focused on protected areas. We use spatially explicit capture-recapture models to estimate leopard density across a mixed-use landscape of protected areas, farmland, and urban areas in the Overberg region of the Western Cape, South Africa. Data from 86 paired camera stations provided 221 independent captures of 25 leopards at 50 camera trap stations with a population density estimate of 0.64 leopards per 100 km2 (95% CI: 0.55-0.73). Elevation, terrain ruggedness, and vegetation productivity were important drivers of leopard density in the landscape, being highest on elevated remnants of natural land outside of protected areas. These results are similar to previous research findings in other parts of the Western Cape, where high-lying natural vegetation was shown to serve as both a refuge and a corridor for leopard movement in otherwise transformed landscapes. Given the low leopard density and the prevalence of transformed land intermixed with patches of more suitable leopard habitat, prioritising and preserving connectivity for leopards is vital in this shared landscape. Ecological corridors should be developed in partnership with private landowners through an inclusive and multifaceted conservation strategy which also incorporates monitoring of and rapid mitigation of emerging threats to leopards.

Complex origins indicate a potential bridgehead introduction of an emerging amphibian invader (Eleutherodactylus planirostris) in China

Identifying the origins of established alien species is important to prevent new introductions in the future. The greenhouse frog (Eleutherodactylus planirostris), native to Cuba, the Bahamas, and the Cayman Islands, has been widely introduced to the Caribbean, North and Central America, Oceania and Asia. This invasive alien amphibian was recently reported in Shenzhen, China, but the potential introduction sources remain poorly understood. Based on phylogenetic analysis using mitochondrial 16S, COI and CYTB sequences, we detected a complex introduction origin of this species, which may be from Hong Kong, China, the Philippines, Panama and Florida, USA, all pointing to a bridgehead introduction. In addition, the nursery trade between the four countries or regions and mainland China from 2011 to 2020 was also significantly higher than other areas with less likelihood of introductions, which supported the molecular results. Our study provides the first genetic evidence of the potential sources of this emerging amphibian invader in mainland China, which may help develop alien species control strategies in the face of growing trade through globalization.

Molecular tracking and prevalence of the red colour morph restricted to a harvested leopard population in South Africa

The red leopard (Panthera pardus) colour morph is a colour variant that occurs only in South Africa, where it is confined to the Central Bushveld bioregion. Red leopards have been spreading over the past 40 years, which raises the speculation that the prevalence of this phenotype is related to low dispersal of young individuals owing to high off‐take in the region. Intensive selective hunting tends to remove large resident male leopards from the breeding population, which gives young male leopards the chance to mate with resident female leopards that are more likely to be their relatives, eventually increasing the frequency of rare genetic variants. To investigate the genetic mechanisms underlying the red coat colour morph in leopards, and whether its prevalence in South Africa relates to an increase in genetic relatedness in the population, we sequenced exons of six coat colour‐associated genes and 20 microsatellite loci in twenty Wild‐type and four red leopards. The results were combined with demographic data available from our study sites. We found that red leopards own a haplotype in homozygosity identified by two SNPs and a 1 bp deletion that causes a frameshift in the tyrosinase‐related protein 1 (TYRP1), a gene known to be involved in the biosynthesis of melanin. Microsatellite analyses indicate clear signs of a population bottleneck and a relatedness of 0.11 among all pairwise relationships, eventually supporting our hypothesis that a rare colour morph in the wild has increased its local frequency due to low natal dispersal, while subject to high human‐induced mortality rate.

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.

Unsustainable anthropogenic mortality disrupts natal dispersal and promotes inbreeding in leopards

Anthropogenic mortality of wildlife is typically inferred from measures of the absolute decline in population numbers. However, increasing evidence suggests that indirect demographic effects including changes to the age, sex, and social structure of populations, as well as the behavior of survivors, can profoundly impact population health and viability. Specifically, anthropogenic mortality of wildlife (especially when unsustainable) and fragmentation of the spatial distribution of individuals (home-ranges) could disrupt natal dispersal mechanisms, with long-term consequences to genetic structure, by compromising outbreeding behavior and gene flow. We investigate this threat in African leopards (Panthera pardus pardus), a polygynous felid with male-biased natal dispersal. Using a combination of spatial (home-range) and genetic (21 polymorphic microsatellites) data from 142 adult leopards, we contrast the structure of two South African populations with markedly different histories of anthropogenically linked mortality. Home-range overlap, parentage assignment, and spatio-genetic autocorrelation together show that historical exploitation of leopards in a recovering protected area has disrupted and reduced subadult male dispersal, thereby facilitating opportunistic male natal philopatry, with sons establishing territories closer to their mothers and sisters. The resultant kin-clustering in males of this historically exploited population is comparable to that of females in a well-protected reserve and has ultimately led to localized inbreeding. Our findings demonstrate novel evidence directly linking unsustainable anthropogenic mortality to inbreeding through disrupted dispersal in a large, solitary felid and expose the genetic consequences underlying this behavioral change. We therefore emphasize the importance of managing and mitigating the effects of unsustainable exploitation on local populations and increasing habitat fragmentation between contiguous protected areas by promoting in situ recovery and providing corridors of suitable habitat that maintain genetic connectivity.

Genetic analyses reveal population structure and recent decline in leopards (Panthera pardus fusca) across the Indian subcontinent

Background

Large carnivores maintain the stability and functioning of ecosystems. Currently, many carnivore species face declining population sizes due to natural and anthropogenic pressures. The leopard, Panthera pardus, is probably the most widely distributed and highly adaptable large felid globally, still persisting in most of its historic range. However, we lack subspecies-level data on country or regional scale on population trends, as ecological monitoring approaches are difficult to apply on such wide-ranging species. We used genetic data from leopards sampled across the Indian subcontinent to investigate population structure and patterns of demographic decline.

Methods

We collected faecal samples from the Terai-Arc landscape of northern India and identified 56 unique individuals using a panel of 13 microsatellite markers. We merged this data with already available 143 leopard individuals and assessed genetic structure at country scale. Subsequently, we investigated the demographic history of each identified subpopulations and compared genetic decline analyses with countrywide local extinction probabilities.

Results

Our genetic analyses revealed four distinct subpopulations corresponding to Western Ghats, Deccan Plateau-Semi Arid, Shivalik and Terai region of the north Indian landscape, each with high genetic variation. Coalescent simulations with microsatellite loci revealed a possibly human-induced 75–90% population decline between ∼120–200 years ago across India. Population-specific estimates of genetic decline are in concordance with ecological estimates of local extinction probabilities in these subpopulations obtained from occupancy modeling of the historic and current distribution of leopards in India.

Conclusions

Our results confirm the population decline of a widely distributed, adaptable large carnivore. We re-iterate the relevance of indirect genetic methods for such species in conjunction with occupancy assessment and recommend that detailed, landscape-level ecological studies on leopard populations are critical to future conservation efforts. Our approaches and inference are relevant to other widely distributed, seemingly unaffected carnivores such as the leopard.

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).

Genetic and spatial characterization of the red fox (Vulpes vulpes) population in the area stretching between the Eastern and Dinaric Alps and its relationship with rabies and canine distemper dynamics

Information on the population dynamics of a reservoir species have been increasingly adopted to understand and eventually predict the dispersal patterns of infectious diseases throughout an area. Although potentially relevant, to date there are no studies which have investigated the genetic structure of the red fox population in relation to infectious disease dynamics. Therefore, we genetically and spatially characterised the red fox population in the area stretching between the Eastern and Dinaric Alps, which has been affected by both distemper and rabies at different time intervals. Red foxes collected from north-eastern Italy, Austria, Slovenia and Croatia between 2006–2012, were studied using a set of 21 microsatellite markers. We confirmed a weak genetic differentiation within the fox population using Bayesian clustering analyses, and we were able to differentiate the fox population into geographically segregated groups. Our finding might be due to the presence of geographical barriers that have likely influenced the distribution of the fox population, limiting in turn gene flow and spread of infectious diseases. Focusing on the Italian red fox population, we observed interesting variations in the prevalence of both diseases among distinct fox clusters, with the previously identified Italy 1 and Italy 2 rabies as well as distemper viruses preferentially affecting different sub-groups identified in the study. Knowledge of the regional-scale population structure can improve understanding of the epidemiology and spread of diseases. Our study paves the way for an integrated approach for disease control coupling pathogen, host and environmental data to inform targeted control programs in the future.

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.

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

Once found throughout Africa and Eurasia, the leopard (Panthera pardus) was recently uplisted from Near Threatened to Vulnerable by the International Union for the Conservation of Nature (IUCN). Historically, more than 50% of the leopard's global range occurred in continental Africa, yet sampling from this part of the species' distribution is only sparsely represented in prior studies examining patterns of genetic variation at the continental or global level. Broad sampling to determine baseline patterns of genetic variation throughout the leopard's historical distribution is important, as these measures are currently used by the IUCN to direct conservation priorities and management plans. By including data from 182 historical museum specimens, faecal samples from ongoing field surveys, and published sequences representing sub-Saharan Africa, we identify previously unrecognized genetic diversity in African leopards. Our mtDNA data indicates high levels of divergence among regional populations and strongly differentiated lineages in West Africa on par with recent studies of other large vertebrates. We provide a reference benchmark of genetic diversity in African leopards against which future monitoring can be compared. These findings emphasize the utility of historical museum collections in understanding the processes that shape present biodiversity. Additionally, we suggest future research to clarify African leopard taxonomy and to differentiate between delineated units requiring monitoring or conservation action.

The importance of correctly identifying the process responsible for spatial genetic structure in Leopard: A response to McManus and Smuts (2016)

Microsatellite analyses suggest that spatial genetic structure among six leopard-sampling sites in southern Africa is the result of isolation by distance.

Low microsatellites used to investigate leopard genetic structure severely restricts the results by Ropiquet et al. (2015) to infer population structure for managers

Low microsatellites used to investigate leopard genetic structure severely restricts the results by Ropiquet et al. (2015) to infer population structure for managers.