GWAS dataset calculating MAF for European dataset
1
0
Entering edit mode
2.2 years ago
v.johnson ▴ 20

How does frequency of existing variant in the combined 1000 genome European data from Variant effector predicator website equate to MAF. I have a SNP data from GWAS data and need to find the MAF for each of the SNPs. Thank you

MAF GWAS • 690 views
ADD COMMENT
0
Entering edit mode
2.2 years ago

Hi Victoria, can you please clarify the format of your data, currently? Per any given variant, MAF can be applied specifically to any population, i.e., MAF in the global population, MAF in the Irish population, MAF among sub-Saharan Arficans, et cetera.

ADD COMMENT
0
Entering edit mode

Dear Kevin, I have a European dataset with SNP data and need to get the MAF for each SNP. I have a lot of data from a GWAM study. Headings are; CHR SNP BP ref_allele alt_allele sample_size case_proportion Thank you

ADD REPLY
0
Entering edit mode

Hi, in that case, you already have enough information to calculate the MAF in your own cohort of patients / samples. However, if you want the MAF for your variants from the large public databases, then you could use Ensembl VEP or ANNOVAR. Both of these will take some time to set up, though. There are other solutions, too.

ADD REPLY
0
Entering edit mode

As in;

If (AF < 0.5) then MAF = AF

if (AF > 0.5) then MAF = 1.0 - AF

ADD REPLY
0
Entering edit mode

More or less, yes! That is, this is your study-specific MAF.

ADD REPLY
1
Entering edit mode

Thank you for your help!

ADD REPLY
0
Entering edit mode

If you find that the AF is > 0.5, though, then you also have to remember to flip the order of the REF and ALT alleles. These situations are quite common and may arise due to this problem: A: Alternate nucleotide is more frequent than reference nucleotide. OMG I'm dizzy.

I would also verify that some of the AFs in your VCF are actually valid by literally counting across all samples... got to be sure in bioinformatics!

ADD REPLY

Login before adding your answer.

Traffic: 864 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6