Estimating inbreeding coefficients from NGS data: impact on genotype calling and allele frequency estimation
- ↵* Corresponding author; email: fgarrettvieira{at}gmail.com
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.
- Received March 7, 2013.
- Accepted August 6, 2013.
- © 2013, Published by Cold Spring Harbor Laboratory Press
This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date (see http://genome.cshlp.org/site/misc/terms.xhtml). After six months, it is available under a Creative Commons License (Attribution-NonCommercial 3.0 Unported), as described at http://creativecommons.org/licenses/by-nc/3.0/.
