Published online before print
September 5, 2006 Genome Research, DOI: 10.1101/gr.5431206
Methods
Inference of population genetic parameters in metagenomics: A clean look at messy data
Philip L.F. Johnson1,3 and
Montgomery Slatkin2
1 Biophysics Graduate Group, University of California, Berkeley, California 94720, USA;
2 Department of Integrative Biology, University of California, Berkeley, California 94720, USA
Metagenomic projects generate short, overlapping fragments of DNA sequence, each deriving from a different individual. We report a new method for inferring the scaled mutation rate, = 2Neu, and the scaled exponential growth rate, R = Ner, from the site-frequency spectrum of these data while accounting for sequencing error via Phred quality scores. After obtaining maximum likelihood parameter estimates for and R, we calculate empirical Bayes quality scores reflecting the posterior probability that each apparently polymorphic site is truly polymorphic; these scores can then be used for other applications such as SNP discovery. For realistic parameter ranges, analytic and simulation results show our estimates to be essentially unbiased with tight confidence intervals. In contrast, choosing an arbitrary quality score cutoff (e.g., trimming reads) and ignoring further quality information during inference yields biased estimates with greater variance. We illustrate the use of our technique on a new project analyzing activated sludge from a lab-scale bioreactor seeded by a wastewater treatment plant.
3 Corresponding author.
E-mail plfjohnson{at}berkeley.edu; fax (510) 643-6264.
Article published online before print. Article and publication date are online at http://www.genome.org/cgi/doi/10.1101/gr.5431206

CiteULike Connotea Del.icio.us Digg Reddit Technorati What's this?
This article has been cited by other articles:

|
 |

|
 |
 
P. L. F. Johnson and M. Slatkin
Accounting for Bias from Sequencing Error in Population Genetic Estimates
Mol. Biol. Evol.,
January 1, 2008;
25(1):
199 - 206.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Chen, R. E. Green, S. Paabo, and M. Slatkin
The Joint Allele-Frequency Spectrum in Closely Related Species
Genetics,
September 1, 2007;
177(1):
387 - 398.
[Abstract]
[Full Text]
[PDF]
|
 |
|
|
|