Matched full-UDG and non-UDG ancient DNA libraries reveal trade-offs in post-mortem damage correction for imputation and kinship inference
- 1Sapiens Life Sciences, Evolution and Medicine Research Center, Kanazawa University
Abstract
Ancient DNA studies increasingly combine full uracil-DNA glycosylase-treated, full-UDG, and non-UDG libraries, but which computational damage correction to apply before imputation and kinship analysis remains unsettled. We compared terminal trimming, base-quality rescaling and known-SNP masking in matched full-UDG and non-UDG libraries from the same extracts of two medieval Mongolian individuals. All methods reduced or reversed the difference in mean alternative-allele fraction between damage-prone and transversion SNPs but retained different proportions of covered sites. In non-UDG libraries, masking and rescaling within five bases of each fragment end produced similar cross-library non-reference discordance, NRD, while retaining 92% and 99.3% of covered sites, respectively. In these 3'-biased libraries, ten-base 3'-only trimming retained more covered sites at lower observed discordance than five-base-per-end symmetric trimming. ancIBD inferred widespread sharing of one identity-by-descent chromosome copy, IBD1, across correction methods. Uncorrected non-UDG data shifted TKGWV2 toward second-degree estimates, whereas corrected mixed-library comparisons supported first-degree relatedness. READv2 and the low rate of opposite-homozygote sites, IBS0, supported a parent-offspring relationship. Mitochondrial data and genetic sex favoured ORT16 as mother and ORT15 as son. No method performed best across all measures, so the choice depends on whether residual damage, site retention or the downstream analysis matters most.
Related
- A Unified Framework for Batch Correction and Imputation in Modern Mass Spectrometry Proteomics
- Archaic ancestry inference in imputed ancient human genomes
- QuImputer: an LD-informed QUBO/Ising formulation of genotype imputation with quantum-hardware execution
- Singular Value Decomposition for the Automated Analysis of NMR Fraction Libraries
- High-throughput variant libraries and machine learning enable robust target-site selection for toehold-switch design in Detectrons
BibTeX
@misc{ravdandorj2026matched,
title = {{Matched full-UDG and non-UDG ancient DNA libraries reveal trade-offs in post-mortem damage correction for imputation and kinship inference}},
author = {O. Ravdandorj and C. Sampildondov and K. Janchiv and T. Gakuhari},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.07.749785},
url = {https://www.biorxiv.org/content/10.64898/2026.09.07.749785v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Matched full-UDG and non-UDG ancient DNA libraries reveal trade-offs in post-mortem damage correction for imputation and kinship inference
AU - O. Ravdandorj
AU - C. Sampildondov
AU - K. Janchiv
AU - T. Gakuhari
PY - 2026
DA - 2026/09/07
PB - bioRxiv
DO - 10.64898/2026.09.07.749785
UR - https://www.biorxiv.org/content/10.64898/2026.09.07.749785v1
AB - Ancient DNA studies increasingly combine full uracil-DNA glycosylase-treated, full-UDG, and non-UDG libraries, but which computational damage correction to apply before imputation and kinship analysis remains unsettled. We compared terminal trimming, base-quality rescaling and known-SNP masking in matched full-UDG and non-UDG libraries from the same extracts of two medieval Mongolian individuals. All methods reduced or reversed the difference in mean alternative-allele fraction between damage-prone and transversion SNPs but retained different proportions of covered sites. In non-UDG libraries, masking and rescaling within five bases of each fragment end produced similar cross-library non-reference discordance, NRD, while retaining 92% and 99.3% of covered sites, respectively. In these 3'-biased libraries, ten-base 3'-only trimming retained more covered sites at lower observed discordance than five-base-per-end symmetric trimming. ancIBD inferred widespread sharing of one identity-by-descent chromosome copy, IBD1, across correction methods. Uncorrected non-UDG data shifted TKGWV2 toward second-degree estimates, whereas corrected mixed-library comparisons supported first-degree relatedness. READv2 and the low rate of opposite-homozygote sites, IBS0, supported a parent-offspring relationship. Mitochondrial data and genetic sex favoured ORT16 as mother and ORT15 as son. No method performed best across all measures, so the choice depends on whether residual damage, site retention or the downstream analysis matters most.
ER -