Dataset Information

Estimating effective population size using RADseq: Effects of SNP selection and sample size.

ABSTRACT: Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.

SUBMITTER: Marandel F

PROVIDER: S-EPMC7042749 | biostudies-literature | 2020 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Estimating effective population size using RADseq: Effects of SNP selection and sample size.

Marandel Florianne F Charrier Grégory G Lamy Jean-Baptiste JB Le Cam Sabrina S Lorance Pascal P Trenkel Verena M VM

Ecology and evolution 20200211 4

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To eval ...[more]

PMID: 32128126

Similar Datasets

Project description:Sample size is a critical aspect of study design in population genomics research, yet few empirical studies have examined the impacts of small sample sizes. We used datasets from eight diverging bird lineages to make pairwise comparisons at different levels of taxonomic divergence (populations, subspecies, and species). Our data are from loci linked to ultraconserved elements and our analyses used one single nucleotide polymorphism per locus. All individuals were genotyped at all loci, effectively doubling sample size for coalescent analyses. We estimated population demographic parameters (effective population size, migration rate, and time since divergence) in a coalescent framework using Diffusion Approximation for Demographic Inference, an allele frequency spectrum method. Using divergence-with-gene-flow models optimized with full datasets, we subsampled at sequentially smaller sample sizes from full datasets of 6-8 diploid individuals per population (with both alleles called) down to 1:1, and then we compared estimates and their changes in accuracy. Accuracy was strongly affected by sample size, with considerable differences among estimated parameters and among lineages. Effective population size parameters (?) tended to be underestimated at low sample sizes (fewer than three diploid individuals per population, or 6:6 haplotypes in coalescent terms). Migration (m) was fairly consistently estimated until <2 individuals per population, and no consistent trend of over-or underestimation was found in either time since divergence (T) or theta (? = 4N ref?). Lineages that were taxonomically recognized above the population level (subspecies and species pairs; that is, deeper divergences) tended to have lower variation in scaled root mean square error of parameter estimation at smaller sample sizes than population-level divergences, and many parameters were estimated accurately down to three diploid individuals per population. Shallower divergence levels (i.e., populations) often required at least five individuals per population for reliable demographic inferences using this approach. Although divergence levels might be unknown at the outset of study design, our results provide a framework for planning appropriate sampling and for interpreting results if smaller sample sizes must be used.

Project description:The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (Ne) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r2=0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of Ne was noted in the very recent generations with a peak value of 19 for Ne at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for FROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (FPed=0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively. The objectives of the present study were to genotype 28 Lundehund using the canine Illumina high density beadchip (Illumina, San Diego, CA, USA) to estimate the effective population size (Ne) from data on linkage disequilibria (LD). We identified runs of homozygosity (ROH) as regions with a local loss of genetic variation. Inbreeding coefficients were calculated on the basis of pedigree data (FPed) and genotype information (FROH, FIS) to compare the results of these different methods for this population. We applied the long-range haplotype test to search for potential selective sweeps.

Dataset Information

Estimating effective population size using RADseq: Effects of SNP selection and sample size.

Publications

Estimating effective population size using RADseq: Effects of SNP selection and sample size.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure