Estimating error models for whole genome sequencing using mixtures of Dirichlet-multinomial distributions

Steven H. Wu, Rachel S. Schwartz, David J. Winter, Donald F. Conrad, Reed Cartwright

Research output: Contribution to journalArticle

4 Scopus citations

Abstract

Motivation: Accurate identification of genotypes is an essential part of the analysis of genomic data, including in identification of sequence polymorphisms, linking mutations with disease and determining mutation rates. Biological and technical processes that adversely affect genotyping include copy-number-variation, paralogous sequences, library preparation, sequencing error and reference-mapping biases, among others. Results: We modeled the read depth for all data as a mixture of Dirichlet-multinomial distributions, resulting in significant improvements over previously used models. In most cases the best model was comprised of two distributions. The major-component distribution is similar to a binomial distribution with low error and low reference bias. The minor-component distribution is overdispersed with higher error and reference bias. We also found that sites fitting the minor component are enriched for copy number variants and low complexity regions, which can produce erroneous genotype calls. By removing sites that do not fit the major component, we can improve the accuracy of genotype calls.

Original languageEnglish (US)
Pages (from-to)2322-2329
Number of pages8
JournalBioinformatics
Volume33
Issue number15
DOIs
StatePublished - Aug 1 2017

    Fingerprint

ASJC Scopus subject areas

  • Statistics and Probability
  • Biochemistry
  • Molecular Biology
  • Computer Science Applications
  • Computational Theory and Mathematics
  • Computational Mathematics

Cite this