Response of an Afro-Palearctic bird migrant to glaciation cycles.

Migration allows animals to exploit spatially separated and seasonally available resources at a continental to global scale. However, responding to global climatic changes might prove challenging, especially for long-distance intercontinental migrants. During glacial periods, when conditions became too harsh for breeding in the north, avian migrants have been hypothesized to retract their distribution to reside within small refugial areas. Here, we present data showing that an Afro-Palearctic migrant continued seasonal migration, largely within Africa, during previous glacial-interglacial cycles with no obvious impact on population size. Using individual migratory track data to hindcast monthly bioclimatic habitat availability maps through the last 120,000 y, we show altered seasonal use of suitable areas through time. Independently derived effective population sizes indicate a growing population through the last 40,000 y. We conclude that the migratory lifestyle enabled adaptation to shifting climate conditions. This indicates that populations of resource-tracking, long-distance migratory species could expand successfully during warming periods in the past, which could also be the case under future climate scenarios.

Comparative genomics provides insights into the aquatic adaptations of mammals.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes.

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Genomic insights into the origin and diversification of late maritime hunter-gatherers from the Chilean Patagonia.

Patagonia was the last region of the Americas reached by humans who entered the continent from Siberia ∼15,000-20,000 y ago. Despite recent genomic approaches to reconstruct the continental evolutionary history, regional characterization of ancient and modern genomes remains understudied. Exploring the genomic diversity within Patagonia is not just a valuable strategy to gain a better understanding of the history and diversification of human populations in the southernmost tip of the Americas, but it would also improve the representation of Native American diversity in global databases of human variation. Here, we present genome data from four modern populations from Central Southern Chile and Patagonia (n = 61) and four ancient maritime individuals from Patagonia (∼1,000 y old). Both the modern and ancient individuals studied in this work have a greater genetic affinity with other modern Native Americans than to any non-American population, showing within South America a clear structure between major geographical regions. Native Patagonian Kawéskar and Yámana showed the highest genetic affinity with the ancient individuals, indicating genetic continuity in the region during the past 1,000 y before present, together with an important agreement between the ethnic affiliation and historical distribution of both groups. Lastly, the ancient maritime individuals were genetically equidistant to a ∼200-y-old terrestrial hunter-gatherer from Tierra del Fuego, which supports a model with an initial separation of a common ancestral group to both maritime populations from a terrestrial population, with a later diversification of the maritime groups.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea.

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of "gene vestiges" in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

Upper Palaeolithic Siberian genome reveals dual ancestry of Native Americans.

The origins of the First Americans remain contentious. Although Native Americans seem to be genetically most closely related to east Asians, there is no consensus with regard to which specific Old World populations they are closest to. Here we sequence the draft genome of an approximately 24,000-year-old individual (MA-1), from Mal'ta in south-central Siberia, to an average depth of 1×. To our knowledge this is the oldest anatomically modern human genome reported to date. The MA-1 mitochondrial genome belongs to haplogroup U, which has also been found at high frequency among Upper Palaeolithic and Mesolithic European hunter-gatherers, and the Y chromosome of MA-1 is basal to modern-day western Eurasians and near the root of most Native American lineages. Similarly, we find autosomal evidence that MA-1 is basal to modern-day western Eurasians and genetically closely related to modern-day Native Americans, with no close affinity to east Asians. This suggests that populations related to contemporary western Eurasians had a more north-easterly distribution 24,000 years ago than commonly thought. Furthermore, we estimate that 14 to 38% of Native American ancestry may originate through gene flow from this ancient population. This is likely to have occurred after the divergence of Native American ancestors from east Asian ancestors, but before the diversification of Native American populations in the New World. Gene flow from the MA-1 lineage into Native American ancestors could explain why several crania from the First Americans have been reported as bearing morphological characteristics that do not resemble those of east Asians. Sequencing of another south-central Siberian, Afontova Gora-2 dating to approximately 17,000 years ago, revealed similar autosomal genetic signatures as MA-1, suggesting that the region was continuously occupied by humans throughout the Last Glacial Maximum. Our findings reveal that western Eurasian genetic signatures in modern-day Native Americans derive not only from post-Columbian admixture, as commonly thought, but also from a mixed ancestry of the First Americans.

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.

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.

Identification of Polynesian mtDNA haplogroups in remains of Botocudo Amerindians from Brazil.

There is a consensus that modern humans arrived in the Americas 15,000-20,000 y ago during the Late Pleistocene, most probably from northeast Asia through Beringia. However, there is still debate about the time of entry and number of migratory waves, including apparent inconsistencies between genetic and morphological data on Paleoamericans. Here we report the identification of mitochondrial sequences belonging to haplogroups characteristic of Polynesians in DNA extracted from ancient skulls of the now extinct Botocudo Indians from Brazil. The identification of these two Polynesian haplogroups was confirmed in independent replications in Brazil and Denmark, ensuring reliability of the data. Parallel analysis of 12 other Botocudo individuals yielded only the well-known Amerindian mtDNA haplogroup C1. Potential scenarios to try to help understand these results are presented and discussed. The findings of this study may be relevant for the understanding of the pre-Columbian and/or post-Columbian peopling of the Americas.

Complete mitochondrial DNA sequence of the endangered giant sable antelope (Hippotragus niger variani): insights into conservation and taxonomy.

The giant sable antelope is one of the most endangered African bovids. Populations of this iconic animal, the national symbol of Angola, were recently rediscovered, after many decades of presumed extinction. Even so, their numbers are scarce and hence conservation plans are essential. However, fundamental information such as its taxonomic position, time of divergence and degree of genetic variation are still lacking. Here, we used a museum preserved horn as a source of DNA to describe, for the first time, the complete mitochondrial genome of the giant sable antelope, and provide insights into its evolutionary history. Reads generated by shotgun sequencing were mapped against the mitochondrial genome of common sable antelope and the nuclear genomes of cow and sheep. Phylogenetic reconstruction and divergence time estimate give support to the monophyly of the giant sable and a maximum divergence time of 170 thousand years to the closest subspecies. About 7% of the nuclear genome was mapped against the reference. The genetic resources reported here are now available for future work in the field of conservation genetics and phylogeny, in this and related species.

Interordinal gene capture, the phylogenetic position of Steller's sea cow based on molecular and morphological data, and the macroevolutionary history of Sirenia.

The recently extinct (ca. 1768) Steller's sea cow (Hydrodamalis gigas) was a large, edentulous North Pacific sirenian. The phylogenetic affinities of this taxon to other members of this clade, living and extinct, are uncertain based on previous morphological and molecular studies. We employed hybridization capture methods and second generation sequencing technology to obtain >30kb of exon sequences from 26 nuclear genes for both H. gigas and Dugong dugon. We also obtained complete coding sequences for the tooth-related enamelin (ENAM) gene. Hybridization probes designed using dugong and manatee sequences were both highly effective in retrieving sequences from H. gigas (mean=98.8% coverage), as were more divergent probes for regions of ENAM (99.0% coverage) that were designed exclusively from a proboscidean (African elephant) and a hyracoid (Cape hyrax). New sequences were combined with available sequences for representatives of all other afrotherian orders. We also expanded a previously published morphological matrix for living and fossil Sirenia by adding both new taxa and nine new postcranial characters. Maximum likelihood and parsimony analyses of the molecular data provide robust support for an association of H. gigas and D. dugon to the exclusion of living trichechids (manatees). Parsimony analyses of the morphological data also support the inclusion of H. gigas in Dugongidae with D. dugon and fossil dugongids. Timetree analyses based on calibration density approaches with hard- and soft-bounded constraints suggest that H. gigas and D. dugon diverged in the Oligocene and that crown sirenians last shared a common ancestor in the Eocene. The coding sequence for the ENAM gene in H. gigas does not contain frameshift mutations or stop codons, but there is a transversion mutation (AG to CG) in the acceptor splice site of intron 2. This disruption in the edentulous Steller's sea cow is consistent with previous studies that have documented inactivating mutations in tooth-specific loci of a variety of edentulous and enamelless vertebrates including birds, turtles, aardvarks, pangolins, xenarthrans, and baleen whales. Further, branch-site dN/dS analyses provide evidence for positive selection in ENAM on the stem dugongid branch where extensive tooth reduction occurred, followed by neutral evolution on the Hydrodamalis branch. Finally, we present a synthetic evolutionary tree for living and fossil sirenians showing several key innovations in the history of this clade including character state changes that parallel those that occurred in the evolutionary history of cetaceans.

Genetic evidence for patrilocal mating behavior among Neandertal groups.

The remains of 12 Neandertal individuals have been found at the El Sidrón site (Asturias, Spain), consisting of six adults, three adolescents, two juveniles, and one infant. Archaeological, paleontological, and geological evidence indicates that these individuals represent all or part of a contemporaneous social group of Neandertals, who died at around the same time and later were buried together as a result of a collapse of an underground karst. We sequenced phylogenetically informative positions of mtDNA hypervariable regions 1 and 2 from each of the remains. Our results show that the 12 individuals stem from three different maternal lineages, accounting for seven, four, and one individual(s), respectively. Using a Y-chromosome assay to confirm the morphological determination of sex for each individual, we found that, although the three adult males carried the same mtDNA lineage, each of the three adult females carried different mtDNA lineages. These findings provide evidence to indicate that Neandertal groups not only were small and characterized by low genetic diversity but also were likely to have practiced patrilocal mating behavior.

Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia.

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Mitochondrial genome diversity and population structure of the giant squid Architeuthis: genetics sheds new light on one of the most enigmatic marine species.

Despite its charismatic appeal to both scientists and the general public, remarkably little is known about the giant squid Architeuthis, one of the largest of the invertebrates. Although specimens of Architeuthis are becoming more readily available owing to the advancement of deep-sea fishing techniques, considerable controversy exists with regard to topics as varied as their taxonomy, biology and even behaviour. In this study, we have characterized the mitochondrial genome (mitogenome) diversity of 43 Architeuthis samples collected from across the range of the species, in order to use genetic information to provide new and otherwise difficult to obtain insights into the life of this animal. The results show no detectable phylogenetic structure at the mitochondrial level and, furthermore, that the level of nucleotide diversity is exceptionally low. These observations are consistent with the hypotheses that there is only one global species of giant squid, Architeuthis dux (Steenstrup, 1857), and that it is highly vagile, possibly dispersing through both a drifting paralarval stage and migration of larger individuals. Demographic history analyses of the genetic data suggest that there has been a recent population expansion or selective sweep, which may explain the low level of genetic diversity.

Ancient DNA analyses exclude humans as the driving force behind late Pleistocene musk ox (Ovibos moschatus) population dynamics.

The causes of the late Pleistocene megafaunal extinctions are poorly understood. Different lines of evidence point to climate change, the arrival of humans, or a combination of these events as the trigger. Although many species went extinct, others, such as caribou and bison, survived to the present. The musk ox has an intermediate story: relatively abundant during the Pleistocene, it is now restricted to Greenland and the Arctic Archipelago. In this study, we use ancient DNA sequences, temporally unbiased summary statistics, and Bayesian analytical techniques to infer musk ox population dynamics throughout the late Pleistocene and Holocene. Our results reveal that musk ox genetic diversity was much higher during the Pleistocene than at present, and has undergone several expansions and contractions over the past 60,000 years. Northeast Siberia was of key importance, as it was the geographic origin of all samples studied and held a large diverse population until local extinction at approximately 45,000 radiocarbon years before present ((14)C YBP). Subsequently, musk ox genetic diversity reincreased at ca. 30,000 (14)C YBP, recontracted at ca. 18,000 (14)C YBP, and finally recovered in the middle Holocene. The arrival of humans into relevant areas of the musk ox range did not affect their mitochondrial diversity, and both musk ox and humans expanded into Greenland concomitantly. Thus, their population dynamics are better explained by a nonanthropogenic cause (for example, environmental change), a hypothesis supported by historic observations on the sensitivity of the species to both climatic warming and fluctuations.

Assessing microplastic exposure of the Critically Endangered Mediterranean monk seal (Monachus monachus) on a remote oceanic island.

Microplastics (<5 mm) are a ubiquitous marine pollutant which are highly bioavailable to marine organisms across all trophic levels. Marine predators are especially vulnerable to microplastic pollution through direct and indirect ingestion (e.g., trophic transfer) due to their high trophic position. In particular, oceanic islands are more susceptible to plastic accumulation, increasing the relative number of microplastics in the environment that are available for consumption. The dynamics of microplastic accumulation in marine predators inhabiting remote islands, however, is sparsely documented. Here we describe microplastic exposure in the Critically Endangered Mediterranean monk seal (Monachus monachus) from the Madeira archipelago (Northeast Atlantic) using scat-based analysis. Microplastics were recovered from 18 scat samples collected between 2014-2021 and were characterized to the polymer level using Fourier-Transform Infrared (u-FTIR) spectroscopy. A total of 390 microplastic particles were recovered, ranging between 0.2-8.6 particles g-1 dry weight (mean 1.84 ± 2.14 particles g-1) consisting mainly of fragments (69 %) of various sizes and polymer composition (e.g., PE, PET, PS). Microplastic prevalence (100 % of samples analysed) was higher than what has been previously recorded using scat-based analysis in other pinniped species. Our results suggest that the levels of microplastic pollution in the coastal food-web in the Madeira archipelago are relatively high, placing higher-trophic level organisms at increased risk of microplastic consumption, including humans. This study provides the first insights into microplastic exposure to Madeira's monk seals that may contribute to future management decisions for the species and their long-term survival.

The half-life of DNA in bone: measuring decay kinetics in 158 dated fossils.

Claims of extreme survival of DNA have emphasized the need for reliable models of DNA degradation through time. By analysing mitochondrial DNA (mtDNA) from 158 radiocarbon-dated bones of the extinct New Zealand moa, we confirm empirically a long-hypothesized exponential decay relationship. The average DNA half-life within this geographically constrained fossil assemblage was estimated to be 521 years for a 242 bp mtDNA sequence, corresponding to a per nucleotide fragmentation rate (k) of 5.50 × 10(-6) per year. With an effective burial temperature of 13.1°C, the rate is almost 400 times slower than predicted from published kinetic data of in vitro DNA depurination at pH 5. Although best described by an exponential model (R(2) = 0.39), considerable sample-to-sample variance in DNA preservation could not be accounted for by geologic age. This variation likely derives from differences in taphonomy and bone diagenesis, which have confounded previous, less spatially constrained attempts to study DNA decay kinetics. Lastly, by calculating DNA fragmentation rates on Illumina HiSeq data, we show that nuclear DNA has degraded at least twice as fast as mtDNA. These results provide a baseline for predicting long-term DNA survival in bone.

The population genomic legacy of the second plague pandemic.

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes.

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.