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.

Million-year-old DNA sheds light on the genomic history of mammoths.

Temporal genomic data hold great potential for studying evolutionary processes such as speciation. However, sampling across speciation events would, in many cases, require genomic time series that stretch well back into the Early Pleistocene subepoch. Although theoretical models suggest that DNA should survive on this timescale1, the oldest genomic data recovered so far are from a horse specimen dated to 780-560 thousand years ago2. Here we report the recovery of genome-wide data from three mammoth specimens dating to the Early and Middle Pleistocene subepochs, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these lineages gave rise to the woolly mammoth and the other represents a previously unrecognized lineage that was ancestral to the first mammoths to colonize North America. Our analyses reveal that the Columbian mammoth of North America traces its ancestry to a Middle Pleistocene hybridization between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were already present one million years ago. These findings highlight the potential of deep-time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes.

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Head to head: the case for fighting behaviour in Megaloceros giganteus using finite-element analysis.

The largest antlers of any known deer species belonged to the extinct giant deer Megaloceros giganteus. It has been argued that their antlers were too large for use in fighting, instead being used only in ritualized displays to attract mates. Here, we used finite-element analysis to test whether the antlers of M. giganteus could have withstood forces generated during fighting. We compared the mechanical performance of antlers in M. giganteus with three extant deer species: red deer (Cervus elaphus), fallow deer (Dama dama) and elk (Alces alces). Von Mises stress results suggest that M. giganteus was capable of withstanding some fighting loads, provided that their antlers interlocked proximally, and that their antlers were best adapted for withstanding loads from twisting rather than pushing actions, as are other deer with palmate antlers. We conclude that fighting in M. giganteus was probably more constrained and predictable than in extant deer.

The role of behaviour in adaptive morphological evolution of African proboscideans.

The fossil record richly illustrates the origin of morphological adaptation through time. However, our understanding of the selective forces responsible in a given case, and the role of behaviour in the process, is hindered by assumptions of synchrony between environmental change, behavioural innovation and morphological response. Here I show, from independent proxy data through a 20-million-year sequence of fossil proboscideans in East Africa, that changes in environment, diet and morphology are often significantly offset chronologically, allowing dissection of the roles of behaviour and different selective drivers. These findings point the way to hypothesis-driven testing of the interplay between habitat change, behaviour and morphological adaptation with the use of independent proxies in the fossil record.

Plant controls on Late Quaternary whole ecosystem structure and function.

Plants and animals influence biomass production and nutrient cycling in terrestrial ecosystems; however, their relative importance remains unclear. We assessed the extent to which mega-herbivore species controlled plant community composition and nutrient cycling, relative to other factors during and after the Late Quaternary extinction event in Britain and Ireland, when two-thirds of the region's mega-herbivore species went extinct. Warmer temperatures, plant-soil and plant-plant interactions, and reduced burning contributed to the expansion of woody plants and declining nitrogen availability in our five study ecosystems. Shrub biomass was consistently one of the strongest predictors of ecosystem change, equalling or exceeding the effects of other biotic and abiotic factors. In contrast, there was relatively little evidence for mega-herbivore control on plant community composition and nitrogen availability. The ability of plants to determine the fate of terrestrial ecosystems during periods of global environmental change may therefore be greater than previously thought.

Growth in fossil and extant deer and implications for body size and life history evolution.

BACKGROUND: Body size variation within clades of mammals is widespread, but the developmental and life-history mechanisms by which this variation is achieved are poorly understood, especially in extinct forms. An illustrative case study is that of the dwarfed morphotypes of Candiacervus from the Pleistocene of Crete versus the giant deer Megaloceros giganteus, both in a clade together with Dama dama among extant species. Histological analyses of long bones and teeth in a phylogenetic context have been shown to provide reliable estimates of growth and life history patterns in extant and extinct mammals. RESULTS: Similarity of bone tissue types across the eight species examined indicates a comparable mode of growth in deer, with long bones mainly possessing primary plexiform fibrolamellar bone. Low absolute growth rates characterize dwarf Candiacervus sp. II and C. ropalophorus compared to Megaloceros giganteus displaying high rates, whereas Dama dama is characterized by intermediate to low growth rates. The lowest recorded rates are those of the Miocene small stem cervid Procervulus praelucidus. Skeletal maturity estimates indicate late attainment in sampled Candiacervus and Procervulus praelucidus. Tooth cementum analysis of first molars of two senile Megaloceros giganteus specimens revealed ages of 16 and 19 years whereas two old dwarf Candiacervus specimens gave ages of 12 and 18 years. CONCLUSIONS: There is a rich histological record of growth across deer species recorded in long bones and teeth, which can be used to understand ontogenetic patterns within species and phylogenetic ones across species. Growth rates sensu Sander & Tückmantel plotted against the anteroposterior bone diameter as a proxy for body mass indicate three groups: one with high growth rates including Megaloceros, Cervus, Alces, and Dama; an intermediate group with Capreolus and Muntiacus; and a group showing low growth rates, including dwarf Candiacervus and Procervulus. Dwarf Candiacervus, in an allometric context, show an extended lifespan compared to other deer of similar body size such as Mazama which has a maximum longevity of 12 years in the wild. Comparison with other clades of mammals reveals that changes in size and life history in evolution have occurred in parallel, with various modes of skeletal tissue modification.

Insular dwarfism in hippos and a model for brain size reduction in Homo floresiensis.

Body size reduction in mammals is usually associated with only moderate brain size reduction, because the brain and sensory organs complete their growth before the rest of the body during ontogeny. On this basis, 'phyletic dwarfs' are predicted to have a greater relative brain size than 'phyletic giants'. However, this trend has been questioned in the special case of dwarfism of mammals on islands. Here we show that the endocranial capacities of extinct dwarf species of hippopotamus from Madagascar are up to 30% smaller than those of a mainland African ancestor scaled to equivalent body mass. These results show that brain size reduction is much greater than predicted from an intraspecific 'late ontogenetic' model of dwarfism in which brain size scales to body size with an exponent of 0.35. The nature of the proportional change or grade shift observed here indicates that selective pressures on brain size are potentially independent of those on body size. This study demonstrates empirically that it is mechanistically possible for dwarf mammals on islands to evolve significantly smaller brains than would be predicted from a model of dwarfing based on the intraspecific scaling of the mainland ancestor. Our findings challenge current understanding of brain-body allometric relationships in mammals and suggest that the process of dwarfism could in principle explain small brain size, a factor relevant to the interpretation of the small-brained hominin found on the Island of Flores, Indonesia.

Serial population extinctions in a small mammal indicate Late Pleistocene ecosystem instability.

The Late Pleistocene global extinction of many terrestrial mammal species has been a subject of intensive scientific study for over a century, yet the relative contributions of environmental changes and the global expansion of humans remain unresolved. A defining component of these extinctions is a bias toward large species, with the majority of small-mammal taxa apparently surviving into the present. Here, we investigate the population-level history of a key tundra-specialist small mammal, the collared lemming (Dicrostonyx torquatus), to explore whether events during the Late Pleistocene had a discernible effect beyond the large mammal fauna. Using ancient DNA techniques to sample across three sites in North-West Europe, we observe a dramatic reduction in genetic diversity in this species over the last 50,000 y. We further identify a series of extinction-recolonization events, indicating a previously unrecognized instability in Late Pleistocene small-mammal populations, which we link with climatic fluctuations. Our results reveal climate-associated, repeated regional extinctions in a keystone prey species across the Late Pleistocene, a pattern likely to have had an impact on the wider steppe-tundra community, and one that is concordant with environmental change as a major force in structuring Late Pleistocene biodiversity.

Faunal record identifies Bering isthmus conditions as constraint to end-Pleistocene migration to the New World.

Human colonization of the New World is generally believed to have entailed migrations from Siberia across the Bering isthmus. However, the limited archaeological record of these migrations means that details of the timing, cause and rate remain cryptic. Here, we have used a combination of ancient DNA, 14C dating, hydrogen and oxygen isotopes, and collagen sequencing to explore the colonization history of one of the few other large mammals to have successfully migrated into the Americas at this time: the North American elk (Cervus elaphus canadensis), also known as wapiti. We identify a long-term occupation of northeast Siberia, far beyond the species's current Old World distribution. Migration into North America occurred at the end of the last glaciation, while the northeast Siberian source population became extinct only within the last 500 years. This finding is congruent with a similar proposed delay in human colonization, inferred from modern human mitochondrial DNA, and suggestions that the Bering isthmus was not traversable during parts of the Late Pleistocene. Our data imply a fundamental constraint in crossing Beringia, placing limits on the age and mode of human settlement in the Americas, and further establish the utility of ancient DNA in palaeontological investigations of species histories.

Discrimination of ivory from extant and extinct elephant species using Raman spectroscopy: A potential non-destructive technique for combating illegal wildlife trade.

The use of elephant ivory as a commodity is a factor in declining elephant populations. Despite recent worldwide elephant ivory trade bans, mammoth ivory trade remains unregulated. This complicates law enforcement efforts, as distinguishing between ivory from extant and extinct species requires costly, destructive and time consuming methods. Elephant and mammoth ivory mainly consists of dentine, a mineralized connective tissue that contains an organic collagenous component and an inorganic component of calcium phosphate minerals, similar in structure to hydroxyapatite crystals. Raman spectroscopy is a non-invasive laser-based technique that has previously been used for the study of bone and mineral chemistry. Ivory and bone have similar biochemical properties, making Raman spectroscopy a promising method for species identification based on ivory. This study aimed to test the hypothesis that it is possible to identify differences in the chemistry of mammoth and elephant ivory using Raman spectroscopy. Mammoth and elephant tusks were obtained from the Natural History Museum in London, UK. Included in this study were eight samples of ivory from Mammuthus primigenius, two samples of carved ivory bangles from Africa (Loxodonta species), and one cross section of a tusk from Elephas maximus. The ivory was scanned using an inVia Raman micro spectrometer equipped with a x50 objective lens and a 785nm laser. Spectra were acquired using line maps and individual spectral points were acquired randomly or at points of interest on all samples. The data was then analysed using principal component analysis (PCA) with use of an in-house MATLAB script. Univariate analysis of peak intensity ratios of phosphate to amide I and III peaks, and carbonate to phosphate peaks showed statistical differences (p<0.0001) in the average peak intensity ratios between Mammuthus primigenius, Loxodonta spp. and Elephas maximus. Full width at half maximum hight (FWHM)analysis of the phosphate peak demonstrated higher crystal maturity of Mammuthus primigenius compared to living elephant species. The results of the study have established that spectra acquired by Raman spectroscopy can be separated into distinct classes through PCA. In conclusion, this study has shown that well-preserved mammoth and elephant ivory has the potential to be characterized using Raman spectroscopy, providing a promising method for species identification. The results of this study will be valuable in developing quick and non-destructive methods for the identification of ivory, which will have direct applications in archaeology and the regulation of international trade.

Holarctic genetic structure and range dynamics in the woolly mammoth.

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32-34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20-15 kyr BP, when a severe population size decline occurred.

Late-glacial recolonization and phylogeography of European red deer (Cervus elaphus L.).

The Pleistocene was an epoch of extreme climatic and environmental changes. How individual species responded to the repeated cycles of warm and cold stages is a major topic of debate. For the European fauna and flora, an expansion-contraction model has been suggested, whereby temperate species were restricted to southern refugia during glacial times and expanded northwards during interglacials, including the present interglacial (Holocene). Here, we test this model on the red deer (Cervus elaphus) a large and highly mobile herbivore, using both modern and ancient mitochondrial DNA from the entire European range of the species over the last c. 40,000 years. Our results indicate that this species was sensitive to the effects of climate change. Prior to the Last Glacial Maximum (LGM) haplogroups restricted today to South-East Europe and Western Asia reached as far west as the UK. During the LGM, red deer was mainly restricted to southern refugia, in Iberia, the Balkans and possibly in Italy and South-Western Asia. At the end of the LGM, red deer expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, Scandinavia, Germany, Poland and Belarus. Ancient DNA data cannot rule out refugial survival of red deer in North-West Europe through the LGM. Had such deer survived, though, they were replaced by deer migrating from Iberia at the end of the glacial. The Balkans served as a separate LGM refugium and were probably connected to Western Asia with genetic exchange between the two areas.

Fluctuating climate and dietary innovation drove ratcheted evolution of proboscidean dental traits.

Identification of the selective forces that shaped adaptive phenotypes generally relies on current habitat and function, but these may differ from the context in which adaptations arose. Moreover, the fixation of adaptive change in a fluctuating environment and the mechanisms of long-term trends are still poorly understood, as is the role of behaviour in triggering these processes. Time series of fossils can provide evidence on these questions, but examples of individual lineages with adequate fossil and proxy data over extended periods are rare. Here, we present new data on proboscidean dental evolution in East Africa over the past 26 million years, tracking temporal patterns of morphological change in relation to proxy evidence of diet, vegetation and climate (aridity). We show that behavioural experimentation in diet is correlated with environmental context, and that major adaptive change in dental traits followed the changes in diet and environment but only after acquisition of functional innovations in the masticatory system. We partition traits by selective agent, showing that the acquisition of high, multiridged molars was primarily a response to an increase in open, arid environments with high dust accumulation, whereas enamel folding was more associated with the amount of grass in the diet. We further show that long-term trends in these features proceeded in a ratchet-like mode, alternating between directional change at times of high selective pressure and stasis when the selective regime reversed. This provides an explanation for morphology adapted to more extreme conditions than current usage (Liem's Paradox). Our study illustrates how, in fossil series with adequate stratigraphic control and proxy data, environmental and behavioural factors can be mapped on to time series of morphological change, illuminating the mode of acquisition of an adaptive complex.

Relationships of Late Pleistocene giant deer as revealed by Sinomegaceros mitogenomes from East Asia.

The giant deer, widespread in northern Eurasia during the Late Pleistocene, have been classified as western Megaloceros and eastern Sinomegaceros through morphological studies. While Megaloceros's evolutionary history has been unveiled through mitogenomes, Sinomegaceros remains molecularly unexplored. Herein, we generated mitogenomes of giant deer from East Asia. We find that, in contrast to the morphological differences between Megaloceros and Sinomegaceros, they are mixed in the mitochondrial phylogeny, and Siberian specimens suggest a range contact or overlap between these two groups. Meanwhile, one deep divergent clade and another surviving until 20.1 thousand years ago (ka) were detected in northeastern China, the latter implying this area as a potential refugium during the Last Glacial Maximum (LGM). Moreover, stable isotope analyses indicate correlations between climate-introduced vegetation changes and giant deer extinction. Our study demonstrates the genetic relationship between eastern and western giant deer and explores the promoters of their extirpation in northern East Asia.

Genomics of adaptive evolution in the woolly mammoth.

Ancient genomes provide a tool to investigate the genetic basis of adaptations in extinct organisms. However, the identification of species-specific fixed genetic variants requires the analysis of genomes from multiple individuals. Moreover, the long-term scale of adaptive evolution coupled with the short-term nature of traditional time series data has made it difficult to assess when different adaptations evolved. Here, we analyze 23 woolly mammoth genomes, including one of the oldest known specimens at 700,000 years old, to identify fixed derived non-synonymous mutations unique to the species and to obtain estimates of when these mutations evolved. We find that at the time of its origin, the woolly mammoth had already acquired a broad spectrum of positively selected genes, including ones associated with hair and skin development, fat storage and metabolism, and immune system function. Our results also suggest that these phenotypes continued to evolve during the last 700,000 years, but through positive selection on different sets of genes. Finally, we also identify additional genes that underwent comparatively recent positive selection, including multiple genes related to skeletal morphology and body size, as well as one gene that may have contributed to the small ear size in Late Quaternary woolly mammoths.

Extreme insular dwarfism evolved in a mammoth.

The insular dwarfism seen in Pleistocene elephants has come to epitomize the island rule; yet our understanding of this phenomenon is hampered by poor taxonomy. For Mediterranean dwarf elephants, where the most extreme cases of insular dwarfism are observed, a key systematic question remains unresolved: are all taxa phyletic dwarfs of a single mainland species Palaeoloxodon antiquus (straight-tusked elephant), or are some referable to Mammuthus (mammoths)? Ancient DNA and geochronological evidence have been used to support a Mammuthus origin for the Cretan 'Palaeoloxodon' creticus, but these studies have been shown to be flawed. On the basis of existing collections and recent field discoveries, we present new, morphological evidence for the taxonomic status of 'P'. creticus, and show that it is indeed a mammoth, most probably derived from Early Pleistocene Mammuthus meridionalis or possibly Late Pliocene Mammuthus rumanus. We also show that Mammuthus creticus is smaller than other known insular dwarf mammoths, and is similar in size to the smallest dwarf Palaeoloxodon species from Sicily and Malta, making it the smallest mammoth species known to have existed. These findings indicate that extreme insular dwarfism has evolved to a similar degree independently in two elephant lineages.

Microsatellite genotyping reveals end-Pleistocene decline in mammoth autosomal genetic variation.

The last glaciation was a dynamic period with strong impact on the demography of many species and populations. In recent years, mitochondrial DNA sequences retrieved from radiocarbon-dated remains have provided novel insights into the history of Late Pleistocene populations. However, genotyping of loci from the nuclear genome may provide enhanced resolution of population-level changes. Here, we use four autosomal microsatellite DNA markers to investigate the demographic history of woolly mammoths (Mammuthus primigenius) in north-eastern Siberia from before 60 000 years ago up until the species' final disappearance c.4000 years ago. We identified two genetic groups, implying a marked temporal genetic differentiation between samples with radiocarbon ages older than 12 thousand radiocarbon years before present (ka) and those younger than 9ka. Simulation-based analysis indicates that this dramatic change in genetic composition, which included a decrease in individual heterozygosity of approximately 30%, was due to a multifold reduction in effective population size. A corresponding reduction in genetic variation was also detected in the mitochondrial DNA, where about 65% of the diversity was lost. We observed no further loss in genetic variation during the Holocene, which suggests a rapid final extinction event.

Multiplex amplification of the mammoth mitochondrial genome and the evolution of Elephantidae.

In studying the genomes of extinct species, two principal limitations are typically the small quantities of endogenous ancient DNA and its degraded condition, even though products of up to 1,600 base pairs (bp) have been amplified in rare cases. Using small overlapping polymerase chain reaction products, longer stretches of sequences or even whole mitochondrial genomes can be reconstructed, but this approach is limited by the number of amplifications that can be performed from rare samples. Thus, even from well-studied Pleistocene species such as mammoths, ground sloths and cave bears, no DNA sequences of more than about 1,000 bp have been reconstructed. Here we report the complete mitochondrial genome sequence of the Pleistocene woolly mammoth Mammuthus primigenius. We used about 200 mg of bone and a new approach that allows the simultaneous retrieval of multiple sequences from small amounts of degraded DNA. Our phylogenetic analyses show that the mammoth was more closely related to the Asian than to the African elephant. However, the divergence of mammoth, African and Asian elephants occurred over a short time, corresponding to only about 7% of the total length of the phylogenetic tree for the three evolutionary lineages.

The earliest record of human activity in northern Europe.

The colonization of Eurasia by early humans is a key event after their spread out of Africa, but the nature, timing and ecological context of the earliest human occupation of northwest Europe is uncertain and has been the subject of intense debate. The southern Caucasus was occupied about 1.8 million years (Myr) ago, whereas human remains from Atapuerca-TD6, Spain (more than 780 kyr ago) and Ceprano, Italy (about 800 kyr ago) show that early Homo had dispersed to the Mediterranean hinterland before the Brunhes-Matuyama magnetic polarity reversal (780 kyr ago). Until now, the earliest uncontested artefacts from northern Europe were much younger, suggesting that humans were unable to colonize northern latitudes until about 500 kyr ago. Here we report flint artefacts from the Cromer Forest-bed Formation at Pakefield (52 degrees N), Suffolk, UK, from an interglacial sequence yielding a diverse range of plant and animal fossils. Event and lithostratigraphy, palaeomagnetism, amino acid geochronology and biostratigraphy indicate that the artefacts date to the early part of the Brunhes Chron (about 700 kyr ago) and thus represent the earliest unequivocal evidence for human presence north of the Alps.