Draft genome of the lowland anoa (Bubalus depressicornis) and comparison with buffalo genome assemblies (Bovidae, Bubalina).

Genomic data for wild species of the genus Bubalus (Asian buffaloes) are still lacking while several whole genomes are currently available for domestic water buffaloes. To address this, we sequenced the genome of a wild endangered dwarf buffalo, the lowland anoa (Bubalus depressicornis), produced a draft genome assembly and made comparison to published buffalo genomes. The lowland anoa genome assembly was 2.56 Gbp long and contained 103,135 contigs, the longest contig being 337.39 kbp long. N50 and L50 values were 38.73 and 19.83 kbp, respectively, mean coverage was 44× and GC content was 41.74%. Two strategies were adopted to evaluate genome completeness: (1) determination of genomic features with de novo and homology-based predictions using annotations of chromosome-level genome assembly of the river buffalo and (2) employment of benchmarking against universal single-copy orthologs (BUSCO). Homology-based predictions identified 94.51% complete and 3.65% partial genomic features. De novo gene predictions identified 32,393 genes, representing 97.14% of the reference's annotated genes, whilst BUSCO search against the mammalian orthologs database identified 71.1% complete, 11.7% fragmented, and 17.2% missing orthologs, indicating a good level of completeness for downstream analyses. Repeat analyses indicated that the lowland anoa genome contains 42.12% of repetitive regions. The genome assembly of the lowland anoa is expected to contribute to comparative genome analyses among bovid species.

Adaptive radiation and speciation in Rhipicephalus ticks: A medley of novel hosts, nested predator-prey food webs, off-host periods and dispersal along temperature variation gradients.

Rhipicephalus are a species-diverse genus of ticks, mainly distributed in the Afrotropics with some species in the Palearctic and Oriental regions. Current taxonomic consensus comprise nine informal species groups/lineages based on immature morphology. This work integrates biogeographic, ecological and molecular lines of evidence to better understand Rhipicephalus evolution. Phylogenetic analysis based on four genes (12S, 16S, 28S-D2 and COI) recovered five distinct clades with nine descendant clades that are generally congruent with current taxonomy, with some exceptions. Historical biogeography is inferred from molecular divergence times, ancestral distribution areas, host-use and climate niches of four phylogenetically significant bioclimatic variables (isothermality, annual, seasonal and diurnal temperature range). Novel hosts enabled host-linked dispersal events into new environments, and ticks exploited new hosts through nested predator-prey connections in food webs. Diversification was further induced by climate niche partitioning along gradients in temperature range during off-host periods. Ancestral climate niche estimates corroborated dispersal events by indicating hypothetical ancestors moved into environments with different annual and seasonal temperature ranges along latitudinal gradients. Host size for immature and adult life stages was important for dispersal and subsequent diversification rates. Clades that utilise large, mobile hosts (ungulates and carnivores) early in development have wider geographic ranges but slower diversification rates, and those utilising small, less mobile hosts (rodents, lagomorphs and afroinsectivores) early in development have smaller ranges but higher diversification rates. These findings suggest diversification is driven by a complex set of factors linked to both host-associations (host size, ranges and mobility) and climate niche partitioning along annual and seasonal temperature range gradients that vary with latitude. Moreover, competitive interactions can reinforce these processes and drive speciation. Off-host periods facilitate adaptive radiation by enabling host switches along nested predator-prey connections in food webs, but at the cost of environmental exposure that partitions niches among dispersing progenitors, disrupting geneflow and driving diversification. As such, the evolution and ecological niches of Rhipicephalus are characterised by trade-offs between on- and off-host periods, and these trade-offs interact with nested predator-prey connections in food webs, host-use at different life stages, as well as gradients in annual and seasonal temperature ranges to drive adaptive radiation and speciation.

Correction: Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

[This corrects the article DOI: 10.1371/journal.pone.0240770.].

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

Characterization of the complete mitochondrial genome of Diplostomum baeri.

Despite Diplostomum baeri (Dubois, 1937) being one of the most widely distributed parasites of freshwater fish, there is no complete mitochondrial (mt) genome currently available. The complicated systematics presented by D. baeri has hampered investigations into the species distributions and infective dynamics of the species. Within this study we obtained complete mt genome sequences of D. baeri and assessed its phylogenetic relationship with other species of Digenea. The complete mitochondrial genome of D. baeri is 14,480 bp in length, containing 36 genes in total. The phylogenetic tree resulting from Bayesian inference of concatenated 12 protein coding gene sequences placed D. baeri alongside published mt genomes of Diplostomidae, with the overall taxonomic placement of the genus being a sister lineage of the order Plagiochiida The characterization of further mitochondrial genomes within the family Diplostomidae will help progress phylogenetic and epidemiological investigations as well as providing a framework for the analysis of diagnostic markers to be used in further monitoring of the parasite worldwide.

Multi-locus phylogeny of the tribe Tragelaphini (Mammalia, Bovidae) and species delimitation in bushbuck: Evidence for chromosomal speciation mediated by interspecific hybridization.

The bushbuck is the most widespread bovid species in Africa. Previous mitochondrial studies have revealed a polyphyletic pattern suggesting the possible existence of two distinct species. To assess this issue, we have sequenced 16 nuclear genes and one mitochondrial fragment (cytochrome b gene + control region) for most species of the tribe Tragelaphini, including seven bushbuck individuals belonging to the two divergent mtDNA haplogroups, Scriptus and Sylvaticus. Our phylogenetic analyses show that the Scriptus lineage is a sister-group of Sylvaticus in the nuclear tree, whereas it is related to Tragelaphus angasii in the mitochondrial tree. This mito-nuclear discordance indicates that the mitochondrial genome of Scriptus was acquired by introgression after one or several past events of hybridization between bushbuck and an extinct species closely related to T. angasii. The division into two bushbuck species is supported by the analyses of nuclear markers and by the karyotype here described for T. scriptus (2n = 57 M/58F), which is strikingly distinct from the one previously found for T. sylvaticus (2n = 33 M/34F). Molecular dating estimates suggest that the two species separated during the Early Pleistocene after an event of interspecific hybridization, which may have mediated massive chromosomal rearrangements in the common ancestor of T. scriptus.

The importance of correctly identifying the process responsible for spatial genetic structure in Leopard: A response to McManus and Smuts (2016).

Microsatellite analyses suggest that spatial genetic structure among six leopard-sampling sites in southern Africa is the result of isolation by distance.

Implications of spatial genetic patterns for conserving African leopards.

The leopard (Panthera pardus) is heavily persecuted in areas where it predates livestock and threatens human well-being. Attempts to resolve human-leopard conflict typically involve translocating problem animals; however, these interventions are rarely informed by genetic studies and can unintentionally compromise the natural spatial genetic structure and diversity, and possibly the long-term persistence, of the species. No significant genetic discontinuities were definable within the southern African leopard population. Analysis of fine-scale genetic data derived from mitochondrial and nuclear DNA revealed that the primary natural process shaping the spatial genetic structure of the species is isolation-by-distance (IBD). The effective gene dispersal (σ) index can inform leopard translocations and is estimated to be 82 km for some South African leopards. The importance of adopting an evidence-based strategy is discussed for supporting the integration of genetic data, spatial planning and social learning institutions so as to promote collaboration between land managers, government agency staff and researchers.

Cryptic speciation and chromosomal repatterning in the South African climbing mice Dendromus (Rodentia, Nesomyidae).

We evaluate the intra- and interspecific diversity in the four South African rodent species of the genus Dendromus. The molecular phylogenetic analysis on twenty-three individuals have been conducted on a combined dataset of nuclear and mitochondrial markers. Moreover, the extent and processes underlying chromosomal variation, have been investigated on three species by mean of G-, C-bands, NORs and Zoo-FISH analysis. The molecular analysis shows the presence of six monophyletic lineages corresponding to D. mesomelas, D. mystacalis and four lineages within D. cfr. melanotis with high divergence values (ranges: 10.6% - 18.3%) that raises the question of the possible presence of cryptic species. The first description of the karyotype for D. mesomelas and D. mystacalis and C- and G- banding for one lineage of D. cfr. melanotis are reported highlighting an extended karyotype reorganization in the genus. Furthermore, the G-banding and Zoo-FISH evidenced an autosome-sex chromosome translocation characterizing all the species and our timing estimates this mutation date back 7.4 mya (Late Miocene). Finally, the molecular clock suggests that cladogenesis took place since the end of Miocene to Plio-Pleistocene, probably due to ecological factors, isolation in refugia followed by differential adaptation to the mesic or dry habitat.

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla, Bovidae).

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships. Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota. The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.

Pattern and timing of diversification of Cetartiodactyla (Mammalia, Laurasiatheria), as revealed by a comprehensive analysis of mitochondrial genomes.

The order Cetartiodactyla includes cetaceans (whales, dolphins and porpoises) that are found in all oceans and seas, as well as in some rivers, and artiodactyls (ruminants, pigs, peccaries, hippos, camels and llamas) that are present on all continents, except Antarctica and until recent invasions, Australia. There are currently 332 recognized cetartiodactyl species, which are classified into 132 genera and 22 families. Most phylogenetic studies have focused on deep relationships, and no comprehensive time-calibrated tree for the group has been published yet. In this study, 128 new complete mitochondrial genomes of Cetartiodactyla were sequenced and aligned with those extracted from nucleotide databases. Our alignment includes 14,902 unambiguously aligned nucleotide characters for 210 taxa, representing 183 species, 107 genera, and all cetartiodactyl families. Our mtDNA data produced a statistically robust tree, which is largely consistent with previous classifications. However, a few taxa were found to be para- or polyphyletic, including the family Balaenopteridae, as well as several genera and species. Accordingly, we propose several taxonomic changes in order to render the classification compatible with our molecular phylogeny. In some cases, the results can be interpreted as possible taxonomic misidentification or evidence for mtDNA introgression. The existence of three new cryptic species of Ruminantia should therefore be confirmed by further analyses using nuclear data. We estimate divergence times using Bayesian relaxed molecular clock models. The deepest nodes appeared very sensitive to prior assumptions leading to unreliable estimates, primarily because of the misleading effects of rate heterogeneity, saturation and divergent outgroups. In addition, we detected that Whippomorpha contains slow-evolving taxa, such as large whales and hippos, as well as fast-evolving taxa, such as river dolphins. Our results nevertheless indicate that the evolutionary history of cetartiodactyls was punctuated by four main phases of rapid radiation during the Cenozoic era: the sudden occurrence of the three extant lineages within Cetartiodactyla (Cetruminantia, Suina and Tylopoda); the basal diversification of Cetacea during the Early Oligocene; and two radiations that involve Cetacea and Pecora, one at the Oligocene/Miocene boundary and the other in the Middle Miocene. In addition, we show that the high species diversity now observed in the families Bovidae and Cervidae accumulated mainly during the Late Miocene and Plio-Pleistocene.

Examination of hemiplasy, homoplasy and phylogenetic discordance in chromosomal evolution of the Bovidae.

Robertsonian chromosomal fusions predominate in shaping the genomes of many species of Bovidae. These and other cytogenetic data (from 52 taxa representing 51 species and 9 tribes of Bovidae) were (i) examined for usefulness in defining phylogenetic relationships and (ii) subsequently mapped to a consensus tree based on mitochondrial and nuclear DNA gene sequences with divergence dates of the corresponding species calculated from cytochrome b sequences. This permitted persistence time estimates for the various rearrangements. The chromosomal data resulted in an unsupported higher-level topology, but with recognition of the monophyly of some genera and tribes within Bovidae. The distribution and temporal spread of character states on the species tree is suggestive of a restricted role for hemiplasy (the retention of an ancestral chromosomal polymorphism through multiple speciation events) and for introgression (resulting from secondary contact among taxa), processes that can potentially lead to phylogenetic discordance. We conclude that the most probable interpretation for these data is that genuine karyotypic homoplasy predominates, but that hemiplasy (and/or introgression) is a realistic hypothesis for the observed patterns of several shared characters in Bovidae.

A paradox revealed: karyotype evolution in the four-horned antelope occurs by tandem fusion (Mammalia, Bovidae, Tetracerus quadricornis).

The four-horned antelope, Tetracerus quadricornis, is a karyotypic novelty in Bovidae since chromosomal evolution in this species is driven by tandem fusions in contradiction to the overwhelming influence of Robertsonian fusions in other species within the family. Using a combination of differential staining and molecular cytogenetic techniques, we provide the first description of the species' karyotype, draw phylogenetic inferences from the cytogenetic data and discuss possible mechanisms underlying the formation of the tandem fusions in this species. We show (a) that pairs 1-6 of Tetracerus correspond to a combination of Bos taurus orthologous chromosomes that are tandemly fused head to tail, (b) the presence of interstitial centromeric satellite DNA at the junctions of orthologous blocks defined by the cross-species painting data and (c) that in some instances, residual telomeric sequences persist at these sites. We conclude that the attendant result of each fusion is an enlarged acrocentric fusion element comprising a single functional centromere and two terminal telomeres that, collectively, led to a reduction of the 2n = 58 bovid ancestral acrocentric chromosomal complement to the 2n = 38 detected in the four-horned antelope.

SuperTRI: A new approach based on branch support analyses of multiple independent data sets for assessing reliability of phylogenetic inferences.

Supermatrix and supertree are two methods for constructing a phylogenetic tree by using multiple data sets. However, these methods are not a panacea, as conflicting signals between data sets can lead to misinterpret the evolutionary history of taxa. In particular, the supermatrix approach is expected to be misleading if the species-tree signal is not dominant after the combination of the data sets. Moreover, most current supertree methods suffer from two limitations: (i) they ignore or misinterpret secondary (non-dominant) phylogenetic signals of the different data sets; and (ii) the logical basis of node robustness measures is unclear. To overcome these limitations, we propose a new approach, called SuperTRI, which is based on the branch support analyses of the independent data sets, and where the reliability of the nodes is assessed using three measures: the supertree Bootstrap percentage and two other values calculated from the separate analyses: the mean branch support (mean Bootstrap percentage or mean posterior probability) and the reproducibility index. The SuperTRI approach is tested on a data matrix including seven genes for 82 taxa of the family Bovidae (Mammalia, Ruminantia), and the results are compared to those found with the supermatrix approach. The phylogenetic analyses of the supermatrix and independent data sets were done using four methods of tree reconstruction: Bayesian inference, maximum likelihood, and unweighted and weighted maximum parsimony. The results indicate, firstly, that the SuperTRI approach shows less sensitivity to the four phylogenetic methods, secondly, that it is more accurate to interpret the relationships among taxa, and thirdly, that interesting conclusions on introgression and radiation can be drawn from the comparisons between SuperTRI and supermatrix analyses.

Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae).

Organisms living at high altitude are exposed to severe environmental stress associated with decreased oxygen pressure, cold temperatures, increased levels of UV radiation, steep slopes, and scarce food supplies, which may have imposed important selective constraints on the evolution of the mitochondrial genome. Within mammals, the tribe Caprini is of particular interest for studying the evolutionary effects of life at high altitude, as most species live in mountain regions, where they developed morphological and physiological adaptations for climbing. In this report, we analyzed the complete mitochondrial genome of 24 ruminants, including 20 species of Caprini. The phylogenetic analyses based on 16,117 nucleotides suggested the existence of a new large clade, here named subtribe Caprina, containing all genera, but Pantholops (Pantholopina), Capricornis, Naemorhedus, and Ovibos (Ovibovina). The alignment of the control region showed that all Caprini have between two and four tandem repeats of ~75 bp in the RS2 region, and that several of these copies emerged from recent and independent duplication events. We proposed therefore that the maintenance of at least two RS2 repeats in the control region of Caprini is positively selected, probably for producing a higher number of D-loop strands 3'-ending at different locations. The analyses of base composition at third-codon positions of protein-coding genes revealed that Caprini have the highest percentage of A nucleotide and the lowest percentage of G nucleotide, a pattern which suggests increased rates of cytosine deamination (C-->T transitions) on the H strand of mtDNA. Two nonexclusive hypotheses related to high-altitude life can explain such a mutational pattern: more severe oxidative stress (ROS) and higher metabolic rates. By comparing the relative rates of nonsynonymous and synonymous substitutions in protein-coding genes, we identified that Caprini have higher levels of adaptive variation in the ATPase complex. In addition, we detected several changes in mitochondrial genes that should be tested for their potential role in mountain adaptation.

Chromosome evolution in the subtribe Bovina (Mammalia, Bovidae): the karyotype of the Cambodian banteng (Bos javanicus birmanicus) suggests that Robertsonian translocations are related to interspecific hybridization.

Three subspecies of banteng (Bos javanicus) have been described: B. j. javanicus in Java, B. j. lowi in Borneo, and B. j. birmanicus in Cambodia, Lao PDR, Myanmar, Thailand and Vietnam. In this paper we provide the first description of the karyotype of the Cambodian banteng. The chromosomal complement of B. j. birmanicus differs from that of B. j. javanicus, which was previously found to be similar to that of cattle, Bos taurus (2n = 60). The Cambodian banteng karyotype has a diploid number of 2n = 56 (FN = 62) and the karyotype consists of 26 pairs of acrocentric chromosomes and two pairs of submetacentric chromosomes. Comparisons with other species of the subtribe Bovina show that the two pairs of bi-armed chromosomes resulted from two centric fusions involving the equivalent of cattle chromosomes 1 and 29, and 2 and 28, respectively. Cross-species fluorescence in-situ hybridization (FISH) with B. taurus whole chromosome paints and satellite DNA I probes was used to identify the chromosomes involved in the translocations, and their orientation. We suggest that Robertsonian translocations (1;29) and (2;28) have been fixed in the common ancestor of Cambodian banteng as a consequence of hybridization with the kouprey (Bos sauveli) during the Pleistocene epoch.

Resolving a zoological mystery: the kouprey is a real species.

The kouprey is a rare and enigmatic forest ox discovered by scientists in Cambodia only in 1937. Numerous morphological hypotheses have been proposed for the origin of the kouprey: that it is a species closely related to banteng and gaur, two other wild oxen of southeast Asia; a morphologically divergent species placed in a separate genus, named Novibos; a wild species linked to aurochs and domestic cattle; a vicariant population of banteng; a feral cattle; or a hybrid of banteng with either zebu cattle, gaur or water buffalo. In a recent paper, which gained a lot of media coverage, Galbreath et al. analysed mitochondrial DNA sequences and concluded that the kouprey never existed as a wild, natural species, and that it was a feral hybrid between banteng and zebu cattle. Here we analyse eight DNA markers-three mitochondrial regions and five nuclear fragments-representing an alignment of 4582 nucleotides for the holotype of the kouprey and all related species. Our results demonstrate that the kouprey is a real and naturally occurring species, and show that Cambodian populations of banteng acquired a mitochondrial genome of kouprey by natural introgressive hybridization during the Pleistocene epoch.

Mitochondrial DNA variability in Giraffa camelopardalis: consequences for taxonomy, phylogeography and conservation of giraffes in West and central Africa.

The giraffe (Giraffa camelopardalis) still survives in four countries of West and central Africa. The populations of Niger and Cameroon are generally assigned to the subspecies peralta, but those of Chad and the Central African Republic are taxonomically problematic, as they are referred to as either peralta, or antiquorum, or congoensis. In this study, a mitochondrial fragment of 1765 nucleotide sites, covering the complete cytochrome b gene, three transfer RNAs and a large part of the control region, was sequenced to assess the relationships between several populations of giraffe. The phylogenetic analyses performed on the 12 identified haplotypes indicate that northern giraffes constitute a natural group, distinct from that of southern giraffes. Surprisingly, the giraffes of Niger are found to be more closely related to the giraffes of East Africa (subspecies rothschildi and reticulata) than to those of central Africa. We conclude therefore that the subspecies peralta contains only the Niger giraffes, whereas the subspecies antiquorum includes all populations living in Cameroon, Chad, the Central African Republic, and southwestern Sudan. We suggest that the ancestor of the Nigerian giraffe dispersed from East to North Africa during the Quaternary period and thereafter migrated to its current Sahelian distribution in West Africa, in response to the development of the Sahara desert. This hypothesis implies that Lake Mega-Chad acted as a strong geographical barrier during the Holocene, preventing any contact between the subspecies peralta and antiquorum. Our study has direct implications for conservation management, as we show that no subspecies peralta is represented in any European zoos, only in Niger, with a small population of less than 200 individuals.

Hybrid origin of the Pliocene ancestor of wild goats.

Recent theories on speciation suggest that interspecific hybridization is an important mechanism for explaining adaptive radiation. According to this view, hybridization can promote the rapid transfer of adaptations between different species; the hybrid population thus invades new habitats and diversifies into a variety of new species. Although hybridization is well accepted as a fairly common mechanism for diversification in plants, its role in the evolution of animals is more controversial, because reduced fitness would typically condemn animal hybrids to an evolutionary dead-end. Here, we examine DNA sequences of four mitochondrial and four nuclear genes selected for resolving phylogenetic relationships between goats, sheep, and their allies. Our analyses provide evidence of strong discordance for the position of Capra between mitochondrial and nuclear phylogenies. We suggest that the common ancestor of wild goats arose from interspecific hybridization, and that the mitochondrial genome of a species better adapted to life at high altitudes was transferred via this route into the common ancestor of Capra. We propose that the acquisition of more efficient mitochondria has conferred a selective advantage on goats, allowing their rapid adaptive radiation during the Plio-Pleistocene epoch. Our study therefore agrees with theories that predict an important role for interspecific hybridization in the evolution and diversification of animal species.

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.