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

Ecological modeling, biogeography, and phenotypic analyses setting the tiger cats' hyperdimensional niches reveal a new species

Recently, the tiger-cat species complex was split into Leopardus tigrinus and Leopardus guttulus, along with other proposed schemes. We performed a detailed analysis integrating ecological modeling, biogeography, and phenotype of the four originally recognized subspecies-tigrinus, oncilla, pardinoides, guttulus-and presented a new multidimensional niche depiction of the species. Species distribution models used > 1400 records from museums and photographs, all checked for species accuracy. Morphological data were obtained from institutional/personal archives. Spotting patterns were established by integrating museum and photographic/camera-trap records. Principal component analysis showed three clearly distinct groups, with the Central American specimens (oncilla) clustering entirely within those of the Andes, namely the pardinoides group of the cloud forests of the southern Central-American and Andean mountain chains (clouded tiger-cat); the tigrinus group of the savannas of the Guiana Shield and central/northeastern Brazil (savanna tiger-cat); and the guttulus group in the lowland forests of the Atlantic Forest domain (Atlantic Forest tiger-cat). This scheme is supported by recent genetic analyses. All species displayed different spotting patterns, with some significant differences in body measurements/proportions. The new distribution presented alarming reductions from the historic range of - 50.4% to - 68.2%. This multidimensional approach revealed a new species of the elusive and threatened tiger-cat complex.

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.

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.

Small spotted bodies with multiple specific mitochondrial DNAs: existence of diverse and differentiated tigrina lineages or species (Leopardus spp: Felidae, Mammalia) throughout Latin America

We analysed two sets of mitochondrial (mt) DNA data from tigrinas (traditionally, Leopardus tigrinus) we sampled in Costa Rica, Venezuela, Colombia, Ecuador, Peru, Bolivia, northwestern and northeastern Argentina and southern Brazil. Additionally, the analysis included some GenBank sequences from southern, central and northeastern Brazil. The first mt set (mt ATP8+mt 16S rRNA with 41 tigrina) revealed the existence of seven different tigrina-like haplogroups. They could represent, at least, 4-6 different tigrina species following the Phylogenetic Species Concept (PSC). In the second mt set (mitogenomics with 18 tigrinas), we detected six different tigrina-like haplogroups. They could represent 4-5 different tigrina species - including a possible full new species, which has gone previously unnoticed to the world of science both morphologic and molecularly. Coat patterns of several of these different tigrinas support the molecular differences. We also detected intense hybridization in many Andean tigrina with margays (Leopardus wiedii) and ocelots (Leopardus pardalis) as well as hybridization of one Bolivian tigrina with Leopardus geoffroyi. Similar hybridization was found for many of the southern Brazilian tigrina (Leopardus guttulus). All of the temporal split estimates for these tigrina haplogroups, together with those of the Leopardus species recognized to date, began in the late Pliocene but mostly occurred during the Pleistocene. In agreement with the existence of multiple species within the traditional L. tigrinus species, we detected strong and significant spatial structure in the two mt data sets. There were clear circular clines. A major part of the analyses detected more genetic resemblance between the Central American + trans Andean Colombian and Ecuadorian tigrina (L. oncilla) with the most geographically distant tigrina from central and southern Brazil (L. guttulus; pure individuals not hybridized with L. geoffroyi). In comparison, the Andean tigrina taxa had intermediate geographical origins but were highly genetically differentiated both from the Central American + trans Andean Colombian-Ecuadorian tigrina and from the central and southern Brazilian tigrina.

Complete mitochondrial genome sequence of Asiatic lion (Panthera leo persica)

The complete mitochondrial genome sequence 17,059 bp of Asiatic lion (Panthera leo persica) has been sequenced with the use of next generation sequencing technology using Ion Torrent ^PGM platform. The complete mitochondrial genome sequence of Asiatic lion consists of 13 protein-coding, 22 tRNA, and two rRNA genes, and 1 control region (CR). The mitochondrial genome is relatively similar to other felid mitochondrial genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases to be typical of those reported for other mammals. The nucleotide composition of Asiatic lion mitogenome shows that there is more A-T% than G-C% on the positive strand as revealed by positive AT and CG skews. The overall base composition is 31.9% of A, 27.2% of C, 14.5% of G, and 26.2% of T. Most of the genes have ATA start codons, except ND1, COX2, ATP8, ATP6, ND4, and ND5 have ATG start codons.

Complete mitochondrial genome sequence of the Indian clouded leopard (Neofelis nebulosa)

The complete mitochondrial genome of sequence 16,859 bp of Indian clouded leopard (Neofelis nebulosa) has been sequenced using next generation sequencing technology Torrent ^PGM platform. The complete mitochondrial genome sequence of clouded leopard consists of 13 protein-coding, 22 tRNA, and two rRNA genes and a control region (CR). The mitochondrial genome is relatively similar to other felid mitochondrial genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases to be typical of those reported for other mammals. The nucleotide composition of the genome shows that there is more A–T% than G–C% on the positive strand as revealed by positive AT and CG skews. The base composition of the mitochondrial genome of clouded leopard is as follows: A, 5362 bp (31.8%); C, 4560bp (27.0%); G, 2475 bp (14.6%); T, 4462 bp (26.4%). Most of the genes have ATG initiation codons, except ND1, ND2, ND3, ND4, ND6, and CYTB (ATA start codon).

Complete mitogenome of Asiatic lion resolves phylogenetic status within Panthera

BACKGROUND

The origin, evolution and speciation of the lion, has been subject of interest, debate and study. The present surviving lions of the genus Panthera comprise of eight sub-species inclusive of Asiatic lion Panthera leo persica of India's Gir forest. Except for the Asiatic lion, the other seven subspecies are found in different parts of Africa. There have been different opinions regarding the phylogenetic status of Panthera leo, as well as classifying lions of different geographic regions into subspecies and races. In the present study, mitogenome sequence of P. leo persica deduced, using Ion Torrent PGM to assess phylogeny and evolution which may play an increasingly important role in conservation biology.

RESULTS

The mtDNA sequence of P. leo persica is 17,057 bp in length with 40.8% GC content. Annotation of mitogenome revealed total 37 genes, including 13 protein coding, 2 rRNA and 22 tRNA. Phylogenetic analysis based on whole mitogenome, suggests Panthera pardus as a neighbouring species to P. leo with species divergence at ~2.96 mya.

CONCLUSION

This work presents first report on complete mitogenome of Panthera leo persica. It sheds light on the phylogenetic and evolutionary status within and across Felidae members. The result compared and evaluated with earlier reports of Felidae shows alteration of phylogenetic status and species evolution. This study may provide information on genetic diversity and population stability.

[Summary of phylogeny in family Felidae of Carnivora]

Felidae (cats) is one of the strict carnivorous groups in the order Carnivora, many of which are most familiar and spectacular to us. They are the top predators in the world. Thirty-six of 37 living cat species are considered as either "endangered" or "threatened". The relationships among species of the family Felidae, which evolved recently and rapidly, are difficult to resolve, and have been the subject of debate. Construction of a reliable Felidae phylogeny will be of evolutionarily significance and conservation value. In this paper, we summarized phylogeny of Felidae, including cytological, morphological and molecular evidence, and pointed out the existing phylogenetic problems. This review is expected to guide future researches of Felidae phylogeny, and to lay a theoretic foundation for the protection of this animal group.

Leopardus braccatus (Carnivora: Felidae)

Leopardus braccatus (Cope, 1889) is a small felid—not much larger than a domestic house cat—commonly called the Pantanal cat. No comprehensive surveys have been conducted to determine how many L. braccatus exist in the wild. It is found in humid, warm grasslands and wooded areas of extreme northwestern Argentina, southwestern and south- and north-central (newly reported ranges) Brazil, Paraguay, and Uruguay. Habitat loss and degradation are considered major threats throughout most of its range. It is uncommon in captivity and museum collections, listed with all Felidae under Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora, and considered “Near Threatened” by the International Union for Conservation of Nature and Natural Resources in the pampas cat group ( = L. colocolo).

Big cat phylogenies, consensus trees, and computational thinking

Phylogenetics seeks to deduce the pattern of relatedness between organisms by using a phylogeny or evolutionary tree. For a given set of organisms or taxa, there may be many evolutionary trees depicting how these organisms evolved from a common ancestor. As a result, consensus trees are a popular approach for summarizing the shared evolutionary relationships in a group of trees. We examine these consensus techniques by studying how the pantherine lineage of cats (clouded leopard, jaguar, leopard, lion, snow leopard, and tiger) evolved, which is hotly debated. While there are many phylogenetic resources that describe consensus trees, there is very little information, written for biologists, regarding the underlying computational techniques for building them. The pantherine cats provide us with a small, relevant example to explore the computational techniques (such as sorting numbers, hashing functions, and traversing trees) for constructing consensus trees. Our hope is that life scientists enjoy peeking under the computational hood of consensus tree construction and share their positive experiences with others in their community.

Highly conserved D-loop-like nuclear mitochondrial sequences (Numts) in tiger (Panthera tigris)

Using oligonucleotide primers designed to match hypervariable segments I (HVS-1) of Panthera tigris mitochondrial DNA (mtDNA), we amplified two different PCR products (500 bp and 287 bp) in the tiger (Panthera tigris), but got only one PCR product (287 bp) in the leopard (Panthera pardus). Sequence analyses indicated that the sequence of 287 bp was a D-loop-like nuclear mitochondrial sequence (Numts), indicating a nuclear transfer that occurred approximately 4.8-17 million years ago in the tiger and 4.6-16 million years ago in the leopard. Although the mtDNA D-loop sequence has a rapid rate of evolution, the 287-bp Numts are highly conserved; they are nearly identical in tiger subspecies and only 1.742% different between tiger and leopard. Thus, such sequences represent molecular 'fossils' that can shed light on evolution of the mitochondrial genome and may be the most appropriate outgroup for phylogenetic analysis. This is also proved by comparing the phylogenetic trees reconstructed using the D-loop sequence of snow leopard and the 287-bp Numts as outgroup.

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.

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.

Phylogeography and genetic ancestry of tigers (Panthera tigris)

Eight traditional subspecies of tiger (Panthera tigris), of which three recently became extinct, are commonly recognized on the basis of geographic isolation and morphological characteristics. To investigate the species' evolutionary history and to establish objective methods for subspecies recognition, voucher specimens of blood, skin, hair, and/or skin biopsies from 134 tigers with verified geographic origins or heritage across the whole distribution range were examined for three molecular markers: (1) 4.0 kb of mitochondrial DNA (mtDNA) sequence; (2) allele variation in the nuclear major histocompatibility complex class II DRB gene; and (3) composite nuclear microsatellite genotypes based on 30 loci. Relatively low genetic variation with mtDNA, DRB, and microsatellite loci was found, but significant population subdivision was nonetheless apparent among five living subspecies. In addition, a distinct partition of the Indochinese subspecies P. t. corbetti into northern Indochinese and Malayan Peninsula populations was discovered. Population genetic structure would suggest recognition of six taxonomic units or subspecies: (1) Amur tiger P. t. altaica; (2) northern Indochinese tiger P. t. corbetti; (3) South China tiger P. t. amoyensis; (4) Malayan tiger P. t. jacksoni, named for the tiger conservationist Peter Jackson; (5) Sumatran tiger P. t. sumatrae; and (6) Bengal tiger P. t. tigris. The proposed South China tiger lineage is tentative due to limited sampling. The age of the most recent common ancestor for tiger mtDNA was estimated to be 72,000–108,000 y, relatively younger than some other Panthera species. A combination of population expansions, reduced gene flow, and genetic drift following the last genetic diminution, and the recent anthropogenic range contraction, have led to the distinct genetic partitions. These results provide an explicit basis for subspecies recognition and will lead to the improved management and conservation of these recently isolated but distinct geographic populations of tigers.

Genetic analysis provides the basis for subspecies recognition among tigers, and will lead to improved conservation strategies for these endangered animals

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.

Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats

Mitochondrial DNA control region structure and variation were determined in the five species of the genus Panthera. Comparative analyses revealed two hypervariable segments, a central conserved region, and the occurrence of size and sequence heteroplasmy. As observed in the domestic cat, but not commonly seen in other animals, two repetitive sequence arrays (RS-2 with an 80-bp motif and RS-3 with a 6-10-bp motif) were identified. The 3' ends of RS-2 and RS-3 were highly conserved among species, suggesting that these motifs have different functional constraints. Control region sequences provided improved phylogenetic resolution grouping the sister taxa lion (Panthera leo) and leopard (Panthera pardus), with the jaguar (Panthera onca).

Phylogeography, population history and conservation genetics of jaguars (Panthera onca, Mammalia, Felidae)

The jaguar (Panthera onca), the largest felid in the American Continent, is currently threatened by habitat loss, fragmentation and human persecution. We have investigated the genetic diversity, population structure and demographic history of jaguars across their geographical range by analysing 715 base pairs of the mitochondrial DNA (mtDNA) control region and 29 microsatellite loci in approximately 40 individuals sampled from Mexico to southern Brazil. Jaguars display low to moderate levels of mtDNA diversity and medium to high levels of microsatellite size variation, and show evidence of a recent demographic expansion. We estimate that extant jaguar mtDNA lineages arose 280 000-510 000 years ago (95% CI 137 000-830 000 years ago), a younger date than suggested by available fossil data. No strong geographical structure was observed, in contrast to previously proposed subspecific partitions. However, major geographical barriers such as the Amazon river and the Darien straits between northern South America and Central America appear to have restricted historical gene flow in this species, producing measurable genetic differentiation. Jaguars could be divided into four incompletely isolated phylogeographic groups, and further sampling may reveal a finer pattern of subdivision or isolation by distance on a regional level. Operational conservation units for this species can be defined on a biome or ecosystem scale, but should take into account the historical barriers to dispersal identified here. Conservation strategies for jaguars should aim to maintain high levels of gene flow over broad geographical areas, possibly through active management of disconnected populations on a regional scale.

Novel gene conversion between X-Y homologues located in the nonrecombining region of the Y chromosome in Felidae (Mammalia)

Genes located on the mammalian Y chromosome outside of the pseudoautosomal region do not recombine with those on the X and are predicted to either undergo selection for male function or gradually degenerate because of an accumulation of deleterious mutations. Here, phylogenetic analyses of X-Y homologues, Zfx and Zfy, among 26 felid species indicate two ancestral episodes of directed genetic exchange (ectopic gene conversion) from X to Y: once during the evolution of pallas cat and once in a common predecessor of ocelot lineage species. Replacement of the more rapidly evolving Y homologue with the evolutionarily constrained X copy may represent a mechanism for adaptive editing of functional genes on the nonrecombining region of the mammalian Y chromosome.

Disparate phylogeographic patterns of molecular genetic variation in four closely related South American small cat species

Tissue specimens from four species of Neotropical small cats (Oncifelis geoffroyi, N = 38; O. guigna, N = 6; Leopardus tigrinus, N = 32; Lynchailurus colocolo, N = 22) collected from throughout their distribution were examined for patterns of DNA sequence variation using three mitochondrial genes, 16S rRNA, ATP8, and NADH-5. Patterns between and among O. guigna and O. geoffroyi individuals were assessed further from size variation at 20 microsatellite loci. Phylogenetic analyses using mitochondrial DNA sequences revealed monophyletic clustering of the four species, plus evidence of natural hybridization between L. tigrinus and L. colocolo in areas of range overlap and discrete population subdivisions reflecting geographical isolation. Several commonly accepted subspecies partitions were affirmed for L. colocolo, but not for O. geoffroyi. The lack of geographical substructure in O. geoffroyi was recapitulated with the microsatellite data, as was the monophyletic clustering of O. guigna and O. geoffroyi individuals. L. tigrinus forms two phylogeographic clusters which correspond to L.t. oncilla (from Costa Rica) and L.t. guttula (from Brazil) and which have mitochondrial DNA (mtDNA) genetic distance estimates comparable to interspecific values between other ocelot lineage species. Using feline-specific calibration rates for mitochondrial DNA mutation rates, we estimated that extant lineages of O. guigna diverged 0.4 million years ago (Ma), compared with 1.7 Ma for L. colocolo, 2.0 Ma for O. geoffroyi, and 3.7 Ma for L. tigrinus.

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.

Phylogenetic reconstruction of the Felidae using 16S rRNA and NADH-5 mitochondrial genes

The Felidae family represents a challenge for molecular phylogenetic reconstruction because it consists of 38 living species that evolved from a relatively recent common ancestor (10-15 million years ago). We have determined mitochondrial DNA sequences from two genes that evolve at relatively rapid evolutionary rates, 16S rRNA (379 bp) and NADH dehydrogenase subunit 5 (NADH-5, 318 bp), from multiple individuals of 35 species. Based on separate and combined gene analyses using minimum evolution, maximum parsimony, and maximum likelihood phylogenetic methods, we recognized eight significant clusters or species clades that likely reflect separate monophyletic evolutionary radiations in the history of this family. The clusters include (1) ocelot lineage, (2) domestic cat lineage, (3) Panthera genus, (4) puma group, (5) Lynx genus, (6) Asian leopard cat group, (7) caracal group, and (8) bay cat group. The results confirm and extend previously hypothesized associations in most cases, but in others, e.g., the bay cat group, suggest novel phylogenetic relationships. The results are compared and evaluated with molecular, cytogenetic, and morphological data to derive a phylogenetic synthesis of field evolutionary history.