Seeing through the eyes of the sabertooth Thylacosmilus atrox (Metatheria, Sparassodonta).

The evolution of mammalian vision is difficult to study because the actual receptor organs-the eyes-are not preserved in the fossil record. Orbital orientation and size are the traditional proxies for inferring aspects of ocular function, such as stereoscopy. Adaptations for good stereopsis have evolved in living predaceous mammals, and it is reasonable to infer that fossil representatives would follow the same pattern. This applies to the sparassodonts, an extinct group of South American hypercarnivores related to marsupials, with one exception. In the sabertooth Thylacosmilus atrox, the bony orbits were notably divergent, like those of a cow or a horse, and thus radically differing from conditions in any other known mammalian predator. Orbital convergence alone, however, does not determine presence of stereopsis; frontation and verticality of the orbits also play a role. We show that the orbits of Thylacosmilus were frontated and verticalized in a way that favored some degree of stereopsis and compensated for limited convergence in orbital orientation. The forcing function behind these morphological tradeoffs was the extraordinary growth of its rootless canines, which affected skull shape in Thylacosmilus in numerous ways, including relative orbital displacement.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments.

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.1 In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.2 Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia3-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,6-12 in paleodietary reconstructions,13-19 and facilitated genomic reconstructions from cave contexts.8,20-22 Murchie et al.6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.24 Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 40,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

Collapse of the mammoth-steppe in central Yukon as revealed by ancient environmental DNA.

The temporal and spatial coarseness of megafaunal fossil records complicates attempts to to disentangle the relative impacts of climate change, ecosystem restructuring, and human activities associated with the Late Quaternary extinctions. Advances in the extraction and identification of ancient DNA that was shed into the environment and preserved for millennia in sediment now provides a way to augment discontinuous palaeontological assemblages. Here, we present a 30,000-year sedimentary ancient DNA (sedaDNA) record derived from loessal permafrost silts in the Klondike region of Yukon, Canada. We observe a substantial turnover in ecosystem composition between 13,500 and 10,000 calendar years ago with the rise of woody shrubs and the disappearance of the mammoth-steppe (steppe-tundra) ecosystem. We also identify a lingering signal of Equus sp. (North American horse) and Mammuthus primigenius (woolly mammoth) at multiple sites persisting thousands of years after their supposed extinction from the fossil record.

No evidence for widespread island extinctions after Pleistocene hominin arrival.

The arrival of modern humans into previously unoccupied island ecosystems is closely linked to widespread extinction, and a key reason cited for Pleistocene megafauna extinction is anthropogenic overhunting. A common assumption based on late Holocene records is that humans always negatively impact insular biotas, which requires an extrapolation of recent human behavior and technology into the archaeological past. Hominins have been on islands since at least the early Pleistocene and Homo sapiens for at least 50 thousand y (ka). Over such lengthy intervals it is scarcely surprising that significant evolutionary, behavioral, and cultural changes occurred. However, the deep-time link between human arrival and island extinctions has never been explored globally. Here, we examine archaeological and paleontological records of all Pleistocene islands with a documented hominin presence to examine whether humans have always been destructive agents. We show that extinctions at a global level cannot be associated with Pleistocene hominin arrival based on current data and are difficult to disentangle from records of environmental change. It is not until the Holocene that large-scale changes in technology, dispersal, demography, and human behavior visibly affect island ecosystems. The extinction acceleration we are currently experiencing is thus not inherent but rather part of a more recent cultural complex.

New radiometric 40Ar-39Ar dates and faunistic analyses refine evolutionary dynamics of Neogene vertebrate assemblages in southern South America.

The vertebrate fossil record of the Pampean Region of Argentina occupies an important place in South American vertebrate paleontology. An abundance of localities has long been the main basis for constructing the chronostratigraphical/geochronological scale for the late Neogene-Quaternary of South America, as well as for understanding major patterns of vertebrate evolution, including the Great American Biotic Interchange. However, few independently-derived dates are available for constraining this record. In this contribution, we present new 40Ar/39Ar dates on escorias (likely the product of meteoric impacts) from the Argentinean Atlantic coast and statistically-based biochronological analyses that help to calibrate Late Miocene-Pliocene Pampean faunal successions. For the type areas of the Montehermosan and Chapadmalalan Ages/Stages, our results delimit their age ranges to 4.7-3.7 Ma and ca. 3.74-3.04 Ma, respectively. Additionally, from Buenos Aires Province, dates of 5.17 Ma and 4.33 Ma were recovered for "Huayquerian" and Montehermosan faunas. This information helps to better calibrate important first appearances of allochthonous taxa in South America, including one of the oldest records for procyonids (7.24-5.95 Ma), cricetids (6.95-5.46 Ma), and tayassuids (> 3.74 Ma, oldest high-confidence record). These results also constrain to ca. 3 Ma the last appearances of the autochthonous sparassodonts, as well as terror birds of large/middle body size in South America. South American faunal turnover during the late Neogene, including Late Pliocene extinctions, is interpreted as a consequence of knock-on effects from global climatic changes and initiation of the icehouse climate regime.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears.

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

American mastodon mitochondrial genomes suggest multiple dispersal events in response to Pleistocene climate oscillations.

Pleistocene glacial-interglacial cycles are correlated with dramatic temperature oscillations. Examining how species responded to these natural fluctuations can provide valuable insights into the impacts of present-day anthropogenic climate change. Here we present a phylogeographic study of the extinct American mastodon (Mammut americanum), based on 35 complete mitochondrial genomes. These data reveal the presence of multiple lineages within this species, including two distinct clades from eastern Beringia. Our molecular date estimates suggest that these clades arose at different times, supporting a pattern of repeated northern expansion and local extirpation in response to glacial cycling. Consistent with this hypothesis, we also note lower levels of genetic diversity among northern mastodons than in endemic clades south of the continental ice sheets. The results of our study highlight the complex relationships between population dispersals and climate change, and can provide testable hypotheses for extant species expected to experience substantial biogeographic impacts from rising temperatures.

Ancient Mitogenomes Reveal the Evolutionary History and Biogeography of Sloths.

Living sloths represent two distinct lineages of small-sized mammals that independently evolved arboreality from terrestrial ancestors. The six extant species are the survivors of an evolutionary radiation marked by the extinction of large terrestrial forms at the end of the Quaternary. Until now, sloth evolutionary history has mainly been reconstructed from phylogenetic analyses of morphological characters. Here, we used ancient DNA methods to successfully sequence 10 extinct sloth mitogenomes encompassing all major lineages. This includes the iconic continental ground sloths Megatherium, Megalonyx, Mylodon, and Nothrotheriops and the smaller endemic Caribbean sloths Parocnus and Acratocnus. Phylogenetic analyses identify eight distinct lineages grouped in three well-supported clades, whose interrelationships are markedly incongruent with the currently accepted morphological topology. We show that recently extinct Caribbean sloths have a single origin but comprise two highly divergent lineages that are not directly related to living two-fingered sloths, which instead group with Mylodon. Moreover, living three-fingered sloths do not represent the sister group to all other sloths but are nested within a clade of extinct ground sloths including Megatherium, Megalonyx, and Nothrotheriops. Molecular dating also reveals that the eight newly recognized sloth families all originated between 36 and 28 million years ago (mya). The early divergence of recently extinct Caribbean sloths around 35 mya is consistent with the debated GAARlandia hypothesis postulating the existence at that time of a biogeographic connection between northern South America and the Greater Antilles. This new molecular phylogeny has major implications for reinterpreting sloth morphological evolution, biogeography, and diversification history.

Palaeoproteomics resolves sloth relationships.

The living tree sloths Choloepus and Bradypus are the only remaining members of Folivora, a major xenarthran radiation that occupied a wide range of habitats in many parts of the western hemisphere during the Cenozoic, including both continents and the West Indies. Ancient DNA evidence has played only a minor role in folivoran systematics, as most sloths lived in places not conducive to genomic preservation. Here we utilize collagen sequence information, both separately and in combination with published mitochondrial DNA evidence, to assess the relationships of tree sloths and their extinct relatives. Results from phylogenetic analysis of these datasets differ substantially from morphology-based concepts: Choloepus groups with Mylodontidae, not Megalonychidae; Bradypus and Megalonyx pair together as megatherioids, while monophyletic Antillean sloths may be sister to all other folivorans. Divergence estimates are consistent with fossil evidence for mid-Cenozoic presence of sloths in the West Indies and an early Miocene radiation in South America.

Ancient DNA of the extinct Jamaican monkey Xenothrix reveals extreme insular change within a morphologically conservative radiation.

The insular Caribbean until recently contained a diverse mammal fauna including four endemic platyrrhine primate species, all of which died out during the Holocene. Previous morphological studies have attempted to establish how these primates are related to fossil and extant platyrrhines, whether they represent ancient or recent colonists, and whether they constitute a monophyletic group. These efforts have generated multiple conflicting hypotheses, from close sister-taxon relationships with several different extant platyrrhines to derivation from a stem platyrrhine lineage outside the extant Neotropical radiation. This diversity of opinion reflects the fact that Caribbean primates were morphologically extremely unusual, displaying numerous autapomorphies and apparently derived conditions present across different platyrrhine clades. Here we report ancient DNA data for an extinct Caribbean primate: a limited-coverage entire mitochondrial genome and seven regions of nuclear genome for the most morphologically derived taxon, the Jamaican monkey Xenothrix mcgregori We demonstrate that Xenothrix is part of the existing platyrrhine radiation rather than a late-surviving stem platyrrhine, despite its unusual adaptations, and falls within the species-rich but morphologically conservative titi monkey clade (Callicebinae) as sister to the newly recognized genus Cheracebus These results are not congruent with previous morphology-based hypotheses and suggest that even morphologically conservative lineages can exhibit phenetic plasticity in novel environments like those found on islands. Xenothrix and Cheracebus diverged ca. 11 Ma, but primates have been present in the Caribbean since 17.5-18.5 Ma, indicating that Caribbean primate diversity was generated by multiple over-water colonizations.

A comprehensive genomic history of extinct and living elephants.

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

A mitogenomic timetree for Darwin's enigmatic South American mammal Macrauchenia patachonica.

The unusual mix of morphological traits displayed by extinct South American native ungulates (SANUs) confounded both Charles Darwin, who first discovered them, and Richard Owen, who tried to resolve their relationships. Here we report an almost complete mitochondrial genome for the litoptern Macrauchenia. Our dated phylogenetic tree places Macrauchenia as sister to Perissodactyla, but close to the radiation of major lineages within Laurasiatheria. This position is consistent with a divergence estimate of ∼66 Ma (95% credibility interval, 56.64-77.83 Ma) obtained for the split between Macrauchenia and other Panperissodactyla. Combined with their morphological distinctiveness, this evidence supports the positioning of Litopterna (possibly in company with other SANU groups) as a separate order within Laurasiatheria. We also show that, when using strict criteria, extinct taxa marked by deep divergence times and a lack of close living relatives may still be amenable to palaeogenomic analysis through iterative mapping against more distant relatives.

Fossil and genomic evidence constrains the timing of bison arrival in North America.

The arrival of bison in North America marks one of the most successful large-mammal dispersals from Asia within the last million years, yet the timing and nature of this event remain poorly determined. Here, we used a combined paleontological and paleogenomic approach to provide a robust timeline for the entry and subsequent evolution of bison within North America. We characterized two fossil-rich localities in Canada's Yukon and identified the oldest well-constrained bison fossil in North America, a 130,000-y-old steppe bison, Bison cf. priscus We extracted and sequenced mitochondrial genomes from both this bison and from the remains of a recently discovered, ∼120,000-y-old giant long-horned bison, Bison latifrons, from Snowmass, Colorado. We analyzed these and 44 other bison mitogenomes with ages that span the Late Pleistocene, and identified two waves of bison dispersal into North America from Asia, the earliest of which occurred ∼195-135 thousand y ago and preceded the morphological diversification of North American bison, and the second of which occurred during the Late Pleistocene, ∼45-21 thousand y ago. This chronological arc establishes that bison first entered North America during the sea level lowstand accompanying marine isotope stage 6, rejecting earlier records of bison in North America. After their invasion, bison rapidly colonized North America during the last interglaciation, spreading from Alaska through continental North America; they have been continuously resident since then.

The evolutionary and phylogeographic history of woolly mammoths: a comprehensive mitogenomic analysis.

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

Evolutionary History of the Nesophontidae, the Last Unplaced Recent Mammal Family.

The mammalian evolutionary tree has lost several major clades through recent human-caused extinctions. This process of historical biodiversity loss has particularly affected tropical island regions such as the Caribbean, an area of great evolutionary diversification but poor molecular preservation. The most enigmatic of the recently extinct endemic Caribbean mammals are the Nesophontidae, a family of morphologically plesiomorphic lipotyphlan insectivores with no consensus on their evolutionary affinities, and which constitute the only major recent mammal clade to lack any molecular information on their phylogenetic placement. Here, we use a palaeogenomic approach to place Nesophontidae within the phylogeny of recent Lipotyphla. We recovered the near-complete mitochondrial genome and sequences for 17 nuclear genes from a ∼750-year-old Hispaniolan Nesophontes specimen, and identify a divergence from their closest living relatives, the Solenodontidae, more than 40 million years ago. Nesophontidae is thus an older distinct lineage than many extant mammalian orders, highlighting not only the role of island systems as "museums" of diversity that preserve ancient lineages, but also the major human-caused loss of evolutionary history.

The phylogenetic affinities of the extinct glyptodonts.

Among the fossils of hitherto unknown mammals that Darwin collected in South America between 1832 and 1833 during the Beagle expedition were examples of the large, heavily armored herbivores later known as glyptodonts. Ever since, glyptodonts have fascinated evolutionary biologists because of their remarkable skeletal adaptations and seemingly isolated phylogenetic position even within their natural group, the cingulate xenarthrans (armadillos and their allies). In possessing a carapace comprised of fused osteoderms, the glyptodonts were clearly related to other cingulates, but their precise phylogenetic position as suggested by morphology remains unresolved. To provide a molecular perspective on this issue, we designed sequence-capture baits using in silico reconstructed ancestral sequences and successfully assembled the complete mitochondrial genome of Doedicurus sp., one of the largest glyptodonts. Our phylogenetic reconstructions establish that glyptodonts are in fact deeply nested within the armadillo crown-group, representing a distinct subfamily (Glyptodontinae) within family Chlamyphoridae. Molecular dating suggests that glyptodonts diverged no earlier than around 35 million years ago, in good agreement with their fossil record. Our results highlight the derived nature of the glyptodont morphotype, one aspect of which is a spectacular increase in body size until their extinction at the end of the last ice age.

Evolutionary Relationships among Extinct and Extant Sloths: The Evidence of Mitogenomes and Retroviruses.

Macroevolutionary trends exhibited by retroviruses are complex and not entirely understood. The sloth endogenized foamy-like retrovirus (SloEFV), which demonstrates incongruence in virus-host evolution among extant sloths (Order Folivora), has not been investigated heretofore in any extinct sloth lineages and its premodern history within folivorans is therefore unknown. Determining retroviral coevolutionary trends requires a robust phylogeny of the viral host, but the highly reduced modern sloth fauna (6 species in 2 genera) does not adequately represent what was once a highly diversified clade (∼100 genera) of placental mammals. At present, the amount of molecular data available for extinct sloth taxa is limited, and analytical results based on these data tend to conflict with phylogenetic inferences made on the basis of morphological studies. To augment the molecular data set, we applied hybridization capture and next-generation Illumina sequencing to two extinct and three extant sloth species to retrieve full mitochondrial genomes (mitogenomes) from the hosts and the polymerase gene of SloEFV. The results produced a fully resolved and well-supported phylogeny that supports dividing crown families into two major clades: 1) The three-toed sloth, Bradypus, and Nothrotheriidae and 2) Megalonychidae, including the two-toed sloth, Choloepus, and Mylodontidae. Our calibrated time tree indicates that the Miocene epoch (23.5 Ma), particularly its earlier part, was an important interval for folivoran diversification. Both extant and extinct sloths demonstrate multiple complex invasions of SloEFV into the ancestral sloth germline followed by subsequent introgressions across different sloth lineages. Thus, sloth mitogenome and SloEFV evolution occurred separately and in parallel among sloths.

Interordinal gene capture, the phylogenetic position of Steller's sea cow based on molecular and morphological data, and the macroevolutionary history of Sirenia.

The recently extinct (ca. 1768) Steller's sea cow (Hydrodamalis gigas) was a large, edentulous North Pacific sirenian. The phylogenetic affinities of this taxon to other members of this clade, living and extinct, are uncertain based on previous morphological and molecular studies. We employed hybridization capture methods and second generation sequencing technology to obtain >30kb of exon sequences from 26 nuclear genes for both H. gigas and Dugong dugon. We also obtained complete coding sequences for the tooth-related enamelin (ENAM) gene. Hybridization probes designed using dugong and manatee sequences were both highly effective in retrieving sequences from H. gigas (mean=98.8% coverage), as were more divergent probes for regions of ENAM (99.0% coverage) that were designed exclusively from a proboscidean (African elephant) and a hyracoid (Cape hyrax). New sequences were combined with available sequences for representatives of all other afrotherian orders. We also expanded a previously published morphological matrix for living and fossil Sirenia by adding both new taxa and nine new postcranial characters. Maximum likelihood and parsimony analyses of the molecular data provide robust support for an association of H. gigas and D. dugon to the exclusion of living trichechids (manatees). Parsimony analyses of the morphological data also support the inclusion of H. gigas in Dugongidae with D. dugon and fossil dugongids. Timetree analyses based on calibration density approaches with hard- and soft-bounded constraints suggest that H. gigas and D. dugon diverged in the Oligocene and that crown sirenians last shared a common ancestor in the Eocene. The coding sequence for the ENAM gene in H. gigas does not contain frameshift mutations or stop codons, but there is a transversion mutation (AG to CG) in the acceptor splice site of intron 2. This disruption in the edentulous Steller's sea cow is consistent with previous studies that have documented inactivating mutations in tooth-specific loci of a variety of edentulous and enamelless vertebrates including birds, turtles, aardvarks, pangolins, xenarthrans, and baleen whales. Further, branch-site dN/dS analyses provide evidence for positive selection in ENAM on the stem dugongid branch where extensive tooth reduction occurred, followed by neutral evolution on the Hydrodamalis branch. Finally, we present a synthetic evolutionary tree for living and fossil sirenians showing several key innovations in the history of this clade including character state changes that parallel those that occurred in the evolutionary history of cetaceans.

Genomic Data from Extinct North American Camelops Revise Camel Evolutionary History.

Recent advances in paleogenomic technologies have enabled an increasingly detailed understanding of the evolutionary relationships of now-extinct mammalian taxa. However, a number of enigmatic Quaternary species have never been characterized with molecular data, often because available fossils are rare or are found in environments that are not optimal for DNA preservation. Here, we analyze paleogenomic data extracted from bones attributed to the late Pleistocene western camel, Camelops cf. hesternus, a species that was distributed across central and western North America until its extinction approximately 13,000 years ago. Despite a modal sequence length of only around 35 base pairs, we reconstructed high-coverage complete mitochondrial genomes and low-coverage partial nuclear genomes for each specimen. We find that Camelops is sister to African and Asian bactrian and dromedary camels, to the exclusion of South American camelids (llamas, guanacos, alpacas, and vicuñas). These results contradict previous morphology-based phylogenetic models for Camelops, which suggest instead a closer relationship between Camelops and the South American camelids. The molecular data imply a Late Miocene divergence of the Camelops clade from lineages that separately gave rise to the extant camels of Eurasia. Our results demonstrate the increasing capacity of modern paleogenomic methods to resolve evolutionary relationships among distantly related lineages.

Ancient proteins resolve the evolutionary history of Darwin's South American ungulates.

No large group of recently extinct placental mammals remains as evolutionarily cryptic as the approximately 280 genera grouped as 'South American native ungulates'. To Charles Darwin, who first collected their remains, they included perhaps the 'strangest animal[s] ever discovered'. Today, much like 180 years ago, it is no clearer whether they had one origin or several, arose before or after the Cretaceous/Palaeogene transition 66.2 million years ago, or are more likely to belong with the elephants and sirenians of superorder Afrotheria than with the euungulates (cattle, horses, and allies) of superorder Laurasiatheria. Morphology-based analyses have proved unconvincing because convergences are pervasive among unrelated ungulate-like placentals. Approaches using ancient DNA have also been unsuccessful, probably because of rapid DNA degradation in semitropical and temperate deposits. Here we apply proteomic analysis to screen bone samples of the Late Quaternary South American native ungulate taxa Toxodon (Notoungulata) and Macrauchenia (Litopterna) for phylogenetically informative protein sequences. For each ungulate, we obtain approximately 90% direct sequence coverage of type I collagen α1- and α2-chains, representing approximately 900 of 1,140 amino-acid residues for each subunit. A phylogeny is estimated from an alignment of these fossil sequences with collagen (I) gene transcripts from available mammalian genomes or mass spectrometrically derived sequence data obtained for this study. The resulting consensus tree agrees well with recent higher-level mammalian phylogenies. Toxodon and Macrauchenia form a monophyletic group whose sister taxon is not Afrotheria or any of its constituent clades as recently claimed, but instead crown Perissodactyla (horses, tapirs, and rhinoceroses). These results are consistent with the origin of at least some South American native ungulates from 'condylarths', a paraphyletic assembly of archaic placentals. With ongoing improvements in instrumentation and analytical procedures, proteomics may produce a revolution in systematics such as that achieved by genomics, but with the possibility of reaching much further back in time.