Cited by 21Forensic and Genetic ResearchGenomics and Phylogenetic StudiesGenetic diversity and population structure

A likelihood method for estimating present-day human contamination in ancient male samples using low-depth X-chromosome data

Open access
Moreno-Mayar, J Víctor, Korneliussen, Thorfinn Sand, Dalal, Jyoti, Renaud, Gabriel, Albrechtsen, Anders, Nielsen, Rasmus, Malaspinas, Anna-Sapfo·2020·Bioinformatics

Published version of a preprint

Abstract

Motivation: The presence of present-day human contaminating DNA fragments is one of the challenges defining ancient DNA (aDNA) research. This is especially relevant to the ancient human DNA field where it is difficult to distinguish endogenous molecules from human contaminants due to their genetic similarity. Recently, with the advent of high-throughput sequencing and new aDNA protocols, hundreds of ancient human genomes have become available. Contamination in those genomes has been measured with computational methods often developed specifically for these empirical studies. Consequently, some of these methods have not been implemented and tested for general use while few are aimed at low-depth nuclear data, a common feature in aDNA datasets. Results: We develop a new X-chromosome-based maximum likelihood method for estimating present-day human contamination in low-depth sequencing data from male individuals. We implement our method for general use, assess its performance under conditions typical of ancient human DNA research, and compare it to previous nuclear data-based methods through extensive simulations. For low-depth data, we show that existing methods can produce unusable estimates or substantially underestimate contamination. In contrast, our method provides accurate estimates for a depth of coverage as low as 0.5× on the X-chromosome when contamination is below 25%. Moreover, our method still yields meaningful estimates in very challenging situations, i.e. when the contaminant and the target come from closely related populations or with increased error rates. With a running time below 5 min, our method is applicable to large scale aDNA genomic studies. Availability and implementation: The method is implemented in C++ and R and is available in github.com/sapfo/contaminationX and popgen.dk/angsd.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1093/bioinformatics/btz660

Cites

35 of 35 references verified0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Moreno-Mayar, J. V., Korneliussen, T. S., Dalal, J., Renaud, G., Albrechtsen, A., Nielsen, R., & Malaspinas, A.-S.
Affiliations
  • SIB Swiss Institute of Bioinformatics (CH)
  • University of Copenhagen (DK)
Corresponding author
  • Moreno-Mayar, J Víctor
  • Malaspinas, Anna-Sapfo
Venue
Bioinformatics
Published
2019-08-26
Volume
36
Issue
3
Pages
828-841
ISSN
  • 1367-4803 (ISSN-L)
  • 1367-4811
OpenAlex
W2970150218
Cited by
21
References
35
Citations by year
2026: 3 · 2025: 5 · 2024: 2 · 2023: 1 · 2022: 2 · 2021: 3 · 2020: 5
FWCI
1.51
Type
journal article
Language
en
Open access
Locations
  • Bioinformatics (publishedVersion)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • PubMed Central (publishedVersion)
  • PubMed (publishedVersion)
  • IRIS (submittedVersion)
  • Research at the University of Copenhagen (University of Copenhagen) (submittedVersion, cc-by-nc-nd)
  • eScholarship (California Digital Library) (submittedVersion)
  • eScholarship (California Digital Library) (submittedVersion, cc-by-nc-nd)
  • SERVAL (Université de Lausanne) (submittedVersion)
  • SERVAL (Université de Lausanne) (submittedVersion, other-oa)
Funders
  • Danish National Research Foundation (DNRF94)
  • European Research Council (679330)
  • Carlsberg Foundation (CF16-0913)
  • NIH (NIH 1R01GM116044-01)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Forensic and Genetic Research
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic diversity and population structure
Keywords
ancient DNA, contamination estimation, X chromosome, low-pass sequencing, maximum likelihood estimation
MeSH terms
DNA, Ancient, Chromosomes, Humans, Male, Likelihood Functions, Sequence Analysis, DNA, High-Throughput Nucleotide Sequencing

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review