This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparison of allele frequencies between population

Hi all,
I would like to examine whether the frequencies of many alleles differ between the African and North European populations.
I downloaded the 1000 genome project wgs vcf file via the FTP site: http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/ALL.wgs.phase3_shapeit2_mvncall_integrated_v5b.20130502.sites.vcf.gz
But the allele frequencies are given for super populations, thus the European population includes Spain and Italy.
Does anyone familiar with this kind of data and could help?

Thanks

snp vcf 1000 genome project population
• 2,029 views • link
updated by chrchang523 chrchang523 • written by elisheva elisheva

1 answer

  1. Download the pedigree file from https://www.internationalgenome.org/faq/can-i-get-phenotype-gender-and-family-relationship-information-samples/ .
  2. Define your European subpopulation, and create a file with just the sample IDs within that subpopulation with e.g. a short shell script.
  3. Compute allele frequencies using just those samples. This can be done efficiently with e.g. plink2 --vcf <VCF file> --keep <sample ID file> --freq.
• 0 views • link
written by chrchang523 chrchang523

Thanks for your response.
It seems that the versions of the vcf file and the pedigree file are different.
I couldn't find the matching version for both files.
Is there is any option to solve this?

• 0 views • link
written by elisheva elisheva

The pedigree file contains some samples that were excluded from the phase 3 VCF, but that's okay, the extra sample IDs will just be ignored by --keep.

• 0 views • link
written by chrchang523 chrchang523

Apparently, there are no sample Ids in the vcf file I posted, and I couldn't find any others.

• 0 views • link
written by elisheva elisheva

Log in to answer this question.

More posts like this