Demographic reconstruction from ancient DNA supports rapid extinction of the great auk

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Ancient DNA is preserved in fish fossils from tropical lake sediments

Tropical freshwater lakes are well known for their high biodiversity, and particularly the East African Great Lakes are renowned for their adaptive radiation of cichlid fishes. While comparative phylogenetic analyses of extant species flocks have revealed patterns and processes of their diversification, little is known about evolutionary trajectories within lineages, the impacts of environmental drivers, or the scope and nature of now-extinct diversity. Time-structured palaeodata from geologically young fossil records, such as fossil counts and particularly ancient DNA (aDNA) data, would help fill this large knowledge gap. High ambient temperatures can be detrimental to the preservation of DNA, but refined methodology now allows data generation even from very poorly preserved samples. Here, we show for the first time that fish fossils from tropical lake sediments yield endogenous aDNA. Despite generally low endogenous content and high sample dropout, the application of high-throughput sequencing and, in some cases, sequence capture allowed taxonomic assignment and phylogenetic placement of 17% of analysed fish fossils to family or tribe level, including remains which are up to 2700 years old or weigh less than 1 mg. The relationship between aDNA degradation and the thermal age of samples is similar to that described for terrestrial samples from cold environments when adjusted for elevated temperature. Success rates and aDNA preservation differed between the investigated lakes Chala, Kivu and Victoria, possibly caused by differences in bottom water oxygenation. Our study demonstrates that the sediment records of tropical lakes can preserve genetic information on rapidly diversifying fish taxa over time scales of millennia.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

The prelude to industrial whaling: identifying the targets of ancient European whaling using zooarchaeology and collagen mass-peptide fingerprinting

Taxonomic identification of whale bones found during archaeological excavations is problematic due to their typically fragmented state. This difficulty limits understanding of both the past spatio-temporal distributions of whale populations and of possible early whaling activities. To overcome this challenge, we performed zooarchaeology by mass spectrometry on an unprecedented 719 archaeological and palaeontological specimens of probable whale bone from Atlantic European contexts, predominantly dating from ca 3500 BCE to the eighteenth century CE. The results show high numbers of Balaenidae (many probably North Atlantic right whale (Eubalaena glacialis)) and grey whale (Eschrichtius robustus) specimens, two taxa no longer present in the eastern North Atlantic. This discovery matches expectations regarding the past utilization of North Atlantic right whales, but was unanticipated for grey whales, which have hitherto rarely been identified in the European zooarchaeological record. Many of these specimens derive from contexts associated with mediaeval cultures frequently linked to whaling: the Basques, northern Spaniards, Normans, Flemish, Frisians, Anglo-Saxons and Scandinavians. This association raises the likelihood that early whaling impacted these taxa, contributing to their extirpation and extinction. Much lower numbers of other large cetacean taxa were identified, suggesting that what are now the most depleted whales were once those most frequently used.

A quiet extirpation of the breeding little auk Alle alle population in Iceland in the shadow of the famous cousin extermination

The little auk Alle alle is an Arctic seabird breeding in the North Atlantic. Its southernmost breeding population in the Low Arctic occurred in Iceland but started to decline in numbers at the beginning of the 20th century ending in extirpation at the end of the same century. Climate warming has been blamed for the disappearance of this population. However, it was also exploited by humans (mainly for eggs). Thus, it is unclear what was the main driver for the population collapse. In this study, we reconstruct population dynamics in relation to changes in environmental conditions, and perform stochastic modelling of population viability considering various scenarios including presence/absence of climate warming and/or egg harvest. We found that extirpation of the studied population was attributed to synergistic effects of both climate warming and human harvest. The simulations revealed that climate warming without harvest would cause a 71% decline in population size but not lead to extirpation. Models with climate warming and egg harvesting resulted in population estimates close to the real data (decrease from 400 individuals in 1903 to 2 in 1996). This is one of the few studies documenting synergistic effect of climate warming and human exploitation on extirpation. A strong harvesting component in the explanation for the decline of the study population emphasizes the continuing need to control commercial harvest of animal species in the face of other pressures such as climate warming.

Ancient DNA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean

Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since at least the 16th century. Here we build upon previous work on ancient DNA of BFT in the Mediterranean by comparing contemporary (2009–2012) specimens with archival (1911–1926) and archaeological (2nd century BCE–15th century CE) specimens that represent population states prior to these two major periods of exploitation, respectively. We successfully genotyped and analysed 259 contemporary and 123 historical (91 archival and 32 archaeological) specimens at 92 SNP loci that were selected for their ability to differentiate contemporary populations or their association with core biological functions. We found no evidence of genetic bottlenecks, inbreeding or population restructuring between temporal sample groups that might explain what has driven catch fluctuations since the 16th century. We also detected a putative adaptive response, involving the cytoskeletal protein synemin which may be related to muscle stress. However, these results require further investigation with more extensive genome-wide data to rule out demographic changes due to overfishing, and other natural and anthropogenic factors, in addition to elucidating the adaptive drivers related to these.

Complex population structure of the Atlantic puffin revealed by whole genome analyses

The factors underlying gene flow and genomic population structure in vagile seabirds are notoriously difficult to understand due to their complex ecology with diverse dispersal barriers and extensive periods at sea. Yet, such understanding is vital for conservation management of seabirds that are globally declining at alarming rates. Here, we elucidate the population structure of the Atlantic puffin (Fratercula arctica) by assembling its reference genome and analyzing genome-wide resequencing data of 72 individuals from 12 colonies. We identify four large, genetically distinct clusters, observe isolation-by-distance between colonies within these clusters, and obtain evidence for a secondary contact zone. These observations disagree with the current taxonomy, and show that a complex set of contemporary biotic factors impede gene flow over different spatial scales. Our results highlight the power of whole genome data to reveal unexpected population structure in vagile marine seabirds and its value for seabird taxonomy, evolution and conservation.

Historical Demographic Processes Dominate Genetic Variation in Ancient Atlantic Cod Mitogenomes

Ancient DNA (aDNA) approaches have been successfully used to infer the long-term impacts of climate change, domestication, and human exploitation in a range of terrestrial species. Nonetheless, studies investigating such impacts using aDNA in marine species are rare. Atlantic cod (Gadus morhua), is an economically important species that has experienced dramatic census population declines during the last century. Here, we investigated 48 ancient mitogenomes from historical specimens obtained from a range of archeological excavations in northern Europe dated up to 6,500 BCE. We compare these mitogenomes to those of 496 modern conspecifics sampled across the North Atlantic Ocean and adjacent seas. Our results confirm earlier observations of high levels of mitogenomic variation and a lack of mutation-drift equilibrium—suggestive of population expansion. Furthermore, our temporal comparison yields no evidence of measurable mitogenomic changes through time. Instead, our results indicate that mitogenomic variation in Atlantic cod reflects past demographic processes driven by major historical events (such as oscillations in sea level) and subsequent gene flow rather than contemporary fluctuations in stock abundance. Our results indicate that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on mitogenomic diversity in a wide-spread marine species with high gene flow such as Atlantic cod. These observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and the importance of genetic variation in implementing conservation strategies. Instead, these observations imply that any measures toward the demographic recovery of Atlantic cod in the eastern Atlantic, will not be constrained by recent loss of historical mitogenomic variation.

Are Skyline Plot-Based Demographic Estimates Overly Dependent on Smoothing Prior Assumptions?

In Bayesian phylogenetics, the coalescent process provides an informative framework for inferring changes in the effective size of a population from a phylogeny (or tree) of sequences sampled from that population. Popular coalescent inference approaches such as the Bayesian Skyline Plot, Skyride, and Skygrid all model these population size changes with a discontinuous, piecewise-constant function but then apply a smoothing prior to ensure that their posterior population size estimates transition gradually with time. These prior distributions implicitly encode extra population size information that is not available from the observed coalescent data or tree. Here, we present a novel statistic, $\Omega$, to quantify and disaggregate the relative contributions of the coalescent data and prior assumptions to the resulting posterior estimate precision. Our statistic also measures the additional mutual information introduced by such priors. Using $\Omega$ we show that, because it is surprisingly easy to overparametrize piecewise-constant population models, common smoothing priors can lead to overconfident and potentially misleading inference, even under robust experimental designs. We propose $\Omega$ as a useful tool for detecting when effective population size estimates are overly reliant on prior assumptions and for improving quantification of the uncertainty in those estimates.[Coalescent processes; effective population size; information theory; phylodynamics; prior assumptions; skyline plots.].

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

Contamination from present-day DNA is a fundamental issue when studying ancient DNA from historical or archaeological material, and quantifying the amount of contamination is essential for downstream analyses. We present AuthentiCT, a command-line tool to estimate the proportion of present-day DNA contamination in ancient DNA datasets generated from single-stranded DNA libraries. The prediction is based solely on the patterns of post-mortem damage observed on ancient DNA sequences. The method has the power to quantify contamination from as few as 10,000 mapped sequences, making it particularly useful for analysing specimens that are poorly preserved or for which little data is available.