Skip to main content
  • Research article
  • Open access
  • Published:

A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics

Abstract

Background

Molecular phylogenetics has provided unprecedented resolution in the ruminant evolutionary tree. However, molecular age estimates using only one or a few (often misapplied) fossil calibration points have produced a diversity of conflicting ages for important evolutionary events within this clade. I here identify 16 fossil calibration points of relevance to the phylogeny of Bovidae and Ruminantia and use these, individually and together, to construct a dated molecular phylogeny through a reanalysis of the full mitochondrial genome of over 100 ruminant species.

Results

The new multi-calibrated tree provides ages that are younger overall than found in previous studies. Among these are young ages for the origin of crown Ruminantia (39.3–28.8 Ma), and crown Bovidae (17.3–15.1 Ma). These are argued to be reasonable hypotheses given that many basal fossils assigned to these taxa may in fact lie on the stem groups leading to the crown clades, thus inflating previous age estimates. Areas of conflict between molecular and fossil dates do persist, however, especially with regard to the base of the rapid Pecoran radiation and the sister relationship of Moschidae to Bovidae. Results of the single-calibrated analyses also show that a very wide range of molecular age estimates are obtainable using different calibration points, and that the choice of calibration point can influence the topology of the resulting tree. Compared to the single-calibrated trees, the multi-calibrated tree exhibits smaller variance in estimated ages and better reflects the fossil record.

Conclusions

The use of a large number of vetted fossil calibration points with soft bounds is promoted as a better approach than using just one or a few calibrations, or relying on internal-congruency metrics to discard good fossil data. This study also highlights the importance of considering morphological and ecological characteristics of clades when delimiting higher taxa. I also illustrate how phylogeographic and paleoenvironmental hypotheses inferred from a tree containing only extant taxa can be problematic without consideration of the fossil record. Incorporating the fossil record of Ruminantia is a necessary step for future analyses aiming to reconstruct the evolutionary history of this clade.

Background

Phylogenetic studies based on genomic data provide an unparalleled tool for the reconstruction of phylogenetic relationships among organisms. However, considering that the overwhelming majority of species that have ever existed on Earth are now extinct, attempts to reconstruct evolutionary history must necessarily consider the fossil record. A key challenge of contemporary phylogenetic methods, therefore, is to combine data from the fossil record with genomic data from living organisms [1]. One of the main contributions of the fossil record to molecular phylogenetics has been the provision of fossil calibrations for molecular rates of evolution that in turn allow for estimates of the timing and rate of cladogenetic events. Molecular age estimation is of particular importance for evolutionary events for which the fossil record is incomplete, ambiguous, or even totally unknown. However, the paleontological data used for molecular rate calibration are often misunderstood and misapplied, resulting in poor (and wildly varying) age estimates [17]. The problem stems in part from the opacity of paleontological literature to non-specialists, as well as to the fact that fossils are not always properly assigned to either the stem group or the crown clade of a taxon [7].

One clade that has been the subject of numerous studies of molecular age estimation is Bovidae, or the clade comprising antelopes, oxen, sheep, and goats. With over 130 living species, Bovidae is the most diverse clade of large mammals alive today, and includes almost half of all living species included in Artiodactyla (= Cetartiodactyla). Bovids are ecologically and geographically diverse, spanning the Rocky Mountains to the African rainforest, and their fossil record is abundant and spans the last 18 million years. Interest in bovid phylogenomics is widespread, covering subjects as diverse as Neogene phylogeography to the reconstruction of human domestication practices in the late Holocene. Age-calibrated trees of Bovidae are also relevant to studies seeking to establish rates of evolution within a larger context of ruminant and mammalian evolution, including, for example, the terrestrial origins of whales e.g. [8, 9]. Though the fossil record of Bovidae and Ruminantia is extensive, difficulties of phylogenetic resolution mean that many topological relationships and phylogeographic reconstructions have come from molecular phylogenies of extant taxa. Similarly, many of the evolutionary divergences within Ruminantia have been dated mainly through molecular age estimates [8, 10, 11]. However, most of these studies have relied on a single or a few (often incorrectly applied) fossil calibration points that are often far removed from many of the nodes being dated. The large number and high diversity of species within ruminants means that molecular rate variations are a real possibility, and that the use of a larger number of fossil calibration points distributed throughout the tree should result in better age estimates that in turn better reflect the evolutionary history of the clade.

In an attempt to test this idea, I here identify 16 fossil calibration points relevant to the phylogeny of Bovidae and Ruminantia, and apply these, together and individually, to a recently published matrix of the complete mitochondrial genome of over 100 ruminant species [8]. Running each of the 16 calibration points alone in independent analyses provides a sense of the degree to which calibration choice can affect molecular age estimates as well as any effects on tree topology. Running all available calibrations together anticipates possible variations in molecular rates, and helps provide new age estimates reflecting a consensus of a larger portion of the fossil record. I also use the results of these analyses to highlight current areas of conflict between genomic and morphological approaches to evolutionary reconstruction, and to promote the incorporation of increasing amounts of paleontological information into molecular phylogenetic approaches.

In what follows, references to Pleistocene refer to the recently redefined Pliocene-Pleistocene boundary at 2.588 Ma (previously at 1.8 Ma) [12]. The term crown clade refers to a node-based clade originating with the last common ancestor of two or more extant species or organisms [13]. The stem group, or stem lineage, is the ancestral lineage leading to a crown group. The branch-based clade that comprises a crown clade plus the stem leading to it is the ‘total’ or ‘pan-’ clade [13]. I use the order name Artiodactyla (rather than Cetartiodactyla) as Cetaceans are nested within Artiodactyla and not its sister group (i.e. whales are artiodactyls too) [9, 14]. I follow long-standing classifications that recognize between about 130 and 145 extant bovid species [1518], rather than a recent treatment proposing 279 species [19].

Results

The 16 fossil calibrations used are listed in Table 1 and full details on each are given as Additional file 1. The maximum clade credibility (MCC) tree summarizing the trees found from the analysis using 16 calibration points is shown in Figure 1 (and Additional file 2: Table S1). The topology and divergence age estimates of this tree differ in a few respects from those of the recent analyses of Artiodactyla (= Cetartiodactyla) by Hassanin et al. [8], the source of the data matrix used in this study. MCC trees from the runs using each calibration point alone and the run using no calibration priors and random starting trees are available (with full metadata) online through TreeBase. Trees resulting from the different calibrated analyses show slight variations in topology (detailed below) indicating that choice of calibration node did have some effect on tree topology. The topology of the analysis using no calibrations and random starting trees is the same as those of some of the calibrated analyses, indicating that use of a starting tree did not have a significant effect on final tree topology.

Figure 1
figure 1

Tree resulting from Bayesian analysis of the full mitochondrial genome of 127 ruminants using 16 fossil calibration points (yellow circles). Calibrations along branches indicate stem calibrations. Node bars represent 95% intervals and node values are posterior probabilities. Names of terminal taxa follow those in ref. [8].

Table 1 Fossil calibration data

The median estimated molecular ages of the different nodes from all the different analyses is shown in Table 2, along with the mean and range of node estimates across all the analyses. Table 2 shows that using only a single calibration point can produce a very wide range of age estimates, with maximum ages being more than double the minimum ages obtainable. In general, older calibration points produced older age estimates, and younger (ingroup) calibration points produced younger age estimates. The use of 16 calibration points results in age estimates intermediate among those of the single-calibration analyses, with older nodes falling toward the mean of the range of estimates, and younger nodes (<10 Ma) falling towards the younger limit of the range of estimates (Figure 2). This may be a result of the age distribution of the 16 calibration points, ten of which are less than 10 Ma in age on the final tree (Figure 1). Node ages in the multicalibrated analysis are also more often associated with smaller error ranges, which was not necessarily to be expected with the use of a larger number of calibration points [55] (Figure 1). Not surprisingly, the multicalibrated tree displays higher rate heterogeneity (in Tracer, ucld.stdev values around 0.22) compared to the single-calibrated analyses (0.14-0.16).

Figure 2
figure 2

Median age estimates by node for each of the 16 single-calibrated analyses and the analysis using all 16 fossil calibrations (data in Table2), covering the tree root (far left) and the 16 calibrated nodes. A large range of ages is obtainable with the use of different calibration points. The multi-calibrated analysis (blue with circles) provides estimates that are intermediate in age for older nodes, and towards the younger end of the range for younger nodes (<10 Ma).

Table 2 Node age estimates (Ma)

Discussion

Fossil calibration points: single, few, or many?

Multiple fossil calibrations have already been argued as preferable to single calibrations, particularly given that rates of molecular evolution are not constant across a phylogeny, and that rate variations may not necessarily be autocorrelated among adjacent branches [5659]. The results here provide further support for the use of multicalibrated analyses and additionally illustrate the very wide range of age estimates that may be obtained when only a single calibration point is used.

Several authors have attempted to selectively choose points from an available pool of fossil calibrations by evaluating the internal congruence of molecular age estimates provided by these calibrations, both against each other and against the fossil record [6064]. While there may be interesting reasons to use the internal congruence of molecular age estimates from multiple calibration points to identify outliers, such an approach is ultimately a poor substitute for the careful selection of calibration points from the paleontological record [2, 6]. It is also potentially circular in reasoning (deriving a calibration point from the fossil record, the estimates from which are then assessed against the fossil record, and these then used to discard a given data point because it apparently doesn’t hold up to the fossil record). Provided a fossil’s age and phylogenetic position have been properly supported as detailed in ref. [2], a low molecular age congruence score should probably not be used to discount it as a molecular rate calibration point. There are numerous reasons why estimates of divergence times derived from molecular data could differ with choice of calibration point, and why molecular and fossil estimates of age might disagree. Among these are rate heterogeneity outside the statistical bounds of the uncorrelated lognormal clock e.g. large population size or generation time changes, [65, 66], node-density effects [67], ancestral polymorphism, and time-dependent rate variations that produce artificial branch length compression at older nodes (e.g. saturation) and extension towards the tips [63, 68, 69]. Ultimately, the use of prior probability distributions and soft bounds for calibration means that the effects of an erroneous fossil calibration within a multiple calibration analysis is limited [70]. Table 2 and Figure 2 show that the use of all 16 calibration points produces a tree that balances aspects of all the calibration points, and is not overly affected by some of the more extremely different points. As might also be expected, the multicalibrated tree based on all 16 calibrations is also more congruent with the predicted ages of the calibrated nodes themselves (lowest sum of percent error of the estimate at each of the 16 nodes with respect to the median age of the probability distribution, Additional file 3: Table S2). Using all available and vetted fossil calibration points therefore seems a better approach than relying on assumptions of molecular clock internal congruence to discard what might be perfectly good information from the fossil record.

New age estimates for the phylogeny of Bovidae and Ruminantia

Many of the calibration points used separately provide very young ages for the entire phylogeny, including ages for certain nodes that are even younger than the minimum age as established from the fossil record (Table 2). This is an interesting result given that molecular rate estimates are very often found to be older than expectations from the fossil record [71]. Molecular estimates of a node’s age are actually ranges defined by probability distributions, and these can vary greatly depending on the type of genomic data used e.g., [63], the choice and number of calibration points (e.g. this study), choice and number of taxa [72], and choice of the many (ever increasing) parameters of a phylogenetic analysis e.g., [73]. Therefore, the molecular age estimates produced by different analyses should be considered as hypotheses open to further testing by the fossil and genomic records. In the discussion that follows I identify node ages using their 95% probability ranges. I refer extensively to the results of the study by Hassanin et al. [8], which used six calibration points from the artiodactyl tree, with a slightly older crown Bovidae calibration among these. Though I used a pruned version of their same genomic dataset, the analysis methods differ greatly between the two studies. References to ages in the study by Hassanin et al. [8] use the mean ± 2 standard deviation values of their BD-Soft model results. References to ages in the study by Meredith et al. [59] use their DNA Soft IR model.

The new multi-calibrated phylogeny (Figure 1), provides divergence age estimates that are younger overall (~1 to 4 Ma) than those of Hassanin et al. [8], from which the current genomic dataset was derived. These authors used a similar but older crown Bovidae calibration in addition to five other calibrations from other parts of the artiodactyl tree, as well as different software for estimation of divergence times. The largest differences are seen with age estimates towards the base of the tree rather than in the ingroup, which might in part result from the use here of a much reduced taxonomic subset, in addition to a larger number of younger calibrations that produce younger dates (Table 2). Lineage-specific rate variations may also have some effect, as two out of six of Hassanin et al.’s calibration points were close to the origin of Cetacea, which as a clade has been shown to have slow rates of molecular evolution relative to other mammals, potentially inflating estimates of divergence times in other parts of the tree [59]. It is worth noting, however, that age estimates from the analysis using only the crown Bovidae calibration (Table 2) are, with the exception of the root node, not greatly different from the age estimates found by Hassanin et al., indicating that the older dates of their study could be replicated here even with a reduced taxonomic subset and a different calibration scheme.

The multicalibrated tree gives a root age (origin of crown Ruminantia) at 39.3–28.8 Ma (95% range, median 33.7 Ma). This is younger than Hassanin et al.’s [8] age of 51.6 ± 9.8 Ma, and a better match for the 41.4–35.1 Ma estimate by Meredith et al. [59]. These younger estimates at first seem too young given that the fossil record of fossil ‘ruminants’ extends back to the middle Eocene, to around 45 Ma [74]. However, the earliest fossil (probably stem) tragulids are in fact only around 33 Ma in age [74, 75]. If Leptomerycidae or Lophiomerycidae are accepted as belonging within crown Ruminantia [9, 74, 76, 77], then that would require the crown clade to be at least 42 Ma in age [74], but these and other fossil ruminant lineages might alternately be found to lie on the stem group. O’Leary et al. [14] also provide support for younger dates. Though they did not sample multiple ruminant taxa, they date the origin of Cetruminantia (ruminants plus hippos and whales) to a minimum divergence of 50 Ma, close to Meredith et al.’s 55 ± 2 Ma, and significantly younger than Hassanin et al.’s estimates of 69.7 ± 9.4 Ma.

The age of origin found here for crown Pecora, and Moschidae in particular, is also young compared to previous sources, including the fossil record. According to my results, crown Pecora originated 23.0–19.4 Ma, and the lineage leading to Moschidae diverged from that leading to Bovidae at 19.3–16.6 Ma. This is a match for Meredith et al.’s [59] estimates of 21.0–18.2 Ma for Pecora and 19.1–16.4 Ma for Bovidae + Moschidae, and significantly younger than Hassanin et al.’s [8] estimates of 27.6 ± 7.6 Ma and 22.4 ± 4.8 Ma. However, if the earliest stem moschids are late Oligocene in age (>23 Ma), this means the origin of crown Pecora must have taken place by this time. The much younger age returned by the analysis is interesting considering that the earliest cervid, moschid, and bovid fossils were used as calibration points, and that, when run separately, these three calibrations provide some of the oldest age estimates among all analyses (Table 2). In contrast, nine of the single-calibrated analyses return very young estimates (<16 Ma median age) for the Moschidae-Bovidae split. While the possibility of a younger (c.16 Ma) origin of crown Bovidae could be considered (below), a 19.3–16.6 Ma split for Moschidae-Bovidae is in direct conflict with our current understanding of the fossil record. Resolving this conflict will not be straightforward. Stretching the branch length leading to Moschidae back into the late Oligocene would create long ghost lineages for the other pecoran families. A more basal phylogenetic placement of Moschidae as the first family to diverge in the pecoran radiation would fit the order of appearance of pecoran families in the fossil record, but is not supported by molecular analyses. An alternate, perhaps more likely, possibility is that the assignment of Dremotherium to the stem group of Moschidae is incorrect and that the oldest stem moschids are the early Miocene Pomelomeryx or the middle to late Miocene Micromeryx[78, 79]. In this case, the divergence of Moschidae from the remaining Pecora could be as young as the molecular data here indicate. Future phylogenetic analyses, particularly on the fossil record of early moschids, may help resolve the conflict.

An early Miocene age for the origin of Pecora, however, does fit the fossil record of Antilocapridae, Giraffidae, Cervidae, and Bovidae, the oldest (stem) fossils of which appear between about 23 and 16 Ma [80, 81]. My estimates for a late Miocene origin of crown Cervidae (12.4–9.0 Ma) match those of Hassanin et al. (10.7 ± 2 Ma) [8] and are slightly older than those of Gilbert et al. (9.6–7.7 Ma) [82]. In separate analyses, Gilbert et al. [82] also obtained much older ages for crown Cervidae, but this result used early Miocene cervids to calibrate the node of origin of crown Cervidae, instead of the origin of the stem group as done here (Additional file 1). A middle to late Miocene origin of crown Cervidae fits the fossil record, in which Miocene species of Euprox, which may lie close to the origin of the crown node, are dated to c.14–11 Ma or later [21, 24].

The young age found here for the origin of crown Bovidae (17.1–15.3 Ma) also contrasts with the older results of most studies to date, but may better fit the fossil record. The oldest fossil bovid is Eotragus noyei from 18.3 Ma in Pakistan [23], but since it is not clear whether this and other Eotragus species belong on the stem lineage or in the crown clade of Bovidae, there is no reason to believe that crown bovids must have appeared by 18 Ma (see Additional file 1). The middle Miocene bovids Pseudoeotragus seegrabensis and Eotragus sansaniens have been proposed as some of the earliest representatives of stem Antilopinae and stem Bovinae, respectively [83, 84]. Other early Bovinae or stem Bovinae include Miocene ‘boselaphines’ such as Kipsigicerus labidotus, Protragocerus spp., Helicoportax spp., and Sivaceros spp, and among early Antilopinae or stem Antilopinae, Caprotragoides potwaricus, Tethytragus spp., Gentrytragus spp., and early forms assigned to ‘Gazella’ [25, 83, 8588]. This pool of taxa indicates that the divergence of Bovinae and Antilopinae probably took place during the middle Miocene. A single site of old age that is well dated and at which several of the above mentioned taxa occur is Fort Ternan in Kenya [86]. The uppermost level bearing mammalian (including bovid) remains in the Fort Ternan “B” beds has an age of 13.9 ± 0.3 Ma [89]. A young (17–15 Ma) age of origin of crown Bovidae is therefore a reasonable hypothesis. An important point here is that the phylogenetic position of the oldest fossil bovid [23] is indeterminate relative to the crown bovid node. Such early bovids may represent stem taxa that predate the last common ancestor of extant bovids, and so therefore should not be used to provide minimum ages for crown Bovidae, as many studies to date have done, therefore potentially overestimating the age of this node.

Bovid phylogeny: topology of the tree using 16 calibrations

The MCC tree resulting from the analysis using all 16 fossil calibration points (Figure 1) and those arising from the 16 single-calibration analyses and uncalibrated analysis are almost all identical. Trees differing slightly in topology to that in Figure 1 were found with 9 of the 16 single-calibrated runs (stem Moschidae, stem Cervidae, stem Hippotragini, crown Hippotragini, crown Tragelaphini, stem Caprini, crown Connochaetes, crown Kobus, stem T. strepsiceros), as well as with the analysis using no calibrations at all (Additional file 1). The main differences observed in these phylogenies involve a rearrangement at the base of the Pecoran radiation (recovery of an Antilocapra + Giraffidae clade and the placement of Cervidae basal to it), and inverted relationships among some of the bovid tribes (Reduncini + Pelea closer to Antilopini, and Cephalophini + Oreotragus closer to the Caprini + Alcelaphini + Hippotragini clade). Some of these relationships are surrounded by poor support values and therefore probably indicate areas of the tree requiring further investigation. These results also indicate that choice of fossil calibration can produce variations in tree topology.

Despite the use of very different methodological approaches, the trees presented here are similar in topology to the maximum likelihood and Bayesian trees of Hassanin et al. [8]. Specifically, Hassanin et al. used 36 nucleotide partitions, and a supertree approach to combine the results of repeated maximum likelihood analyses on overlapping subsets of the matrix. These authors also performed a Bayesian analysis, but using different software and a different evolutionary model (CAT - GTR + Gamma 4). My tree differs from the Bayesian tree of Hassanin et al. (appendix 5 in ref. [8]) and agrees with their maximum likelihood tree (Figure one in ref. [8]) in three main aspects: 1) a resolved basal branching order of Tragulidae - > Antilocapra - > Giraffidae - > Cervidae - > Moschidae (all PP = 1) rather than a 3-way polytomy with Antilocapra + Giraffidae, Cervidae, and Moschidae + Bovidae cf. ref. [59]; 2) a sister relationship between Alcelaphini and Hippotragini (PP = 0.84), rather than Alcelaphini and Caprini; 3) placement of Reduncini + Pelea as closer to Caprini + Alcelaphini + Hippotragini (PP = 0.71) rather than Antilopini; and 4) placement of Cephalophini + Oreotragus as sister to Antilopini (PP = 0.72), rather than basal to Caprini + Alcelaphini + Hippotragini (there are variations on 1, 3, and 4 however, in some of the single-calibration runs mentioned above). Naturally, all these relationships would benefit from further testing, especially at nodes with low posterior probability values (PP <0.9). Topological differences towards the base of the tree especially might also be affected by taxon sampling, as the current study only considers a subset of the taxa analyzed by Hassanin et al..

The current MCC tree (Figure 1) differs from the maximum likelihood tree presented by Hassanin et al. [8] and agrees with their Bayesian tree and also ref. [90] in the placement of the saola (Pseudoryx) as the sister taxon to Bovini rather than within Bovini (PP = 1), and also in the placement of the rhebok (Pelea) as the sister taxon to Reduncini and not within Reduncini (PP = 1). These are significant results, as they provide support for the monophyly of Bovini and Reduncini as traditionally defined on the basis of major morphological differences, and highlight an area of conflict between morphological and genomic approaches to taxonomic classification.

Genomic vs morphological taxonomy

Over the last 20 years, bovid taxonomy has received much needed revision and clarification as a result of molecular phylogenetic investigations. However, conflicts between genomic and phenotypic (e.g. morphological, behavioral, ecological) delimitations of species and higher taxa persist in different parts of the bovid tree. Presumably on the basis of low molecular distances (short molecular branch lengths), some workers have proposed expanding the composition of traditionally defined Reduncini, Bovini, and Caprini (formerly Caprinae) to include their nearest sister taxa, the saola (Pseudoryx) with Bovini, the rhebok (Pelea) with Reduncini, and the chiru (Pantholops) with Caprini e.g. Figure 1 in ref. [8]. This is in conflict with traditional classifications that distinguish these tribes from their sister taxa on the basis of large suites of advanced phenotypic characters [9193]. In all three cases, there is an apparent lack of correlation between the rate and amount of molecular and morphological evolution that has taken place e.g. [94]. Perhaps homoplasy, whether of nucleotide substitutions or as a result of the convergence or reversal of morphological characters, could be further investigated as a potential cause for the discrepancies between molecular and morphological distances at these nodes.

Since there are no agreed upon rules to govern the delimitation of monophyletic higher taxa, a consensus may be difficult to achieve, but at minimum, independent lines of evidence should be considered together. Higher taxa are probably only necessary in so far as they are useful to both specialists and non-specialists alike. In the case of Bovini for example, it is much more practical and intuitive for paleontologists, ecologists, zoologists, and lay people to restrict composition of this clade to the buffaloes and oxen, which are easily recognized bovids of large size with a significant suite of skeletal, dietary, and behavioral adaptations, and broadly similar ecological niches (Figure 3). The saola, in contrast, is a small forest dweller lacking most if not all of the derived features characterizing Bovini, so much so that its skull was found to resemble caprins such as Capricornis and Naemorhedus[95], but see [96]. The fact that DNA differences between the saola and Bovini are few should probably not be used as the sole criterion for the delimitation of Bovini, otherwise the tribal name risks losing much of its coherency and common utility. Granted, expanding the composition of Bovini to inclide the saola would highlight the surprising genomic similarities among all these animals, but it would do so in spite of large phenotype-level differences, which clearly the genomic sequences alone do not predict. The same case applies to Pantholops with respect to Caprini, and Pelea with respect to Reduncini. The latter case is particularly worthy of further investigation because though the branch length (and resulting time estimate) between the divergence of Pelea and the origin of Reduncini is very short (Figure 1), and the support for reduncin monophyly is poor (PP = 0.71), reduncin fossils of Kobus porrecticornis as old as 9.3 Ma are recorded in the Siwaliks of Pakistan [29]. Further study of the Siwalik reduncins with the intent of establishing their relationships to living Redunca, Kobus, and Pelea could therefore go a long way toward resolving the unstable molecular arrangements at the base of Reduncini.

Figure 3
figure 3

A large number of morphological, ecological, and behavioral characteristics distinguish Bovini (A-C) from their sister taxon the saola (D, Pseudoryx nghetinhensis ). A, European bison (Bos (Bison) bonasus), B, African buffalo (Syncerus caffer), and C, domestic cattle (Bos taurus). All images from Wikimedia Commons, user credits: Michael Gäbler (A), Haplochromis (B), DFoidl (C), Silviculture (D).

Conflict is also present with regard to taxonomic revisions made at the species level. Relationships based on mitochondrial genomes benefit from independent confirmation by nuclear genes e.g. [97]. For example, the mitochondrial sister relationships of Capra sibirica to Hemitragus jemlahicus and Bos (Bison) bonasus to Bos taurus, which result in the polyphyly of Capra and the bisons, are contradicted by both nuclear genomic and morphological analyses that both support the monophyly of these two taxa [97102]. Analyses such as those presented here or by Hassanin et al. [8] similarly support previous findings [103] regarding mitochondrial polyphyly within the bushbuck (here labelled T. scriptus 1 and T. scriptus 2 + 3 in Figure 1). However, these mitochondrial relationships contradict morphological evidence, and, in light of the examples of the bisons and Capra, the elevation of T. scriptus sylvaticus to the level of a separate species (T. sylvaticus) unrelated to T. scriptus should await further confirmation from the nuclear genome [103]. The possibility of incomplete lineage sorting, or mitochondrial lineage introgression among geographically proximate tragelaphin species, as probably occurred within Caprini and Bovini, has not yet been ruled out.

Neogene diversity dynamics and biogeography of Bovidae

It is tempting to reconstruct the timing and mode of interesting evolutionary, ecological, and biogeographic events on a well-resolved tree of extant taxa. However, doing so in the absence of the fossil record can be highly misleading, as a lineage’s characteristics (e.g. rates of extinction, geographic distribution, habitat preferences) may have changed significantly through the course of its evolution. I provide a few examples of evolutionary hypotheses developed from the dated phylogeny in Figure 1, and compare these to the fossil record.

The bovid subclades (tribes) Bovini, Tragelaphini, Cephalophini, Reduncini, Antilopini, Alcelaphini, Hippotragini, and Caprini comprise some 120 of ~140 living bovids species. Rather than being entirely taxonomic constructs, these bovid tribes comprise large and distinct ecological guilds, each distinguished by suites of unique morphological and behavioral characters that reflect distinct dietary and habitat preferences. These tribal-ecological associations have proved so compelling that paleontologists regularly use the proportions of bovid tribes at fossil sites as an indicator of the range of vegetational habitats found at that site [104107]. The age of origin of the main ecologically distinct clades (or ‘guilds’) of bovids is therefore of interest from the perspective of the development of modern vegetational biomes, and particularly in Africa, where bovid diversity mostly resides today. According to the dated tree (Figure 1), the origins of Reduncini, Caprini, Alcelaphini, Hippotragini, Cephalophini, and Tragelaphini all took place between 10.1 and 5.4 Ma (maximal 95% interval, or just 9.2–6.3 Ma using median ages), which is a short interval of time considering the 15 to 18 million year history of Bovidae. However in almost all cases, the stem branches leading to these originations are long, stretching back to the middle Miocene, which makes it possible that their late Miocene ‘radiation’ is actually an artifact of differential (non-random) extinction of older lineages. Given that the majority of middle and early late Miocene fossil bovids have not been well studied from a phylogenetic perspective, and given limitations of testing for deep extinction using only extant species, evaluating punctuated (radiation) versus gradualistic evolution (differential extinction) scenarios for the origins of the major ecological subclades (tribes) of Bovidae will benefit greatly from further work on the Miocene fossil record.

The origins of Bovini and Antilopini are older than those of the other bovid tribes (Figure 1). Antilopini is a highly diverse clade that is complexly distributed in terms of geography, such that any single calibrated phylogeny of this clade is full of evolutionary, biogeographic, and ecological hypotheses that could be tested with the fossil record. For example, the springbuck (Antidorcas) and the gerenuk (Litocranius) are African antilopins adapted to and living in relatively arid habitats. They have totally non-overlapping ranges, the springbuck being restricted to southwestern Africa and the gerenuk to the African Horn. The 8.8–6.4 Ma divergence between these two arid-adapted species might be interpreted as evidence for the late Miocene presence of arid corridors linking the Somali and South West arid zones, much as has been proposed for times in the Pliocene and Pleistocene [108]. The problem with such a scenario is that the restriction of Antidorcas to southwestern Africa appears to be a relatively recent phenomenon — the genus is recorded from the late Pliocene of Chad [109] and from the late Pliocene to the Pleistocene in eastern Africa, up until 2.5 Ma or younger in the Afar [110, 111], the early Pleistocene in the Turkana Basin [51, 112], and the middle Pleistocene in northern Tanzania [50, 51]. While no fossil record is known for the gerenuk [51], the wide range of Antidorcas throughout most of its history means Pliocene and Miocene phylogeographic inferences based on its modern day geographic range must be made with caution.

In another example, the blackbuck (Antilope cervicapra) is historically restricted to southern Asia (Pakistan, Nepal, India, and Bangladesh), but is overwhelmingly bracketed by African species (Figure 1) see also ref. [113]. The age of divergence of the blackbuck lineage from the remaining African Antilopini might therefore be taken to indicate a dispersal event from Africa into southern Asia sometime around 7.4–5.2 Ma. Critically missing from the phylogeny, however, is a large number of spiral-horned antilopin species from the late Miocene Greco-Iranian region, specifically Prostrepsiceros spp., from which the blackbuck may be descended [51, 114], and which indicate the blackbuck’s origins are probably to be found in Eurasia.

In a similar position, the Reduncini + Pelea clade is today entirely African in distribution, and phylogeographic inferences based on a tree of extant taxa might suggest that the origins of the Reduncini are to be found in Africa as well. The fact is, however, that the earliest reduncins known are from Pakistan [25, 28, 29, 115], raising the possibility that the morphological and ecological characteristics of this African clade may have their evolutionary roots in southern Asia, and not Africa.

In Bovini, a deep divergence of African (Syncerus) and Asian (Bubalus) buffaloes around 7 Ma (8.5–5.7 Ma 95%) might reflect findings from the fossil record. Around 8–7 Ma is the age of hypothesized dispersal of Bovini from southern Asia into Africa in the late Miocene based on the fossil record [26], and 7 Ma is also the age of the oldest Bovini currently known from Africa [116, 117]. Furthermore, there is evidence that faunal connections between southern Asia and Africa were stronger in the Miocene than today, becoming progressively weaker over time, and especially by the end of the late Miocene [115, 118]. The disruption of gene flow between Bubalus and Syncerus at around 7 Ma might also be corroborated by evidence for the evolution of Bubalus from Proamphibos and Hemibos in southern Asia, and Syncerus from Ugandax in Africa since the middle Pliocene or earlier [25, 51, 119], and may therefore reflect the strengthening of dispersal barriers between what may once have been relatively permeable Afrotropical and Indo-Malayan faunal zones [115, 120]. These hypotheses would of course benefit from further paleontological investigation, particularly of the spotty late Miocene African record.

Conclusions

In a recent paper, Bibi and Vrba [7] proposed that molecular phylogenetic studies of Bovidae could be advanced with the use of better calibration practices. I have here identified sixteen fossil calibration points that can be used as molecular age calibrations or reference points for phylogenetic analyses of Bovidae, Ruminantia, and Artiodactyla. I have shown that the choice of calibration points used in an analysis can significantly affect the resulting molecular age estimates for nodes all across the tree, as well as produce slight variations in topology. As such, incorporation of all vetted calibration data available is argued to be a better approach than using just one or a few calibrations. An analysis of the complete mitochondrial genome using a large number of fossil calibration points here produces younger ages for Bovidae (and Ruminantia), which may be a better fit to the fossil record than the older estimates of previous studies. The resulting tree has also been used to highlight areas of conflict between phylogeographic and evolutionary inferences derived from trees containing only extant taxa and those derived from the fossil record. Further development of nuclear genomic datasets and further incorporation of the fossil record are required to address central questions on rates of genomic versus morphological evolution, gene and species trees, molecular and fossil ages, molecular and morphological differences, and the recognition and timing of ecological radiations and dispersal events in Ruminantia. It is hoped that the inclusion of fossil taxa into a combined morphological-molecular analysis, along with the use of multiple molecular clock calibration points and perhaps different methodological approaches e.g. [73, 121] can bring us closer to a more complete evolutionary history of ruminants.

Methods

I used a subset of the matrix of complete mitochondrial genome matrix presented by Hassanin et al. [8], comprising 112 extant bovids, two moschids, six cervids, three giraffids, one antilocaprid, and two tragulids (these counts include in a few cases more than one individual of the same species). The matrix and analysis parameters were assembled in Beauti v.1.7.4 [122] and the analysis was run in BEAST 1.7.4 [122] using the online computational service Bioportal [123].

The mitochondrial dataset was analyzed as a single partition, using a lognormal relaxed clock with uncorrelated rates [122], the GTR + G + I model as determined by both AIC and BIC criteria in jModelTest 2.1.1 [124, 125], with 4 gamma categories, a Yule process speciation tree prior, and empirical base frequencies. A starting tree was used based on that of Hassanin et al.’s Bayesian analysis (Appendix 5 in ref. [8]), pruned down to the 127 taxonomic units used in this analysis, with polytomies arbitrarily resolved, and with branch lengths adjusted to fit prior age distributions (required for BEAST to run with so many priors). An independent analysis of the same dataset but including no calibrations and using a random starting tree was also run to test the effect of a starting tree on the final topology.

The 16 fossil calibrations used are listed in Table 1 and full details on each are given as Additional file 1 following the recommendations of Parham et al. [2]. Of the 16 points, three provide ages for nodes toward the base of the tree, and 13 provide ages for ingroup nodes (i.e. nodes inside Bovidae). These 16 calibrations are by no means exhaustive, and further study and research of the bovid record may find others, and even refine the ones currently used. While covering large parts of the bovid tree, there are conspicuous areas with no calibrations, namely Cephalophini, Antilopini, and Caprini. Cephalophini have nearly no fossil record at all [51]. Antilopini and Caprini on the other hand have a rich fossil record, but there is much ambiguity in the phylogenetic determination of most of the material. Further work should help clarify relationships, especially at the shallower nodes of both of these clades.

Several calibration points were considered but not used because in the current taxonomic sample they would point to the same node as another already used and older calibration point (e.g., stem Kobus, a duplicate of crown Reduncini; stem Alcelaphini, a duplicate of stem Hippotragini; stem Connochaetes, a duplicate of crown Alcelaphini; stem Taurotragus, duplicate of stem T. strepsiceros). Other excluded potential calibrations come from the fossil records of bisons, Capra, and Tragelaphus scriptus, given that all three taxa are found to be polyphyletic on the basis of mitochondrial genomic analysis (but not nuclear genes, discussed in text). In an isolated case, separate analyses (not shown) showed that use of a stem Reduncini calibration produced aberrant results, such as age probability distributions (error bars) much larger than in any other analysis, inclusion of Pelea within Reduncini (not found in any other analysis here), and a reversed branching order of Cervidae and Antilocapra with the nodes underlying these clades very weakly supported (found in a few other analyses, detailed in text). Including this data point in the multi-calibrated analysis affected the final topology much more so than with other calibration points, and this point was a result excluded.

All calibration points use a probability distribution with soft bounds [126]. The oldest fossil record in a clade provides a minimum age of origination for that clade [3, 71]. In cases where the phylogenetic position of the fossil is known with some certainty, a lognormal distribution is used with a 2.5% probability quantile (youngest bound) set at the age of the fossil. If the phylogenetic position of the fossil is not precisely known relative to the clade in question, a normal distribution is used allowing for the actual node age to be equally younger or older than the fossil. Maximal (97.5% probability quantile) node ages are difficult to determine from the fossil record and this setting is to a large degree unavoidably arbitrary. However, the use of multiple calibrations and soft bounds means accurately determining maximal bounds is less crucial than with a single calibration or hard bounds [61, 70]. Case-by-case assessments are given in the Additional file 1, but in general I chose conservatively large probability distributions. The use of soft bounds meant the analysis could still date a node to an age older or younger than the 95% probability interval.

No monophyletic constraints were used so as not to constrain tree topology. Tree root height was assigned a uniform distribution between 25 and 65 Ma with an initial starting value of 45. All analyses comprised 2 independent runs of 20,000,000 generations each, sampling every 1000 generations, and burn-in set to 2000 trees (10%) for each run. Examination of the resulting log files in Tracer [122] indicated convergence of all parameters in the independent runs and combined effective sample sizes greater than 95 for all values (and in most cases well above 200). Trees files from each run were combined using LogCombiner [122] and a maximum clade credibility (MCC) tree was created using Tree Annotator [122] with a burn-in of 4000 trees, using median heights and a posterior probability limit of 0.5.

Availability of supporting data

Full details on the sixteen fossil calibration points are given as Additional file 1 accompanying this manuscript online. The datasets (matrix and trees resulting from the phylogenetic analyses) supporting the results of this article are available in the TreeBase repository, http://purl.org/phylo/treebase/phylows/study/TB2:S14132.

References

  1. Donoghue PCJ, Benton MJ: Rocks and clocks: calibrating the tree of life using fossils and molecules. Trends Ecol Evol. 2007, 22: 424-431. 10.1016/j.tree.2007.05.005.

    PubMed  Google Scholar 

  2. Parham JF, Donoghue PCJ, Bell CJ, Calway TD, Head JJ, Holroyd PA, Inoue JG, Irmis RB, Joyce WG, Ksepka DT: Best practices for justifying fossil calibrations. Syst Biol. 2012, 61: 346-359. 10.1093/sysbio/syr107.

    PubMed Central  PubMed  Google Scholar 

  3. Reisz RR, Müller J: Molecular timescales and the fossil record: a paleontological perspective. Trends Genet. 2004, 20: 237-241. 10.1016/j.tig.2004.03.007.

    CAS  PubMed  Google Scholar 

  4. Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision. Trends Genet. 2004, 20: 80-86. 10.1016/j.tig.2003.12.003.

    CAS  PubMed  Google Scholar 

  5. Müller J, Reisz RR: Four well-constrained calibration points from the vertebrate fossil record for molecular clock estimates. Bioessays. 2005, 27: 1069-1075. 10.1002/bies.20286.

    PubMed  Google Scholar 

  6. Parham JF, Irmis RB: Caveats on the use of fossil calibrations for molecular dating: a comment on Near et al. Am Nat. 2008, 171: 132-136. 10.1086/524198.

    PubMed  Google Scholar 

  7. Bibi F, Vrba ES: Unraveling bovin phylogeny: accomplishments and challenges. BMC Biol. 2010, 8: 50-10.1186/1741-7007-8-50.

    PubMed Central  PubMed  Google Scholar 

  8. Hassanin A, Delsuc F, Ropiquet A, Hammer C, Jansen van Vuuren B, Matthee C, Ruiz-Garcia M, Catzeflis F, Areskoug V, Nguyen TT, Couloux A: Pattern and timing of diversification of Cetartiodactyla (Mammalia, Laurasiatheria), as revealed by a comprehensive analysis of mitochondrial genomes. C R Biol. 2012, 335: 32-50. 10.1016/j.crvi.2011.11.002.

    PubMed  Google Scholar 

  9. Spaulding M, O’Leary MA, Gatesy J: Relationships of Cetacea (Artiodactyla) among mammals: increased taxon sampling alters interpretations of key fossils and character evolution. PLoS ONE. 2009, 4: e7062-10.1371/journal.pone.0007062.

    PubMed Central  PubMed  Google Scholar 

  10. Hassanin A, Douzery EJP: Molecular and morphological phylogenies of Ruminantia and the alternative position of the Moschidae. Syst Biol. 2003, 52: 206-228. 10.1080/10635150390192726.

    PubMed  Google Scholar 

  11. Guha S, Goyal SP, Kashyap VK: Molecular phylogeny of musk deer: a genomic view with mitochondrial 16S rRNA and cytochrome b gene. Mol Phylogen Evol. 2007, 42: 585-597. 10.1016/j.ympev.2006.06.020.

    CAS  Google Scholar 

  12. Gibbard PL, Head MJ, Walker MJ: Formal ratification of the quaternary system/period and the pleistocene series/epoch with a base at 2.58 Ma. J Quat Sci. 2010, 25: 96-102. 10.1002/jqs.1338.

    Google Scholar 

  13. de Queiroz K: Toward an integrated system of clade names. Syst Biol. 2007, 56: 956-974. 10.1080/10635150701656378.

    PubMed  Google Scholar 

  14. O’Leary MA, Bloch JI, Flynn JJ, Gaudin TJ, Giallombardo A, Giannini NP, Goldberg SL, Kraatz BP, Luo Z-X, Meng J, et al: The placental mammal ancestor and the post-K-Pg radiation of placentals. Science. 2013, 339: 662-667. 10.1126/science.1229237.

    PubMed  Google Scholar 

  15. Nowak RM: Walker’s Mammals of the World, vol. 2. 1999, Baltimore, USA: Johns Hopkins University Press, 6

    Google Scholar 

  16. Antelopes, Deer, and Relatives. Edited by: Vrba ES, Schaller G. 2000, New Haven: Yale University Press

    Google Scholar 

  17. IUCN red list of threatened species. 2012, [http://www.iucnredlist.org. Accessed April 2013], .2

  18. Heller R, Frandsen P, Lorenzen E, Siegismund H: Are there really twice as many bovid species as we thought?. Syst Biol. 2013, 10.1093/sysbio/syt004.

    Google Scholar 

  19. Groves C, Grubb P: Ungulate Taxonomy. 2011, Baltimore: Johns Hopkins Univ

    Google Scholar 

  20. Janis CM, Scott KM: The interrelationships of higher ruminant families with special emphasis on the members of the Cervoidea. Am Mus Novit. 1987, 2893: 1-85.

    Google Scholar 

  21. Gentry AW, Rössner GE, Heizmann EPJ: Suborder Ruminantia. The Miocene Land Mammals of Europe. Edited by: Rössner GE, Heissig K. 1999, Munich: Verlag Friedrich Pfeil, 225-258.

    Google Scholar 

  22. Gradstein FM, Ogg JG, Smith AG: A Geologic Time Scale 2004. 2004, Cambridge: Cambridge University Press

    Google Scholar 

  23. Solounias N, Barry JC, Bernor RL, Lindsay EH, Raza SM: The oldest bovid from the Siwaliks, Pakistan. J Vertebr Paleontol. 1995, 15: 806-814. 10.1080/02724634.1995.10011263.

    Google Scholar 

  24. Agustí J, Cabrera L, Garces M, Krijgsman W, Oms O, Pares JM: A calibrated mammal scale for the Neogene of Western Europe. State of the art. Earth-Sci Rev. 2001, 52: 247-260. 10.1016/S0012-8252(00)00025-8.

    Google Scholar 

  25. Pilgrim GE: The fossil Bovidae of India. Palaeontol Indica. 1939, NS 26: 1-356.

    Google Scholar 

  26. Bibi F: Origin, paleoecology, and paleobiogeography of early Bovini. Palaeogeogr Palaeoclimatol Palaeoecol. 2007, 248: 60-72. 10.1016/j.palaeo.2006.11.009.

    Google Scholar 

  27. Bibi F: Ph.D. thesis. Evolution, Systematics, and Paleoecology of Bovinae (Mammalia: Artiodactyla) from the Late Miocene to the Recent. 2009, Yale University, Geology & Geophysics

    Google Scholar 

  28. Barry JC, Morgan ME, Flynn LJ, Pilbeam D, Behrensmeyer AK, Raza SM, Khan IA, Badgley C, Hicks J, Kelley J: Faunal and environmental change in the late Miocene Siwaliks of northern Pakistan. Paleobiology. 2002, 28 (S2): 1-71.

    Google Scholar 

  29. Badgley C, Barry JC, Morgan ME, Nelson SV, Behrensmeyer AK, Cerling TE, Pilbeam D: Ecological changes in miocene mammalian record show impact of prolonged climatic forcing. Proc Natl Acad Sci. 2008, 105: 12145-12149. 10.1073/pnas.0805592105.

    CAS  PubMed Central  PubMed  Google Scholar 

  30. Alcalá L, Morales J: A primitive caprine from the Upper Vallesian of La Roma 2 (Alfambra, Teruel, Aragón, Spain). Comptes Rendus de l’Académie des Sciences-Series IIA-Earth and Planetary Science. 1997, 324: 947-953.

    Google Scholar 

  31. Van Dam JA, Alcalá L, Zarza AA, Calvo JP, Garcés M, Krijgsman W: The upper miocene mammal record from the Teruel-Alfambra region (Spain). The MN system and continental stage/age concepts discussed. J Vertebr Paleontol. 2001, 21: 367-385. 10.1671/0272-4634(2001)021[0367:TUMMRF]2.0.CO;2.

    Google Scholar 

  32. Geraads D, Blondel C, Likius A, Mackaye HT, Vignaud P, Brunet M: New Hippotragini (Bovidae, Mammalia) from the late miocene of Toros-Menalla (Chad). J Vertebr Paleontol. 2008, 28: 231-242. 10.1671/0272-4634(2008)28[231:NHBMFT]2.0.CO;2.

    Google Scholar 

  33. Lebatard AE, Bourles DL, Duringer P, Jolivet M, Brauchert R, Carcaillet J, Schuster M, Arnaud N, Monie P, Lihoreau F, et al: Cosmogenic nuclide dating of Sahelanthropus tchadensis and Australopithecus bahrelghazali: mio-pliocene hominids from Chad. Proc Natl Acad Sci. 2008, 105: 3226-3231. 10.1073/pnas.0708015105.

    CAS  PubMed Central  PubMed  Google Scholar 

  34. Thomas H: Les bovidés du Miocène supérieur des couches de Mpesida et de la formation de Lukeino (district de Baringo, Kenya). Proceedings of the 8th Panafrican Congress of Prehistory and Quaternary Studies (Nairobi 1977). Edited by: Leakey REF, Ogot BA. 1980, Nairobi, 82-91.

    Google Scholar 

  35. Haile-Selassie Y, Vrba ES, Bibi F: Ardipithecus kadabba: Late Miocene Evidence from the Middle Awash, Ethiopia. Bovidae. Edited by: Haile-Selassie Y, WoldeGabriel G. 2009, Berkeley: University of California Press, 277-330.

    Google Scholar 

  36. WoldeGabriel G, Haile-Selassie Y, Renne PR, Hart WK, Ambrose SH, Asfaw B, Heiken G, White T: Geology and palaeontology of the late miocene middle Awash valley, Afar rift, Ethiopia. Nature. 2001, 412: 175-178. 10.1038/35084058.

    CAS  PubMed  Google Scholar 

  37. Deino AL, Tauxe L, Monaghan M, Hill A: Ar-40/Ar-30 geochronology and paleomagnetic stratigraphy of the Lukeino and lower Chemeron formations at Tabarin and Kapcheberek, Tugen Hills, Kenya. J Hum Evol. 2002, 42: 117-140. 10.1006/jhev.2001.0521.

    PubMed  Google Scholar 

  38. Renne PR, WoldeGabriel G, Hart WK, Heiken G, White TD: Chronostratigraphy of the miocene-pliocene Sagantole formation, middle Awash valley, Afar rift, Ethiopia. Geol Soc Am Bull. 1999, 111: 869-885. 10.1130/0016-7606(1999)111<0869:COTMPS>2.3.CO;2.

    CAS  Google Scholar 

  39. Gentry AW: Fossil Bovidae (Mammalia) from Langebaanweg, South Africa. Ann S Afr Mus. 1980, 79: 213-337.

    Google Scholar 

  40. Vrba ES: New fossils of Alcelaphini and Caprinae (Bovidae; Mammalia) from Awash, Ethiopia, and phylogenetic analysis of Alcelaphini. Palaeontol Afr. 1997, 34: 127-198.

    Google Scholar 

  41. Hendey QB: Geological succession at Langebaanweg, Cape province, and global events of the late tertiary. S Afr J Sci. 1981, 77: 33-38.

    Google Scholar 

  42. Vrba ES, Gatesy J: New antelope fossils from Awash, Ethiopia, and phylogenetic analysis of Hippotragini (Bovidae, Mammalia). Palaeontol Afr. 1994, 31: 55-72.

    Google Scholar 

  43. Gentry AW: Bovidae. Paleontology and Geology of Laetoli: Human Evolution in Context Volume 2. Edited by: Harrison T. 2011, New York: Springer, 363-465.

    Google Scholar 

  44. Deino AL: 40Ar/39Ar dating of Laetoli, Tanzania. Paleontology and Geology of Laetoli: Human Evolution in Context. Edited by: Harrison T. 2011, Dordrecht: Springer, 77-97.

    Google Scholar 

  45. Bibi F: Tragelaphus nakuae: Evolutionary change, biochronology, and turnover in the African plio-pleistocene. Zool J Linn Soc. 2011, 162: 699-711. 10.1111/j.1096-3642.2010.00691.x.

    Google Scholar 

  46. WoldeGabriel G, Endale T, White TD, Thouveny N, Hart WK, Renne PR, Asfaw B: The role of Tephra studies in African paleoanthropology as exemplified by the Sidi Hakoma Tuff. J Afr Earth Sci. 2013, 77: 41-58.

    Google Scholar 

  47. Reed K, Bibi F: Fossil Tragelaphini (Artiodactyla: Bovidae) from the Hadar formation, Afar Regional State, Ethiopia. J Mamm Evol. 2011, 1: 57-69.

    Google Scholar 

  48. Gentry AW: The Bovidae of the Omo Group deposits, Ethiopia (French and American collections). Les faunes Plio-Pléistocènes de la basse Vallée de l’Omo (Ethiopie); I: Perissodactyles-Artiodactyles (Bovidae). Volume 1985. Edited by: Coppens Y, Howell FC. 1985, Paris, France: CNRS, 119-191. Cahiers de Paléontologie

    Google Scholar 

  49. Feibel CS, Brown FH, Mcdougall I: Stratigraphic context of fossil hominids from the Omo group deposits - northern Turkana Basin, Kenya and Ethiopia. Am J Phys Anthropol. 1989, 78: 595-622. 10.1002/ajpa.1330780412.

    CAS  PubMed  Google Scholar 

  50. Gentry AW, Gentry A: Fossil Bovidae (Mammalia) of Olduvai Gorge, Tanzania; Part I. Bulletin of the British Museum Natural History Geology Series. 1978, 29: 289-446.

    Google Scholar 

  51. Gentry AW: Bovidae. Cenozoic Mammals of Africa. Edited by: Werdelin L, Sanders WJ. 2010, Berkeley: University of California Press, 747-803.

    Google Scholar 

  52. Hay RL: Geology of the Olduvai Gorge: A Study of Sedimentation in a Semiarid Basin. 1976, Berkeley: University of California Press

    Google Scholar 

  53. Deino AL: 40Ar/39Ar dating of Bed I, Olduvai Gorge, Tanzania, and the chronology of early pleistocene climate change. J Hum Evol. 2012, 63: 251-273. 10.1016/j.jhevol.2012.05.004.

    PubMed  Google Scholar 

  54. Clark JD, de Heinzelin J, Schick KD, Hart WK, White TD, WoldeGabriel G, Walter RC, Suwa G, Asfaw B, Vrba E: African Homo erectus: old radiometric ages and young Oldowan assemblages in the middle Awash valley, Ethiopia. Science. 1994, 264: 1907-1910. 10.1126/science.8009220.

    CAS  PubMed  Google Scholar 

  55. Pulquerio MJ, Nichols RA: Dates from the molecular clock: how wrong can we be?. Trends Ecol Evol. 2007, 22: 180-184. 10.1016/j.tree.2006.11.013.

    PubMed  Google Scholar 

  56. Drummond AJ, Ho SY, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.

    PubMed Central  PubMed  Google Scholar 

  57. Lanfear R, Welch JJ, Bromham L: Watching the clock: studying variation in rates of molecular evolution between species. Trends Ecol Evol. 2010, 25: 495-503. 10.1016/j.tree.2010.06.007.

    PubMed  Google Scholar 

  58. Ho SY: An examination of phylogenetic models of substitution rate variation among lineages. Biol Lett. 2009, 5: 421-424. 10.1098/rsbl.2008.0729.

    PubMed Central  PubMed  Google Scholar 

  59. Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simão TLL, Stadler T, et al: Impacts of the cretaceous terrestrial revolution and KPg extinction on mammal diversification. Science. 2011, 334: 521-524. 10.1126/science.1211028.

    CAS  PubMed  Google Scholar 

  60. Near TJ, Meylan PA, Shaffer HB: Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am Nat. 2005, 165: 137-146. 10.1086/427734.

    PubMed  Google Scholar 

  61. Sanders KL, Lee MS: Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol Lett. 2007, 3: 275-279. 10.1098/rsbl.2007.0063.

    CAS  PubMed Central  PubMed  Google Scholar 

  62. Marshall CR: A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am Nat. 2008, 171: 726-742. 10.1086/587523.

    PubMed  Google Scholar 

  63. Lukoschek V, Keogh JS, Avise JC: Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches. Syst Biol. 2012, 61: 22-43. 10.1093/sysbio/syr075.

    PubMed  Google Scholar 

  64. Rutschmann F, Eriksson T, Salim KA, Conti E: Assessing calibration uncertainty in molecular dating: the assignment of fossils to alternative calibration points. Syst Biol. 2007, 56: 591-608. 10.1080/10635150701491156.

    CAS  PubMed  Google Scholar 

  65. Ayala FJ: Molecular clock mirages. Bioessays. 1999, 21: 71-75. 10.1002/(SICI)1521-1878(199901)21:1<71::AID-BIES9>3.0.CO;2-B.

    CAS  PubMed  Google Scholar 

  66. Wilke T, Schultheiss R, Albrecht C: As time goes by: a simple fool’s guide to molecular clock approaches in invertebrates. Am Malacol Bull. 2009, 27: 25-45. 10.4003/006.027.0203.

    Google Scholar 

  67. Hugall AF, Lee MS: The likelihood node density effect and consequences for evolutionary studies of molecular rates. Evolution. 2007, 61: 2293-2307. 10.1111/j.1558-5646.2007.00188.x.

    PubMed  Google Scholar 

  68. Ho SY, Phillips MJ, Cooper A, Drummond AJ: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol. 2005, 22: 1561-1568. 10.1093/molbev/msi145.

    CAS  PubMed  Google Scholar 

  69. Phillips MJ: Branch-length estimation bias misleads molecular dating for a vertebrate mitochondrial phylogeny. Gene. 2009, 441: 132-140. 10.1016/j.gene.2008.08.017.

    CAS  PubMed  Google Scholar 

  70. Ho SY, Phillips MJ: Accounting for calibration uncertainty in phylogenetic estimation of evolutionary divergence times. Syst Biol. 2009, 58: 367-380. 10.1093/sysbio/syp035.

    PubMed  Google Scholar 

  71. Benton MJ, Ayala FJ: Dating the tree of life. Science. 2003, 300: 1698-1700. 10.1126/science.1077795.

    CAS  PubMed  Google Scholar 

  72. Linder HP, Hardy CR, Rutschmann F: Taxon sampling effects in molecular clock dating: an example from the African Restionaceae. Mol Phylogen Evol. 2005, 35: 569-582. 10.1016/j.ympev.2004.12.006.

    CAS  Google Scholar 

  73. Ronquist F, Klopfstein S, Vilhelmsen L, Schulmeister S, Murray DL, Rasnitsyn AP: A total-evidence approach to dating with fossils, applied to the early radiation of the Hymenoptera. Syst Biol. 2012, 61: 973-999. 10.1093/sysbio/sys058.

    PubMed Central  PubMed  Google Scholar 

  74. Métais G, Vislobokova I: Basal ruminants. The Evolution of Artiodactyls. Edited by: Prothero DR, Foss SE. 2007, Baltimore: Johns Hopkins University Press, 189-212.

    Google Scholar 

  75. Métais G, Chaimanee Y, Jaeger JJ, Ducrocq S: New remains of primitive ruminants from Thailand; evidence of the early evolution of the Ruminantia in Asia. Zool Scr. 2001, 30: 231-248. 10.1046/j.0300-3256.2001.00071.x.

    Google Scholar 

  76. Janis CM, Scott KM: The phylogeny of the Ruminantia (Artiodactyla, Mammalia). The Phylogeny and Classification of Tetrapods; Volume 2, Mammals. Volume 35b. Edited by: Benton MJ. 1988, London-New York: International: Academic Press [for the] Systematics Association, 273-282. Systematics Association Special Volume

    Google Scholar 

  77. O’Leary MA, Gatesy J: Impact of increased character sampling on the phylogeny of Cetartiodactyla (Mammalia): combined analysis including fossils. Cladistics. 2008, 24: 397-442. 10.1111/j.1096-0031.2007.00187.x.

    Google Scholar 

  78. Vislobokova I, Lavrov A: The earliest musk deer of the genus Moschus and their significance in clarifying of evolution and relationships of the family Moschidae. Paleontol J. 2009, 43: 326-338. 10.1134/S0031030109030125.

    Google Scholar 

  79. Prothero DR: Family Moschidae. The Evolution of Artiodactyls. Edited by: Prothero D, Foss S. 2007, Baltimore: Johns Hopkins University Press, 221-226.

    Google Scholar 

  80. Barry JC, Cote S, MacLatchy L, Lindsay EH, Kityo R, Rajpar AR: Oligocene and Early Miocene Ruminants (Mammalia, Artiodactyla) from Pakistan and Uganda. Palaeontol Electron. 2005, 8: 22A-29p.

    Google Scholar 

  81. Prothero DR: Foss SE (Eds.): The Evolution of Artiodactyls. 2007, Baltimore: Johns Hopkins University,

    Google Scholar 

  82. Gilbert C, Ropiquet A, Hassanin A: Mitochondrial and nuclear phylogenies of Cervidae (Mammalia, Ruminantia): systematics, morphology, and biogeography. Mol Phylogen Evol. 2006, 40: 101-117. 10.1016/j.ympev.2006.02.017.

    CAS  Google Scholar 

  83. Azanza B, Morales J: Tethytragus nov. gen et Gentrytragus nov. gen. Deux nouveaux Bovidés (Artiodactyla, Mammalia) du Miocène moyen. Proceedings, Koninklijke Nederlandse Akademie Van Wetenschappen. 1994, B 97: 249-282.

    Google Scholar 

  84. Made J: The bovid Pseudoeotragus seegrabensis nov. gen., nov. sp. from the Aragonia (Miocene) of Seegraben near Leoben (Austria). Proceedings Koninklijke Nederlandse Akademie Van Wetenschappen Ser B. 1989, 92: 215-240.

    Google Scholar 

  85. Thomas H: Les Bovidés anté-hipparions des Siwaliks inférieurs (Plateau du Potwar, Pakistan). Mémoires de la Société Géologique de France, Nouvelle Série. 1984, 145: 1-68.

    Google Scholar 

  86. Gentry AW: The Bovidae (Mammalia) of the Fort Ternan fossil fauna. Fossil Vertebrates of Africa, Vol 2. 1970, London: Academic Press, 243-323.

    Google Scholar 

  87. Köhler M: Boviden des turkischen Miozäns (Känozoikum und Braunkohlen der Türkei). Paleontologia y Evolució. 1987, 21: 133-246.

    Google Scholar 

  88. Gentry AW: The miocene differentiation of old world Pecora (Mammalia). Hist Biol. 1994, 7: 115-158. 10.1080/10292389409380449.

    Google Scholar 

  89. Shipman P, Walker A, Van Couvering JA, Hooker PJ, Miller JA: The Fort Ternan hominoid site, Kenya: geology, age, taphonomy and paleoecology. J Hum Evol. 1981, 10: 49-72. 10.1016/S0047-2484(81)80025-5.

    Google Scholar 

  90. Hassanin A, An J, Ropiquet A, Nguyen TT, Couloux A: Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: application to the tribe Bovini (Cetartiodactyla, Bovidae). Mol Phylogen Evol. 2013, 766-775.

    Google Scholar 

  91. Vrba ES, Schaller G: Phylogeny of Bovidae based on behavior, glands, skulls, and postcrania. Antelopes, Deer, and Relatives. Edited by: Vrba ES, Schaller G. 2000, New Haven: Yale University Press, 203-222.

    Google Scholar 

  92. Groves C, Schaller G: The phylogeny and biogeography of the newly discovered anamite artiodactyls. Antelopes, Deer, and Relatives. Edited by: Vrba ES, Schaller G. 2000, New Haven: Yale University Press, 261-282.

    Google Scholar 

  93. Vrba ES, Vaisnys JR, Gatesy JE, Desalle R, Wei KY: Analysis of pedomorphosis using allometric characters – the example of Reduncini antelopes (Bovidae, Mammalia). Syst Biol. 1994, 43: 92-116.

    Google Scholar 

  94. Omland KE: Correlated rates of molecular and morphological evolution. Evolution. 1997, 51: 1381-1393. 10.2307/2411190.

    Google Scholar 

  95. Thomas H: Anatomie crânienne et relations phylogénétiques du nouveau bovidé (Pseudoryx nghetinensis) découvert dans la cordillère annamitique au Vietnam. Mammalia. 1994, 58: 453-481.

    Google Scholar 

  96. Gatesy J, Arctander P: Hidden morphological support for the phylogenetic placement of Pseudoryx nghetinhensis with bovine bovids: a combined analysis of gross anatomical evidence and DNA sequences from five genes. Syst Biol. 2000, 49: 515-538. 10.1080/10635159950127376.

    CAS  PubMed  Google Scholar 

  97. Ropiquet A, Hassanin A: Molecular evidence for the polyphyly of the genus Hemitragus (Mammalia, Bovidae). Mol Phylogen Evol. 2005, 36: 154-168. 10.1016/j.ympev.2005.01.002.

    CAS  Google Scholar 

  98. Verkaar ELC, Nijman IJ, Beeke M, Hanekamp E, Lenstra JA: Maternal and paternal lineages in cross-breeding bovine species. Has wisent a hybrid origin?. Mol Biol Evol. 2004, 21: 1165-1170. 10.1093/molbev/msh064.

    CAS  PubMed  Google Scholar 

  99. Nijman IJ, Van Boxtel DCJ, Van Cann LM, Marnoch Y, Cuppen E, Lenstra JA: Phylogeny of Y chromosomes from bovine species. Cladistics. 2008, 24: 723-726. 10.1111/j.1096-0031.2008.00201.x.

    Google Scholar 

  100. Buntjer JB, Otsen M, Nijman IJ, Kuiper MTR, Lenstra JA: Phylogeny of bovine species based on AFLP fingerprinting. Heredity. 2002, 88: 46-51. 10.1038/sj.hdy.6800007.

    CAS  PubMed  Google Scholar 

  101. Pidancier N, Jordan S, Luikart G, Taberlet P: Evolutionary history of the genus Capra (Mammalia, Artiodactyla): discordance between mitochondrial DNA and Y-chromosome phylogenies. Mol Phylogen Evol. 2006, 40: 739-749. 10.1016/j.ympev.2006.04.002.

    CAS  Google Scholar 

  102. Bibi F, Fack F, Vrba ES: A new fossil caprin and a combined molecular and morphological bayesian analysis of Caprini (Mammalia: Bovidae). J Evol Biol. 2012, 25: 1843-1854. 10.1111/j.1420-9101.2012.02572.x.

    CAS  PubMed  Google Scholar 

  103. Moodley Y, Bruford MW, Bleidorn C, Wronski T, Apio A, Plath M: Analysis of mitochondrial DNA data reveals non-monophyly in the bushbuck (Tragelaphus scriptus) complex. Mamm Biol. 2009, 74: 418ñ422-

    Google Scholar 

  104. Vrba ES: The significance of bovid remains as indicators of environment and predation patterns. Fossils in the Making: Vertebrate Taphonomy and Paleoecology. Edited by: Behrensmeyer AK, Hill A. 1980, Chicago: Univ. of Chicago Press, 247-271.

    Google Scholar 

  105. Shipman P, Harris JM: Evolutionary history of the ìrobustî australopithecines New York: Aldine de Gruyter. Habitat preference and paleoecology of Australopithecus boisei in Eastern Africa. 1988, 343ñ381-

    Google Scholar 

  106. Bobe R, Eck GG: Responses of African bovids to pliocene climatic change. Paleobiology. 2001, 27: 1-47.

    Google Scholar 

  107. White TD, Ambrose SH, Suwa G, Su DF, DeGusta D, Bernor RL, Boisserie J-R, Brunet M, Delson E, Frost S, et al: Macrovertebrate paleontology and the pliocene habitat of Ardipithecus ramidus. Science. 2009, 326: 50-56. 10.1126/science.1179475.

    Google Scholar 

  108. Grubb P, Sandrock O, Kullmer O, Kaiser TM, Schrenk F: Relationships between eastern and southern African mammal faunas. African Biogeography, Climate Change and Early Hominid Evolution. Edited by: Bromage TG, Schrenk F. 2000, New York: Oxford University Press, 253-267.

    Google Scholar 

  109. Geraads D, Brunet M, Mackaye HT, Vignaud P: Pliocene Bovidae (Mammalia) from the Koro Toro Australopithecine sites, Chad. J Vertebr Paleontol. 2001, 21: 335-346. 10.1671/0272-4634(2001)021[0335:PBMFTK]2.0.CO;2.

    Google Scholar 

  110. de Heinzelin J, Clark JD, White T, Hart W, Renne P, WoldeGabriel G, Beyene Y, Vrba ES: Environment and behavior of 2.5-million-year-old Bouri hominids. Science. 1999, 284: 625-629. 10.1126/science.284.5414.625.

    CAS  PubMed  Google Scholar 

  111. Gilbert WH: Bovidae. Homo erectus: Pleistocene Evidence from the Middle Awash, Ethiopia. Edited by: Gilbert WH, Asfaw B. 2008, Berkeley: Univ. of California Press, 45-94.

    Google Scholar 

  112. Gentry AW: Bovidae of the Omo group deposits. Earliest man and environments in the Lake Rudolf Basin; stratigraphy, paleoecology, and evolution. Edited by: Coppens Y, Howell FC, Isaac GL, Leakey REF. 1976, Chicago, Ill., United States: Univ. Chicago Press, 275-292.

    Google Scholar 

  113. Bärmann EV, Rössner GE, Wörheide G: A revised phylogeny of Antilopini (Bovidae, Artiodactyla) using combined mitochondrial and nuclear genes. Mol Phylogen Evol. 2013, 484-493.

    Google Scholar 

  114. Gentry AW: Ruminantia (Artiodactyla). Geology and Paleontology of the Miocene Sinap Formation, Turkey. Edited by: Fortelius M, Kappelman J, Sen S, Bernor RL. 2003, New York: Columbia University Press

    Google Scholar 

  115. Bibi F: Mio-pliocene faunal exchanges and African biogeography: the record of fossil bovids. PLoS ONE. 2011, 6: e16688-10.1371/journal.pone.0016688.

    CAS  PubMed Central  PubMed  Google Scholar 

  116. Harris JM: Bovidae from the Lothagam succession. Lothagam: The Dawn of Humanity in Eastern Africa. Edited by: Leakey MG, Harris JM. 2003, New York: Columbia University Press, 531-579.

    Google Scholar 

  117. Vignaud P, Duringer P, Mackaye HT, Likius A, Blondel C, Boisserie JR, de Bonis L, Eisenmann V, Etienne ME, Geraads D, et al: Geology and palaeontology of the upper miocene Toros-Menalla hominid locality, Chad. Nature. 2002, 418: 152-155. 10.1038/nature00880.

    CAS  PubMed  Google Scholar 

  118. Thomas H: Les Bovidae (Artiodactyla; Mammalia) du Miocène du sous-continent indien, de la péninsule arabique et de l’Afrique: biostratigraphie, biogéographie et écologie. Palaeogeogr, Palaeoclimatol, Palaeoecol. 1984, 45: 251-299. 10.1016/0031-0182(84)90009-9.

    Google Scholar 

  119. Gentry AW: A new bovine (Bovidae, Artiodactyla) from the Hadar formation, Ethiopia. Trans R Soc S Afr. 2006, 61: 41-50. 10.1080/00359190609519952.

    Google Scholar 

  120. Bibi F, Hill A, Beech M, Yasin W: Late Miocene fossils from the Baynunah Formation, United Arab Emirates: summary of a decade of new work. Fossil Mammals of Asia: Neogene Biostratigraphy and Chronology. Edited by: Wang X, Flynn LJ, Fortelius M. 2013, New York: Columbia Univ. Press, 583-594.

    Google Scholar 

  121. Pyron RA: Divergence time estimation using fossils as terminal Taxa and the origins of Lissamphibia. Syst Biol. 2011, 60: 466-10.1093/sysbio/syr047.

    PubMed  Google Scholar 

  122. Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.

    CAS  PubMed Central  PubMed  Google Scholar 

  123. Kumar S, Skjaeveland A, Orr R, Enger P, Ruden T, Mevik B-H, Burki F, Botnen A, Shalchian-Tabrizi K: AIR: a batch-oriented web program package for construction of supermatrices ready for phylogenomic analyses. BMC Bioinformatics. 2009, 10: 357-10.1186/1471-2105-10-357.

    PubMed Central  PubMed  Google Scholar 

  124. Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012, 9: 772-772.

    CAS  PubMed Central  PubMed  Google Scholar 

  125. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.

    PubMed  Google Scholar 

  126. Yang Z, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006, 23: 212-226.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

I would like to thank T. Guiher, R. MacPhee, J. Müller and P. Velazco for discussions and comments that greatly helped this study and improved a previous draft of the manuscript. The editor and three anonymous reviewers also provided helpful comments and improvements. G. Rössner provided information on cervids. This study was initiated toward the end of a fellowship at the Museum für Naturkunde, funded by a Leibniz-DAAD fellowship. Support at the AMNH comes from a Gerstner Scholarship.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Faysal Bibi.

Additional information

Competing interests

The author declares that he has no competing interests.

Electronic supplementary material

Additional file 1: Supplementary text with details of fossil calibrations used.(PDF 247 KB)

12862_2013_2402_MOESM2_ESM.xls

Additional file 2: Table S1: Molecular divergence age estimates for ruminant clades from the analysis using all 16 fossil calibrations. (XLS 8 KB)

12862_2013_2402_MOESM3_ESM.xls

Additional file 3: Table S2: Percent difference between the molecular estimate and the median fossil age for the root and 16 calibrated nodes. These references are then summed at far right for each tree, and below for each node. Not surprisingly, the tree using all 16 calibrations also provides estimates with the lowest difference from the expected median age. The nodes with the highest summed difference scores are the Alcelaphus buselaphus crown and the Kobus ellipsiprymnus + K. leche crown nodes. (XLS 16 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Bibi, F. A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics. BMC Evol Biol 13, 166 (2013). https://doi.org/10.1186/1471-2148-13-166

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2148-13-166

Keywords