CADD-SV plugin not working

[kittysher]

CADD-SV plugin not working

[nakib103]

Hello @kittysher ,

Thanks for your query!

I can confirm that I can re-produce the issue. We are looking more into it and let you know once we have more information.

Best regards,
Nakib

[nakib103]

Hi @kittysher,

Thanks for your patience. Here is why your data is not getting annotated -
CADD plugin matches SV using position and variant type as these are data available from the CADD annotation file. Currently supported type is INS, DUP, and DEL. Check the CADD-SV website for details -
https://kircherlab.bihealth.org/download/CADD-SV/v1.1/column-annotations.txt (See Type column)

Looking more at you input file -

• the first line have a typo - there is a new line after the REF allele. Even after the fix it won't get annotated as it is not an SV and won't get any match from the annotation file.
• the second line variant type is BND and so would not get annotated.
• the third line I don't see any matching variant in the CADD annotation file -

$ tabix -D https://kircherlab.bihealth.org/download/CADD-SV/v1.1/1000G_phase3_SVs.tsv.gz 1:43221167-43221167 | wc -l
0

Hope it helps answer your question. Let me know if you have any further query.

Best regards
Nakib

[kittysher]

Hi Nakib,

Thank you for getting back to me about this issue.
Sorry about the typo in the first line of the input file I gave - that was due to me typing it incorrectly here, the VCF itself doesn't have this. On this page, it says that breakend annotation is available in VEP v110:
https://www.ensembl.info/2023/07/21/cool-stuff-ensembl-vep-can-do-enhanced-structural-variant-annotation/
Please could you let me know if this is the case?

Many thanks
Kitty

[nakib103]

Hello @kittysher,

Yes, that is the case. That is why you will see the breakend variant getting annotated in the output with things such as transcript affected and consequence on that transcript -

manta_MantaBND:1693:0:1:0:0:0:0	chr1:19031896	[chr3:122974046[C	ENSG00000082684	ENST00000650207	Transcript	feature_truncation ...

But annotation from CADD plugin is different matter. Plugins are extension to Ensembl VEP and often depends on files from external resources. As I pointed out CADD annotation file (which is an external resource) currently supports INS, DUP, and DEL. From https://kircherlab.bihealth.org/download/CADD-SV/v1.1/column-annotations.txt file -

Chrom	Chromosome
Start	Start Coordinates
End	End Coordinates
Type	DELetion, INSertion, DUPlication. <- see here
....

You can learn more about Ensembl VEP plugins from here - https://www.ensembl.org/info/docs/tools/vep/script/vep_plugins.html

Best regards,
Nakib

[kittysher]

Hi Nakib,

OK thank you for clarifying the difference between the VEP and CADD annotations.
I am still not getting the CADD VEP plugin to work though - the output is empty for every variant in my VCF which I don't think is correct. Please could you explain what this error message means? Is there some incompatibility between my VCF and the reference CADD SV file?

WARNING:` Plugin 'CADD' went wrong:
-------------------- EXCEPTION --------------------
MSG: Missing reference or alternate sequence

Many thanks
Kitty

[nakib103]

Hi @kittysher,

The warning is given when the alternate allele type is not supported. In the example you provided it was the BND type that was not supported and is generating the warning.

It could be there are no variants in your input files that gets matched with CADD SV file that are of supported types. Can you provide any other variants that is not annotated but should be?

Best regards,
Nakib

[nakib103]

Hello @kittysher,
I am closing this ticket. If you have further problem feel free to open another ticket.

Best regards,
Nakib