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.

A high coverage Mesolithic aurochs genome and effective leveraging of ancient cattle genomes using whole genome imputation

A high coverage Mesolithic aurochs genome and effective leveraging of ancient cattle genomes using whole genome imputation

The genomic natural history of the aurochs

The Genomic Legacy of Aurochs hybridization in ancient and modern Iberian Cattle

Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent

The Efficient Phasing and Imputation Pipeline of Low‐coverage Whole Genome Sequencing Data Using a High‐quality and Publicly Available Reference Panel in Cattle

The Construction of a Haplotype Reference Panel Using Extremely Low Coverage Whole Genome Sequences and Its Application in Genome-Wide Association Studies and Genomic Prediction in Duroc Pigs.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

Evaluation of ancient DNA imputation: a simulation study

An empirical evaluation of genotype imputation of ancient DNA

Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history

Ancestral alleles defined for 70 million cattle variants using a population-based likelihood ratio test

DNA methylation analysis to differentiate reference, breed, and parent-of-origin effects in the bovine pangenome era

Genomic evaluation of hybridization in historic and modern North American Bison (Bison bison)

Evolution and legacy of East Asian aurochs

Paging through history: parchment as a reservoir of ancient DNA for next generation sequencing

Genome-wide association analyses of multiple traits in Duroc pigs using low-coverage whole-genome sequencing strategy

Epigenomics and Genotype-Phenotype Association Analyses Reveal Conserved Genetic Architecture of Complex Traits in Cattle and Human

Assessment of the performance of different imputation methods for low-coverage sequencing in Holstein cattle

A Computational Framework for Improving Genetic Variants Identification from 5,061 Sheep Sequencing Data.

Recalibrating Equus Evolution Using the Genome Sequence of an Early Middle Pleistocene Horse.