SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
From VCF to Fst AdrianP Bioinformatics 19 03-08-2016 08:59 AM
NexteraXT MiSeq Bowtie2 strange insert size distribution M4TTN Bioinformatics 63 04-08-2015 06:51 AM
custom ssDNA library with strange size distribution of qPCR products pyridine Sample Prep / Library Generation 4 07-17-2013 03:18 AM
Periodic variation in nucleotide distribution along the read & other strange things lionelguy Illumina/Solexa 16 08-03-2011 09:04 AM
[Galaxy] Strange QC Nucleotides Distribution Chart zippered_ohio Bioinformatics 5 06-28-2011 11:43 PM

Reply
 
Thread Tools
Old 10-06-2015, 09:13 AM   #1
sacketlc
Junior Member
 
Location: Washington DC

Join Date: Feb 2015
Posts: 1
Default Strange Fst distribution in BayeScan

I have used BayeScan (2.1) successfully several times on various subsets of a large SNP dataset that I am currently working on. Recently, I tried to use it to infer outliers between several population pairs using my whole dataset (57,000 SNPs), and came up with a peculiar result in all 9 population pairs. The plot of posterior density for Fst looks as expected, except with a mean Fst ~ 0.4 - 0.6, depending on the population pair (not ~0.05 as I know them to be). If I plot a frequency distribution of Fst values inferred by BayeScan, the distribution is the opposite of what I would expect -- the highest frequency is Fst at these high values around 0.5, with a tail on the left for low Fst values. These are bird populations that we know are differentiated but still experience some gene flow, and are geographically quite close. Calculating Fst in vcftools or Genepop gives a mean and distribution more in line with expectations -- mean ~0.05 with a tail to the right of higher Fst values.

Has anyone else experienced something similar? I have been picking my brain to try to figure out what is causing this to happen, but don't have any hypotheses yet. The current dataset has less missing data than subsets I have previously used, so I wouldn't expect that to be an explanation, and I used the same approach as previously. I used PGDSpider to convert from vcf format to BayeScan format.

I would appreciate any thoughts!
sacketlc is offline   Reply With Quote
Old 10-06-2015, 09:35 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,059
Default

For reference

cross-posted on Biostars: https://www.biostars.org/p/160738/
GenoMax is offline   Reply With Quote
Reply

Tags
bayescan, fst outliers, snp

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 05:41 PM.


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