Abstract:Ancient genomic data is becoming increasingly available thanks to recent advances in high-throughput sequencing technologies. Yet, post-mortem degradation of endogenous ancient DNA often results in low depth of coverage and subsequently high levels of genotype missingness and uncertainty. Genotype imputation is a potential strategy for increasing the information available in ancient DNA samples and thus improving the power of downstream population genetic analyses. However, the performance of genotype imputation on ancient genomes under different conditions has not yet been fully explored, with all previous work primarily using an empirical approach of downsampling high coverage paleogenomes. While these studies have provided invaluable insights into best practices for imputation, they rely on a fairly limited number of existing high coverage samples with significant temporal and geographical biases. As an alternative, we used a coalescent simulation approach to generate genomes with characteristics of ancient DNA in order to more systematically evaluate the performance of two popular imputation software, BEAGLE and GLIMPSE, under variable divergence times between the target sample and reference haplotypes, as well as different depths of coverage and reference sample size. Our results suggest that for genomes with coverage <=0.1x imputation performance is poor regardless of the strategy employed. Beyond 0.1x coverage imputation is generally improved as the size of the reference panel increases, and imputation accuracy decreases with increasing divergence between target and reference populations. It may thus be preferable to compile a smaller set of less diverged reference samples than a larger more highly diverged dataset. In addition, the imputation accuracy may plateau beyond some level of divergence between the reference and target populations. While accuracy at common variants is similar regardless of divergence time, rarer variants are better imputed on less diverged target samples. Furthermore, both imputation software, but particularly GLIMPSE, overestimate high genotype probability calls, especially at low coverages. Our results provide insight into optimal strategies for ancient genotype imputation under a wide set of scenarios, complementing previous empirical studies based on imputing downsampled high-coverage ancient genomes.

Deep estimation of the intensity and timing of natural selection from ancient genomes

Predicting functional consequences of recent natural selection in Britain

Modeling the spatiotemporal spread of beneficial alleles using ancient genomes

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

Estimating temporally variable selection intensity from ancient DNA Data

Inferring the model and onset of natural selection under varying population size from the site frequency spectrum and haplotype structure

Estimating temporally variable selection intensity from ancient DNA data with the flexibility of modelling linkage and epistasis

Inference of natural selection from ancient DNA

Bayesian Inference of Natural Selection from Allele Frequency Time Series

Inferring the timing and strength of natural selection and gene migration in the evolution of chicken from ancient DNA data

Fast and Accurate Estimation of Selection Coefficients and Allele Histories from Ancient and Modern DNA

A novel expectation-maximization approach to infer general diploid selection from time-series genetic data

A novel spectral method for inferring general diploid selection from time series genetic data

Evaluation of ancient DNA imputation: a simulation study

Accurate inference of population history in the presence of background selection

A Chronological Atlas of Natural Selection in the Human Genome during the Past Half-million Years

Natural selection acting on complex traits hampers the predictive accuracy of polygenic scores in ancient samples

Testing for Ancient Selection Using Cross-population Allele Frequency Differentiation

Natural Selection Exerted by Historical Coronavirus Epidemic(s): Comparative Genetic Analysis in China Kadoorie Biobank and UK Biobank

The Selection Landscape and Genetic Legacy of Ancient Eurasians

Inferring Allele Frequency Trajectories from Ancient DNA Indicates That Selection on a Chicken Gene Coincided with Changes in Medieval Husbandry Practices