Population genomics based on low coverage sequencing: how low should we go?
- 22 November 2012
- journal article
- research article
- Published by Wiley in Molecular Ecology
- Vol. 22 (11), 3028-3035
- https://doi.org/10.1111/mec.12105
Abstract
Research in molecular ecology is now often based on large numbers of DNA sequence reads. Given a time and financial budget for DNA sequencing, the question arises as to how to allocate the finite number of sequence reads among three dimensions: (i) sequencing individual nucleotide positions repeatedly and achieving high confidence in the true genotype of individuals, (ii) sampling larger numbers of individuals from a population, and (iii) sampling a larger fraction of the genome. Leaving aside the question of what fraction of the genome to sample, we analyze the trade-off between repeatedly sequencing the same nucleotide position (coverage depth) and the number of individuals in the sample. We review simple Bayesian models for allele frequencies and utilize these in the analysis of how to obtain maximal information about population genetic parameters. The models indicate that sampling larger numbers of individuals, at the expense of coverage depth per nucleotide position, provides more information about population parameters. Dividing the sequencing effort maximally among individuals and obtaining approximately one read per locus and individual (1 × coverage) yields the most information about a population. Some analyses require genetic parameters for individuals, in which case Bayesian population models also support inference from lower coverage sequence data than are required for simple likelihood models. Low coverage sequencing is not only sufficient to support inference, but it is optimal to design studies to utilize low coverage because they will yield highly accurate and precise parameter estimates based on more individuals or sites in the genome.Keywords
This publication has 41 references indexed in Scilit:
- Reconstructing Native American population historyNature, 2012
- Extremely low-coverage sequencing and imputation increases power for genome-wide association studiesNature Genetics, 2012
- The genomic basis of adaptive evolution in threespine sticklebacksNature, 2012
- Genomic islands of divergence in hybridizing Heliconius butterflies identified by large-scale targeted sequencingPhilosophical Transactions Of The Royal Society B-Biological Sciences, 2012
- Genomics of isolation in hybridsPhilosophical Transactions Of The Royal Society B-Biological Sciences, 2012
- Multiplexed shotgun genotyping for rapid and efficient genetic mappingGenome Research, 2011
- The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing dataGenome Research, 2010
- Isolation with Migration Models for More Than Two PopulationsMolecular Biology and Evolution, 2009
- The Sequence Alignment/Map format and SAMtoolsBioinformatics, 2009
- Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencingNature Biotechnology, 2009