Harnessing ancient genomes to study the history of human adaptation

Assessing the impact of post-mortem damage and contamination on imputation performance in ancient DNA

Low-coverage imputation is becoming ever more present in ancient DNA (aDNA) studies. Imputation pipelines commonly used for present-day genomes have been shown to yield accurate results when applied to ancient genomes. However, post-mortem damage (PMD), in the form of C-to-T substitutions at the reads termini, and contamination with DNA from closely related species can potentially affect imputation performance in aDNA. In this study, we evaluated imputation performance (i) when using a genotype caller designed for aDNA, ATLAS, compared to bcftools, and (ii) when contamination is present. We evaluated imputation performance with principal component analyses and by calculating imputation error rates. With a particular focus on differently imputed sites, we found that using ATLAS prior to imputation substantially improved imputed genotypes for a very damaged ancient genome (42% PMD). Trimming the ends of the sequencing reads led to similar improvements in imputation accuracy. For the remaining genomes, ATLAS brought limited gains. Finally, to examine the effect of contamination on imputation, we added various amounts of reads from two present-day genomes to a previously downsampled high-coverage ancient genome. We observed that imputation accuracy drastically decreased for contamination rates above 5%. In conclusion, we recommend (i) accounting for PMD by either trimming sequencing reads or using a genotype caller such as ATLAS before imputing highly damaged genomes and (ii) only imputing genomes containing up to 5% of contamination.

Predicting Functional Consequences of Recent Natural Selection in Britain

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are noncoding, so we expect that a large proportion of the phenotypic effects of selection will also act through noncoding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation [JTI]) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4,500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (false discovery rate [FDR] < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-single nucleotide polymorphism (SNP) selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally, we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection.

Genetic history of Cambridgeshire before and after the Black Death

The extent of the devastation of the Black Death pandemic (1346-1353) on European populations is known from documentary sources and its bacterial source illuminated by studies of ancient pathogen DNA. What has remained less understood is the effect of the pandemic on human mobility and genetic diversity at the local scale. Here, we report 275 ancient genomes, including 109 with coverage >0.1×, from later medieval and postmedieval Cambridgeshire of individuals buried before and after the Black Death. Consistent with the function of the institutions, we found a lack of close relatives among the friars and the inmates of the hospital in contrast to their abundance in general urban and rural parish communities. While we detect long-term shifts in local genetic ancestry in Cambridgeshire, we find no evidence of major changes in genetic ancestry nor higher differentiation of immune loci between cohorts living before and after the Black Death.

Unveiling recent and ongoing adaptive selection in human populations

Genome-wide scans for signals of selection have become a routine part of the analysis of population genomic variation datasets and have resulted in compelling evidence of selection during recent human evolution. This Essay spotlights methodological innovations that have enabled the detection of selection over very recent timescales, even in contemporary human populations. By harnessing large-scale genomic and phenotypic datasets, these new methods use different strategies to uncover connections between genotype, phenotype, and fitness. This Essay outlines the rationale and key findings of each strategy, discusses challenges in interpretation, and describes opportunities to improve detection and understanding of ongoing selection in human populations.

Allelic bias when performing in-solution enrichment of ancient human DNA

In-solution hybridisation enrichment of genetic variation is a valuable methodology in human paleogenomics. It allows enrichment of endogenous DNA by targeting genetic markers that are comparable between sequencing libraries. Many studies have used the 1240k reagent-which enriches 1,237,207 genome-wide SNPs-since 2015, though access was restricted. In 2021, Twist Biosciences and Daicel Arbor Biosciences independently released commercial kits that enabled all researchers to perform enrichments for the same 1240 k SNPs. We used the Daicel Arbor Biosciences Prime Plus kit to enrich 132 ancient samples from three continents. We identified a systematic assay bias that increases genetic similarity between enriched samples and that cannot be explained by batch effects. We present the impact of the bias on population genetics inferences (e.g. Principal Components Analysis, ƒ-statistics) and genetic relatedness (READ). We compare the Prime Plus bias to that previously reported of the legacy 1240k enrichment assay. In ƒ-statistics, we find that all Prime-Plus-generated data exhibit artefactual excess shared drift, such that within-continent relationships cannot be correctly determined. The bias is more subtle in READ, though interpretation of the results can still be misleading in specific contexts. We expect the bias may affect analyses we have not yet tested. Our observations support previously reported concerns for the integration of different data types in paleogenomics. We also caution that technological solutions to generate 1240k data necessitate a thorough validation process before their adoption in the paleogenomic community.

Predicting functional consequences of recent natural selection in Britain

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are non-coding, so we expect that a large proportion of the phenotypic effects of selection will also act through non-coding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation; JTI) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (FDR < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-SNP selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection.

Going local with ancient DNA: A review of human histories from regional perspectives

Ancient DNA (aDNA) has added a wealth of information about our species' history, including insights on genetic origins, migrations and gene flow, genetic admixture, and health and disease. Much early work has focused on continental-level questions, leaving many regional questions, especially those relevant to the Global South, comparatively underexplored. A few success stories of aDNA studies from smaller laboratories involve more local aspects of human histories and health in the Americas, Africa, Asia, and Oceania. In this Review, we cover some of these contributions by synthesizing finer-scale questions of importance to the archaeogenetics field, as well as to Indigenous and Descendant communities. We further highlight the potential of aDNA to uncover past histories in regions where colonialism has neglected the oral histories of oppressed peoples.

Methods for Assessing Population Relationships and History Using Genomic Data

Genetic data contain a record of our evolutionary history. The availability of large-scale datasets of human populations from various geographic areas and timescales, coupled with advances in the computational methods to analyze these data, has transformed our ability to use genetic data to learn about our evolutionary past. Here, we review some of the widely used statistical methods to explore and characterize population relationships and history using genomic data. We describe the intuition behind commonly used approaches, their interpretation, and important limitations. For illustration, we apply some of these techniques to genome-wide autosomal data from 929 individuals representing 53 worldwide populations that are part of the Human Genome Diversity Project. Finally, we discuss the new frontiers in genomic methods to learn about population history. In sum, this review highlights the power (and limitations) of DNA to infer features of human evolutionary history, complementing the knowledge gleaned from other disciplines, such as archaeology, anthropology, and linguistics.

Ancient dolphin genomes reveal rapid repeated adaptation to coastal waters

Parallel evolution provides strong evidence of adaptation by natural selection due to local environmental variation. Yet, the chronology, and mode of the process of parallel evolution remains debated. Here, we harness the temporal resolution of paleogenomics to address these long-standing questions, by comparing genomes originating from the mid-Holocene (8610-5626 years before present, BP) to contemporary pairs of coastal-pelagic ecotypes of bottlenose dolphin. We find that the affinity of ancient samples to coastal populations increases as the age of the samples decreases. We assess the youngest genome (5626 years BP) at sites previously inferred to be under parallel selection to coastal habitats and find it contained coastal-associated genotypes. Thus, coastal-associated variants rose to detectable frequencies close to the emergence of coastal habitat. Admixture graph analyses reveal a reticulate evolutionary history between pelagic and coastal populations, sharing standing genetic variation that facilitated rapid adaptation to newly emerged coastal habitats.

Nutrigenomics in the context of evolution

Nutrigenomics describes the interaction between nutrients and our genome. Since the origin of our species most of these nutrient-gene communication pathways have not changed. However, our genome experienced over the past 50,000 years a number of evolutionary pressures, which are based on the migration to new environments concerning geography and climate, the transition from hunter-gatherers to farmers including the zoonotic transfer of many pathogenic microbes and the rather recent change of societies to a preferentially sedentary lifestyle and the dominance of Western diet. Human populations responded to these challenges not only by specific anthropometric adaptations, such as skin color and body stature, but also through diversity in dietary intake and different resistance to complex diseases like the metabolic syndrome, cancer and immune disorders. The genetic basis of this adaptation process has been investigated by whole genome genotyping and sequencing including that of DNA extracted from ancient bones. In addition to genomic changes, also the programming of epigenomes in pre- and postnatal phases of life has an important contribution to the response to environmental changes. Thus, insight into the variation of our (epi)genome in the context of our individual's risk for developing complex diseases, helps to understand the evolutionary basis how and why we become ill. This review will discuss the relation of diet, modern environment and our (epi)genome including aspects of redox biology. This has numerous implications for the interpretation of the risks for disease and their prevention.

Modeling the spatiotemporal spread of beneficial alleles using ancient genomes

Ancient genome sequencing technologies now provide the opportunity to study natural selection in unprecedented detail. Rather than making inferences from indirect footprints left by selection in present-day genomes, we can directly observe whether a given allele was present or absent in a particular region of the world at almost any period of human history within the last 10,000 years. Methods for studying selection using ancient genomes often rely on partitioning individuals into discrete time periods or regions of the world. However, a complete understanding of natural selection requires more nuanced statistical methods which can explicitly model allele frequency changes in a continuum across space and time. Here we introduce a method for inferring the spread of a beneficial allele across a landscape using two-dimensional partial differential equations. Unlike previous approaches, our framework can handle time-stamped ancient samples, as well as genotype likelihoods and pseudohaploid sequences from low-coverage genomes. We apply the method to a panel of published ancient West Eurasian genomes to produce dynamic maps showcasing the inferred spread of candidate beneficial alleles over time and space. We also provide estimates for the strength of selection and diffusion rate for each of these alleles. Finally, we highlight possible avenues of improvement for accurately tracing the spread of beneficial alleles in more complex scenarios.

CONGA: Copy number variation genotyping in ancient genomes and low-coverage sequencing data

To date, ancient genome analyses have been largely confined to the study of single nucleotide polymorphisms (SNPs). Copy number variants (CNVs) are a major contributor of disease and of evolutionary adaptation, but identifying CNVs in ancient shotgun-sequenced genomes is hampered by typical low genome coverage (<1×) and short fragments (<80 bps), precluding standard CNV detection software to be effectively applied to ancient genomes. Here we present CONGA, tailored for genotyping CNVs at low coverage. Simulations and down-sampling experiments suggest that CONGA can genotype deletions >1 kbps with F-scores >0.75 at ≥1×, and distinguish between heterozygous and homozygous states. We used CONGA to genotype 10,002 outgroup-ascertained deletions across a heterogenous set of 71 ancient human genomes spanning the last 50,000 years, produced using variable experimental protocols. A fraction of these (21/71) display divergent deletion profiles unrelated to their population origin, but attributable to technical factors such as coverage and read length. The majority of the sample (50/71), despite originating from nine different laboratories and having coverages ranging from 0.44×-26× (median 4×) and average read lengths 52-121 bps (median 69), exhibit coherent deletion frequencies. Across these 50 genomes, inter-individual genetic diversity measured using SNPs and CONGA-genotyped deletions are highly correlated. CONGA-genotyped deletions also display purifying selection signatures, as expected. CONGA thus paves the way for systematic CNV analyses in ancient genomes, despite the technical challenges posed by low and variable genome coverage.

Parallel evolution in human populations: A biocultural perspective

Parallel evolution-where different populations evolve similar traits in response to similar environments-has been a topic of growing interest to biologists and biological anthropologists for decades. Parallel evolution occurs in human populations thanks to myriad biological and cultural mechanisms that permit humans to survive and thrive in diverse environments worldwide. Because humans shape and are shaped by their environments, biocultural approaches that emphasize the interconnections between biology and culture are key to understanding parallel evolution in human populations as well as the nuances of human biological variation and adaptation. In this review, we discuss how biocultural theory has been and can be applied to studies of parallel evolution and adaptation more broadly. We illustrate this through four examples of parallel evolution in humans: malaria resistance, lactase persistence, cold tolerance, and high-altitude adaptation.

Achieving improved accuracy for imputation of ancient DNA

MOTIVATION

Genotype imputation has the potential to increase the amount of information that can be gained from the often limited biological material available in ancient samples. As many widely used tools have been developed with modern data in mind, their design is not necessarily reflective of the requirements in studies of ancient DNA. Here, we investigate if an imputation method based on the full probabilistic Li and Stephens model of haplotype frequencies might be beneficial for the particular challenges posed by ancient data.

RESULTS

We present an implementation called prophaser and compare imputation performance to two alternative pipelines that have been used in the ancient DNA community based on the Beagle software. Considering empirical ancient data downsampled to lower coverages as well as present-day samples with artificially thinned genotypes, we show that the proposed method is advantageous at lower coverages, where it yields improved accuracy and ability to capture rare variation. The software prophaser is optimized for running in a massively parallel manner and achieved reasonable runtimes on the experiments performed when executed on a GPU.

AVAILABILITY AND IMPLEMENTATION

The C++ code for prophaser is available in the GitHub repository https://github.com/scicompuu/prophaser.

SUPPLEMENTARY INFORMATION

Supplementary information is available at Bioinformatics online.

A Late Pleistocene human genome from Southwest China

Southern East Asia is the dispersal center regarding the prehistoric settlement and migrations of modern humans in Asia-Pacific regions. However, the settlement pattern and population structure of paleolithic humans in this region remain elusive, and ancient DNA can provide direct information. Here, we sequenced the genome of a Late Pleistocene hominin (MZR), dated ∼14.0 thousand years ago from Red Deer Cave located in Southwest China, which was previously reported possessing mosaic features of modern and archaic hominins. MZR is the first Late Pleistocene genome from southern East Asia. Our results indicate that MZR is a modern human who represents an early diversified lineage in East Asia. The mtDNA of MZR belongs to an extinct basal lineage of the M9 haplogroup, reflecting a rich matrilineal diversity in southern East Asia during the Late Pleistocene. Combined with the published data, we detected clear genetic stratification in ancient southern populations of East/Southeast Asia and some degree of south-versus-north divergency during the Late Pleistocene, and MZR was identified as a southern East Asian who exhibits genetic continuity to present day populations. Markedly, MZR is linked deeply to the East Asian ancestry that contributed to First Americans.

The untapped potential of macrofossils in ancient plant DNA research

The rapid development of ancient DNA analysis in the last decades has induced a paradigm shift in ecology and evolution. Driven by a combination of breakthroughs in DNA isolation techniques, high-throughput sequencing, and bioinformatics, ancient genome-scale data for a rapidly growing variety of taxa are now available, allowing researchers to directly observe demographic and evolutionary processes over time. However, the vast majority of paleogenomic studies still focus on human or animal remains. In this article, we make the case for a vast untapped resource of ancient plant material that is ideally suited for paleogenomic analyses: plant remains, such as needles, leaves, wood, seeds, or fruits, that are deposited in natural archives, such as lake sediments, permafrost, or even ice caves. Such plant remains are commonly found in large numbers and in stratigraphic sequence through time and have so far been used primarily to reconstruct past local species presences and abundances. However, they are also unique repositories of genetic information with the potential to revolutionize the fields of ecology and evolution by directly studying microevolutionary processes over time. Here, we give an overview of the current state-of-the-art, address important challenges, and highlight new research avenues to inspire future research.

Ergothioneine: A Stress Vitamin with Antiaging, Vascular, and Neuroprotective Roles?

Significance: Ergothioneine (ET) is an unusual sulfur-containing amino acid derived from histidine, acquired predominantly from food. Its depletion is associated with deleterious consequences in response to stress stimuli in cell culture models, prompting us to classify it as a vitamin in 2010, which was later supported by in vivo studies. ET is obtained from a variety of foods and is taken up by a selective transporter. ET possesses antioxidant and anti-inflammatory properties that confer cytoprotection. ET crosses the blood-brain barrier and has been reported to have beneficial effects in the brain. In this study, we discuss the cytoprotective and neuroprotective properties of ET, which may be harnessed for combating neurodegeneration and decline during aging. Recent Advances: The designation of ET as a stress vitamin is gaining momentum, opening a new field of investigation involving small molecules that are essential for optimal physiological functioning and maintenance of health span. Critical Issues: Although ET was discovered more than a century ago, its physiological functions are still being elucidated, especially in the brain. As ET is present in most foods, toxicity associated with its deprivation has been difficult to assess. Future Directions: Using genetically engineered cells and mice, it may now be possible to elucidate roles of ET. This coupled with advances in genomics and metabolomics may lead to identification of ET function. As ET is a stable antioxidant with anti-inflammatory properties, whose levels decline during aging, supplementing ET in the diet or consuming an ET-rich diet may prove beneficial. Antioxid. Redox Signal. 36, 1306-1317.

Ethical Guidance in Human Paleogenomics: New and Ongoing Perspectives

Over the past two decades, the study of ancient genomes from Ancestral humans, or human paleogenomic research, has expanded rapidly in both scale and scope. Ethical discourse has subsequently emerged to address issues of social responsibility and scientific robusticity in conducting research. Here, we highlight and contextualize the primary sources of professional ethical guidance aimed at paleogenomic researchers. We describe the tension among existing guidelines, while addressing core issues such as consent, destructive research methods, and data access and management. Currently, there is a dissonance between guidelines that focus on scientific outcomes and those that hold scientists accountable to stakeholder communities, such as descendants. Thus, we provide additional tools to navigate the complexities of ancient DNA research while centering engagement with stakeholder communities in the scientific process. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Performance of innovative nanomaterials for bone remains consolidation and effect on 14C dating and on palaeogenetic analysis

An innovative protocol for the consolidation of ancient bone remains based on the use of nanometric HydroxyAPatite (HAP) was set up and tested through a multidisciplinary approach. A new protocol for the synthesis of HAP nanoparticles was developed, and the composition of the obtained nanomaterial was investigated through Fourier Transform Infrared Spectroscopy (FTIR) and X-Ray Diffraction (XRD); sizes, shape and morphology of the synthesized particles were studied by Scanning Electron Microscopy (SEM). The consolidation performance was evaluated by testing the new nanomaterial on degraded ancient bone findings. An increase of the mineral density and of the micro-hardness of the bone were observed. The new consolidation method was also tested to assess possible effects on the palaeogenetic analysis and radiocarbon dating on the treated bones. The consolidation treatment does not introduce any contaminations that could affect radiocarbon dating and has no general detrimental impact on the genetic characterization of the skeletal remains. This consolidation procedure represents a more compatible conservation tool with respect to traditional procedures: it has been shown that the treatment is effective, easily-applicable and compatible with post-consolidation analysis.

Predicting Archaic Hominin Phenotypes from Genomic Data

Ancient DNA provides a powerful window into the biology of extant and extinct species, including humans' closest relatives: Denisovans and Neanderthals. Here, we review what is known about archaic hominin phenotypes from genomic data and how those inferences have been made. We contend that understanding the influence of variants on lower-level molecular phenotypes-such as gene expression and protein function-is a promising approach to using ancient DNA to learn about archaic hominin traits. Molecular phenotypes have simpler genetic architectures than organism-level complex phenotypes, and this approach enables moving beyond association studies by proposing hypotheses about the effects of archaic variants that are testable in model systems. The major challenge to understanding archaic hominin phenotypes is broadening our ability to accurately map genotypes to phenotypes, but ongoing advances ensure that there will be much more to learn about archaic hominin phenotypes from their genomes. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A geometric relationship of F 2 , F 3 and F 4 -statistics with principal component analysis

Principal component analysis (PCA) and F -statistics sensu Patterson are two of the most widely used population genetic tools to study human genetic variation. Here, I derive explicit connections between the two approaches and show that these two methods are closely related. F -statistics have a simple geometrical interpretation in the context of PCA, and orthogonal projections are a key concept to establish this link. I show that for any pair of populations, any population that is admixed as determined by an F 3 -statistic will lie inside a circle on a PCA plot. Furthermore, the F 4 -statistic is closely related to an angle measurement, and will be zero if the differences between pairs of populations intersect at a right angle in PCA space. I illustrate my results on two examples, one of Western Eurasian, and one of global human diversity. In both examples, I find that the first few PCs are sufficient to approximate most F -statistics, and that PCA plots are effective at predicting F -statistics. Thus, while F -statistics are commonly understood in terms of discrete populations, the geometric perspective illustrates that they can be viewed in a framework of populations that vary in a more continuous manner.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Tracing the Evolution of Human Gene Regulation and Its Association with Shifts in Environment

As humans populated the world, they adapted to many varying environmental factors, including climate, diet, and pathogens. Because many of these adaptations were mediated by multiple noncoding variants with small effects on gene regulation, it has been difficult to link genomic signals of selection to specific genes, and to describe the regulatory response to selection. To overcome this challenge, we adapted PrediXcan, a machine learning method for imputing gene regulation from genotype data, to analyze low-coverage ancient human DNA (aDNA). First, we used simulated genomes to benchmark strategies for adapting PrediXcan to increase robustness to incomplete data. Applying the resulting models to 490 ancient Eurasians, we found that genes with the strongest divergent regulation among ancient populations with hunter-gatherer, pastoralist, and agricultural lifestyles are enriched for metabolic and immune functions. Next, we explored the contribution of divergent gene regulation to two traits with strong evidence of recent adaptation: dietary metabolism and skin pigmentation. We found enrichment for divergent regulation among genes proposed to be involved in diet-related local adaptation, and the predicted effects on regulation often suggest explanations for known signals of selection, for example, at FADS1, GPX1, and LEPR. In contrast, skin pigmentation genes show little regulatory change over a 38,000-year time series of 2,999 ancient Europeans, suggesting that adaptation mainly involved large-effect coding variants. This work demonstrates that combining aDNA with present-day genomes is informative about the biological differences among ancient populations, the role of gene regulation in adaptation, and the relationship between genetic diversity and complex traits.

A Dual Role for Behavior in Evolution and Shaping Organismal Selective Environments

The hypothesis that evolved behaviors play a determining role in facilitating and impeding the evolution of other traits has been discussed for more than 100 years with little consensus beyond an agreement that the ideas are theoretically plausible in accord with the Modern Synthesis. Many recent reviews of the genomic, epigenetic, and developmental mechanisms underpinning major behavioral transitions show how facultative expression of novel behaviors can lead to the evolution of obligate behaviors and structures that enhance behavioral function. Phylogenetic and genomic studies indicate that behavioral traits are generally evolutionarily more labile than other traits and that they help shape selective environments on the latter traits. Adaptive decision-making to encounter resources and avoid stress sources requires specific sensory inputs, which behaviorally shape selective environments by determining those features of the external world that are biologically relevant. These recent findings support the hypothesis of a dual role for behavior in evolution and are consistent with current evolutionary theory.

Similarity-based analysis of allele frequency distribution among multiple populations identifies adaptive genomic structural variants

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Spatially explicit paleogenomic simulations support cohabitation with limited admixture between Bronze Age Central European populations

The Bronze Age is a complex period of social, cultural and economic changes. Recent paleogenomic studies have documented a large and rapid genetic change in early Bronze Age populations from Central Europe. However, the detailed demographic and genetic processes involved in this change are still debated. Here we have used spatially explicit simulations of genomic components to better characterize the demographic and migratory conditions that may have led to this change. We investigated various scenarios representing the expansion of pastoralists from the Pontic steppe, potentially linked to the Yamnaya cultural complex, and their interactions with local populations in Central Europe, considering various eco-evolutionary factors, such as population admixture, competition and long-distance dispersal. Our results do not support direct competition but rather the cohabitation of pastoralists and farmers in Central Europe, with limited gene flow between populations. They also suggest occasional long-distance migrations accompanying the expansion of pastoralists and a demographic decline in both populations following their initial contact. These results link recent archaeological and paleogenomic observations and move further the debate of genomic changes during the early Bronze Age.

COMBINING ANCIENT DNA AND RADIOCARBON DATING DATA TO INCREASE CHRONOLOGICAL ACCURACY

This paper examines how ancient DNA data can enhance radiocarbon dating. Because there is a limit to the number of years that can separate the dates of death of related individuals, the ability to identify relatives through ancient DNA analysis can serve as a constraint on radiocarbon date range estimates. To determine the number of years that can separate related individuals, we modeled maximums derived from biological extremes of human reproduction and death ages and compiled data from historic and genealogical death records. We used these data to jointly study the date ranges of a global dataset of individuals that have been radiocarbon dated and for which ancient DNA analysis identified at least one relative. We found that many of these individuals could have their date uncertainties reduced by building in date of death separation constraints. We examined possible reasons for date discrepancies of related individuals, such as dating of different skeletal elements or wiggles in the radiocarbon curve. We also developed a program, refinedate, which researchers can download and use to help refine the radiocarbon date distributions of related individuals. Our research demonstrates that when combined, radiocarbon dating and ancient DNA analysis can provide a refined and richer view of the past.

Paleogenomics of the prehistory of Europe: human migrations, domestication and disease

A substantial portion of ancient DNA research has been centred on understanding European populations' origin and evolution. A rchaeological evidence has already shown that the peopling of Europe involved an intricate pattern of demic and/or cultural diffusion since the Upper Palaeolithic, which became more evident during the Neolithic and Bronze Age periods. However, ancient DNA data has been crucial in determining if cultural changes occurred due to the movement of ideas or people. With the advent of next-generation sequencing and population-based paleogenomic research, ancient DNA studies have been directed not only at the study of continental human migrations, but also to the detailed analysis of particular archaeological sites, the processes of domestication, or the spread of disease during prehistoric times. With this vast paleogenomic effort added to a proper archaeological contextualisation of results, a deeper understanding of Europe's peopling is starting to emanate.

Mitochondrial genetic variation in human bioenergetics, adaptation, and adult disease

OBJECTIVES

Mitochondria are critical for the survival of eukaryotic organisms due to their ability to produce cellular energy, which drives virtually all aspects of host biology. However, the effects of mitochondrial DNA (mtDNA) variation in relation to disease etiology and adaptation within contemporary global human populations remains incompletely understood.

METHODS

To develop a more holistic understanding of the role of mtDNA diversity in human adaptation, health, and disease, we investigated mitochondrial biology and bioenergetics. More specifically, we synthesized details from studies of mitochondrial function and variation in the context of haplogroup background, climatic adaptation, and oxidative disease.

RESULTS

The majority of studies show that mtDNA variation arose during modern human dispersal around the world. Some of these variants appear to have been positively selected for their adaptiveness in colder climates, with these sequence changes having implications for tissue-specific function and thermogenic capacity. In addition, many variants modulating energy production are also associated with damaging metabolic byproducts and mitochondrial dysfunction, which, in turn, are implicated in the onset and severity of several different adult mitochondrial diseases. Thus, mtDNA variation that governs bioenergetics, metabolism, and thermoregulation may potentially have adverse consequences for human health, depending on the genetic background and context in which it occurs.

CONCLUSIONS

Our review suggests that the mitochondrial research field would benefit from independently replicating mtDNA haplogroup-phenotype associations across global populations, incorporating potentially confounding environmental, demographic, and disease covariates into studies of mtDNA variation, and extending association-based studies to include analyses of complete mitogenomes and assays of mitochondrial function.

Ancient DNA analysis of rare genetic bone disorders

OBJECTIVE

Review of the current advancements in the field of paleogenetics that provide new opportunities in studying the evolution of rare genetic bone diseases.

MATERIAL AND METHODS

Based on cases from the literature, the genetics of rare bone diseases will be introduced and the main methodological issues will be addressed, focusing on the opportunities presented by the application of aDNA analyses in the field of paleopathology.

RESULTS

Medical literature provides large datasets on the genes responsible for rare bone disorders. These genes, subdivided in functional categories, display important future targets when analyzing rare genetic bone disorders in ancient human remains.

CONCLUSIONS

Knowledge on both phenotype and genotype is required to study rare diseases in ancient human remains.

SIGNIFICANCE

The proposed interdisciplinary research will provide new insight into the occurrence and spread of genetic risk factors in the past and will help in the diagnostics of these rare and often neglected diseases.

LIMITATIONS

The current limitations in ancient DNA research and targeting the disease-causing specific mutations (e.g., somatic or germline).

SUGGESTIONS FOR FURTHER RESEARCH

Methodological advancements and candidate gene lists provide the optimal basis for future interdisciplinary studies of rare genetic bone disorders in ancient human remains.

A global carbon and nitrogen isotope perspective on modern and ancient human diet

Stable carbon and nitrogen isotope analyses are widely used to infer diet and mobility in ancient and modern human populations, potentially providing a means to situate humans in global food webs. We collated 13,666 globally distributed analyses of ancient and modern human collagen and keratin samples. We converted all data to a common "Modern Diet Equivalent" reference frame to enable direct comparison among modern human diets, human diets prior to the advent of industrial agriculture, and the natural environment. This approach reveals a broad diet prior to industrialized agriculture and continued in modern subsistence populations, consistent with the human ability to consume opportunistically as extreme omnivores within complex natural food webs and across multiple trophic levels in every terrestrial and many marine ecosystems on the planet. In stark contrast, isotope dietary breadth across modern nonsubsistence populations has compressed by two-thirds as a result of the rise of industrialized agriculture and animal husbandry practices and the globalization of food distribution networks.

Inferring archaic introgression from hominin genetic data

Questions surrounding the timing, extent, and evolutionary consequences of archaic admixture into human populations have a long history in evolutionary anthropology. More recently, advances in human genetics, particularly in the field of ancient DNA, have shed new light on the question of whether or not Homo sapiens interbred with other hominin groups. By the late 1990s, published genetic work had largely concluded that archaic groups made no lasting genetic contribution to modern humans; less than a decade later, this conclusion was reversed following the successful DNA sequencing of an ancient Neanderthal. This reversal of consensus is noteworthy, but the reasoning behind it is not widely understood across all academic communities. There remains a communication gap between population geneticists and paleoanthropologists. In this review, we endeavor to bridge this gap by outlining how technological advancements, new statistical methods, and notable controversies ultimately led to the current consensus.

GGVD: A goat genome variation database for tracking the dynamic evolutionary process of selective signatures and ancient introgressions

Understanding the evolutionary history and adaptive process depends on the knowledge that we can acquire from both ancient and modern genomic data. With the availability of a deluge of whole-genome sequencing data from ancient and modern goat samples, a user-friendly database making efficient reuse of these important resources is needed. Here, we use the genomes of 208 modern domestic goats, 24 bezoars, 46 wild ibexes, and 82 ancient goats to present a comprehensive goat genome variation database (GGVD). GGVD hosts a total of ∼41.44 million SNPs, ∼5.14 million indels, 6,193 selected loci, and 112 introgression regions. Users can freely visualize the frequency of genomic variations in geographical maps, selective sweeps in interactive tables, Manhattan plots, or line charts, as well as the heatmap patterns of the SNP genotype. Ancient data can be shown in haplotypes to track the state of genetic variants of selection and introgression events in the early, middle, and late stages. For facilitating access to sequence features, the UCSC Genome Browser, BLAT, BLAST, LiftOver, and pcadapt are also integrated into GGVD. GGVD will be a convenient tool for population genetic studies and molecular marker designing in goat breeding programs, and it is publicly available at http://animal.nwsuaf.edu.cn/GoatVar.

Scanning the human genome for "signatures" of positive selection: Transformative opportunities and ethical obligations

The relationship history of evolutionary anthropology and genetics is complex. At best, genetics is a beautifully integrative part of the discipline. Yet this integration has also been fraught, with punctuated, disruptive challenges to dogma, periodic reluctance by some members of the field to embrace results from analyses of genetic data, and occasional over-assertions of genetic definitiveness by geneticists. At worst, evolutionary genetics has been a tool for reinforcing racism and colonialism. While a number of genetics/genomics papers have disproportionately impacted evolutionary anthropology, here we highlight the 2002 presentation of an elegantly powerful approach for identifying "signatures" of past positive selection from haplotype-based patterns of genetic variation. Together with technological advances in genotyping methods, this article transformed our field by facilitating genome-wide "scans" for signatures of past positive selection in human populations. This approach helped researchers test longstanding evolutionary anthropology hypotheses while simultaneously providing opportunities to develop entirely new ones. Genome-wide scans for signatures of positive selection have since been conducted in diverse worldwide populations, with striking findings of local adaptation and convergent evolution. Yet there are ethical considerations with respect to the ubiquity of these studies and the cross-application of the genome-wide scan approach to existing datasets, which we also discuss.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

An efficient pipeline for ancient DNA mapping and recovery of endogenous ancient DNA from whole-genome sequencing data

Ancient DNA research has developed rapidly over the past few decades due to improvements in PCR and next-generation sequencing (NGS) technologies, but challenges still exist. One major challenge in relation to ancient DNA research is to recover genuine endogenous ancient DNA sequences from raw sequencing data. This is often difficult due to degradation of ancient DNA and high levels of contamination, especially homologous contamination that has extremely similar genetic background with that of the real ancient DNA. In this study, we collected whole-genome sequencing (WGS) data from 6 ancient samples to compare different mapping algorithms. To further explore more effective methods to separate endogenous DNA from homologous contaminations, we attempted to recover reads based on ancient DNA specific characteristics of deamination, depurination, and DNA fragmentation with different parameters. We propose a quick and improved pipeline for separating endogenous ancient DNA while simultaneously decreasing homologous contaminations to very low proportions. Our goal in this research was to develop useful recommendations for ancient DNA mapping and for separation of endogenous DNA to facilitate future studies of ancient DNA.

Adaptation and co-adaptation of skin pigmentation and vitamin D genes in native Americans

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest.

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records. This article is part of the theme issue 'Cross-disciplinary approaches to prehistoric demography'.

The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics

The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.

Evaluating genotype imputation pipeline for ultra-low coverage ancient genomes

Although ancient DNA data have become increasingly more important in studies about past populations, it is often not feasible or practical to obtain high coverage genomes from poorly preserved samples. While methods of accurate genotype imputation from > 1 × coverage data have recently become a routine, a large proportion of ancient samples remain unusable for downstream analyses due to their low coverage. Here, we evaluate a two-step pipeline for the imputation of common variants in ancient genomes at 0.05-1 × coverage. We use the genotype likelihood input mode in Beagle and filter for confident genotypes as the input to impute missing genotypes. This procedure, when tested on ancient genomes, outperforms a single-step imputation from genotype likelihoods, suggesting that current genotype callers do not fully account for errors in ancient sequences and additional quality controls can be beneficial. We compared the effect of various genotype likelihood calling methods, post-calling, pre-imputation and post-imputation filters, different reference panels, as well as different imputation tools. In a Neolithic Hungarian genome, we obtain ~ 90% imputation accuracy for heterozygous common variants at coverage 0.05 × and > 97% accuracy at coverage 0.5 ×. We show that imputation can mitigate, though not eliminate reference bias in ultra-low coverage ancient genomes.

A systematic investigation of human DNA preservation in medieval skeletons

Ancient DNA (aDNA) analyses necessitate the destructive sampling of archaeological material. Currently, the cochlea, part of the osseous inner ear located inside the petrous pyramid, is the most sought after skeletal element for molecular analyses of ancient humans as it has been shown to yield high amounts of endogenous DNA. However, destructive sampling of the petrous pyramid may not always be possible, particularly in cases where preservation of skeletal morphology is of top priority. To investigate alternatives, we present a survey of human aDNA preservation for each of ten skeletal elements in a skeletal collection from Medieval Germany. Through comparison of human DNA content and quality we confirm best performance of the petrous pyramid and identify seven additional sampling locations across four skeletal elements that yield adequate aDNA for most applications in human palaeogenetics. Our study provides a better perspective on DNA preservation across the human skeleton and takes a further step toward the more responsible use of ancient materials in human aDNA studies.

The Impact of Ancient Genome Studies in Archaeology

The study of ancient genomes has burgeoned at an incredible rate in the last decade. The result is a shift in archaeological narratives, bringing with it a fierce debate on the place of genetics in anthropological research. Archaeogenomics has challenged and scrutinized fundamental themes of anthropological research, including human origins, movement of ancient and modern populations, the role of social organization in shaping material culture, and the relationship between culture, language, and ancestry. Moreover, the discussion has inevitably invoked new debates on indigenous rights, ownership of ancient materials, inclusion in the scientific process, and even the meaning of what it is to be a human. We argue that the broad and seemingly daunting ethical, methodological, and theoretical challenges posed by archaeogenomics, in fact, represent the very cutting edge of social science research. Here, we provide a general review of the field by introducing the contemporary discussion points and summarizing methodological and ethical concerns, while highlighting the exciting possibilities of ancient genome studies in archaeology from an anthropological perspective.

Skin colour and vitamin D: An update

Homo sapiens evolved in East Africa and had dark skin, hair, and eyes, in order to protect against deleterious consequences of intensive UV radiation at equatorial latitudes. Intensive skin pigmentation was thought to bear the risk of inefficient vitamin D3 synthesis in the skin. This initiated the hypothesis that within the past 75 000 years, in which humans migrated to higher latitudes in Asia and Europe, the need for vitamin D3 synthesis served as an evolutionary driver for skin lightening. In this review, we summarize the recent archeogenomic reconstruction of population admixture in Europe and demonstrate that skin lightening happened as late as 5000 years ago through immigration of lighter pigmented populations from western Anatolia and the Russian steppe but not primarily via evolutionary pressure for vitamin D3 synthesis. We show that variations in genes encoding for proteins being responsible for the transport, metabolism and signalling of vitamin D provide alternative mechanisms of adaptation to a life in northern latitudes without suffering from consequences of vitamin D deficiency. This includes hypotheses explaining differences in the vitamin D status and response index of European populations.

The Pastoral Origin of Semiotically Functional Tonal Organization of Music

This paper presents a new line of inquiry into when and how music as a semiotic system was born. Eleven principal expressive aspects of music each contains specific structural patterns whose configuration signifies a certain affective state. This distinguishes the tonal organization of music from the phonetic and prosodic organization of natural languages and animal communication. The question of music’s origin can therefore be answered by establishing the point in human history at which all eleven expressive aspects might have been abstracted from the instinct-driven primate calls and used to express human psycho-emotional states. Etic analysis of acoustic parameters is the prime means of cross-examination of the typical patterns of expression of the basic emotions in human music versus animal vocal communication. A new method of such analysis is proposed here. Formation of such expressive aspects as meter, tempo, melodic intervals, and articulation can be explained by the influence of bipedal locomotion, breathing cycle, and heartbeat, long before Homo sapiens. However, two aspects, rhythm and melodic contour, most crucial for music as we know it, lack proxies in the Paleolithic lifestyle. The available ethnographic and developmental data leads one to believe that rhythmic and directional patterns of melody became involved in conveying emotion-related information in the process of frequent switching from one call-type to another within the limited repertory of calls. Such calls are usually adopted for the ongoing caretaking of human youngsters and domestic animals. The efficacy of rhythm and pitch contour in affective communication must have been spontaneously discovered in new important cultural activities. The most likely scenario for music to have become fully semiotically functional and to have spread wide enough to avoid extinctions is the formation of cross-specific communication between humans and domesticated animals during the Neolithic demographic explosion and the subsequent cultural revolution. Changes in distance during such communication must have promoted the integration between different expressive aspects and generated the basic musical grammar. The model of such communication can be found in the surviving tradition of Scandinavian pastoral music - kulning. This article discusses the most likely ways in which such music evolved.