English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

hapCon: Estimating contamination of ancient genomes by copying from reference haplotypes

MPS-Authors
/persons/resource/persons275566

Huang,  Yilei
Department of Archaeogenetics, Max Planck Institute for Evolutionary Anthropology, Max Planck Society;
The Leipzig School of Human Origins (IMPRS), Max Planck Institute for Evolutionary Anthropology, Max Planck Society;

/persons/resource/persons264741

Ringbauer,  Harald       
Department of Archaeogenetics, Max Planck Institute for Evolutionary Anthropology, Max Planck Society;

External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)

Huang_hapCon_Bioinf_2022.pdf
(Publisher version), 5MB

Supplementary Material (public)

Huang_hapCon_Bioinf_Suppl_2022.pdf
(Supplementary material), 11MB

Citation

Huang, Y., & Ringbauer, H. (2022). hapCon: Estimating contamination of ancient genomes by copying from reference haplotypes. Bioinformatics, 38(15): btac390, pp. 3768-3777. doi:10.1093/bioinformatics/btac390.


Cite as: https://hdl.handle.net/21.11116/0000-000A-A869-3
Abstract
Human ancient DNA (aDNA) studies have surged in recent years, revolutionizing the study of the human past. Typically, aDNA is preserved poorly, making such data prone to contamination from other human DNA. Therefore, it is important to rule out substantial contamination before proceeding to downstream analysis. As most aDNA samples can only be sequenced to low coverages (<1× average depth), computational methods that can robustly estimate contamination in the low coverage regime are needed. However, the ultra low-coverage regime (0.1× and below) remains a challenging task for existing approaches.We present a new method to estimate contamination in aDNA for male modern humans. It utilizes a Li&Stephens haplotype copying model for haploid X chromosomes, with mismatches modeled as errors or contamination. We assessed this new approach, hapCon, on simulated and down-sampled empirical aDNA data. Our experiments demonstrate that hapCon outperforms a commonly used tool for estimating male X contamination (ANGSD), with substantially lower variance and narrower confidence intervals, especially in the low coverage regime. We found that hapCon provides useful contamination estimates for coverages as low as 0.1× for SNP capture data (1240k) and 0.02× for whole genome sequencing data, substantially extending the coverage limit of previous male X chromosome-based contamination estimation methods. Our experiments demonstrate that hapCon has little bias for contamination up to 25–30\% as long as the contaminating source is specified within continental genetic variation, and that its application range extends to human aDNA as old as ∼45 000 and various global ancestries.We make hapCon available as part of a python package (hapROH), which is available at the Python Package Index (https://pypi.org/project/hapROH) and can be installed via pip. The documentation provides example use cases as blueprints for custom applications (https://haproh.readthedocs.io/en/latest/hapCon.html). The program can analyze either BAM files or pileup files produced with samtools. An implementation of our software (hapCon) using Python and C is deposited at https://github.com/hyl317/hapROH.Supplementary data are available at Bioinformatics online.