Hello everyone!
I am trying to get the gnomad allele frequencies of variants in the VCF files (I share the screenshot of one example which allele frequency I want to have). I used the following protocols and operations:
protocols="gnomad211_genome,gnomad211_exome,clinvar_20221231,dbnsfp42a,avsnp150,refGene" operations="f,f,f,f,f,g"
I have two sets of AF values, I believe one retrieved from gnomad211_genome database, one retrieved from gnomad211_exome database. Here are the columns I get after the annotation: (2 of the each one with the same name: 2 AF, 2 AF_popmax, 2 AF_male... etc.)
AF AF_popmax AF_male AF_female AF_raw AF_afr AF_sas AF_amr AF_eas AF_nfe AF_fin AF_asj AF_oth non_topmed_AF_popmax non_neuro_AF_popmax non_cancer_AF_popmax controls_AF_popmax
But any of these columns have the value that I want (which can be seen in the screenshot I share) How can I get gnomad allele frequencies of the variants?
Also, another weird thing is that some variants have an AF values that are too much low, when I check them on gnomad database, they don't exist. How can I get some frequency values if the variants are not in the gnomad database, should I trust these too low values?
Thank you all!
Notice that the allele frequency is 0.2391
before we dive in too deep let's confirm we are on the same page that gnomad 2.1.1 (in your text) is in hg37. gnomad 3.1.2 (in your screenshot) in is hg38. annovar may be lifting the 2.1.1 over to hg38 but it won't match 3.1.2.
Thanks for pointing out. Let me explain:
I had VCF files with gnomad 2.1.1 variants (chromosome positions). I performed a liftover with Pickard liftover tool which means I converted the chromosome positions of the variants from hg19 to hg38. And with these new VCF files (containing hg38 chromosome positions) I performed annotation with ANNOVAR. But my supervisor suggested me to use specifically protocols="gnomad211_genome" and "gnomad211_exome" for the annotation because they are studied more and can include more individuals. Annotation went well, but I have this kind of variants with too little frequency values. In the screenshot I just gave an example, and the fact that it is hg38 or hg19 does not matter because I tried to look for these variants (the ones with too little frequency values) on gnomad, both for hg19 and hg38 chromosome positions and I could not find, they don't exist. My questions are briefly:
You performed a liftover on your VCF and then used a liftover-gnomad211 for the annotation? Why not simply use the regular gnomad211 with your hg19 callset? I use the gnomad liftover because my calls are done using hg38 and gnomAD 3.x is not as robust as 2.x yet, but you have no reason to do a double liftover.
Thank you so much, I decided to move forward as you suggested. Initially, I thought using the most up-to-date databases and positions would provide a better quality annotation, but as you said it was not healthy to do the liftover and work on that.
There are plenty of variants that don't exist in gnomAD - only 759M SNPs (chr-pos-ref-alt) out of a possible 9B. Your question is a bit confusing to me because I don't know how you are getting any frequency values for variants that are not in gnomAD. They should just be NULL.
Yes, I was confused as well. I think the same. Maybe Annovar uses some other resources besides gnomad, but it does not make sense. Also, I get these columns with the annotation I performed (I provide an example). Can you tell me what they stand for:
I am asking this, because it seems that these columns represent classical annotation columns Chr Start End Ref Alt ... But they appear in the end of the annotation as an addition. Chr Start positions differ from the first column.
Here are the protocols I used:
Thank you.
0.60625 is 97/160 so if you have 80 individuals in your VCF that's an internal allele frequency, not gnomAD