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

Worldwide Late Pleistocene and Early Holocene population declines in extant megafauna are associated with Homo sapiens expansion rather than climate change

The worldwide extinction of megafauna during the Late Pleistocene and Early Holocene is evident from the fossil record, with dominant theories suggesting a climate, human or combined impact cause. Consequently, two disparate scenarios are possible for the surviving megafauna during this time period - they could have declined due to similar pressures, or increased in population size due to reductions in competition or other biotic pressures. We therefore infer population histories of 139 extant megafauna species using genomic data which reveal population declines in 91% of species throughout the Quaternary period, with larger species experiencing the strongest decreases. Declines become ubiquitous 32-76 kya across all landmasses, a pattern better explained by worldwide Homo sapiens expansion than by changes in climate. We estimate that, in consequence, total megafauna abundance, biomass, and energy turnover decreased by 92-95% over the past 50,000 years, implying major human-driven ecosystem restructuring at a global scale.

Ancient DNA of narrow-headed vole reveal common features of the Late Pleistocene population dynamics in cold-adapted small mammals

The narrow-headed vole, collared lemming and common vole were the most abundant small mammal species across the Eurasian Late Pleistocene steppe-tundra environment. Previous ancient DNA studies of the collared lemming and common vole have revealed dynamic population histories shaped by climatic fluctuations. To investigate the extent to which species with similar adaptations share common evolutionary histories, we generated a dataset comprised the mitochondrial genomes of 139 ancient and 6 modern narrow-headed voles from several sites across Europe and northwestern Asia covering approximately the last 100 thousand years (kyr). We inferred Bayesian time-aware phylogenies using 11 radiocarbon-dated samples to calibrate the molecular clock. Divergence of the main mtDNA lineages across the three species occurred during marine isotope stages (MIS) 7 and MIS 5, suggesting a common response of species adapted to open habitat during interglacials. We identified several time-structured mtDNA lineages in European narrow-headed vole, suggesting lineage turnover. The timing of some of these turnovers was synchronous across the three species, allowing us to identify the main drivers of the Late Pleistocene dynamics of steppe- and cold-adapted species.

Mammal extinction facilitated biome shift and human population change during the last glacial termination in East-Central Europe

The study of local extinction times, together with the associated environmental and human population changes in the last glacial termination, provides insights into the causes of mega- and microfauna extinctions. In East-Central (EC) Europe, groups of Palaeolithic humans were present throughout the last glacial maximum, but disappeared suddenly around 15,200 cal BP. In this study cave sediment profiles dated using radiocarbon techniques and a large set of mammal bones dated directly by AMS ¹⁴C were used to determine local extinction times. These were, in turn, compared to changes in the total megafauna population of EC Europe derived from coprophilous fungi, the Epigravettian population decline, quantitative climate models, pollen and plant macrofossil inferred climate, as well as to biome reconstructions. The results suggest that the population size of large herbivores decreased in the area after 17,700 cal BP, when temperate tree abundance and warm continental steppe cover both increased in the lowlands. Boreal forest expansion started around 16,200 cal BP. Cave sediments show the decline of narrow-headed vole and arctic lemming populations specifically associated with a tundra environment at the same time and the expansion of the common vole, an inhabitant of steppes. The last dated appearance of arctic lemming was at ~ 16,640 cal BP, while that of the narrow-headed vole at ~ 13,340, and the estimated extinction time of woolly mammoth was either at 13,830 (GRIWM) or 15,210 (PHASE), and reindeer at 11,860 (GRIWM) or 12,550 cal BP (PHASE). The population decline of the large herbivore fauna slightly preceded changes in terrestrial vegetation, and likely facilitated it via a reduction in the intensity of grazing and the concomitant accumulation of plant biomass. Furthermore, it is possible to conclude that the Late Epigravettian population had high degree of quarry-fidelity; they left the basin when these mammals vanished.

Museum specimens tell the history of rhinoceroses

Rhinoceroses are facing an uncertain future due to human-induced population-size reductions. In this issue of Cell, Liu et al. analyze rhinoceros genomes to study their past and find support for an early divergence between African and Asian species, signatures of past hybridization, and evidence of loss of genetic diversity.

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.

A new rhinoceros clade from the Pleistocene of Asia sheds light on mammal dispersals to the Philippines

Rhinoceroses are among the most endangered mammalian species today. Their past diversity is well documented from the Eocene onward, although their evolutionary history is far from being fully understood. Here, we elucidate the systematic affinities of a Pleistocene rhinoceros species represented by a partial skeleton from 709 ± 68 kya archaeological deposits in Luzon Island, Philippines. We perform a comprehensive phylogenetic analysis, including all living species and a wide array of extinct rhinocerotid species. We confirm the early split between Elasmotheriinae and Rhinocerotinae at c. 35.5 Mya and constrain the divergence between recent Asian and African rhinoceroses at c. 24 Mya, with contrasting phenotypic evolutionary rates in Diceroti and Rhinoceroti. Dental features reveal the existence of an unsuspected Asian Pleistocene clade, referred to as Nesorhinus gen. nov.. It includes the rhinoceros from the Philippines and another extinct species from Taiwan, N. hayasakai. Nesorhinus is the sister-group to a cluster comprising Dicerorhinus and Rhinoceros. Our phylogenetic results strongly suggest an island-hopping dispersal for Nesorhinus, from the Asian mainland towards Luzon via Taiwan by the Late Miocene or later, and Pleistocene dispersals for representatives of Rhinoceros. Nesorhinus philippinensis would be the first perissodactyl species supporting the island-rule hypothesis, with decreased body weight and limb-bone robustness.

Reevaluating the timing of Neanderthal disappearance in Northwest Europe

Elucidating when Neanderthal populations disappeared from Eurasia is a key question in paleoanthropology, and Belgium is one of the key regions for studying the Middle to Upper Paleolithic transition. Previous radiocarbon dating placed the Spy Neanderthals among the latest surviving Neanderthals in Northwest Europe with reported dates as young as 23,880 ± 240 B.P. (OxA-8912). Questions were raised, however, regarding the reliability of these dates. Soil contamination and carbon-based conservation products are known to cause problems during the radiocarbon dating of bulk collagen samples. Employing a compound-specific approach that is today the most efficient in removing contamination and ancient genomic analysis, we demonstrate here that previous dates produced on Neanderthal specimens from Spy were inaccurately young by up to 10,000 y due to the presence of unremoved contamination. Our compound-specific radiocarbon dates on the Neanderthals from Spy and those from Engis and Fonds-de-Forêt demonstrate that they disappeared from Northwest Europe at 44,200 to 40,600 cal B.P. (at 95.4% probability), much earlier than previously suggested. Our data contribute significantly to refining models for Neanderthal disappearance in Europe and, more broadly, show that chronometric models regarding the appearance or disappearance of animal or hominin groups should be based only on radiocarbon dates obtained using robust pretreatment methods.

Functional traits of the world's late Quaternary large-bodied avian and mammalian herbivores

Prehistoric and recent extinctions of large-bodied terrestrial herbivores had significant and lasting impacts on Earth's ecosystems due to the loss of their distinct trait combinations. The world's surviving large-bodied avian and mammalian herbivores remain among the most threatened taxa. As such, a greater understanding of the ecological impacts of large herbivore losses is increasingly important. However, comprehensive and ecologically-relevant trait datasets for extinct and extant herbivores are lacking. Here, we present HerbiTraits, a comprehensive functional trait dataset for all late Quaternary terrestrial avian and mammalian herbivores ≥10 kg (545 species). HerbiTraits includes key traits that influence how herbivores interact with ecosystems, namely body mass, diet, fermentation type, habitat use, and limb morphology. Trait data were compiled from 557 sources and comprise the best available knowledge on late Quaternary large-bodied herbivores. HerbiTraits provides a tool for the analysis of herbivore functional diversity both past and present and its effects on Earth's ecosystems.

Bone need not remain an elephant in the room for radiocarbon dating

Radiocarbon (¹⁴C) analysis of skeletal remains by accelerator mass spectrometry is an essential tool in multiple branches of science. However, bone ¹⁴C dating results can be inconsistent and not comparable due to disparate laboratory pretreatment protocols that remove contamination. And, pretreatments are rarely discussed or reported by end-users, making it an 'elephant in the room' for Quaternary scientists. Through a questionnaire survey, I quantified consensus on the reliability of collagen pretreatments for ¹⁴C dating across 132 experts (25 countries). I discovered that while more than 95% of the audience was wary of contamination and would avoid gelatinization alone (minimum pretreatment used by most ¹⁴C facilities), 52% asked laboratories to choose the pretreatment method for them, and 58% could not rank the reliability of at least one pretreatment. Ultrafiltration was highly popular, and purification by XAD resins seemed restricted to American researchers. Isolating and dating the amino acid hydroxyproline was perceived as the most reliable pretreatment, but is expensive, time-consuming and not widely available. Solid evidence supports that only molecular-level dating accommodates all known bone contaminants and guarantees complete removal of humic and fulvic acids and conservation substances, with three key areas of progress: (i) innovation and more funded research is required to develop affordable analytical chemistry that can handle low-mass samples of collagen amino acids, (ii) a certification agency overseeing dating-quality control is needed to enhance methodological reproducibility and dating accuracy among laboratories, and (iii) more cross-disciplinary work with better ¹⁴C reporting etiquette will promote the integration of ¹⁴C dating across disciplines. Those developments could conclude long-standing debates based on low-accuracy data used to build chronologies for animal domestications, human/megafauna extirpations and migrations, archaeology, palaeoecology, palaeontology and palaeoclimate models.

The origin of Rhinocerotoidea and phylogeny of Ceratomorpha (Mammalia, Perissodactyla)

Rhinoceroses have been considered to have originated from tapiroids in the middle Eocene; however, the transition remains controversial, and the first unequivocal rhinocerotoids appeared about 4 Ma later than the earliest tapiroids of the Early Eocene. Here we describe 5 genera and 6 new species of rhinoceroses recently discovered from the early Eocene to the early middle Eocene deposits of the Erlian Basin of Inner Mongolia, China. These new materials represent the earliest members of rhinocerotoids, forstercooperiids, and/or hyrachyids, and bridge the evolutionary gap between the early Eocene ceratomorphs and middle Eocene rhinocerotoids. The phylogenetic analyses using parsimony and Bayesian inference methods support their affinities with rhinocerotoids, and also illuminate the phylogenetic relationships and biogeography of Ceratomorpha, although some discrepancies are present between the two criteria. The nearly contemporary occurrence of various rhinocerotoids indicates that the divergence of different rhinocerotoid groups occurred no later than the late early Eocene, which is soon after the split between the rhinocerotoids and the tapiroids in the early early Eocene. However, the Bayesian tip-dating estimate suggests that the divergence of different ceratomorph groups occurred in the middle Paleocene.

Heavy reliance on plants for Romanian cave bears evidenced by amino acid nitrogen isotope analysis

Heavy reliance on plants is rare in Carnivora and mostly limited to relatively small species in subtropical settings. The feeding behaviors of extinct cave bears living during Pleistocene cold periods at middle latitudes have been intensely studied using various approaches including isotopic analyses of fossil collagen. In contrast to cave bears from all other regions in Europe, some individuals from Romania show exceptionally high δ¹⁵N values that might be indicative of meat consumption. Herbivory on plants with high δ¹⁵N values cannot be ruled out based on this method, however. Here we apply an approach using the δ¹⁵N values of individual amino acids from collagen that offsets the baseline δ¹⁵N variation among environments. The analysis yielded strong signals of reliance on plants for Romanian cave bears based on the δ¹⁵N values of glutamate and phenylalanine. These results could suggest that the high variability in bulk collagen δ¹⁵N values observed among cave bears in Romania reflects niche partitioning but in a general trophic context of herbivory.

Preparative HPLC Separation of Underivatized Amino Acids for Isotopic Analysis

Single-compound analysis of stable or radioactive isotopes has found application in a number of fields ranging from archaeology to forensics. Often, the most difficult part of these analyses is the development of a method for isolating the compound(s) of interest, which can derive from a wide range of sample types including the hair, nails, and bone.Here we describe three complementary preparative HPLC techniques suitable for separating and isolating amino acids from bone collagen and hair keratin. Using preparative reversed-phase, ion-pair, or mixed-mode chromatography in aqueous carbon-free mobile phases, or those from which carbon can easily be removed, underivatized single amino acids can be isolated and further analyzed using mass spectrometric techniques.

Recent mitochondrial lineage extinction in the critically endangered Javan rhinoceros

The Javan rhinoceros (Rhinoceros sondaicus) is one of five extant rhinoceros species and among the rarest large mammals on Earth. Once widespread across Southeast Asia, it is now on the verge of extinction, with only one wild population remaining (estimated at ~60 individuals) on the island of Java, Indonesia. To assess the past genetic diversity of the female lineage of R. sondaicus, we generated mitochondrial genome data from eight museum specimens dating back to the 19th century, before the range of the Javan rhinoceros was dramatically reduced, for comparison against mitochondrial DNA (mtDNA) sequences of current R. sondaicus and other rhinoceros species. We succeeded in reconstructing five full and three partial ancient mitogenomes from the eight samples. We used BEAST to assess the phylogenetic relationship of the five extant rhinoceros species and the historical samples. The results show that the oldest and most diverse mtDNA lineages of R. sondaicus are found in historical samples, indicating a significant reduction of mtDNA diversity in modern Javan rhinos. We anticipate that the newly sequenced data will represent a useful resource for improving our understanding of evolutionary history of this species, should future studies be able to increase the available dataset. We hope this information may help in conservation efforts for this species.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Micro Methods for Megafauna: Novel Approaches to Late Quaternary Extinctions and Their Contributions to Faunal Conservation in the Anthropocene

Drivers of Late Quaternary megafaunal extinctions are relevant to modern conservation policy in a world of growing human population density, climate change, and faunal decline. Traditional debates tend toward global solutions, blaming either dramatic climate change or dispersals of Homo sapiens to new regions. Inherent limitations to archaeological and paleontological data sets often require reliance on scant, poorly resolved lines of evidence. However, recent developments in scientific technologies allow for more local, context-specific approaches. In the present article, we highlight how developments in five such methodologies (radiocarbon approaches, stable isotope analysis, ancient DNA, ancient proteomics, microscopy) have helped drive detailed analysis of specific megafaunal species, their particular ecological settings, and responses to new competitors or predators, climate change, and other external phenomena. The detailed case studies of faunal community composition, extinction chronologies, and demographic trends enabled by these methods examine megafaunal extinctions at scales appropriate for practical understanding of threats against particular species in their habitats today.

Introducing Trait Networks to Elucidate the Fluidity of Organismal Evolution Using Palaeontological Data

Explaining the evolution of animals requires ecological, developmental, paleontological, and phylogenetic considerations because organismal traits are affected by complex evolutionary processes. Modeling a plurality of processes, operating at distinct time-scales on potentially interdependent traits, can benefit from approaches that are complementary treatments to phylogenetics. Here, we developed an inclusive network approach, implemented in the command line software ComponentGrapher, and analyzed trait co-occurrence of rhinocerotoid mammals. We identified stable, unstable, and pivotal traits, as well as traits contributing to complexes, that may follow to a common developmental regulation, that point to an early implementation of the postcranial Bauplan among rhinocerotoids. Strikingly, most identified traits are highly dissociable, used repeatedly in distinct combinations and in different taxa, which usually do not form clades. Therefore, the genes encoding these traits are likely recruited into novel gene regulation networks during the course of evolution. Our evo-systemic framework, generalizable to other evolved organizations, supports a pluralistic modeling of organismal evolution, including trees and networks.

Palaeogenomics of the Hydrocarbon Producing Microalga Botryococcus braunii

Botryococcus braunii is a colonial microalga that appears early in the fossil record and is a sensitive proxy of environmental and hydroclimatic conditions. Palaeozoic Botryococcus fossils which contribute up to 90% of oil shales and approximately 1% of crude oil, co-localise with diagnostic geolipids from the degradation of source-signature hydrocarbons. However more recent Holocene sediments demonstrate no such association. Consequently, Botryococcus are identified in younger sediments by morphology alone, where potential misclassifications could lead to inaccurate paleoenvironmental reconstructions. Here we show that a combination of flow cytometry and ancient DNA (aDNA) sequencing can unambiguously identify Botryococcus microfossils in Holocene sediments with hitherto unparalleled accuracy and rapidity. The application of aDNA sequencing to microfossils offers a far-reaching opportunity for understanding environmental change in the recent geological record. When allied with other high-resolution palaeoenvironmental information such as aDNA sequencing of humans and megafauna, aDNA from microfossils may allow a deeper and more precise understanding of past environments, ecologies and migrations.

Compound-specific radiocarbon dating and mitochondrial DNA analysis of the Pleistocene hominin from Salkhit Mongolia

A skullcap found in the Salkhit Valley in northeast Mongolia is, to our knowledge, the only Pleistocene hominin fossil found in the country. It was initially described as an individual with possible archaic affinities, but its ancestry has been debated since the discovery. Here, we determine the age of the Salkhit skull by compound-specific radiocarbon dating of hydroxyproline to 34,950-33,900 Cal. BP (at 95% probability), placing the Salkhit individual in the Early Upper Paleolithic period. We reconstruct the complete mitochondrial genome (mtDNA) of the specimen. It falls within a group of modern human mtDNAs (haplogroup N) that is widespread in Eurasia today. The results now place the specimen into its proper chronometric and biological context and allow us to begin integrating it with other evidence for the human occupation of this region during the Paleolithic, as well as wider Pleistocene sequences across Eurasia.