Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
RAD-seq SNPs: pooled population library prep? a.cardilini Sample Prep / Library Generation 2 06-12-2012 05:43 PM
Postdoc position in population genomics, CSIRO, Canberra habm Academic/Non-Profit Jobs 0 05-03-2012 10:43 PM
Bioinformatics position available: Population genetics and genomics pbonnen Academic/Non-Profit Jobs 0 07-20-2010 02:17 PM

Thread Tools
Old 11-25-2013, 03:47 PM   #1
Location: Austin, TX

Join Date: Sep 2011
Posts: 20
Post Permutation Strategies for RAD population genomics

I have a set of RAD data sampled from many individuals from several populations and I would like to calculate various population genetics statistics continuously along the genome with the ultimate goal of identifying regions of the genome that may be under selection.

I currently have a vcf file containing all genotyped bases including invariant sites, and I have developed a few scripts to calculate the stats I'm interested in from that file. Calculating the statistics the first time is easy, but the literature suggests that this is not enough. GWAS literature typically permutes the case/control (in my situation, "population") labels, generates a bootstrap sample via sampling with replacement, and recalculates the statistic, then repeats this process many times to produce a smooth distribution for this particular statistic under neutral expectations.

This is pretty much what Stacks already does very nicely. However, for several reasons, Stacks won't work with my Type II B RAD data.

And I cannot use the otherwise excellent vcftools package to calculate some of my statistics such as pi or Tajima's D, because it assumes that any base that is omitted from the vcf (as in, not genotyped) is invariant, which is most definitely not true for my data. Additionally, vcftools is more of a first-pass tool, and while I could probably work it into a pipeline, it doesn't do permutation/bootstrapping/etc by itself.

I am hoping that someone may have encountered similar problems using GWAS approaches for population genetics/population genomics and may have some advice for me.


1) Are there any existing tools besides Stacks and BayeScan (model-based approach for Fst) that do what I'm trying to do already?

2) Is this permute-bootstrap-recalc per locus strategy the best approach? Basically, my strategy is to develop a null distribution per each locus; Hohenlohe et al 2010, at least for pi, heterozygosity, Fst, and rho, generated a genome-wide null distribution via calculating the stats once per locus, then sampling those values from across the genome, weighting for sample size, to develop the genome-wide null distribution.

As a third option, someone suggested that simply testing for autocorrelation/LD among my pop gen statistics may be enough (i.e., "do I expect to see these three SNPs with a value greater than X as neighbors, or could this happen by chance?"). I'm not convinced; also, this approach necessitates defining whatever "X" cutoff you're using in some manner.

Which approach is more robust? I think that the per-locus approach may be more conservative, and may also produce fewer false positives... but I am also terrible at statistics. Thoughts? Alternatives?

3) Regardless of strategy, how many times should the resampling process be repeated? Until the null distribution "looks smooth enough"? How can this be tested?

4) What are the benefits of using a sliding window approach? I understand the kernel smoothing idea where one weights each point's value by both the value and distance of the statistics at neighboring SNPs, and I like it -- but is it appropriate to use if my null distribution is locus-specific and not genome-wide?

Thank you for any hints, comments or advice you can give me!
rcapper is offline   Reply With Quote
Old 11-25-2013, 10:34 PM   #2
Junior Member
Location: Eugene, OR

Join Date: Nov 2013
Posts: 1
Default 2bRAD, Stacks, and Fst

Hi Roxana,

What about the 2bRAD data makes it not useable in Stacks? I know that the demultiplexing/cleaning program in Stacks (process_radtags) does not yet support 2bRAD data, but I assume you have cleaned and demultiplexed your data already. Are there other factors that make Stacks unable to process it?

I don't have any 2bRAD data myself but would like to see it supported by our software.


catchen is offline   Reply With Quote

genomics, population genetics, rad, stacks, vcftools

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:51 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO