Question: alternative splicing analysis
1
gravatar for marina-orlova
6.4 years ago by
Russian Federation
marina-orlova80 wrote:

Hello

I have WIG file and need to find genes regulated by the protein which is alternatively spliced (to be precise - the type of alternative splicing should be a5ss - alternative donor sites)

Please advise what programs should I use for this purpose.

 

ADD COMMENTlink modified 4.3 years ago by Biostar ♦♦ 20 • written 6.4 years ago by marina-orlova80
1

Alternative splicing analysis tools require BAM files so that they can use both read abundance across exons and reads spanning splice junctions to identify isoforms. WIG files are for visualization purpose. I don't know of any tools that takes wig files as input and predicts alternative splicing events. Rather, WIG files are often produced as output files from RNA-seq analysis.   

ADD REPLYlink modified 6.4 years ago • written 6.4 years ago by Ashutosh Pandey12k

Thank you for your answer. I found also sam and fastq files for input and chip-seq, hope this is enough. Could you please advise what programs should I use now?

ADD REPLYlink written 6.4 years ago by marina-orlova80
1

You can use MISO (http://genes.mit.edu/burgelab/miso/) that will take the SAM files as input. You can also use Tophat and Cufflinks that will take fastq files as input. I would suggest you to use Tophat+Cufflinks.  

ADD REPLYlink written 6.4 years ago by Ashutosh Pandey12k

Thank you a lot!

ADD REPLYlink written 6.4 years ago by marina-orlova80
4
gravatar for Charles Warden
6.4 years ago by
Charles Warden7.9k
Duarte, CA
Charles Warden7.9k wrote:

MATS is my favorite tool for splicing event analysis

As mentioned in the previous comment, MISO is another option. In both cases, you need .bam files for RNA-Seq data (neither uses ChIP-Seq data, which is probably not a reliable tool for analyzing splicing events).

ADD COMMENTlink modified 13 months ago by _r_am31k • written 6.4 years ago by Charles Warden7.9k

Thank you very much for the answer. I have both CHiP-seq and RNA-seq data. At first I should find peaks which are situated in introns and 3-UTR regions, and then check if the genes where peaks are situated are alternatavely spliced (a5ss). Could you please advise how to connect these two parts of the task?

ADD REPLYlink written 6.4 years ago by marina-orlova80

I think you have the right idea.  I don't know of a specific program to accomplish this all at once.  For the ChIP-Seq peak calling, I assume you can use something standard like MACS.  If you have a list of interesting coordinates, you could use bedtools to look for overlaps between peaks and regions

For an alternative 5` splice site, I would probably look for peaks around promoters where the ChIP-seq binding is upstream of the start of transcription, but maybe you have something different in mind.

ADD REPLYlink written 6.4 years ago by Charles Warden7.9k

Since this got mentioned in a more recent thread, I think I should probably say that I currently like QoRTS + JunctionSeq

ADD REPLYlink written 21 months ago by Charles Warden7.9k
1
gravatar for eddie.im
6.4 years ago by
eddie.im130
Brazil
eddie.im130 wrote:

If you have RNA-seq data you can use TopHat, however it uses fastq/bam/sam files: http://bioinformatics.oxfordjournals.org/content/25/9/1105.short
Or if you already know R language, you can use the Bioconductor package Rsubread, using the "subjunc" function.

ADD COMMENTlink written 6.4 years ago by eddie.im130

Thank you! I will try using both

ADD REPLYlink written 6.4 years ago by marina-orlova80
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: 951 users visited in the last hour