[ensembl-dev] Annotating variants against custom transcripts
Anja Thormann
anja at ebi.ac.uk
Thu Jan 18 15:11:00 GMT 2018
Hi Shawn,
the Bio::DB::HTS module should have been installed by INSTALL.pl. Did you set the DYLD_LIBRARY_PATH environment variable as asked during the installation?
You can also install Bio::DB::HTS manually: https://github.com/Ensembl/Bio-DB-HTS <https://github.com/Ensembl/Bio-DB-HTS>.
Anja
> On 18 Jan 2018, at 14:46, Shawn Yost <yostshawn at gmail.com> wrote:
>
> Hi,
> I've installed the latest version of vep and I'm not getting a new
> error. This error occurred also during the install/test step of vep:
>
> -------------------- EXCEPTION --------------------
> MSG: ERROR: Cannot use format gff without Bio::DB::HTS::Tabix module installed
>
> STACK Bio::EnsEMBL::VEP::AnnotationSource::File::new
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/AnnotationSource/File.pm:162
> STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all_custom
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:228
> STACK Bio::EnsEMBL::VEP::AnnotationSourceAdaptor::get_all
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/AnnotationSourceAdaptor.pm:93
> STACK Bio::EnsEMBL::VEP::BaseRunner::get_all_AnnotationSources
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/BaseRunner.pm:175
> STACK Bio::EnsEMBL::VEP::Runner::init
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/Runner.pm:123
> STACK Bio::EnsEMBL::VEP::Runner::next_output_line
> /mnt/scratch/DGE/GENSUSC/syost/vep/ensembl-vep/modules/Bio/EnsEMBL/VEP/Runner.pm:356
> STACK toplevel ./t/Runner.t:739
> Date (localtime) = Thu Jan 18 14:26:33 2018
> Ensembl API version = 91
>
>
>
> I'm using a conda env to install/run VEP. Inside of the
> 'ensembl-vep/' directory Tabix does exist:
>
>> ls -la Bio/DB/HTS/Tabix*
> -rw-rw---- 1 syost cancgene 6721 Jan 18 14:24 Bio/DB/HTS/Tabix.pm
>
> Bio/DB/HTS/Tabix:
> total 12
> drwxrwx--- 2 syost cancgene 4096 Jan 18 14:24 .
> drwxrwx--- 5 syost cancgene 4096 Jan 18 14:24 ..
> -rw-rw---- 1 syost cancgene 2974 Jan 18 14:24 Iterator.pm
>
>
> I've also installed Tabix separately and it works. Do you have any
> suggestions? Why isn't it recognizing the Tabix.pm?
>
>
> Thanks,
> Shawn
>
>
>
>
> On Tue, Jan 16, 2018 at 2:21 PM, Anja Thormann <anja at ebi.ac.uk> wrote:
>> Hi Shawn,
>>
>> I noticed that you are not using the supported VEP code. You can install the
>> new code by following the instructions here:
>> http://www.ensembl.org/info/docs/tools/vep/script/vep_download.html
>>
>> The new VEP code supports annotations against a GFF or GTF file:
>> http://www.ensembl.org/info/docs/tools/vep/script/vep_options.html#other
>>
>> Anja
>>
>>
>> On 16 Jan 2018, at 14:06, Shawn Yost <yostshawn at gmail.com> wrote:
>>
>> Hi,
>> I would like to annotate my VCF file against a custom transcript database.
>> I've created both a GFF v3 file and a GTF file (see below) and I have been
>> unsuccessful in getting VEP to annotate against these transcripts. The
>> examples below are before using sort + bgzip + tabix (so that is not the
>> problem).
>>
>> I'm currently running VEP v85.
>>
>> The command I used was:
>> variant_effect_predictor.pl -i IN.vcf --custom test.gff.gz,,gff -fasta
>> test.fa --cache -o OUT -dir /path/to/cache --hgvs --cache_version 75
>> --offline --force_overwrite
>>
>> In the outputted file I only see ENSTs and can't find the transcripts I
>> inputted along with them. The same thing occurs if I run --custom
>> test.gtf.gz,,gtf. If I change the command to --custom
>> test.gtf.gz,,gtf,overlap it will tell me if it overlaps the inputted
>> transcript but it doesn't annotate against the transcript.
>>
>>
>> Is there a problem with the command options I am using? Is there a problem
>> with the inputted files? How do I get VEP to annotate the variant against
>> my custom transcript (i.e. 412662 2:73613056 A GENE1
>> TRANSCRIPT1 Transcript synonymous_variant ....)?
>>
>>
>>
>> Example output:
>> 241004 2:73613032-73613049 - ENSG00000116127 ENST00000264448
>> Transcript inframe_deletion 147-164 36-53 12-18
>> LEEEEEE/L ctGGAGGAGGAGGAGGAGGAg/ctg -
>> IMPACT=MODERATE;STRAND=1;HGVSc=ENST00000264448.6:c.36_53delNNNNNNNNNNNNNNNNNN;HGVSp=ENSP00000264448.6:p.Glu23_Glu28del
>> 412662 2:73613056 A ENSG00000116127 ENST00000377715 Transcript
>> synonymous_variant 171 60 20 E gaG/gaA -
>> IMPACT=LOW;STRAND=1;HGVSc=ENST00000377715.1:c.60N>A;HGVSp=ENST00000377715.1:c.60N>A(p.%3D)
>> 412662 2:73613056 A ENSG00000116127 ENST00000409009 Transcript
>> synonymous_variant 171 60 20 E gaG/gaA -
>> IMPACT=LOW;STRAND=1;HGVSc=ENST00000409009.1:c.60N>A;HGVSp=ENST00000409009.1:c.60N>A(p.%3D)
>> 412662 2:73613056 A ENSG00000116127 ENST00000264448 Transcript
>> synonymous_variant 171 60 20 E gaG/gaA -
>> IMPACT=LOW;STRAND=1;HGVSc=ENST00000264448.6:c.60N>A;HGVSp=ENST00000264448.6:c.60N>A(p.%3D)
>> 402364 2:73613066-73613071 - ENSG00000116127 ENST00000377715
>> Transcript inframe_deletion 181-186 70-75 24-25 EE/-
>> GAGGAA/- -
>> IMPACT=MODERATE;STRAND=1;HGVSc=ENST00000377715.1:c.70_75delNNNNNN;HGVSp=ENSP00000366944.1:p.Glu27_Glu28del
>>
>>
>>
>> GFF v3 file:
>> 15 . transcript 74701625 74726300 . -
>> . ID=TRANSCRIPT1;Alias=10741;Name=SEMA7A
>> 15 . exon 74726082 74726300 . - .
>> ID=EXON37A10411.1;Parent=TRANSCRIPT1
>> 15 . exon 74711142 74711293 . - .
>> ID=EXON37A10411.2;Parent=TRANSCRIPT1
>> 15 . exon 74710609 74710650 . - .
>> ID=EXON37A10411.3;Parent=TRANSCRIPT1
>> 15 . exon 74710218 74710310 . - .
>> ID=EXON37A10411.4;Parent=TRANSCRIPT1
>> 15 . exon 74709932 74710016 . - .
>> ID=EXON37A10411.5;Parent=TRANSCRIPT1
>> 15 . exon 74709676 74709786 . - .
>> ID=EXON37A10411.6;Parent=TRANSCRIPT1
>> 15 . exon 74708916 74709055 . - .
>> ID=EXON37A10411.7;Parent=TRANSCRIPT1
>> 15 . exon 74708142 74708326 . - .
>> ID=EXON37A10411.8;Parent=TRANSCRIPT1
>> 15 . exon 74707179 74707287 . - .
>> ID=EXON37A10411.9;Parent=TRANSCRIPT1
>> 15 . exon 74706888 74707086 . - .
>> ID=EXON37A10411.10;Parent=TRANSCRIPT1
>> 15 . exon 74704226 74704353 . - .
>> ID=EXON37A10411.11;Parent=TRANSCRIPT1
>> 15 . exon 74703897 74704051 . - .
>> ID=EXON37A10411.12;Parent=TRANSCRIPT1
>> 15 . exon 74703636 74703697 . - .
>> ID=EXON37A10411.13;Parent=TRANSCRIPT1
>> 15 . exon 74701625 74703326 . - .
>> ID=EXON37A10411.14;Parent=TRANSCRIPT1
>> 15 . CDS 74726082 74726259 . - 0
>> ID=CDS37A10411.1;Parent=TRANSCRIPT1
>> 15 . CDS 74711142 74711293 . - 2
>> ID=CDS37A10411.2;Parent=TRANSCRIPT1
>> 15 . CDS 74710609 74710650 . - 0
>> ID=CDS37A10411.3;Parent=TRANSCRIPT1
>> 15 . CDS 74710218 74710310 . - 0
>> ID=CDS37A10411.4;Parent=TRANSCRIPT1
>> 15 . CDS 74709932 74710016 . - 0
>> ID=CDS37A10411.5;Parent=TRANSCRIPT1
>> 15 . CDS 74709676 74709786 . - 2
>> ID=CDS37A10411.6;Parent=TRANSCRIPT1
>> 15 . CDS 74708916 74709055 . - 2
>> ID=CDS37A10411.7;Parent=TRANSCRIPT1
>> 15 . CDS 74708142 74708326 . - 0
>> ID=CDS37A10411.8;Parent=TRANSCRIPT1
>> 15 . CDS 74707179 74707287 . - 1
>> ID=CDS37A10411.9;Parent=TRANSCRIPT1
>> 15 . CDS 74706888 74707086 . - 0
>> ID=CDS37A10411.10;Parent=TRANSCRIPT1
>> 15 . CDS 74704226 74704353 . - 2
>> ID=CDS37A10411.11;Parent=TRANSCRIPT1
>> 15 . CDS 74703897 74704051 . - 0
>> ID=CDS37A10411.12;Parent=TRANSCRIPT1
>> 15 . CDS 74703636 74703697 . - 1
>> ID=CDS37A10411.13;Parent=TRANSCRIPT1
>> 15 . CDS 74702965 74703326 . - 2
>> ID=CDS37A10411.14;Parent=TRANSCRIPT1
>>
>>
>> GTF file:
>> TRANSCRIPT1 15 - 74701624 74726300 74702964
>> 74726259 14
>> 74701624,74703635,74703896,74704225,74706887,74707178,74708141,74708915,74709675,74709931,74710217,74710608,74711141,74726081,
>> 74703326,74703697,74704051,74704353,74707086,74707287,74708326,74709055,74709786,74710016,74710310,74710650,74711293,74726300,
>> 0 HGNC:10741 cmplcmpl 1,2,0,1,0,2,0,1,1,0,0,0,1,0,
>> TRANSCRIPT2 17 - 8108048 8113944 8108188 8113542 9
>> 8108048,8108533,8109808,8110067,8110493,8110888,8111055,8113494,8113847,
>> 8108362,8108708,8109957,8110206,8110685,8110943,8111158,8113567,8113944,
>> 0 HGNC:11390 cmpl cmpl0,2,0,2,2,1,0,0,-1,
>> TRANSCRIPT3 11 - 73711325 73720282 73712456
>> 73718087 7
>> 73711325,73714871,73715528,73716774,73717213,73717961,73720022,
>> 73712571,73715052,73715630,73716978,73717424,73718182,73720282, 0
>> HGNC:12519 cmplcmpl 2,1,1,1,0,0,-1,
>> TRANSCRIPT4 11 - 123594634 123612391 123596704
>> 123601596 9
>> 123594634,123598183,123598840,123599833,123600322,123601194,123610824,123611124,123612256,
>> 123597699,123598303,123598970,123599922,123600533,123601693,123610900,123611244,123612391,
>> 0 HGNC:12994 cmpl cmpl 1,1,0,1,0,0,-1,-1,-1,
>> TRANSCRIPT5 19 - 51994480 52005043 51994894
>> 52004987 8
>> 51994480,52000133,52000602,52001271,52002423,52002691,52003173,52004560,
>> 51995083,52000230,52000699,52001541,52002471,52002970,52003554,52005043,
>> 0 HGNC:15482 cmpl cmpl 0,2,1,1,1,1,1,0,
>> TRANSCRIPT6 1 + 228194722 228248972 228194829
>> 228247166 4 228194722,228210367,228238356,228246686,
>> 228194900,228210609,228238622,228248972, 0 HGNC:15983 cmpl
>> cmpl 0,2,1,0,
>>
>>
>>
>>
>>
>>
>>
>> Thank you for your help,
>> Shawn
>>
>> _______________________________________________
>> Dev mailing list Dev at ensembl.org
>> Posting guidelines and subscribe/unsubscribe info:
>> http://lists.ensembl.org/mailman/listinfo/dev
>> Ensembl Blog: http://www.ensembl.info/
>>
>>
>>
>> _______________________________________________
>> Dev mailing list Dev at ensembl.org
>> Posting guidelines and subscribe/unsubscribe info:
>> http://lists.ensembl.org/mailman/listinfo/dev
>> Ensembl Blog: http://www.ensembl.info/
>>
> _______________________________________________
> Dev mailing list Dev at ensembl.org
> Posting guidelines and subscribe/unsubscribe info: http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20180118/71ed7149/attachment.html>
More information about the Dev
mailing list