A cost-effective blood DNA methylation-based age estimation method in domestic cats, Tsushima leopard cats (Prionailurus bengalensis euptilurus) and Panthera species, using targeted bisulphite sequencing and machine learning models

Individual age can be used to design more efficient and suitable management plans in both in situ and ex situ conservation programmes for targeted wildlife species. DNA methylation is a promising marker of epigenetic ageing that can accurately estimate age from small amounts of biological material, which can be collected in a minimally invasive manner. In this study, we sequenced five targeted genetic regions and used 8-23 selected CpG sites to build age estimation models using machine learning methods at only about $3-7 per sample. Blood samples of seven Felidae species were used, ranging from small to big, and domestic to endangered species: domestic cats (Felis catus, 139 samples), Tsushima leopard cats (Prionailurus bengalensis euptilurus, 84 samples) and five Panthera species (96 samples). The models achieved satisfactory accuracy, with the mean absolute error of the most accurate models recorded at 1.966, 1.348 and 1.552 years in domestic cats, Tsushima leopard cats and Panthera spp. respectively. We developed the models in domestic cats and Tsushima leopard cats, which were applicable to individuals regardless of health conditions; therefore, these models are applicable to samples collected from individuals with diverse characteristics, which is often the case in conservation. We also showed the possibility of developing universal age estimation models for the five Panthera spp. using only two of the five genetic regions. We do not recommend building a common age estimation model for all the target species using our markers, because of the degraded performance of models that included all species.

The genome of the black-footed cat: Revealing a rich natural history and urgent conservation priorities for small felids

Habitat degradation and loss of genetic diversity are common threats faced by almost all of today's wild cats. Big cats, such as tigers and lions, are of great concern and have received considerable conservation attention through policies and international actions. However, knowledge of and conservation actions for small wild cats are lagging considerably behind. The black-footed cat, Felis nigripes, one of the smallest felid species, is experiencing increasing threats with a rapid reduction in population size. However, there is a lack of genetic information to assist in developing effective conservation actions. A de novo assembly of a high-quality chromosome-level reference genome of the black-footed cat was made, and comparative genomics and population genomics analyses were carried out. These analyses revealed that the most significant genetic changes in the evolution of the black-footed cat are the rapid evolution of sensory and metabolic-related genes, reflecting genetic adaptations to its characteristic nocturnal hunting and a high metabolic rate. Genomes of the black-footed cat exhibit a high level of inbreeding, especially for signals of recent inbreeding events, which suggest that they may have experienced severe genetic isolation caused by habitat fragmentation. More importantly, inbreeding associated with two deleterious mutated genes may exacerbate the risk of amyloidosis, the dominant disease that causes mortality of about 70% of captive individuals. Our research provides comprehensive documentation of the evolutionary history of the black-footed cat and suggests that there is an urgent need to investigate genomic variations of small felids worldwide to support effective conservation actions.

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.

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.

Simulating impacts of rapid forest loss on population size, connectivity and genetic diversity of Sunda clouded leopards (Neofelis diardi) in Borneo

Habitat loss is the greatest threat to biodiversity in Borneo, and to anticipate and combat its effects it is important to predict the pattern of loss and its consequences. Borneo is a region of extremely high biodiversity from which forest is being lost faster than in any other. The little-known Sunda clouded leopard (Neofelis diardi) is the top predator in Borneo and is likely to depend critically on habitat connectivity that is currently being rapidly lost to deforestation. We modeled the effects of landscape fragmentation on population size, genetic diversity and population connectivity for the Sunda clouded leopard across the entirety of Borneo. We modelled the impacts of land use change between the years 2000, 2010 and projected forwards to 2020. We found substantial reductions across all metrics between 2000 and 2010: the proportion of landscape connected by dispersal fell by approximately 12.5% and the largest patch size declined by around 15.1%, leading to a predicted 11.4% decline in clouded leopard numbers. We also predict that these trends will accelerate greatly towards 2020, with the percentage of the landscape being connected by dispersal falling by about 57.8%, the largest patch size falling by around 62.8% and the predicted clouded leopard population falling by 62.5% between 2010 and 2020. We predicted that these large declines in clouded leopard population size and connectivity will also substantially reduce the genetic diversity of the remaining clouded leopard population.

Microsatellite DNA analysis reveals lower than expected genetic diversity in the threatened leopard cat (Prionailurus bengalensis) in South Korea

To optimize conservation efforts, it is necessary to determine the risk of extinction by collecting reliable population information for a given species. We developed eight novel, polymorphic microsatellite markers and used these markers in conjunction with twelve existing markers to measure genetic diversity of South Korean populations of leopard cat (Prionailurus bengalensis), a species for which population size and habitat area data are unknown in the country, to assess its conservation status. The average number of alleles and the observed heterozygosity of the species were 3.8 and 0.41, respectively, and microsatellite diversity was lower than the average genetic diversity of 57 populations of 12 other felid species, and lower than that of other mammal populations occurring in South Korea, including the raccoon dog (Nyctereutes procyonoides), water deer (Hydropotes inermis), and endangered long-tailed goral (Naemorhedus caudatus). Furthermore, analysis of genetic structure in the national leopard cat population showed no clear genetic differentiation, suggesting that it is not necessary to divide the South Korean leopard cat population into multiple management units for the purposes of conservation. These results indicate that the genetic diversity of the leopard cat in South Korea is unexpectedly low, and that the risk of local extinction is, as a result, substantial. Thus, it is necessary to begin appropriate conservation efforts at a national level to conserve the leopard cat population in South Korea.

Comparing hair-morphology and molecular methods to identify fecal samples from Neotropical felids

To avoid certain problems encountered with more-traditional and invasive methods in behavioral-ecology studies of mammalian predators, such as felids, molecular approaches have been employed to identify feces found in the field. However, this method requires a complete molecular biology laboratory, and usually also requires very fresh fecal samples to avoid DNA degradation. Both conditions are normally absent in the field. To address these difficulties, identification based on morphological characters (length, color, banding, scales and medullar patterns) of hairs found in feces could be employed as an alternative. In this study we constructed a morphological identification key for guard hairs of eight Neotropical felids (jaguar, oncilla, Geoffroy's cat, margay, ocelot, Pampas cat, puma and jaguarundi) and compared its efficiency to that of a molecular identification method, using the ATP6 region as a marker. For this molecular approach, we simulated some field conditions by postponing sample-conservation procedures. A blind test of the identification key obtained a nearly 70% overall success rate, which we considered equivalent to or better than the results of some molecular methods (probably due to DNA degradation) found in other studies. The jaguar, puma and jaguarundi could be unequivocally discriminated from any other Neotropical felid. On a scale ranging from inadequate to excellent, the key proved poor only for the margay, with only 30% of its hairs successfully identified using this key; and have intermediate success rates for the remaining species, the oncilla, Geoffroy's cat, ocelot and Pampas cat, were intermediate. Complementary information about the known distributions of felid populations may be necessary to substantially improve the results obtained with the key. Our own molecular results were even better, since all blind-tested samples were correctly identified. Part of these identifications were made from samples kept in suboptimal conditions, with some samples remaining outdoors for up to seven days, simulating conditions in the field. It appears that both methods can be used, depending on the available laboratory facilities and on the expected results.

[Reproductive Biology and a Genome Resource Bank of Felidae]

The main achievements in applying modern reproductive technologies to the banking of the genetic resources of the Felidae family are reviewed. The classification of felids at the level of species and subspecies is revised in the light of recent molecular data. Special emphasis is made on such mainstream technologies as semen collection and cryopreservation followed by artificial insemination, as well as on in vitro maturation and fertilization of oocytes combined with the culture of in vitro-derived felid embryos.

Synergistic roles of climate warming and human occupation in Patagonian megafaunal extinctions during the Last Deglaciation

The causes of Late Pleistocene megafaunal extinctions (60,000 to 11,650 years ago, hereafter 60 to 11.65 ka) remain contentious, with major phases coinciding with both human arrival and climate change around the world. The Americas provide a unique opportunity to disentangle these factors as human colonization took place over a narrow time frame (~15 to 14.6 ka) but during contrasting temperature trends across each continent. Unfortunately, limited data sets in South America have so far precluded detailed comparison. We analyze genetic and radiocarbon data from 89 and 71 Patagonian megafaunal bones, respectively, more than doubling the high-quality Pleistocene megafaunal radiocarbon data sets from the region. We identify a narrow megafaunal extinction phase 12,280 ± 110 years ago, some 1 to 3 thousand years after initial human presence in the area. Although humans arrived immediately prior to a cold phase, the Antarctic Cold Reversal stadial, megafaunal extinctions did not occur until the stadial finished and the subsequent warming phase commenced some 1 to 3 thousand years later. The increased resolution provided by the Patagonian material reveals that the sequence of climate and extinction events in North and South America were temporally inverted, but in both cases, megafaunal extinctions did not occur until human presence and climate warming coincided. Overall, metapopulation processes involving subpopulation connectivity on a continental scale appear to have been critical for megafaunal species survival of both climate change and human impacts.

Using a Novel Absolute Ontogenetic Age Determination Technique to Calculate the Timing of Tooth Eruption in the Saber-Toothed Cat, Smilodon fatalis

Despite the superb fossil record of the saber-toothed cat, Smilodon fatalis, ontogenetic age determination for this and other ancient species remains a challenge. The present study utilizes a new technique, a combination of data from stable oxygen isotope analyses and micro-computed tomography, to establish the eruption rate for the permanent upper canines in Smilodon fatalis. The results imply an eruption rate of 6.0 millimeters per month, which is similar to a previously published average enamel growth rate of the S. fatalis upper canines (5.8 millimeters per month). Utilizing the upper canine growth rate, the upper canine eruption rate, and a previously published tooth replacement sequence, this study calculates absolute ontogenetic age ranges of tooth development and eruption in S. fatalis. The timing of tooth eruption is compared between S. fatalis and several extant conical-toothed felids, such as the African lion (Panthera leo). Results suggest that the permanent dentition of S. fatalis, except for the upper canines, was fully erupted by 14 to 22 months, and that the upper canines finished erupting at about 34 to 41 months. Based on these developmental age calculations, S. fatalis individuals less than 4 to 7 months of age were not typically preserved at Rancho La Brea. On the whole, S. fatalis appears to have had delayed dental development compared to dental development in similar-sized extant felids. This technique for absolute ontogenetic age determination can be replicated in other ancient species, including non-saber-toothed taxa, as long as the timing of growth initiation and growth rate can be determined for a specific feature, such as a tooth, and that growth period overlaps with the development of the other features under investigation.

Evolutionary allometry of the thoracolumbar centra in felids and bovids

Mammals have evolved a remarkable range of body sizes, yet their overall body plan remains unaltered. One challenge of evolutionary biology is to understand the mechanisms by which this size diversity is achieved, and how the mechanical challenges associated with changing body size are overcome. Despite the importance of the axial skeleton in body support and locomotion, and much interest in the allometry of the appendicular skeleton, little is known about vertebral allometry outside primates. This study compares evolutionary allometry of the thoracolumbar centra in two families of quadrupedal running mammals: Felidae and Bovidae. I test the hypothesis that, as size increases, the thoracolumbar region will resist increasing loads by becoming a) craniocaudally shorter, and b) larger in cross-sectional area, particularly in the sagittal plane. Length, width, and height of the thoracolumbar centra of 23 felid and 34 bovid species were taken. Thoracic, prediaphragmatic, lumbar, and postdiaphragmatic lengths were calculated, and diameters were compared at three equivalent positions: the midthoracic, the diaphragmatic and the midlumbar vertebra. Allometric slopes were calculated using a reduced major axis regression, on both raw and independent contrasts data. Slopes and elevations were compared using an ANCOVA. As size increases the thoracolumbar centra become more robust, showing preferential reinforcement in the sagittal plane. There was less allometric shortening of the thoracic than the lumbar region, perhaps reflecting constraints due to its connection with the respiratory apparatus. The thoracic region was more robust in bovids than felids, whereas the lumbar region was longer and more robust in felids than bovids. Elongation of lumbar centra increases the outlever of sagittal bending at intervertebral joints, increasing the total pelvic displacement during dorsomobile running. Both locomotor specializations and functional regionalization of the axial skeleton appear to have influenced its response to increasing size.

Evolution in the sabre-tooth cat, Smilodon fatalis, in response to Pleistocene climate change

The late Pleistocene was a time of environmental change, culminating in an extinction event. Few fossil localities record a temporal series of carnivore fossil populations from this interesting interval as well as Rancho La Brea (RLB). We analysed mandibles of Smilodon fatalis from RLB using 2-D geometric morphometrics to examine whether, and how, mandibular shape changes through time. Smilodon fatalis shows mandibular evolution with oscillations between a small, ancestral-type morph in pits 77 (≈37 Kybp) and 2051 (≈26 Kybp), a larger, more derived morph in pits 91 (≈28 Kybp) and 61-67 (≈13.6 Kybp), and an intermediate morph from pit 13 (≈17.7 Kybp). These oscillations end in pit 61-67, with greatest body size, and are estimated to have its widest gape and lowest bite force. Additionally, variation is lowest in pit 61-67, which was deposited concurrent with the Bølling–Allerød warming event, which may have important implications for the timing or conditions during the extinction event. Contra to a temporal Bergmann's rule, such rapid warming events appear to be correlated with larger, derived, morphologies whereas static, cooler, climates correlate with gracile, ancestral morphologies.

Comparative assessment of genetic and morphological variation at an extensive hybrid zone between two wild cats in southern Brazil

Increased attention towards the Neotropical cats Leopardus guttulus and L. geoffroyi was prompted after genetic studies identified the occurrence of extensive hybridization between them at their geographic contact zone in southern Brazil. This is a region where two biomes intersect, each of which is associated with one of the hybridizing species (Atlantic Forest with L. guttulus and Pampas with L. geoffroyi). In this study, we conducted in-depth analyses of multiple molecular markers aiming to characterize the magnitude and spatial structure of this hybrid zone. We also performed a morphological assessment of these species, aiming to test their phenotypic differentiation at the contact zone, as well as the correlation between morphological features and the admixture status of the individuals. We found strong evidence for extensive and complex hybridization, with at least 40% of the individuals sampled in Rio Grande do Sul state (southernmost Brazil) identified as hybrids resulting from post-F1 generations. Despite such a high level of hybridization, samples collected in this state still comprised two recognizable clusters (genetically and morphologically). Genetically pure individuals were sampled mainly in regions farther from the contact zone, while hybrids concentrated in a central region (exactly at the interface between the two biomes). The morphological data set also revealed a strong spatial structure, which was correlated with the molecular results but displayed an even more marked separation between the clusters. Hybrids often did not present intermediate body sizes and could not be clearly distinguished morphologically from the parental forms. This observation suggests that some selective pressure may be acting on the hybrids, limiting their dispersal away from the hybrid zone and perhaps favoring genomic combinations that maintain adaptive phenotypic features of one or the other parental species.

The dynamic proliferation of CanSINEs mirrors the complex evolution of Feliforms

BACKGROUND

Repetitive short interspersed elements (SINEs) are retrotransposons ubiquitous in mammalian genomes and are highly informative markers to identify species and phylogenetic associations. Of these, SINEs unique to the order Carnivora (CanSINEs) yield novel insights on genome evolution in domestic dogs and cats, but less is known about their role in related carnivores. In particular, genome-wide assessment of CanSINE evolution has yet to be completed across the Feliformia (cat-like) suborder of Carnivora. Within Feliformia, the cat family Felidae is composed of 37 species and numerous subspecies organized into eight monophyletic lineages that likely arose 10 million years ago. Using the Felidae family as a reference phylogeny, along with representative taxa from other families of Feliformia, the origin, proliferation and evolution of CanSINEs within the suborder were assessed.

RESULTS

We identified 93 novel intergenic CanSINE loci in Feliformia. Sequence analyses separated Feliform CanSINEs into two subfamilies, each characterized by distinct RNA polymerase binding motifs and phylogenetic associations. Subfamily I CanSINEs arose early within Feliformia but are no longer under active proliferation. Subfamily II loci are more recent, exclusive to Felidae and show evidence for adaptation to extant RNA polymerase activity. Further, presence/absence distributions of CanSINE loci are largely congruent with taxonomic expectations within Feliformia and the less resolved nodes in the Felidae reference phylogeny present equally ambiguous CanSINE data. SINEs are thought to be nearly impervious to excision from the genome. However, we observed a nearly complete excision of a CanSINEs locus in puma (Puma concolor). In addition, we found that CanSINE proliferation in Felidae frequently targeted existing CanSINE loci for insertion sites, resulting in tandem arrays.

CONCLUSIONS

We demonstrate the existence of at least two SINE families within the Feliformia suborder, one of which is actively involved in insertional mutagenesis. We find SINEs are powerful markers of speciation and conclude that the few inconsistencies with expected patterns of speciation likely represent incomplete lineage sorting, species hybridization and SINE-mediated genome rearrangement.

Loss of genetic diversity among ocelots in the United States during the 20th century linked to human induced population reductions

Ocelots (Leopardus pardalis) in the United States currently exhibit low levels of genetic diversity. One hypothesis for this observation is that habitat fragmentation, resulting from human induced changes in the landscape during the 20(th) century, created island populations with highly reduced gene flow and increased genetic drift and inbreeding. In an effort to investigate this, we used a portion of the mitochondrial control region and 11 autosomal microsatellite loci to examine historical levels of genetic diversity and infer temporal changes in ocelot populations between 1853 and 2005. Levels of genetic diversity were higher in historical ocelot populations than in extant populations from Texas. The earliest documented loss of mitochondrial haplotype diversity occurred at Laguna Atascosa National Wildlife Refuge. The second extant population inhabiting private lands in Willacy County retained higher levels of genetic diversity through the 1990s, but subsequently lost diversity over the next decade. A similar pattern was observed for autosomal microsatellite loci. This supports the argument that low levels of genetic diversity in Texas are related to human induced population reductions and fragmentation, both of which threaten the remaining ocelots in the United States. At this time, the best means of mitigating the continued erosion of genetic variation are translocation of individuals either from larger populations in Mexico to Texas, or between the Texas populations.

Species detection and identification in sexual organisms using population genetic theory and DNA sequences

Phylogenetic trees of DNA sequences of a group of specimens may include clades of two kinds: those produced by stochastic processes (random genetic drift) within a species, and clades that represent different species. The ratio of the mean pairwise sequence difference between a pair of clades (K) to the mean pairwise sequence difference within a clade (θ) can be used to determine whether the clades are samples from different species (K/θ ≥ 4) or the same species (K/θ<4) with probability ≥ 0.95. Previously I applied this criterion to delimit species of asexual organisms. Here I use data from the literature to show how it can also be applied to delimit sexual species using four groups of sexual organisms as examples: ravens, spotted leopards, sea butterflies, and liverworts. Mitochondrial or chloroplast genes are used because these segregate earlier during speciation than most nuclear genes and hence detect earlier stages of speciation. In several cases the K/θ ratio was greater than 4, confirming the original authors' intuition that the clades were sufficiently different to be assigned to different species. But the K/θ ratio split each of two liverwort species into two evolutionary species, and showed that support for the distinction between the common and Chihuahuan raven species is weak. I also discuss some possible sources of error in using the K/θ ratio; the most significant one would be cases where males migrate between different populations but females do not, making the use of maternally inherited organelle genes problematic. The K/θ ratio must be used with some caution, like all other methods for species delimitation. Nevertheless, it is a simple theory-based quantitative method for using DNA sequences to make rigorous decisions about species delimitation in sexual as well as asexual eukaryotes.

Convergence and divergence in the evolution of cat skulls: temporal and spatial patterns of morphological diversity

BACKGROUND

Studies of biological shape evolution are greatly enhanced when framed in a phylogenetic perspective. Inclusion of fossils amplifies the scope of macroevolutionary research, offers a deep-time perspective on tempo and mode of radiations, and elucidates life-trait changes. We explore the evolution of skull shape in felids (cats) through morphometric analyses of linear variables, phylogenetic comparative methods, and a new cladistic study of saber-toothed cats.

METHODOLOGY/PRINCIPAL FINDINGS

A new phylogenetic analysis supports the monophyly of saber-toothed cats (Machairodontinae) exclusive of Felinae and some basal felids, but does not support the monophyly of various saber-toothed tribes and genera. We quantified skull shape variation in 34 extant and 18 extinct species using size-adjusted linear variables. These distinguish taxonomic group membership with high accuracy. Patterns of morphospace occupation are consistent with previous analyses, for example, in showing a size gradient along the primary axis of shape variation and a separation between large and small-medium cats. By combining the new phylogeny with a molecular tree of extant Felinae, we built a chronophylomorphospace (a phylogeny superimposed onto a two-dimensional morphospace through time). The evolutionary history of cats was characterized by two major episodes of morphological divergence, one marking the separation between saber-toothed and modern cats, the other marking the split between large and small-medium cats.

CONCLUSIONS/SIGNIFICANCE

Ancestors of large cats in the 'Panthera' lineage tend to occupy, at a much later stage, morphospace regions previously occupied by saber-toothed cats. The latter radiated out into new morphospace regions peripheral to those of extant large cats. The separation between large and small-medium cats was marked by considerable morphologically divergent trajectories early in feline evolution. A chronophylomorphospace has wider applications in reconstructing temporal transitions across two-dimensional trait spaces, can be used in ecophenotypical and functional diversity studies, and may reveal novel patterns of morphospace occupation.

Mitogenomic analysis of the genus Panthera

The complete sequences of the mitochondrial DNA genomes of Panthera tigris, Panthera pardus, and Panthera uncia were determined using the polymerase chain reaction method. The lengths of the complete mitochondrial DNA sequences of the three species were 16990, 16964, and 16773 bp, respectively. Each of the three mitochondrial DNA genomes included 13 protein-coding genes, 22 tRNA, two rRNA, one O(L)R, and one control region. The structures of the genomes were highly similar to those of Felis catus, Acinonyx jubatus, and Neofelis nebulosa. The phylogenies of the genus Panthera were inferred from two combined mitochondrial sequence data sets and the complete mitochondrial genome sequences, by MP (maximum parsimony), ML (maximum likelihood), and Bayesian analysis. The results showed that Panthera was composed of Panthera leo, P. uncia, P. pardus, Panthera onca, P. tigris, and N. nebulosa, which was included as the most basal member. The phylogeny within Panthera genus was N. nebulosa (P. tigris (P. onca (P. pardus, (P. leo, P. uncia)))). The divergence times for Panthera genus were estimated based on the ML branch lengths and four well-established calibration points. The results showed that at about 11.3 MYA, the Panthera genus separated from other felid species and then evolved into the several species of the genus. In detail, N. nebulosa was estimated to be founded about 8.66 MYA, P. tigris about 6.55 MYA, P. uncia about 4.63 MYA, and P. pardus about 4.35 MYA. All these estimated times were older than those estimated from the fossil records. The divergence event, evolutionary process, speciation, and distribution pattern of P. uncia, a species endemic to the central Asia with core habitats on the Qinghai-Tibetan Plateau and surrounding highlands, mostly correlated with the geological tectonic events and intensive climate shifts that happened at 8, 3.6, 2.5, and 1.7 MYA on the plateau during the late Cenozoic period.

Why the leopard got its spots: relating pattern development to ecology in felids

A complete explanation of the diversity of animal colour patterns requires an understanding of both the developmental mechanisms generating them and their adaptive value. However, only two previous studies, which involved computer-generated evolving prey, have attempted to make this link. This study examines variation in the camouflage patterns displayed on the flanks of many felids. After controlling for the effects of shared ancestry using a fully resolved molecular phylogeny, this study shows how phenotypes from plausible felid coat pattern generation mechanisms relate to ecology. We found that likelihood of patterning and pattern attributes, such as complexity and irregularity, were related to felids' habitats, arboreality and nocturnality. Our analysis also indicates that disruptive selection is a likely explanation for the prevalence of melanistic forms in Felidae. Furthermore, we show that there is little phylogenetic signal in the visual appearance of felid patterning, indicating that camouflage adapts to ecology over relatively short time scales. Our method could be applied to any taxon with colour patterns that can reasonably be matched to reaction-diffusion and similar models, where the kinetics of the reaction between two or more initially randomly dispersed morphogens determines the outcome of pattern development.

Influence of ancient glacial periods on the Andean fauna: the case of the pampas cat (Leopardus colocolo)

BACKGROUND

While numerous studies revealed the major role of environmental changes of the Quaternary on the evolution of biodiversity, research on the influence of that period on current South-American fauna is scarce and have usually focused on lowland regions. In this study, the genetic structure of the pampas cat (Leopardus colocolo), a widely distributed felid, was determined and linked to ancient climate fluctuations on the Andean region.

RESULTS

Using both mitochondrial sequences and nuclear microsatellites, we inferred the existence of at least four groups of populations in the central Andes, while other three localities, with little sample sizes (n = 3), presented differences in only one of these markers. The distribution of these groups is correlated to latitude, with a central area characterized by admixture of numerous mitochondrial clades. This suggests colonization from at least three glacial refuges and a contact zone between 20 degrees and 23 degrees S following a glaciation event. The similar coalescence times of the mitochondrial haplotypes indicated that the major clades split approximately one million years ago, likely during the Pre-Pastonian glacial period (0.80 - 1.30 MYA), followed by a demographic expansion in every clade during the Aftonian interglacial period (0.45 - 0.62 MYA). Interestingly, this structure roughly corresponds to the current recognised distribution of morphological subspecies.

CONCLUSION

The four groups of populations identified here must be considered different management units, and we propose the three localities showing differences in only mtDNA or ncDNA as provisional management units. The results revealed the influence of ancient climate fluctuations on the evolutionary history of this species. It is expected that the other species of land vertebrates with a smaller or similar mobility have been affected in the same manner by the glacial and interglacial periods in the central Andes.

[Felines: an alternative in genetic toxicology studies?]

The micronuclei (MN) test carry out in peripheral blood is fast, simple, economic and it is used to detect genotoxic environmental agents. MN are fragments of chromosomes or complete chromosomes remaining in the cytoplasm after cell division, which increase when organisms are exposed to genotoxic agents. Therefore, species with the highest values of spontaneous micronucleated erythrocytes (MNE) are the most suitable to be potentials biomonitor of micronucleogenic agents, using a drop of blood. Nine species of Felines that present spontaneous MNE in peripheral blood are shown. From these species, the cat has been previously proven, with positive results and also lion (Panthera leo), yaguaroundi (Felis yagoaroundi), lynx (Lynx ruffus), jaguar (Panthera onca), puma (Puma concolor), tiger (Panthera tigris), ocelote (Felis padalis) and leopard (Panthera pardus) display spontaneous MNE, and with this characteristic this Family can be propose like a potential group to be used in toxicogenetic studies.

Development of Y chromosome intraspecific polymorphic markers in the Felidae

Y chromosome haplotyping based on microsatellites and single nucleotide polymorphisms (SNPs) has proved to be a powerful tool for population genetic studies of humans. However, the promise of the approach is hampered in the majority of nonhuman mammals by the lack of Y-specific polymorphic markers. We were able to identify new male-specific polymorphisms in the domestic cat Felis catus and 6 additional Felidae species with a combination of molecular genetic and cytogenetic approaches including 1) identifying domestic cat male-specific microsatellites from markers generated from a male cat microsatellite-enriched genomic library, a flow-sorted Y cosmid library, or a Y-specific cat bacteria artificial chromosome (BAC) clone, (2) constructing microsatellite-enriched libraries from flow-sorted Y chromosomes isolated directly from focal wildcat species, and (3) screening Y chromosome conserved anchored tagged sequences primers in Felidae species. Forty-one male-specific microsatellites were identified, but only 6 were single-copy loci, consistent with the repetitive nature of the Y chromosome. Nucleotide diversity (pi) of Y-linked intron sequences (2.1 kbp) was in the range of 0 (tiger) to 9.95 x 10(-4) (marbled cat), and the number of SNPs ranged from none in the tiger to 7 in the Asian leopard cat. The Y haplotyping system described here, consisting of 4 introns (SMCY3, SMCY7, UTY11, and DBY7) and 1 polymorphic microsatellite (SMCY-STR), represents the first available markers for tracking intraspecific male lineage polymorphisms in Felidae species and promises to provide significant insights to evolutionary and population genetic studies of the species.

New taxonomy and the origin of species

The trend for elevating known subspecies to the status of species on the basis of inappropriate evidence can potentially divert important conservation funds away from other species.

Genetic diversity of microsatellite loci in Leopardus pardalis, Leopardus wiedii and Leopardus tigrinus

The microsatellite loci FCA045, FCA077, FCA008, and FCA096 are highly variable molecular markers which were used to determine the genetic diversity in 148 captive Leopardus sp. The PCR-amplified products of microsatellite loci were characterized in ABI Prism 310 Genetic Analyzer. Allele numbers, heterozygosity, polymorphism information content, exclusive allele number, and shared alleles were calculated. Sixty-five alleles were found and their sizes ranged from 116 to 216 bp in four microsatellite loci. The heterozygosity ranged from 0.36 to 0.81 in Leopardus pardalis, 0.57 to 0.67 in L. tigrinus and 0.80 to 0.92 in L. wiedii. The polymorphism information content was from 0.80 to 0.88 in L. pardalis, 0.76 to 0.88 in L. tigrinus and 0.77 to 0.90 in L. wiedii. The margay (L. wiedii) showed the highest index of polymorphism among the three species in this study. These results imply that microsatellite DNA markers can help in the study of the genetic diversity of Leopardus specimens.

Evolution of the male-determining gene SRY within the cat family Felidae

In most placental mammals, SRY is a single-copy gene located on the Y chromosome and is the trigger for male sex determination during embryonic development. Here, we present comparative genomic analyses of SRY (705 bp) along with the adjacent noncoding 5' flank (997 bp) and 3' flank (948 bp) in 36 species of the cat family Felidae. Phylogenetic analyses indicate that the noncoding genomic flanks and SRY closely track species divergence. However, several inconsistencies are observed in SRY. Overall, the gene exhibits purifying selection to maintain function (omega = 0.815) yet SRY is under positive selection in two of the eight felid lineages. SRY has low numbers of nucleotide substitutions, yet most encode amino acid changes between species, and four different species have significantly altered SRY due to insertion/deletions. Moreover, fixation of nonsynonymous substitutions between sister taxa is not consistent and may occur rapidly, as in the case of domestic cat, or not at all over long periods of time, as observed within the Panthera lineage. The former resembles positive selection during speciation, and the latter purifying selection to maintain function. Thus, SRY evolution in cats likely reflects the different phylogeographic histories, selection pressures, and patterns of speciation in modern felids.

Molecular evidence for species-level distinctions in clouded leopards

Among the 37 living species of Felidae, the clouded leopard (Neofelis nebulosa) is generally classified as a monotypic genus basal to the Panthera lineage of great cats. This secretive, mid-sized (16-23 kg) carnivore, now severely endangered, is traditionally subdivided into four southeast Asian subspecies (Figure 1A). We used molecular genetic methods to re-evaluate subspecies partitions and to quantify patterns of population genetic variation among 109 clouded leopards of known geographic origin (Figure 1A, Tables S1 ans S2 in the Supplemental Data available online). We found strong phylogeographic monophyly and large genetic distances between N. n. nebulosa (mainland) and N. n. diardi (Borneo; n = 3 individuals) with mtDNA (771 bp), nuclear DNA (3100 bp), and 51 microsatellite loci. Thirty-six fixed mitochondrial and nuclear nucleotide differences and 20 microsatellite loci with nonoverlapping allele-size ranges distinguished N. n. nebulosa from N. n. diardi. Along with fixed subspecies-specific chromosomal differences, this degree of differentiation is equivalent to, or greater than, comparable measures among five recognized Panthera species (lion, tiger, leopard, jaguar, and snow leopard). These distinctions increase the urgency of clouded leopard conservation efforts, and if affirmed by morphological analysis and wider sampling of N. n. diardi in Borneo and Sumatra, would support reclassification of N. n. diardi as a new species (Neofelis diardi).

The geography of mammalian speciation: mixed signals from phylogenies and range maps

The importance of geographic isolation in speciation has been debated since the 19th century. Since the beginning of the 20th century, the consensus has been that most speciation involves divergence in allopatry. This consensus was based largely on decades of observations by naturalists and verbal arguments against speciation without isolation. Recent attempts to quantify the importance of allopatric versus sympatric speciation using comparative methods called "age-range correlation" (ARC) suggest that allopatric speciation is more common than sympatric speciation. However, very few taxa have been studied and there are concerns about the adequacy of the methods. We propose methodological improvements including changes in the way overlap between clades is quantified and Monte Carlo methods to test the null hypothesis of no relationship between phylogenetic relatedness and geographic range overlap. We analyze 14 clades of mammals, chosen because of the availability of data and the consensus among mammalogists that speciation is routinely allopatric. Although data from a few clades clearly indicate allopatric speciation, divergence with gene flow is plausible in others and many results are inconclusive. The relative rarity of significant correlations between phylogenetic distance and range overlap may have three distinct causes: (1) post-speciation range changes, (2) relative rarity of range overlap, and (3) a mixture of geographic modes of speciation. Our results support skepticism about ARC's power for inferring the biogeography of speciation. Yet, even if few clades provide clear signals, meta-analytic approaches such as ARC may set bounds on the prevalence of alternative modes of speciation.

Follicle-stimulating hormone receptor in felids: intra- and interspecies variation

Assisted reproductive technologies are increasingly applied to support breeding efforts for many endangered felids. To explain the highly variable responses among felids to exogenous gonadotropins (FSH, eCG), we analyzed a 567bp fragment spanning a hyper-variable region of the FSH receptor in the domestic cat (catFSHR) and nine wild felid species/subspecies (felFSHR). Phylogenetic analysis indicated that the newly sequenced felFSHRs, together with the bear FSHR, belong to the carnivore group closely related to the ungulate clade. Within Felidae, genetic distances were 0.0089 +/- 0.0018 for nucleotide and 0.0183 +/- 0.0044 for amino acid (aa) sequences. In pairwise comparisons among catFSHR and all new felFSHRs, similarity ranged from 98.6 to 99.5% for nucleotides and from 97.4 to 98.9% for aa. Besides interspecies variability, intraspecies variation was also detected on both the cDNA and the protein level. There were no indications for an expression of tissue-specific isoforms of FSHR in testis and ovary.

The late Miocene radiation of modern Felidae: a genetic assessment

Modern felid species descend from relatively recent (<11 million years ago) divergence and speciation events that produced successful predatory carnivores worldwide but that have confounded taxonomic classifications. A highly resolved molecular phylogeny with divergence dates for all living cat species, derived from autosomal, X-linked, Y-linked, and mitochondrial gene segments (22,789 base pairs) and 16 fossil calibrations define eight principal lineages produced through at least 10 intercontinental migrations facilitated by sea-level fluctuations. A ghost lineage analysis indicates that available felid fossils underestimate (i.e., unrepresented basal branch length) first occurrence by an average of 76%, revealing a low representation of felid lineages in paleontological remains. The phylogenetic performance of distinct gene classes showed that Y-chromosome segments are appreciably more informative than mitochondrial DNA, X-linked, or autosomal genes in resolving the rapid Felidae species radiation.

Big cat genomics

Advances in population and quantitative genomics, aided by the computational algorithms that employ genetic theory and practice, are now being applied to biological questions that surround free-ranging species not traditionally suitable for genetic enquiry. Here we review how applications of molecular genetic tools have been used to describe the natural history, present status, and future disposition of wild cat species. Insight into phylogenetic hierarchy, demographic contractions, geographic population substructure, behavioral ecology, and infectious diseases have revealed strategies for survival and adaptation of these fascinating predators. Conservation, stabilization, and management of the big cats are important areas that derive benefit from the genome resources expanded and applied to highly successful species, imperiled by an expanding human population.

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.

Phylogenetic assessment of introns and SINEs within the Y chromosome using the cat family felidae as a species tree

The cat family Felidae was used as a species tree to assess the phylogenetic performance of genes, and their embedded SINE elements, within the nonrecombining region of the Y chromosome (NRY). Genomic segments from single-copy X-Y homologs SMCY, UBE1Y, and ZFY (3,604 bp) were amplified in 36 species of cat. These genes are located within the X-degenerate region of the NRY and are thought to be molecular "fossils" that ceased conventional recombination with the X chromosome early within the placental mammal evolution. The pattern and tempo of evolution at these three genes is significant in light of the recent, rapid evolution of the family over approximately 12 Myr and provides exceptional support for each of the eight recognized felid lineages, as well as clear diagnostic substitutions identifying nearly all species. Bootstrap support and Bayesian posterior probabilities are uniformly high for defining each of the eight monophyletic lineages. Further, the preferential use of specific target-site motifs facilitating SINE insertion is empirically supported by sequence analyses of SINEs embedded within the three genes. Target-site insertion is thought to explain the contradiction between intron phylogeny and results of the SMCY SINE phylogeny that unites distantly related species. Overall, our data suggest X-degenerate genes within the NRY are singularly powerful markers and offer a valuable patrilineal perspective in species evolution.