Question: How to identify the ancestral population to which an individual belongs knowing its allele frequencies?
gravatar for cyril-cros
4.7 years ago by
cyril-cros900 wrote:

Hi, I have been given as part as a homework a set of genotyped human SNP for many individuals; I am supposed to identify their ancestral populations (disclaimer: as a bonus question I won't submit and which is likely an example from the readings I have to do, so I am not asking for a key in hand solution). Using STRUCTURE, I can get an estimate of the number of ancestral populations of my sample and their allele frequencies.

My SNPs are coded by 0 or 1 (ancestral or derived states), but I know their labels (rsXXXXX). Some SNP are found in one state at a much higher frequency in some population, which makes them informative: I can can manually crosscheck my frequency estimates with the population frequencies found in HapMap.

For example, with SNP rs924201 , populations A and C present the derived state in around 50% of cases and B the ancestral in more than 80%. HapMap tells me that ~80% of Africans share a same variant of this SNP , but only roughly 50% of European and Asians. I can guess that B would likely be from Africa, given the fact that my sample is rather large.

Is there a way to fit my imputed allele frequencies to real world frequencies? Like, finding all markers with noticeably different frequencies between my predicted groups and matching them to real world populations. I would need a way to retrieve the frequencies of the most common variants for each of the HapMap populations. I could maybe then try some kind of lasso on the frequencies of the most common allele, in order to know which real world population is the closest of my predicted ancestral population.

Thanks for your advice!

snp structure infer ancestry • 2.0k views
ADD COMMENTlink modified 4.6 years ago by stolarek.ir650 • written 4.7 years ago by cyril-cros900
gravatar for
4.6 years ago by
stolarek.ir650 wrote:

Read on. STRUCTURE is all you need for homework. You could also use somewhat visual way with PCA (eigensoft)

ADD COMMENTlink written 4.6 years ago by stolarek.ir650
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1698 users visited in the last hour