Phylogenetic relationships of artiodactyls and cetaceans as deduced from the comparison of cytochrome b and 12S rRNA mitochondrial sequences

A data set of complete mitochondrial cytochrome b and 12S rDNA sequences is presented here for 17 representatives of Artiodactyla and Cetacea, together with potential outgroups (two Perissodactyla, two Carnivora, two Tethytheria, four Rodentia, and two Marsupialia). We include seven sequences not previously published from Hippopotamidae (Ancodonta) and Camelidae (Tylopoda), yielding a total of nearly 2.1 kb for both genes combined. Distance and parsimony analyses of each gene indicate that 11 clades are well supported, including the artiodactyl taxa Pecora, Ruminantia (with low 12S rRNA support), Tylopoda, Suina, and Ancodonta, as well as Cetacea, Perissodactyla, Carnivora, Tethytheria, Muridae, and Caviomorpha. Neither the cytochrome b nor the 12S rDNA genes resolve the relationships between these major clades. The combined analysis of the two genes suggests a monophyletic Cetacea +Artiodactyla clade (defined as "Cetartiodactyla"), whereas Perissodactyla, Carnivora, and Tethytheria fall outside this clade. Perissodactyla could represent the sister taxon of Cetartiodactyla, as deduced from resampling studies among outgroup lineages. Cetartiodactyla includes five major lineages: Ruminantia, Tylopoda, Suina, Ancodonta, and Cetacea, among which the phylogenetic relationships are not resolved. Thus, Suiformes do not appear to be monophyletic, justifying their split into the Suina and Ancodonta infraorders. An association between Cetacea and Hippopotamidae is supported by the cytochrome b gene but not by the 12S rRNA gene. Calculation of divergence dates suggests that the Cetartiodactyla could have diverged from other Ferungulata about 60 MYA.      

The complete mitochondrial DNA sequence of the greater Indian rhinoceros, Rhinoceros unicornis, and the Phylogenetic relationship among Carnivora, Perissodactyla, and Artiodactyla (+ Cetacea)

The sequence (16,829 nt) of the complete mitochondrial genome of the greater Indian rhinoceros, Rhinoceros unicornis, was determined. Like other perissodactyls studied (horse and donkey) the rhinoceros demonstrates length variation (heteroplasmy) associated with different numbers of repetitive motifs in the control region. The 16,829-nt variety of the molecule includes 36 identical control region motifs. The evolution of individual peptide-coding genes was examined by comparison with a distantly related perissodactyl, the horse, and the relationships among the orders Carnivora, Perissodactyla, and Artiodactyla (+ Cetacea) were examined on the basis of concatenated sequences of 12 mitochondrial peptide-coding genes. The phylogenetic analyses grouped Carnivora, Perissodactyla, and Artiodactyla (+ Cetacea) into a superordinal clade and within this clade a sister group relationship was recognized between Carnivora and Perissodactyla to the exclusion of Artiodactyla (+ Cetacea) . On the basis of the molecular difference between the rhinoceros and the horse and by applying as a reference to Artiodactyl/Cetacean divergence set at 60 million years ago (MYA), the evolutionary divergence between the families Rhinocerotidae and Equidae was dated to approximately 50 MYA.      

Characterization and phylogenetic utility of the mammalian protamine p1 gene

We sequenced the protamine P1 gene (ca. 450 bp) from 20 bats (order Chiroptera) and the flying lemur (order Dermoptera). We compared these sequences with published sequences from 19 other mammals representing seven orders (Artiodactyla, Carnivora, Cetacea, Perissodactyla, Primates, Proboscidea, and Rodentia) to assess structure, base compositional bias, and phylogenetic utility. Approximately 80% of second codon positions were guanine, resulting in protamine proteins containing a high frequency of arginine residues. Our data indicate that codon usage for arginine differs among higher mammalian taxa. Parsimony analysis of 40 species representing nine orders produced a well-resolved tree in which most nodes were supported strongly, except at the lowest taxonomic levels (e.g., within Artiodactyla and Vespertilionidae). These data support monophyly of several taxa proposed by morphologic and molecular studies (all nine orders: Laurasiatheria, Cetartiodactytla, Yangochiroptera, Noctilionoidea, Rhinolophoidea, Vespertilionoidea, Phyllostomidae, Natalidae, and Vespertilionidae) and, in agreement with recent molecular studies, reject monophyly of Archonta, Volitantia, and Microchiroptera. Bats were sister to a clade containing Perissodactyla, Carnivora, and Cetartiodactyla, and, although not unequivocally, rhinolophoid bats (traditional microchiropterans) were sister to megachiropterans. Sequences of the protamine P1 gene are useful for resolving relationships at and above the familial level in bats, and generally within and among mammalian orders, but with some drawbacks. The coding and intervening sequences are small, producing few phylogenetically informative characters, and aligning the intron is difficult, even among closely related families. Given these caveats, the protamine P1 gene may be important to future systematic studies because its functional and evolutionary constraints differ from other genes currently used in systematic studies.     

A molecular perspective on mammalian evolution from the gene encoding interphotoreceptor retinoid binding protein, with convincing evidence for bat monophyly.

The evolutionary relationships of the various orders of placental mammals remain an issue of uncertainty and controversy. Molecular studies of mammalian phylogeny at the DNA level that include more than just a few orders are still relatively meager. Here we report results on mammalian phylogeny deduced from the coding sequence of the single-copy nuclear gene for the interphotoreceptor retinoid binding protein (IRBP). Analysis of 13 species representing eight eutherian orders and one marsupial yielded results that falsify the hypothesis that megachiropteran bats are "flying primates," only convergently resembling microchiropteran bats. Instead, in agreement with more traditional views, as well as those from other recent molecular studies, the results strongly support a monophyletic Chiroptera (micro- and megabats grouped together). The IRBP results also offer some rare molecular support for the Glires concept, in which rodents and lagomorphs form a superordinal grouping. Also in congruence with other recent molecular evidence, IRBP sequences do not support the view of a superorder Archonta that includes Chiroptera along with Dermoptera (flying lemur), Scandentia (tree shrew), and Primates. IRBP was not however, without its shortcomings as a molecular phylogenetic system: high levels of homoplasy, evident in the marsupial outgroup, did not allow us to properly root the tree, and several of the higher level eutherian clades were only weakly supported (e.g., a Carnivora/Chiroptera clade and an Artiodactyla/Carnivora/Chiroptera clade). We suggest that these shortcomings may be diminished as the phylogenetic density of the data set is increased.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

The phylogenetic position of the Talpidae within eutheria based on analysis of complete mitochondrial sequences

The complete mitochondrial (mt) genome of the mole Talpa europaea was sequenced and included in phylogenetic analyses together with another lipotyphlan (insectivore) species, the hedgehog Erinaceus europaeus, and 22 other eutherian species plus three outgroup taxa (two marsupials and a monotreme). The phylogenetic analyses reconstructed a sister group relationship between the mole and fruit bat Artibeus jamaicensis (order Chiroptera). The Talpa/Artibeus clade constitutes a sister clade of the cetferungulates, a clade including Cetacea, Artiodactyla, Perissodactyla, and Carnivora. A monophyletic relationship between the hedgehog and the mole was significantly rejected by maximum parsimony and maximum likelihood. Consistent with current systematic schemes, analyses of complete cytochrome b genes including the shrew Sorex araneus (family Soricidae) revealed a close relationship between Talpidae and Soricidae. The analyses of complete mtDNAs, along with the findings of other insectivore studies, challenge the maintenance of the order Lipotyphla as a taxonomic unit and support the elevation of the Soricomorpha (with the families Talpidae and Soricidae and possibly also the Solenodontidae and Tenrecidae) to the level of an order, as previously proposed in some morphological studies.     

Phylogenetic assessment of molecular and morphological data for eutherian mammals

The interordinal relationships of eutherian (placental) mammals were evaluated by a phylogenetic analysis of four published data sets (three sequences and one morphological). The nature and degree of support and conflict for particular groups were assessed by separate bootstrap and homogeneity tests, which were followed by combined analyses of the sequence and morphological data. Between orders, strong support (i.e., > or = 95% bootstrap scores) was found for a paraphyletic Artiodactyla (relative to Cetacea) and a monophyletic Cetartiodactyla (Artiodactyla and Cetacea) and Paenungulata (Hyracoidea, Proboscidea, and Sirenia). In turn, some reasonable to strong evidence (> or = 85%) was obtained for Hyracoidea with Sirenia, Dermoptera with Scandentia, Glires (Lagomorpha with Rodentia), and Afrotheria (Amblysomus, Macroscelidea, Paenungulata, and Tubulidentata). Otherwise, no other interordinal clades were supported at these reasonable to strong levels. This overall lack of resolution for eutherian interordinal clusters agrees with other studies that suggest further progress will continue to be slow and difficult. Further resolution will require the integration of more recently published data, the continued sampling of taxa and characters, and the use of more powerful methods of data analysis.     

PHYLOGENY OF PRIMATES AND OTHER EUTHERIAN ORDERS: A CLADISTIC ANALYSIS USING AMINO ACID AND NUCLEOTIDE SEQUENCE DATA

Abstract— Genealogical reconstructions carried out by the parsimony method on protein amino acid and DNA nucleotide sequence data are providing fresh evidence on cladistic branching patterns at taxonomic levels from the classes of Vertebrata and orders of Eutheria to the genera of Hominoidea. Minimum length trees constructed from amino acid sequence data group Mammalia with Archosauria (i.e., Aves plus Crocodilia), Amniota with Amphibia, and Tetrapoda with Teleostei. Within Mammalia, Edentata and Paenungulata (e.g., Proboscidea) appear as the most anciently separated from other eutherians. Another superordinal eutherian clade consists of Artiodactyla, Cetacea, and Perissodactyla. A third consistently contains Primates, Lagomorpha, and Tupaia. The cladistic positions of such orders as Carnivora, Chiroptera, Insectivora, and Rodentia are not well resolved by the currently still sparse body of sequence data. However, recent dramatic progress in the technology of gene cloning and nucleotide sequencing has opened the way for so enlarging the body of sequence data that it should become possible to solve almost any problem concerning the phylogenetic systematice of extant mammals. An example is provided by hominoid genera. Minimum length trees constructed from mitochondrial DNA nucleotide sequence data very strongly group Pan, Homo , and Gorilla into Homininae and then join Homininae and Ponginae (pongo) into Hominidae as the sister family of Hylobatidae (Hylobates). Resolution of the hominine trichotomy into two dichotomous branchings should be forthcoming as kilobase sequencing of nuclear genes progresses.     

A new phylogenetic marker,apolipoprotein B,provides compelling evidence for eutherian relationships

Higher-level relationships within, and the root of Placentalia, remain contentious issues. Resolution of the placental tree is important to the choice of mammalian genome projects and model organisms, as well as for understanding the biogeography of the eutherian radiation. We present phylogenetic analyses of 63 species representing all extant eutherian mammal orders for a new molecular phylogenetic marker, a 1.3kb portion of exon 26 of the apolipoprotein B (APOB) gene. In addition, we analyzed a multigene concatenation that included APOB sequences and a previously published data set (Murphy et al., 2001b) of three mitochondrial and 19 nuclear genes, resulting in an alignment of over 17kb for 42 placentals and two marsupials. Due to computational difficulties, previous maximum likelihood analyses of large, multigene concatenations for placental mammals have used quartet puzzling, less complex models of sequence evolution, or phylogenetic constraints to approximate a full maximum likelihood bootstrap. Here, we utilize a Unix load sharing facility to perform maximum likelihood bootstrap analyses for both the APOB and concatenated data sets with a GTR+Gamma+I model of sequence evolution, tree-bisection and reconnection branch-swapping, and no phylogenetic constraints. Maximum likelihood and Bayesian analyses of both data sets provide support for the superordinal clades Boreoeutheria, Euarchontoglires, Laurasiatheria, Xenarthra, Afrotheria, and Ostentoria (pangolins+carnivores), as well as for the monophyly of the orders Eulipotyphla, Primates, and Rodentia, all of which have recently been questioned. Both data sets recovered an association of Hippopotamidae and Cetacea within Cetartiodactyla, as well as hedgehog and shrew within Eulipotyphla. APOB showed strong support for an association of tarsier and Anthropoidea within Primates. Parsimony, maximum likelihood and Bayesian analyses with both data sets placed Afrotheria at the base of the placental radiation. Statistical tests that employed APOB to examine a priori hypotheses for the root of the placental tree rejected rooting on myomorphs and hedgehog, but did not discriminate between rooting at the base of Afrotheria, at the base of Xenarthra, or between Atlantogenata (Xenarthra+Afrotheria) and Boreoeutheria. An orthologous deletion of 363bp in the aligned APOB sequences proved phylogenetically informative for the grouping of the order Carnivora with the order Pholidota into the superordinal clade Ostentoria. A smaller deletion of 237-246bp was diagnostic of the superordinal clade Afrotheria.     

Mammalian mitochondrial DNA evolution: A comparison of the cytochrome b and cytochrome c oxidase II genes

The evolution of two mitochondrial genes, cytochrome b and cytochrome c oxidase subunit II, was examined in several eutherian mammal orders, with special emphasis on the orders Artiodactyla and Rodentia. When analyzed using both maximum parsimony, with either equal or unequal character weighting, and neighbor joining, neither gene performed with a high degree of consistency in terms of the phylogenetic hypotheses supported. The phylogenetic inconsistencies observed for both these genes may be the result of several factors including differences in the rate of nucleotide substitution among particular lineages (especially between orders), base composition bias, transition/transversion bias, differences in codon usage, and different constraints and levels of homoplasy associated with first, second, and third codon positions. We discuss the implications of these findings for the molecular systematics of mammals, especially as they relate to recent hypotheses concerning the polyphyly of the order Rodentia, relationships among the Artiodactyla, and various interordinal relationships.Correspondence to: R.L. Honeycutt     

Mammalian evolution and the interphotoreceptor retinoid binding protein (IRBP) gene: Convincing evidence for several superordinal clades

Phylogenetic relationships of 25 mammalian species representing 17 of the 18 eutherian orders were examined using DNA sequences from a 1.2-kb region of the 5′ end of exon 1 of the single-copy nuclear gene known as interphotoreceptor retinoid binding protein (IRBP). A wide variety of methods of analysis of the DNA sequence, and of the translated products, all supported a five-order clade consisting of elephant shrew (Macroscelidea)/aardvark (Tubulidentata)/and the paenungulates (hyracoids, sirenians, and elephants), with bootstrap support in all cases of 100%. The Paenungulata was also strongly supported by these IRBP data. In the majority of analyses this monophyletic five-order grouping was the first branch off the tree after the Edentata. These results are highly congruent with two other recent sources of molecular data. Another superordinal grouping, with similar 100% bootstrap support in all of the same wide-ranging types of analyses, was Artiodactyla/Cetacea. Other superordinal affinities, suggested by the analyses, but with less convincing support, included a Perissodactyla/Artiodactyla/Cetacea clade, an Insectivora/Chiroptera clade, and Glires (an association of rodents and lagomorphs). Correspondence to: M. J. Stanhope     

The eutherian stapedial artery: character analysis and implications for superordinal relationships

Evidence from outgroups, ontogeny, neontology, and fossils is used to distinguish primitive and derived character states for the major components of the eutherian stapedial artery in 17 modern orders. Derived states support the following higher-level phylogenetic hypotheses: (1) Microchiroptera and Megachiroptera are monophyletic+ADs- and (2) within Ungulata, Tubulidentata is the outgroup to the remaining modern orders, followed in succession by Artiodactyla and then Cetacea. Three branches of the stapedial artery, the a. diploetica magna, ramus temporalis, and ramus posterior, all but neglected in previous syntheses, are shown to be primitive for Eutheria and Amniota.     

Basal metabolic rates in mammals: taxonomic differences in the allometry of BMR and body mass

No single equation adequately describes the allometric relation between body mass and BMR for mammals. Least squares regression of log-transformed data for 248 eutherian species results in a line with a slope (-0.30) significantly different from that of Kleiber's line (-0.25). Interordinal comparisons of least squares regressions of log-transformed BMR and mass suggest that the Insectivora have a significantly steeper slope to their allometric relationship than do most other orders, while the non-insectivore orders are statistically homogeneous with respect to slope. With respect to elevation, Edentata have the lowest BMRs; Marsupialia, Primates and Chiroptera are indistinguishable from each other but above the edentates; Primates, Chiroptera, Rodentia, Lagomorpha and Carnivora form the next highest homogeneous grouping; and Artiodactyla have the highest BMRs, significantly greater than all but Lagomorpha and Carnivora. Analysis of intraordinal variation within the Rodentia suggests significant heterogeneity among families in BMR-mass allometry.     

Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data

Extensive phylogenetic analyses of the updated sequence data of mammalian mitochondrial genomes were carried out using the maximum likelihood method in order to resolve deep branchings in eutherian evolution. The divergence times in the mammalian tree were estimated by a relaxed molecular clock of the mitochondrial proteins calibrated with multiple references. A Chiroptera/Eulipotyphla (i.e. bat/mole) clade and a close relationship of this clade to Fereuungulata (Carnivora+Perissodactyla+Cetartiodactyla) were reconfirmed with high statistical significance. However, a support for a monophyly of Fereuungulata relative to the Chiroptera/Eulipotyphla clade was fragile, and we suggest that the three branchings among Carnivora, Perissodactyla, Cetartiodactyla and Chiroptera/Eulipotyphla occurred successively in a short time period, estimated to be approximately 77Myr BP. The Chiroptera/Eulipotyphla divergence was estimated to roughly coincide with the Cretaceous-Tertiary boundary (65Myr BP). The monophyly of Rodentia, the Lagomorpha/Rodentia clade (traditionally called Glires), and the Afrotheria/Xenarthra clade were preferred over alternative relationships, but the supports of these clades were not strong enough to exclude other possibilities. Although several super-order taxa of eutherians were strongly supported by the analyses of the mitochondrial genome data, the branching order in the deepest part of the eutherian tree remained ambiguous from the data presently available.     

Model dependence of the phylogenetic inference: relationship among carnivores, Perissodactyls and cetartiodactyls as inferred from mitochondrial genome sequences

Some previous analysis of mitochondrial proteins strongly support the Carnivora/Perissodactyla grouping excluding Cetartiodactyla (Artiodactyla + Cetacea) as an outgroup, but the support of the hypothesis remains equivocal from the analysis of several nuclear-encoded proteins. In order to evaluate the strength of the support by mitochondrial proteins, phylogenetic relationship among Carnivora, Perissodactyla, and Cetartiodactyla was estimated with the ML method by using the updated data set of the 12 mitochondrial proteins with several alternative models. The analyses demonstrate that the phylogenetic inference depends on the model used in the ML analysis; i.e., whether the site-heterogeneity is taken into account and whether the rate parameters are estimated for each individual proteins or for the concatenated sequences. Although the analysis of concatenated sequences strongly supports the Carnivora/Perissodactyla grouping, the total evaluation of the separate analyses of individual proteins, which approximates the data better than the concatenated analysis, gives only ambiguous results, and therefore it is concluded that more data are needed to resolve this trichotomy.     

Phylogenetic position of cetaceans relative to artiodactyls: reanalysis of mitochondrial and nuclear sequences

By a maximum likelihood analysis of mitochondrial DNA sequences, we examine Graur and Higgins' hypothesis of the Ruminantia/Cetacea clade with Suiformes as an outgroup. Graur and Higgins analyzed these sequences by the neighbor-joining and parsimony methods, as well as by the maximum likelihood method under the assumption that the substitution rate is the same for all sites. The Ruminantia/Suiformes clade assumed by the traditional taxonomy was rejected strongly by this analysis and the Ruminantia/Cetacea clade was supported. Adoption of a more realistic model distinguishing among rates at different codon positions in the maximum likelihood analysis of the same data, however, grossly reduces the significance level on the Graur-Higgins hypothesis. Thus, although the Ruminantia/Suiformes grouping is indeed least likely from Graur and Higgins' data set of mitochondrial DNA, this traditional tree cannot be rejected with statistical significance under the new analysis, and more data are needed to settle the issue. In the same way, we examine Irwin and Arnason's suggestion of the Hippopotamus/Cetacea clade by using cytochrome b and hemoglobins alpha and beta, and it turn out that their suggestion is also fragile. This analysis demonstrates the importance of selecting an appropriate model among the alternatives in the maximum likelihood analysis and of using many different genes from many relevant species in order to make reliable phylogenetic inferences.      

c-myc gene sequences and the phylogeny of bats and other eutherian mammals

The complete protein-coding sequences of the c-myc proto-oncogene were determined for five species of four new orders of eutherian (placental) mammals. These newly obtained sequences were aligned to each other and to other available orthologs for the phylogenetic estimation of eutherian interordinal relationships. Several measures of sequence difference and base composition were first calculated to assess the major evolutionary properties of the three codon positions and two protein-coding exons of the gene. On the basis of these calculations, different parsimony, distance, and maximum likelihood approaches were adopted, with the most sophisticated involving the separate, then combined, likelihood analyses of the third codon positions of exon 2 versus all other sites. These phylogenetic approaches provided clear support for the grouping of Chiroptera (bats) with Artiodactyla (ruminants, camels, and pigs) and Carnivora (cats, dogs, and their allies), an interordinal arrangement that receives strong corroboration from other lines of evidence including complete mitochondrial DNA sequences. In contrast, these analyses failed to provide strong to reasonable support for any other interordinal group. This study concludes with specific recommendations about sampling and other strategies for maximizing the phylogenetic contributions of the c-myc gene to the continued resolution of the eutherian ordinal tree.     

Mitogenomic analyses of eutherian relationships

Reasonably correct phylogenies are fundamental to the testing of evolutionary hypotheses. Here, we present phylogenetic findings based on analyses of 67 complete mammalian mitochondrial (mt) genomes. The analyses, irrespective of whether they were performed at the amino acid (aa) level or on nucleotides (nt) of first and second codon positions, placed Erinaceomorpha (hedgehogs and their kin) as the sister group of remaining eutherians. Thus, the analyses separated Erinaceomorpha from other traditional lipotyphlans (e.g., tenrecs, moles, and shrews), making traditional Lipotyphla polyphyletic. Both the aa and nt data sets identified the two order-rich eutherian clades, the Cetferungulata (comprising Pholidota, Carnivora, Perissodactyla, Artiodactyla, and Cetacea) and the African clade (Tenrecomorpha, Macroscelidea, Tubulidentata, Hyracoidea, Proboscidea, and Sirenia). The study corroborated recent findings that have identified a sister-group relationship between Anthropoidea and Dermoptera (flying lemurs), thereby making our own order, Primates, a paraphyletic assembly. Molecular estimates using paleontologically well-established calibration points, placed the origin of most eutherian orders in Cretaceous times, 70-100 million years before present (MYBP). The same estimates place all primate divergences much earlier than traditionally believed. For example, the divergence between Homo and Pan is estimated to have taken place approximately 10 MYBP, a dating consistent with recent findings in primate paleontology.     

The position of Cetacea within mammalia: phylogenetic analysis of morphological data from extinct and extant taxa

Knowledge of the phylogenetic position of the order Cetacea (whales, dolphins, and porpoises) within Mammalia is of central importance to evolutionary biologists studying the transformations of biological form and function that accompanied the shift from fully terrestrial to fully aquatic life in this clade. Phylogenies based on molecular data and those based on morphological data both place cetaceans among ungulates but are incongruent in other respects. Morphologists argue that cetaceans are most closely related to mesonychians, an extinct group of terrestrial ungulates. They have disagreed, however, as to whether Perissodactyla (odd-toed ungulates) or Artiodactyla (even-toed ungulates) is the extant clade most closely related to Cetacea, and have long maintained that each of these orders is monophyletic. The great majority of molecule-based phylogenies show, by contrast, not only that artiodactyls are the closest extant relatives of Cetacea, but also that Artiodactyla is paraphyletic unless cetaceans are nested within it, often as the sister group of hippopotamids. We tested morphological evidence for several hypotheses concerning the sister taxon relationships of Cetacea in a maximum parsimony analysis of 123 morphological characters from 10 extant and 30 extinct taxa. We advocate treating certain multistate characters as ordered because such a procedure incorporates information about hierarchical morphological transformation. In all most-parsimonious trees, whether multistate characters are ordered or unordered, Artiodactyla is the extant sister taxon of Cetacea. With certain multistate characters ordered, the extinct clade Mesonychia (Mesonychidae + Hapalodectidae) is the sister taxon of Cetacea, and Artiodactyla is monophyletic. When all fossils are removed from the analysis, Artiodactyla is paraphyletic with Cetacea nested inside, indicating that inclusion of mesonychians and other extinct stem taxa in a phylogenetic analysis of the ungulate clade is integral to the recovery of artiodactyl monophyly. Phylogenies derived from molecular data alone may risk recovering inconsistent branches because of an inability to sample extinct clades, which by a conservative estimate, amount to 89% of the ingroup. Addition of data from recently described astragali attributed to cetaceans does not overturn artiodactyl monophyly.