Saturday, 20 October 2007

genetics - How to compare SNP from genotyping results for multiple people with a known phenotype?

Each rsid identifies a unique SNP in the genome. Thus there should not be any entry in the files that have the same rsid but different chromosomes and/or positions. If you do find this, it is likely that you have data from different versions of the assembled genomes.



To find genetic variants associated, i.e. correlated, with your trait, you need to focus on the genotype at each of the SNPs. For example, SNP rs3094315 is an A/G polymorphism, i.e. an individual can either have the genotype AA, AG, or GG. To find if it is associated with tongue rolling, you would count the number of AA, AG, and GG individuals (alternatively you could sum the number of A and G alleles) in the group that can roll their tongue and compare those numbers to the group that cannot. You would then repeat this for each SNP that you want to test. For most SNPs, there will be no difference in proportion.



That being said, this is not a modest undertaking as this a cornerstone of human and statistical genetics. There are many issues that could strongly bias your results (e.g. population stratification, genotyping error, linkage disequilibrium, multiple testing). To learn more, you can read about the Cochrane-Armitage Trend Test and the PLINK software. If you are really serious I recommend Applied Statistical Genetics with R by Andrea S. Foulkes.

No comments:

Post a Comment