get allele frequency (MAF) value from snp

Question

[Barista]

Hi all! :)

I have a list of snp ids, to which I need their correspodning MAF (minor allele frequency) values.

I tried to do it in R using biomaRt but it seems that it is not possible- correct me if I am wrong.

The only way I managed to do it is in Python using the API from ensembl, but since the list is over half a milion of snps, it takes so much time and I always get connection errors or something, so I cannot make this finish even once :(

Please help me and give me some tips, I will highly appreciate that, because it is my first time working on such a task.

[Pierre Lindenbaum]

download a VCF containing the frequencies (dbsnp, etc....) and use bcftools to get the data in a set of regions

[Ben Moore]

In addition to Pierre's answer, the Ensembl Variant Effect Predictor (VEP) might be useful here as you can retrieve allele frequencies (as well a slots of other data) for a custom list of variants:

https://www.ensembl.org/info/docs/tools/vep/index.html

[Barista]

Where can I find such a VCF file? I am currently trying to use dbsnp API in Python but with no success so far.

[Pierre Lindenbaum]

gnomad see the download section : https://gnomad.broadinstitute.org/downloads

[Barista]

Does it really take so long to download one file? I got ..teen hours for one chromosome file (and I want several). Am I doing something wrong, is there a faster way? Highly appreciate your help in advance!