junctions usage

fadista · 07-09-2012, 08:11 AM

Dear all,

I would like to know if there is a tool out there that can handle the junctions file coming out from Cufflinks and treat them as for instance the DEXSeq Bioconductor package, but instead of exon usage it would do junctions usage analysis (or transcript isoform deconvolution).

Thanks.

sudders · 11-28-2012, 03:22 AM

Hi, I've done this using DEXSeq, I just used DEXSeq and treated the junciton counts as exon counts and then did DEXSeq as usual.

I found that the results as WAY more conservative than DEXSeq on exon counts (on exon counts, differentail genes: 426, on junction counts differential genes: 50), but that when I looked at the resuts(in for e.g. IGV) I believed far more of them .

I found that the problem with doing DEXSeq analysis on exon read counts is that almost all my differentially used exons were 3' UTRs becuase the were the longest. Using junctions gets around this because all junctions are the same 'length'.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
samtools psl2sam.pl usage?	ehlin	Bioinformatics	3	09-17-2013 02:38 PM
Tophat memory usage during "Searching for junctions via segment mapping"	biznatch	RNA Sequencing	9	02-18-2013 10:47 AM
Missing Junctions in Tophat! ( after providing a known junctions & gene models files)	avi	Bioinformatics	2	08-02-2011 05:21 PM
dwgsim usage	gprakhar	Bioinformatics	18	07-12-2011 09:59 AM
Memory usage	Seta	General	2	04-05-2011 11:44 AM

07-09-2012, 08:11 AM	#1
fadista Member Location: Malmö Join Date: Sep 2008 Posts: 37	junctions usage Dear all, I would like to know if there is a tool out there that can handle the junctions file coming out from Cufflinks and treat them as for instance the DEXSeq Bioconductor package, but instead of exon usage it would do junctions usage analysis (or transcript isoform deconvolution). Thanks.

11-28-2012, 03:22 AM	#2
sudders Member Location: Sheffield, UK Join Date: Dec 2011 Posts: 32	Hi, I've done this using DEXSeq, I just used DEXSeq and treated the junciton counts as exon counts and then did DEXSeq as usual. I found that the results as WAY more conservative than DEXSeq on exon counts (on exon counts, differentail genes: 426, on junction counts differential genes: 50), but that when I looked at the resuts(in for e.g. IGV) I believed far more of them . I found that the problem with doing DEXSeq analysis on exon read counts is that almost all my differentially used exons were 3' UTRs becuase the were the longest. Using junctions gets around this because all junctions are the same 'length'. Last edited by sudders; 11-28-2012 at 03:25 AM.