SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data

Rasmus Nielsen; Thorfinn Sand Korneliussen; Anders Albrechtsen; Yingrui Li; Jun Wang

doi:10.1371/journal.pone.0037558

SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data

Rasmus Nielsen, Thorfinn Sand Korneliussen, Anders Albrechtsen, Yingrui Li, Jun Wang

186 Citations (Scopus)

3231 Downloads (Pure)

Abstract

We present a statistical framework for estimation and application of sample allele frequency spectra from New-Generation Sequencing (NGS) data. In this method, we first estimate the allele frequency spectrum using maximum likelihood. In contrast to previous methods, the likelihood function is calculated using a dynamic programming algorithm and numerically optimized using analytical derivatives. We then use a bayesian method for estimating the sample allele frequency in a single site, and show how the method can be used for genotype calling and SNP calling. We also show how the method can be extended to various other cases including cases with deviations from Hardy-Weinberg equilibrium. We evaluate the statistical properties of the methods using simulations and by application to a real data set.

Original language	English
Article number	e37558
Journal	PLoS ONE
Volume	7
Issue number	7
Number of pages	11
ISSN	1932-6203
DOIs	https://doi.org/10.1371/journal.pone.0037558
Publication status	Published - 24 Jul 2012

Access to Document

10.1371/journal.pone.0037558Licence: CC BY

Nielsen_2012_SNP_callingFinal published version, 1.25 MBLicence: CC BY

Cite this

@article{9115b976bf424b1f9a1de08605d7fcbd,

title = "SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data",

abstract = "We present a statistical framework for estimation and application of sample allele frequency spectra from New-Generation Sequencing (NGS) data. In this method, we first estimate the allele frequency spectrum using maximum likelihood. In contrast to previous methods, the likelihood function is calculated using a dynamic programming algorithm and numerically optimized using analytical derivatives. We then use a bayesian method for estimating the sample allele frequency in a single site, and show how the method can be used for genotype calling and SNP calling. We also show how the method can be extended to various other cases including cases with deviations from Hardy-Weinberg equilibrium. We evaluate the statistical properties of the methods using simulations and by application to a real data set.",

author = "Rasmus Nielsen and Korneliussen, {Thorfinn Sand} and Anders Albrechtsen and Yingrui Li and Jun Wang",

note = "e37558",

year = "2012",

month = jul,

day = "24",

doi = "10.1371/journal.pone.0037558",

language = "English",

volume = "7",

journal = "PLoS Computational Biology",

issn = "1932-6203",

publisher = "Public Library of Science",

number = "7",

}

TY - JOUR

T1 - SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data

AU - Nielsen, Rasmus

AU - Korneliussen, Thorfinn Sand

AU - Albrechtsen, Anders

AU - Li, Yingrui

AU - Wang, Jun

N1 - e37558

PY - 2012/7/24

Y1 - 2012/7/24

N2 - We present a statistical framework for estimation and application of sample allele frequency spectra from New-Generation Sequencing (NGS) data. In this method, we first estimate the allele frequency spectrum using maximum likelihood. In contrast to previous methods, the likelihood function is calculated using a dynamic programming algorithm and numerically optimized using analytical derivatives. We then use a bayesian method for estimating the sample allele frequency in a single site, and show how the method can be used for genotype calling and SNP calling. We also show how the method can be extended to various other cases including cases with deviations from Hardy-Weinberg equilibrium. We evaluate the statistical properties of the methods using simulations and by application to a real data set.

AB - We present a statistical framework for estimation and application of sample allele frequency spectra from New-Generation Sequencing (NGS) data. In this method, we first estimate the allele frequency spectrum using maximum likelihood. In contrast to previous methods, the likelihood function is calculated using a dynamic programming algorithm and numerically optimized using analytical derivatives. We then use a bayesian method for estimating the sample allele frequency in a single site, and show how the method can be used for genotype calling and SNP calling. We also show how the method can be extended to various other cases including cases with deviations from Hardy-Weinberg equilibrium. We evaluate the statistical properties of the methods using simulations and by application to a real data set.

U2 - 10.1371/journal.pone.0037558

DO - 10.1371/journal.pone.0037558

M3 - Journal article

C2 - 22911679

SN - 1932-6203

VL - 7

JO - PLoS Computational Biology

JF - PLoS Computational Biology

IS - 7

M1 - e37558

ER -

SNP calling, genotype calling, and sample allele frequency estimation from new-generation sequencing data

Abstract

Access to Document

Fingerprint

Cite this