Side note: Canonical isn't always what people may be looking for. For example, the longest transcript for ESR1 and BRCA1 genes are not the most popular ones. Canonical is mostly but not always what we are looking for. Canonical translates to "longest", not "isoform 1 on NCBI". For that, you may want to use the MANE + RefSeq Select.
Annotations vep : keep high impact in any transcript + canonical
Hello all,
I use vep to annotate my vcf. I would like to know if there is a way, if a variants as a high impact in a non canonical transcript, to keep it and the annotations in the canonical. I don't want to annotate any impact in any transcript to avoid huge amount of data, but i don't want to miss a high impact in alternative transcript.
• 4,271 views
•
link
updated
by
Ram
•
written
by
quentin54520 •
•
written
by
quentin54520 •
1 answer
Hi quentin54520,
I would advise annotating your variants with the VEP without any filters, then using the filter_vep tool to filter your VEP output: https://www.ensembl.org/info/docs/tools/vep/script/vep_filter.html#filter_run
You can write custom filters to return lines that refer to either the canonical transcript or have a particular IMPACT rating.
This is very true, Ram. Here is the Ensembl documentation on canonical transcript selection: https://www.ensembl.org/info/genome/genebuild/canonical.html
Thank you, Ben. Another super-useful feature is being able to assign priority to flag-pick-allele. Use the option to flag an entry based on one's own priorities: https://useast.ensembl.org/info/docs/tools/vep/script/vep_options.html#opt_pick_order
The option to flag something helps us get VEP's take but ask VEP to not throw anything away on its own discretion.
Yes, indeed it is very difficult to know what a canonical transcript is, there can be several definitions depending on the people and the context. I work on rare diseases, and I have to annotate human genomes. Usually I use the --pick option but I'm always afraid to miss a potentially interesting variant. But on the other hand I don't want to end up with ten annotations per variant, there are already so many variants. Do you think there would be a better solution than the --pick option? Using the --mane option and then filtering to keep only what has the mane flag?
• 0 views
•
link
written
by
quentin54520 •
Use --flag options instead of plain --pick and get as much annotation info as you can. Perform filtering later once you have determined the optimal filters to remove noise while retaining important information. That's the most advice we can give you without getting into the particulars of your requirements.
Log in to answer this question.
More posts like this
-
Annotate VCF via VEP for gnomAD annotation considering FILTER column
written by asalimih •
-
How to choose the right annotation from multiple annotation of single variant from snpeff result
written by sonsunjirachote •
-
VEP and loftee. How to get the GT value
written by Manuel •
-
Vcf stat about consequences
written by quentin54520 •
-
Vcf : remove variant absent in affected sample
written by quentin54520 •
-
keep samples that carry non-reference allele in list of variants snpeff/snpsift
written by curious
-
Vep : MNP phased genotype
written by quentin54520 •
-
How does VEP decide on canonical transcripts and is there a list?
written by hagen.klett •
-
How to select predefine transcript when annotate VCF
written by Paul
-
Wrong predict consequence in variants using VEP Ensembl
written by theodorobiotec •
Please note:
I've fixed these in your post and title.