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 Paleogenomic Reconstruction of the Deep Population History of the Andes.

There are many unanswered questions about the population history of the Central and South Central Andes, particularly regarding the impact of large-scale societies, such as the Moche, Wari, Tiwanaku, and Inca. We assembled genome-wide data on 89 individuals dating from ∼9,000-500 years ago (BP), with a particular focus on the period of the rise and fall of state societies. Today's genetic structure began to develop by 5,800 BP, followed by bi-directional gene flow between the North and South Highlands, and between the Highlands and Coast. We detect minimal admixture among neighboring groups between ∼2,000-500 BP, although we do detect cosmopolitanism (people of diverse ancestries living side-by-side) in the heartlands of the Tiwanaku and Inca polities. We also highlight cases of long-range mobility connecting the Andes to Argentina and the Northwest Andes to the Amazon Basin. VIDEO ABSTRACT.

Reconstructing the Deep Population History of Central and South America.

We report genome-wide ancient DNA from 49 individuals forming four parallel time transects in Belize, Brazil, the Central Andes, and the Southern Cone, each dating to at least ∼9,000 years ago. The common ancestral population radiated rapidly from just one of the two early branches that contributed to Native Americans today. We document two previously unappreciated streams of gene flow between North and South America. One affected the Central Andes by ∼4,200 years ago, while the other explains an affinity between the oldest North American genome associated with the Clovis culture and the oldest Central and South Americans from Chile, Brazil, and Belize. However, this was not the primary source for later South Americans, as the other ancient individuals derive from lineages without specific affinity to the Clovis-associated genome, suggesting a population replacement that began at least 9,000 years ago and was followed by substantial population continuity in multiple regions.

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.

Global impacts of an extreme solar particle event under different geomagnetic field strengths.

Solar particle events (SPEs) are short-lived bursts of high-energy particles from the solar atmosphere and are widely recognized as posing significant economic risks to modern society. Most SPEs are relatively weak and have minor impacts on the Earth's environment, but historic records contain much stronger SPEs which have the potential to alter atmospheric chemistry, impacting climate and biological life. The impacts of such strong SPEs would be far more severe when the Earth's protective geomagnetic field is weak, such as during past geomagnetic excursions or reversals. Here, we model the impacts of an extreme SPE under different geomagnetic field strengths, focusing on changes in atmospheric chemistry and surface radiation using the atmosphere-ocean-chemistry-climate model SOCOL3-MPIOM and the radiation transfer model LibRadtran. Under current geomagnetic conditions, an extreme SPE would increase NOx concentrations in the polar stratosphere and mesosphere, causing reductions in extratropical stratospheric ozone lasting for about a year. In contrast, with no geomagnetic field, there would be a substantial increase in NOx throughout the entire atmosphere, resulting in severe stratospheric ozone depletion for several years. The resulting ground-level ultraviolet (UV) radiation would remain elevated for up to 6 y, leading to increases in UV index up to 20 to 25% and solar-induced DNA damage rates by 40 to 50%. The potential evolutionary impacts of past extreme SPEs remain an important question, while the risks they pose to human health in modern conditions continue to be underestimated.

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.

The role of genetic selection and climatic factors in the dispersal of anatomically modern humans out of Africa.

The evolutionarily recent dispersal of anatomically modern humans (AMH) out of Africa (OoA) and across Eurasia provides a unique opportunity to examine the impacts of genetic selection as humans adapted to multiple new environments. Analysis of ancient Eurasian genomic datasets (~1,000 to 45,000 y old) reveals signatures of strong selection, including at least 57 hard sweeps after the initial AMH movement OoA, which have been obscured in modern populations by extensive admixture during the Holocene. The spatiotemporal patterns of these hard sweeps provide a means to reconstruct early AMH population dispersals OoA. We identify a previously unsuspected extended period of genetic adaptation lasting ~30,000 y, potentially in the Arabian Peninsula area, prior to a major Neandertal genetic introgression and subsequent rapid dispersal across Eurasia as far as Australia. Consistent functional targets of selection initiated during this period, which we term the Arabian Standstill, include loci involved in the regulation of fat storage, neural development, skin physiology, and cilia function. Similar adaptive signatures are also evident in introgressed archaic hominin loci and modern Arctic human groups, and we suggest that this signal represents selection for cold adaptation. Surprisingly, many of the candidate selected loci across these groups appear to directly interact and coordinately regulate biological processes, with a number associated with major modern diseases including the ciliopathies, metabolic syndrome, and neurodegenerative disorders. This expands the potential for ancestral human adaptation to directly impact modern diseases, providing a platform for evolutionary medicine.

Systematic benchmark of ancient DNA read mapping.

The current standard practice for assembling individual genomes involves mapping millions of short DNA sequences (also known as DNA 'reads') against a pre-constructed reference genome. Mapping vast amounts of short reads in a timely manner is a computationally challenging task that inevitably produces artefacts, including biases against alleles not found in the reference genome. This reference bias and other mapping artefacts are expected to be exacerbated in ancient DNA (aDNA) studies, which rely on the analysis of low quantities of damaged and very short DNA fragments (~30-80 bp). Nevertheless, the current gold-standard mapping strategies for aDNA studies have effectively remained unchanged for nearly a decade, during which time new software has emerged. In this study, we used simulated aDNA reads from three different human populations to benchmark the performance of 30 distinct mapping strategies implemented across four different read mapping software-BWA-aln, BWA-mem, NovoAlign and Bowtie2-and quantified the impact of reference bias in downstream population genetic analyses. We show that specific NovoAlign, BWA-aln and BWA-mem parameterizations achieve high mapping precision with low levels of reference bias, particularly after filtering out reads with low mapping qualities. However, unbiased NovoAlign results required the use of an IUPAC reference genome. While relevant only to aDNA projects where reference population data are available, the benefit of using an IUPAC reference demonstrates the value of incorporating population genetic information into the aDNA mapping process, echoing recent results based on graph genome representations.

Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus.

Recent genomic data have revealed multiple interactions between Neanderthals and modern humans, but there is currently little genetic evidence regarding Neanderthal behaviour, diet, or disease. Here we describe the shotgun-sequencing of ancient DNA from five specimens of Neanderthal calcified dental plaque (calculus) and the characterization of regional differences in Neanderthal ecology. At Spy cave, Belgium, Neanderthal diet was heavily meat based and included woolly rhinoceros and wild sheep (mouflon), characteristic of a steppe environment. In contrast, no meat was detected in the diet of Neanderthals from El Sidrón cave, Spain, and dietary components of mushrooms, pine nuts, and moss reflected forest gathering. Differences in diet were also linked to an overall shift in the oral bacterial community (microbiota) and suggested that meat consumption contributed to substantial variation within Neanderthal microbiota. Evidence for self-medication was detected in an El Sidrón Neanderthal with a dental abscess and a chronic gastrointestinal pathogen (Enterocytozoon bieneusi). Metagenomic data from this individual also contained a nearly complete genome of the archaeal commensal Methanobrevibacter oralis (10.2× depth of coverage)-the oldest draft microbial genome generated to date, at around 48,000 years old. DNA preserved within dental calculus represents a notable source of information about the behaviour and health of ancient hominin specimens, as well as a unique system that is useful for the study of long-term microbial evolution.

Aboriginal mitogenomes reveal 50,000 years of regionalism in Australia.

Aboriginal Australians represent one of the longest continuous cultural complexes known. Archaeological evidence indicates that Australia and New Guinea were initially settled approximately 50 thousand years ago (ka); however, little is known about the processes underlying the enormous linguistic and phenotypic diversity within Australia. Here we report 111 mitochondrial genomes (mitogenomes) from historical Aboriginal Australian hair samples, whose origins enable us to reconstruct Australian phylogeographic history before European settlement. Marked geographic patterns and deep splits across the major mitochondrial haplogroups imply that the settlement of Australia comprised a single, rapid migration along the east and west coasts that reached southern Australia by 49-45 ka. After continent-wide colonization, strong regional patterns developed and these have survived despite substantial climatic and cultural change during the late Pleistocene and Holocene epochs. Remarkably, we find evidence for the continuous presence of populations in discrete geographic areas dating back to around 50 ka, in agreement with the notable Aboriginal Australian cultural attachment to their country.

Extreme hydrothermal conditions at an active plate-bounding fault.

Temperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre. At temperatures above 300-450 degrees Celsius, usually found at depths greater than 10-15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional-mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults.

Early Last Interglacial ocean warming drove substantial ice mass loss from Antarctica.

The future response of the Antarctic ice sheet to rising temperatures remains highly uncertain. A useful period for assessing the sensitivity of Antarctica to warming is the Last Interglacial (LIG) (129 to 116 ky), which experienced warmer polar temperatures and higher global mean sea level (GMSL) (+6 to 9 m) relative to present day. LIG sea level cannot be fully explained by Greenland Ice Sheet melt (∼2 m), ocean thermal expansion, and melting mountain glaciers (∼1 m), suggesting substantial Antarctic mass loss was initiated by warming of Southern Ocean waters, resulting from a weakening Atlantic meridional overturning circulation in response to North Atlantic surface freshening. Here, we report a blue-ice record of ice sheet and environmental change from the Weddell Sea Embayment at the periphery of the marine-based West Antarctic Ice Sheet (WAIS), which is underlain by major methane hydrate reserves. Constrained by a widespread volcanic horizon and supported by ancient microbial DNA analyses, we provide evidence for substantial mass loss across the Weddell Sea Embayment during the LIG, most likely driven by ocean warming and associated with destabilization of subglacial hydrates. Ice sheet modeling supports this interpretation and suggests that millennial-scale warming of the Southern Ocean could have triggered a multimeter rise in global sea levels. Our data indicate that Antarctica is highly vulnerable to projected increases in ocean temperatures and may drive ice-climate feedbacks that further amplify warming.

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.

Using hominin introgression to trace modern human dispersals.

The dispersal of anatomically modern human populations out of Africa and across much of the rest of the world around 55 to 50 thousand years before present (ka) is recorded genetically by the multiple hominin groups they met and interbred with along the way, including the Neandertals and Denisovans. The signatures of these introgression events remain preserved in the genomes of modern-day populations, and provide a powerful record of the sequence and timing of these early migrations, with Asia proving a particularly complex area. At least 3 different hominin groups appear to have been involved in Asia, of which only the Denisovans are currently known. Several interbreeding events are inferred to have taken place east of Wallace's Line, consistent with archaeological evidence of widespread and early hominin presence in the area. However, archaeological and fossil evidence indicates archaic hominins had not spread as far as the Sahul continent (New Guinea, Australia, and Tasmania), where recent genetic evidence remains enigmatic.

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.

Genome-wide patterns of selection in 230 ancient Eurasians.

Ancient DNA makes it possible to observe natural selection directly by analysing samples from populations before, during and after adaptation events. Here we report a genome-wide scan for selection using ancient DNA, capitalizing on the largest ancient DNA data set yet assembled: 230 West Eurasians who lived between 6500 and 300 bc, including 163 with newly reported data. The new samples include, to our knowledge, the first genome-wide ancient DNA from Anatolian Neolithic farmers, whose genetic material we obtained by extracting from petrous bones, and who we show were members of the population that was the source of Europe's first farmers. We also report a transect of the steppe region in Samara between 5600 and 300 bc, which allows us to identify admixture into the steppe from at least two external sources. We detect selection at loci associated with diet, pigmentation and immunity, and two independent episodes of selection on height.

Massive migration from the steppe was a source for Indo-European languages in Europe.

We generated genome-wide data from 69 Europeans who lived between 8,000-3,000 years ago by enriching ancient DNA libraries for a target set of almost 400,000 polymorphisms. Enrichment of these positions decreases the sequencing required for genome-wide ancient DNA analysis by a median of around 250-fold, allowing us to study an order of magnitude more individuals than previous studies and to obtain new insights about the past. We show that the populations of Western and Far Eastern Europe followed opposite trajectories between 8,000-5,000 years ago. At the beginning of the Neolithic period in Europe, ∼8,000-7,000 years ago, closely related groups of early farmers appeared in Germany, Hungary and Spain, different from indigenous hunter-gatherers, whereas Russia was inhabited by a distinctive population of hunter-gatherers with high affinity to a ∼24,000-year-old Siberian. By ∼6,000-5,000 years ago, farmers throughout much of Europe had more hunter-gatherer ancestry than their predecessors, but in Russia, the Yamnaya steppe herders of this time were descended not only from the preceding eastern European hunter-gatherers, but also from a population of Near Eastern ancestry. Western and Eastern Europe came into contact ∼4,500 years ago, as the Late Neolithic Corded Ware people from Germany traced ∼75% of their ancestry to the Yamnaya, documenting a massive migration into the heartland of Europe from its eastern periphery. This steppe ancestry persisted in all sampled central Europeans until at least ∼3,000 years ago, and is ubiquitous in present-day Europeans. These results provide support for a steppe origin of at least some of the Indo-European languages of Europe.

Coprolites reveal ecological interactions lost with the extinction of New Zealand birds.

Over the past 50,000 y, biotic extinctions and declines have left a legacy of vacant niches and broken ecological interactions across global terrestrial ecosystems. Reconstructing the natural, unmodified ecosystems that preceded these events relies on high-resolution analyses of paleoecological deposits. Coprolites are a source of uniquely detailed information about trophic interactions and the behaviors, gut parasite communities, and microbiotas of prehistoric animal species. Such insights are critical for understanding the legacy effects of extinctions on ecosystems, and can help guide contemporary conservation and ecosystem restoration efforts. Here we use high-throughput sequencing (HTS) of ancient eukaryotic DNA from coprolites to reconstruct aspects of the biology and ecology of four species of extinct moa and the critically endangered kakapo parrot from New Zealand (NZ). Importantly, we provide evidence that moa and prehistoric kakapo consumed ectomycorrhizal fungi, suggesting these birds played a role in dispersing fungi that are key to NZ's natural forest ecosystems. We also provide the first DNA-based evidence that moa frequently supplemented their broad diets with ferns and mosses. Finally, we also find parasite taxa that provide insight into moa behavior, and present data supporting the hypothesis of coextinction between moa and several parasite species. Our study demonstrates that HTS sequencing of coprolites provides a powerful tool for resolving key aspects of ancient ecosystems and may rapidly provide information not obtainable by conventional paleoecological techniques, such as fossil analyses.

When did Homo sapiens first reach Southeast Asia and Sahul?

Anatomically modern humans (Homo sapiens, AMH) began spreading across Eurasia from Africa and adjacent Southwest Asia about 50,000-55,000 years ago (ca 50-55 ka). Some have argued that human genetic, fossil, and archaeological data indicate one or more prior dispersals, possibly as early as 120 ka. A recently reported age estimate of 65 ka for Madjedbebe, an archaeological site in northern Sahul (Pleistocene Australia-New Guinea), if correct, offers what might be the strongest support yet presented for a pre-55-ka African AMH exodus. We review evidence for AMH arrival on an arc spanning South China through Sahul and then evaluate data from Madjedbebe. We find that an age estimate of >50 ka for this site is unlikely to be valid. While AMH may have moved far beyond Africa well before 50-55 ka, data from the region of interest offered in support of this idea are not compelling.

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.