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.

Benchmarking Imputed Low Coverage Genomes in a Human Population Genetics Context

Evaluation of genotype imputation using Glimpse tools on low coverage ancient DNA

Comprehensive Assessment of Genotype Imputation Performance.

Evaluation of ancient DNA imputation: a simulation study

Imputation accuracy across global human populations

Benchmarking for genotyping and imputation using degraded DNA for forensic applications across diverse populations

Rapid and accurate genotype imputation from low coverage short read, long read, and cell free DNA sequence

Genotype-Imputation Accuracy Across Worldwide Human Populations

A Novel Efficient Algorithm for Common Variants Genotyping from Low-Coverage Sequencing Data

An empirical evaluation of genotype imputation of ancient DNA

A New Genotype Imputation Method with Tolerance to High Missing Rate and Rare Variants

Performance of Genotype Imputation for Low Frequency and Rare Variants from the 1000 Genomes

A Comprehensive SNP and Indel Imputability Database

FastImpute: A Baseline for Open-source, Reference-Free Genotype Imputation Methods -- A Case Study in PRS313

On Combining Reference Data to Improve Imputation Accuracy

Estimating heterozygosity from a low-coverage genome sequence, leveraging data from other individuals sequenced at the same sites

Improved imputation of low-frequency and rare variants using the UK10K haplotype reference panel

Improving Power of Association Tests Using Multiple Sets of Imputed Genotypes from Distributed Reference Panels

MaCH-admix: Genotype Imputation for Admixed Populations.

Integer programming framework for pangenome-based genome inference

Performance of Genotype Imputation for Rare Variants Identified in Exons and Flanking Regions of Genes