The late Miocene radiation of modern Felidae: a genetic assessment

Evolutionary insights into Felidae iris color through ancestral state reconstruction

Few studies have explored eye (iris) color evolution beyond humans and domesticated animals. Felids exhibit significant eye color diversity, unlike their brown-eyed relatives, making them an ideal model to study the evolution of eye color in natural populations. Through machine learning analysis of public photographs, five felid eye colors were identified: brown, green, yellow, gray, and blue. The presence or absence of these colors was reconstructed on a phylogeny, as well as their specific quantitative shades. The ancestral felid population likely had brown-eyed and gray-eyed individuals, the latter color being pivotal for the diversification of eye color seen in modern felids. Additionally, yellow eyes are highly associated with and may be necessary for, the evolution of round pupils in felids. These findings enhance the understanding of eye color evolution, and the methods presented in this work are widely applicable and will facilitate future research into the phylogenetic reconstruction of color beyond irises.

Spermatogonial stem cell technologies: applications from human medicine to wildlife conservation†

Spermatogonial stem cell (SSC) technologies that are currently under clinical development to reverse human infertility hold the potential to be adapted and applied for the conservation of endangered and vulnerable wildlife species. The biobanking of testis tissue containing SSCs from wildlife species, aligned with that occurring in pediatric human patients, could facilitate strategies to improve the genetic diversity and fitness of endangered populations. Approaches to utilize these SSCs could include spermatogonial transplantation or testis tissue grafting into a donor animal of the same or a closely related species, or in vitro spermatogenesis paired with assisted reproduction approaches. The primary roadblock to progress in this field is a lack of fundamental knowledge of SSC biology in non-model species. Herein, we review the current understanding of molecular mechanisms controlling SSC function in laboratory rodents and humans, and given our particular interest in the conservation of Australian marsupials, use a subset of these species as a case-study to demonstrate gaps-in-knowledge that are common to wildlife. Additionally, we review progress in the development and application of SSC technologies in fertility clinics and consider the translation potential of these techniques for species conservation pipelines.

Feline eye-inspired artificial vision for enhanced camouflage breaking under diverse light conditions

Biologically inspired artificial vision research has led to innovative robotic vision systems with low optical aberration, wide field of view, and compact form factor. However, challenges persist in object detection and recognition against complex backgrounds and varied lighting. Inspired by the feline eye, which features a vertically elongated pupil and tapetum lucidum, this study introduces an artificial vision system designed for superior object detection and recognition in a monocular framework. Using a slit-like elliptical aperture and a patterned metal reflector beneath a hemispherical silicon photodiode array, the system reduces excessive light and enhances photosensitivity. This design achieves clear focus under bright light and enhanced sensitivity in dim conditions. Theoretical and experimental analyses demonstrate the system's ability to filter redundant information and detect camouflaged objects in diverse lighting, representing a substantial advancement in monocular camera technology and the potential of biomimicry in optical innovations.

Mitogenomic analysis of a late Pleistocene jaguar from North America

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, United States, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

Toward telomere-to-telomere cat genomes for precision medicine and conservation biology

Genomic data from species of the cat family Felidae promise to stimulate veterinary and human medical advances, and clarify the coherence of genome organization. We describe how interspecies hybrids have been instrumental in the genetic analysis of cats, from the first genetic maps to propelling cat genomes toward the T2T standard set by the human genome project. Genotype-to-phenotype mapping in cat models has revealed dozens of health-related genetic variants, the molecular basis for mammalian pigmentation and patterning, and species-specific adaptations. Improved genomic surveillance of natural and captive populations across the cat family tree will increase our understanding of the genetic architecture of traits, population dynamics, and guide a future of genome-enabled biodiversity conservation.

Evolutionary patterns of cat-like carnivorans unveil drivers of the sabertooth morphology

The sabertooth morphology stands as a classic case of convergence, manifesting recurrently across various vertebrate groups, prominently within two carnivorans clades: felids and nimravids. Nonetheless, the evolutionary mechanisms driving these recurring phenotypes remain insufficiently understood, lacking a robust phylogenetic and spatiotemporal framework. We reconstruct the tempo and mode of craniomandibular evolution of Felidae and Nimravidae and evaluate the strength of the dichotomy between conical and saber-toothed species, as well as within saber-toothed morphotypes. To do so, we investigate morphological variation, convergence, phenotypic integration, and evolutionary rates, employing a comprehensive dataset of nearly 200 3D models encompassing mandibles and crania from both extinct and extant feline-like carnivorans, spanning their entire evolutionary timeline. Our results reject the hypothesis of a distinctive sabertooth morphology, revealing instead a continuous spectrum of feline-like phenotypes in both the cranium and mandible, with sporadic instances of unequivocal convergence. Disparity peaked at the end of the Miocene and is usually higher in clades containing taxa with extreme sabertoothed adaptations. We show that taxa with saberteeth exhibit a lower degree of craniomandibular integration, allowing to exhibit a greater range of phenotypes. Those same groups usually show a burst of morphological evolutionary rate at the beginning of their evolutionary history. Consequently, we propose that a reduced degree of integration coupled with rapid evolutionary rates emerge as key components in the development of a sabertooth morphology in multiple clades.

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.

Is there a loophole in Dollo's law? A DevoEvo perspective on irreversibility (of felid dentition)

There is a longstanding interest in whether the loss of complex characters is reversible (so-called "Dollo's law"). Reevolution has been suggested for numerous traits but among the first was Kurtén, who proposed that the presence of the second lower molar (M2 ) of the Eurasian lynx (Lynx lynx) was a violation of Dollo's law because all other Felids lack M2 . While an early and often cited example for the reevolution of a complex trait, Kurtén and Werdelin used an ad hoc parsimony argument to support their. Here I revisit the evidence that M2 reevolved lynx using explicit parsimony and maximum likelihood models of character evolution and find strong evidence that Kurtén and Werdelin were correct-M2 reevolved in E. lynx. Next, I explore the developmental mechanisms which may explain this violation of Dollo's law and suggest that the reevolution of lost complex traits may arise from the reevolution of cis-regulatory elements and protein-protein interactions, which have a longer half-life after silencing that protein coding genes. Finally, I present a developmental model to explain the reevolution M2 in E. lynx, which suggest that the developmental programs required for the establishment of serially homologous characters may never really be lost so long as a single instance of the character remains-thus the gain and loss and regain of serially homologous characters, such mammalian molars, may be developmentally and evolutionarily "simple."

Extensive Phylogenomic Discordance and the Complex Evolutionary History of the Neotropical Cat Genus Leopardus

Even in the genomics era, the phylogeny of Neotropical small felids comprised in the genus Leopardus remains contentious. We used whole-genome resequencing data to construct a time-calibrated consensus phylogeny of this group, quantify phylogenomic discordance, test for interspecies introgression, and assess patterns of genetic diversity and demographic history. We infer that the Leopardus radiation started in the Early Pliocene as an initial speciation burst, followed by another in its subgenus Oncifelis during the Early Pleistocene. Our findings challenge the long-held notion that ocelot (Leopardus pardalis) and margay (L. wiedii) are sister species and instead indicate that margay is most closely related to the enigmatic Andean cat (L. jacobita), whose whole-genome data are reported here for the first time. In addition, we found that the newly sampled Andean tiger cat (L. tigrinus pardinoides) population from Colombia associates closely with Central American tiger cats (L. tigrinus oncilla). Genealogical discordance was largely attributable to incomplete lineage sorting, yet was augmented by strong gene flow between ocelot and the ancestral branch of Oncifelis, as well as between Geoffroy's cat (L. geoffroyi) and southern tiger cat (L. guttulus). Contrasting demographic trajectories have led to disparate levels of current genomic diversity, with a nearly tenfold difference in heterozygosity between Andean cat and ocelot, spanning the entire range of variability found in extant felids. Our analyses improved our understanding of the speciation history and diversity patterns in this felid radiation, and highlight the benefits to phylogenomic inference of embracing the many heterogeneous signals scattered across the genome.

Reproductive evaluation of Luisa, the last jaguar of the Caatinga

The in situ population of jaguars in the Caatinga is less than 250 individuals, subdivided into five subpopulations, and is classified as endangered regarding its risk of extinction. Luisa, a 15-year-old female weighing 36 kg, was the last known ex situ jaguar from this biome. Her reproductive evaluation is detailed in this manuscript. Luisa was subjected to both a clinical and laparoscopic evaluation of her reproductive system. After 45 days of reproductive investigation, she died unexpectedly, and skin fragments were taken to establish the postmortem fibroblast lineage. At the clinical evaluation, Luisa had small, undeveloped mammary gland and a small vulva, characteristic of a nulliparous female, with no mammary gland nodules, edema, or abnormal masses. By laparoscopy, normal-appearing bladder and bowel loops were observed, as were uterine horns with standard color, shape, and length with no striae. Ovaries and uterine horns seem free of fibrinous adhesions. Both ovaries showed a yellowish color, a fibrous consistency, a decreased size (atrophied), and no follicles, hemorrhagic corpus, corpus luteum, luteal scars, or other abnormal structures. We may assume that this jaguar female was infertile based on Luisa's mature age and the absence of birthing or ovarian activity signs. The harsh conditions of the Caatinga biome, which included low food availability and frequent conflicts with humans, may have impacted both the pregnancy and lactation of Luisa's mother and her development after birth.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

Ancient DNA reveals genetic admixture in China during tiger evolution

The tiger (Panthera tigris) is a charismatic megafauna species that originated and diversified in Asia and probably experienced population contraction and expansion during the Pleistocene, resulting in low genetic diversity of modern tigers. However, little is known about patterns of genomic diversity in ancient populations. Here we generated whole-genome sequences from ancient or historical (100-10,000 yr old) specimens collected across mainland Asia, including a 10,600-yr-old Russian Far East specimen (RUSA21, 8× coverage) plus six ancient mitogenomes, 14 South China tigers (0.1-12×) and three Caspian tigers (4-8×). Admixture analysis showed that RUSA21 clustered within modern Northeast Asian phylogroups and partially derived from an extinct Late Pleistocene lineage. While some of the 8,000-10,000-yr-old Russian Far East mitogenomes are basal to all tigers, one 2,000-yr-old specimen resembles present Amur tigers. Phylogenomic analyses suggested that the Caspian tiger probably dispersed from an ancestral Northeast Asian population and experienced gene flow from southern Bengal tigers. Lastly, genome-wide monophyly supported the South China tiger as a distinct subspecies, albeit with mitochondrial paraphyly, hence resolving its longstanding taxonomic controversy. The distribution of mitochondrial haplogroups corroborated by biogeographical modelling suggested that Southwest China was a Late Pleistocene refugium for a relic basal lineage. As suitable habitat returned, admixture between divergent lineages of South China tigers took place in Eastern China, promoting the evolution of other northern subspecies. Altogether, our analysis of ancient genomes sheds light on the evolutionary history of tigers and supports the existence of nine modern subspecies.

Updated list of the mammals of Costa Rica, with notes on recent taxonomic changes

Although Costa Rica occupies a mere 0.03% of the Earths land area, it nevertheless has recorded within its borders approximately 5% of the global diversity of mammals, thus making it one of the worlds megadiverse countries. Over the past ten years, 22 species have been added to the countrys inventory, bringing the total number known as here documented to 271; Chiroptera account for ten of these, having grown to 124 from 114; rodents have increased by eight species, from 47 to 55, with the caveat that we include three invasive species of Muridae that have gone feral. In contrast, the number of orders has decreased by one, by Artiodactyla incorporating the former Cetacea. Notes are provided for all taxonomic novelties since the last update. Since the first taxonomic compendium of the mammals of Costa Rica in 1869, the number of known species has grown by approximately 1.22 species year-1 (R2 = 0.96). Since 1983 however, this growth rate has been 1.64 species year-1 (R2 = 0.98). Despite this strong growth, an asymptote in the number of known species has not been reached. Conservation remains a primary need: over 60% of the countrys mammal species show population trends that are decreasing (13%), unknown (37%), or not assessed (11%), based on IUCN criteria. These analyses suggest that much remains to be known regarding the number of mammal species living in Costa Rica, but also that much more remains to be done to safeguard Costa Ricas exceptional biodiversity heritage.

How genomic insights into the evolutionary history of clouded leopards inform their conservation

Clouded leopards (Neofelis spp.), a morphologically and ecologically distinct lineage of big cats, are severely threatened by habitat loss and fragmentation, targeted hunting, and other human activities. The long-held poor understanding of their genetics and evolution has undermined the effectiveness of conservation actions. Here, we report a comprehensive investigation of the whole genomes, population genetics, and adaptive evolution of Neofelis. Our results indicate the genus Neofelis arose during the Pleistocene, coinciding with glacial-induced climate changes to the distributions of savannas and rainforests, and signatures of natural selection associated with genes functioning in tooth, pigmentation, and tail development, associated with clouded leopards' unique adaptations. Our study highlights high-altitude adaptation as the main factor driving nontaxonomic population differentiation in Neofelis nebulosa. Population declines and inbreeding have led to reduced genetic diversity and the accumulation of deleterious variation that likely affect reproduction of clouded leopards, highlighting the urgent need for effective conservation efforts.

Bioacoustics supports genus identification in piranhasa)

In different teleost species, sound production can utilize specific coding schemes to avoid confusion between species during communication. Piranhas are vocal Neotropical fishes, and both Pygocentrus and Serrasalmus produce similar pulsed sounds using the same sound-producing mechanism. In this study, we analysed the sounds of three Pygocentrus and nine Serrasalmus species to determine whether sounds can be used to discriminate piranha species at both the species and genus levels. Our analysis of temporal and frequency data supports the idea that the sounds of Serrasalmus and Pygocentrus species are species specific, and that different acoustic features can be used to differentiate taxa at the genus level. Specifically, the sounds of Serrasalmus species are shorter, louder, and have a shorter pulse period (as determined after correction for standard length). This suggests that sounds can be used to support taxonomy at the genus level as well as the species level.

"Spotting" Mycobacterium bovis infection in leopards (Panthera pardus) - novel application of diagnostic tools

Background

Mycobacterium bovis (M. bovis) is the causative agent of animal tuberculosis (TB) which poses a threat to many of South Africa's most iconic wildlife species, including leopards (Panthera pardus). Due to limited tests for wildlife, the development of accurate ante-mortem tests for TB diagnosis in African big cat populations is urgently required. The aim of this study was to evaluate currently available immunological assays for their ability to detect M. bovis infection in leopards.

Methods

Leopard whole blood (n=19) was stimulated using the QuantiFERON Gold Plus In-Tube System (QFT) to evaluate cytokine gene expression and protein production, along with serological assays. The GeneXpert® MTB/RIF Ultra (GXU®) qPCR assay, mycobacterial culture, and speciation by genomic regions of difference PCR, was used to confirm M. bovis infection in leopards.

Results

Mycobacterium bovis infection was confirmed in six leopards and individuals that were tuberculin skin test (TST) negative were used for comparison. The GXU® assay was positive using all available tissue homogenates (n=5) from M. bovis culture positive animals. Mycobacterium bovis culture-confirmed leopards had greater antigen-specific responses, in the QFT interferon gamma release assay, CXCL9 and CXCL10 gene expression assays, compared to TST-negative individuals. One M. bovis culture-confirmed leopard had detectable antibodies using the DPP® Vet TB assay.

Conclusion

Preliminary results demonstrated that immunoassays and TST may be potential tools to identify M. bovis-infected leopards. The GXU® assay provided rapid direct detection of infected leopards. Further studies should aim to improve TB diagnosis in wild felids, which will facilitate disease surveillance and screening.

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.

Morphological and Genetics Support for a Hitherto Undescribed Spotted Cat Species (Genus Leopardus; Felidae, Carnivora) from the Southern Colombian Andes

In 1989, a skin of a small spotted cat, from the Galeras Volcano in southern Colombia (Nariño Department), was donated to the Instituto Alexander von Humboldt (identification, ID 5857) at Villa de Leyva (Boyacá Department, Colombia). Although originally classified as Leopardus tigrinus, its distinctiveness merits a new taxonomic designation. The skin is distinct from all known L. tigrinus holotypes as well as from other Leopardus species. Analysis of the complete mitochondrial genomes from 44 felid specimens (including 18 L. tigrinus and all the current known species of the genus Leopardus), the mtND5 gene from 84 felid specimens (including 30 L. tigrinus and all the species of the genus Leopardus), and six nuclear DNA microsatellites (113 felid specimens of all the current known species of the genus Leopardus) indicate that this specimen does not belong to any previously recognized Leopardus taxon. The mtND5 gene suggests this new lineage (the Nariño cat as we name it) is a sister taxon of Leopardus colocola. The mitogenomic and nuclear DNA microsatellite analyses suggest that this new lineage is the sister taxon to a clade formed by Central American and trans-Andean L. tigrinus + (Leopardus geoffroyi + Leopardus guigna). The temporal split between the ancestor of this new possible species and the most recent ancestor within Leopardus was dated to 1.2-1.9 million years ago. We consider that this new unique lineage is a new species, and we propose the scientific name Leopardus narinensis.

Chromosome-level, nanopore-only genome and allele-specific DNA methylation of Pallas's cat, Otocolobus manul

Pallas's cat, or the manul cat (Otocolobus manul), is a small felid native to the grasslands and steppes of central Asia. Population strongholds in Mongolia and China face growing challenges from climate change, habitat fragmentation, poaching, and other sources. These threats, combined with O. manul's zoo collection popularity and value in evolutionary biology, necessitate improvement of species genomic resources. We used standalone nanopore sequencing to assemble a 2.5 Gb, 61-contig nuclear assembly and 17097 bp mitogenome for O. manul. The primary nuclear assembly had 56× sequencing coverage, a contig N50 of 118 Mb, and a 94.7% BUSCO completeness score for Carnivora-specific genes. High genome collinearity within Felidae permitted alignment-based scaffolding onto the fishing cat (Prionailurus viverrinus) reference genome. Manul contigs spanned all 19 felid chromosomes with an inferred total gap length of less than 400 kilobases. Modified basecalling and variant phasing produced an alternate pseudohaplotype assembly and allele-specific DNA methylation calls; 61 differentially methylated regions were identified between haplotypes. Nearest features included classical imprinted genes, non-coding RNAs, and putative novel imprinted loci. The assembled mitogenome successfully resolved existing discordance between Felinae nuclear and mtDNA phylogenies. All assembly drafts were generated from 158 Gb of sequence using seven minION flow cells.

cDNA Cloning of Feline PIWIL1 and Evaluation of Expression in the Testis of the Domestic Cat

The PIWI clade of Argonaute proteins is essential for spermatogenesis in all species examined to date. This protein family binds specific classes of small non-coding RNAs known as PIWI-interacting RNAs (piRNAs) which together form piRNA-induced silencing complexes (piRISCs) that are recruited to specific RNA targets through sequence complementarity. These complexes facilitate gene silencing through endonuclease activity and guided recruitment of epigenetic silencing factors. PIWI proteins and piRNAs have been found to play multiple roles in the testis including the maintenance of genomic integrity through transposon silencing and facilitating the turnover of coding RNAs during spermatogenesis. In the present study, we report the first characterization of PIWIL1 in the male domestic cat, a mammalian system predicted to express four PIWI family members. Multiple transcript variants of PIWIL1 were cloned from feline testes cDNA. One isoform shows high homology to PIWIL1 from other mammals, however, the other has characteristics of a "slicer null" isoform, lacking the domain required for endonuclease activity. Expression of PIWIL1 in the male cat appears limited to the testis and correlates with sexual maturity. RNA-immunoprecipitation revealed that feline PIWIL1 binds small RNAs with an average size of 29 nt. Together, these data suggest that the domestic cat has two PIWIL1 isoforms expressed in the mature testis, at least one of which interacts with piRNAs.

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

Narrowly defined taxa on a global scale: The phylogeny and taxonomy of the genera Catriona and Tenellia (Nudibranchia, Trinchesiidae) favours fine-scale taxonomic differentiation and dissolution of the "lumpers & splitters" dilemma

By applying morphological and molecular data on two genera of the nudibranch molluscs it is shown that the tension between taxonomic practice and evolutionary processes persists. A review of the related genera Catriona and Tenellia is used to demonstrate that the fine-scale taxonomic differentiation is an important tool in the integration of morphological and molecular data. This is highlighted by the hidden species problem and provides strong argument that the genus must be kept as a maximally narrowly-defined entity. Otherwise, we are forced to compare a highly disparate species under the putatively lumped name "Tenellia". We demonstrate this in the present study by applying a suite of delimitation methods and describing a new species of Tenellia from the Baltic Sea. The new species possesses fine-scale morphological distinguishing features, which were not investigated before. The true, narrowly defined genus Tenellia represents a peculiar taxon with a clearly expressed paedomorphic characters and predominantly brackish-water habitats. The phylogenetically related genus Catriona, of which three new species are described here, clearly demonstrates different features. A lumping decision to call many morphologically and evolutionary different taxa as "Tenellia" will downgrade the taxonomic and phylogenetic resolution of the entire family Trinchesiidae to just a single genus. The dissolution of the dilemma of "lumpers & splitters", which still significantly affects taxonomy, will further help to make systematics a true evolutionary discipline.

Confidence intervals in molecular dating by maximum likelihood

Molecular dating has been widely used to infer the times of past evolutionary events using molecular sequences. This paper describes three bootstrap methods to infer confidence intervals under a penalized likelihood framework. The basic idea is to use data pseudoreplicates to infer uncertainty in the branch lengths of a phylogeny reconstructed with molecular sequences. The three specific bootstrap methods are nonparametric (direct tree bootstrapping), semiparametric (rate smoothing), and parametric (Poisson simulation). Our extensive simulation study showed that the three methods perform generally well under a simple strict clock model of molecular evolution; however, the results were less positive with data simulated using an uncorrelated or a correlated relaxed clock model. Several factors impacted, possibly in interaction, the performance of the confidence intervals. Increasing the number of calibration points had a positive effect, as well as increasing the sequence length or the number of sequences although both latter effects depended on the model of evolution. A case study is presented with a molecular phylogeny of the Felidae (Mammalia: Carnivora). A comparison was made with a Bayesian analysis: the results were very close in terms of confidence intervals and there was no marked tendency for an approach to produce younger or older bounds compared to the other.

Whole genome analysis of clouded leopard species reveals an ancient divergence and distinct demographic histories

Similar to other apex predator species, populations of mainland (Neofelis nebulosa) and Sunda (Neofelis diardi) clouded leopards are declining. Understanding their patterns of genetic variation can provide critical insights on past genetic erosion and a baseline for understanding their long-term conservation needs. As a step toward this goal, we present draft genome assemblies for the two clouded leopard species to quantify their phylogenetic divergence, genome-wide diversity, and historical population trends. We estimate that the two species diverged 5.1 Mya, much earlier than previous estimates of 1.41 Mya and 2.86 Mya, suggesting they separated when Sundaland was becoming increasingly isolated from mainland Southeast Asia. The Sunda clouded leopard displays a distinct and reduced effective population size trajectory, consistent with a lower genome-wide heterozygosity and SNP density, relative to the mainland clouded leopard. Our results provide new insights into the evolutionary history and genetic health of this unique lineage of felids.

Genome report: chromosome-level draft assemblies of the snow leopard, African leopard, and tiger (Panthera uncia, Panthera pardus pardus, and Panthera tigris)

The big cats (genus Panthera) represent some of the most popular and charismatic species on the planet. Although some reference genomes are available for this clade, few are at the chromosome level, inhibiting high-resolution genomic studies. We assembled genomes from 3 members of the genus, the tiger (Panthera tigris), the snow leopard (Panthera uncia), and the African leopard (Panthera pardus pardus), at chromosome or near-chromosome level. We used a combination of short- and long-read technologies, as well as proximity ligation data from Hi-C technology, to achieve high continuity and contiguity for each individual. We hope that these genomes will aid in further evolutionary and conservation research of this iconic group of mammals.

Genetics of randomly bred cats support the cradle of cat domestication being in the Near East

Cat domestication likely initiated as a symbiotic relationship between wildcats (Felis silvestris subspecies) and the peoples of developing agrarian societies in the Fertile Crescent. As humans transitioned from hunter-gatherers to farmers ~12,000 years ago, bold wildcats likely capitalized on increased prey density (i.e., rodents). Humans benefited from the cats' predation on these vermin. To refine the site(s) of cat domestication, over 1000 random-bred cats of primarily Eurasian descent were genotyped for single-nucleotide variants and short tandem repeats. The overall cat population structure suggested a single worldwide population with significant isolation by the distance of peripheral subpopulations. The cat population heterozygosity decreased as genetic distance from the proposed cat progenitor's (F.s. lybica) natural habitat increased. Domestic cat origins are focused in the eastern Mediterranean Basin, spreading to nearby islands, and southernly via the Levantine coast into the Nile Valley. Cat population diversity supports the migration patterns of humans and other symbiotic species.

A Reference Genome Assembly of the Bobcat, Lynx rufus

The bobcat (Lynx rufus) is a medium-sized carnivore well adapted to various environments and an indicator species for landscape connectivity. It is one of the 4 species within the extant Lynx genus in the family Felidae. Because of its broad geographic distribution and central role in food webs, the bobcat is important for conservation. Here we present a high-quality de novo genome assembly of a male bobcat located in Mendocino County, CA, as part of the California Conservation Genomics Project (CCGP). The assembly was generated using the standard CCGP pipeline from a combination of Omni-C and HiFi technologies. The primary assembly comprises 76 scaffolds spanning 2.4 Gb, represented by a scaffold N50 of 142 Mb, a contig N50 of 66.2 Mb, and a BUSCO completeness score of 95.90%. The bobcat genome will be an important resource for the effective management and conservation of this species and comparative genomics exploration.

Role of feeding specialization in taste receptor loss: insights from sweet and umami receptor evolution in Carnivora

Controversy and misunderstanding surround the role of feeding specialization in taste receptor loss in vertebrates. We refined and tested the hypothesis that this loss is caused by feeding specializations. Specifically, feeding specializations were proposed to trigger time-dependent process of taste receptor loss through deprivation of benefit of using the receptor's gustatory function. We propose that this process may be accelerated by abiotic environmental conditions or decelerated/stopped because of extragustatory functions of the receptor's protein(s). As test case we used evolution of the sweet (TAS1R2+TAS1R3) and umami (TAS1R1+TAS1R3) receptors in Carnivora (dogs, cats, and kin). We predicted these receptors' absence/presence using data on presence/absence of inactivating mutations in these receptors' genes and data from behavioral sweet/umami preference tests. We identified 20 evolutionary events of sweet (11) or umami (9) receptor loss. These events affected species with feeding specializations predicted to favor sweet/umami receptor loss (27 and 22 species, respectively). All species with feeding habits predicted to favor sweet/umami receptor retention (11 and 24, respectively) were found to retain that receptor. Six species retained the sweet (5) or umami (1) receptor despite feeding specialization predicted to favor loss of that receptor, which can be explained by the time dependence of sweet/umami receptor loss process and the possible decelerating effect of TAS1R extragustatory functions so that the sweet/umami receptor process is ongoing in these species. Our findings support the idea that feeding specialization leads to taste receptor loss and is the main if not only triggering factor for evolutionary loss of taste receptors.

Collaborative behaviour and coalitions in male jaguars (Panthera onca)—evidence and comparison with other felids

Abstract

Most large felids are classified as solitary species, with only lions (Panthera leo) and cheetahs (Acinonyx jubatus) exhibiting social, collaborative behaviours. Herein, we present evidence of the formation of male coalitions by jaguars (Panthera onca), based on data from five studies conducted with camera trapping, GPS telemetry, and direct observations in the Venezuelan Llanos and Brazilian Pantanal. Out of 7062 male records obtained with camera traps or visual observations, we detected 105 cases of male-male interactions, of which we classified 18 as aggression, nine as tolerance, 70 as cooperation/coalition, and eight as unidentified. In two studies, two male jaguars formed stable coalitions lasting over 7 years each. In the Llanos, each coalition male paired and mated with several females. For male jaguar coalitions, we documented similar behaviours as recorded earlier in lions or cheetahs, which included patrolling and marking territory together, invading territories of other males, collaborative chasing and killing other jaguars, and sharing prey. However, different from lions or cheetahs, associated male jaguars spent less time together, did not cooperate with females, and did not hunt cooperatively together. Our analysis of literature suggested that male jaguar coalitions were more likely to form when females had small home range size, a proxy of females’ concentration, while in lions, the male group size was directly correlated with the female group size. Similarly, locally concentrated access to females may drive formation of male coalitions in cheetahs. We conclude that high biomass and aggregation of prey are likely drivers of sociality in felids.

Significance statement

The division into social and solitary species in large felids has so far seemed unambiguous, with only lions and cheetahs classified as social species, in which male coalitions also occurred. Our data show that, under certain conditions, male coalitions may also form in jaguar populations. Factors that drive formation of male coalitions in lions and cheetahs, but not in other species of large cats, have not been clear until now. Our analyses indicate that in jaguars, lions, and cheetahs, the concentration of females likely plays the most important role. In jaguars, the probability of male coalition occurrence is highest in populations with the smallest mean female home range size (and thus likely high local density of females), while in lions, male group size is most strongly correlated with female group size.

Determining Significant Correlation Between Pairs of Extant Characters in a Small Parsimony Framework

When studying the evolutionary relationships among a set of species, the principle of parsimony states that a relationship involving the fewest number of evolutionary events is likely the correct one. Due to its simplicity, this principle was formalized in the context of computational evolutionary biology decades ago by, for example, Fitch and Sankoff. Because the parsimony framework does not require a model of evolution, unlike maximum likelihood or Bayesian approaches, it is often a good starting point when no reasonable estimate of such a model is available. In this work, we devise a method for determining if pairs of discrete characters are significantly correlated across all most parsimonious reconstructions, given a set of species on these characters, and an evolutionary tree. The first step of this method is to use Sankoff's algorithm to compute all most parsimonious assignments of ancestral states (of each character) to the internal nodes of the phylogeny. Correlation between a pair of evolutionary events (e.g., absent to present) for a pair of characters is then determined by the (co-) occurrence patterns between the sets of their respective ancestral assignments. The probability of obtaining a correlation this extreme (or more) under a null hypothesis where the events happen randomly on the evolutionary tree is then used to assess the significance of this correlation. We implement this method: parcours (PARsimonious CO-occURrenceS) and use it to identify significantly correlated evolution among vocalizations and morphological characters in the Felidae family.

Molecular phylogenies map to biogeography better than morphological ones

Phylogenetic relationships are inferred principally from two classes of data: morphological and molecular. Currently, most phylogenies of extant taxa are inferred from molecules and when morphological and molecular trees conflict the latter are often preferred. Although supported by simulations, the superiority of molecular trees has rarely been assessed empirically. Here we test phylogenetic accuracy using two independent data sources: biogeographic distributions and fossil first occurrences. For 48 pairs of morphological and molecular trees we show that, on average, molecular trees provide a better fit to biogeographic data than their morphological counterparts and that biogeographic congruence increases over research time. We find no significant differences in stratigraphic congruence between morphological and molecular trees. These results have implications for understanding the distribution of homoplasy in morphological data sets, the utility of morphology as a test of molecular hypotheses and the implications of analysing fossil groups for which molecular data are unavailable.

Using biogeographical and phylogenetic data, it is shown that molecular trees fit species geographical data better than trees inferred from morphology, and that these differences are not simply due to better tree resolution.

Genomic signatures of divergent ecological strategies in a recent radiation of Neotropical wild cats

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, e.g. a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

Purging of deleterious burden in the endangered Iberian lynx

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.

PWM2Vec: An Efficient Embedding Approach for Viral Host Specification from Coronavirus Spike Sequences

The study of host specificity has important connections to the question about the origin of SARS-CoV-2 in humans which led to the COVID-19 pandemic-an important open question. There are speculations that bats are a possible origin. Likewise, there are many closely related (corona)viruses, such as SARS, which was found to be transmitted through civets. The study of the different hosts which can be potential carriers and transmitters of deadly viruses to humans is crucial to understanding, mitigating, and preventing current and future pandemics. In coronaviruses, the surface (S) protein, or spike protein, is important in determining host specificity, since it is the point of contact between the virus and the host cell membrane. In this paper, we classify the hosts of over five thousand coronaviruses from their spike protein sequences, segregating them into clusters of distinct hosts among birds, bats, camels, swine, humans, and weasels, to name a few. We propose a feature embedding based on the well-known position weight matrix (PWM), which we call PWM2Vec, and we use it to generate feature vectors from the spike protein sequences of these coronaviruses. While our embedding is inspired by the success of PWMs in biological applications, such as determining protein function and identifying transcription factor binding sites, we are the first (to the best of our knowledge) to use PWMs from viral sequences to generate fixed-length feature vector representations, and use them in the context of host classification. The results on real world data show that when using PWM2Vec, machine learning classifiers are able to perform comparably to the baseline models in terms of predictive performance and runtime-in some cases, the performance is better. We also measure the importance of different amino acids using information gain to show the amino acids which are important for predicting the host of a given coronavirus. Finally, we perform some statistical analyses on these results to show that our embedding is more compact than the embeddings of the baseline models.

Effects of reduced mobility on trabecular bone density in captive big cats

Bone responds to elevated mechanical loading by increasing in mass and density. Therefore, wild animals should exhibit greater skeletal mass and density than captive conspecifics. This expectation is pertinent to testing bone functional adaptation theories and to comparative studies, which commonly use skeletal remains that combine zoo and wild-caught specimens. Conservationists are also interested in the effects of captivity on bone morphology as it may influence rewilding success. We compared trabecular bone volume fraction (BVF) between wild and captive mountain lions, cheetahs, leopards and jaguars. We found significantly greater BVF in wild than in captive felids. Effects of captivity were more marked in the humerus than in the femur. A ratio of humeral/femoral BVF was also lower in captive animals and showed a positive relationship to home range size in wild animals. Results are consistent with greater forelimb than hindlimb loading during terrestrial travel, and possibly reduced loading of the forelimb associated with lack of predatory behaviour in captive animals. Thus, captivity among felids has general effects on BVF in the postcranial skeleton and location-specific effects related to limb use. Caution should be exercised when identifying skeletal specimens for use in comparative studies and when rearing animals for conservation purposes.

Phylogenetic analyses of Eurasian lynx (Lynx lynx Linnaeus, 1758) including new mitochondrial DNA sequences from Iran

The Eurasian lynx (Lynx lynx) is one of the widespread felids in Eurasia; however, relatively little is known about the Asian subspecies, and especially the Iranian populations, which comprise the most southwestern part of its range. The current study aimed to assess the phylogenetic status of Iranian populations relative to other populations of Eurasia, by sequencing a 613 bp fragment of the mitochondrial control region. In total, 44 haplotypes were recorded from 83 sequences throughout Eurasia, two of which were found in Iran. The haplotype (H1) is dominant in all Iranian lynx populations and identical to specimens from SW Russia and central China. The second haplotype (H2) is unique and was recorded only from Ghazvin Province in the central Alborz Mountains. Both haplotypes occur in Ghazvin Province. The phylogenetic tree and a median-joining network identified four clades (i.e., East, West 1, West 2, and South). These results are congruent with previous studies and suggest that Eurasian lynx was restricted to the southern part of its range during the glacial maxima and expanded from there to East Asia and to Europe during several independent re-colonization events. The Caucasus region most like plays an important role as a refugium during glacial cycles.

The First Case of Cytauxzoon spp. in Russia: The Parasite Conquers Eurasia

Simple Summary

Cytauxzoon felis was first described more than 40 years ago in the US, later (in the 21st century) similar pathogens were detected in Eurasian and Iberian lynxes, European wildcats and domestic cats in Southern and Central Europe. Our findings have shown the previously unrecorded presence of this parasite in Russia (50 km from Moscow). We described the crucial decrease in the number of leukocytes and erythrocytes, as well as in hemoglobin concentration, throughout the captive serval’s disease, and their increase during the animal’s recovery over six months. Molecular genetic methods allowed us to detect and describe this parasite in four cat species in captivity. The analysis showed high genetic variability and high occurrence of the parasites, which suggests their presence in free-ranging domestic cats and wild felids in Russia.

Abstract

Over the last two decades, Cytauxzoon spp. has been conquering Eurasia, although this fact has only been brought to light through recent more intensive research after the discovery of C. manul in Pallas’ cat. In Europe, Cytauxzoon was detected mainly in southern countries and later in central Europe. This pathogen has now been found in Russia for the first time (50 km from Moscow), this being the most northern sighting in Eurasia. A captive serval (Leptailurus serval) was found to be infected. Hematological analysis showed a crucial decrease in the number of leukocytes and erythrocytes, as well as in hemoglobin concentration. Genetic analysis confirmed the presence of Cytauxzoon spp. in serval blood at the beginning of the disease period. The identical pathogen was found in one bobcat at the same breeding center. Two other haplotypes of Cytauxzoon spp. were obtained from domestic cats at the same location, identical to the samples from Italy. One new haplotype, which was sequenced for the first time, was found in 7/7 investigated Amur wildcats (100%). The high occurrence and diversity of these pathogens suggest that they are present in free-ranging domestic cats and wild felids in Russia, and may be considered a potential threat to the endangered species. Current research shows that the genetic diversity of this pathogen may be even higher than it was suggested previously. Further genetic research is necessary to describe the diversity and phylogeny of this pathogen in Eurasia.

Comparative phylogeography among eight Neotropical wild cat species: no single evolutionary pattern

The felid species of South America are thought to have arrived on the continent during the Great American Biotic Interchange (GABI) in the Pleistocene. However, molecular and palaeontological data do not agree on how this event affected speciation in felids. Here, we determine both the number of colonization events and the period when felines first migrated from North America to South America. In addition, we evaluate whether similar evolutionary events could have affected the eight Neotropical cat species in their levels of genetic diversity, spatial genetic structure and demographic changes. We analysed four concatenated mitochondrial genes of the jaguar, ocelot, margay, tigrina, pampas cat, Andean cat, puma and jaguarundi. The samples were representative of a wide distribution of these species in Central and South America. Our analysis suggests either three or four colonization events from North America to South America over the past 3 Myr, followed by subsequent speciation events and the attainment of high or very high genetic diversity levels for seven of the species. No unique evolutionary process was detected for any of the current Neotropical cat species.

Recurrent erosion of COA1/MITRAC15 exemplifies conditional gene dispensability in oxidative phosphorylation

Skeletal muscle fibers rely upon either oxidative phosphorylation or the glycolytic pathway with much less reliance on oxidative phosphorylation to achieve muscular contractions that power mechanical movements. Species with energy-intensive adaptive traits that require sudden bursts of energy have a greater dependency on glycolytic fibers. Glycolytic fibers have decreased reliance on OXPHOS and lower mitochondrial content compared to oxidative fibers. Hence, we hypothesized that gene loss might have occurred within the OXPHOS pathway in lineages that largely depend on glycolytic fibers. The protein encoded by the COA1/MITRAC15 gene with conserved orthologs found in budding yeast to humans promotes mitochondrial translation. We show that gene disrupting mutations have accumulated within the COA1 gene in the cheetah, several species of galliform birds, and rodents. The genomic region containing COA1 is a well-established evolutionary breakpoint region in mammals. Careful inspection of genome assemblies of closely related species of rodents and marsupials suggests two independent COA1 gene loss events co-occurring with chromosomal rearrangements. Besides recurrent gene loss events, we document changes in COA1 exon structure in primates and felids. The detailed evolutionary history presented in this study reveals the intricate link between skeletal muscle fiber composition and the occasional dispensability of the chaperone-like role of the COA1 gene.

Rapid Macrosatellite Evolution Promotes X-Linked Hybrid Male Sterility in a Feline Interspecies Cross

The sterility or inviability of hybrid offspring produced from an interspecific mating result from incompatibilities between parental genotypes that are thought to result from divergence of loci involved in epistatic interactions. However, attributes contributing to the rapid evolution of these regions also complicates their assembly, thus discovery of candidate hybrid sterility loci is difficult and has been restricted to a small number of model systems. Here we reported rapid interspecific divergence at the DXZ4 macrosatellite locus in an interspecific cross between two closely related mammalian species: the domestic cat (Felis silvestris catus) and the Jungle cat (Felis chaus). DXZ4 is an interesting candidate due to its structural complexity, copy number variability, and described role in the critical yet complex biological process of X-chromosome inactivation. However, the full structure of DXZ4 was absent or incomplete in nearly every available mammalian genome assembly given its repetitive complexity. We compared highly continuous genomes for three cat species, each containing a complete DXZ4 locus, and discovered that the felid DXZ4 locus differs substantially from the human ortholog, and that it varies in copy number between cat species. Additionally, we reported expression, methylation, and structural conformation profiles of DXZ4 and the X chromosome during stages of spermatogenesis that have been previously associated with hybrid male sterility. Collectively, these findings suggest a new role for DXZ4 in male meiosis and a mechanism for feline interspecific incompatibility through rapid satellite divergence.