Phylogenetic studies of pantherine cats (Felidae) based on multiple genes, with novel application of nuclear beta-fibrinogen intron 7 to carnivores

The pantherine lineage of the cat family Felidae (order: Carnivora) includes five big cats of genus Panthera and a great many midsized cats known worldwide. Presumably because of their recent and rapid radiation, the evolutionary relationship among pantherines remains ambiguous. We provide an independent assessment of the evolutionary history of pantherine lineage using two complete mitochondrial (mt) genes (ND2 and ND4) and the nuclear beta-fibrinogen intron 7 gene, whose utility in carnivoran phylogeny was first explored. The available four mt (ND5, cytb, 12S, and 16SrRNA) and two nuclear (IRBP and TTR) sequence loci were also combined to reconstruct phylogeny of 14 closely related cat species. Our analyses of combined mt data (six genes; approximately 3750 bp) and combined mt and nuclear data (nine genes; approximately 6500 bp) obtained identical tree topologies, which were well-resolved and strongly supported for almost all nodes. Monophyly of Panthera genus in pantherine lineage was confirmed and interspecific affinities within this genus revealed a novel branching pattern, with P. tigris diverging first in Panthera genus, followed by P. onca, P. leo, and last two sister species P. pardus and P. uncia. In addition, close association of Neofelis nebulosa to Panthera, the phylogenetic redefinition of Otocolobus manul within the domestic cat group, and the relatedness of Acinonyx jubatus and Puma concolor were all important findings in the resulting phylogenies. The potential utilities of nine different genes for phylogenetic resolution of closely related pantherine species were also evaluated, with special interest in that of the novel nuclear beta-fibrinogen intron 7.     

Patterns of Y and X chromosome DNA sequence divergence during the Felidae radiation.

The 37 species of modern cats have evolved from approximately eight phylogenetic lineages within the past 10 to 15 million years. The Felidae family has been described with multiple measures of morphologic and molecular evolutionary methods that serve as a framework for tracking gene divergence during brief evolutionary periods. In this report, we compare the mode and tempo of evolution of noncoding sequences of a large intron within Zfy (783 bp) and Zfx (854 bp), homologous genes located on the felid Y and X chromosomes, respectively. Zfy sequence variation evolves at about twice the rate of Zfx, and both gene intron sequences track feline hierarchical topologies accurately. As homoplasies are infrequent in patterns of nucleotide substitution, the Y chromosome sequence displays a remarkable degree of phylogenetic consistency among cat species and provides a highly informative glimpse of divergence of sex chromosome sequences in Felidae.     

Resolution of recent radiations within three evolutionary lineages of felidae using mitochondrial restriction fragment length polymorphism variation

Patterns of mitochondrial restriction fragment length polymorphism (RFLP) variation were used to resolve more recent relationships among the species of the Felidae ocelot lineage, domestic cat lineage, and pantherine lineage. Twenty-five of 28 restriction enzymes revealed site variation in at least 1 of 21 cat species. The ocelot lineage was resolved into three separate sistertaxa groups: Geoffroy's cat (Oncifelis geoffroyi) and kodkod (O. guigna), ocelot (Leopardus pardalis) and margay (L. wiedii), and pampas cat (Lynchailurus colocolo) and most of the tigrina samples (Leopardus tigrina). Within the domestic cat lineage, domestic cat (Felis catus), European wild cat (F. silvestris), and African wild cat (F. libyca) formed a monophyletic trichotomy, which was joined with sand cat (F. margarita) to a common ancestor. Jungle cat (F. chaus) and black-footed cat (F. nigripes) mtDNAs diverged earlier than those of the other domestic cat lineage species and are less closely related. Within the pantherine lineage, phylogenetic analysis identified two distinct groups, uniting lion (P. leo) with leopard (P. pardus) and tiger (P. tigris) with snow leopard (P. uncia).To whom correspondence should be addressed.     

猫科动物的分类——分子系统学框架和化石证据

猫科动物的演化关系传统上一直是富有争议的,这很大程度上造成了猫科在属一级的分类众说纷纭。新兴的分子生物学技术为我们理解猫科动物的系统发育提供了新的视角,本文回顾了近些年来的猫科分子系统学研究,以Johnson等（2006）和李钢等（2016）的研究结果为核心,同时结合确切的化石记录进行佐证,推断了猫科的演化历史。我们发现,分子系统学提出的现代猫科动物于晚中新世开始适应辐射、先后分化出了8个支系等观点,与化石证据基本吻合;但是分子系统学对部分支系起源的推断与化石证据存在差异。结合化石记录的分析,我们认为现代猫科8个支系中除狞猫支系和虎猫支系以外,其余均最有可能起源于亚洲,现代猫科动物在演化过程中至少经历了30次洲际迁徙。结合演化历史和形态学研究,我们遵循邱占祥等（2004）的观点,将现存的猫科动物全部归于猫亚科（Felinae）,并将现代猫亚科动物划分为15个属40个物种。     

食肉目猫科物种的系统发育学研究概述

猫科动物(Felidae)是食肉目中肉食性最强的一科, 其中许多成员是人们最熟悉、最引人注目的动物, 也是各地的顶级食肉动物。目前37个现存猫科物种中有36个已经被列为濒危和稀有对象。食肉目猫科物种的进化历史是一个快速辐射和较近时期发生的物种形成事件, 使得猫科物种之间系统发育关系的重建非常困难, 一直处于广泛争论的状态。构建可靠的猫科系统发育关系, 具有重要的进化理论意义和保护生物学价值。文章对猫科物种的系统发育学研究进展, 包括来自于形态学特征、细胞学和分子生物学方面的证据做简要概述, 并提出目前研究中存在的问题。以期对今后猫科物种的系统发育方面的进一步研究工作具有指导意义, 并为该类群的生物多样性资源保护提供科学依据。     

Mitochondrial introgressions into the nuclear genome of the domestic cat

Translocation of mtDNA into the nuclear genome, also referred to as numt, was first reported in the domestic cat (Felis catus) by Lopez et al. (1994). The Lopez-numt consisted of a translocation of 7.9 kbp of mtDNA that inserted into the domestic cat chromosome D2 around 1.8 million years ago. More than a decade later, the release of the domestic cat whole-genome shotgun sequences (1.9x coverage) provides the resource to obtain more comprehensive insight into the extent of mtDNA transfer over time in the domestic cat genome. MegaBLAST searches revealed that the cat genome harbors a wide variety of numts (298 320 bp), one-third of which likely correspond to the Lopez-numt tandem repeat, whereas the remaining numts are probably derived from multiple independent insertions, which in some cases were followed by segmental duplication after insertion in the nucleus. Numts were detected across most cat chromosomes, but the number of numts assigned to chromosomes is underestimated due to the relatively high number of numt sequences with insufficient flanking sequence to map. The catalog of cat numts provides a valuable resource for future studies in Felidae species, including its use as a tool to avoid numt contaminations that may confound population genetics and phylogenetic studies.     

Albumin evolution and phylogenetic relationships among Greek rodents of the families Arvicolidae and Muridae

The phylogenetic relationships of seven rodent species Microtus atticus, M. thomasi, M. epiroticus (family Arvicolidae) and Mus domesticus, Rattus norvegicus, Apodemus flavicollis and A. mystacinus (family Muridae) have been studied. In order to define these relationships we study the albumin evolution using the micro-complement fixation test (MC'F). No phylogenetic (immunological) distance between M. atticus and M. thomasi was found, a fact which confirms from the biochemical point of view the opinion that the former taxon is a synonym of the latter one. A molecular time scale relating MC'F immunological distances and geological time was established based on the assumption of a rate of 100 amino acid substitutions per–20 million years. The time of divergence between M. epiroticus and M. thomasi was estimated to be 0.5–0.6 million years ago (Pleistocene). Such a recent divergence corroborates the opinion based on morphological and protein electrophoretic criteria according to which Terricola (formerly Pitymys ) must be considered as a subgenus of the genus Microtus and not as a distinct genus Pitymys , as previously had been accepted. Apodemus flavicollis and A. mystacinus were separated about 0.65–0.8 million years ago (Pleistocene). The Rattus norvegicus lineage was separated–12.5 million years ago (end of Miocene), shortly before the Mus and Apodemus divergence. Our data indicate that the common ancestor of Arvicolidae and Muridae lived–25 million years ago (early Miocene). All these results are in agreement with paleontological and some recent DNA-DNA hybridization and electrophoretic data.     

Evolution of acid phosphatase-1 in the genus drosophila. Immunological studies

Summary The enzyme acid phosphatase-1 was partially purified from 10 Drosophila species. Four antisera were produced and the ten enzymes were reacted against each serum. The method used to quantitate the reactions involved the electrophoretic separation of antigen-antibody complexes from uncomplexed enzyme, followed by densitometry of the free enzyme. Immunological distances were used to obtain correlation coefficients for all pairwise combinations of the 10 species. From these correlation coefficients, a dendrogram was constructed which is very similar to one diagramming the presumed phylogenetic relationships of the ten species. In addition, the data indicate acid phosphatase-1 has evolved at different rates in different lineages within the genus. A preliminary estimate of the unit evolutionary period for this enzyme is 3.25 million years. The method of determining immunological distances which was used in this study is compared to the method of microcomplement fixation in theDiscussion.     

Exchanges of short polymorphic DNA segments predating speciation in feline major histocompatibility complex class I genes

Sequence comparisons of 14 distinct MHC class I cDNA clones isolated from species representing the three major taxonomic lineages of Felidae (domestic cat lineage, ocelot lineage, and pantherine lineage) revealed that feline MHC class I alleles have highly mosaic structures with short polymorphic sequence motifs that are rearranged between alleles of individual MHC loci, between MHC class I genes within cat species, and between homologous MHC loci in different species. The pattern of sequence variation in felids supports the role of the following factors in production and maintenance of MHC variation: (1) gradual spontaneous mutation; (2) selective pressure to conserve certain residues but also to vary in hypervariable regions, notably residues that functionally participate in antigen recognition and presentation; and (3) recombination-mediated gene exchange between alleles and between related genes. The overall amount of genetic variation observed among MHC class I genes in the Felidae family is no greater than the amount of variation within any outbred cat species (i.e., domestic cat, ocelot). The occurrence of equivalent levels of polymorphism plus the simultaneous persistence of the same sequence motifs in divergent feline species suggest that most MHC class I nucleotide site polymorphism predated species divergences. Ancient polymorphisms have been transmitted through the speciation events and modern feline MHC class I alleles were derived by recombinational exchange of polymorphic sequence motifs. Moreover, some of these sequence motifs were found in other mammalian MHC class I genes, such as classical human HLA-B5, nonclassical human HLA-E class I genes, and bovine class I genes. These results raise the prospect of an ancient origin for some motifs, although the possibility of convergence in parallel mammalian radiations cannot be excluded. Correspondence to: N. Yuhki     

Molecular Phylogeny of Mitochondrial Cytochrome b and 12S rRNA Sequences in the Felidae: Ocelot and Domestic Cat Lineages

Molecular phylogeny of the cat family Felidae is derived using two mitochondrial genes, cytochrome b and 12S rRNA. Phylogenetic methods of weighted maximum parsimony and minimum evolution estimated by neighbor-joining are employed to reconstruct topologies among 20 extant felid species. Sequence analyses of 363 bp of cytochrome b and 376 bp of the 12S rRNA genes yielded average pair-wise similarity values between felids ranging from 94 to 99% and from 85 to 99%, respectively. Phylogenetic reconstruction supports more recent, intralineage associations but fails to completely resolve interlineage relationships. Both genes produce a monophyletic group ofFelisspecies but vary in the placement of the pallas cat. The ocelot lineage represents an early divergence within the Felidae, with strong associations between ocelot and margay, Geoffroy's cat and kodkod, and pampas cat and tigrina. Implications of the relative recency of felid evolution, presence of ancestral polymorphisms, and influence of outgroups in placement of the topological root are discussed.     

Evolution of arginine kinase within the genus Drosophila

The rate of evolution and extant levels of polymorphism for arginine kinase were examined in species of the genus Drosophila. Surveys of 11 species for electrophoretic variation revealed an average heterozygosity of 0.003. Using antisera prepared against arginine kinase purified from a representative of each of three major subdivisions of the genus, immunological distances were measured by microcomplement fixation (MC'F) assay for 20 species of Drosophila. These data are consistent with the broad outline of the phylogenetic relationships within the genus suggested by other kinds of data. The unit evolutionary period for Drosophila arginine kinase was estimated to be 59 million years (MY). This contrasts with an estimate of 30 MY as the unit evolutionary period for vertebrate creatine kinase. It is suggested that this difference in evolutionary rate arises because a single gene encodes Drosophila arginine kinase, whereas at least three different genes encode the various forms of vertebrate creatine kinase.     

Evidence suggesting preadaptation to domestication throughout the small Felidae

One obstacle in the development of a coherent theoretical framework for the process of animal domestication is the rarity of domestication events in human history. It is unclear whether: (1) many species are suitable for domestication, the limiting factor being the requirement of people for new domestic animals; or (2) very few species are preadapted for domestication. Comparisons between 16 species and subspecies of small cats (Felidae) kept in zoos indicated that affiliative behaviour towards people, an important preadaptation to domestication, is widely, if patchily, distributed throughout this taxon, and is not concentrated in species closely related to the domestic cat, Felis silvestris catus . The highest proportion of individuals showing affiliative behaviour was found in the ocelot lineage, which is estimated to have diverged from the rest of the Felidae between 5 and 13 Mya. The domestication of F. silvestris alone among felids is therefore likely to have been the result of a specific set of human cultural events and requirements in the Egyptian New Kingdom, rather than the consequence of a unique tendency to tame-ness in this subspecies.© 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 75 , 2002, 361–366.     

Cross-Species Transmission and Differential Fate of an Endogenous Retrovirus in Three Mammal Lineages

Endogenous retroviruses (ERVs) arise from retroviruses chromosomally integrated in the host germline. ERVs are common in vertebrate genomes and provide a valuable fossil record of past retroviral infections to investigate the biology and evolution of retroviruses over a deep time scale, including cross-species transmission events. Here we took advantage of a catalog of ERVs we recently produced for the bat Myotis lucifugus to seek evidence for infiltration of these retroviruses in other mammalian species (>100) currently represented in the genome sequence database. We provide multiple lines of evidence for the cross-ordinal transmission of a gammaretrovirus endogenized independently in the lineages of vespertilionid bats, felid cats and pangolin ~13–25 million years ago. Following its initial introduction, the ERV amplified extensively in parallel in both bat and cat lineages, generating hundreds of species-specific insertions throughout evolution. However, despite being derived from the same viral species, phylogenetic and selection analyses suggest that the ERV experienced different amplification dynamics in the two mammalian lineages. In the cat lineage, the ERV appears to have expanded primarily by retrotransposition of a single proviral progenitor that lost infectious capacity shortly after endogenization. In the bat lineage, the ERV followed a more complex path of germline invasion characterized by both retrotransposition and multiple infection events. The results also suggest that some of the bat ERVs have maintained infectious capacity for extended period of time and may be still infectious today. This study provides one of the most rigorously documented cases of cross-ordinal transmission of a mammalian retrovirus. It also illustrates how the same retrovirus species has transitioned multiple times from an infectious pathogen to a genomic parasite (i.e. retrotransposon), yet experiencing different invasion dynamics in different mammalian hosts.     

Directional selection in the evolution of elongated upper canines in clouded leopards and sabre‐toothed cats

Extremely developed or specialized traits such as the elongated upper canines of extinct sabre‐toothed cats are often not analogous to those of any extant species, which limits our understanding of their evolutionary cause. However, an extant species may have undergone directional selection for a similar extreme phenotype. Among living felids, the clouded leopard, Neofelis nebulosa, has exceptionally long upper canines for its body size. We hypothesized that directional selection generated the elongated upper canines of clouded leopards in a manner similar to the process in extinct sabre‐toothed cats. To test this, we developed an approach that compared the effect of directional selection among lineages in a phylogeny using a simulation of trait evolution and approximate Bayesian computation. This approach was applied to analyse the evolution of upper canine length in the Felidae phylogeny. Our analyses consistently showed directional selection favouring longer upper canines in the clouded leopard lineage and a lineage leading to the sabre‐toothed cat with the longest upper canines, Smilodon. Most of our analyses detected an effect of directional selection for longer upper canines in the lineage leading to another sabre‐toothed cat, Homotherium, although this selection may have occurred exclusively in the primitive species. In all the analyses, the clouded leopard and Smilodon lineages showed comparable directional selection. This implies that clouded leopards share a selection advantage with sabre‐toothed cats in having elongated upper canines.     

Molecular evolution of mitochondrial 12S RNA and cytochrome b sequences in the pantherine lineage of Felidae

DNA sequence comparisons of two mitochondrial DNA genes were used to infer phylogenetic relationships among 17 Felidae species, notably 15 in the previously described pantherine lineage. The polymerase chain reaction (PCR) was used to generate sequences of 358 base pairs of the mitochondrial 12S RNA gene and 289 base pairs of the cytochrome b protein coding gene. DNA sequences were compared within and between 17 felid and five nonfelid carnivore species. Evolutionary trees were constructed using phenetic, cladistic, and maximum likelihood algorithms. The combined results suggested several phylogenetic relationships including (1) the recognition of a recently evolved monophyletic genus Panthera consisting of Panthera leo, P. pardus, P. onca, P. uncia, P. tigris, and Neofelis nebulosa; (2) the recent common ancestry of Acinonyx jubatus, the African cheetah, and Puma concolor, the American puma; and (3) two golden cat species, Profelis temmincki and Profelis aurata, are not sister species, and the latter is strongly associated with Caracal caracal. These data add to the growing database of vertebrate mtDNA sequences and, given the relatively recent divergence among the felids represented here (1-10 Myr), allow 12S and cytochrome b sequence evolution to be addressed over a time scale different from those addressed in most work on vertebrate mtDNA.      

Evolution of a major drug metabolizing enzyme defect in the domestic cat and other felidae: phylogenetic timing and the role of hypercarnivory

The domestic cat (Felis catus) shows remarkable sensitivity to the adverse effects of phenolic drugs, including acetaminophen and aspirin, as well as structurally-related toxicants found in the diet and environment. This idiosyncrasy results from pseudogenization of the gene encoding UDP-glucuronosyltransferase (UGT) 1A6, the major species-conserved phenol detoxification enzyme. Here, we established the phylogenetic timing of disruptive UGT1A6 mutations and explored the hypothesis that gene inactivation in cats was enabled by minimal exposure to plant-derived toxicants. Fixation of the UGT1A6 pseudogene was estimated to have occurred between 35 and 11 million years ago with all extant Felidae having dysfunctional UGT1A6. Out of 22 additional taxa sampled, representative of most Carnivora families, only brown hyena (Parahyaena brunnea) and northern elephant seal (Mirounga angustirostris) showed inactivating UGT1A6 mutations. A comprehensive literature review of the natural diet of the sampled taxa indicated that all species with defective UGT1A6 were hypercarnivores (>70% dietary animal matter). Furthermore those species with UGT1A6 defects showed evidence for reduced amino acid constraint (increased dN/dS ratios approaching the neutral selection value of 1.0) as compared with species with intact UGT1A6. In contrast, there was no evidence for reduced amino acid constraint for these same species within UGT1A1, the gene encoding the enzyme responsible for detoxification of endogenously generated bilirubin. Our results provide the first evidence suggesting that diet may have played a permissive role in the devolution of a mammalian drug metabolizing enzyme. Further work is needed to establish whether these preliminary findings can be generalized to all Carnivora.     

Skull morphology and functionality of extant Felidae (Mammalia: Carnivora): a phylogenetic and evolutionary perspective

Felids morphology and ecological role as hypercarnivores are quite constant, despite considerable body size variation among species. Skull morphological and functional features of 34 extant cat species were evaluated under a phylogenetic framework of the Felidae. Twenty skull measurements were analysed through Principal Component Analysis to assess the species morphofunctional spaces. Force indexes were obtained from static equilibrium equations to infer jaw mechanics. Correlations between morphological, functional, and ecological traits were tested by phylogenetically independent contrasts. In spite of the general cat‐like pattern, specific features on the skulls allowed differentiation among groups. Acinonyx jubatus, for instance, showed a shorter and shallower temporal fossa than other big cats, and their bite functionality is marked by an increased contribution of the masseteric system. A morphofunctional dichotomy between Neotropical and Eurasian/African small cats was detected, and is associated with the major transversal axes of the skulls. According to the contrast analyses, the skull size is correlated with the bite force and prey size, but it is uncorrelated with the variations on jaw mechanics (from temporalis or masseter muscle optimizations). Also, there was no correlation between functional differences on jaw muscles and the ratio of prey weight to cat weight. The efficiency of the jaw apparatus among cats is quite consistent; therefore, the different evolutionary trends of jaw mechanics seem to be caused by the casuistic fixation of phenotypical variations, rather than by specific adaptative selections. © 2010 The Linnean Society of London, Zoological Journal of the Linnean Society, 2011, 161 , 414–462.     

Sympatric Asian felid phylogeography reveals a major Indochinese–Sundaic divergence

The dynamic geological and climatological history of Southeast Asia has spawned a complex array of ecosystems and 12 of the 37 known cat species, making it the most felid‐rich region in the world. To examine the evolutionary histories of these poorly studied fauna, we compared phylogeography of six species (leopard cat Prionailurus bengalensis, fishing cat P. viverrinus, Asiatic golden cat Pardofelis temminckii, marbled cat P. marmorata, tiger Panthera tigris and leopard P. pardus) by sequencing over 5 kb of DNA each from 445 specimens at multiple loci of mtDNA, Y and X chromosomes. All species except the leopard displayed significant phylogenetic partitions between Indochina and Sundaland, with the central Thai–Malay Peninsula serving as the biogeographic boundary. Concordant mtDNA and nuclear DNA genealogies revealed deep Indochinese–Sundaic divergences around 2 MYA in both P. bengalensis and P. marmorata comparable to previously described interspecific distances within Felidae. The divergence coincided with serial sea level rises during the late Pliocene and early Pleistocene, and was probably reinforced by repeated isolation events associated with environmental changes throughout the Pleistocene. Indochinese–Sundaic differentiations within P. tigris and P. temminckii were more recent at 72–108 and 250–1570 kya, respectively. Overall, these results illuminate unexpected, deep vicariance events in Southeast Asian felids and provide compelling evidence of species‐level distinction between the Indochinese and Sundaic populations in the leopard cat and marbled cat. Broader sampling and further molecular and morphometric analyses of these species will be instrumental in defining conservation units and effectively preserving Southeast Asian biodiversity.     

Expression of an onco-developmental antegen among avian species.

Rabbit antisera directed against an onco-developmental antigen on chicken red blood cells have been serologically dissected through specific adsorptions. It is now possible to detect 13 antigenic determinants with the fractionated antisera. The onco-developmental antigen referred to as chicken fetal-leukemic antigen (CFA) is fetal-specific in the white Leghorn chicken, being present on the embryonic but not adult peripheral red blood cells of non-being present on the embryonic but not adult peripheral red blood cells of non-leukemic birds. However, one or more of the onco-developmental antigenic determinants have been detected on adult peripheral red blood cells of non-Gallus avian species, as well as on red blood cells from two adult chicken varieties. For phylogenetic purposes, red blood cells from avian species were characterized for their combinations of CFA determinants. Comparisons among species revealed specific patterns of antigenic expression within phylogenetic groups. Several CFA determinants were restricted in their occurrence to species within a single family, and one determinant was found in all cases where CFA was expressed. The distribution of CFA determinants was used to determine immunological distances among four Galliform species. These distances agreed with the immunological relationships established using different serological markers.     

Evolutionary analysis of a large mtDNA translocation (numt) into the nuclear genome of the Panthera genus species

Translocation of cymtDNA into the nuclear genome, also referred to as numt, has been reported in many species, including several closely related to the domestic cat (Felis catus). We describe the recent transposition of 12,536 bp of the 17 kb mitochondrial genome into the nucleus of the common ancestor of the five Panthera genus species: tiger, P. tigris; snow leopard, P. uncia; jaguar, P. onca; leopard, P. pardus; and lion, P. leo. This nuclear integration, representing 74% of the mitochondrial genome, is one of the largest to be reported in eukaryotes. The Panthera genus numt differs from the numt previously described in the Felis genus in: (1) chromosomal location (F2-telomeric region vs. D2-centromeric region), (2) gene make up (from the ND5 to the ATP8 vs. from the CR to the COII), (3) size (12.5 vs. 7.9 kb), and (4) structure (single monomer vs. tandemly repeated in Felis). These distinctions indicate that the origin of this large numt fragment in the nuclear genome of the Panthera species is an independent insertion from that of the domestic cat lineage, which has been further supported by phylogenetic analyses. The tiger cymtDNA shared around 90% sequence identity with the homologous numt sequence, suggesting an origin for the Panthera numt at around 3.5 million years ago, prior to the radiation of the five extant Panthera species.