USE OF SOLENSIA (FRUNEVETMAB) IN NONDOMESTIC FELIDS (PANTHERA TIGRIS, LYNX RUFUS, PANTHERA UNCIA) GUIDED BY PROTEIN SEQUENCE ANALYSIS.

Osteoarthritis and degenerative joint disease affect many species of nondomestic felid and negatively impact quality of life in managed care settings. Previously, pain management options were limited because of the frequency of comorbidities such as renal disease and common difficulty encountered in medication administration. SolensiaTM (frunevetmab) is a felinized monoclonal antibody that binds to nerve growth factor in arthritic joints, thereby inhibiting pain response cascades, and produces marked improvement in clinical signs and quality of life in domestic cats diagnosed with osteoarthritis. Protein sequence analysis was performed to inform the application of frunevetmab in four individuals from three nondomestic felid species (Panthera tigris, Lynx rufus, and Panthera uncia) diagnosed with degenerative joint disease to predict safety, contraindications, and likely response to treatment. Patients were then treated utilizing doses extrapolated from domestic cats and following manufacturer guidelines for administration, and clinical response was evaluated over a minimum 2-mon period. Significant improvement was noted in clinical signs in all four animals, resulting in marked improvement in mobility, lameness, activity level, demeanor, and overall quality of life. Frunevetmab presents an excellent adjunctive therapeutic alternative for nondomestic felids suffering from degenerative joint disease and may complement or decrease the need for other conventional therapeutic regimens.

Different coexistence patterns between apex carnivores and mesocarnivores based on temporal, spatial, and dietary niche partitioning analysis in Qilian Mountain National Park, China.

Carnivores play key roles in maintaining ecosystem structure and function as well as ecological processes. Understanding how sympatric species coexist in natural ecosystems is a central research topic in community ecology and biodiversity conservation. In this study, we explored intra- and interspecific niche partitioning along spatial, temporal, and dietary niche partitioning between apex carnivores (wolf Canis lupus, snow leopard Panthera uncia, Eurasian lynx Lynx lynx) and mesocarnivores (Pallas's cat Otocolobus manul, red fox Vulpes vulpes, Tibetan fox Vulpes ferrilata) in Qilian Mountain National Park, China, using camera trapping data and DNA metabarcoding sequencing data. Our study showed that apex carnivore species had more overlap temporally (coefficients of interspecific overlap ranging from 0.661 to 0.900) or trophically (Pianka's index ranging from 0.458 to 0.892), mesocarnivore species had high dietary overlap with each other (Pianka's index ranging from 0.945 to 0.997), and apex carnivore and mesocarnivore species had high temporal overlap (coefficients of interspecific overlap ranging from 0.497 to 0.855). Large dietary overlap was observed between wolf and snow leopard (Pianka's index = 0.892) and Pallas's cat and Tibetan fox (Pianka's index = 0.997), suggesting the potential for increased resource competition for these species pairs. We concluded that spatial niche partitioning is likely to key driver in facilitating the coexistence of apex carnivore species, while spatial and temporal niche partitioning likely facilitate the coexistence of mesocarnivore species, and spatial and dietary niche partitioning facilitate the coexistence between apex and mesocarnivore species. Our findings consider partitioning across temporal, spatial, and dietary dimensions while examining diverse coexistence patterns of carnivore species in Qilian Mountain National Park, China. These findings will contribute substantially to current understanding of carnivore guilds and effective conservation management in fragile alpine ecosystems.

Impact of carbon neutralization policy on the suitable habitat distribution of the North China leopard.

The Chinese government has introduced a carbon neutral policy to cope with the rapid changes in the global climate. It is not clear what impact this policy will have on wildlife. Therefore, this study analyzed the suitable habitat distribution of China's unique leopard subspecies in northern Shaanxi, and simulated the potential suitable habitat distribution under different carbon emission scenarios at two time points of future carbon peak and carbon neutralization. We found that in the future SSPs 126 scenario, the suitable habitat area and the number of suitable habitat patches of North China leopard will continue to increase. With the increase of carbon emissions, it is expected that the suitable habitat of North China leopard will continue to be fragmented and shifted. When the annual average temperature is lower than 8 °C, the precipitation seasonality is 80-90 mm and the precipitation of the warmest quarter is greater than 260 mm, the probability of occurrence of North China leopard is higher. The increase in carbon emissions will lead to the reduction, migration, and fragmentation of the suitable habitat distribution of the North China leopard. Carbon neutrality policies can protect suitable wild habitats. In the future, the impact of carbon neutrality policies on future wildlife habitat protection should be carried out in depth to effectively promote the construction of wildlife protection projects.

Home range, habitat use and capture-release of translocated leopards in Gir landscape, Gujarat, India.

Understanding the spatial ecology of translocated leopards (Panthera pardus fusca) is crucial for their conservation and the effective assessment of conflict management strategies. We investigated the home range and habitat preferences of five radio-collared leopards (n = 5; 2 males; 3 females) in the Gir landscape. Additionally, we examined the usefulness of the capture-release strategy for these animals. We assessed home range and habitat selection using kernel density estimation (at 95% and 50% levels) and compositional analysis. Our findings revealed that leopards exhibited distinct patterns of movement, often returning to their original capture site or nearby locations or exploring new areas within 3 to 25 days, covering distances ranging from 48 to 260 km. The average home range (95% FK) was estimated at 103.96±36.37 (SE) km2, with a core area usage (50% FK) of 21.38±5.95 km2. Seasonally, we observed the largest home ranges during summer and the smallest during winter. Males exhibited larger home ranges (95% FK, 151±64.28 km2) compared to females (56.18±14.22 km2). The habitat analysis indicated that agricultural areas were consistently preferred in the multi-use landscape at the 2nd order habitat selection level. Additionally, habitat around water bodies was highly favoured at the 3rd order, with distinct variations in habitat selection observed during day and night. This study highlights the significance of riverine and scrubland habitats, as leopards exhibited strong preferences for these habitats within their home ranges. We emphasize the importance of conserving natural habitat patches, particularly those surrounding water bodies. We also report on the characteristics of the capture-release strategy and provide our observations indicating no escalated aggression by leopards' post-release. In conclusion, this study evaluates widely employed approaches to conflict mitigation and suggests the continuous review and assessment of management strategies for mitigating human-leopard conflicts.

Identifying risk zones and landscape features that affect common leopard depredation.

Human-wildlife conflict (HWC) is a pressing issue worldwide but varies by species over time and place. One of the most prevalent forms of HWC in the mid-hills of Nepal is human-common-leopard conflict (HLC). Leopard attacks, especially in forested areas, can severely impact villagers and their livestock. Information on HLC in the Gorkha district was scarce, thus making it an ideal location to identify high-risk zones and landscape variables associated with such events. Registered cases were collected and reviewed from the Division Forest Office (DFO) during 2019-2021. Claims from DFO records were confirmed with herders and villagers via eight focus group discussions. To enhance modeling success, researchers identified a total of 163 leopard attack locations on livestock, ensuring a minimum distance of at least 100 meters between locations. Using maximum entropy (MaxEnt) and considering 13 environmental variables, we mapped common leopard attack risk zones. True Skill Statistics (TSS) and area under receiver-operator curve (AUC) were used to evaluate and validate the Output. Furthermore, 10 replications, 1,000 maximum iterations, and 1000 background points were employed during modeling. The average AUC value for the model, which was 0.726 ± 0.021, revealed good accuracy. The model performed well, as indicated by a TSS value of 0.61 ± 0.03. Of the total research area (27.92 km2), about 74% was designated as a low-risk area, 19% as a medium-risk area, and 7% as a high-risk area. Of the 13 environmental variables, distance to water (25.2%) was the most significant predictor of risk, followed by distance to road (16.2%) and elevation (10.7%). According to response curves, the risk of common leopard is highest in the areas between 1.5 to 2 km distances from the water sources, followed by the closest distance from a road and an elevation of 700 to 800 m. Results suggest that managers and local governments should employ intervention strategies immediately to safeguard rural livelihoods in high-risk areas. Improvements include better design of livestock corrals, insurance, and total compensation of livestock losses. Settlements near roads and water sources should improve the design and construction of pens and cages to prevent livestock loss. More studies on the characteristics of victims are suggested to enhance understanding of common leopard attacks, in addition to landscape variables. Such information can be helpful in formulating the best management practices.

Leopard density and determinants of space use in a farming landscape in South Africa.

Protected areas are traditionally the foundation of conservation strategy, but land not formally protected is of particular importance for the conservation of large carnivores because of their typically wide-ranging nature. In South Africa, leopard (Panthera pardus) population decreases are thought to be occurring in areas of human development and intense negative interactions, but research is biased towards protected areas, with quantitative information on population sizes and trends in non-protected areas severely lacking. Using Spatially Explicit Capture-Recapture and occupancy techniques including 10 environmental and anthropogenic covariates, we analysed camera trap data from commercial farmland in South Africa where negative human-wildlife interactions are reported to be high. Our findings demonstrate that leopards persist at a moderate density (2.21 /100 km2) and exhibit signs of avoidance from areas where lethal control measures are implemented. This suggests leopards have the potential to navigate mixed mosaic landscapes effectively, enhancing their chances of long-term survival and coexistence with humans. Mixed mosaics of agriculture that include crops, game and livestock farming should be encouraged and, providing lethal control is not ubiquitous in the landscape, chains of safer spaces should permit vital landscape connectivity. However, continuing to promote non-lethal mitigation techniques remains vital.

Whole-genome analysis of Escherichia coli isolated from wild Amur tiger (Panthera tigris altaica) and North China leopard (Panthera pardus japonensis).

Background: Escherichia coli is an important intestinal flora, of which pathogenic E. coli is capable of causing many enteric and extra-intestinal diseases. Antibiotics are essential for the treatment of bacterial infections caused by pathogenic E. coli; however, with the widespread use of antibiotics, drug resistance in E. coli has become particularly serious, posing a global threat to human, animal, and environmental health. While the drug resistance and pathogenicity of E. coli carried by tigers and leopards in captivity have been studied intensively in recent years, there is an extreme lack of information on E. coli in these top predators in the wild environment. Methods: Whole genome sequencing data of 32 E. coli strains collected from the feces of wild Amur tiger (Panthera tigris altaica, n = 24) and North China leopard (Panthera pardus japonensis, n = 8) were analyzed in this article. The multi-locus sequence types, serotypes, virulence and resistance genotypes, plasmid replicon types, and core genomic SNPs phylogeny of these isolates were studied. Additionally, antimicrobial susceptibility testing (AST) was performed on these E. coli isolates. Results: Among the E. coli isolates studied, 18 different sequence types were identified, with ST939 (21.9%), ST10 (15.6%), and ST3246 (9.4%) being the most prevalent. A total of 111 virulence genes were detected, averaging about 54 virulence genes per sample. They contribute to invasion, adherence, immune evasion, efflux pump, toxin, motility, stress adaption, and other virulence-related functions of E. coli. Sixty-eight AMR genes and point mutations were identified. Among the detected resistance genes, those belonging to the efflux pump family were the most abundant. Thirty-two E. coli isolates showed the highest rate of resistance to tetracycline (14/32; 43.8%), followed by imipenem (4/32; 12.5%), ciprofloxacin (3/32; 9.4%), doxycycline (2/32; 6.3%), and norfloxacin (1/32; 3.1%). Conclusions: Our results suggest that E. coli isolates carried by wild Amur tigers and North China leopards have potential pathogenicity and drug resistance.

Species discrimination from hair using ATR-FTIR spectroscopy: Application in wildlife forensics.

Hair is a commonly encountered trace evidence in wildlife crimes involving mammals and can be used for species identification which is essential for subsequent judicial proceedings. This proof of concept study aims, to distinguish the black guard hair of three wild cat species belonging to the genus Panthera i.e. Royal Bengal Tiger (Panthera tigris tigris), Indian Leopard (Panthera pardus fusca), and Snow Leopard (Panthera uncia) using a rapid and non-destructive ATR-FTIR spectroscopic technique in combination with chemometrics. A training dataset including 72 black guard hair samples of three species (24 samples from each species) was used to construct chemometric models. A PLS2-DA model successfully classified these three species into distinct classes with R-Square values of 0.9985 (calibration) and 0.8989 (validation). VIP score was also computed, and a new PLS2DA-V model was constructed using variables with a VIP score ≥ 1. External validation was performed using a validation dataset including 18 black guard hair samples (6 samples per species) to validate the constructed PLS2-DA model. It was observed that PLS2-DA model provides greater accuracy and precision compared to the PLS2DA-V model during cross-validation and external validation. The developed PLS2-DA model was also successful in differentiating human and non-human hair with R-Square values of 0.99 and 0.91 for calibration and validation, respectively. Apart from this, a blind test was also carried out using 10 unknown hair samples which were correctly classified into their respective classes providing 100 % accuracy. This study highlights the advantages of ATR-FTIR spectroscopy associated with PLS-DA for differentiation and identification of the Royal Bengal Tiger, Indian Leopard, and Snow Leopard hairs in a rapid, accurate, eco-friendly, and non-destructive way.

Turning the tide on big cat trade: Expert opinion on trends and conservation lessons from the Republic of Korea.

Unsustainable trade in big cats affects all species in the genus, Panthera, and is one of the foremost threats to their conservation. To provide further insight into the impact of policy interventions intended to address this issue, we examine the case study of the Republic of Korea (South Korea), which in the early 1990s was one of the world's largest importers of tiger (Panthera tigris) bone and a major manufacturer of tiger-derived medicinal products. In 1993, South Korea became a Party to the Convention on International Trade in Endangered Species (CITES) and introduced a ban on commercial trade in CITES Appendix I-listed big cats a year later. We used an expert-based questionnaire survey and an exploration of the CITES trade database to investigate what has since happened to big cat trade in South Korea. Expert opinion suggested that big cat trade has likely substantially reduced since the early 1990s, as a result of the trade ban and broad socioeconomic changes. However, illegal trade has not been eradicated entirely and we were able to confirm that products reportedly derived from big cats were still publicly available for sale on a range of Korean online marketplaces, sometimes openly. The items most commonly reported by respondents from post-1994 trade and supported by expert-led evidence were tiger and leopard (Panthera pardus) skins and tiger bone wine. Although South Korea may provide a useful case study of a historically significant consumer country for tiger which has made strong progress in addressing unsustainable levels of big cat trade within a short period of time, there remains a need to address recalcitrant small-scale, illegal trade. We also recommend further investigation regarding reports of South Korean nationals being involved in illegal trade in tiger-derived products in Southeast Asia.

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.

NONDOMESTIC FELID ABC BLOOD PHENOTYPING, GENOTYPING, AND CROSSMATCHING.

Based upon previous clinical experience with domestic cats (Felis catus), the ability to assess ABC blood types and blood (in-)compatibilities of nondomestic felids, and adequately consider and plan for blood transfusions, may be important. Although nondomestic felids appear to have an ABC blood group system similar to domestic cats, typing with point-of-care kits and by CMAH genotyping for domestic cats have not been reported. In this study, 162 blood samples from 18 different nondomestic felid species (cheetah [Acinonyx jubatus, n = 42], lion [Panthera leo, n = 33], tiger [Panthera tigris, n = 23], Canada lynx [Lynx canadensis, n = 11], snow leopard [Uncia uncia, n = 10], puma [Puma concolor, n = 7], clouded leopard [Neofelis nebulosa, n = 6], serval [Leptailurus serval, n = 5], jaguar [Panthera onca, n = 5], fishing cat [Prionailurus viverrinus, n = 4], Pallas cat [Felis manul, n = 3], bobcat [Lynx rufus, n = 3], ocelot [Leopardus pardalis, n = 3], black footed cat [Felis nigripes, n = 2], leopard [Panthera pardus, n = 2], African wildcat [Felis lybica, n = 1], caracal [Caracal caracal, n = 1], and sand cat [Felis margarita, n = 1]) were ABC blood typed by laboratory and point-of-care tests, genotyped for four known CMAH variants for type B and type C (AB) phenotypes, and crossmatched with one another and domestic type A cats. Traditional tube typing identified blood type A (n = 106), type B (n = 8), type C (n = 43), and no discernible ABC type (n = 4). Several discrepancies were found between point-of-care and traditional typing test results. None of the tested felids possessed the four CMAH variants responsible for type B and C (AB) in domestic cats. Crossmatch incompatibilities (≥2+ agglutination) were identified within and between nondomestic felid species and beyond ABC incompatibilities. Of 26 crossmatches performed between domestic cats and various nondomestic felids, only 7 (27%) were compatible. In conclusion, point-of-care typing kits and CMAH genotyping, successfully used in domestic cats, may not identify the correct ABC blood type in nondomestic felids. Prior crossmatching is recommended to increase the likelihood of compatible transfusions between any nondomestic felids.

Preserving the spots: Jaguar (Panthera onca) distribution and priority conservation areas in Colombia.

The jaguar (Panthera onca) is a charismatic species considered Vulnerable in Colombia but yet largely unknown in the country. The species is mostly threatened by the continuous decline in its habitats, mostly derived from deforestation and habitat loss, additional to hunting and conflicts with humans. Thus, the future of jaguars in Colombia depends on protecting and recovering existing habitats. The aims of this study were to 1) evaluate jaguar distribution and identify the remnant patches of habitat in Colombia, 2) define an ecological connectivity network within the country, and 3) propose a priority areas portfolio for the conservation and recovery of jaguars. We used a presence background model for estimating species potential distribution and subsequently identified remaining habitat patches across the country based on land cover and species-specific ecological attributes. We then created an ecological connectivity network based on circuit theory and following a multi-criteria approach identified jaguar priority areas for conservation (JPCA) and recovery (JPRA). Jaguar potential distribution comprises 1103122.43 km2, from which 56.71% maintain suitable patches of potential habitat. We identified 960 corridors between remnant patches of natural or semi-natural vegetation. Based on the criteria, JPCAs with greater importance were identified in each of the five Colombian regions. JPRAs were located mainly towards the Andean and Caribbean regions. These JPCAs and JPRAs could serve as a guide for designing and implementing management strategies for the long-term conservation and recovery of the species in Colombia.

NEW RECORDS OF HELMINTHS OF THE JAGUAR IN MEXICO, WITH AN UPDATED LIST OF SPECIES IN THE AMERICAS.

An inventory of parasites infecting the jaguar (Panthera onca) across its distribution range is relevant for the conservation of this threatened big cat. In this study, we report the occurrence of helminths in a jaguar from Mexico using morphological techniques (cleared and stained mounts and scanning electron microscopy) and partial sequences of the 28S ribosomal RNA (28S rRNA) gene and the cytochrome c oxidase 1 mitochondrial (COI) gene. We also provide an updated list of helminth species reported in jaguars in the Americas. Three helminth taxa are identified in the jaguar examined from Mexico: Toxocara cati, Physaloptera sp., and Taenia sp. The new 28S rRNA sequences of To. cati, Physaloptera sp., and Taenia sp. and the COI sequence of Taenia sp. corroborate the identity of the helminths isolated from this host. One hundred and twenty-nine records of helminths parasitizing jaguars from 49 studies up to May 2023 were identified in the Americas. In most of these studies (73.6%), helminths were identified using coproparasitological techniques. Sixteen helminths (7 nematodes, 5 cestodes, 3 acanthocephalans, and 1 trematode) were identified at the species level in free-ranging and captive jaguars. The study demonstrates the value of an integrative taxonomy approach to increase the accuracy of parasite identification in wildlife, especially when helminth specimens are scarce or poorly fixed.

Dietary type (carnivore, herbivore and omnivore) and animal species modulate the nutritional metabolome of terrestrial species.

Ecometabolomics could be implemented as a powerful tool in molecular ecology studies, but it is necessary to know the baseline of certain metabolites and understand how different traits could affect the metabolome of the animals. Therefore, the main objective of this study was to provide values for the nutritional metabolome profile of different diet groups and animal species, as well as to study the differences in the metabolomic profile due to the effect of diet type and species. To achieve this goal, blood samples were taken from healthy animals (n = 43) of different species: lion (Panthera leo), jaguar (Panthera onca), chimpanzee (Pan troglodytes), bison (Bison bison), gazelle (Gazella cuvieri) and fallow deer (Dama dama), and with different types of diet (carnivore, herbivore and omnivore). Each blood sample was analysed to determine nutritional metabolites. The main results this study provides are the nutritional metabolic profile of these animals based on the type of diet and the animal species. A significant effect of the dietary type was found on nutritional metabolite levels, with those metabolites related to protein metabolism (total protein and creatine) being higher in carnivores. There is also an effect of the species on nutritional metabolites, observing a metabolome differentiation between lion and jaguar. In the case of herbivores, bison showed higher levels of uric acid and cholesterol, and lower urea levels than gazelle and fallow deer. More molecular ecology studies are needed to further the knowledge of the metabolism of these animals.

Ancestry dynamics and trait selection in a designer cat breed.

The Bengal cat breed was developed from intercrosses between the Asian leopard cat, Prionailurus bengalensis, and the domestic cat, Felis catus, with a last common ancestor approximately 6 million years ago. Predicted to derive ∼94% of their genome from domestic cats, regions of the leopard cat genome are thought to account for the unique pelage traits and ornate color patterns of the Bengal breed, which are similar to those of ocelots and jaguars. We explore ancestry distribution and selection signatures in the Bengal breed by using reduced representation and whole-genome sequencing from 947 cats. The mean proportion of leopard cat DNA in the Bengal breed is 3.48%, lower than predicted from breed history, and is broadly distributed, covering 93% of the Bengal genome. Overall, leopard cat introgressions do not show strong signatures of selection across the Bengal breed. However, two popular color traits in Bengal cats, charcoal and pheomelanin intensity, are explained by selection of leopard cat genes whose expression is reduced in a domestic cat background, consistent with genetic incompatibility resulting from hybridization. We characterize several selective sweeps in the Bengal genome that harbor candidate genes for pelage and color pattern and that are associated with domestic, rather than leopard, cat haplotypes. We identify the molecular and phenotypic basis of one selective sweep as reduced expression of the Fgfr2 gene, which underlies glitter, a trait desired by breeders that affects hair texture and light reflectivity.

Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus).

The anthropogenic impact on wildlife is ever increasing. With shrinking habitats, wild populations are being pushed to co-exist in proximity to humans leading to an increased threat of infectious diseases. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immune system (IIS) is the body's first line of defense against pathogens. High variability in IIS genes, like toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in vulnerable species for conservation. Large predators are threatened globally including leopards and cheetahs, both listed as 'vulnerable' by IUCN. To examine IIS diversity in these sympatric species, we used next-generation-sequencing to compare selected TLR genes in African leopards and cheetahs. Despite differences, both species show some TLR haplotype similarity. Historic cheetahs from all subspecies exhibit greater genetic diversity than modern Southern African cheetahs. The diversity in investigated TLR genes is lower in modern Southern African cheetahs than in African leopards. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs. However, this may not yet impact the health of this cheetah subspecies.

Positive selection in gamete interaction proteins in Carnivora.

The absence of robust interspecific isolation barriers among pantherines, including the iconic South American jaguar (Panthera onca), led us to study molecular evolution of typically rapidly evolving reproductive proteins within this subfamily and related groups. In this study, we delved into the evolutionary forces acting on the zona pellucida (ZP) gamete interaction protein family and the sperm-oocyte fusion protein pair IZUMO1-JUNO across the Carnivora order, distinguishing between Caniformia and Feliformia suborders and anticipating few significant diversifying changes in the Pantherinae subfamily. A chromosome-resolved jaguar genome assembly facilitated coding sequences, enabling the reconstruction of protein evolutionary histories. Examining sequence variability across more than 30 Carnivora species revealed that Feliformia exhibited significantly lower diversity compared to its sister taxa, Caniformia. Molecular evolution analyses of ZP2 and ZP3, subunits directly involved in sperm-recognition, unveiled diversifying positive selection in Feliformia, Caniformia and Pantherinae, although no significant changes were linked to sperm binding. Structural cross-linking ZP subunits, ZP4 and ZP1 exhibited lower levels or complete absence of positive selection. Notably, the fusion protein IZUMO1 displayed prominent positive selection signatures and sites in basal lineages of both Caniformia and Feliformia, extending along the Caniformia subtree but absent in Pantherinae. Conversely, JUNO did not exhibit any positive selection signatures across tested lineages and clades. Eight Caniformia-specific positive selected sites in IZUMO1 were detected within two JUNO-interaction clusters. Our findings provide for the first time insights into the evolutionary trajectories of ZP proteins and the IZUMO1-JUNO gamete interaction pair within the Carnivora order.