CADD can not get any scores

[zyxNo1]

CADD can not get any scores

[zyxNo1]

YES!!

Use correct path to specify --dir_plugins, then the CADD scores are annotated without any warnings.

Many thanks for your time and consideration. @olaaustine

[olaaustine]

Hi @zyxNo1,
Thank you very much for opening this query.
To debug this issue, can you share more examples of your input variants, what version of VEP plugin you are using and also what version of CADD file you are using?
Also in your command, I can see you are using --cache, can you tell us what cache version you are using also
Thank you
OIa.

[zyxNo1]

Thanks for your kindly reply.
The cache version is consistent with VEP, also being 105. The CADD file used is the latest one, which corresponds to version 111.
I have attached partial input and output and hope they can be helpful.

input.vcf.gz

variant_effect_output.txt

[olaaustine]

Hi @zyxNo1,
Thank you for your response, I have not been able to recreate the issue using that input file, I was able to get CADD scores.
To understand whats happening, can you confirm that you have tabix installed in your path and you also have the tabix files of ~/resource/ensembl-vep/cadd/whole_genome_SNVs.tsv.gz and ~/resource/ensembl-vep/cadd/gnomad.genomes.r4.0.indel.tsv.gz in the same directory, if it is not, you can download it from the same place you downloaded those files.

Secondly, can you you run the command above like this
vep -i ~/WGS.bwa.dedup-NC_T_1_vs_NC_N_1-MuTect2.vcf --plugin CADD,~/resource/ensembl-vep/cadd/whole_genome_SNVs.tsv.gz,~/resource/ensembl-vep/cadd/gnomad.genomes.r4.0.indel.tsv.gz --cache --force_overwrite --fork 10 --dir_cache <directory where your cache is> --offline
Let us know if this works.
Thank you
Ola.

[zyxNo1]

I have just comfired that tabix tool is in the system path (version 1.11), and four required files are located in the same directory:
gnomad.genomes.r4.0.indel.tsv.gz gnomad.genomes.r4.0.indel.tsv.gz.tbi whole_genome_SNVs.tsv.gz whole_genome_SNVs.tsv.gz.tbi

Then I attempted to run VEP with the following command:
vep -i ~/input.vcf --plugin CADD,~/resource/ensembl-vep/cadd/whole_genome_SNVs.tsv.gz,~/resource/ensembl-vep/cadd/gnomad.genomes.r4.0.indel.tsv.gz --cache --force_overwrite --fork 10 --dir_cache ~/.vep/ --offline

But still got no scores.

Moreover, I can confirm that VEP’s basic annotation functionality is working as expected, as well as “custom” annotations like ClinVar information. The plugin CADD ouput generates a notably large warning file of uninitialized values; while input.vcf is 2.5 M but the warning file is about 45 M. I don't know if this information can be of any help.

[olaaustine]

Hi @zyxNo1,
Please can you specify the --dir_plugins to use also and the --assembly flag.
Can you let us know if this resolves the problem?
Thank you
Ola