Morphometrical relationships between South-east Asian deer (Cervidae, tribe Cervini): evolutionary and biogeographic implications

Antler tine homologies and cervid systematics: A review of past and present controversies with special emphasis on Elaphurus davidianus

Antlers are the most conspicuous trait of cervids and have been used in the past to establish a classification of their fossil and living representatives. Since the availability of molecular data, morphological characters have generally become less important for phylogenetic reconstructions. In recent years, however, the appreciation of morphological characters has increased, and they are now more frequently used in addition to molecular data to reconstruct the evolutionary history of cervids. A persistent challenge when using antler traits in deer systematics is finding a consensus on the homology of structures. Here, we review early and recent attempts to homologize antler structures and objections to this approach, compare and evaluate recent advances on antler homologies, and critically discuss these different views in order to offer a basis for further scientific exchange on the topic. We further present some developmental aspects of antler branching patterns and discuss their potential for reconstructing cervid systematics. The use of heterogeneous data for reconstructing phylogenies has resulted in partly conflicting hypotheses on the systematic position of certain cervid species, on which we also elaborate here. We address current discussions on the use of different molecular markers in cervid systematics and the question whether antler morphology and molecular data can provide a consistent picture on the evolutionary history of cervids. In this context, special attention is given to the antler morphology and the systematic position of the enigmatic Pere David's deer (Elaphurus davidianus).

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Genetic analysis reveals a distinct lineage of hog deer (Axis porcinus) in Kratie province, Cambodia

The hog deer (Axis porcinus) is an endangered cervid with drastic population declines. There are two recognized subspecies of hog deer: A. p. porcinus, ranging from Punjab Province in Pakistan, Nepal and the Northern part of India to Myanmar, and A. p. annamiticus found in Indo-China, Thailand, Laos, Cambodia and Vietnam. The current geographic range of A. p. annamiticus is still ambiguous. We analyzed variation in the mitochondrial DNA control region (mtDNA CR) to investigate the intra-species structure, differentiation, and demographic history of hog deer from Cambodia (Kratie Province), which we compared with the populations from India and Thailand. We also generated divergence time estimates using a concatenated dataset of complete Cyt b and partial CR. The CR data showed that Cambodian hog deer are genetically differentiated from the mainland Indian and Thai populations, forming a distinct basal clade. The time of divergence indicates that the Cambodian lineage split from the other two hog deer lineages around 0.51 Mya, during the Late Pleistocene. The results also suggest strong phylogeographic structure among hog deer: lineage A extends from Terai Arc (foothills of the Himalayas) to Assam, India (A. p. porcinus), lineage B from Manipur, India to Thailand (A. p. annamiticus), and lineage C is only known from Kratie Province, Cambodia. Lineage A exhibited a higher level of genetic diversity than lineages B and C, with recent demographic stability. Thus, the hog deer population in Kratie Province appears to be a distinct lineage that should be treated as an evolutionarily significant unit.

The island rule explains consistent patterns of body size evolution in terrestrial vertebrates

Island faunas can be characterized by gigantism in small animals and dwarfism in large animals, but the extent to which this so-called 'island rule' provides a general explanation for evolutionary trajectories on islands remains contentious. Here we use a phylogenetic meta-analysis to assess patterns and drivers of body size evolution across a global sample of paired island-mainland populations of terrestrial vertebrates. We show that 'island rule' effects are widespread in mammals, birds and reptiles, but less evident in amphibians, which mostly tend towards gigantism. We also found that the magnitude of insular dwarfism and gigantism is mediated by climate as well as island size and isolation, with more pronounced effects in smaller, more remote islands for mammals and reptiles. We conclude that the island rule is pervasive across vertebrates, but that the implications for body size evolution are nuanced and depend on an array of context-dependent ecological pressures and environmental conditions.

Unexpectedly high levels of lineage diversity in Sundaland puddle frogs (Dicroglossidae: Occidozyga Kuhl and van Hasselt, 1822)

One of the most urgent contemporary tasks for taxonomists and evolutionary biologists is to estimate the number of species on earth. Recording alpha diversity is crucial for protecting biodiversity, especially in areas of elevated species richness, which coincide geographically with increased anthropogenic environmental pressures - the world's so-called biodiversity hotspots. Although the distribution of Puddle frogs of the genus Occidozyga in South and Southeast Asia includes five biodiversity hotspots, the available data on phylogeny, species diversity, and biogeography are surprisingly patchy. Samples analyzed in this study were collected throughout Southeast Asia, with a primary focus on Sundaland and the Philippines. A mitochondrial gene region comprising ~ 2000 bp of 12S and 16S rRNA with intervening tRNA Valine and three nuclear loci (BDNF, NTF3, POMC) were analyzed to obtain a robust, time-calibrated phylogenetic hypothesis. We found a surprisingly high level of genetic diversity within Occidozyga, based on uncorrected p-distance values corroborated by species delimitation analyses. This extensive genetic diversity revealed 29 evolutionary lineages, defined by the > 5% uncorrected p-distance criterion for the 16S rRNA gene, suggesting that species diversity in this clade of phenotypically homogeneous forms probably has been underestimated. The comparison with results of other anuran groups leads to the assumption that anuran species diversity could still be substantially underestimated in Southeast Asia in general. Many genetically divergent lineages of frogs are phenotypically similar, indicating a tendency towards extensive morphological conservatism. We present a biogeographic reconstruction of the colonization of Sundaland and nearby islands which, together with our temporal framework, suggests that lineage diversification centered on the landmasses of the northern Sunda Shelf. This remarkably genetically structured group of amphibians could represent an exceptional case for future studies of geographical structure and diversification in a widespread anuran clade spanning some of the most pronounced geographical barriers on the planet (e.g., Wallace's Line). Studies considering gene flow, morphology, ecological and bioacoustic data are needed to answer these questions and to test whether observed diversity of Puddle frog lineages warrants taxonomic recognition.

Functional traits of the world's late Quaternary large-bodied avian and mammalian herbivores

Prehistoric and recent extinctions of large-bodied terrestrial herbivores had significant and lasting impacts on Earth's ecosystems due to the loss of their distinct trait combinations. The world's surviving large-bodied avian and mammalian herbivores remain among the most threatened taxa. As such, a greater understanding of the ecological impacts of large herbivore losses is increasingly important. However, comprehensive and ecologically-relevant trait datasets for extinct and extant herbivores are lacking. Here, we present HerbiTraits, a comprehensive functional trait dataset for all late Quaternary terrestrial avian and mammalian herbivores ≥10 kg (545 species). HerbiTraits includes key traits that influence how herbivores interact with ecosystems, namely body mass, diet, fermentation type, habitat use, and limb morphology. Trait data were compiled from 557 sources and comprise the best available knowledge on late Quaternary large-bodied herbivores. HerbiTraits provides a tool for the analysis of herbivore functional diversity both past and present and its effects on Earth's ecosystems.

A review of the genus Rusa in the indo-malayan archipelago and conservation efforts

Genus Rusa, belonging to the deer family Cervidae is native to the Indo-Malaya Archipelago (IMA). However, detailed information on the Rusa genus in the IMA is limited. This review provides comprehensive information on the Rusa genus in the IMA including, threats and conservation efforts. There are four species of deer in Rusa genus, which is Sambar deer (Rusa unicolor), Javan deer (Rusa timorensis), Visayan spotted deer (Rusa alfredi) and Philippine deer (Rusa marianna). Despite their wide distribution in the South Asian and Southeast Asian regions, they are under serious threats. Some conservation efforts that are being done to protect and conserve them among others are; (1) facilities protection, (2) habitat enrichment programme, (3) Ex-situ conservation, (4) legislations, and (5) captive breeding. Conservation through genetics is also an important step in conserving these species. Recommendations for conservation of the genus are also discussed; 1. maintenance of ecosystem. 2. more effective monitoring system on the existing protected area. 3. ex-situ conservation, and 4. habitat monitoring.

The systematics of the Cervidae: a total evidence approach

Systematic relationships of cervids have been controversial for decades. Despite new input from molecular systematics, consensus could only be partially reached. The initial, gross (sub) classification based on morphology and comparative anatomy was mostly supported by molecular data. The rich fossil record of cervids has never been extensively tested in phylogenetic frameworks concerning potential systematic relationships of fossil cervids to extant cervids. The aim of this work was to investigate the systematic relationships of extant and fossil cervids using molecular and morphological characters and make implications about their evolutionary history based on the phylogenetic reconstructions. To achieve these objectives, molecular data were compiled consisting of five nuclear markers and the complete mitochondrial genome of 50 extant and one fossil cervids. Several analyses using different data partitions, taxon sampling, partitioning schemes, and optimality criteria were undertaken. In addition, the most extensive morphological character matrix for such a broad cervid taxon sampling was compiled including 168 cranial and dental characters of 41 extant and 29 fossil cervids. The morphological and molecular data were analysed in a combined approach and other comprehensive phylogenetic reconstructions. The results showed that most Miocene cervids were more closely related to each other than to any other cervids. They were often positioned between the outgroup and all other cervids or as the sister taxon to Muntiacini. Two Miocene cervids were frequently placed within Muntiacini. Plio- and Pleistocene cervids could often be affiliated to Cervini, Odocoileini or Capreolini. The phylogenetic analyses provide new insights into the evolutionary history of cervids. Several fossil cervids could be successfully related to living representatives, confirming previously assumed affiliations based on comparative morphology and introducing new hypotheses. New systematic relationships were observed, some uncertainties persisted and resolving systematics within certain taxa remained challenging.

Demographic and genetic structure of a severely fragmented population of the endangered hog deer (Axis porcinus) in the Indo-Burma biodiversity hotspot

The population of the globally endangered hog deer (Axis porcinus) has declined severely across its geographic range. Intensive monitoring of its demographic and genetic status is necessary. We examined the demographic and genetic structure of a small hog deer population in Keibul Lamjao National Park (KLNP), located on the western fringe of the Indo-Burma biodiversity hotspot for conservation planning. The distribution pattern of hog deer in the Park was derived based on the presence/absence of faecal pellets in 1 km × 1 km grids. We used double-observer distance sampling method to derive the hog deer abundance and population structure and compared with previous data to derive the population trend. We determined the genetic diversity of the population through microsatellite screening and bottleneck detection. The overall pellet density was 0.34 ± 0.02 pellets km^-2 restricted to only 22.34 ± 0.20 km² area of the park. The estimated density of the deer in the park was 1.82–4.32 individuals km^-2. The population showed a declining trend from 2006–08 (p < 0.05, R² = 0.916) with 8% annum^-1 and an increasing trend from 2003–2018 (p < 0.05, R² = 0.9304) with 10% annum^-1. The adult male-to-female ratio and fawn-to-doe ratio were 36.2 ± 1.9 males per 100 females and 16.5 ± 0.4 fawns per 100 females, respectively. The molecular examination suggested that the mean number of alleles at 23 loci was 2.70 ± 0.18, the observed heterozygosity (Ho) ranged from 0.26 to 0.63 (mean 0.42 ± 0.02), the expected heterozygosity (He) ranged from 0.23 to 0.73 (χ = 0.51 ± 0.03), and the polymorphic information content (PIC) ranged from 0.2 to 0.67 (χ = 0.43 ± 0.03) indicating a moderate level of genetic diversity. Although no bottleneck in the population was observed, the loss of genetic diversity may affect the evolutionary potential of the species at the site by limiting the selection flexibility. Conservation planning coupled with scientific management regime will help in the long term persistence of the population in the region.

Variability, morphometrics, and co-variation of the os lacrimale in Cervidae

In Ruminantia, the lacrimal bone forms a considerable part of the facial skeleton, and the morphology of its facial facet is highly variable when compared to other mammals. In this study, we quantify the species-specific variability in size and shape of the lacrimal facial facet in species of Cervidae (deer) and relate it to systematics and various aspects of their ecology and behavior. We sampled 143 skull specimens from 10 genera; 12 Moschus and 3 Tragulus specimens were used as outgroups. We find that size and shape of the lacrimal facial facet allow differentiating most species analyzed here, except for Mazama gouazoubira and Capreolus capreolus. Size and shape of the lacrimal facial facet vary widely across Cervidae regardless of their systematic relationships, ecology or behavior. Thus, we could not detect a unique signature of adaptational criteria in lacrimal morphology. Our data indicate that the lacrimal facial facet scales allometrically with skull size, in particular, the lacrimojugal length scales positively and the lacrimomaxillar length scales negatively. However, correlation analyses did not reveal any differences in the integration of the lacrimal bone with any specific skull module in any of the species compared. Lastly, we could not ascertain any correlation between the size and position of the preorbital depression with the size and shape of the lacrimal facial facet. We conclude that the lacrimal facial facet is highly flexible and may rapidly adjust to its surrounding bones. Its allometric growth appears to be an example of exaptation: changes in size and shape in the context of the increase of the skull length provide lacrimal contacts, in particular, a lacrimojugal one, which may serve to reduce mechanical loads resulting from increasingly larger antlers in large cervids.

Genetic analysis of endangered hog deer (Axis porcinus) reveals two distinct lineages from the Indian subcontinent

The hog deer (Axis porcinus) is threatened by habitat alteration, fragmentation, and poaching, which have led to a drastic decline of its wild population. Two subspecies of A. porcinus have been described from its distribution range. A. p. porcinus is reported to occur from Pakistan along the Himalayan foothills through Nepal, India and Myanmar, and A. p. annamiticus is found in Thailand, Indo-China, Laos, Cambodia, and Vietnam. However, the current distribution range of A. p. annamiticus is still unclear. We used the partial control region (CR) of mitochondrial DNA (mtDNA) and seven microsatellite loci to investigate the intra-species structure, differentiation, and demographic history of hog deer populations from three landscapes, the Terai Arc, Northeast, and Indo-Burma (Keibul Lamjao National Park (KLNP), Manipur, India) landscapes. We also carried out divergence time estimation using the complete mitogenome. The level of variation was ~4%, and the time of divergence of the KLNP population and the other Indian populations was about 0.22 Mya, i.e., during the last glaciation periods of the Late Pleistocene/Early Holocene. The KLNP haplotypes of the control region were shared with the Southeast Asian subspecies, A. p. annamiticus. The results of the investigations of the microsatellite loci supported the mtDNA results unambiguously. Two genetically distinct lineages are found in India: one is found from the Terai Arc to Assam (A. p. porcinus) and the other in Manipur (A. p. annamiticus). The genetic diversity in KLNP was low and exhibited a higher degree of genetic differentiation compared with major Indian populations. The Bayesian skyline plots indicated that after a long phase of historic demographic stability, the populations of both the lineages of hog deer suffered pronounced declines during the period from ~800 years BP to 5000 years BP. In summary, our finding provided evidence that the KLNP population is probably a prime, isolated and sustaining stock of A. p. annamiticus and should be managed as evolutionarily significant units (ESUs).

The modular organization of roe deer (Capreolus capreolus) body during ontogeny: the effects of sex and habitat

Background

As a small artiodactyl, the roe deer (Capreolus capreolus L.) is characterized by biological plasticity and great adaptability demonstrated by their survival under a wide variety of environmental conditions. In order to depict patterns of phenotypic variation of roe deer body this study aims to quantify variation during ontogenetic development and determine how sex-specific reproductive investment and non-uniform habitat differences relate to phenotypic variation and do these differential investments mold the patterns of phenotypic variation through modular organisation.

Results

Patterns of phenotypic correlation among body traits change during the ontogeny of roe deer, with differential influence of sex and habitat type. Modularity was found to be a feature of closed habitats with trunk+forelimbs+hindlimbs as the best supported integration/modularity hypothesis for both sexes. The indices of integration and evolvability vary with habitat type, age and sex where increased integration is followed by decreased evolvability.

Conclusion

This is the first study that quantifies patterns of correlation in the roe deer body and finds pronounced changes in correlation structure during ontogeny affected by sex and habitat type. The correlation structure of the roe deer body is developmentally written over the course of ontogeny but we do not exclude the influence of function on ontogenetic changes. Modularity arises with the onset of reproduction (subadults not being modular) and is differentially expressed in males and females from different habitats. Both adult males and females show modularity in primordial, closed habitats. Overall, all these findings are important as they provide support to the idea that modularity can evolve at the population level and change fast within a species.

Human-mediated introduction of introgressed deer across Wallace's line: Historical biogeography of Rusa unicolor and R. timorensis

In this study we compared the phylogeographic patterns of two Rusa species, Rusa unicolor and Rusa timorensis, in order to understand what drove and maintained differentiation between these two geographically and genetically close species and investigated the route of introduction of individuals to the islands outside of the Sunda Shelf. We analyzed full mitogenomes from 56 archival samples from the distribution areas of the two species and 18 microsatellite loci in a subset of 16 individuals to generate the phylogeographic patterns of both species. Bayesian inference with fossil calibration was used to estimate the age of each species and major divergence events. Our results indicated that the split between the two species took place during the Pleistocene, ~1.8 Mya, possibly driven by adaptations of R. timorensis to the drier climate found on Java compared to the other islands of Sundaland. Although both markers identified two well-differentiated clades, there was a largely discrepant pattern between mitochondrial and nuclear markers. While nDNA separated the individuals into the two species, largely in agreement with their museum label, mtDNA revealed that all R. timorensis sampled to the east of the Sunda shelf carried haplotypes from R. unicolor and one Rusa unicolor from South Sumatra carried a R. timorensis haplotype. Our results show that hybridization occurred between these two sister species in Sundaland during the Late Pleistocene and resulted in human-mediated introduction of hybrid descendants in all islands outside Sundaland.

The novel primers for mammal species identification-based mitochondrial cytochrome b sequence: implication for reserved wild animals in Thailand and endangered mammal species in Southeast Asia

We presented the powerful techniques for species identification using the short amplicon of mitochondrial cytochrome b gene sequence. Two faecal samples and one single hair sample of the Asian tapir were tested using the new cytochrome b primers. The results showed a high sequence similarity with the mainland Asian tapir group. The comparative sequence analysis of the reserved wild mammals in Thailand and the other endangered mammal species from Southeast Asia comprehensibly verified the potential of our novel primers. The forward and reverse primers were 94.2 and 93.2%, respectively, by the average value of the sequence identity among 77 species sequences, and the overall mean distance was 35.9%. This development technique could provide rapid, simple, and reliable tools for species confirmation. Especially, it could recognize the problematic biological specimens contained less DNA material from illegal products and assist with wildlife crime investigation of threatened species and related forensic casework.

The Middle Pleistocene vertebrate fauna from Khok Sung (Nakhon Ratchasima, Thailand): biochronological and paleobiogeographical implications

The fluviatile terrace deposits of Khok Sung, Nakhon Ratchasima province, have yielded more than one thousand fossils, making this the richest Pleistocene vertebrate fauna of Thailand. The excellent preservation of the specimens allows precise characterization of the faunal composition. The mammalian fauna consists of fifteen species in thirteen genera, including a primate, a canid, a hyaenid, proboscideans, rhinoceroses, a suid, cervids, and bovids. Most species correspond to living taxa but globally (Stegodon cf. orientalis) and locally (Crocuta crocuta ultima, Rhinoceros unicornis, Sus barbatus, and Axis axis) extinct taxa were also present. The identification of Axis axis in Khok Sung, a chital currently restricted to the Indian Subcontinent, represents the first record of the species in Southeast Asia. Three reptilian taxa: Crocodylus cf. siamensis, Python sp., and Varanus sp., are also identified. Faunal correlations with other Southeast Asian sites suggest a late Middle to early Late Pleistocene age for the Khok Sung assemblage. However, the Khok Sung mammalian fauna is most similar to that of Thum Wiman Nakin, dated to older than 169 ka. The Khok Sung large mammal assemblage mostly comprises mainland Southeast Asian taxa that migrated to Java during the latest Middle Pleistocene, supporting the hypothesis that Thailand was a biogeographic pathway for the Sino-Malayan migration event from South China to Java.

Systematic relationships of five newly sequenced cervid species

Cervid phylogenetics has been puzzling researchers for over 150 years. In recent decades, molecular systematics has provided new input for both the support and revision of the previous results from comparative anatomy but has led to only partial consensus. Despite all of the efforts to reach taxon-wide species sampling over the last two decades, a number of cervid species still lack molecular data because they are difficult to access in the wild. By extracting ancient DNA from museum specimens, in this study, we obtained partial mitochondrial cytochrome b gene sequences for Mazama bricenii, Mazama chunyi, Muntiacus atherodes, Pudu mephistophiles, and Rusa marianna, including three holotypes. These new sequences were used to enrich the existing mitochondrial DNA alignments and yielded the most taxonomically complete data set for cervids to date. Phylogenetic analyses provide new insights into the evolutionary history of these five species. However, systematic uncertainties within Muntiacus persist and resolving phylogenetic relationships within Pudu and Mazama remain challenging.

First records of preorbital gland opening in rare wild barasingha (Rucervus duvaucelii) in social contexts may help to explain this phenomenon in cervids

The opening of the preorbital gland in deer serves as a visual communication and has been linked to a wide variety of behavioural situations. As we reported recently, all previous long-term studies on preorbital gland opening were carried out on only one species, whereas case reports on six other rarely studied species have shown associations with different behaviours, thus greatly increasing our overall understanding of the real function of this gland in animal visual communication. Here, we report for the first time preorbital gland opening in the barasingha (Rucervus duvaucelii) in social contexts as observed in a wild population in Kanha National Park, India. We observed this behaviour in two different contexts: agonistic and sexual. Moreover, our record of preorbital gland opening during copulation is the first one amongst cervids. Our findings of preorbital gland opening in both contexts in wild barasingha indicate that the gland was opened only when the individual was highly excited. We suggest that preorbital gland opening may be an important behavioural indicator of an individual involved in a serious intraspecific interaction, and thus a useful tool with which to distinguish between playful and serious behaviours, especially in agonistic and sexual situations.

Forest fruit production is higher on Sumatra than on Borneo

Background

Various studies have shown that the population densities of a number of forest vertebrates, such as orangutans, are higher on Sumatra than Borneo, and that several species exhibit smaller body sizes on Borneo than Sumatra and mainland Southeast Asia. It has been suggested that differences in forest fruit productivity between the islands can explain these patterns. Here we present a large-scale comparison of forest fruit production between the islands to test this hypothesis.

Methodology/Principal Findings

Data on fruit production were collated from Sumatran and Bornean sites. At six sites we assessed fruit production in three forest types: riverine, peat swamp and dryland forests. We compared fruit production using time-series models during different periods of overall fruit production and in different tree size classes. We examined overall island differences and differences specifically for fruiting period and tree size class. The results of these analyses indicate that overall the Sumatran forests are more productive than those on Borneo. This difference remains when each of the three forest types (dryland, riverine, and peat) are examined separately. The difference also holds over most tree sizes and fruiting periods.

Conclusions/Significance

Our results provide strong support for the hypothesis that forest fruit productivity is higher on Sumatra than Borneo. This difference is most likely the result of the overall younger and more volcanic soils on Sumatra than Borneo. These results contribute to our understanding of the determinants of faunal density and the evolution of body size on both islands.

Mitochondrial and nuclear phylogenies of Cervidae (Mammalia, Ruminantia): Systematics, morphology, and biogeography

The family Cervidae includes 40 species of deer distributed throughout the northern hemisphere, as well as in South America and Southeast Asia. Here, we examine the phylogeny of this family by analyzing two mitochondrial protein-coding genes and two nuclear introns for 25 species of deer representing most of the taxonomic diversity of the family. Our results provide strong support for intergeneric relationships. To reconcile taxonomy and phylogeny, we propose a new classification where the family Cervidae is divided in two subfamilies and five tribes. The subfamily Cervinae is composed of two tribes: the tribe Cervini groups the genera Cervus, Axis, Dama, and Rucervus, with the Père David's deer (Elaphurus davidianus) included in the genus Cervus, and the swamp deer (Cervus duvauceli) placed in the genus Rucervus; the tribe Muntiacini contains Muntiacus and Elaphodus. The subfamily Capreolinae consists of the tribes Capreolini (Capreolus and Hydropotes), Alceini (Alces), and Odocoileini (Rangifer + American genera). Deer endemic to the New World fall in two biogeographic lineages: the first one groups Odocoileus and Mazama americana and is distributed in North, Central, and South America, whereas the second one is composed of South American species only and includes Mazama gouazoubira. This implies that the genus Mazama is not a valid taxon. Molecular dating suggests that the family originated and radiated in central Asia during the Late Miocene, and that Odocoileini dispersed to North America during the Miocene/Pliocene boundary, and underwent an adaptive radiation in South America after their Pliocene dispersal across the Isthmus of Panama. Our phylogenetic inferences show that the evolution of secondary sexual characters (antlers, tusk-like upper canines, and body size) has been strongly influenced by changes in habitat and behaviour.

Evolution and phylogeny of old world deer

The phylogenetic pattern and timing of the radiation of Old World deer was determined based on the complete mitochondrial cytochrome b gene from 33 Cervinae taxa. Using rooted and unrooted phylogenies derived from distinct theoretical approaches, strong support was achieved for monophyly of the Old World deer with muntjacs as sister group as well as for the divergence of at least three distinct genera: Rucervus, Dama, and Cervus. The latter clade comprises what have previously been regarded as the genera or subgenera Panolia, Rusa, Cervus, Sika, and probably Przewalskium. Our data also consistently confirmed paraphyly of nominate C. elaphus and did not support the monophyly of Axis. We used these molecular phylogenies to assess the homoplastic evolution of morphological, geographical, ecological, and selected behavioural character state differences within the Cervinae. Reliable fossil calibrations, large molecular data sets, and improved dating methods are shaping a molecular time scale for the evolutionary radiation of Old World deer that occurred at the Miocene/Pliocene transition and is largely compatible with existing palaeontological evidence. Using node ages estimated from sequence data, we estimated an average per-lineage diversification rate of 0.51+/-0.1 species per million years (my) over roughly the last 6 mya.