Estimating inbreeding coefficients from NGS data: impact on genotype calling and allele frequency estimation

Research output: Contribution to journalJournal articleResearchpeer-review

68 Citations (Scopus)

Abstract

Most methods for Next-Generation Sequencing (NGS) data analyses incorporate information regarding allele frequencies using the assumption of Hardy-Weinberg Equilibrium (HWE) as a prior. However, many organisms including domesticated, partially selfing or with asexual life cycles show strong deviations from HWE. For such species, and specially for low coverage data, it is necessary to obtain estimates of inbreeding coefficients (F) for each individual beforecalling genotypes. Here, we present two methods for estimating inbreeding coefficients from NGS data based on an Expectation-Maximization (EM) algorithm. We assess the impact of taking inbreeding into account when calling genotypes or estimating the Site Frequency Spectrum (SFS), and demonstrate a marked increase in accuracy on low coverage highly inbred samples. We demonstrate the applicability and efficacy of these methods in both simulated and real datasets.
Original languageEnglish
JournalGenome Research
Volume23
Pages (from-to)1852-1861
Number of pages10
ISSN1088-9051
DOIs
Publication statusPublished - 2013

Cite this