Question: Structural Variations (Sv) Annotation Workflow
5
gravatar for Doctoroots
7.3 years ago by
Doctoroots780
Doctoroots780 wrote:

I want to know if there is some kind of pipeline i can use in order to annotate a given list of structural variants detected by NGS.

specifically, i just ran Breakdancer on a paired end sequencing data and recieved the output specifying the different types of predicted variations (insertions / deletions / inversions / translocations) example :

1 1413910 9+0- 1 1452528 0+10- DEL 38636 99

1 1588571 0+17- 1 1653965 17+0- ITX 64960 99

the columns are : chr1, position1, supportingreads(forward/reverse), chr2, position2, supportingreads(forward/reverse), variationtype, variationsize, quality

i basically want to annotate these variants in regards to gene name, exon/intron, transcription binding sites, ncRNA (miRNA) and any other known feature regarding the genetic region affected by them.

is there a known tool for this? if not, i would appreciate any suggestions on how to perform such a task (for example some sort of outline for a workflow i can build)

thanks.

annotation workflow structural • 4.9k views
ADD COMMENTlink modified 7.3 years ago by Raymond301140 • written 7.3 years ago by Doctoroots780
1

How did you end up annotating your SV data. I find myself in the same kind of situation.

ADD REPLYlink written 6.5 years ago by William4.4k

SV annotation (with OMIM, DGV, 1000g, haploinsufficiency, TAD, ... and also with your own in-house information) can be easily automated !

You can look at this post describing the annotSV tool: Annotation for SV and CNV

ADD REPLYlink modified 10 months ago • written 10 months ago by LGMgeo90
1
gravatar for Bioscientist
7.2 years ago by
Bioscientist1.7k
Bioscientist1.7k wrote:

Bedtools may help you

ADD COMMENTlink written 7.2 years ago by Bioscientist1.7k
1
gravatar for Raymond301
7.2 years ago by
Raymond301140
Rochester, MN
Raymond301140 wrote:

There happens to be a complete pipeline that results in IGV visualization as well as flat files for tertiary analysis.

Look at this publication:

http://bioinformatics.oxfordjournals.org/content/early/2011/11/14/bioinformatics.btr612.full.pdf

Can be downloaded Here:

http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm

--> This is a front to back full pipeline. However you can see that there are option capabilities to run from any step. Either all, just alignment, just variants, just annotation.

** It's important to note that this pipeline is based on using .BAMS and .vcf as standard file types for modular employment. You'll need to refer to pac-bio or some other source for converters.

ADD COMMENTlink written 7.2 years ago by Raymond301140

*** Note that there is also a Cloud deployment in case you don't have the computational resources to run this pipeline.

ADD REPLYlink written 7.2 years ago by Raymond301140

hi raymond, although the tool looks useful for other purposes i couldnt see how it handles structural variants such as translocations and inversions

ADD REPLYlink written 7.2 years ago by Doctoroots780

You are correct. Upon closer inspection translocations & inversions are not readily annotated through this software. I apologize.

ADD REPLYlink written 7.2 years ago by Raymond301140
0
gravatar for Robert Ernst
7.2 years ago by
Robert Ernst60
Rotterdam, The Netherlands.
Robert Ernst60 wrote:

I have basically the same question. I have to create a SV pipeline as a project on my master. We also have to include some sort of annotation. For now we annotate the output via a GFF file. The GFF contains all the positions of genes, exons, introns etc. We compare these positions with our SV list.

Furthermore we compare SV calls from different software tools to see which SV are supported by multiple SV callers.

I am really wondering what other people's thoughts are about annotating SV's.

ADD COMMENTlink written 7.2 years ago by Robert Ernst60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1486 users visited in the last hour