Molecular and morphological phylogenies of ruminantia and the alternative position of the moschidae

The ruminants constitute the largest group of ungulates, with >190 species, and its distribution is widespread throughout all continents except Australia and Antarctica. Six families are traditionally recognized within the suborder Ruminantia: Antilocapridae (pronghorns), Bovidae (cattle, sheep, and antelopes), Cervidae (deer), Giraffidae (giraffes and okapis), Moschidae (musk deer), and Tragulidae (chevrotains). The interrelationships of the families have been an area of controversy among morphology, palaeontology, and molecular studies, and almost all possible evolutionary scenarios have been proposed in the literature. We analyzed a large DNA data set (5,322 nucleotides) for 23 species including both mitochondrial (cytochrome b, 12S ribosomal RNA (rRNA), and 16S rRNA) and nuclear (kappa-casein, cytochrome P-450, lactoferrin, and alpha-lactalbumin) markers. Our results show that the family Tragulidae occupies a basal position with respect to all other ruminant families, confirming the traditional view that separates Tragulina and Pecora. Within the pecorans, Antilocapridae and Giraffidae emerge first, and the families Bovidae, Moschidae, and Cervidae are allied, with the unexpected placement of Moschus close to bovids rather than to cervids. We used these molecular results to assess the homoplastic evolution of morphological characters within the Ruminantia. A Bayesian relaxed molecular clock approach based on the continuous autocorrelation of evolutionary rates along branches was applied to estimate the divergence ages between the major clades of ruminants. The evolutionary radiation of Pecora occurred at the Early/Late Oligocene transition, and Pecoran families diversified and dispersed rapidly during the Early and Middle Miocene. We propose a biogeographic scenario to explain the extraordinary expansion of this group during the Cenozoic era.     

Phylogenetic relationships among deer in China derived from mitochondrial DNA cytochromeb sequences

The phylogenetic relationships of Cervidae and Moschidae were examined using partial sequence data of mitochondrial DNA (mtDNA) cytochromeb. Ten new sequences were obtained for six species of Cervidae and Moschidae, and aligned with those previously reported for other deer species. Our results demonstrated that the phylogenetic status of the taxa inferred from molecular data was congruent with taxonomy based on morphological studies. Cervidae formed a monophyletic group that consists of four subfamilies: Cervinae, Muntiacinae, Hydropotinae, and Odocoileinae. Moschidae occurred at the base of the Cervidae clade. On the basis of molecular clocks for genetic distance, the divergence time of mtDNA haplotypes within the subfamily Cervinae, among subfamilies in Cervidae, and between Moschidae and Cervidae was estimated to date 2–7 MYA, 6–10 MYA and 8–13 MYA, respectively.     

Molecular phylogeny of deer (Cervidae: Artiodactyla)

Sequences of mitochondrial genes 12S and 16S rRNA (2 445 bp) and the region of the nuclear beta-spectrin gene (828 bp) were analyzed in members of the family Cervidae and in other artiodactyls. Several molecular synapomorphies characteristic both of Cervidae and musk deer have been found. According to our data, Cervidae is a sister clade to Bovidae, which are very close to Moschidae. The family Giraffidae is exterior to this common clade, while Antilocapridae occupies a more basal position. The family Cervidae proper splits into three clades including the genera Cervus and Muntiacus (1), Capreolus, Hydropotes, Alces (2), and Rangifer, Odocoileus, and the remaining genera (3). In general, our phylogenetic reconstructions conform to the results of earlier molecular genetic studies, but substantially differ from the traditional taxonomy of Ruminantia.     

The tribal radiation of the family Bovidae (Artiodactyla) and the evolution of the mitochondrial cytochrome b gene

The nucleotide sequence of the complete mitochondrial cytochrome b gene has been determined and compared for 51 species of the family Bovidae and 10 potential pecoran and tragulid outgroups. A detailed saturation analysis at each codon position relative to the maximum parsimony procedure indicates that all transitions on third codon positions do not accumulate in a similar fashion: C-T are more saturated than A-G substitutions. The same trend is observed for second positions but not for first positions where A-G and C-T transitions exhibit roughly the same levels of saturation. Maximum parsimony reconstructions were weighted according to these observations. Maximum parsimony, maximum likelihood, and distance phylogenetic reconstructions all depict a major split within Bovidae. The subfamily Bovinae includes four multifurcating tribes and subtribes: Boselaphini, Tragelaphini, cattle-Bovini (Bos and Bison), and buffalo-Bovini (Bubalus and Syncerus). Its sister group is the subfamily Antilopinae, i.e., all non-Bovinae taxa, represented by seven lineages: Antilopini (including Saiga), Caprini sensu lato (i. e., Caprinae including Pantholops), Hippotragini, Alcelaphini, Reduncini (including Pelea), Aepyceros possibly linked to Neotragus, and Cephalophini possibly linked to Oreotragus (the suni and the klipspringer being members of a polyphyletic Neotragini). These various tribes and major lineages were produced by two noteworthy explosive radiations, which occurred simultaneously between 12.0 and 15.3 MY (Middle Miocene) in the subfamilies Bovinae and Antilopinae.     

Molecular Phylogeny of Deer (Cervidae: Artiodactyla)

Sequences of mitochondrial genes 12S and 16S rRNA (2445 bp) and the region of the nuclear betaspectrin gene (828 bp) were analyzed in members of the family Cervidae and in some other artiodactyls. Several molecular synapomorphies characteristic both of Cervidae and musk deer have been found. According to our data, Cervidae is a sister clade to Bovidae, which are very close to Moschidae. The family Giraffidae is exterior to this common clade, while Antilocapridae occupies a more basal position. The family Cervidae proper splits into three clades including the genera Cervus and Muntiacus (1), Capreolus, Hydropotes, Alces (2), and Rangifer, Odocoileus, and the remaining genera (3). In general, our phylogenetic reconstructions conform to the results of earlier molecular genetic studies, but substantially differ from the traditional taxonomy of Ruminantia.__________Translated from Genetika, Vol. 41, No. 7, 2005, pp. 910–918.Original Russian Text Copyright © 2005 by Kuznetsova, Kholodova, Danilkin.     

The melanocyte-stimulating hormone receptor (MC1-R) gene as a tool in evolutionary studies of artiodactyles

The complete coding region of the melanocyte-stimulating hormone receptor (MC1-R) gene was characterized in species belonging to the two families Bovidae and Cervidae; cattle (Bos taurus), sheep (Ovis aries), goat (Capra hircus), muskox (Ovibos moschatus), roe deer (Capreolus capreolus), reindeer (Rangifer tarandus), moose (Alces alces), red deer (Cervus elaphus) and fallow deer (Dama dama). This well conserved gene is a central regulator of mammalian coat colour. Examination of the interspecies variability revealed a 5.3-6.8% divergence between the Cervidae and Bovidae families, whereas the divergence within the families were 1.0-3.1% and 1.2-4.6%, respectively. Complete identity was found when two subspecies of reindeer, Eurasian tundra reindeer (R.t. tarandus) and Svalbard reindeer (R.t. platvrhynehus), were analyzed. An rooted phylogenetic tree based on Bovidae and Cervidae MC1-R DNA sequences was in complete agreement with current taxonomy, and was supported by bootstrapping analysis. Due to different frequencies of silent vs. replacement mutations, the amino acid based phylogenetic tree contains several dissimilarities when compared to the DNA based phylogenetic tree.     

A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla)

Despite the biological and economic importance of the Cetartiodactyla, the phylogeny of this clade remains controversial. Using the supertree approach of matrix representation with parsimony, we present the first phylogeny to include all 290 extant species of the Cetacea (whales and dolphins) and Artiodactyla (even-toed hoofed mammals). At the family-level, the supertree is fully resolved. For example, the relationships among the Ruminantia appear as (((Cervidae, Moschidae) Bovidae) (Giraffidae, Antilocapridae) Tragulidae). However, due to either lack of phylogenetic study or contradictory information, polytomies occur within the clades Sus, Muntiacus, Cervus, Delphinidae, Ziphiidae and Bovidae. Complete species-level phylogenies are necessary for both illustrating and analysing biological, geographical and ecological patterns in an evolutionary framework. The present species-level tree of the Cetartiodactyla provides the first opportunity to examine comparative hypotheses across entirely aquatic and terrestrial species within a single mammalian order.     

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.     

Phylogeny and evolution of Cervidae based on complete mitochondrial genomes

Mitochondrial DNA sequences can be used to estimate phylogenetic relationships among animal taxa and for molecular phylogenetic evolution analysis. With the development of sequencing technology, more and more mitochondrial sequences have been made available in public databases, including whole mitochondrial DNA sequences. These data have been used for phylogenetic analysis of animal species, and for studies of evolutionary processes. We made phylogenetic analyses of 19 species of Cervidae, with Bos taurus as the outgroup. We used neighbor joining, maximum likelihood, maximum parsimony, and Bayesian inference methods on whole mitochondrial genome sequences. The consensus phylogenetic trees supported monophyly of the family Cervidae; it was divided into two subfamilies, Plesiometacarpalia and Telemetacarpalia, and four tribes, Cervinae, Muntiacinae, Hydropotinae, and Odocoileinae. The divergence times in these families were estimated by phylogenetic analysis using the Bayesian method with a relaxed molecular clock method; the results were consistent with those of previous studies. We concluded that the evolutionary structure of the family Cervidae can be reconstructed by phylogenetic analysis based on whole mitochondrial genomes; this method could be used broadly in phylogenetic evolutionary analysis of animal taxa.     

Phylogenetic position of the saola (Pseudoryx nghetinhensis) inferred from cytogenetic analysis of eleven species of Bovidae

Previous morphological and molecular analyses failed to resolve the phylogenetic position of the critically endangered saola (Pseudoryx nghetinhensis) with respect to its placement in Bovina (cattle, bison, and yak) or Bubalina (Asian and African buffaloes). In the present study, G- and C-banding, Ag-staining and FISH with 28S and telomeric probes was undertaken for 17 bovid species. An analysis of these data allowed us to identify 49 structural rearrangements that included autosomes, gonosomes and 17 different NOR sites. The combined data set was subjected to a cladistic analysis aimed at: (i) providing new insights on phylogenetic relationships of the saola and other species within the subfamily Bovinae, and (ii) testing the suitability of different classes of chromosomal characters for phylogenetic reconstruction of the family Bovidae. The study revealed that nucleolar organizing regions (NORs) are phylogenetically informative. It was shown that at least one, or sometimes two of these characters punctuate divergences that include nodes that are the most basal in the tree, to those that are the most recent. In this context, the shared presence of three NORs in saola and species of Syncerus and Bubalus strongly suggests the saola's placement within the subtribe Bubalina. This contrasts with Robertsonian rearrangements which are informative only at the generic level. These findings suggest that NORs are an important and frequently overlooked source of additional phylogenetic information within the Bovidae that may also have applicability at higher taxonomic levels, possibly even for Pecora.     

Mitochondrial 12S rDNA sequence relationships suggest that the enigmatic bovid "Linh Duong" Pseudonovibos spiralis is closely related to buffalo

The phylogenetic position and taxonomic status of the recently described Southeast Asian endemic bovid Pseudonovibos spiralis were deduced from nearly complete 12S mitochondrial rDNA sequences of this species and Bubalus bubalis alongside 26 sequences of Bovidae from GenBank using Cervus elaphus (Cervidae) as outgroup. Maximum-likelihood analyses performed by PUZZLE and fastDNAml nested P. spiralis at the base of the subtribe buffalo Bovini, suggesting the close relationship of this enigmatic species with buffalo and enabling its distinction into the separate     

New insights into the karyotypic relationships of Chinese muntjac (Muntiacus reevesi), forest musk deer (Moschus berezovskii) and gayal (Bos frontalis)

To investigate the karyotypic relationships between Chinese muntjac (Muntiacus reevesi), forest musk deer (Moschus berezovskii) and gayal (Bos frontalis), a complete set of Chinese muntjac chromosome-specific painting probes has been assigned to G-banded chromosomes of these three species. Sixteen autosomal probes (i.e. 6-10, 12-22) of the Chinese muntjac each delineated one pair of conserved segments in the forest musk deer and gayal, respectively. The remaining six autosomal probes (1-5, and 11) each delineated two to five pairs of conserved segments. In total, the 22 autosomal painting probes of Chinese muntjac delineated 33 and 34 conserved chromosomal segments in the genomes of forest musk deer and gayal, respectively. The combined analysis of comparative chromosome painting and G-band comparison reveals that most interspecific homologous segments show a high degree of conservation in G-banding patterns. Eleven chromosome fissions and five chromosome fusions differentiate the karyotypes of Chinese muntjac and forest musk deer; twelve chromosome fissions and six fusions are required to convert the Chinese muntjac karyotype to that of gayal; one chromosome fission and one fusion separate the forest musk deer and gayal. The musk deer has retained a highly conserved karyotype that closely resembles the proposed ancestral pecoran karyotype but shares none of the rearrangements characteristic for the Cervidae and Bovidae. Our results substantiate that chromosomes 1-5 and 11 of Chinese muntjac originated through exclusive centromere-to-telomere fusions of ancestral acrocentric chromosomes.     

Evolutionary affinities of the enigmatic saola (Pseudoryx nghetinhensis) in the context of the molecular phylogeny of Bovidae.

To elucidate the systematic status of the enigmatic saola (Pseudoryx nghetinhensis), a new bovid genus recently discovered in Vietnam, and to investigate phylogenetic relationships within the family Bovidae, four distinct DNA markers were sequenced. Complete mitochondrial cytochrome b (1143 bp) and 12S rRNA (956 bp) genes and non-coding regions from the nuclear genes for aromatase cytochrome P-450 (199 bp) and lactoferrin (338 bp) have been compared for 25 bovid species and three Cervidae and Antilocapridae outgroups. Independent and/or combined analyses of the four nucleotide matrices through maximum parsimony and maximum-likelihood methods indicated that Bovidae consists of two major lineages, i.e. Bovinac which contains the tribes Bovini, Boselaphini and Tragelaphini, and Antilopinae which encompasses all other bovids. Within Bovinae, the tribe Bovini is divided into buffalo Bovini (Bubalus and Syncerus) and cattle Bovini (Bos and Bison) and Tragelaphini are possibly related to Boselaphini. Pseudoryx is shown to be (i) robustly nested within Bovinae; (ii) strongly associated with Bovini; and (iii) tentatively sharing a sister-group relationship with cattle Bovini. Within Antilopinae, three robust clades are in evidence: (i) Hippotragus and Damaliscus are linked to Ovis; (ii) Aepyceros joins Neotragus; and (iii) Cephalophus clusters with Oreotragus.     

A comparative study of mammalian diversification pattern

Although mammals have long been regarded as a successful radiation, the diversification pattern among the clades is still poorly known. Higher-level phylogenies are conflicting and comprehensive comparative analyses are still lacking. Using a recently published supermatrix encompassing nearly all extant mammalian families and a novel comparative likelihood approach (MEDUSA), the diversification pattern of mammalian groups was examined. Both order- and family-level phylogenetic analyses revealed the rapid radiation of Boreoeutheria and Euaustralidelphia in the early mammalian history. The observation of a diversification burst within Boreoeutheria at approximately 100 My supports the Long Fuse model in elucidating placental diversification progress, and the rapid radiation of Euaustralidelphia suggests an important role of biogeographic dispersal events in triggering early Australian marsupial rapid radiation. Diversification analyses based on family-level diversity tree revealed seven additional clades with exceptional diversification rate shifts, six of which represent accelerations in net diversification rate as compared to the background pattern. The shifts gave origin to the clades Muridae+Cricetidae, Bovidae+Moschidae+Cervidae, Simiiformes, Echimyidae, Odontoceti (excluding Physeteridae+Kogiidae+Platanistidae), Macropodidae, and Vespertilionidae. Moderate to high extinction rates from background and boreoeutherian diversification patterns indicate the important role of turnovers in shaping the heterogeneous taxonomic richness observed among extant mammalian groups. Furthermore, the present results emphasize the key role of extinction on erasing unusual diversification signals, and suggest that further studies are needed to clarify the historical radiation of some mammalian groups for which MEDUSA did not detect exceptional diversification rates.     

The Levator Veli Palatini Muscle in Artiodactyls—A Comparative Ontogenetic Study

It is accepted in the literature that the levator veli palatini muscle of artiodactyls originates at the ectotympanic bone, a statement based on macroscopic dissection of adult specimens. The study of 34 histological serial sections of fetal heads of 23 species of artiodactyls revealed that this generalization must be modified: in the studied Camelidae, Suidae, Hippopotamidae, Giraffidae, and in some Bovidae (namely Tragelaphus and Antidorcas) the primary attachment of this functionally important muscle is at the tendinous intersection with the tensor veli palatini. Primary ontogenetic attachments are considered as relevant for defining homologies. By outgroup comparison (Felis and Diceros), this structural connection (character state 1) is also found in the Scrotifera—and hence may be considered as plesiomorphic for the Artiodactyla and its subunits. Only in the Tragulidae, Cervidae, Moschidae, and some Bovidae is a secondary attachment at the ectotympanic observed, which is interpreted as apomorphic for these taxa; possibly this character state 2 developed homoplastically several times. Bovidae show a mixed distribution of this character: Tragelaphus, Aepyceros, and Antidorcas show only a connection of the levator with the tensor veli; in Neotragus, Raphicerus, and Sylvicapra there exists an additional insertion at the ectotympanic; only Bos, Cephalophus, Damaliscus, and Ovis have a primary origin at the ectotympanic. It can be demonstrated in late fetal Sus that a secondary insertion of the levator veli at the ectotympanic is established during ontogeny; in a late fetal Ovis a secondary contact with the tensor veli is realized. The interpretation of this character distribution depends not only on an intrinsic polarity (‘Lesrichtung’), but also on the assumed character state of the groundplan of the common ancestor of the Bovidae. The anatomical observations are documented by photographs of relevant histological sections. The character states are mapped on a simplified and synoptic cladogram of extant artiodactyls; their pattern of evolutionary transformation as well as their relevance for the phylogenetic systematics of this mammalian order are discussed.     

A complete estimate of the phylogenetic relationships in Ruminantia: a dated species-level supertree of the extant ruminants

This paper presents the first complete estimate of the phylogenetic relationships among all 197 species of extant and recently extinct ruminants combining morphological, ethological and molecular information. The composite tree is derived by applying matrix representation using parsimony analysis to 164 previous partial estimates, and is remarkably well resolved, containing 159 nodes (> 80 % of the potential nodes in the completely resolved phylogeny). Bremer decay index has been used to indicate the degree of certainty associated with each clade. The ages of over 80% of the clades in the tree have been estimated from information in the literature. The supertree for Ruminantia illustrates which areas of ruminant phylogeny are still only roughly known because of taxa with controversial relationships (e.g. Odocoileini, Antilopinae) or not studied in great detail (e.g. Muntiacus). It supports the monophyly of the ruminant families and Pecora. According to this analysis Antilocapridae and Giraffidae constitute the superfamily Giraffoidea, which is the sister group of a clade clustering Bovoidea and Cervoidea. The position of several taxa whose systematic positions have remained controversial in the past (Saiga, Pelea, Aepycerus, Pantholops, Ammotragus, Pseudois) is unambiguously established. Nevertheless, the position of Neotragus and Oreotragus within the original radiation of the non-bovine bovids remains unresolved in the present analysis. It also shows that six successive rapid cladogenesis events occurred within the infraorder Pecora during the Oligocene to middle Pliocene, which coincided with periods of global climatic change. Finally, the presented supertree will be a useful framework for comparative and evolutionary biologists interested in studies involving the ruminants.     

Evolutionary histories of highly repeated DNA families among the artiodactyla (mammalia)

Six highly repeated DNA families were analyzed using Southern blotting and fluorescence in situ hybridization in a comparative study of 46 species of artiodactyls belonging to seven of the eight extant taxonomic families. Two of the repeats, the dispersed bovine-Pst family and the localized 1.715 component, were found to have the broadest taxonomic distributions, being present in all pecoran ruminants (Giraffidae, Cervidae, Antilocapridae, and Bovidae), indicating that these repeats may be 25–40 million years old. Different 1.715 restriction patterns were observed in different taxonomic families, indicating that independent concerted evolution events have homogenized different motifs in different lineages. The other four satellite arrays were restricted to the Bovini and sometimes to the related Boselaphini and Tragelaphini. Results reveal that among the two compound satellites studied, the two components of the 1.711a originated simultaneously, whereas the two components of the 1.711b originated at two different historical times, perhaps as many as 15 million years apart. Systematic conclusions support the monophyly of the infraorder Pecora, the monophyly of the subfamily Bovinae (containing the Boselaphini, Bovini, and Tragelaphini), an inability to resolve any interrelationships among the other tribes of bovids, paraphyly of the genus Bos with respect to Bison, and a lack of molecular variation among two morphologically and ecologically distinct subspecies of African buffaloes (Syncerus caffer cafer and S. c. nanus). Cytogenetically, a reduction in diploid chromosome numbers through centric fusion in derived karyotypes is accompanied by a loss of centromeric satellite DNA. The nilgai karyotype contains an apparent dicentric chromosome as evidenced by the sites of 1.715 hybridization. Telomeric sequences have been translocated to the centromeres without concomitant chromosomal rearrangement in Thompson's gazelle. Received: 18 June 1995 / Accepted: 1 September 1995     

ALLOZYME DIVERGENCE WITHIN THE BOVIDAE

We describe a phylogeny of the Bovidae based on 40 allozyme loci in 27 species, representing 10 of the 14 bovid tribes described by Vrba (1985). Giraffe represented a related family (Giraffidae). A phenogram was derived using the unweighted pair-group method with arithmetic means (UPGMA), based on Nei's genetic distances (ND) between species. A tree was also derived using the neighbor-joining technique, also based on ND. To provide a cladistic interpretation, the data were analyzed by a maximum parsimony method (phylogenetic analysis using parsimony, PAUP). We found marked divergence within the Bovidae, consistent with the appearance of the family in the early Miocene. Unexpectedly, the most divergent species was the impala, which occupied a basal position in all trees. Species in the tribe Alcelaphini were the most derived taxa in all trees. These patterns conflict strongly with the previous taxonomic alliance, based on immuno-distance and anatomical evidence, of the impala as a sister group of the Alcelaphini. All trees agreed that tribes described by Vrba (1985) are monophyletic, except the Neotragini, which was polyphyletic, with suni occupying a long branch by itself. The dikdik and klipspringer were consistently placed as sister taxa to species in the Antilopini. Three tribes (Aepycerotini, Tragelaphini and Cephalophini), whose fossils have not been found outside Africa, were basal in all trees, suggesting that bovids originated in Africa. Nodes connecting the remaining tribes were closely clustered, a pattern that agrees with fossil evidence of rapid divergence within the Bovidae in the mid-Miocene (about 15 mybp). The allozyme data suggested a second phase of rapid divergence within tribes during the Plio-Pleistocene, a pattern that also agrees with fossil evidence. Rates of bovid divergence have therefore been far from constant. However, the clustering of nodes imparts considerable uncertainty to the branching order leading to the derived tribes, and to a lesser extent, species within tribes. The classical division of the Bovidae into the Boodontia and Aegeodontia does not agree with the phylogenetic grouping of tribes presented in this analysis. However, the maximum parsimony tree derived using ‘local’ branch swapping clustered all grazing species into a derived, monophyletic group, suggesting that grazing may have evolved only once in bovid evolution.     

Sexually and Developmentally Dimorphic Grooming: A Comparative Survey of the Ungulata

Grooming is effective in removing fitness‐compromising ectoparasites, such as ticks. Tick‐removal grooming is regulated both by a central control mechanism (programmed grooming model) and by cutaneous stimulation from tick bites (stimulus‐driven model). The programmed grooming model predicts that organismic factors that impact the cost‐benefit ratio of grooming (e.g., sex and body size) will influence the rate of grooming. The ‘vigilance principle’ predicts that breeding males of sexually dimorphic species will groom less than females to maintain high levels of vigilance for rival males and oestrus females. The intraspecific body size principle predicts that juveniles will engage in more frequent grooming than larger‐bodied adults to compensate for the higher cost of tick infestation for smaller animals. To examine the generalizability of these predictions, we surveyed the grooming rate of 53 species/subspecies of ungulates (primarily Bovidae and Cervidae) in tick‐free zoological parks in which stimulus‐driven grooming was controlled for. Matched‐pair analysis supported both predictions. Males delivered fewer oral grooming episodes per hour than females in 85% of sexually dimorphic species (n = 40), but only 56% of monomorphic species (n = 11) exhibited sexually dimorphic grooming, which was not different from random. Juvenile oral episodes per hour exceeded that of adult females in 74% of surveyed species (n = 36). As predicted by the body size principle, the grooming rate of juveniles declined as juveniles grew larger and less vulnerable to tick infestation. In separate analyses of Bovidae and Cervidae to control for common ancestry, both families supported sexually dimorphic grooming, but developmentally dimorphic grooming was supported only for Bovidae. Our results indicate that sexually dimorphic grooming is widespread in the ungulate species surveyed, suggesting that programmed grooming evolved at least as early as the common ancestor to the Artiodactyla.