Dataset Information

Assessing the effect of sequencing depth and sample size in population genetics inferences.

ABSTRACT: Next-Generation Sequencing (NGS) technologies have dramatically revolutionised research in many fields of genetics. The ability to sequence many individuals from one or multiple populations at a genomic scale has greatly enhanced population genetics studies and made it a data-driven discipline. Recently, researchers have proposed statistical modelling to address genotyping uncertainty associated with NGS data. However, an ongoing debate is whether it is more beneficial to increase the number of sequenced individuals or the per-sample sequencing depth for estimating genetic variation. Through extensive simulations, I assessed the accuracy of estimating nucleotide diversity, detecting polymorphic sites, and predicting population structure under different experimental scenarios. Results show that the greatest accuracy for estimating population genetics parameters is achieved by employing a large sample size, despite single individuals being sequenced at low depth. Under some circumstances, the minimum sequencing depth for obtaining accurate estimates of allele frequencies and to identify polymorphic sites is [Formula: see text], where both alleles are more likely to have been sequenced. On the other hand, inferences of population structure are more accurate at very large sample sizes, even with extremely low sequencing depth. This all points to the conclusion that under various experimental scenarios, in cost-limited population genetics studies, large sample sizes at low sequencing depth are desirable to achieve high accuracy. These findings will help researchers design their experimental set-ups and guide further investigation on the effect of protocol design for genetic research.

SUBMITTER: Fumagalli M

PROVIDER: S-EPMC3832539 | biostudies-literature | 2013

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Assessing the effect of sequencing depth and sample size in population genetics inferences.

Fumagalli Matteo M

PloS one 20131118 11

Next-Generation Sequencing (NGS) technologies have dramatically revolutionised research in many fields of genetics. The ability to sequence many individuals from one or multiple populations at a genomic scale has greatly enhanced population genetics studies and made it a data-driven discipline. Recently, researchers have proposed statistical modelling to address genotyping uncertainty associated with NGS data. However, an ongoing debate is whether it is more beneficial to increase the number of ...[more]

PMID: 24260275

Similar Datasets

Project description:In landscape genetics, it is largely unknown how choices regarding sampling density and study area size impact inferences upon which habitat features impede vs. facilitate gene flow. While it is recommended that sampling locations be spaced no further apart than the average individual's dispersal distance, for low-mobility species, this could lead to a challenging number of sampling locations, or an unrepresentative study area. We assessed the effects of sampling density and study area size on landscape genetic inferences for a dispersal-limited amphibian, Plethodon mississippi, via analysis of nested datasets. Microsatellite-based genetic distances among individuals were divided into three datasets representing sparse sampling across a large study area, dense sampling across a small study area, or sparse sampling across the same small study area. These datasets were a proxy for gene flow (i.e., the response variable) in maximum-likelihood population effects models that assessed the nature and strength of their relationship with each of five land-use classes (i.e., potential predictor variables). Comparisons of outcomes were based on the rank order of effect, sign of effect (i.e., gene flow resistance vs. facilitation), spatial scale of effect, and functional relationship with gene flow. The best-fit model for each dataset had the same sign of effect for hardwood forests, manmade structures, and pine forests, indicating the impacts of these land-use classes on dispersal and gene flow in P. mississippi are robust to sampling scheme. Contrasting sampling densities led to a different inferred functional relationship between agricultural areas and gene flow. Study area size appeared to influence the scale of effect of manmade structures and the sign of effect of pine forests. Our findings provided evidence for an influence of sampling density, study area size, and sampling effort upon inferences. Accordingly, we recommend iterative subsampling of empirical datasets and continued investigation into the sensitivities of landscape genetic analyses using simulations.

Project description:Eradication or local extinction of the human parasite Schistosoma mansoni is a goal for many control programs. Population genetic analyses are helping to evaluate and guide these efforts, yet what to sample, how to sample and how densely to sample is not well established. We determined the S. mansoni allele frequency profile of nearly all infected inhabitants in two small Brazilian communities and created sub-samples representing 5-50% of all detected human infections (infrapopulations). Samples were selected at random with replacement, and each size class was replicated 100 times. Mean pairwise differentiation for all infrapopulations (Di) and the variance effective population size (Ne) were calculated for each sample. Prior to community-wide treatment, the true mean Di was moderate (0.095-0.123) and Ne large (>30,000). Most samples of <50% of those infected produced estimates outside of 5% of the true value. For estimates within 10%, sample sizes of >15% of all infrapopulations were required. At the 3 year follow-up after treatment, the Di increased and Ne was reduced by >15 fold. At this time sampling of >30-45% was needed to achieve the same accuracy. Following a second treatment and 4 years from baseline, the Di further increased and Ne decreased with little change in the sampling effort required. Extensive sampling is required for accurate estimates of these important population parameters. Characteristics such as population census size, infection prevalence, the community's treatment history and the degree of infrapopulation differentiation should be taken into account. The intensity of infection was weakly correlated with the ability of a single infrapopulation to represent the component population (Dic), indicating a tendency toward random acquisition of parasite genotypes. This also suggests that targeted sampling from those most heavily infected will better represent the genetic diversity of the whole community than a random sample of infrapopulations.

Dataset Information

Assessing the effect of sequencing depth and sample size in population genetics inferences.

Publications

Assessing the effect of sequencing depth and sample size in population genetics inferences.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets