Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

Ancient genomes reveal over two thousand years of dingo population structure.

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge.

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Genetic evidence for post-glacial expansion from a southern refugium in the eastern moa (Emeus crassus).

Cycles of glacial expansion and contraction throughout the Pleistocene drove increases and decreases, respectively, in the geographical range and population size of many animal species. Genetic data have revealed that during glacial maxima the distribution of many Eurasian animals was restricted to small refugial areas, from which species expanded to reoccupy parts of their former range as the climate warmed. It has been suggested that the extinct eastern moa (Emeus crassus)-a large, flightless bird from New Zealand-behaved analogously during glacial maxima, possibly surviving only in a restricted area of lowland habitat in the southern South Island of New Zealand during the Last Glacial Maximum (LGM). However, previous studies have lacked the power and geographical sampling to explicitly test this hypothesis using genetic data. Here we analyse 46 ancient mitochondrial genomes from Late Pleistocene and Holocene bones of the eastern moa from across their post-LGM distribution. Our results are consistent with a post-LGM increase in the population size and genetic diversity of eastern moa. We also demonstrate that genetic diversity was higher in eastern moa from the southern extent of their range, supporting the hypothesis that they expanded from a single glacial refugium following the LGM.

Widespread male sex bias in mammal fossil and museum collections.

A recent study of mammoth subfossil remains has demonstrated the potential of using relatively low-coverage high-throughput DNA sequencing to genetically sex specimens, revealing a strong male-biased sex ratio [P. Pecnerová et al., Curr. Biol. 27, 3505-3510.e3 (2017)]. Similar patterns were predicted for steppe bison, based on their analogous female herd-based structure. We genetically sexed subfossil remains of 186 Holarctic bison (Bison spp.), and also 91 brown bears (Ursus arctos), which are not female herd-based, and found that ∼75% of both groups were male, very close to the ratio observed in mammoths (72%). This large deviation from a 1:1 ratio was unexpected, but we found no evidence for sex differences with respect to DNA preservation, sample age, material type, or overall spatial distribution. We further examined ratios of male and female specimens from 4 large museum mammal collections and found a strong male bias, observable in almost all mammalian orders. We suggest that, in mammals at least, 1) wider male geographic ranges can lead to considerably increased chances of detection in fossil studies, and 2) sexual dimorphic behavior or appearance can facilitate a considerable sex bias in fossil and modern collections, on a previously unacknowledged scale. This finding has major implications for a wide range of studies of fossil and museum material.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Mitogenomes Uncover Extinct Penguin Taxa and Reveal Island Formation as a Key Driver of Speciation.

The emergence of islands has been linked to spectacular radiations of diverse organisms. Although penguins spend much of their lives at sea, they rely on land for nesting, and a high proportion of extant species are endemic to geologically young islands. Islands may thus have been crucial to the evolutionary diversification of penguins. We test this hypothesis using a fossil-calibrated phylogeny of mitochondrial genomes (mitogenomes) from all extant and recently extinct penguin taxa. Our temporal analysis demonstrates that numerous recent island-endemic penguin taxa diverged following the formation of their islands during the Plio-Pleistocene, including the Galápagos (Galápagos Islands), northern rockhopper (Gough Island), erect-crested (Antipodes Islands), Snares crested (Snares) and royal (Macquarie Island) penguins. Our analysis also reveals two new recently extinct island-endemic penguin taxa from New Zealand's Chatham Islands: Eudyptes warhami sp. nov. and a dwarf subspecies of the yellow-eyed penguin, Megadyptes antipodes richdalei ssp. nov. Eudyptes warhami diverged from the Antipodes Islands erect-crested penguin between 1.1 and 2.5 Ma, shortly after the emergence of the Chatham Islands (∼3 Ma). This new finding of recently evolved taxa on this young archipelago provides further evidence that the radiation of penguins over the last 5 Ma has been linked to island emergence. Mitogenomic analyses of all penguin species, and the discovery of two new extinct penguin taxa, highlight the importance of island formation in the diversification of penguins, as well as the extent to which anthropogenic extinctions have affected island-endemic taxa across the Southern Hemisphere's isolated archipelagos.

Reconstructing the Evolution of Giant Extinct Kangaroos: Comparing the Utility of DNA, Morphology, and Total Evidence.

Combined "total evidence" analysis of molecular and morphological data offers the opportunity to objectively merge fossils into the tree of life, and challenges the primacy of solely DNA based phylogenetic and dating inference, even among modern taxa. To investigate the relative utility of DNA, morphology, and total evidence for evolutionary inference, we sequenced the first near-complete mitochondrial genomes from extinct Australian megafauna: a 40-50 thousand year old giant short-faced kangaroo (Simosthenurus occidentalis) and giant wallaby (Protemnodon anak). We analyzed the ancient DNA and fossil data alongside comparable data from extant species to infer phylogeny, divergence times, and ancestral body mass among macropods (kangaroos and wallabies). Our results confirm a close relationship between Protemnodon and the iconic kangaroo genus complex "Macropus", and unite the giant Simothenurus with the hare-sized Lagostrophus fasciatus (banded hare-wallaby), suggesting that the latter is the closest living link to the once diverse sthenurine kangaroo radiation. We find that large body size evolved multiple times among kangaroos, coincident with expansion of open woodland habitats beginning in the Late Miocene. In addition, our results suggest that morphological data mislead macropod phylogeny reconstruction and in turn can distort total evidence estimation of divergence dates. However, a novel result with potentially broad application is that the accuracy and precision of reconstructing ancestral body mass was improved by tracing body mass on morphological branch lengths. This is likely due to positive allometric correlation between morphological and body size variation-a relationship that may be masked or even misleadingly inverted with the temporal or molecular branch lengths that typically underpin ancestral body size reconstruction. Our study supports complementary roles for DNA and morphology in evolutionary inference, and opens a new window into the evolution of Australia's unique marsupial fauna.

Oligocene divergence of frogmouth birds (Podargidae) across Wallace's Line.

Wallace's Line demarcates the transition between the differentiated regional faunas of Asia and Australia. However, while patterns of biotic differentiation across these two continental landmasses and the intervening island groups (Wallacea) have been extensively studied, patterns of long-term dispersal and diversification across this region are less well understood. Frogmouths (Aves: Podargidae) are a relictual family of large nocturnal birds represented by three extant genera occurring, respectively, in Asia, 'Sahul' (Australia and New Guinea) and the Solomon Islands, thus spanning Wallace's Line. We used new mitochondrial genomes from each of the extant frogmouth genera to estimate the timeline of frogmouth evolution and dispersal across Wallace's Line. Our results suggest that the three genera diverged and dispersed during the mid-Cenozoic between approximately 30 and 40 Mya. These divergences are among the oldest inferred for any trans-Wallacean vertebrate lineage. In addition, our results reveal that the monotypic Solomons frogmouth (Rigidipenna inexpectata) is one of the most phylogenetically divergent endemic bird lineages in the southwest Pacific. We suggest that the contemporary distribution of exceptionally deep divergences among extant frogmouth lineages may be explained by colonization of, and subsequent long-term persistence on, island arcs in the southwest Pacific during the Oligocene. These island arcs may have provided a pathway for biotic dispersal out of both Asia and Australia that preceded the formation of extensive emergent landmasses in Wallacea by at least 10 million years.

Molecular resolution to a morphological controversy: The case of North American fossil muskoxen Bootherium and Symbos.

The musk ox (Ovibos moschatus) is the only surviving member of a group of Pleistocene North American musk ox genera (Praeovibos, Ovibos, Bootherium, Euceratherium, and Soergelia) whose taxonomy is uncertain. The helmeted musk ox (Bootherium bombifrons) and the woodland musk ox (Symbos cavifrons) have been synonymised as male and female forms of a single Nearctic species found from Alaska, in the north, to Texas, in the south. However, this reclassification has not been tested using molecular data, despite the potential to use ancient DNA to examine these late Pleistocene taxa. In the present study, we sequenced mitochondrial genomes from seven subfossil musk ox specimens (originally identified as Bootherium and/or Symbos), allowing us to evaluate the identity of these muskoxen, explore their phylogeography, and estimate the timeline for their evolution. We also used nuclear genomic data to determine the sex of six of our seven samples. Ultimately, our molecular data support the synonymisation of the North American muskoxen Bootherium and Symbos.

Molecular phylogenetics supports the origin of an endemic Balearic shrew lineage (Nesiotites) coincident with the Messinian Salinity Crisis.

The red-toothed shrews (Soricinae) are the most widespread subfamily of shrews, distributed from northern South America to North America and Eurasia. Within this subfamily, the tribe Nectogalini includes the fossil species Nesiotites hidalgo recorded from the Late Pleistocene to Holocene of the Balearic Islands (Western Mediterranean). Although there is a consensus about the close relationship between the extinct red-toothed shrew genera Nesiotites and Asoriculus based on morphology, molecular data are necessary to further evaluate the phylogenetic relationships of the Balearic fossils. We obtained a near complete mitochondrial genome of N. hidalgo, allowing the first molecular phylogenetic analysis of this species. Analyses based on 15,167 bp of the mitochondrial genome placed N. hidalgo as close relative to the extant Himalayan shrew (Soriculus nigrescens), and a combined analysis using molecular and morphological data confirm that N. hidalgo and Asoriculus gibberodon are sister-taxa with S. nigrescens as the immediate outgroup. Molecular clock and divergence estimates suggest that the split between N. hidalgo and its closest living relative occurred around 6.44 Ma, which is in agreement with the previously proposed colonisation of the Balearic Islands from mainland Europe by nectogaline shrews during the Messinian Salinity Crisis (5.97-5.33 My ago). Our results highlight that it is possible to retrieve genetic data from extinct small mammals from marginal environments for DNA preservation. Additional finds from the fossil record of Soricinae from the Eurasian Late Miocene/Early Pliocene are needed to shed further light on the still confusing taxonomy and paleobiogeography of this clade.

Genetic diversity and drivers of dwarfism in extinct island emu populations.

Australia's iconic emu (Dromaius novaehollandiae novaehollandiae) is the only living representative of its genus, but fossil evidence and reports from early European explorers suggest that three island forms (at least two of which were dwarfs) became extinct during the nineteenth century. While one of these-the King Island emu-has been found to be conspecific with Australian mainland emus, little is known about how the other two forms-Kangaroo Island and Tasmanian emus-relate to the others, or even the size of Tasmanian emus. We present a comprehensive genetic and morphological analysis of Dromaius diversity, including data from one of the few definitively genuine Tasmanian emu specimens known. Our genetic analyses suggest that all the island populations represent sub-populations of mainland Dnovaehollandiae Further, the size of island emus and those on the mainland appears to scale linearly with island size but not time since isolation, suggesting that island size-and presumably concomitant limitations on resource availability-may be a more important driver of dwarfism in island emus, though its precise contribution to emu dwarfism remains to be confirmed.

The origin and phylogenetic relationships of the New Zealand ravens.

The relationships of the extinct New Zealand ravens (Corvus spp.) are poorly understood. We sequenced the mitogenomes of the two currently recognised species and found they were sister-taxa to a clade comprising the Australian raven, little raven, and forest raven (C.coronoides, C. mellori and C. tasmanicus respectively). The divergence between the New Zealand ravens and Australian raven clade occurred in the latest Pliocene, which coincides with the onset of glacial deforestation. We also found that the divergence between the two putative New Zealand species C. antipodum and C. moriorum probably occurred in the late Pleistocene making their separation as species untenable. Consequently, we consider Corax antipodum (Forbes, 1893) to be a subspecies of Corvus moriorum Forbes, 1892. We re-examine the osteological evidence that led 19th century researchers to assign the New Zealand taxa to a separate genus, and re-assess these features in light of our new phylogenetic hypotheses. Like previous researchers, we conclude that the morphology of the palate of C. moriorum is unique among the genus Corvus, and suggest this may be an adaptation for a specialist diet.

Ancient DNA analysis of the extinct North American flat-headed peccary (Platygonus compressus).

The geographical range of extant peccaries extends from the southwestern United States through Central America and into northern Argentina. However, from the Miocene until the Pleistocene now-extinct peccary species inhabited the entirety of North America. Relationships among the living and extinct species have long been contentious. Similarly, how and when peccaries moved from North to South America is unclear. The North American flat-headed peccary (Platygonus compressus) became extinct at the end of the Pleistocene and is one of the most abundant subfossil taxa found in North America, yet despite this extensive fossil record its phylogenetic position has not been resolved. This study is the first to present DNA data from the flat-headed peccary and full mitochondrial genome sequences of all the extant peccary species. We performed a molecular phylogenetic analysis to determine the relationships among ancient and extant peccary species. Our results suggested that the flat-headed peccary is sister-taxon to a clade comprising the extant peccary species. Divergence date estimates from our molecular dating analyses suggest that if extant peccary diversification occurred in South America then their common ancestor must have dispersed from North America to South America well before the establishment of the Isthmus of Panama. We also investigated the genetic diversity of the flat-headed peccary by performing a preliminary population study on specimens from Sheriden Cave, Ohio. Flat-headed peccaries from Sheriden Cave appear to be genetically diverse and show no signature of population decline prior to extinction. Including additional extinct Pleistocene peccary species in future phylogenetic analyses will further clarify peccary evolution.

Late pleistocene Australian marsupial DNA clarifies the affinities of extinct megafaunal kangaroos and wallabies.

Understanding the evolution of Australia's extinct marsupial megafauna has been hindered by a relatively incomplete fossil record and convergent or highly specialized morphology, which confound phylogenetic analyses. Further, the harsh Australian climate and early date of most megafaunal extinctions (39-52 ka) means that the vast majority of fossil remains are unsuitable for ancient DNA analyses. Here, we apply cross-species DNA capture to fossils from relatively high latitude, high altitude caves in Tasmania. Using low-stringency hybridization and high-throughput sequencing, we were able to retrieve mitochondrial sequences from two extinct megafaunal macropodid species. The two specimens, Simosthenurus occidentalis (giant short-faced kangaroo) and Protemnodon anak (giant wallaby), have been radiocarbon dated to 46-50 and 40-45 ka, respectively. This is significantly older than any Australian fossil that has previously yielded DNA sequence information. Processing the raw sequence data from these samples posed a bioinformatic challenge due to the poor preservation of DNA. We explored several approaches in order to maximize the signal-to-noise ratio in retained sequencing reads. Our findings demonstrate the critical importance of adopting stringent processing criteria when distant outgroups are used as references for mapping highly fragmented DNA. Based on the most stringent nucleotide data sets (879 bp for S. occidentalis and 2,383 bp for P. anak), total-evidence phylogenetic analyses confirm that macropodids consist of three primary lineages: Sthenurines such as Simosthenurus (extinct short-faced kangaroos), the macropodines (all other wallabies and kangaroos), and the enigmatic living banded hare-wallaby Lagostrophus fasciatus (Lagostrophinae). Protemnodon emerges as a close relative of Macropus (large living kangaroos), a position not supported by recent morphological phylogenetic analyses.

Ancient DNA from the extinct South American giant glyptodont Doedicurus sp. (Xenarthra: Glyptodontidae) reveals that glyptodonts evolved from Eocene armadillos.

Glyptodonts were giant (some of them up to ~2400 kg), heavily armoured relatives of living armadillos, which became extinct during the Late Pleistocene/early Holocene alongside much of the South American megafauna. Although glyptodonts were an important component of Cenozoic South American faunas, their early evolution and phylogenetic affinities within the order Cingulata (armoured New World placental mammals) remain controversial. In this study, we used hybridization enrichment and high-throughput sequencing to obtain a partial mitochondrial genome from Doedicurus sp., the largest (1.5 m tall, and 4 m long) and one of the last surviving glyptodonts. Our molecular phylogenetic analyses revealed that glyptodonts fall within the diversity of living armadillos. Reanalysis of morphological data using a molecular 'backbone constraint' revealed several morphological characters that supported a close relationship between glyptodonts and the tiny extant fairy armadillos (Chlamyphorinae). This is surprising as these taxa are among the most derived cingulates: glyptodonts were generally large-bodied and heavily armoured, while the fairy armadillos are tiny (~9-17 cm) and adapted for burrowing. Calibration of our phylogeny with the first appearance of glyptodonts in the Eocene resulted in a more precise timeline for xenarthran evolution. The osteological novelties of glyptodonts and their specialization for grazing appear to have evolved rapidly during the Late Eocene to Early Miocene, coincident with global temperature decreases and a shift from wet closed forest towards drier open woodland and grassland across much of South America. This environmental change may have driven the evolution of glyptodonts, culminating in the bizarre giant forms of the Pleistocene.

Ancient mitochondrial genomes clarify the evolutionary history of New Zealand's enigmatic acanthisittid wrens.

The New Zealand acanthisittid wrens are the sister-taxon to all other "perching birds" (Passeriformes) and - including recently extinct species - represent the most diverse endemic passerine family in New Zealand. Consequently, they are important for understanding both the early evolution of Passeriformes and the New Zealand biota. However, five of the seven species have become extinct since the arrival of humans in New Zealand, complicating evolutionary analyses. The results of morphological analyses have been largely equivocal, and no comprehensive genetic analysis of Acanthisittidae has been undertaken. We present novel mitochondrial genome sequences from four acanthisittid species (three extinct, one extant), allowing us to resolve the phylogeny and revise the taxonomy of acanthisittids. Reanalysis of morphological data in light of our genetic results confirms a close relationship between the extant rifleman (Acanthisitta chloris) and an extinct Miocene wren (Kuiornis indicator), making Kuiornis a useful calibration point for molecular dating of passerines. Our molecular dating analyses reveal that the stout-legged wrens (Pachyplichas) diverged relatively recently from a more gracile (Xenicus-like) ancestor. Further, our results suggest a possible Early Oligocene origin of the basal Lyall's wren (Traversia) lineage, which would imply that Acanthisittidae survived the Oligocene marine inundation of New Zealand and therefore that the inundation was not complete.

Ancient mitochondrial DNA reveals convergent evolution of giant short-faced bears (Tremarctinae) in North and South America.

The Tremarctinae are a subfamily of bears endemic to the New World, including two of the largest terrestrial mammalian carnivores that have ever lived: the giant, short-faced bears Arctodus simus from North America and Arctotherium angustidens from South America (greater than or equal to 1000 kg). Arctotherium angustidens became extinct during the Early Pleistocene, whereas Arctodus simus went extinct at the very end of the Pleistocene. The only living tremarctine is the spectacled bear (Tremarctos ornatus), a largely herbivorous bear that is today only found in South America. The relationships among the spectacled bears (Tremarctos), South American short-faced bears (Arctotherium) and North American short-faced bears (Arctodus) remain uncertain. In this study, we sequenced a mitochondrial genome from an Arctotherium femur preserved in a Chilean cave. Our molecular phylogenetic analyses revealed that the South American short-faced bears were more closely related to the extant South American spectacled bear than to the North American short-faced bears. This result suggests striking convergent evolution of giant forms in the two groups of short-faced bears (Arctodus and Arctotherium), potentially as an adaptation to dominate competition for megafaunal carcasses.