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.

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.

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.

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.

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.

Estimating the dwarfing rate of an extinct Sicilian elephant.

Evolution on islands, together with the often extreme phenotypic changes associated with it, has attracted much interest from evolutionary biologists. However, measuring the rate of change of phenotypic traits of extinct animals can be challenging, in part due to the incompleteness of the fossil record. Here, we use combined molecular and fossil evidence to define the minimum and maximum rate of dwarfing in an extinct Mediterranean dwarf elephant from Puntali Cave (Sicily).1 Despite the challenges associated with recovering ancient DNA from warm climates,2 we successfully retrieved a mitogenome from a sample with an estimated age between 175,500 and 50,000 years. Our results suggest that this specific Sicilian elephant lineage evolved from one of the largest terrestrial mammals that ever lived3 to an island species weighing less than 20% of its original mass with an estimated mass reduction between 0.74 and 200.95 kg and height reduction between 0.15 and 41.49 mm per generation. We show that combining ancient DNA with paleontological and geochronological evidence can constrain the timing of phenotypic changes with greater accuracy than could be achieved using any source of evidence in isolation.

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.

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.

Pre-extinction Demographic Stability and Genomic Signatures of Adaptation in the Woolly Rhinoceros.

Ancient DNA has significantly improved our understanding of the evolution and population history of extinct megafauna. However, few studies have used complete ancient genomes to examine species responses to climate change prior to extinction. The woolly rhinoceros (Coelodonta antiquitatis) was a cold-adapted megaherbivore widely distributed across northern Eurasia during the Late Pleistocene and became extinct approximately 14 thousand years before present (ka BP). While humans and climate change have been proposed as potential causes of extinction [1-3], knowledge is limited on how the woolly rhinoceros was impacted by human arrival and climatic fluctuations [2]. Here, we use one complete nuclear genome and 14 mitogenomes to investigate the demographic history of woolly rhinoceros leading up to its extinction. Unlike other northern megafauna, the effective population size of woolly rhinoceros likely increased at 29.7 ka BP and subsequently remained stable until close to the species' extinction. Analysis of the nuclear genome from a ∼18.5-ka-old specimen did not indicate any increased inbreeding or reduced genetic diversity, suggesting that the population size remained steady for more than 13 ka following the arrival of humans [4]. The population contraction leading to extinction of the woolly rhinoceros may have thus been sudden and mostly driven by rapid warming in the Bølling-Allerød interstadial. Furthermore, we identify woolly rhinoceros-specific adaptations to arctic climate, similar to those of the woolly mammoth. This study highlights how species respond differently to climatic fluctuations and further illustrates the potential of palaeogenomics to study the evolutionary history of extinct species.

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.

Evolution and extinction of the giant rhinoceros Elasmotherium sibiricum sheds light on late Quaternary megafaunal extinctions.

Understanding extinction events requires an unbiased record of the chronology and ecology of victims and survivors. The rhinoceros Elasmotherium sibiricum, known as the 'Siberian unicorn', was believed to have gone extinct around 200,000 years ago-well before the late Quaternary megafaunal extinction event. However, no absolute dating, genetic analysis or quantitative ecological assessment of this species has been undertaken. Here, we show, by accelerator mass spectrometry radiocarbon dating of 23 individuals, including cross-validation by compound-specific analysis, that E. sibiricum survived in Eastern Europe and Central Asia until at least 39,000 years ago, corroborating a wave of megafaunal turnover before the Last Glacial Maximum in Eurasia, in addition to the better-known late-glacial event. Stable isotope data indicate a dry steppe niche for E. sibiricum and, together with morphology, a highly specialized diet that probably contributed to its extinction. We further demonstrate, with DNA sequencing data, a very deep phylogenetic split between the subfamilies Elasmotheriinae and Rhinocerotinae that includes all the living rhinoceroses, settling a debate based on fossil evidence and confirming that the two lineages had diverged by the Eocene. As the last surviving member of the Elasmotheriinae, the demise of the 'Siberian unicorn' marked the extinction of this subfamily.

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.

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.

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.

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.

Millennial climatic fluctuations are key to the structure of last glacial ecosystems.

Whereas fossil evidence indicates extensive treeless vegetation and diverse grazing megafauna in Europe and northern Asia during the last glacial, experiments combining vegetation models and climate models have to-date simulated widespread persistence of trees. Resolving this conflict is key to understanding both last glacial ecosystems and extinction of most of the mega-herbivores. Using a dynamic vegetation model (DVM) we explored the implications of the differing climatic conditions generated by a general circulation model (GCM) in "normal" and "hosing" experiments. Whilst the former approximate interstadial conditions, the latter, designed to mimic Heinrich Events, approximate stadial conditions. The "hosing" experiments gave simulated European vegetation much closer in composition to that inferred from fossil evidence than did the "normal" experiments. Given the short duration of interstadials, and the rate at which forest cover expanded during the late-glacial and early Holocene, our results demonstrate the importance of millennial variability in determining the character of last glacial ecosystems.

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.

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.