Email updates

Keep up to date with the latest news and content from Frontiers in Zoology and BioMed Central.

Open Access Highly Accessed Research

Clouded leopard phylogeny revisited: support for species recognition and population division between Borneo and Sumatra

Andreas Wilting1*, Valerie A Buckley-Beason2, Heike Feldhaar3, Jürgen Gadau4, Stephen J O'Brien2 and K Eduard Linsenmair1

  • * Corresponding author: Andreas Wilting a.wilting@gmx.de

  • † Equal contributors

Author Affiliations

1 University of Würzburg, Biocentre, Department of Animal Ecology and Tropical Biology, 97074 Würzburg, Germany

2 Laboratory of Genomic Diversity, National Cancer Institute, Frederick, Maryland, USA

3 University of Würzburg, Biocentre, Department of Behavioural Physiology and Sociobiology, 97074 Würzburg, Germany

4 School of Life Science, Arizona State University, Tempe, AZ 85287, USA

For all author emails, please log on.

Frontiers in Zoology 2007, 4:15  doi:10.1186/1742-9994-4-15


The electronic version of this article is the complete one and can be found online at: http://www.frontiersinzoology.com/content/4/1/15


Received:18 March 2007
Accepted:29 May 2007
Published:29 May 2007

© 2007 Wilting et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The clouded leopard (Neofelis nebulosa) is one of the least known cat species and depletion of their forested habitats puts it under heavy pressure. Recently reclassification of Bornean clouded leopards (N. nebulosa diardi) to species level (N.diardi) was suggested based on molecular and morphological evidence. Since the genetic results were based solely on three Bornean samples we re-evaluated this partition using additional samples of Bornean clouded leopards (N = 7) and we were also able to include specimens from Sumatra (N = 3), which were lacking in previous analysis.

Results

We found strong support for the distinction between N. nebulosa and N. diardi based on three fragments of mtDNA (900 bp) and 18 microsatellites. Forty-one fixed mitochondrial nucleotide differences and non-overlapping allele sizes in 8 of 18 microsatellite loci distinguished N. nebulosa and N. diardi. This is equivalent to the genetic divergence among recognized species in the genus Panthera. Sumatran clouded leopards clustered with specimens from Borneo, suggesting that Sumatran individuals also belong to N. diardi. Additionally, a significant population subdivision was apparent among N. diardi from Sumatra and Borneo based on mtDNA and microsatellite data.

Conclusion

Referring to their origin on two Sunda Islands we propose to give N. diardi the common name "Sundaland clouded leopard". The reduced gene flow between Borneo and Sumatra might suggest the recognition of two subspecies of N. diardi. Based on this reclassification of clouded leopards not only species, but also the populations on Borneo and Sumatra should be managed separately and a higher priority should be placed to protect the different populations from extinction.

Background

Clouded leopards have drawn less attention in the past compared to other species of the Panthera lineage. This medium-sized cat (11–25 kg) primarily inhabits densely vegetated habitats and remote areas in South-East Asian subtropical and tropical forests [1,2]. As a consequence of their secretive behaviour [3,4] very little is known about their behaviour and status in the wild and very few studies have focused on this species [5-7]. Likewise, their taxonomy and phylogenetic status remained obscure for a long time [8-11]. Phylogenetic studies consistently reveal that clouded leopards separated from the Panthera lineage approximately 6 million years ago (MYA) [12-14]. Only recently the phylogenetic relationships among clouded leopards were assessed [15]. Based on mtDNA, nuclear DNA sequences, microsatellites, fixed chromosomal differences and morphological variation, the reclassification of Bornean clouded leopards (N. nebulosa diardi) to species level N. diardi was suggested [15,16]. The small sample size of only three Bornean samples in the phylogenetic study reduces confidence in this reclassification of N. diardi. Therefore, we present here the phylogenetic analysis with additional individuals from Borneo. Furthermore, we examine the phylogenetic relationships among clouded leopards more closely by including Sumatran samples. Sumatran individuals are necessary for the recognition of a new species on the Sunda Islands, but until now have not been investigated genetically. The goal of this study is to provide a better understanding of the systematic classification of one of the most threatened cat species in Asia. This classification is of utmost importance for conservation and management purposes [17].

Results

Mitochondrial DNA analysis

Analysis of the phylogenetic relationship among clouded leopards based on a concatenated fragment of 900 bp of mitochondrial DNA comprising sequences from three mitochondrial gene segments (185 bp ATPase-8, 286 bp Cyt b and 426 bp control region) revealed 54 variable sites among clouded leopards. Of those 39 were fixed nucleotide differences between specimens from the islands of Borneo/Sumatra (N. diardi) and specimens from the mainland (N. nebulosa nebulosa). There were 12 nucleotide differences in the Cyt-b gene (4.2 %), 10 in the ATPase-8 gene (5.4 %) and 17 in the control region (3.9 %). In comparison, three Panthera species (P. pardus, P.tigris and P. onca) were separated by 38–52 nucleotide differences in the same fragment. Four fixed nucleotide differences distinguished the individuals from Sumatra and Borneo. Altogether we found 13 haplotypes among all the clouded leopards examined. N. nebulosa nebulosa (N = 59 GenBank accession no. in Table 2. see Figure 1) had six haplotypes (NEB 1–6), N. diardi (Borneo, N = 7) had 5 haplotypes (DIB 1–5) and in N. diardi (Sumatra, N = 3) we found 2 haplotypes (DIS 1 & 2).

Table 1. Specimens sampled in this study

thumbnailFigure 1. Table 2: Haplotypes and variable sites in combined analysis of 900 bp of mtDNA sequences.

Phylogenetic analysis of mtDNA haplotypes using minimum evolution (ME), maximum parsimony (MP) and maximum likelihood (ML) approaches produced congruent topologies that strongly support the reciprocally monophyletic status of N.nebulosa and N. diardi with high bootstrap values (100 % ME/MP, 98 % ML) (Figure 2A). Furthermore, Sumatran individuals were separated from Bornean clouded leopards supported by high bootstrap values; except for the ME approach (46% ME, 85 % MP, and 81 % ML).

thumbnailFigure 2. Phylogenetic relationships among clouded leopards from mtDNA haplotypes and a Bayesian analysis of microsatellite genotypes. (A) Phylogenetic relationships based on minimum evolution (ME) among the clouded leopard mtDNA haplotypes from the concatenated 900 bp mitochondrial sequences comprising Cyt-b (286 bp), ATPase-8 (186 bp) and control region (426 bp) gene fragments. Panthera samples and Felis catus samples were taken as outgroups. Trees constructed with PAUP [56] obtained under maximum parsimony (MP) and maximum likelihood (ML) criteria have identical topologies. Numbers above the branches represent bootstrap support (1000 replicates) for each three methods (ME/MP/ML); only those with > 80 % are shown. Numbers in parentheses represent the number of individuals sharing the same haplotype. We used Kimura 2-parameter distance with neighbor-joining (NJ) algorithm followed by tree-bisection reconnection branch-swapping procedure (TBR) for the ME analysis. MP trees were constructed using a heuristic search, with a random addition of taxa and TBR branch swapping. The ML approach was performed using a HKY85 model [60]. Haplotype codes are shown in Table 1. NEB 1 – 5, DIB 1 and 2, and outgroups have been described previously [15]. * Sequences of only two mtDNA genes (ATPase-8 and Cyt b) were included. (B) Greyscale bars are from the Bayesian admixture analysis of the microsatellite analysis. Different colours stand for different genetic groups. ID codes in are shown in Table 1.

A nuclear pseudogene insertion of cytoplasmic mtDNA often referred as numt have been reported in numerous organisms [reviewed in [18]] including the Panthera genus and the domestic cat [19-21]. This can complicate the analysis of mtDNA, because numt can coamplify with the mitochondrial genes. We cannot exclude the presence of numts in clouded leopards. However, it is very unlikely that we sequenced numt instead of mtDNA, because all our protein coding sequences had open reading frames [22]. Furthermore other nuclear genes did not amplify with our ancient DNA, even if we tried to sequence very short fragments. Therefore, we conclude that most likely we did not encounter any evidence of numt.

Microsatellite analysis

Composite genotypes from 18 felid-specific microsatellite loci [23,24] were obtained from 11 clouded leopard samples (four mainland specimens NNE 1 – 4, four Bornean specimens NDB 1 – 4 and three Sumatran specimens NDS 1 – 3), two domestic cats (Felis catus; FCA 1 – 2) and one leopard (PPA 1). Of the 10 loci, which were used before in clouded leopards, Buckley-Beason et al. [15] showed that six microsatellite loci (FCA 82, FCA 105, FCA 132, FCA 144, FCA 261, FCA 310) had non-overlapping allele sizes between N. nebulosa (mainland) and N. diardi (Borneo). For only one of these six loci, FCA 261, non-overlapping allele-sizes could not be confirmed by wider sampling of Bornean and Sumatran specimens. Out of eight microsatellite loci, which were not tested in clouded leopards before, three (FCA 23, FCA 43 and HDZ 859) did not overlap in allele sizes. Overall expected heterozygosity HE in clouded leopards ranged from 0.488 in the Bornean population to 0.652 in N. nebulosa and exceeds the observed heterozygosity H0 in all three populations, with Bornean specimens having the lowest observed heterozygosity H0 with 0.361 (Table 3).

Table 3. Measure of genetic variance among 18 microsatellite loci in clouded leopards.

Neighbor-joining analysis of individual clouded leopard genotypes based on two microsatellite genetic distance estimators (Dps & Dkf) produced concordant topologies; both trees support the species distinction among clouded leopards (Figure 3). Individuals from Borneo and Sumatra form a monophyletic clade with 100 % (Dps & Dkf) bootstrap support (BS) distinguishing them from mainland specimens and the outgroups. The microsatellite analysis lends further support to the phylogeographic subdivision observed in mtDNA analysis between Borneo and Sumatra individuals, however with lower bootstrap support values (Sumatra clade with 69 % Dps/61%Dkf BS, and Borneo clade with 50/46 % BS).

thumbnailFigure 3. Phylogenetic relationships among the individual clouded leopards from composite microsatellite genotypes of 18 loci. One Panthera pardus sample and two Felis catus samples were included as outgroups. Branches of the same greyscale represent individuals from the same geographical region. Trees are based on the proportion of shared alleles (Dps) and kinship coefficient (Dkf) genetic distances with 1 – (kf/ps) option in MICROSAT [65] produced identical topologies. Dps tree is shown here. Bootstrap values over 70 % are shown on the divergence node (Dps/Dkf). ID codes are shown in Table 1.

Population substructures

To evaluate population distinctiveness we tested the microsatellite data using a Bayesian algorithm as implemented in BAPS version 4.14 [25]. To estimate hidden substructures, within clouded leopards BAPS suggested to partition clouded leopard into five populations (p > 0.99) three different mainland populations, Sumatra and Borneo. However mtDNA results and previous analysis based on a broader sampling, including more microsatellite loci found no substructure among mainland populations [15]. Therefore we assume that the population partition within our four mainland individuals is a result of the small sampling size and thus set a maximum of three populations. BAPS then grouped all mainland, Sumatra, and Borneo individuals together, respectively (p = 1). We confirmed the individual assignments to each population with admixture analysis in BAPS [26,27]. In this scenario all individuals were assigned to three unique clusters N. nebulosa,N. diardi (Borneo) and N. diardi (Sumatra) with very high admixture coefficients (all q > 0.98, except for one Sumatran individual (NDS1) q > 0.95) (Figure 2B). No individual had a Bayesian p-value less than 0.05 (all p > 0.91) revealing no evidence of admixed background and none or very low degrees of past and present gene flow between all three populations.

MtDNA sequence differences were used to estimate the divergence time of N.nebulosa and N. diardi and the time of origin of Sumatran clouded leopards. Using a calibration of 6.37 MYA for the divergence of clouded leopards from the Panthera lineage [14], N. diardi diverged from N. nebulosa about 2.86 MYA (95 % CI of 1.71–4.02 MYA) and Sumatran and Bornean clouded leopards diverged about 437,000 years ago (95 % CI of 30,000 – 845,000 years ago).

Discussion

Our results strongly support the reclassification of extant clouded leopards into two distinct species N. nebulosa and N. diardi. Based on mtDNA data clouded leopards on the islands of Sumatra and Borneo have been reproductively isolated from the mainland species since middle to late Pliocene (~2.86 MYA). We are aware that recent studies showed that mtDNA is least robust in node resolution [14], which might lead to an overestimate of our calculated length of time. However previously estimated divergence time of 1.41 MYA between N. nebulosa and N. diardi [15] is still within the same range (1 – 3 MYA) as species level distinctions across Panthera [14]. The wider sampling of Bornean individuals dispels the doubt that the high distinction between N. nebulosa and N. diardi described before was a consequence of an inadequate sample size [15]. Referring to their origin on two Sunda Islands, we would propose to give N. diardi the common name "Sundaland clouded leopard".

Furthermore, the consistent results of mtDNA and microsatellite data provide an evidence for a reduced gene flow between Sumatra and Borneo. The trend in the reduction of observed heterozygosity might indicate a Wahlund effect, within Neofelis, as well as within N. diardi. Due to the low number of samples from each population within the sampling this result should be considered as preliminary, although it supports a hypothesised disconnection of gene flow between clouded leopard populations on the islands of Borneo and Sumatra. On the basis of our genetic results we recommend the recognition of two distinct subspecies of N. diardi (Figure 4), considering previous criteria for the designation of subspecies [28,29]. We estimate that since the middle to late Pleistocene Bornean and Sumatran clouded leopards were most likely isolated from each other and unable to move freely between islands. However, it has to be considered that our sample size of Sumatran specimens is small and further molecular genetic and morphological studies will be needed to confirm our findings. If these studies will support our data, we suggest giving clouded leopards on Borneo and Sumatra the scientific names N. diardi borneensis and N.diardi sumatrensis, referring to their origin.

thumbnailFigure 4. Sampling locations and suggested new classification of clouded leopards. The new classification is based on this molecular analysis and data obtained from Buckley-Beason et al. [15]. Numbers before the slash indicate the number of samples from this study, those after the slash indicate the number of samples that were included in the mtDNA analysis obtained from Buckley-Beason et al. [15].

Our results have several implications for the conservation and management strategies for modern clouded leopards. N. nebulosa and N. diardi should be managed separately, and treated as different species as suggested before [15,16]. Furthermore, on the basis of the population division between Bornean and Sumatran clouded leopards we suggest that these populations should also be managed separately. Considering previous criteria for the designation of evolutionarily significant units (ESU) [30-35], Bornean and Sumatran clouded leopards should be treated as different conservation units with separate management plans. Both populations are on the one hand reproductively isolated from each other and monophyletic for mtDNA and microsatellites and on the other hand Bornean and Sumatran individuals represent an important component of the evolutionary legacy of the species. To our knowledge there have been no translocations between clouded leopards from Borneo and Sumatra so far.

The continued depletion of tropical rainforests and fragmentation of natural habitats in Borneo and Sumatra put the reclassified species N diardi under severe pressure of extinction [7,36]. Therefore a higher priority should be placed on effective conservation of Sundaland clouded leopards and their shrinking habitats. Further research is urgently needed to reveal the distribution and status of different species in situ, because smaller distribution ranges associated with reduced gene pools of the reclassified species puts clouded leopards under a greater risk of extinction. The current IUCN Red Data category "vulnerable" [37] might underestimate the threat these cats are facing. To date only for Sabah (north-eastern part of Borneo) preliminary rough population estimates exist, demonstrating the severe threat of Bornean clouded leopards [7].

The inclusion of Sumatran samples reveals a very interesting aspect of clouded leopards' evolutionary history. At present, shallow seas separate Sumatra, Borneo, Java and Malay Peninsula. These areas were connected several times via land bridges during periods of low sea levels in the late Pliocene and Pleistocene [e. g. [38-40]]. Therefore many authors have characterized Sundaland as a geographical unit across which species should have been able to move freely during glacial periods until 10,000 years ago when higher sea-levels started to separate the islands [41-44]. Based on the presented data we cannot support the hypothesis of one geographical unit for the clouded leopard, because Bornean and Sumatran animals were reproductively isolated from the mainland individuals even during glaciation periods with accompanying low sea levels and postulated land bridges. In contrast, clouded leopards were only able to move throughout the exposed shelf between Borneo and Sumatra at least once during the early or middle Pleistocene. Therefore, despite the late Pleistocene existence of land connections [39,40], present-day distribution patterns, exemplarily shown here for the clouded leopard, indicate that dispersal was restricted and there appear to have been considerable barriers for animal migrations [40,45,46]. The wider sampling of this study in addition to previous studies [15,16] have provided a better insight into phylogeographic history of one of the least known cat species in South-East Asia.

Conclusion

The results we present in this paper are of importance for the understanding of clouded leopard phylogeny and for their in situ and ex situ conservation and management. The different species N. nebulosa and N. diardi, as well as the two distinct populations on Borneo and Sumatra, should be managed separately, to protect the genetic diversity upon which future evolutionary potential depends. Our paper gives a good example of the importance of taxonomic splitting for conservation, because the two species and the distinct populations on Borneo and Sumatra face a much greater risk of extinction, due to smaller distribution ranges, than previously assessed based on the former classification.

The wider genetic sampling contributed to the understanding of the evolutionary history of clouded leopards. The long isolation of N. diardi and N. nebulosa revealed that the clouded leopard, comparable to other forest-dwelling species [47-51], had a deep history of vicariant evolution. Our results show that geographical or ecological barriers must have existed during Pleistocene glaciation periods hindering migrations of clouded leopards. This raises the question about the phylogeographic history of other taxa in the Sunda shelf and further research will be needed for a better understanding of evolutionary processes in this region

Methods

Samples and DNA extraction

We sampled four wild-born specimens from the islands of Borneo, three ancient samples from Sumatra and four recent animals from known geographic origin on the mainland (Table 1) in addition to 58 clouded leopard individuals described previously [15]. Panthera species and domestic cat individuals as outgroups for mitochondrial DNA (mtDNA) analysis were specified before [15] and for outgroup comparison of the microsatellite analysis one leopard (Panthera pardus) and two domestic cats (Felis catus) were sampled (Table 1).

DNA was extracted from fecal samples using QiAmp Stool Mini Kit (Qiagen, Hilden, Germany) and from serum and whole blood using QiAmp DNeasy (Qiagen, Hilden, Germany). For the extraction of DNA from historical museum samples hide and dry tissue samples were cut into small pieces using a sterile scalpel followed by a standard proteinase K digestion with an extended incubation interval at 56°C for up to 72 hours [40] until only little solid tissue remained. A standard phenol/chloroform extraction procedure was then used for DNA extraction [52]. DNA was suspended in 50 μl of double distilled H2O.

Mitochondrial DNA analysis

We used a 426 bp portion of a central conserved region within the D-loop of the control region [53] in addition to two mtDNA (ATPase-8 and Cyt-b) genes [15]. The control region primers were modified by Janecka JE (pers. comm.) [(F) CTC AAC TAT CCG AAA GAG CTT] and [(R) CCT GTG GAA CAT TAG GAA TT]. In total we amplified 900 bp of mtDNA.

PCR reactions were performed in a final volume of 25 μl containing 2.5 μl MolTaq 10x PCR Buffer, 2 mM MgCl2, 0.2 mM dNTPs, 1 μM of each primer, 2 units of MolTaq polymerase (Molzym GmbH, Bremen, Germany) and 2 μl of genomic DNA. PCR reactions were performed in an Eppendorf Mastercycler (Eppendorf GmbH, Wesseling-Berzdorf, Germany), with an initial denaturation step at 95°C for 3 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 45 s, elongation at 72°C for 45 s, and were completed with a final elongation step at 72°C for 10 min. PCR products were purified by ultra filtration through Montage TM filter devices (Millipore GmbH, Schwalbach, Germany) and sent to Seqlab (Seqlab Laboratories, Göttingen, Germany) for sequencing. Sequences were edited, assembled and aligned using ClustalW [54] implemented in BioEdit (Version 7.0.5.2) [55], before being exported to PAUP (Version 4.0b10) [56] for phylogenetic analysis. Sequences from each of the mtDNA fragments were concatenated into a 900 bp sequences, because results from separate analysis of each gene fragment showed identical topologies [57] and mitochondrial genes usually do not recombine [58]. Phylogenetic relationships among haplotypes were estimated using minimum evolution (ME), maximum likelihood (ML) and maximum parsimony (MP) [56,59]. We used Kimura 2-parameter distance with neighbor-joining (NJ) algorithm followed by tree-bisection reconnection branch-swapping procedure (TBR) for the ME analysis. MP trees were conducted using a heuristic search, with 10 random taxon addition replicates and TBR branch swapping. The ML approach was performed using the HKY85 model [60]. Each phylogenetic tree was rooted with the domestic cat sequence. Reliability of all trees was tested with bootstrap values by 1000 replicates of heuristic search and TBR branch swapping.

The approximate age of separation between N. nebulosa and N. diardi was estimated using LINTREE [61]. A neighbor-joining tree [59] was generated with Kimura 2-parameter γ-corrected distances [62] using the combined 900 bp mtDNA sequence. The molecular clock test implemented in LINTREE [61] showed that the sequences did not deviate significantly from the rate constancy test (p > 0.05). The coalescence point between clouded leopards and the Panthera genus, being 6.37 MYA based upon a comprehensive analysis of nuclear gene sequences and multiple fossil dates [14], was chosen to be the calibration point for this study. We used a range of two standard errors to calculate a 95 % confidence interval.

Microsatellite markers

We used 10 felid dinucleotide microsatellite primers (FCA 8, FCA 45, FCA 77, FCA82, FCA 105, FCA 126, FCA 132, FCA 144, FCA 261, FCA 310) [23], which were already used in a previous study on clouded leopards [15]. Those microsatellites are located on 8 felid autosomes and all of them are at least 5 centimorgans apart from each other [23]. In addition to those ten microsatellite loci of known allele size ranges for clouded leopards, we applied 8 microsatellites of unknown allele sizes for clouded leopards, FCA 23 and FCA 43 [23], and HDZ 3, HDZ 57, HDZ 64, HDZ 89, HDZ817, HDZ 859 [24]. PCR amplifications were performed in a final reaction volume of 10 μl utilizing described methods [23,24]. The IR-dye-labeled PCR products were diluted and analyzed on a LI-COR 4300 DNA-Analyser (LI-COR Bioscience GmbH, Bad Homburg, Germany). Data were collected and analyzed using Saga Generation 2 (Version 3.2.1).

To test their performance as population genetic markers, all microsatellites were tested for deviations from linkage disequilibrium (LD) using GENEPOP on the web version 3.4 [63]. Measures of microsatellite genetic variation in terms of observed and expected heterozygosities were estimated with Arlequin 3.1 [64].

Pairwise genetic distance among clouded leopards and the two outgroup species was estimated with two microsatellite genetic distance estimators: the proportion of shared alleles (Dps) and the kinship coefficient (Dkf) with the [1 – ps/kf] option in MICROSAT [65]. Phylogenetic NJ-trees were constructed from the Dps and Dkf distance matrixes using NEIGHBOR (included in PHYLIP version 3.66) [66]. Bootstrap values for 1000 bootstrap replicates in MICROSAT were calculated using CONSENSE TREE (included in PHYLIP version 3.66) [66]. Trees were drawn using the program TREEVIEW (version 1.6.6) [67].

Population structure analysis

A Bayesian clustering method as implemented in BAPS [25-27] was used to infer population structure based on multilocus microsatellite genotype data. This program estimates the hidden population substructure by testing whether the allele frequencies between populations are significantly different. A major advantage compared to most other methods is that the number of populations is treated here as an unknown parameter that can be estimated from the dataset. We performed 10 independent runs of clustering of individuals with the microsatellite genotypes to ensure homogenous results. In all ten runs we obtained similar results (data not shown). After the clustering of individuals by their allele frequencies the results were used to perform an admixture analysis. We used 500 iterations and a number of 1000 reference individuals per population each with 20 iterations. The estimated admixture coefficient for an individual in each cluster q (maximum = 1) was used as a measure of correct assignments. The Bayesian p-value tells the proportion of reference individuals simulated from the population in which the individual was originally clustered having the admixture coefficient to the cluster smaller than or equal to the individual [27]. Individuals having p-values larger than 0.05 are by default considered as having "non-significant" evidence for admixture [27].

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

AW collected the samples for this study and carried out the molecular work at the University of Würzburg including genotyping of microsatellite markers, and mtDNA sequencing. AW also performed all sequence and microsatellite analysis and drafted the manuscript. VABB carried out the molecular work at the National Cancer Institute as well as gave very valuable comments on the manuscript. HF mentored the study at the University of Würzburg and helped to analyse the data. HF also has contributed to the preparation of the manuscript. JG participated in the study design. SJO and KEL coordinated the study at the National Cancer Institute and at the University of Würzburg, respectively, and contributed to the final revisions of the manuscript.

All authors read and approved the final manuscript.

Acknowledgements

We thank D. Hesse for comments on earlier drafts of this manuscript. Special thanks also go to all institutions listed in Table 1 that supplied the biological specimens this work is based upon. We would also like to thank CITES agencies for issuing the specific permits to the University of Würzburg and the National Cancer Institute. Point Defiance Zoo and Aquarium and Duisburg Zoo provided financial support for data collection.

References

  1. Brehm AE: Brehm's illustrirtes Thierleben – Eine allgemeine Kunde des Thierreichs, 2. Auflage. Hildburghausen: Bibliographisches Institut; 1876:408-409. OpenURL

  2. Corbett GB, Hill JE: Family Felidae – Cats. In The Mammals of the Indo-Malayan Region: A systematic review. Oxford University Press for British Museum (Natural History), London; 1992:219-228. OpenURL

  3. Banks E: A popular account of the mammals of Borneo.

    J Malayan Brit Roy Asiatic Soc 1931, 9:1-139. OpenURL

  4. Pocock RI: Family Felidae. In Fauna of British India (Mammals). Taylor and Francis, London; 1939:190-331. OpenURL

  5. Austin SC, Tewes ME: Ecology of the clouded leopard in Khao Yai National Park, Thailand.

    Cat News/IUCN SSC 1999, 31:17-18. OpenURL

  6. Grassman LI Jr, Tewes ME, Silvy NJ, Kreetiyutanont K: Ecology of three sympatric felids in a mixed evergreen forest in North-Central Thailand.

    J Mamm 2005, 86:29-38. Publisher Full Text OpenURL

  7. Wilting A, Fischer F, Abu Bakar S, Linsenmair KE: Clouded leopards, the secretive top-carnivore of South-East Asian rainforests: their distribution, status and conservation needs in Sabah, Malaysia.

    BMC Ecol 2006, 6:16. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  8. Pocock RI: The classification of existing Felidae.

    Annals and Magazine of Natural History 1917, 20:329-350. OpenURL

  9. Haltenorth T: Die verwandtschaftliche Stellung der Großkatzen zueinander.

    Zeitschrift fuer Saeugetierkunde 1936, 11:32-105. OpenURL

  10. Hemmer H: Zur systematischen Stellung des Nebelparders, Neofelis nebulosa (Griffith, 1821).

    Saeugetierkundliche Mitteilungen Muenchen 1964, 12:13-17. OpenURL

  11. Hemmer H: Studien zur Ethologie des Nebelparders Neofelis nebulosa (Griffith 1821) und des Irbis Uncia uncia (Schreber 1775).

    Veroeffentlichungen der Zoologischen Staatssammlung Muenchen 1968, 12:155-247. OpenURL

  12. Yu L, Zhang YP: Phylogenetic studies of pantherine cats (Felidae) based on multiple genes, with novel application of nuclear beta-fibrinogen intron 7 to carnivores.

    Mol Phy Evol 2005, 35:483-495. Publisher Full Text OpenURL

  13. Johnson WE, O'Brien SJ: Phylogenetic reconstruction of the Felidae using 16S rRNA and NADH-5 mitochondrial genes.

    J Mol Evol 1997, 44(Suppl 1):98-116. Publisher Full Text OpenURL

  14. Johnson WE, Eizirik E, Pecon-Slattery J, Murphy WJ, Antunes A, Teeling E, O'Brien SJ: The late Miocene radiation of modern Felidae: a genetic assessment.

    Science 2006, 311:73-77. PubMed Abstract | Publisher Full Text OpenURL

  15. Buckley-Beason VA, Johnson WE, Nash WG, Stanyon R, Menninger JC, Driscoll CA, Howard JG, Bush M, Page JE, Roelke ME, Stone G, Martelli PP, Wen C, Ling L, Duraisingam RK, Lam PV, O'Brien SJ: Molecular evidence for species-level distinctions in clouded leopards.

    Curr Biol 2006, 16:2371-2376. PubMed Abstract | Publisher Full Text OpenURL

  16. Kitchener AC, Beaumont MA, Richardson D: Geographical variation in the clouded leopard, Neofelis nebulosa, reveals two species.

    Curr Biol 2006, 16:2377-2383. PubMed Abstract | Publisher Full Text OpenURL

  17. Mace GM: The role of taxonomy in species conservation.

    Phil Trans Roy Soc B 2004, 359:711-719. Publisher Full Text OpenURL

  18. Bensasson D, Zhang DX, Hartl DL, Hewitt GM: Mitochondrial pseudogenes: evolution's misplaced witnesses.

    Trends Ecol Evol 2001, 16:314-321. PubMed Abstract | Publisher Full Text OpenURL

  19. Lopez JV, Cevario S, O'Brien SJ: Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (Numt) in the nuclear genome.

    Genomics 1996, 33:229-246. PubMed Abstract | Publisher Full Text OpenURL

  20. Lopez JV, Yuhki N, Masuda R, Modi W, O'Brien SJ: Numt, a recent transfer and tandem amplification of mitochondrial DNA to the nuclear genome of the domestic cat.

    J Mol Evol 1994, 39:174-190. PubMed Abstract OpenURL

  21. Kim JH, Antunes A, Luo SJ, Menninger J, Nash WG, O'Brien SJ, Johnson WA: Evolutionary analysis of a large mtDNA translocation (numt) into the nuclear genome of the Panthera genus species.

    Gene 2006, 366:292-302. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Collura RV, Stewart C: Insertions and duplications of mtDNA in the nuclear genomes of Old World monkeys and hominoids.

    Nature 1995, 378:485-489. PubMed Abstract | Publisher Full Text OpenURL

  23. Menotti-Raymond M, David VA, Lyons LA, Schäffer AA, Tomlin JF, Hutton MK, O'Brien SJ: A genetic linkage map of microsatellites in the domestic cat (Felis catus).

    Genomics 1999, 57:9-23. PubMed Abstract | Publisher Full Text OpenURL

  24. Williamson JE, Huebinger RM, Sommer JA, Louis EE, Barber RC: Development and cross-species amplification of 18 microsatellite markers in the Sumatran tiger (Panthera tigris sumatrae).

    Mol Ecol Notes 2002, 2:110-112. Publisher Full Text OpenURL

  25. Corander J, Waldmann P, Sillanpaa MJ: Bayesian analysis of genetic differentiation between populations.

    Genetics 2003, 163:367-374. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Corander J, Marttinen P: Bayesian identification of admixture events using multilocus molecular markers.

    Mol Ecol 2006, 15:2833-2843. PubMed Abstract | Publisher Full Text OpenURL

  27. Corander J, Marttinen P, Sirén J, Tang J: BAPS: BayesianAnalysis of Population Structure, Manual v. 4.14. Department of Mathematics and statistics, University of Helsinki; 2006. OpenURL

  28. Avise JC, Ball RM: Principles of genealogical concordance in species concepts and biological taxonomy.

    Oxf Surv Evol Biol 1990, 7:45-67. OpenURL

  29. O'Brien SJ, Mayr E: Bureaucratic mischief: recognizing endangered species and subspecies.

    Science 1991, 251:1187-1188. Publisher Full Text OpenURL

  30. Avise JC: Molecular markers, natural history and evolution. Chapman & Hall, New York; 1994.

  31. Crandall KA, Bininda-Emonds ORP, Mace GM, Wayne RK: Considering evolutionary processes in conservation biology.

    Trends Ecol Evol 2000, 17:390-395. OpenURL

  32. Fraser DJ, Bernatchez L: Adaptive evolutionary conservation: towards a unified concept for defining conservation units.

    Mol Ecol 2001, 10:2741-2752. PubMed Abstract | Publisher Full Text OpenURL

  33. Ryder OA: Species conservation and systematics: the dilemma of subspecies.

    Trends Ecol Evol 1986, 1:9-10. Publisher Full Text OpenURL

  34. Waples RS: Pacific salmon, Oncorhynchus spp. & the definition of 'species' under the endangered species act.

    Mar Fish Reviews 1991, 53:11-22. OpenURL

  35. Waples RS: Evolutionary significant units and the conservation of biological diversity under the endangered species act. In Evolution and the aquatic ecosystem: Defining unique units in population conservation. Edited by Nielson JL, Powers GA. American Fisheries Society, Bethesda, Maryland; 1995:8-27. OpenURL

  36. Santiapillai C, Ashby KR: The Clouded leopard in Sumatra.

    Oryx 1988, 22:44-45. OpenURL

  37. IUCN The World Conservation Union: Neofelis nebulosa – status accessed 2002.

    2006 IUCN Red List of Threatened Species 2006. OpenURL

  38. Tjia HD: The Sunda shelf, Southeast Asia.

    Zeitschrift fuer Geomorphologie 1980, 24:405-427. OpenURL

  39. Heaney LR: A synopsis of climatic and vegetational change in southeast Asia.

    Climatic Change 1991, 19:53-61. Publisher Full Text OpenURL

  40. Voris HK: Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations.

    J Biogeography 2000, 27:1153-1167. Publisher Full Text OpenURL

  41. Heaney LR: Zoogeographic evidence for Middle and Late Pleistocene land bridges to the Philippine islands.

    Modern Quaternary Research SE Asia 1985, 9:127-143. OpenURL

  42. Koopman KF: Distributional patterns of the Indo-Malayan bats (Mammalia: Chiroptera).

    Am Mus Nov 1989, 2942:1-19. OpenURL

  43. Ruedi M: Phylogenetic evolution and biogeography of Southeast Asian shrews (genus Crocidura : Soricidae).

    Biol J Linn Soc 1996, 58:197-219. Publisher Full Text OpenURL

  44. Karns DR, O'Bannon A, Voris HK, Weigt LA: Biogeographical implications of mitochondrial DNA variation in the Bockadam snake (Cerberus rynchops, Serpentes: Homalopsinae) in Southeast Asia.

    J Biogeography 2000, 27:391-401. Publisher Full Text OpenURL

  45. Hewitt G: The genetic legacy of the Quaternary ice ages.

    Nature 2000, 405:907-913. PubMed Abstract | Publisher Full Text OpenURL

  46. Inger RF, Voris HK: The biogeographical relations of the frogs and snakes of Sundaland.

    J Biogeography 2001, 28:863-891. Publisher Full Text OpenURL

  47. Gorog AJ, Sinaga MH, Engstrom MD: Vicariance or dispersal? Historical biogeography of three Sunda shelf murine rodents (Maxomys surifer, Leopoldamys sabanus and Maxomys whiteheadi).

    Biol J Linn Soc 2004, 81:91-109. Publisher Full Text OpenURL

  48. Brandon-Jones D: The Asian Colobinae (Mammalia: Cercopithecidae) as indicators of Quaternary climatic change.

    Biol J Linn Soc 1996, 59:327-350. Publisher Full Text OpenURL

  49. Steiper ME: Population history, biogeography, and taxonomy of orangutans (Genus: Pongo) based on a population genetic meta-analysis of multiple loci.

    J Hum Evol 2006, 50:509-522. PubMed Abstract | Publisher Full Text OpenURL

  50. Warren KS, Verschoor EJ, Langenhuijzen S, Heryanto Swan RA, Vigilant L, Heeney JL: Speciation and intrasubspecific variation of Bornean orangutans Pongo pygmaeus pygmaeus.

    Mol Biol Evol 2001, 18:471-480. OpenURL

  51. Kanthaswamy S, Kurushima JD, Smith DG: Inferring Pongo conservation units: a perspective based on microsatellite and mitochondrial DNA analyses.

    Primates 2006, 47:310-321. PubMed Abstract | Publisher Full Text OpenURL

  52. Sambrook J, Fitsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, Cold Spring Harbor Press; 1989. OpenURL

  53. Jae-Heup K, Eizirik E, O'Brien SJ, Johnson WE: Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats.

    Mitochondrion 2001, 1:279-292. PubMed Abstract | Publisher Full Text OpenURL

  54. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.

    Nucleic Acids Res 1994, 22:4673-4680. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.

    Nucleic Acids Symp Ser 1999, 41:95-98. OpenURL

  56. Swofford DL: PAUP* phylogenetic analysis using parsimony and other methods. Sunderland, MA: Sinauer Associates; 2001. OpenURL

  57. Huelsenbeck JP, Bull JJ, Cunningham CW: Combining data in phylogenetic analysis.

    Trends Ecol Evol 1996, 11:152-157. Publisher Full Text OpenURL

  58. Eyre-Walker A, Awadalla P: Does human mtDNA recombine?

    J Mol Evol 2001, 51:430-435. Publisher Full Text OpenURL

  59. Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees.

    Mol Biol Evol 1987, 4:406-425. PubMed Abstract | Publisher Full Text OpenURL

  60. Hasegawa M, Kishino HYT: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.

    J Mol Evol 1985, 22:160-174. PubMed Abstract | Publisher Full Text OpenURL

  61. Takezaki N, Rzhetsky A, Nei M: Phylogenetic test of the molecular clock and linearized trees.

    Mol Biol Evol 1995, 12:823-833. PubMed Abstract | Publisher Full Text OpenURL

  62. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences.

    J Mol Evol 1980, 16:112-120. Publisher Full Text OpenURL

  63. Raymond M, Rousset F: GENEPOP (Version 1.2): A population genetic software for exact test and ecumenicism. [http://wbiomed.curtin.edu.au/genepop/] webcite

    J Heredity 1995, 86:248-249.

    Version 3.4

    OpenURL

  64. Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis.

    Evol Bioinformatics online 2005, 1:47-50. OpenURL

  65. Minch E, Ruiz-Linares A, Goldstein DB: MICROSAT. [http://hpgl.stanford.edu/projects/microsat/] webcite

    1995.

    Accessed 13 November 2006

  66. Felsenstein J: PHYLIP (Phylogeny inference package) version 3.66. Distributed by the author. Department of Genome Sciences, University of Washington, Seattle.; 2006. OpenURL

  67. Page RDM: Tree View: An application to display phylogenetic trees on personal computers.

    Comput Appl Biosc 1996, 12:357-358. OpenURL