I have some ADMIXTURE results that I am struggling to plot, but I am interested in determining what SNPs are going into which components. I know that ADMIXTURE itself does not specify which components correspond to which SNPs. I was wondering if anyone has any suggestions on how I could do this.
I am looking into ELAI, LAMP, and HAPMIX, but I am open to suggestions.
Thanks!
1 answer
You don't really have SNPs associated with a particular ancestral component but rather, an estimate of allele frequency in each ancestral component. That's the P file. See the manual page 3:
https://www.genetics.ucla.edu/software/admixture/admixture-manual.pdf
In that case, many of the programs you listed are a really good place to start, essentially you are asking about local ancestry estimation, for which ELAI, LAMP, HAPMIX, chromosome painter etc. are good alternatives.
There also programs that use local ancestry to increase power for association testing if that is of interest to you, in this case maybe start by reading http://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1002325
The same author wrote a fairly good review in 2013 .. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3556814/
I remember Sriram Sankararaman did something along those lines for the Altai Neanderthal paper. Have a look at this maybe:
He essentially used an HMM to identify these chunks in phased genomes.
•
written
by
Gabriel R.
Ah then have a look at ChromoPainter and fineStructure
•
written
by
Gabriel R.
Log in to answer this question.
-
Optimizing mutect/sarek for somatic substitution variant calling
written by ramiro.barrantes •
-
Upset plot in ComplexHeatmap
written by Omics data mining
-
Plotting the results (p-values) from various conditional GWAS
written by sam.far.7 •
-
Where can I upload SNP array data for public accessibility?
written by devenvyas
-
Creating a tree from pairwise FST values for ADMIXTURE components
written by devenvyas
-
Generating an ADMIXTURE plot that stacks the different values of K using R
written by devenvyas
-
Perplexing ADMIXTURE results - one component present at K =10 disappears at 11, reappears at 12
written by devenvyas
-
Looking for a quick way of consistently order components between different ADMIXTURE analyses
written by devenvyas
-
Setting dimensions of barplot area in R (for STRUCTURE-like plot)
written by devenvyas
-
Help making ADMIXTURE output digestible
written by devenvyas
