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.

Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears.

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans.

52,000 years of woolly rhinoceros population dynamics reveal extinction mechanisms.

The extinction of the woolly rhinoceros (Coelodonta antiquitatis) at the onset of the Holocene remains an enigma, with conflicting evidence regarding its cause and spatiotemporal dynamics. This partly reflects challenges in determining demographic responses of late Quaternary megafauna to climatic and anthropogenic causal drivers with available genetic and paleontological techniques. Here, we show that elucidating mechanisms of ancient extinctions can benefit from a detailed understanding of fine-scale metapopulation dynamics, operating over many millennia. Using an abundant fossil record, ancient DNA, and high-resolution simulation models, we untangle the ecological mechanisms and causal drivers that are likely to have been integral in the decline and later extinction of the woolly rhinoceros. Our 52,000-y reconstruction of distribution-wide metapopulation dynamics supports a pathway to extinction that began long before the Holocene, when the combination of cooling temperatures and low but sustained hunting by humans trapped woolly rhinoceroses in suboptimal habitats along the southern edge of their range. Modeling indicates that this ecological trap intensified after the end of the last ice age, preventing colonization of newly formed suitable habitats, weakening stabilizing metapopulation processes, triggering the extinction of the woolly rhinoceros in the early Holocene. Our findings suggest that fragmentation and resultant metapopulation dynamics should be explicitly considered in explanations of late Quaternary megafauna extinctions, sending a clarion call to the fragility of the remaining large-bodied grazers restricted to disjunct fragments of poor-quality habitat due to anthropogenic environmental change.

Elucidating the sustainability of 700 y of Inuvialuit beluga whale hunting in the Mackenzie River Delta, Northwest Territories, Canada.

Beluga whales play a critical role in the subsistence economies and cultural heritage of Indigenous communities across the Arctic, yet the effects of Indigenous hunting on beluga whales remain unknown. Here, we integrate paleogenomics, genetic simulations, and stable δ13C and δ15N isotope analysis to investigate 700 y of beluga subsistence hunting in the Mackenzie Delta area of northwestern Canada. Genetic identification of the zooarchaeological remains, which is based on radiocarbon dating, span three time periods (1290 to 1440 CE; 1450 to 1650 CE; 1800 to 1870 CE), indicates shifts across time in the sex ratio of the harvested belugas. The equal number of females and males harvested in 1450 to 1650 CE versus more males harvested in the two other time periods may reflect changes in hunting practices or temporal shifts in beluga availability. We find temporal shifts and sex-based differences in δ13C of the harvested belugas across time, suggesting historical adaptability in the foraging ecology of the whales. We uncovered distinct mitochondrial diversity unique to the Mackenzie Delta belugas, but found no changes in nuclear genomic diversity nor any substructuring across time. Our findings indicate the genomic stability and continuity of the Mackenzie Delta beluga population across the 700 y surveyed, indicating the impact of Inuvialuit subsistence harvests on the genetic diversity of contemporary beluga individuals has been negligible.

Late Pleistocene polar bear genomes reveal the timing of allele fixation in key genes associated with Arctic adaptation.

The polar bear (Ursus maritimus) occupies a relatively narrow ecological niche, with many traits adapted for cold temperatures, movement across snow, ice and open water, and for consuming highly lipid-dense prey species. The divergence of polar bears from brown bears (Ursus arctos) and their adaptation to their Arctic lifestyle is a well-known example of rapid evolution. Previous research investigating whole genomes uncovered twelve key genes that are highly differentiated between polar and brown bears, show signatures of selection in the polar bear lineage, and are associated with polar bear adaptation to the Arctic environment. Further research suggested fixed derived alleles in these genes arose from selection on both standing variation and de novo mutations in the evolution of polar bears. Here, we reevaluate these findings based on a larger and geographically more representative dataset of 119 polar bears and 135 brown bears, and assess the timing of derived allele fixation in polar bears by incorporating the genomes of two Late Pleistocene individuals (aged 130-100,000 years old and 100-70,000 years old). In contrast with previous results, we found no evidence of derived alleles fixed in present-day polar bears within the key genes arising from de novo mutation. Most derived alleles fixed in present-day polar bears were also fixed in the Late Pleistocene polar bears, suggesting selection occurred prior to 70,000 years ago. However, some derived alleles fixed in present-day polar bears were not fixed in the two Late Pleistocene polar bears, including at sites within APOB, LYST, and TTN. These three genes are associated with cardiovascular function, metabolism, and pigmentation, suggesting selection may have acted on different loci at different times.

The origins and diversification of Holarctic brown bear populations inferred from genomes of past and present populations.

The brown bear (Ursus arctos) is one of the survivors of the Late Quaternary megafauna extinctions. However, despite being widely distributed across the Holarctic, brown bears have experienced extensive range reductions, and even extirpations in some geographical regions. Previous research efforts using genetic data have provided valuable insights into their evolutionary history. However, most studies have been limited to contemporary individuals or mitochondrial DNA, limiting insights into population processes that preceded the present. Here, we present genomic data from two Late Pleistocene brown bears from Honshu, Japan and eastern Siberia, and combine them with published contemporary and ancient genomes from across the Holarctic range of brown bears to investigate the evolutionary relationships among brown bear populations through time and space. By including genomic data from Late Pleistocene and Holocene individuals sampled outside the current distribution range, we uncover diversity not present in contemporary populations. Notably, although contemporary individuals display geographically structured populations most likely driven by isolation-by-distance, this pattern varies among the ancient samples across different regions. The inclusion of ancient brown bears in our analysis provides novel insights into the evolutionary history of brown bears and contributes to understanding the populations and diversity lost during the Late Quaternary.

Topographic barriers drive the pronounced genetic subdivision of a range-limited fossorial rodent.

Due to their limited dispersal ability, fossorial species with predominantly belowground activity usually show increased levels of population subdivision across relatively small spatial scales. This may be exacerbated in harsh mountain ecosystems, where landscape geomorphology limits species' dispersal ability and leads to small effective population sizes, making species relatively vulnerable to environmental change. To better understand the environmental drivers of species' population subdivision in remote mountain ecosystems, particularly in understudied high-elevation systems in Africa, we studied the giant root-rat (Tachyoryctes macrocephalus), a fossorial rodent confined to the afro-alpine ecosystem of the Bale Mountains in Ethiopia. Using mitochondrial and low-coverage nuclear genomes, we investigated 77 giant root-rat individuals sampled from nine localities across its entire ~1000 km2 range. Our data revealed a distinct division into a northern and southern group, with no signs of gene flow, and higher nuclear genetic diversity in the south. Landscape genetic analyses of the mitochondrial and nuclear genomes indicated that population subdivision was driven by slope and elevation differences of up to 500 m across escarpments separating the north and south, potentially reinforced by glaciation of the south during the Late Pleistocene (~42,000-16,000 years ago). Despite this landscape-scale subdivision between the north and south, weak geographic structuring of sampling localities within regions indicated gene flow across distances of at least 16 km at the local scale, suggesting high, aboveground mobility for relatively long distances. Our study highlights that despite the potential for local-scale gene flow in fossorial species, topographic barriers can result in pronounced genetic subdivision. These factors can reduce genetic variability, which should be considered when developing conservation strategies.

A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea.

The impact of post-divergence gene flow in speciation has been documented across a range of taxa in recent years, and may have been especially widespread in highly mobile, wide-ranging marine species, such as cetaceans. Here, we studied individual genomes from nine species across the three families of the toothed whale superfamily Delphinoidea (Delphinidae, Phocoenidae and Monodontidae). To investigate the role of post-divergence gene flow in the speciation process, we used a multifaceted approach, including (i) phylogenomics, (ii) the distribution of shared derived alleles and (iii) demographic inference. We found the divergence of lineages within Delphinoidea did not follow a process of pure bifurcation, but was much more complex. Sliding-window phylogenomics reveal a high prevalence of discordant topologies within the superfamily, with further analyses indicating these discordances arose due to both incomplete lineage sorting and gene flow. D-statistics and f-branch analyses supported gene flow between members of Delphinoidea, with the vast majority of gene flow occurring as ancient interfamilial events. Demographic analyses provided evidence that introgressive gene flow has likely ceased between all species pairs tested, despite reports of contemporary interspecific hybrids. Our study provides the first steps towards resolving the large complexity of speciation within Delphinoidea; we reveal the prevalence of ancient interfamilial gene flow events prior to the diversification of each family, and suggest that contemporary hybridisation events may be disadvantageous, as hybrid individuals do not appear to contribute to the parental species' gene pools.

Faster ocean warming threatens richest areas of marine biodiversity.

The vulnerability of marine biodiversity to accelerated rates of climatic change is poorly understood. By developing a new method for identifying extreme oceanic warming events during Earth's most recent deglaciation, and comparing these to 21st century projections, we show that future rates of ocean warming will disproportionately affect the most speciose marine communities, potentially threatening biodiversity in more than 70% of current-day global hotspots of marine species richness. The persistence of these richest areas of marine biodiversity will require many species to move well beyond the biogeographic realm where they are endemic, at rates of redistribution not previously seen. Our approach for quantifying exposure of biodiversity to past and future rates of oceanic warming provides new context and scalable information for deriving and strengthening conservation actions to safeguard marine biodiversity under climate change.

Circumpolar phylogeography and demographic history of beluga whales reflect past climatic fluctuations.

Several Arctic marine mammal species are predicted to be negatively impacted by rapid sea ice loss associated with ongoing ocean warming. However, consequences for Arctic whales remain uncertain. To investigate how Arctic whales responded to past climatic fluctuations, we analysed 206 mitochondrial genomes from beluga whales (Delphinapterus leucas) sampled across their circumpolar range, and four nuclear genomes, covering both the Atlantic and the Pacific Arctic region. We found four well-differentiated mitochondrial lineages, which were established before the onset of the last glacial expansion ~110 thousand years ago. Our findings suggested these lineages diverged in allopatry, reflecting isolation of populations during glacial periods when the Arctic sea-shelf was covered by multiyear sea ice. Subsequent population expansion and secondary contact between the Atlantic and Pacific Oceans shaped the current geographic distribution of lineages, and may have facilitated mitochondrial introgression. Our demographic reconstructions based on both mitochondrial and nuclear genomes showed markedly lower population sizes during the Last Glacial Maximum (LGM) compared to the preceding Eemian and current Holocene interglacial periods. Habitat modelling similarly revealed less suitable habitat during the LGM (glacial) than at present (interglacial). Together, our findings suggested the association between climate, population size, and available habitat in belugas. Forecasts for year 2100 showed that beluga habitat will decrease and shift northwards as oceans continue to warm, putatively leading to population declines in some beluga populations. Finally, we identified vulnerable populations which, if extirpated as a consequence of ocean warming, will lead to a substantial decline of species-wide haplotype diversity.

Unraveling elephant-shrews: Phylogenetic relationships and unexpected introgression among giant sengis.

Giant sengis, or elephant-shrews (Macroscelidea; Macroscelididae; Rhynchocyon), are small-bodied mammals found in central and eastern African forests. Studies have provided contrasting views of the extent and direction of introgression among species. We generated full mitochondrial genomes, and compiled publically available mtDNA 12S and nuclear vWF sequences from Rhynchocyon cirnei, R. petersi and R. udzungwensis that had not previously been analyzed in concert, to elucidate the phylogenetic and population-specific context of potential introgression. Our spatially and phylogenetically broad sampling across species revealed substantial, unidirectional mitochondrial introgression of the R. petersi lineage into R. cirnei reichardi and R. udzungwensis, and from R. udzungwensis into R. c. reichardi. All introgression was highly localized and found only in the eastern Udzungwa Mountains forests in Tanzania. The nuclear data showed another pattern, with R. petersi haplotypes in R. cirnei cirnei and R. c. reichardi. No individuals showed both mitochondrial and nuclear introgression. Our results suggest higher levels of hybridization among giant sengi species than previously recognized, but also highlight the need for further genome-wide analysis and increased spatial sampling to clarify the many aspects of diversification and introgression in this group.

Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role.

BACKGROUND: Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage. RESULTS: We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage. CONCLUSIONS: Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.

Influence of past climate change on phylogeography and demographic history of narwhals, Monodon monoceros.

The Arctic is warming at an unprecedented rate, with unknown consequences for endemic fauna. However, Earth has experienced severe climatic oscillations in the past, and understanding how species responded to them might provide insight into their resilience to near-future climatic predictions. Little is known about the responses of Arctic marine mammals to past climatic shifts, but narwhals (Monodon monoceros) are considered one of the endemic Arctic species most vulnerable to environmental change. Here, we analyse 121 complete mitochondrial genomes from narwhals sampled across their range and use them in combination with species distribution models to elucidate the influence of past and ongoing climatic shifts on their population structure and demographic history. We find low levels of genetic diversity and limited geographic structuring of genetic clades. We show that narwhals experienced a long-term low effective population size, which increased after the Last Glacial Maximum, when the amount of suitable habitat expanded. Similar post-glacial habitat release has been a key driver of population size expansion of other polar marine predators. Our analyses indicate that habitat availability has been critical to the success of narwhals, raising concerns for their fate in an increasingly warming Arctic.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

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.

Persistent organic pollutants, skull size and bone density of polar bears (Ursus maritimus) from East Greenland 1892-2015 and Svalbard 1964-2004.

We investigated skull size (condylobasal length; CBL) and bone mineral density (BMD) in polar bears (Ursus maritimus) from East Greenland (n = 307) and Svalbard (n = 173) sampled during the period 1892-2015 in East Greenland and 1964-2004 at Svalbard. Adult males from East Greenland showed a continuous decrease in BMD from 1892 to 2015 (linear regression: p < 0.01) indicating that adult male skulls collected in the early pre-pollution period had the highest BMD. A similar decrease in BMD over time was not found for the East Greenland adult females. However, there was a non-significant trend that the skull size of adult East Greenland females was negatively correlated with collection year 1892-2015 (linear regression: p = 0.06). No temporal change was found for BMD or skull size in Svalbard polar bears (ANOVA: all p > 0.05) nor was there any significant difference in BMD between Svalbard and East Greenland subpopulations. Skull size was larger in polar bears from Svalbard than from East Greenland (two-way ANOVA: p = 0.003). T-scores reflecting risk of osteoporosis showed that adult males from both East Greenland and Svalbard are at risk of developing osteopenia. Finally, when correcting for age and sex, BMD in East Greenland polar bears increased with increasing concentrations of persistent organic pollutants (POPs) i.e. ΣPCB (polychlorinated biphenyls), ΣHCH (hexachlorohexane), HCB (hexachlorobenzene) and ΣPBDE (polybrominated diphenyl ethers) while skull size increased with ΣHCH concentrations all in the period 1999-2014 (multiple linear regression: all p < 0.05, n = 175). The results suggest that environmental changes over time, including exposure to POPs, may affect bone density and size of polar bears.

Species-specific responses of Late Quaternary megafauna to climate and humans.

Despite decades of research, the roles of climate and humans in driving the dramatic extinctions of large-bodied mammals during the Late Quaternary period remain contentious. Here we use ancient DNA, species distribution models and the human fossil record to elucidate how climate and humans shaped the demographic history of woolly rhinoceros, woolly mammoth, wild horse, reindeer, bison and musk ox. We show that climate has been a major driver of population change over the past 50,000 years. However, each species responds differently to the effects of climatic shifts, habitat redistribution and human encroachment. Although climate change alone can explain the extinction of some species, such as Eurasian musk ox and woolly rhinoceros, a combination of climatic and anthropogenic effects appears to be responsible for the extinction of others, including Eurasian steppe bison and wild horse. We find no genetic signature or any distinctive range dynamics distinguishing extinct from surviving species, emphasizing the challenges associated with predicting future responses of extant mammals to climate and human-mediated habitat change.

Extinct New Zealand megafauna were not in decline before human colonization.

The extinction of New Zealand's moa (Aves: Dinornithiformes) followed the arrival of humans in the late 13th century and was the final event of the prehistoric Late Quaternary megafauna extinctions. Determining the state of the moa populations in the pre-extinction period is fundamental to understanding the causes of the event. We sampled 281 moa individuals and combined radiocarbon dating with ancient DNA analyses to help resolve the extinction debate and gain insights into moa biology. The samples, which were predominantly from the last 4,000 years preceding the extinction, represent four sympatric moa species excavated from five adjacent fossil deposits. We characterized the moa assemblage using mitochondrial DNA and nuclear microsatellite markers developed specifically for moa. Although genetic diversity differed significantly among the four species, we found that the millennia preceding the extinction were characterized by a remarkable degree of genetic stability in all species, with no loss of heterozygosity and no shifts in allele frequencies over time. The extinction event itself was too rapid to be manifested in the moa gene pools. Contradicting previous claims of a decline in moa before Polynesian settlement in New Zealand, our findings indicate that the populations were large and stable before suddenly disappearing. This interpretation is supported by approximate Bayesian computation analyses. Our analyses consolidate the disappearance of moa as the most rapid, human-facilitated megafauna extinction documented to date.

Ancient human genome sequence of an extinct Palaeo-Eskimo.

We report here the genome sequence of an ancient human. Obtained from approximately 4,000-year-old permafrost-preserved hair, the genome represents a male individual from the first known culture to settle in Greenland. Sequenced to an average depth of 20x, we recover 79% of the diploid genome, an amount close to the practical limit of current sequencing technologies. We identify 353,151 high-confidence single-nucleotide polymorphisms (SNPs), of which 6.8% have not been reported previously. We estimate raw read contamination to be no higher than 0.8%. We use functional SNP assessment to assign possible phenotypic characteristics of the individual that belonged to a culture whose location has yielded only trace human remains. We compare the high-confidence SNPs to those of contemporary populations to find the populations most closely related to the individual. This provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of that giving rise to the modern Native Americans and Inuit.

Phylogenetic estimation of timescales using ancient DNA: the effects of temporal sampling scheme and uncertainty in sample ages.

In recent years, ancient DNA has increasingly been used for estimating molecular timescales, particularly in studies of substitution rates and demographic histories. Molecular clocks can be calibrated using temporal information from ancient DNA sequences. This information comes from the ages of the ancient samples, which can be estimated by radiocarbon dating the source material or by dating the layers in which the material was deposited. Both methods involve sources of uncertainty. The performance of bayesian phylogenetic inference depends on the information content of the data set, which includes variation in the DNA sequences and the structure of the sample ages. Various sources of estimation error can reduce our ability to estimate rates and timescales accurately and precisely. We investigated the impact of sample-dating uncertainties on the estimation of evolutionary timescale parameters using the software BEAST. Our analyses involved 11 published data sets and focused on estimates of substitution rate and root age. We show that, provided that samples have been accurately dated and have a broad temporal span, it might be unnecessary to account for sample-dating uncertainty in Bayesian phylogenetic analyses of ancient DNA. We also investigated the sample size and temporal span of the ancient DNA sequences needed to estimate phylogenetic timescales reliably. Our results show that the range of sample ages plays a crucial role in determining the quality of the results but that accurate and precise phylogenetic estimates of timescales can be made even with only a few ancient sequences. These findings have important practical consequences for studies of molecular rates, timescales, and population dynamics.